diff --git a/src/thermophysicalModels/solidThermo/solidThermo/heSolidThermo.C b/src/thermophysicalModels/solidThermo/solidThermo/heSolidThermo.C
index b98b178c55b2ed2c4c907091b8df045a7268fd2a..e6c354e7dc4fb6ae059d90acba00e23e5a6cc9c1 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(this->mu_.dimensions(), Zero);
+    this->psi_ == dimensionedScalar(this->psi_.dimensions(), Zero);
 }
 
 
@@ -176,6 +178,8 @@ heSolidThermo
     heThermo<BasicSolidThermo, MixtureType>(mesh, dict, phaseName)
 {
     calculate();
+    this->mu_ == dimensionedScalar(this->mu_.dimensions(), Zero);
+    this->psi_ == dimensionedScalar(this->psi_.dimensions(), Zero);
 }