CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Radiation transfer equation problem (https://www.cfd-online.com/Forums/openfoam-solving/199487-radiation-transfer-equation-problem.html)

kane March 7, 2018 08:52

Radiation transfer equation problem
 
Hello!!

I need some help in setting up radiation transfer equation.

As you know, there is rte equation in radiativeIntensityRay.C

Code:

forAll(ILambda_, lambdaI)
    {
        const volScalarField& k = dom_.aLambda(lambdaI);

        const surfaceScalarField Ji(dAve_ & mesh_.Sf());

        fvScalarMatrix IiEq
        (
            fvm::div(Ji, ILambda_[lambdaI], "div(Ji,Ii_h)")
          + fvm::Sp(k*omega_, ILambda_[lambdaI])
        ==
            1.0/constant::mathematical::pi*omega_
          *(
                // Remove aDisp from k
                (k - absorptionEmission_.aDisp(lambdaI))
              *blackBody_.bLambda(lambdaI)

              + absorptionEmission_.E(lambdaI)/4
            )
        );

I know that these codes mean 'Radiation transfer equation'.
http://www.cfd-online.com/Forums/att...1&d=1380536642


But I don't know what each term means.
Could you explain meaning of each term in detail?
Or any links which explain radiation transfer equation will be fine.

My final goal is setting up radiation transfer equation including scattering.

Can you please help me in this regards?

Best regards

kane March 9, 2018 00:07

greeting all.

I have tried to undestand radiation transfer equation.
I understand second term and third term in rte.

But I can't understand the first term and the final additional source.

Would you explain the first term and final additional source??

Additionally, why the source(1.0/constant::mathematical::pi*omega_) does exist?

Thanks in advance for your help

kane March 12, 2018 23:36

Nobody know??

why the term(1.0/constant::mathematical::pi*omega) is added??


All times are GMT -4. The time now is 02:54.