CFD Online Discussion Forums - radiation source term

Hi foamers,

I hope you are doing well.

I have a doubt about the calculation of radiation source term in openfoam. The source term is calculated in radiationModel.C using the function Sh(fluidThermo& thermo).

My confusion is about the term Ru(). It is calculated in radiation models, e.g. fvDOM, P1.

Code:

Foam::tmp<Foam::DimensionedField<Foam::scalar, Foam::volMesh> >

Foam::radiation::fvDOM(or P1)::Ru() const

{



    const DimensionedField<scalar, volMesh>& G =

        G_.dimensionedInternalField();

    const DimensionedField<scalar, volMesh> E =

        absorptionEmission_->ECont()().dimensionedInternalField();

    const DimensionedField<scalar, volMesh> a =

        a_.dimensionedInternalField();



    return  a*G - E;

}

I don't understand why it returns only dimensionedInternalField(). Why not use volScalarField and return a_*G_?
Any ideas? Thank you in advance.

Cheers,
Xu