ENH: Added include file to calculate the mechanical energy

applications/solvers/compressible/rhoPimpleFoam/createFields.H

@@ -69,8 +69,7 @@ volScalarField dpdt
     dimensionedScalar("dpdt", p.dimensions()/dimTime, 0)
 );
 
-Info<< "Creating field kinetic energy K\n" << endl;
-volScalarField K("K", 0.5*magSqr(U));
+#include "createK.H"
 
 #include "createMRF.H"
 

src/finiteVolume/cfdTools/general/include/createK.H

+Info<< "Creating field kinetic energy K\n" << endl;
+volScalarField K("K", 0.5*magSqr(U));
+
+if (U.nOldTimes())
+{
+    volVectorField* Uold = &U.oldTime();
+    volScalarField* Kold = &K.oldTime();
+    *Kold == 0.5*magSqr(*Uold);
+
+    while (Uold->nOldTimes())
+    {
+        Uold = &Uold->oldTime();
+        Kold = &Kold->oldTime();
+        *Kold == 0.5*magSqr(*Uold);
+    }
+}