diff --git a/src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/turbulentHeatFluxTemperature/turbulentHeatFluxTemperatureFvPatchScalarField.C b/src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/turbulentHeatFluxTemperature/turbulentHeatFluxTemperatureFvPatchScalarField.C index 15632c9026..40af0b6b9f 100644 --- a/src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/turbulentHeatFluxTemperature/turbulentHeatFluxTemperatureFvPatchScalarField.C +++ b/src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/turbulentHeatFluxTemperature/turbulentHeatFluxTemperatureFvPatchScalarField.C @@ -27,6 +27,7 @@ License #include "addToRunTimeSelectionTable.H" #include "fvPatchFieldMapper.H" #include "volFields.H" +#include "RASModel.H" // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // @@ -173,20 +174,21 @@ void turbulentHeatFluxTemperatureFvPatchScalarField::updateCoeffs() const scalarField& alphaEffp = patch().lookupPatchField(alphaEffName_); - const scalarField& Cpp = - patch().lookupPatchField(CpName_); + // retrieve (constant) specific heat capacity from transport dictionary + const RASModel& rasModel = db().lookupObject("RASProperties"); + const scalar Cp0(readScalar(rasModel.transport().lookup("Cp0"))); switch (heatSource_) { case hsPower: { const scalar Ap = gSum(patch().magSf()); - gradient() = q_/(Ap*Cpp*alphaEffp); + gradient() = q_/(Ap*Cp0*alphaEffp); break; } case hsFlux: { - gradient() = q_/(Cpp*alphaEffp); + gradient() = q_/(Cp0*alphaEffp); break; } default: @@ -215,7 +217,6 @@ void turbulentHeatFluxTemperatureFvPatchScalarField::write(Ostream& os) const << token::END_STATEMENT << nl; q_.writeEntry("q", os); os.writeKeyword("alphaEff") << alphaEffName_ << token::END_STATEMENT << nl; - os.writeKeyword("Cp") << CpName_ << token::END_STATEMENT << nl; writeEntry("value", os); }