Radiation Model

julienh

Hi,
I am trying to program a radiation model using the Gibb's method (spherical harmonics)

I used the "reactingFoam" solver, where I had a file IEqn.H with the following radiation equation :

{

fvScalarMatrix IEqn

(

fvm::laplacian(I)
- (I/(4*boltzmannCoeff)
+ pow(T,4))*(12*pow(knu,2 *boltzmannCoeff)

==
m

);

IEqn.relax();

IEqn.solve();

thermo->correct();

}

where I is the radiativ intensity I wanted to get from this equation

After the compilation, when I execute reactingFoam, I get this message :


--> FOAM FATAL ERROR : gradientInternalCoeffs cannot be called for a calculatedFvPatchField.

You are probably trying to solve for a field with a calculated or default boundary conditions.

From function calculatedFvPatchField<type>::gradientInternalCoef fs() const in file fields/fvPatchFields/basicFvPatchFields/calculated/calculatedFvPatchField.H at line 174.


In order to solve the I equation, I need of course some boundary conditions. I though the file I in the "0" directory of the case with the initial and boundary conditions was enough but it is apparently not the case.

I try to see in the code an equation file like UEqn.H or hEqn.H but I did not see how the boundary equations are coded

Could somebody give me some advices please ?

Thank you in advance

Julienh