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; }