ENH: Updated pressure-work term creation for moving meshes

src/dynamicFvMesh/dynamicFvMesh/dynamicFvMesh.H

@@ -40,6 +40,8 @@ SourceFiles
 #include "autoPtr.H"
 #include "runTimeSelectionTables.H"
 
+#define MOVING_MESH
+
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
 namespace Foam

src/engine/engineMesh/engineMesh/engineMesh.H

@@ -40,6 +40,8 @@ SourceFiles
 #include "autoPtr.H"
 #include "runTimeSelectionTables.H"
 
+#define MOVING_MESH
+
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
 namespace Foam

src/finiteVolume/cfdTools/compressible/createDpdt.H

@@ -1,6 +1,25 @@
-Info<< "Creating field dpdt\n" << endl;
+#if defined(MOVING_MESH)
-volScalarField dpdt
+    Info<< "Creating field dpdt for moving meshes\n" << endl;
-(
+
+    // Note: set to read if present and auto write to simplify dpdt correction
+    // by meshPhi
+    volScalarField dpdt
+    (
+        IOobject
+        (
+            "dpdt",
+            runTime.timeName(),
+            mesh,
+            IOobject::READ_IF_PRESENT,
+            IOobject::AUTO_WRITE
+        ),
+        fvc::ddt(p)
+    );
+#else
+    Info<< "Creating field dpdt\n" << endl;
+
+    volScalarField dpdt
+    (
         IOobject
         (
             "dpdt",
@@ -10,9 +29,11 @@ volScalarField dpdt
             IOobject::NO_WRITE
         ),
         fvc::ddt(p)
-);
+    );
+#endif
 
 if (!thermo.dpdt())
 {
     dpdt == dimensionedScalar("0", dpdt.dimensions(), 0);
+    dpdt.writeOpt() = IOobject::NO_WRITE;
 }