
[Sponsors] 
April 15, 2014, 17:10 
Lack of knowledge simulating radiation in solids

#1 
Senior Member
Alex
Join Date: Oct 2013
Posts: 286
Rep Power: 13 
Hello everybody!
I'm trying to solve a problem using chtMultiRegionSimpleFoam involving radiation between some solids surrounded by air as a nonparticipating media. The case is explained in this post: Problems in my first attempt with chtMultiRegionSimpleFoam. There you can see a picture of my geometry. I am using fvDOM for the calculation of the radiation, but I am a bit confused because there is no multiregion tutorial using fvDOM for the radiation calculation. Thus, I am trying to set it up by following the multiRegionHeaterRadiation and the hotRoomRadiationFvDOM tutorials. Now I need some help and advices to know If I am doing the radiation set up right. 1) Since air is a nonparticipating media I used this radiationProperties file (taken from the hotRoomRadiationFvDOM tutorial and adapted to my case): Code:
radiationModel fvDOM; fvDOMCoeffs { nPhi 3; // azimuthal angles in PI/2 on XY.(from Y to X) nTheta 5; // polar angles in PI (from Z to XY plane) convergence 1e3; // convergence criteria for radiation iteration maxIter 5; // maximum number of iterations } // Number of flow iterations per radiation iteration solverFreq 10; absorptionEmissionModel constantAbsorptionEmission; constantAbsorptionEmissionCoeffs { absorptivity absorptivity [ 0 1 0 0 0 0 0 ] 0.0; emissivity emissivity [ 0 1 0 0 0 0 0 ] 0.0; E E [ 1 1 3 0 0 0 0 ] 0; } scatterModel none; 2) On the other hand, the solid regions set up is what is confusing me more. This is the radiationProperties file I used for one of the solid regions (Also taken from the tutorial and adapted): Code:
radiationModel opaqueSolid; absorptionEmissionModel constantAbsorptionEmission; constantAbsorptionEmissionCoeffs { absorptivity absorptivity [ 0 1 0 0 0 0 0 ] 0.0; //opaque emissivity emissivity [ 0 1 0 0 0 0 0 ] 0.5; E E [ 1 1 3 0 0 0 0 ] 0; } scatterModel none; Then, according to Kirchhoff's law we can say that (only valid for black body, in other cases they are equal but different to 1) (so I maybe should define the same value for both coefficients). But what I don't really understand is that opacity is defined by absorptivity (alpha) in the tutorial, when, as per my knowledge in matters of radiation, opacity is defined by the transmissivity (thau), being 0 when the solid is opaque. Once I found out this little inconsistency (at least to my eyes, I guess I must have missed something or I am mistaken in something) I sought all along the english and spanish literature to see if it was a matter of a bad translation from spanish (the language I studied it from originally) to english, but I couldn't find anything clear. 3) Finally, the last question I have is about the units of both coefficients, that is, m⁻đ. Shouldn't it be dimensionless? As you can see, I'm a little lost with this simulation and I really need a little help to see what I am doing wrong. Any hint or advice will be welvcome and much appreciated. Many thanks in advance!! Alex
__________________
Web site where I present my Final Thesis developed with OpenFOAM: foamingtime.wordpress.com The case I talk about in this site was solved with chtMultiRegionSimpleFoam solver and involves radiation. If you are interested in these matters, you are invited to come in! I'm newbie in OpenFOAM's world and not an Englishspeaking, so if I make any mistake a correction will be welcome! Last edited by zfaraday; April 16, 2014 at 07:27. Reason: Edited a mistake, clarified in green within brackets 

April 20, 2014, 14:48 

#2 
Senior Member
Alex
Join Date: Oct 2013
Posts: 286
Rep Power: 13 
I don't know if my explanation was correct enough or not, but i really need a little help to understand the radiation set up process in order to solve my problem. Any hint will be much appreciated!
Alex
__________________
Web site where I present my Final Thesis developed with OpenFOAM: foamingtime.wordpress.com The case I talk about in this site was solved with chtMultiRegionSimpleFoam solver and involves radiation. If you are interested in these matters, you are invited to come in! I'm newbie in OpenFOAM's world and not an Englishspeaking, so if I make any mistake a correction will be welcome! 

May 9, 2014, 06:03 

#3 
Senior Member
Join Date: Oct 2013
Posts: 321
Rep Power: 6 
I'm currently also looking at radiation modelling with fvDOM, using a participating plasma with a multigroup approach (wideBandAbsorptionEmissionModel). I can't help you too much yet I'm afraid, but absorptivity and emissivity (the properties of the medium), are in m⁻1 because it's the factor in exponential decay of intensity absorped in a unit length.


May 10, 2014, 14:39 

#4 
Senior Member
Alex
Join Date: Oct 2013
Posts: 286
Rep Power: 13 
Thank you so much chriss85 for your replay, that clarifies a little. Then, according to what you said I understant that if I want to work in a non participating media I have to turn these coefficients to 0, haven't I? Or, otherwise, switch absorptionEmissionModel to none.
Thanks again. Alex
__________________
Web site where I present my Final Thesis developed with OpenFOAM: foamingtime.wordpress.com The case I talk about in this site was solved with chtMultiRegionSimpleFoam solver and involves radiation. If you are interested in these matters, you are invited to come in! I'm newbie in OpenFOAM's world and not an Englishspeaking, so if I make any mistake a correction will be welcome! 

May 12, 2014, 03:15 

#5 
Senior Member
Join Date: Oct 2013
Posts: 321
Rep Power: 6 
Yes, that's correct I think. If your radiation comes from surfaces only and the medium between doesn't absorb or emit, you can turn these values to 0.


May 12, 2014, 13:34 
radiation model in solids

#6  
Senior Member
Alex
Join Date: Oct 2013
Posts: 286
Rep Power: 13 
Thank you for your help chriss85, now I understand a little bit more the matter of radiation. But I still keep having a doubt about it, specifically this is about the radiation in solids as I asked in the second question on the first post of this thread. Here comes the question again just in case you can help me understanding it...
Quote:
__________________
Web site where I present my Final Thesis developed with OpenFOAM: foamingtime.wordpress.com The case I talk about in this site was solved with chtMultiRegionSimpleFoam solver and involves radiation. If you are interested in these matters, you are invited to come in! I'm newbie in OpenFOAM's world and not an Englishspeaking, so if I make any mistake a correction will be welcome! 

May 13, 2014, 03:00 

#7 
Senior Member
Join Date: Oct 2013
Posts: 321
Rep Power: 6 
I haven't worked with the opaqueSolid model yet, so I can't really tell you anything about it. I suggest that you also take a look into the code of it, to see what equations are being solved. Some comments may also help you.
Regarding the absorptivity, I agree that it should not be 0 in an opaque solid, but rather a very high value instead, because the absorption and the exponential decay of intensity will happen on a very short length scale. I think transmittivity can only be defined per length, while transmittance is the absolute value related to a specific length. Don't quote me on that though, it's been a while for me. Reflectivity is only useful on surfaces, not in solids. If you look at a solid of a fixed size, then your equation saying alpha + rho + tau = 1 is correct, but it refers to the absolute ratios of intensity (i.e. reflectance, transmittance, and absorption), not on absorptivity/emissivity, which are per length. 

May 19, 2014, 11:47 

#8 
Senior Member
Alex
Join Date: Oct 2013
Posts: 286
Rep Power: 13 
Thank you criss85 for your advice. Now I have taken a look into the code (my knowledge in C++ is more than weak, though) another doubt came out to my mind. This is what the solver solves for the solid region (solveSolid.H):
Code:
fvScalarMatrix hEqn (  fvm::laplacian(betav*alpha, h, "laplacian(alpha,h)") == fvOptions(rho, h) ); hEqn.relax(); fvOptions.constrain(hEqn); hEqn.solve(); fvOptions.correct(h); Code:
fvScalarMatrix EEqn ( fvm::div(phi, he) + ( he.name() == "e" ? fvc::div(phi, volScalarField("Ekp", 0.5*magSqr(U) + p/rho)) : fvc::div(phi, volScalarField("K", 0.5*magSqr(U))) )  fvm::laplacian(turb.alphaEff(), he) == rad.Sh(thermo) + fvOptions(rho, he) ); Code:
Description Radiation for solid opaque solids  does nothing to energy equation source terms (returns zeros) but creates absorptionEmissionModel and scatterModel. Thanks in advance! Alex
__________________
Web site where I present my Final Thesis developed with OpenFOAM: foamingtime.wordpress.com The case I talk about in this site was solved with chtMultiRegionSimpleFoam solver and involves radiation. If you are interested in these matters, you are invited to come in! I'm newbie in OpenFOAM's world and not an Englishspeaking, so if I make any mistake a correction will be welcome! 

January 20, 2015, 07:01 

#9 
New Member
Pedro
Join Date: Dec 2014
Posts: 9
Rep Power: 2 
Hi Alex,
Where you able to discover exactly how to model radiation between the two solids? The tutorial seems to be only applicable in the case where the solid is not absorbing any of the radiation, is this correct? My question more in a modelling sense is: Letīs say that you have two bodies if you model them both as opaque and you take a radiation model for air. Does this mean that air will absorb all the radiation? 

January 20, 2015, 08:48 

#10  
Senior Member
Alex
Join Date: Oct 2013
Posts: 286
Rep Power: 13 
Hello Pedro,
Quote:
Quote:
Quote:
Hope it helps. Alex
__________________
Web site where I present my Final Thesis developed with OpenFOAM: foamingtime.wordpress.com The case I talk about in this site was solved with chtMultiRegionSimpleFoam solver and involves radiation. If you are interested in these matters, you are invited to come in! I'm newbie in OpenFOAM's world and not an Englishspeaking, so if I make any mistake a correction will be welcome! 

Tags 
fvdom, radiation 
Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Modeling both radiation and convection on surfaces  Ansys Transient Thermal R13  s.mishra  ANSYS  0  March 31, 2012 04:12 
radiation between solids in CFX  groonm  Main CFD Forum  0  September 15, 2011 17:28 
Multiphase radiation in CFX  Luk_Fiz  CFX  5  March 17, 2009 04:51 
Thermal radiation with solids rotating  luis  CFX  0  March 16, 2008 22:04 
code for simulating radiation heat transfer  Seaspace  Main CFD Forum  1  October 3, 2002 13:16 