[mabinty]

dear all!!

i would like to implement the P1 radiation model (used in the buoyantSimplRadiationFoam solver) into the chtMultiRegionFoam solver. i m using the 1.5 version and started first by implementing the P1-model into the transient buoyantFoam solver which i think works well as a parameter study with respect to the absorption/emission coefficients showed plausible results. Now i m about to realize the implementation into the chtMultiRegionFoam solver, but the same parameter study indicates that the energy balance is decoupled from the radiative transfer equation (plausible behavior of the radiative flux but no changes in temperature distribution). my implementation looks like this:

1.) create a radiationModel "rad" in creatFluidFields.H with the temperature field thermof[i].T().

2.) the object "rad" is given to the function solveEnthalpyEquation, done at the top of hEqn.H and solveEnthalpyEquation.C.

3.) add the radiation source term to the enthalpy equation in solveEnthalpyEquation.C in the form rad.Sh(thermo).

i ll keep working on it and report, and would greatly appreciate any comments!! thx in advance!!

aram