|
[Sponsors] |
September 15, 2016, 13:29 |
Radiation in chtMultiRegionSimpleFoam
|
#1 |
New Member
Mauricio Aguilar
Join Date: Jul 2016
Posts: 15
Rep Power: 10 |
Hello Foamers,
I have a couple of questions to those that have experience with radiation modeling. I have set up a simple case to validate my simulation, based on the multiregionheaterRadiation tutorial. I have a couple of parallel solid plates with a non participating fluid (air) the middle. Plate 1 has a fixed T of 294K. Plate 2 has a fixed T of 273. Both measure .1 by .1 m. And the fluid thickness is of .01m The radiation model in the Air is viewFactors and in the plates is opaqueSolid. My interpretation of the tutorial setup is that the Qr solved by the viewFactors radiation model is taken from the solid radiation inherited from each plate. Question 1: in the /0/Air/viewFactorsField i can see values greater than 1. Why is this happening? Analytically calculating the view factor for the entire plate I get a value of 0.82 using the equation in Incropera's Fundamentals of Mass and Heat Transfer. Shouldn't the greatest value of a middle cell whose radiation is completely intercepted be limited to 1. Question 2: When I set the emissivity of both solids to 1 and their absorptivity to 0 (of an opaque body) I get a net radiative flux , Qr or q'', that corresponds to the analytical solution. Simulations Qr = 106Wm^-2; Analytical Qr = 108 Wm^-2. However if I set the emissivity on the cold plate to 0.5 the Qr is halved. Taking the net radiative flux (Qr) to be equal to the radiosity (J) minus the Irradiation (G) , Qr = J -G , the results should be different as a change in emissivity in the cold plate should only reduce the Irradiation on the hot side and thus increase the magnitude of heat flux on that plate. How should I set up these boundary conditions so that I can have solids with different emissivities? I'll keep trying to figure things out in the mean time but I hope someone will know how to approach this. kind regards, Mauricio |
|
September 16, 2016, 14:28 |
|
#2 |
New Member
Mauricio Aguilar
Join Date: Jul 2016
Posts: 15
Rep Power: 10 |
Question 2:
Total brain fart and I apologize for it. I forgot that I was dealing with opaque solids and non absorbed energy bounces back(increasing the returning G value ) and satisfies the conservation of energy. Question 1: Still uncertain about the view factor matrix and will increase the number of coarse faces to evaluate differences. |
|
September 17, 2016, 07:37 |
|
#3 | |
Senior Member
Sergei
Join Date: Dec 2009
Posts: 261
Rep Power: 21 |
Quote:
|
||
September 19, 2016, 07:57 |
|
#4 |
New Member
Mauricio Aguilar
Join Date: Jul 2016
Posts: 15
Rep Power: 10 |
Sergei,
I just read your thread and the error % you are getting is huge! My error is limited to a 10% discrepancy by comparison to my expected results. If I cannot solve this by increasing the number of coarse faces I may just parse the viewfactor matrix and replace any value>1 with 1. |
|
September 19, 2016, 10:24 |
|
#5 |
Senior Member
Join Date: Sep 2013
Posts: 353
Rep Power: 21 |
Code:
viewFactorCoeffs { smoothing true; //Smooth view factor matrix (use when in a close surface //to force Sum(Fij = 1) ... } Did you use this option? |
|
September 23, 2016, 12:00 |
|
#6 | |
Senior Member
Sergei
Join Date: Dec 2009
Posts: 261
Rep Power: 21 |
Quote:
Yes. And, in fact, it didn't make any difference. |
||
September 29, 2016, 13:27 |
|
#7 |
New Member
Mauricio Aguilar
Join Date: Jul 2016
Posts: 15
Rep Power: 10 |
I've been messing around with radiation lately and I found that values greater than 1 in the ViewFactorField dictionary must get parsed later on, because I wrote a simple program to replace values greater than 1 into 1 and it didn't make much difference.
In the setup I've been working with involving parallel plates at 373 and 273 K with 0.07 emissivity, when spaced at 0.01m the viewFactors for the hot and cold plates revolved around 8, changed them to 1 and got like a 0.02% difference in the resulting heat flux. However, when the viewFactorGen fails so does the conservation of energy. With very small gaps, as in the previous case I get an unbalanced energy transfer between all surfaces. As the gap distance increases and the values in the viewFactorField approach 1 or less the net radiative transfer approaches 0. Although my residuals show convergence this post-processing on the heat flux clearly shows otherwise. |
|
Tags |
chtmultiregionsimpefoam, emissivity, radiation, viewfactor |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
chtMultiRegionSimpleFoam: Thermal Conduction + Surface-To-Surface Radiation | Zeppo | OpenFOAM Running, Solving & CFD | 16 | May 18, 2017 19:04 |
Monte Carlo Simulation: H-Energy is not convergating & high Incident Radiation | volleyHC | CFX | 5 | April 3, 2016 06:41 |
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 |