ENH: Adding option mesh flux correction to solidFoam

applications/solvers/heatTransfer/solidFoam/createFields.H

@@ -90,5 +90,24 @@ else
     );
 }
 
+// Consider mesh flux to correct for mesh deformation
+bool meshFluxCorr(false);
+if (mesh.solutionDict().found("SIMPLE"))
+{
+    meshFluxCorr =
+        mesh.solutionDict().subDict("SIMPLE").getOrDefault<bool>
+        (
+            "meshFluxCorrection", false
+        );
+}
+else if (mesh.solutionDict().found("PIMPLE"))
+{
+    meshFluxCorr =
+        mesh.solutionDict().subDict("PIMPLE").getOrDefault<bool>
+        (
+            "meshFluxCorrection", false
+        );
+}
+
 #include "createRadiationModel.H"
 #include "createFvOptions.H"

applications/solvers/heatTransfer/solidFoam/hEqn.H

@@ -11,14 +11,14 @@
         fvOptions(rho, h)
     );
 
-    if (mesh.changing())
+    if (meshFluxCorr)
     {
         surfaceScalarField phihMesh
         (
             fvc::interpolate(betav*rho*h)*mesh.phi()
         );
 
-        hEqn -= fvc::div(phihMesh);
+       hEqn -= fvc::div(phihMesh);
     }
 
     hEqn.relax();

tutorials/heatTransfer/solidFoam/movingCone/system/fvSolution

@@ -46,6 +46,7 @@ PIMPLE
     nOuterCorrectors    2;
     nCorrectors         1;
     nNonOrthogonalCorrectors 0;
+    meshFluxCorrection  true;
 }
 
 relaxationFactors