[kuria]

Hey Foamers

I am looking at the reactingFoam solver in OpenFOAM 8 and trying to understand how the energy and species equations are solved.

In energy equation:
fvScalarMatrix EEqn
(
fvm::ddt(rho, he) + mvConvection->fvmDiv(phi, he)
+ fvc::ddt(rho, K) + fvc::div(phi, K)
+ (
he.name() == "e"
? fvc::div
(
fvc::absolute(phi/fvc::interpolate(rho), U),
p,
"div(phiv,p)"
)
: -dpdt
)
+ thermophysicalTransport->divq(he)
==
reaction->Qdot()
+ fvOptions(rho, he)
);

EEqn.relax();

fvOptions.constrain(EEqn);

EEqn.solve();

thermophysicalTransport->divq(he) is calculated for the turbulent part as
template<class TurbulenceThermophysicalTransportModel>
tmp<fvScalarMatrix>
eddyDiffusivity<TurbulenceThermophysicalTransportM odel>::divq
(
volScalarField& he
) const
{
return -fvm::laplacian(this->alpha()*this->alphaEff(), he);
}

Similarly for species transport equation, the laplacian term is calcualted based on
template<class TurbulenceThermophysicalTransportModel>
tmp<fvScalarMatrix>
eddyDiffusivity<TurbulenceThermophysicalTransportM odel>::divj
(
volScalarField& Yi
) const
{
return -fvm::laplacian(this->alpha()*this->DEff(Yi), Yi);
}


I understood until here, but I could not figure out how both DEff and alphaEff are calculated.

Has someone looked into how these are calculated ?