CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Problems with chtMultiRegionSimpleFoam and radiation with fvDOM

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Search this Thread Display Modes
Old   March 7, 2015, 18:20
Default Problems with chtMultiRegionSimpleFoam and radiation with fvDOM
Senior Member
Join Date: Oct 2013
Posts: 337
Rep Power: 21
zfaraday will become famous soon enough
Dear Foamers,

I'm working in a case where a sample is heated inside a cylindrical heater by radiation. Convection should be neglected since the heater is placed within an enclosure in vaccum conditions. The heater is made of a very thin metal (0.1mm thickness) in which an electrical current is injected to make it generate the required thermal power to heat up the sample. Below you can see my geometry.

You can see that I generated it as a cylindrical mesh using blockMesh because all the sistem has cylindrical symmetry.

As I don't need to compute anything in the air but the energy transference (only by radiation) I used frozenFlow option in the SIMPLE subdictionary within fvSolution file, so that no velocity nor pressure are computed to speed up the simulation.

As I said, the heater is really thin (0.1mm) and for this reason at the very beginning I tried to create it as a baffle with not much success. I created the shields as a baffle too (2mm thickness) but it all gave me a lot of trouble and finally I decided to create the thinner regions as solid elements. To simulate the power generated by the heater I used the fvOptions framework.

I ran some simulations and the results I got, when it comes to temperature distributions, are quite satisfying: the ranges of temperatures are very similar to the ones obtained with another software for heat transfer. Maybe the distributions are not exactly the same, but the mean values are quite close compared with the results I first obtained when I created the thinner regions as baffles... That made me a little glad!

However, not totally satisfied with that, I used som function objects and wallHeatFlux utility in order to compute all heat fluxes in the air region (some fluxes in the solids also) to check heat balance. And here comes the problem... Even though my case has converged (the solver is not iterating for any region for the last few time steps), no balance is gotten when it comes to energy.

I computed radiative power in all patches belonging to the air region as
Thus, if I sum all of the values gotten from all the patches I should get a value of 0 (energy cannot be created nor destroyed!) if I'm not wrong... Well, then I don't get any 0 but a value around 6-9 depending on the case. May anyone explain me why is it happening to me?

Not only the balance in the air is not balanced, but also is unbalanced for the solid regions, some of them are more or less balanced, but others have a magical gain or a loss of energy. Here you can see the balance for the sample as an example (extracted from the wallHeatFlux utility):
Wall heat fluxes [W]
sample_to_air 0.05049047
sample_to_support -0.0600329
Another point I don't understand is the following: If I set the value of the generation in the fvOptions file as, say, 50W, am I not supposed to get the same value when I compute the power in the post-processing stage? I ask that because it's not happening to me. As an example I share some numbers with you. I set the generation in one of the tests I've done as 16W in the fvOptions. After the simulation I got these values with wallHeatFlux and function objects:
  • wallHeatFlux
    Wall heat fluxes [W]
    heater_to_air -17.56851
  • function object to compute radiative power
    Expression Qr_air_patches on air_to_heater:  sum=-17.56811

Why is not the same value of power I set as a generation?

At this stage, I'm a little lost because I ran a few simulations using different values but I always got an unbalanced system. I don't know if the problem may be related to the geometry and the thin regions I need to use for this simulation or maybe it is something regarding the fvDOM specification (since I'm not very experienced in the use of this method). I would appreciate any word you can give me regarding this matter since I can't go further on my own...

Many thanks in advance.

Best regards,

Attached Images
File Type: jpg heater.jpg (19.3 KB, 28 views)
File Type: jpg heater_thickness.jpg (36.9 KB, 24 views)
Web site where I present my Master's Thesis:

The case I talk about in this site was solved with chtMultiRegionSimpleFoam solver and involves radiation. Some basic tutorials are also resolved step by step in the web. If you are interested in these matters, you are invited to come in!

Last edited by zfaraday; March 7, 2015 at 18:28. Reason: extra picture of the heater attached
zfaraday is offline   Reply With Quote

Old   February 3, 2018, 20:28
Join Date: Jun 2015
Posts: 34
Rep Power: 10
jaydeep is on a distinguished road
too little too late from me
jaydeep is offline   Reply With Quote


chtmultiregionsimplefoam, fvdom, radiation, wallheatflux

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Problems using fvdom radiation model and convection heat transfer jrsilvio_ver OpenFOAM Running, Solving & CFD 3 January 22, 2015 14:01
Radiation Modeling Chris89 CFX 20 August 14, 2014 08:51
GETVAR Error in Multiband Monte Carlo Radiation Simulation with Directional Source silvan CFX 3 June 16, 2014 10:49
Modeling both radiation and convection on surfaces - Ansys Transient Thermal R13 s.mishra ANSYS 0 March 31, 2012 05:12
radiation model fvDOM thomasduerr OpenFOAM Running, Solving & CFD 0 February 25, 2010 04:11

All times are GMT -4. The time now is 01:53.