CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

ExternalWallHeatFluxTemperature with radiation + convection not working

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 30, 2024, 23:40
Question ExternalWallHeatFluxTemperature with radiation + convection not working
  #1
New Member
 
Theo Firelli
Join Date: Aug 2024
Posts: 4
Rep Power: 2
Engineers_Firelli is on a distinguished road
Hello OpenFOAM users,


I'm new on the forum. I'm trying to simulate a hot gas flow through ducts in OpenFOAM. I need to setup a boundary condition that combines:
* a heat transfer coefficient + outer temperature (representing the convection between the external wall and the ambient)
* an emissivity at the outer wall + outer temperature (representing the radiation from the externall wall to the ambient)


(* perhaps it would be nice also to include the conduction in the wall, but this may be calculated analytically and included in the heat transfer coefficient mentioned above)


I have struggled a lot with externalWallHeatFluxTemperature. I am simulating a simple case with a cylindrical tube (1 inlet, 1 outlet, 1 wall). The inlet is hot fluid (specified uniform velocity), and at the outlet there is a defined pressure (1 bar).



I am using externalWallHeatFluxTemperature for the wall.

walls
{
type externalWallHeatFluxTemperature;
mode coefficient;

//kappaMethod fluidThermo;
Ta constant 298.15;
h uniform 0.01;

emisssivity 0.95;
//thicknessLayers List<scalar> 1(0.1);
//kappaLayers List<scalar> 1(5.0);
//qr none;
//relaxation 0.3;
//qrRelaxation 0.3;
value $internalField;
}




I tried many different combinations of the parameters, but it seems that setting "emissivity" has no effect. Even, if I check the results file at some timestep different from "0", the emissivity does not even show up !


walls
{
type externalTemperature;
h uniform 0.01;
Ta
{
type constant;
value 298.15;
}
refValue uniform 298.15;
refGradient uniform 0;
valueFraction uniform 4.18458e-06;
value nonuniform List<scalar>
}


I looked into the code for externalWallHeatFluxTemperature for mode coefficient and I saw the coefficient there is defined as I would expect: h_rad = (Tp^4 - Ta^4)/(Tp-Ta) is added to the h defined by the user.



However I cannot make it function to compute radiation as well. It only uses the h that I provide.




I also tried to include a radiation model in my simulation (P1) (although I actually don't want to simulate the radiation inside the fluid domain, JUST to calculate the radiation from the external wall of the duct to the ambient), but it didn't change this behaviour.


I checked different related posts on the forum, but none seemed to give a good solution to this.



Has anyone had success in implementing convection + radiation (+conduction) using externalWallHeatFluxTemperature ? Did it work as you had expected ? Was the radiation contribution to the heat loss from the wall computed correctly ?



I am using OpenFOAM 11, the solver is fluid (and I noticed the same behaviour with CHTMultiRegion - now called foamMultiRun).




Thanks in advance.
Theo
Engineers_Firelli is offline   Reply With Quote

Old   September 3, 2024, 12:26
Default With the ESI versions (tested 2306 and 2406) externalWallHeatFluxTemperature works !!
  #2
New Member
 
Theo Firelli
Join Date: Aug 2024
Posts: 4
Rep Power: 2
Engineers_Firelli is on a distinguished road
I have tested the same case with OpenFOAM 2306 and 2406, and it worked exactly as expected. Even setting h = 0.000001 but with emissivity 0.99, a large amount of heat was released from the flow of hot gas, so the gas cooled down significantly.


simpleCase_externalWall_Radiation.zip

Screenshot_tube_longitudinal_section.jpg

In OpenFOAM 11 (solver "fluid"), the emissivity did not make any difference and in the same case as mentioned above, the gas stayed hot throughout the entire duct !


wall
{
type externalWallHeatFluxTemperature;
mode coefficient;
Ta constant 290.0;
h uniform 1e-06;
emissivity 0.99;
kappaMethod fluidThermo;
value $internalField;
}




Any idea where this difference is coming from ?

Last edited by Engineers_Firelli; September 3, 2024 at 12:35. Reason: added case files
Engineers_Firelli is offline   Reply With Quote

Reply

Tags
boundary condition, emissivity, externalwallheatfluxtemp, radiation

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
Radiation in semi-transparent media with surface-to-surface model? mpeppels CFX 11 August 22, 2019 08:30
Heat transfer: Sim. convection incl. radiation from power consumption on circuitboard maett OpenFOAM Running, Solving & CFD 1 December 1, 2017 13:40
Simulation of natural convection and radiation of a radiant heater in a room ami236 FLUENT 0 July 5, 2017 22:51
Heat Transfer - Convection and Radiation dyku OpenFOAM Running, Solving & CFD 3 September 3, 2015 00:22
CFD code for Convection + Radiation Rodrigo Lucianetti Main CFD Forum 4 December 23, 1998 16:33


All times are GMT -4. The time now is 18:16.