3 alphat.correctBoundaryConditions();
dimensionedScalar rhoCpRef("rhoCpRef", dimDensity *dimEnergy/dimMass/dimTemperature, 1.0)
autoPtr< radiation::radiationModel > radiation(radiation::radiationModel::New(T))
compressible::turbulenceModel & turbulence
volScalarField alphaEff("alphaEff", turbulence->nu()/Pr+alphat)
fvScalarMatrix TEqn(fvm::ddt(T)+fvm::div(phi, T) - fvm::laplacian(alphaEff, T)==radiation->ST(rhoCpRef, T)+fvOptions(T))
dimensionedScalar TRef("TRef", dimTemperature, laminarTransport)
dimensionedScalar beta("beta", dimless/dimTemperature, laminarTransport)
dimensionedScalar Prt("Prt", dimless, laminarTransport)
dimensionedScalar Pr("Pr", dimless, laminarTransport)