diff --git a/src/thermophysicalModels/radiation/derivedFvPatchFields/greyDiffusiveRadiation/greyDiffusiveRadiationMixedFvPatchScalarField.C b/src/thermophysicalModels/radiation/derivedFvPatchFields/greyDiffusiveRadiation/greyDiffusiveRadiationMixedFvPatchScalarField.C index e7daa58aaf..1f9a6ac63c 100644 --- a/src/thermophysicalModels/radiation/derivedFvPatchFields/greyDiffusiveRadiation/greyDiffusiveRadiationMixedFvPatchScalarField.C +++ b/src/thermophysicalModels/radiation/derivedFvPatchFields/greyDiffusiveRadiation/greyDiffusiveRadiationMixedFvPatchScalarField.C @@ -151,7 +151,10 @@ updateCoeffs() const scalarField& Tp = patch().lookupPatchField(TName_); - const fvDOM& dom = db().lookupObject("radiationProperties"); + const radiationModel& radiation = + db().lookupObject("radiationProperties"); + + const fvDOM& dom(refCast(radiation)); const label patchI = patch().index(); diff --git a/src/thermophysicalModels/radiation/derivedFvPatchFields/wideBandDiffusiveRadiation/wideBandDiffusiveRadiationMixedFvPatchScalarField.C b/src/thermophysicalModels/radiation/derivedFvPatchFields/wideBandDiffusiveRadiation/wideBandDiffusiveRadiationMixedFvPatchScalarField.C index eb65337aca..eec4d22021 100644 --- a/src/thermophysicalModels/radiation/derivedFvPatchFields/wideBandDiffusiveRadiation/wideBandDiffusiveRadiationMixedFvPatchScalarField.C +++ b/src/thermophysicalModels/radiation/derivedFvPatchFields/wideBandDiffusiveRadiation/wideBandDiffusiveRadiationMixedFvPatchScalarField.C @@ -147,7 +147,10 @@ updateCoeffs() return; } - const fvDOM& dom = db().lookupObject("radiationProperties"); + const radiationModel& radiation = + db().lookupObject("radiationProperties"); + + const fvDOM& dom(refCast(radiation)); const label patchI = patch().index(); diff --git a/src/thermophysicalModels/radiation/radiationModel/fvDOM/radiativeIntensityRay/radiativeIntensityRay.C b/src/thermophysicalModels/radiation/radiationModel/fvDOM/radiativeIntensityRay/radiativeIntensityRay.C index 96090a3ef2..74ec95827c 100644 --- a/src/thermophysicalModels/radiation/radiationModel/fvDOM/radiativeIntensityRay/radiativeIntensityRay.C +++ b/src/thermophysicalModels/radiation/radiationModel/fvDOM/radiativeIntensityRay/radiativeIntensityRay.C @@ -182,7 +182,7 @@ Foam::scalar Foam::radiation::radiativeIntensityRay::correct() fvScalarMatrix IiEq ( - fvm::div(Ji, ILambda_[lambdaI], " div(Ji,Ii_h)") + fvm::div(Ji, ILambda_[lambdaI], "div(Ji,Ii_h)") + fvm::Sp(k*omega_, ILambda_[lambdaI]) == 1.0/Foam::mathematicalConstant::pi