diff --git a/src/turbulenceModels/compressible/turbulenceModel/derivedFvPatchFields/externalWallHeatFluxTemperature/externalWallHeatFluxTemperatureFvPatchScalarField.C b/src/turbulenceModels/compressible/turbulenceModel/derivedFvPatchFields/externalWallHeatFluxTemperature/externalWallHeatFluxTemperatureFvPatchScalarField.C index 8446d56e1bd756bb5ee1881bc253200b8554e0ad..47bc69aeb6d8dea60c07cf3156a799b6c925e49f 100644 --- a/src/turbulenceModels/compressible/turbulenceModel/derivedFvPatchFields/externalWallHeatFluxTemperature/externalWallHeatFluxTemperatureFvPatchScalarField.C +++ b/src/turbulenceModels/compressible/turbulenceModel/derivedFvPatchFields/externalWallHeatFluxTemperature/externalWallHeatFluxTemperatureFvPatchScalarField.C @@ -229,7 +229,8 @@ void Foam::externalWallHeatFluxTemperatureFvPatchScalarField::updateCoeffs() scalarField q(size(), 0.0); const scalarField Tc(patchInternalField()); - const scalarField KWall(kappa(*this)); + const scalarField Tp(*this); + const scalarField KWall(kappa(Tp)); const scalarField KDelta(KWall*patch().deltaCoeffs()); switch (mode_) @@ -241,7 +242,7 @@ void Foam::externalWallHeatFluxTemperatureFvPatchScalarField::updateCoeffs() } case fixedHeatTransferCoeff: { - q = (Ta_ - Tc)*h_; + q = (Ta_ - Tp)*h_; break; } default: