[OPFO]

I am using OpenFOAM-v2006 and the icoreactingmultiphaseinterfoam solver with LaserDTRM radiation model and the constantScatter model.

According to [1], in the constantScatter model parameters of radiation properties, the user can input sigma and C constants that ultimately affect a parameter called SigmaEff(). In the P1 radiation model, SigmaEff() is used to alter mesh values to affect diffusion of properties such as T.

// Construct diffusion
const volScalarField gamma
(
IOobject
(
"gammaRad",
G_.mesh().time().timeName(),
G_.mesh(),
IOobject::NO_READ,
IOobject::NO_WRITE
),
1.0/(3.0*a_ + sigmaEff + a0)
);

// Solve G transport equation
solve
(
fvm::laplacian(gamma, G_)
- fvm::Sp(a_, G_)
==
- 4.0*(e_*physicoChemical::sigma*pow4(T_)) - E_
);

In LaserDTRM, I cannot find where SigmaEff() is incorporated into the solver as in the example of P1. For various sigma and C values that I plug in, the simulation is unchanged which would suggest that the scatter model is not actually incorporated into the solver. I have also tried varying the size of my mesh elements but with no effect.

Can someone help explain how the constantScatter model works for the LaserDTRM radiation model?

Thanks in advance!

[1] http://www.tfd.chalmers.se/~hani/kur...Foam_final.pdf

[mkh_sakr]

I also face the same concern about the scatter model. What I think is, the scatter model is not applicable to the icoReactingMultiPhaseInterFoam solver. It could be defined in the radiationProperties file, due to the solver includes the models from the radiation model.