36 Foam::radiation::radiativeIntensityRay::radiativeIntensityRay
42 const scalar deltaPhi,
43 const scalar deltaTheta,
52 absorptionEmission_(absorptionEmission),
53 blackBody_(blackBody),
59 mesh_.time().timeName(),
72 mesh_.time().timeName(),
85 mesh_.time().timeName(),
98 mesh_.time().timeName(),
119 omega_ = 2.0*sinTheta*
Foam::sin(deltaTheta/2.0)*deltaPhi;
120 d_ =
vector(sinTheta*sinPhi, sinTheta*cosPhi, cosTheta);
141 intensityPrefix +
"_" +
name(rayId) +
"_" +
name(lambdaI),
142 mesh_.time().timeName(),
149 if (IHeader.headerOk())
160 if (!IDefaultPtr.
valid())
169 mesh_.time().timeName(),
181 noReadHeader.
readOpt() = IOobject::NO_READ;
204 Qr_.boundaryField() = 0.0;
216 fvm::div(Ji, ILambda_[lambdaI],
"div(Ji,Ii_h)")
217 +
fvm::Sp(k*omega_, ILambda_[lambdaI])
221 k*blackBody_.bLambda(lambdaI)
222 + absorptionEmission_.ECont(lambdaI)/4
234 maxResidual =
max(eqnResidual, maxResidual);
247 I_ += absorptionEmission_.addIntensity(lambdaI, ILambda_[lambdaI]);