Foam::tmp<Foam::volScalarField>
Foam::radiation::cloudAbsorptionEmission::EDisp(const label bandI) const
{
tmp<volScalarField> tE
(
new volScalarField
(
IOobject
(
"E",
mesh_.time().timeName(),
mesh_,
IOobject::NO_READ,
IOobject::NO_WRITE,
false
),
mesh_,
dimensionedScalar("E", dimMass/dimLength/pow3(dimTime), 0.0)
)
);
forAll(cloudNames_, i)
{
const thermoCloud& tc
(
mesh_.objectRegistry::lookupObject<thermoCloud>(cloudNames_[i])
);
tE() += tc.Ep();
}
return tE;
}
Foam::tmp<Foam::volScalarField> Foam::MaxwellViscoelastic::E(scalar t) const
{
scalar tau = eta_.value()/k_.value();
tmp<volScalarField> tE
(
new volScalarField
(
IOobject
(
"E",
mesh().time().timeName(),
mesh(),
IOobject::NO_READ,
IOobject::NO_WRITE
),
mesh(),
k_*exp(-t/tau),
zeroGradientFvPatchScalarField::typeName
)
);
if (t < 0)
{
tE().internalField() = 0.0;
tE().correctBoundaryConditions();
}
return tE;
}
Foam::tmp<Foam::volScalarField>
Foam::radiation::constantAbsorptionEmission::ECont(const label bandI) const
{
tmp<volScalarField> tE
(
new volScalarField
(
IOobject
(
"E",
mesh_.time().timeName(),
mesh_,
IOobject::NO_READ,
IOobject::NO_WRITE,
false
),
mesh_,
E_
)
);
return tE;
}
