compressibleInterFoam: Corrected handling for rAU for restart of moving-mesh cases

applications/solvers/multiphase/compressibleInterFoam/compressibleInterFoam.C

@@ -123,14 +123,7 @@ int main(int argc, char *argv[])
 
                     if (correctPhi)
                     {
-                        // Calculate absolute flux
-                        // from the mapped surface velocity
-                        phi = mesh.Sf() & Uf();
-
                         #include "correctPhi.H"
-
-                        // Make the fluxes relative to the mesh motion
-                        fvc::makeRelative(phi, U);
                     }
 
                     mixture.correct();

applications/solvers/multiphase/compressibleInterFoam/correctPhi.H

@@ -1,3 +1,7 @@
+// Calculate absolute flux
+// from the mapped surface velocity
+phi = mesh.Sf() & Uf();
+
 CorrectPhi
 (
     U,
@@ -10,3 +14,6 @@ CorrectPhi
 );
 
 #include "continuityErrs.H"
+
+// Make the fluxes relative to the mesh motion
+fvc::makeRelative(phi, U);

applications/solvers/multiphase/compressibleInterFoam/createFieldRefs.H

@@ -11,3 +11,20 @@ const volScalarField& psi1 = mixture.thermo1().psi();
 const volScalarField& psi2 = mixture.thermo2().psi();
 
 tmp<volScalarField> rAU;
+
+if (correctPhi)
+{
+    rAU = new volScalarField
+    (
+        IOobject
+        (
+            "rAU",
+            runTime.timeName(),
+            mesh,
+            IOobject::READ_IF_PRESENT,
+            IOobject::AUTO_WRITE
+        ),
+        mesh,
+        dimensionedScalar(dimTime/dimDensity, 1)
+    );
+}

applications/solvers/multiphase/interFoam/correctPhi.H

@@ -1,3 +1,7 @@
+// Calculate absolute flux
+// from the mapped surface velocity
+phi = mesh.Sf() & Uf();
+
 if (divU.valid())
 {
     CorrectPhi
@@ -26,3 +30,6 @@ else
 }
 
 #include "continuityErrs.H"
+
+// Make the flux relative to the mesh motion
+fvc::makeRelative(phi, U);

applications/solvers/multiphase/interFoam/interFoam.C

@@ -135,14 +135,7 @@ int main(int argc, char *argv[])
 
                     if (correctPhi)
                     {
-                        // Calculate absolute flux
-                        // from the mapped surface velocity
-                        phi = mesh.Sf() & Uf();
-
                         #include "correctPhi.H"
-
-                        // Make the flux relative to the mesh motion
-                        fvc::makeRelative(phi, U);
                     }
 
                     mixture.correct();