![]() |
Difference between IDefault and G (radiation)
Hi,
Sorry but I am a new OpenFOAM user so maybe it is an easy question. When the radiation is activated, we must have two files of boundary condition: G and IDefault. I understand that G gives the Marshak boundary condition to the walls (usually). Then, IDefault gives greydiffusivewalls (in tutorials). I do not understand that. It appears to be redundant to me, we already give a radiation boundary condition. So I just would like to know the meaning of this two files. I think there is something that I do not understand. Thanks. |
Nobody knows ?
|
G is used for P1 model, IDefault for fvDOM model and Qr for viewFactor model.
|
Thank you for the response.
|
Hi, from my understanding, G is incident radiation and iDefault is radiation intensity ray which calculates radiation intensity.
G= Qin + Qem so if you think about it this way, G is the total heat that's been absorbed and emitted within a media and IDefault which i believe is used for calculating Qin, which stands for incident radiation heat flux. not too sure but you can easily find out using the search box here, its pretty effective. http://www.openfoam.com/documentatio...b34c63f9ffab05 That's IDefault https://github.com/OpenFOAM/OpenFOAM...IntensityRay.C |
Hello,
Actually, Qr=Qin-Qem where Qem is emissivity*Sigma*T⁴ I have checked it with an easy case. Qin is the total incident radiative heat flux field. G is used with the P1 model. But, even if I use the fvDOM method, OpenFoam calculates G. G is incident radiation. Actually, we have Qr=-1/(3a)*grad(G) (if we do not have scattering). So I do not understand one thing. Why do we call G the incident radiation ? Qin and G are different but their name seem to be the same for me. Both are in W/mē and it is what the wall receives. If someone can explain what is the physical difference between G and Qin, I would be grateful. |
Hello Fabien,
hmm interesting: well looking at the source code, it shows that Qr, Qin and Qem are independently calculated: Quote:
to answer your question specifically G accounts for the black body radiation and is the summation of all radiation contribution from all rays and can be calculated my multiplying the intensity of the ray by its solid angle. so once G is updated, Qr, is calculated from G. and I as you clearly stated in your last reply, the Qem and Qin are then calculated. However, i would reconfirm this with an expert like Bruno santos if i were you |
Well I think you are right concerning the calculations of each parameter.
But my last question is what is the difference between Qin and G ? They have similar names but it is not the same thing. What is the power received by the wall for example, Qin or G ? Thank you for your responses. |
they are completely different!!!!
if you try to run a simple case, you will realize that Qin is only calculated at the boundary and G is calculated on the boundaries and throughout the region itself. Qin is the energy generated at the surface. for example, if you have a black body under the sun, Qin will be the heat flux on the intercepting surface, which will be conducted by the regions bellow the surface. in practical terms, Qin is calculated based on absorptivity value. G, like I said earlier, is: Quote:
I will confirm all these details with an expert like Bruno Santos and a few others. |
Hi,
i am wondering where you got this from? is it based on your results? Quote:
|
Hello,
here: https://www.sharcnet.ca/Software/Flu...ug/node577.htm It is the user guide of fluent. We use the same boundary condition (Marshak) in OpenFOAM so I guess G is the same for both softwares. |
1 Attachment(s)
In my opinion the calculation of G in the fvDOM model is wrong. At our institute we developed our own and improved view factor model and work with it quite a lot. We also calculate and write the G-field there and so I compared both models in a test case. I found out that the G-field differs exactly by a factor of 4. I digged in the source code and found the calculation of G:
Code:
G_ += IRay_[rayI].I()*IRay_[rayI].omega(); Do you agree with me and should I write a bug report? I never did this so far |
|
|
Exactly. It is due to the fact that G is used for the energy balance of a single ray and not of a physical surface. The programmer redefined it a little, but for me knowing G from my lectures and books, it was really confusing.
|
All times are GMT -4. The time now is 01:29. |