[makaveli_lcf]

I noticed an error in turbulentHeatFluxTemperatureFvPatchScalarField class for temperature gradient calculation:

00146 void turbulentHeatFluxTemperatureFvPatchScalarField::updateCoeffs()
00147 {
00148     if (updated())
00149     {
00150         return;
00151     }
00152 
00153     const scalarField& alphaEffp =
00154         patch().lookupPatchField<volScalarField, scalar>(alphaEffName_);
00155 
00156     const scalarField& Cpp =
00157         patch().lookupPatchField<volScalarField, scalar>(CpName_);
00158 
 00159     gradient() = q_/(Cpp*alphaEffp);
00160 
00161     fixedGradientFvPatchScalarField::updateCoeffs();
00162 }

From my point of view and based on Fourier's law heat flux q is calculated as:

q = lambda * grad (T) = alpha * rho * Cp * grad (T)

where
lambda - thermal conductivity
alpha - thermal diffusivity
Cp - cpecific heat capacity
rho - density

Thus temperature gradient should be

grad (T) = q / (alpha * rho * Cp)

Please correct me if I'm wrong!