diff --git a/src/thermophysicalModels/radiationModels/submodels/absorptionEmissionModel/greyMeanAbsorptionEmission/greyMeanAbsorptionEmission.C b/src/thermophysicalModels/radiationModels/submodels/absorptionEmissionModel/greyMeanAbsorptionEmission/greyMeanAbsorptionEmission.C index 3edfd4ce38..08dd1e3764 100644 --- a/src/thermophysicalModels/radiationModels/submodels/absorptionEmissionModel/greyMeanAbsorptionEmission/greyMeanAbsorptionEmission.C +++ b/src/thermophysicalModels/radiationModels/submodels/absorptionEmissionModel/greyMeanAbsorptionEmission/greyMeanAbsorptionEmission.C @@ -2,7 +2,7 @@ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | - \\ / A nd | Copyright (C) 2011-2012 OpenFOAM Foundation + \\ / A nd | Copyright (C) 2011-2013 OpenFOAM Foundation \\/ M anipulation | ------------------------------------------------------------------------------- License @@ -312,7 +312,30 @@ Foam::radiation::greyMeanAbsorptionEmission::ECont(const label bandI) const { const volScalarField& dQ = mesh_.lookupObject("dQ"); - E().internalField() = EhrrCoeff_*dQ; + + if (dQ.dimensions() == dimEnergy/dimTime) + { + E().internalField() = EhrrCoeff_*dQ/mesh_.V(); + } + else if (dQ.dimensions() == dimEnergy/dimTime/dimVolume) + { + E().internalField() = EhrrCoeff_*dQ; + } + else + { + if (debug) + { + WarningIn + ( + "tmp" + "radiation::greyMeanAbsorptionEmission::ECont" + "(" + "const label" + ") const" + ) + << "Incompatible dimensions for dQ field" << endl; + } + } } return E;