diff --git a/src/thermophysicalModels/solidThermo/solidThermo/heSolidThermo.C b/src/thermophysicalModels/solidThermo/solidThermo/heSolidThermo.C
index f99ac1ba9b2da73998bde98857ab18ed84257093..0cdfeec3d678a442679d2cbc5adf1bbf016fdf49 100644
--- a/src/thermophysicalModels/solidThermo/solidThermo/heSolidThermo.C
+++ b/src/thermophysicalModels/solidThermo/solidThermo/heSolidThermo.C
@@ -161,6 +161,8 @@ heSolidThermo
     heThermo<BasicSolidThermo, MixtureType>(mesh, phaseName)
 {
     calculate();
+    this->mu_ == dimensionedScalar("zero", this->mu_.dimensions(), 0.0);
+    this->psi_ == dimensionedScalar("zero", this->psi_.dimensions(), 0.0);
 }
 
 
@@ -176,6 +178,8 @@ heSolidThermo
     heThermo<BasicSolidThermo, MixtureType>(mesh, dict, phaseName)
 {
     calculate();
+    this->mu_ == dimensionedScalar("zero", this->mu_.dimensions(), 0.0);
+    this->psi_ == dimensionedScalar("zero", this->psi_.dimensions(), 0.0);
 }