[upuli]

Dear members
I am using OpenFoam version 4.0 and want to add radiation heat transfer in to my temperature equation. below is the equation

{
alphat = turbulence->nut()/Prt;
alphat.correctBoundaryConditions();

volScalarField alphaEff("alphaEff", turbulence->nu()/Pr + alphat);

fvScalarMatrix TEqn
(
fvm:T)
+ fvm::div(phi, T)
- fvm::laplacian(alphaEff, T)
==
radiation->ST(rhoCpRef, T)
+ fvOptions(T)
);

TEqn.relax();

fvOptions.constrain(TEqn);

TEqn.solve();

radiation->correct();

fvOptions.correct(T);

rhok = 1.0 - beta*(T - TRef);
}
I included the
radiation.H in the solver and canged the make/options file also.

But when I compile the solver it gives errors that rhoCpRef was not declared.

I will be thankfull if you can help me to solve this error.

rgds

K.U.C.Perera

[DustExplosion]

Hi Upuli,

I am not sure about your specific problem, but you may find this reference useful for implementing radiation in an OF solver: http://www.tfd.chalmers.se/~hani/kur...Foam_final.pdf