diff --git a/src/thermoTools/derivedFvPatchFields/externalWallHeatFluxTemperature/externalWallHeatFluxTemperatureFvPatchScalarField.C b/src/thermoTools/derivedFvPatchFields/externalWallHeatFluxTemperature/externalWallHeatFluxTemperatureFvPatchScalarField.C
index 9b186b78d0408e180344b6f67ec6dfda2ac4f83e..67f622f2ccf417b65799a3a8018513c906c72839 100644
--- a/src/thermoTools/derivedFvPatchFields/externalWallHeatFluxTemperature/externalWallHeatFluxTemperatureFvPatchScalarField.C
+++ b/src/thermoTools/derivedFvPatchFields/externalWallHeatFluxTemperature/externalWallHeatFluxTemperatureFvPatchScalarField.C
@@ -331,7 +331,7 @@ void Foam::externalWallHeatFluxTemperatureFvPatchScalarField::updateCoeffs()
                 Q_->value(this->db().time().timeOutputValue());
 
             refGrad() = (heatPower/gSum(patch().magSf()) + qr)/kappa(Tp);
-            refValue() = 0;
+            refValue() = 293.15;  // prevents FPE, no impact on condition
             valueFraction() = 0;
 
             break;
@@ -342,7 +342,7 @@ void Foam::externalWallHeatFluxTemperatureFvPatchScalarField::updateCoeffs()
                 q_->value(this->db().time().timeOutputValue());
 
             refGrad() = (heatFlux + qr)/kappa(Tp);
-            refValue() = 0;
+            refValue() = 293.15;  // prevents FPE, no impact on condition
             valueFraction() = 0;
 
             break;