void lcsRhoEffPluginFunction::doEvaluation()
{
autoPtr<volScalarField> prhoEff;
// pick up the first fitting class
castAndCall(prhoEff,volScalarField,basicKinematicCloud,kinematicCloud,rhoEff());
castAndCall(prhoEff,volScalarField,swakFluidThermoCloudType,thermoCloud,rhoEff());
castAndCall(prhoEff,volScalarField,basicReactingCloud,reactingCloud,rhoEff());
castAndCall(prhoEff,volScalarField,basicReactingMultiphaseCloud,reactingMultiphaseCloud,rhoEff());
noCloudFound(prhoEff);
result().setObjectResult(prhoEff);
}
tmp<Foam::volScalarField>
continuousGasKEpsilon<BasicTurbulenceModel>::nuEff() const
{
volScalarField blend
(
max
(
min
(
(this->alpha_ - scalar(0.5))/(alphaInversion_ - 0.5),
scalar(1)
),
scalar(0)
)
);
return tmp<volScalarField>
(
new volScalarField
(
IOobject::groupName("nuEff", this->U_.group()),
blend*this->nut_
+ (1.0 - blend)*rhoEff()*nutEff_/this->transport().rho()
+ this->nu()
)
);
}
