applications/solvers/heatTransfer: Added -postProcess option

See also commit cc455173ff

applications/solvers/heatTransfer/buoyantBoussinesqPimpleFoam/Make/options

@@ -1,4 +1,5 @@
 EXE_INC = \
+    -I. \
     -I../buoyantBoussinesqSimpleFoam \
     -I$(LIB_SRC)/TurbulenceModels/turbulenceModels/lnInclude \
     -I$(LIB_SRC)/TurbulenceModels/incompressible/lnInclude \

applications/solvers/heatTransfer/buoyantBoussinesqPimpleFoam/buoyantBoussinesqPimpleFoam.C

@@ -56,19 +56,18 @@ Description
 
 int main(int argc, char *argv[])
 {
+    #include "postProcess.H"
+
     #include "setRootCase.H"
     #include "createTime.H"
     #include "createMesh.H"
-
+    #include "createControl.H"
-    pimpleControl pimple(mesh);
-
     #include "createFields.H"
-    #include "createIncompressibleRadiationModel.H"
     #include "createFvOptions.H"
-    #include "initContinuityErrs.H"
     #include "createTimeControls.H"
     #include "CourantNo.H"
     #include "setInitialDeltaT.H"
+    #include "initContinuityErrs.H"
 
     turbulence->validate();
 

applications/solvers/heatTransfer/buoyantBoussinesqPimpleFoam/createFields.H

@@ -122,3 +122,4 @@ if (p_rgh.needReference())
 mesh.setFluxRequired(p_rgh.name());
 
 #include "createMRF.H"
+#include "createIncompressibleRadiationModel.H"

applications/solvers/heatTransfer/buoyantBoussinesqSimpleFoam/buoyantBoussinesqSimpleFoam.C

@@ -56,14 +56,13 @@ Description
 
 int main(int argc, char *argv[])
 {
+    #include "postProcess.H"
+
     #include "setRootCase.H"
     #include "createTime.H"
     #include "createMesh.H"
-
+    #include "createControl.H"
-    simpleControl simple(mesh);
-
     #include "createFields.H"
-    #include "createIncompressibleRadiationModel.H"
     #include "createFvOptions.H"
     #include "initContinuityErrs.H"
 

applications/solvers/heatTransfer/buoyantBoussinesqSimpleFoam/createFields.H

@@ -122,3 +122,4 @@ if (p_rgh.needReference())
 mesh.setFluxRequired(p_rgh.name());
 
 #include "createMRF.H"
+#include "createIncompressibleRadiationModel.H"

applications/solvers/heatTransfer/buoyantPimpleFoam/buoyantPimpleFoam.C

@@ -44,15 +44,15 @@ Description
 
 int main(int argc, char *argv[])
 {
+    #include "postProcess.H"
+
     #include "setRootCase.H"
     #include "createTime.H"
     #include "createMesh.H"
-
+    #include "createControl.H"
-    pimpleControl pimple(mesh);
-
     #include "createFields.H"
+    #include "createFieldRefs.H"
     #include "createFvOptions.H"
-    #include "createRadiationModel.H"
     #include "initContinuityErrs.H"
     #include "createTimeControls.H"
     #include "compressibleCourantNo.H"

applications/solvers/heatTransfer/buoyantPimpleFoam/createFieldRefs.H

@@ -0,0 +1 @@
+const volScalarField& psi = thermo.psi();

applications/solvers/heatTransfer/buoyantPimpleFoam/createFields.H

@@ -18,8 +18,6 @@ volScalarField rho
 );
 
 volScalarField& p = thermo.p();
-const volScalarField& psi = thermo.psi();
-
 
 Info<< "Reading field U\n" << endl;
 volVectorField U
@@ -92,3 +90,4 @@ Info<< "Creating field kinetic energy K\n" << endl;
 volScalarField K("K", 0.5*magSqr(U));
 
 #include "createMRF.H"
+#include "createRadiationModel.H"

applications/solvers/heatTransfer/buoyantSimpleFoam/buoyantSimpleFoam.C

@@ -41,15 +41,15 @@ Description
 
 int main(int argc, char *argv[])
 {
+    #include "postProcess.H"
+
     #include "setRootCase.H"
     #include "createTime.H"
     #include "createMesh.H"
-
+    #include "createControl.H"
-    simpleControl simple(mesh);
-
     #include "createFields.H"
+    #include "createFieldRefs.H"
     #include "createFvOptions.H"
-    #include "createRadiationModel.H"
     #include "initContinuityErrs.H"
 
     turbulence->validate();

applications/solvers/heatTransfer/buoyantSimpleFoam/createFieldRefs.H

@@ -0,0 +1 @@
+const volScalarField& psi = thermo.psi();

applications/solvers/heatTransfer/buoyantSimpleFoam/createFields.H

@@ -18,7 +18,6 @@ volScalarField rho
 );
 
 volScalarField& p = thermo.p();
-const volScalarField& psi = thermo.psi();
 
 Info<< "Reading field U\n" << endl;
 volVectorField U
@@ -88,3 +87,4 @@ dimensionedScalar initialMass = fvc::domainIntegrate(rho);
 dimensionedScalar totalVolume = sum(mesh.V());
 
 #include "createMRF.H"
+#include "createRadiationModel.H"

applications/solvers/heatTransfer/chtMultiRegionFoam/Make/options

@@ -1,5 +1,6 @@
 
 EXE_INC = \
+    -I. \
     -I./fluid \
     -I./solid \
     -I./porousFluid \

applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionFoam.C

@@ -50,22 +50,17 @@ Description
 
 int main(int argc, char *argv[])
 {
+    #define NO_CONTROL
+    #define CREATE_MESH createMeshesPostProcess.H
+    #include "postProcess.H"
+
     #include "setRootCase.H"
     #include "createTime.H"
-
+    #include "createMeshes.H"
-    regionProperties rp(runTime);
+    #include "createFields.H"
-
-    #include "createFluidMeshes.H"
-    #include "createSolidMeshes.H"
-
-    #include "createFluidFields.H"
-    #include "createSolidFields.H"
-
     #include "initContinuityErrs.H"
     #include "createTimeControls.H"
     #include "readSolidTimeControls.H"
-
-
     #include "compressibleMultiRegionCourantNo.H"
     #include "solidRegionDiffusionNo.H"
     #include "setInitialMultiRegionDeltaT.H"

applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/Make/options

@@ -1,4 +1,6 @@
 EXE_INC = \
+    -I. \
+    -I.. \
     -Ifluid \
     -Isolid \
     -I../solid \

applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/chtMultiRegionSimpleFoam.C

@@ -43,17 +43,14 @@ Description
 
 int main(int argc, char *argv[])
 {
+    #define NO_CONTROL
+    #define CREATE_MESH createMeshesPostProcess.H
+    #include "postProcess.H"
+
     #include "setRootCase.H"
     #include "createTime.H"
-
+    #include "createMeshes.H"
-    regionProperties rp(runTime);
+    #include "createFields.H"
-
-    #include "createFluidMeshes.H"
-    #include "createSolidMeshes.H"
-
-    #include "createFluidFields.H"
-    #include "createSolidFields.H"
-
     #include "initContinuityErrs.H"
 
 

applications/solvers/heatTransfer/chtMultiRegionFoam/createFields.H

@@ -0,0 +1,2 @@
+#include "createFluidFields.H"
+#include "createSolidFields.H"

applications/solvers/heatTransfer/chtMultiRegionFoam/createMeshes.H

@@ -0,0 +1,4 @@
+regionProperties rp(runTime);
+
+#include "createFluidMeshes.H"
+#include "createSolidMeshes.H"

applications/solvers/heatTransfer/chtMultiRegionFoam/createMeshesPostProcess.H

@@ -0,0 +1,2 @@
+#include "createMeshes.H"
+fvMesh& mesh = fluidRegions[0];

applications/solvers/heatTransfer/thermoFoam/createFields.H

@@ -1,54 +1,56 @@
-    Info<< "Reading thermophysical properties\n" << endl;
+Info<< "Reading thermophysical properties\n" << endl;
 
-    autoPtr<rhoThermo> pThermo(rhoThermo::New(mesh));
+autoPtr<rhoThermo> pThermo(rhoThermo::New(mesh));
-    rhoThermo& thermo = pThermo();
+rhoThermo& thermo = pThermo();
-    thermo.validate(args.executable(), "h", "e");
+thermo.validate(args.executable(), "h", "e");
 
-    volScalarField rho
+volScalarField rho
+(
+    IOobject
     (
-        IOobject
+        "rho",
-        (
+        runTime.timeName(),
-            "rho",
-            runTime.timeName(),
-            mesh,
-            IOobject::NO_READ,
-            IOobject::NO_WRITE
-        ),
-        thermo.rho()
-    );
-
-    volScalarField& p = thermo.p();
-
-    Info<< "Reading field U\n" << endl;
-    volVectorField U
-    (
-        IOobject
-        (
-            "U",
-            runTime.timeName(),
-            mesh,
-            IOobject::MUST_READ,
-            IOobject::AUTO_WRITE
-        ),
-        mesh
-    );
-
-    #include "compressibleCreatePhi.H"
-
-    #include "setAlphaEff.H"
-
-    Info<< "Creating field dpdt\n" << endl;
-    volScalarField dpdt
-    (
-        IOobject
-        (
-            "dpdt",
-            runTime.timeName(),
-            mesh
-        ),
         mesh,
-        dimensionedScalar("dpdt", p.dimensions()/dimTime, 0)
+        IOobject::NO_READ,
-    );
+        IOobject::NO_WRITE
+    ),
+    thermo.rho()
+);
 
-    Info<< "Creating field kinetic energy K\n" << endl;
+volScalarField& p = thermo.p();
-    volScalarField K("K", 0.5*magSqr(U));
+
+Info<< "Reading field U\n" << endl;
+volVectorField U
+(
+    IOobject
+    (
+        "U",
+        runTime.timeName(),
+        mesh,
+        IOobject::MUST_READ,
+        IOobject::AUTO_WRITE
+    ),
+    mesh
+);
+
+#include "compressibleCreatePhi.H"
+
+#include "setAlphaEff.H"
+
+Info<< "Creating field dpdt\n" << endl;
+volScalarField dpdt
+(
+    IOobject
+    (
+        "dpdt",
+        runTime.timeName(),
+        mesh
+    ),
+    mesh,
+    dimensionedScalar("dpdt", p.dimensions()/dimTime, 0)
+);
+
+Info<< "Creating field kinetic energy K\n" << endl;
+volScalarField K("K", 0.5*magSqr(U));
+
+#include "createRadiationModel.H"

applications/solvers/heatTransfer/thermoFoam/setAlphaEff.H

@@ -1,49 +1,47 @@
-    Info<< "Creating turbulence model\n" << endl;
+Info<< "Creating turbulence model\n" << endl;
-    tmp<volScalarField> talphaEff;
+tmp<volScalarField> talphaEff;
 
-    IOobject turbulencePropertiesHeader
+IOobject turbulencePropertiesHeader
+(
+    "turbulenceProperties",
+    runTime.constant(),
+    mesh,
+    IOobject::MUST_READ_IF_MODIFIED,
+    IOobject::NO_WRITE,
+    false
+);
+
+if (turbulencePropertiesHeader.headerOk())
+{
+    autoPtr<compressible::turbulenceModel> turbulence
     (
-        "turbulenceProperties",
+        compressible::turbulenceModel::New
-        runTime.constant(),
+        (
-        mesh,
+            rho,
-        IOobject::MUST_READ_IF_MODIFIED,
+            U,
-        IOobject::NO_WRITE,
+            phi,
-        false
+            thermo
+        )
     );
 
-    if (turbulencePropertiesHeader.headerOk())
+    talphaEff = turbulence->alphaEff();
-    {
+}
-        autoPtr<compressible::turbulenceModel> turbulence
+else
+{
+    talphaEff = tmp<volScalarField>
+    (
+        new volScalarField
         (
-            compressible::turbulenceModel::New
+            IOobject
             (
-                rho,
+                "alphaEff",
-                U,
+                runTime.timeName(),
-                phi,
-                thermo
-            )
-        );
-
-        talphaEff = turbulence->alphaEff();
-    }
-    else
-    {
-        talphaEff = tmp<volScalarField>
-        (
-            new volScalarField
-            (
-                IOobject
-                (
-                    "alphaEff",
-                    runTime.timeName(),
-                    mesh,
-                    IOobject::NO_READ,
-                    IOobject::NO_WRITE
-                ),
                 mesh,
-                dimensionedScalar("0", dimMass/dimLength/dimTime, 0.0)
+                IOobject::NO_READ,
-            )
+                IOobject::NO_WRITE
-        );
+            ),
-    }
+            mesh,
-
+            dimensionedScalar("0", dimMass/dimLength/dimTime, 0.0)
-    const volScalarField& alphaEff = talphaEff();
+        )
+    );
+}

applications/solvers/heatTransfer/thermoFoam/thermoFoam.C

@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2013-2015 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2013-2016 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -43,13 +43,16 @@ Description
 
 int main(int argc, char *argv[])
 {
-    #include "setRootCase.H"
+    #define NO_CONTROL
+    #include "postProcess.H"
 
+    #include "setRootCase.H"
     #include "createTime.H"
     #include "createMesh.H"
     #include "createFields.H"
     #include "createFvOptions.H"
-    #include "createRadiationModel.H"
+
+    const volScalarField& alphaEff = talphaEff();
 
     // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //