void lcsEnthalpySourcePluginFunction::doEvaluation()
{
typedef DimensionedField<scalar,volMesh> dimScalarField;
autoPtr<dimScalarField> pSh;
// pick up the first fitting class
castAndCall(pSh,dimScalarField,basicThermoCloud,thermoCloud,hsTrans());
castAndCall(pSh,dimScalarField,basicReactingCloud,reactingCloud,hsTrans());
castAndCall(pSh,dimScalarField,basicReactingMultiphaseCloud,reactingMultiphaseCloud,hsTrans());
noCloudFound(pSh);
const dimScalarField &Sh=pSh();
autoPtr<volScalarField> pSource(
new volScalarField(
IOobject(
cloudName()+"EnthalpySource",
mesh().time().timeName(),
mesh(),
IOobject::NO_READ,
IOobject::NO_WRITE
),
mesh(),
Sh.dimensions()/(dimTime*dimVolume),
"zeroGradient"
)
);
pSource->internalField()=Sh.field()/(mesh().V()*mesh().time().deltaT().value());
result().setObjectResult(pSource);
}
void lcsMassSourcePluginFunction::doEvaluation()
{
typedef DimensionedField<scalar,volMesh> dimScalarField;
autoPtr<dimScalarField> pSrho;
// pick up the first fitting class
castAndCall(pSrho,dimScalarField,basicReactingCloud,reactingCloud,Srho());
castAndCall(pSrho,dimScalarField,basicReactingMultiphaseCloud,reactingMultiphaseCloud,Srho());
noCloudFound(pSrho);
const dimScalarField &Srho=pSrho();
autoPtr<volScalarField> pSource(
new volScalarField(
IOobject(
cloudName()+"MassSource",
mesh().time().timeName(),
mesh(),
IOobject::NO_READ,
IOobject::NO_WRITE
),
mesh(),
Srho.dimensions(),
"zeroGradient"
)
);
pSource->internalField()=Srho.field();
result().setObjectResult(pSource);
}
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);
}