interPhaseChangeFoam, multiphaseInterFoam, potentialFreeSurfaceFoam: Added -postProcess option

See also commit cc455173ff

applications/solvers/multiphase/interPhaseChangeFoam/createControls.H

@@ -0,0 +1 @@
+pimpleControl pimple(mesh);

applications/solvers/multiphase/interPhaseChangeFoam/createFields.H

@@ -38,7 +38,6 @@ volScalarField& alpha2(mixture->alpha2());
 
 const dimensionedScalar& rho1 = mixture->rho1();
 const dimensionedScalar& rho2 = mixture->rho2();
-const dimensionedScalar& pSat = mixture->pSat();
 
 
 // Need to store rho for ddt(rho, U)

applications/solvers/multiphase/interPhaseChangeFoam/interPhaseChangeDyMFoam/interPhaseChangeDyMFoam.C

@@ -56,13 +56,14 @@ Description
 
 int main(int argc, char *argv[])
 {
+    #include "postProcess.H"
+
     #include "setRootCase.H"
     #include "createTime.H"
     #include "createDynamicFvMesh.H"
-
+    #include "createControls.H"
-    pimpleControl pimple(mesh);
+    #include "createTimeControls.H"
-
+    #include "../interFoam/interDyMFoam/createDyMControls.H"
-    #include "../interFoam/interDyMFoam/createControls.H"
     #include "initContinuityErrs.H"
     #include "createFields.H"
     #include "createFvOptions.H"

applications/solvers/multiphase/interPhaseChangeFoam/interPhaseChangeDyMFoam/pEqn.H

@@ -38,7 +38,8 @@
         fvScalarMatrix p_rghEqn
         (
             fvc::div(phiHbyA) - fvm::laplacian(rAUf, p_rgh)
-          - (vDotvP - vDotcP)*(pSat - rho*gh) + fvm::Sp(vDotvP - vDotcP, p_rgh)
+          - (vDotvP - vDotcP)*(mixture->pSat() - rho*gh)
+          + fvm::Sp(vDotvP - vDotcP, p_rgh)
         );
 
         p_rghEqn.setReference(pRefCell, pRefValue);

applications/solvers/multiphase/interPhaseChangeFoam/interPhaseChangeFoam.C

@@ -53,12 +53,12 @@ Description
 
 int main(int argc, char *argv[])
 {
+    #include "postProcess.H"
+
     #include "setRootCase.H"
     #include "createTime.H"
     #include "createMesh.H"
-
+    #include "createControls.H"
-    pimpleControl pimple(mesh);
-
     #include "createFields.H"
     #include "createFvOptions.H"
     #include "createTimeControls.H"

applications/solvers/multiphase/interPhaseChangeFoam/pEqn.H

@@ -32,7 +32,8 @@
         fvScalarMatrix p_rghEqn
         (
             fvc::div(phiHbyA) - fvm::laplacian(rAUf, p_rgh)
-          - (vDotvP - vDotcP)*(pSat - rho*gh) + fvm::Sp(vDotvP - vDotcP, p_rgh)
+          - (vDotvP - vDotcP)*(mixture->pSat() - rho*gh)
+          + fvm::Sp(vDotvP - vDotcP, p_rgh)
         );
 
         p_rghEqn.setReference(pRefCell, pRefValue);

applications/solvers/multiphase/multiphaseInterFoam/Make/options

@@ -1,4 +1,5 @@
 EXE_INC = \
+    -I. \
     -I../interFoam \
     -ImultiphaseMixture/lnInclude \
     -I$(LIB_SRC)/transportModels \

applications/solvers/multiphase/multiphaseInterFoam/createFields.H

@@ -44,8 +44,6 @@ volScalarField rho
 );
 rho.oldTime();
 
-const surfaceScalarField& rhoPhi(mixture.rhoPhi());
-
 // Construct incompressible turbulence model
 autoPtr<incompressible::turbulenceModel> turbulence
 (

applications/solvers/multiphase/multiphaseInterFoam/multiphaseInterDyMFoam/multiphaseInterDyMFoam.C

@@ -44,14 +44,15 @@ Description
 
 int main(int argc, char *argv[])
 {
+    #include "postProcess.H"
+
     #include "setRootCase.H"
     #include "createTime.H"
     #include "createDynamicFvMesh.H"
     #include "initContinuityErrs.H"
-
-    pimpleControl pimple(mesh);
-
     #include "createControls.H"
+    #include "createTimeControls.H"
+    #include "createDyMControls.H"
     #include "createFields.H"
     #include "createMRF.H"
     #include "createFvOptions.H"
@@ -75,6 +76,8 @@ int main(int argc, char *argv[])
     #include "CourantNo.H"
     #include "setInitialDeltaT.H"
 
+    const surfaceScalarField& rhoPhi(mixture.rhoPhi());
+
     turbulence->validate();
 
     // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

applications/solvers/multiphase/multiphaseInterFoam/multiphaseInterFoam.C

@@ -43,12 +43,12 @@ Description
 
 int main(int argc, char *argv[])
 {
+    #include "postProcess.H"
+
     #include "setRootCase.H"
     #include "createTime.H"
     #include "createMesh.H"
-
+    #include "createControls.H"
-    pimpleControl pimple(mesh);
-
     #include "initContinuityErrs.H"
     #include "createFields.H"
     #include "createMRF.H"
@@ -60,6 +60,8 @@ int main(int argc, char *argv[])
 
     turbulence->validate();
 
+    const surfaceScalarField& rhoPhi(mixture.rhoPhi());
+
     // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
     Info<< "\nStarting time loop\n" << endl;

applications/solvers/multiphase/potentialFreeSurfaceFoam/createControls.H

@@ -0,0 +1 @@
+pimpleControl pimple(mesh);

applications/solvers/multiphase/potentialFreeSurfaceFoam/potentialFreeSurfaceDyMFoam/potentialFreeSurfaceDyMFoam.C

@@ -49,14 +49,15 @@ Description
 
 int main(int argc, char *argv[])
 {
+    #include "postProcess.H"
+
     #include "setRootCase.H"
     #include "createTime.H"
     #include "createDynamicFvMesh.H"
     #include "initContinuityErrs.H"
-
-    pimpleControl pimple(mesh);
-
     #include "createControls.H"
+    #include "createTimeControls.H"
+    #include "createDyMControls.H"
     #include "createFields.H"
     #include "createMRF.H"
     #include "createFvOptions.H"

applications/solvers/multiphase/potentialFreeSurfaceFoam/potentialFreeSurfaceFoam.C

@@ -44,13 +44,12 @@ Description
 
 int main(int argc, char *argv[])
 {
-    #include "setRootCase.H"
+    #include "postProcess.H"
 
+    #include "setRootCase.H"
     #include "createTime.H"
     #include "createMesh.H"
-
+    #include "createControls.H"
-    pimpleControl pimple(mesh);
-
     #include "createTimeControls.H"
     #include "createFields.H"
     #include "createMRF.H"