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 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 lcsSpeciesSourcePluginFunction::doEvaluation()
{
wordList speciesList;
label speciesIndex=getIndex(speciesList);
if(speciesList.size()==0) {
listAvailableClouds(Info);
FatalErrorIn("lcsSpeciesSourcePluginFunction::doEvaluation()")
<< "No species list found (probably wrong cloud type)"
<< endl
<< "Supported cloud types (typeid): " \
<< supportedTypes() << endl \
<< exit(FatalError);
}
if(speciesIndex<0) {
listAvailableClouds(Info);
FatalErrorIn("lcsSpeciesSourcePluginFunction::doEvaluation()")
<< "Species " << speciesName_ << " not found in gas composition "
<< speciesList
<< endl
<< "Supported cloud types (typeid): " \
<< supportedTypes() << endl \
<< exit(FatalError);
}
autoPtr<dimScalarField> pSrho=internalEvaluate(speciesIndex);;
noCloudFound(pSrho);
const dimScalarField &Srho=pSrho();
autoPtr<volScalarField> pSource(
new volScalarField(
IOobject(
cloudName()+"SpeciesSource",
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);
}
