CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Simple radiation validation problem (http://www.cfd-online.com/Forums/openfoam-solving/122297-simple-radiation-validation-problem.html)

Logan Page August 16, 2013 08:36

Simple radiation validation problem
 
5 Attachment(s)
Hi All

I am trying to solve a very simple radiation case in order to validate OpenFoam's radiation model/s, so far with little success.

The case I'm trying to solve is a 2D square enclosure (1m x 1m):
  • left wall: uniform T at 100K (black body - emission)
  • right, top and bottom walls: uniform T at 0 K (black body - absorption)
  • no convection or conduction (vacuum)
The analytical results for this are quite easy to calculate and are as follows:
  • Qr from the left wall: 5.67 W/m2
  • Qr into right wall: 2.3474 W/m2
    • 5.67 * ViewFactor = 5.67 * 0.414
  • Qr into top wall: 1.661 W/m2
    • 5.67 * ViewFactor = 5.67 * 0.293
  • Qr into bottom wall: 1.661 W/m2
    • 5.67 * ViewFactor = 5.67 * 0.293
For this I have tried the P1 and fvDOM radiation models, I have tried using the buoyanySimpleFoam solver as well as a custom solver I created and all give me more or less the same incorrect solution.

I have attached the following:

Qr on the leftWall should be uniform across the whole surface equal to 5.67, which I am able to get with the custom solver, but not with the buoyantSimpleFoam solver (for some reason the thermal conductivity in the energy equation has an effect on Qr).
Qr on the rightWall should give a area-weighted sum of 2.3474, in paraFoam I get an area-weighted sum of 0.6103. I am unable to get the correct Qr for the right, top and bottom surfaces.

I hoping that this is the case of me missing something fundamental so any comments / advice would be useful :D

Logan Page August 19, 2013 06:59

Can someone please explain to me how and where these coefficients (in the radiationProperties file) play a role in the radiation models:

Code:

absorptionEmissionModel constantAbsorptionEmission;

constantAbsorptionEmissionCoeffs
{
  absorptivity    absorptivity    [ 0 -1 0 0 0 0 0 ]  0;
  emissivity      emissivity      [ 0 -1 0 0 0 0 0 ]  0;
  E              E              [ 1 -1 -3 0 0 0 0 ] 0;
}

Because if I change the absorptivity and emissivity to 0 (as shown above) i get the exact analytical solution for this test case.

Zeppo September 7, 2013 13:23

Quote:

Originally Posted by Logan Page (Post 446604)
Because if I change the absorptivity and emissivity to 0 (as shown above) i get the exact analytical solution for this test case.

Have you tried out ViewFactors radiation model?

ARTem September 30, 2013 02:31

Hello, Logan Page.
The spherical harmonic approx. method (P1) gives me 8.47 [W/m2] and 2.79 [W/m2] for q_left and q_right correspondingly.
The discrete ordinates method gives me 5.67 [W/m2] and 2.31 [W/m2] for q_left and q_right correspondingly. This result is very close to theoretical one.
I haven't used viewFactors model yet.

If you're still interesting in this stuff, I can share test case here.

olivierG September 30, 2013 03:52

hello,

You can not use P1 model with an empty cavity, since P1 model is for optical thick media (a*l >1), with a= absorptivity and l = carac. length
DO and viewfactor method should work however.

regards,
olivier

Logan Page September 30, 2013 06:33

1 Attachment(s)
Hi

Thanks for the feedback.

I was able to get the theoretical results using the DO method by setting "absorptionEmissionModel" to "none" in the "radiationProperties" file.

I was also able to figure out the theory and 90% of the implementation thereof in OpenFOAM for the DO method through the use of the book by M. Modest (Radiative Heat Transfer, 3rd Edition)


However for the life of me I cannot figure out why there is an additional source radiation term implemented for the DO method in OpenFoam.
For a participating, non-scattering, medium the governing RTE is given by:
http://www.cfd-online.com/Forums/att...1&d=1380536642

However in OpenFOAM there is an additional source term that has been added:
Code:

            IiEq =
            (
                fvm::div(Ji, ILambda_[lambdaI], "div(Ji,Ii_h)")  //<-- first term
              + fvm::Sp(k*omega_, ILambda_[lambdaI])              //<-- second term
            ==
                1.0/constant::mathematical::pi*omega_
              * (
                    k*blackBody_.bLambda(lambdaI)                //<-- third term
                  + absorptionEmission_.ECont(lambdaI)/4          //<-- additional source ??
                )
            );

The first 3 terms implemented in OpenFOAM are 100% correct according to the theory, however I cannot find any theory / literature as to why there would be an additional source term in the RTE.

For a "constantAbsorptionEmission" model "ECont(lambdaI)" is simply the "E" value specified by the user in the "radiationProperties" file. For a "greyMeanAbsorptionEmission" model "ECont(lambdaI)" is "EhrrCoeff * dQ" where again "EhrrCoeff" is specified by the user in the "radiationProperties" file.

wc34071209 February 18, 2014 14:13

Quote:

Originally Posted by ARTem (Post 454193)
Hello, Logan Page.
The spherical harmonic approx. method (P1) gives me 8.47 [W/m2] and 2.79 [W/m2] for q_left and q_right correspondingly.
The discrete ordinates method gives me 5.67 [W/m2] and 2.31 [W/m2] for q_left and q_right correspondingly. This result is very close to theoretical one.
I haven't used viewFactors model yet.

If you're still interesting in this stuff, I can share test case here.

Hello,

could you kindly share your test case ? I am facing the same problem and I want to learn a little about how to set up radiation from your case.

thank you very much.

scintilla March 10, 2014 19:05

Quote:

Originally Posted by Logan Page (Post 454232)
Hi

Thanks for the feedback.

I was able to get the theoretical results using the DO method by setting "absorptionEmissionModel" to "none" in the "radiationProperties" file.

I was also able to figure out the theory and 90% of the implementation thereof in OpenFOAM for the DO method through the use of the book by M. Modest (Radiative Heat Transfer, 3rd Edition)


However for the life of me I cannot figure out why there is an additional source radiation term implemented for the DO method in OpenFoam.
For a participating, non-scattering, medium the governing RTE is given by:
http://www.cfd-online.com/Forums/att...1&d=1380536642

However in OpenFOAM there is an additional source term that has been added:
Code:

            IiEq =
            (
                fvm::div(Ji, ILambda_[lambdaI], "div(Ji,Ii_h)")  //<-- first term
              + fvm::Sp(k*omega_, ILambda_[lambdaI])              //<-- second term
            ==
                1.0/constant::mathematical::pi*omega_
              * (
                    k*blackBody_.bLambda(lambdaI)                //<-- third term
                  + absorptionEmission_.ECont(lambdaI)/4          //<-- additional source ??
                )
            );

The first 3 terms implemented in OpenFOAM are 100% correct according to the theory, however I cannot find any theory / literature as to why there would be an additional source term in the RTE.

For a "constantAbsorptionEmission" model "ECont(lambdaI)" is simply the "E" value specified by the user in the "radiationProperties" file. For a "greyMeanAbsorptionEmission" model "ECont(lambdaI)" is "EhrrCoeff * dQ" where again "EhrrCoeff" is specified by the user in the "radiationProperties" file.

Hi Logan,

I was a bit confused by this when I first saw it too -- as you correctly say, only the first three terms appear in the theoretical RTE. However this is a numerical implementation and I believe the reason for the E parameter in the absorptionEmissionModel is to allow the user the choice of defining the emissive power of the gas directly (W/m3), rather than via the gray gas relation (emissivity * planck function). This would be useful if one wanted to apply a radiative source term that does not vary in proportion to the Planck function.

S

Majed March 10, 2014 21:32

I Need Help
 
1 Attachment(s)
Hello all,

I am new to this forum. I really need your help to solve my problems in running openfoam. attached is my geometry. There are 5 rows of vanes in front of an intake channel. I have produced my geometry, but I think it is not correct. Does any body know, how I can define the water surface in my geometry? I want to consider it as rigid lid. Then my flow in the channel in turbulent, and I want to use k-epsilon turbulent model for the simulations.

Any help is really appreciated.

jrsilvio_ver October 11, 2014 16:28

Dear Friends,
In the example of validation discussed, as you would to include the ViewFactor model? What files should include?
I appreciate everyone's attention.

wc34071209 October 13, 2014 04:49

Quote:

Originally Posted by jrsilvio_ver (Post 513906)
Dear Friends,
In the example of validation discussed, as you would to include the ViewFactor model? What files should include?
I appreciate everyone's attention.

Hi,

I am eager to know as well.

Amit-1911 March 2, 2015 06:37

Dear Logan page,

I am new to the OpenFOAM and hence may sound stupid!
Could you pls tell me how did you manage to set up your problem with no convection and conduction ( Conduction i can understand as no solid body present).
and if I want to consider convection with radiation what necessary changes i have to make??

Thanks in advance!

Regards,
Amit Dhage


All times are GMT -4. The time now is 23:19.