From 22a5cee8a4a9ab7a93550c27305650ada7e1d3cb Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Daniel=20Jasi=C5=84ski?= <daniel.jasinski@gmail.com>
Date: Thu, 29 Oct 2015 10:23:47 +0100
Subject: [PATCH] buoyantBoussinesqSimpleFoam added radiation

---
 .../solvers/heatTransfer/buoyantBoussinesqSimpleFoam/Make/options    | 2 ++
 applications/solvers/heatTransfer/buoyantBoussinesqSimpleFoam/TEqn.H | 5 ++++-
 .../buoyantBoussinesqSimpleFoam/buoyantBoussinesqSimpleFoam.C        | 2 ++
 3 files changed, 8 insertions(+), 1 deletion(-)

diff --git a/applications/solvers/heatTransfer/buoyantBoussinesqSimpleFoam/Make/options b/applications/solvers/heatTransfer/buoyantBoussinesqSimpleFoam/Make/options
index 48d8383..5f9e323 100644
--- a/applications/solvers/heatTransfer/buoyantBoussinesqSimpleFoam/Make/options
+++ b/applications/solvers/heatTransfer/buoyantBoussinesqSimpleFoam/Make/options
@@ -3,6 +3,7 @@ EXE_INC = \
     -I$(LIB_SRC)/TurbulenceModels/incompressible/lnInclude \
     -I$(LIB_SRC)/transportModels \
     -I$(LIB_SRC)/transportModels/incompressible/singlePhaseTransportModel \
+    -I$(LIB_SRC)/thermophysicalModels/radiation/lnInclude \
     -I$(LIB_SRC)/finiteVolume/lnInclude \
     -I$(LIB_SRC)/sampling/lnInclude \
     -I$(LIB_SRC)/meshTools/lnInclude \
@@ -12,6 +13,7 @@ EXE_LIBS = \
     -lturbulenceModels \
     -lincompressibleTurbulenceModels \
     -lincompressibleTransportModels \
+    -lradiationModels \
     -lfiniteVolume \
     -lsampling \
     -lmeshTools \
diff --git a/applications/solvers/heatTransfer/buoyantBoussinesqSimpleFoam/TEqn.H b/applications/solvers/heatTransfer/buoyantBoussinesqSimpleFoam/TEqn.H
index c495e58..fec770c 100644
--- a/applications/solvers/heatTransfer/buoyantBoussinesqSimpleFoam/TEqn.H
+++ b/applications/solvers/heatTransfer/buoyantBoussinesqSimpleFoam/TEqn.H
@@ -9,7 +9,8 @@
         fvm::div(phi, T)
       - fvm::laplacian(alphaEff, T)
      ==
-        fvOptions(T)
+        radiation->ST(rhoCpRef, T)
+      + fvOptions(T)
     );
 
     TEqn.relax();
@@ -18,6 +19,8 @@
 
     TEqn.solve();
 
+    radiation->correct();
+
     fvOptions.correct(T);
 
     rhok = 1.0 - beta*(T - TRef);
diff --git a/applications/solvers/heatTransfer/buoyantBoussinesqSimpleFoam/buoyantBoussinesqSimpleFoam.C b/applications/solvers/heatTransfer/buoyantBoussinesqSimpleFoam/buoyantBoussinesqSimpleFoam.C
index 6dff329..75f0c23 100644
--- a/applications/solvers/heatTransfer/buoyantBoussinesqSimpleFoam/buoyantBoussinesqSimpleFoam.C
+++ b/applications/solvers/heatTransfer/buoyantBoussinesqSimpleFoam/buoyantBoussinesqSimpleFoam.C
@@ -48,6 +48,7 @@ Description
 #include "fvCFD.H"
 #include "singlePhaseTransportModel.H"
 #include "turbulentTransportModel.H"
+#include "radiationModel.H"
 #include "fvIOoptionList.H"
 #include "simpleControl.H"
 #include "fixedFluxPressureFvPatchScalarField.H"
@@ -63,6 +64,7 @@ int main(int argc, char *argv[])
     simpleControl simple(mesh);
 
     #include "createFields.H"
+    #include "createIncompressibleRadiationModel.H"
     #include "createMRF.H"
     #include "createFvOptions.H"
     #include "initContinuityErrs.H"
-- 
1.9.1

From bc22199c03d6c08511e2963c110ea0dc5abfab11 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Daniel=20Jasi=C5=84ski?= <daniel.jasinski@gmail.com>
Date: Thu, 29 Oct 2015 18:16:11 +0100
Subject: [PATCH] Reading incompressible radiation model properties consistent
 with Boussinesq readTransportProperties

---
 .../radiation/include/createIncompressibleRadiationModel.H            | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/src/thermophysicalModels/radiation/include/createIncompressibleRadiationModel.H b/src/thermophysicalModels/radiation/include/createIncompressibleRadiationModel.H
index ab7466d..13a8dc1 100644
--- a/src/thermophysicalModels/radiation/include/createIncompressibleRadiationModel.H
+++ b/src/thermophysicalModels/radiation/include/createIncompressibleRadiationModel.H
@@ -25,8 +25,8 @@
             )
         );
 
-        dimensionedScalar rhoRef(transportProperties.lookup("rhoRef"));
-        dimensionedScalar CpRef(transportProperties.lookup("CpRef"));
+        dimensionedScalar rhoRef("rhoRef",dimDensity,transportProperties);
+        dimensionedScalar CpRef("CpRef",dimSpecificHeatCapacity,transportProperties);
 
         rhoCpRef = rhoRef*CpRef;
     }
-- 
1.9.1

From 48c96a955a9391adcc238f6084346e1fdf6eab9f Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Daniel=20Jasi=C5=84ski?= <daniel.jasinski@gmail.com>
Date: Thu, 29 Oct 2015 13:45:41 +0100
Subject: [PATCH] Added radiative heat flux to wallHeatFlux utility

---
 .../wall/wallHeatFlux/createFields.H               | 16 ++++++++++++++++
 .../wall/wallHeatFlux/wallHeatFlux.C               | 22 ++++++++++++++--------
 2 files changed, 30 insertions(+), 8 deletions(-)

diff --git a/applications/utilities/postProcessing/wall/wallHeatFlux/createFields.H b/applications/utilities/postProcessing/wall/wallHeatFlux/createFields.H
index da7585b..02043ef 100644
--- a/applications/utilities/postProcessing/wall/wallHeatFlux/createFields.H
+++ b/applications/utilities/postProcessing/wall/wallHeatFlux/createFields.H
@@ -54,3 +54,19 @@ if (isA<fluidThermo>(thermo()))
         refCast<const fluidThermo>(thermo())
     );
 }
+
+//Read radiative heat flux if available
+volScalarField Qr
+(
+    IOobject
+    (
+        "Qr",
+        runTime.timeName(),
+        mesh,
+        IOobject::READ_IF_PRESENT,
+        IOobject::NO_WRITE
+    ),
+    mesh,
+    dimensionedScalar("Qr", dimMass/pow3(dimTime), 0.0)
+);
+
diff --git a/applications/utilities/postProcessing/wall/wallHeatFlux/wallHeatFlux.C b/applications/utilities/postProcessing/wall/wallHeatFlux/wallHeatFlux.C
index 8c80f0c..876affa 100644
--- a/applications/utilities/postProcessing/wall/wallHeatFlux/wallHeatFlux.C
+++ b/applications/utilities/postProcessing/wall/wallHeatFlux/wallHeatFlux.C
@@ -70,19 +70,25 @@ int main(int argc, char *argv[])
         const surfaceScalarField::GeometricBoundaryField& patchHeatFlux =
             heatFlux.boundaryField();
 
+        const volScalarField::GeometricBoundaryField& patchRadHeatFlux =
+            Qr.boundaryField();
+
+        const surfaceScalarField::GeometricBoundaryField& magSf =
+            mesh.magSf().boundaryField();
+
         Info<< "\nWall heat fluxes [W]" << endl;
         forAll(patchHeatFlux, patchi)
         {
             if (isA<wallFvPatch>(mesh.boundary()[patchi]))
             {
-                Info<< mesh.boundary()[patchi].name()
-                    << " "
-                    << gSum
-                       (
-                           mesh.magSf().boundaryField()[patchi]
-                          *patchHeatFlux[patchi]
-                       )
-                    << endl;
+                scalar convFlux = gSum(magSf[patchi]*patchHeatFlux[patchi]);
+
+                scalar radFlux = gSum(magSf[patchi]*patchRadHeatFlux[patchi]);
+
+                Info<< mesh.boundary()[patchi].name() << endl
+                    << "    convective: " << convFlux << endl
+                    << "    radiative:  " << radFlux << endl
+                    << "    total:      " << convFlux + radFlux << endl;
             }
         }
         Info<< endl;
-- 
1.9.1

From a54aefff4f5b80347d07e39668cdb6ebfd857ee9 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Daniel=20Jasi=C5=84ski?= <daniel.jasinski@gmail.com>
Date: Thu, 29 Oct 2015 15:08:43 +0100
Subject: [PATCH] basicThermo API for optional read

---
 .../basic/basicThermo/basicThermo.C                | 38 ++++++++++++++++++++++
 .../basic/basicThermo/basicThermo.H                | 12 ++++++-
 .../basic/basicThermo/basicThermoTemplates.C       |  9 +++--
 3 files changed, 55 insertions(+), 4 deletions(-)

diff --git a/src/thermophysicalModels/basic/basicThermo/basicThermo.C b/src/thermophysicalModels/basic/basicThermo/basicThermo.C
index e7264b0..e780379 100644
--- a/src/thermophysicalModels/basic/basicThermo/basicThermo.C
+++ b/src/thermophysicalModels/basic/basicThermo/basicThermo.C
@@ -276,6 +276,44 @@ Foam::autoPtr<Foam::basicThermo> Foam::basicThermo::New
     return New<basicThermo>(mesh, phaseName);
 }
 
+Foam::autoPtr<Foam::basicThermo> Foam::basicThermo::NewIfPresent
+(
+    const fvMesh& mesh,
+    const word& phaseName,
+    bool printInfo
+)
+{
+    IOobject thermoIO
+    (
+        phasePropertyName(dictName, phaseName),
+        mesh.time().constant(),
+        mesh,
+        IOobject::MUST_READ_IF_MODIFIED,
+        IOobject::NO_WRITE,
+        false
+    );
+
+    if(!thermoIO.headerOk())
+    {
+        return autoPtr<basicThermo>();
+    }
+
+    IOdictionary thermoDict
+    (
+        thermoIO
+    );
+
+    basicThermo::fvMeshConstructorTable::iterator cstrIter =
+        lookupThermo<basicThermo, basicThermo::fvMeshConstructorTable>
+        (
+            thermoDict,
+            basicThermo::fvMeshConstructorTablePtr_,
+            printInfo
+        );
+
+    return autoPtr<basicThermo>(cstrIter()(mesh, phaseName));
+}
+
 
 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
 
diff --git a/src/thermophysicalModels/basic/basicThermo/basicThermo.H b/src/thermophysicalModels/basic/basicThermo/basicThermo.H
index 9d4fc58..600e295 100644
--- a/src/thermophysicalModels/basic/basicThermo/basicThermo.H
+++ b/src/thermophysicalModels/basic/basicThermo/basicThermo.H
@@ -140,7 +140,8 @@ public:
         static typename Table::iterator lookupThermo
         (
             const dictionary& thermoDict,
-            Table* tablePtr
+            Table* tablePtr,
+            bool printInfo=true
         );
 
         //- Generic New for each of the related thermodynamics packages
@@ -167,6 +168,15 @@ public:
             const word& phaseName=word::null
         );
 
+        //- Specialisation of the Generic New for basicThermo
+        //- Returns empty pointer if dictionary not available
+        static autoPtr<basicThermo> NewIfPresent
+        (
+            const fvMesh&,
+            const word& phaseName=word::null,
+            bool printInfo=true
+        );
+
 
     //- Destructor
     virtual ~basicThermo();
diff --git a/src/thermophysicalModels/basic/basicThermo/basicThermoTemplates.C b/src/thermophysicalModels/basic/basicThermo/basicThermoTemplates.C
index 4521750..334ee89 100644
--- a/src/thermophysicalModels/basic/basicThermo/basicThermoTemplates.C
+++ b/src/thermophysicalModels/basic/basicThermo/basicThermoTemplates.C
@@ -31,7 +31,8 @@ template<class Thermo, class Table>
 typename Table::iterator Foam::basicThermo::lookupThermo
 (
     const dictionary& thermoDict,
-    Table* tablePtr
+    Table* tablePtr,
+    bool printInfo
 )
 {
     word thermoTypeName;
@@ -40,7 +41,8 @@ typename Table::iterator Foam::basicThermo::lookupThermo
     {
         const dictionary& thermoTypeDict(thermoDict.subDict("thermoType"));
 
-        Info<< "Selecting thermodynamics package " << thermoTypeDict << endl;
+        if(printInfo)
+            Info<< "Selecting thermodynamics package " << thermoTypeDict << endl;
 
         const int nCmpt = 7;
         const char* cmptNames[nCmpt] =
@@ -114,7 +116,8 @@ typename Table::iterator Foam::basicThermo::lookupThermo
     {
         thermoTypeName = word(thermoDict.lookup("thermoType"));
 
-        Info<< "Selecting thermodynamics package " << thermoTypeName << endl;
+        if(printInfo)
+            Info<< "Selecting thermodynamics package " << thermoTypeName << endl;
 
         typename Table::iterator cstrIter = tablePtr->find(thermoTypeName);
 
-- 
1.9.1

From 9e61f4095f51b7856ad7c68d2e4babdc42b76077 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Daniel=20Jasi=C5=84ski?= <daniel.jasinski@gmail.com>
Date: Thu, 29 Oct 2015 18:50:24 +0100
Subject: [PATCH] wallHeatFlux switches to Boussinesq when thermo package not
 found

---
 .../postProcessing/wall/wallHeatFlux/Make/options  |   5 +
 .../wall/wallHeatFlux/createFields.H               | 141 +++++++++++++++++----
 .../wall/wallHeatFlux/wallHeatFlux.C               |  44 +++++--
 3 files changed, 157 insertions(+), 33 deletions(-)

diff --git a/applications/utilities/postProcessing/wall/wallHeatFlux/Make/options b/applications/utilities/postProcessing/wall/wallHeatFlux/Make/options
index fdcd3d7..b134db7 100644
--- a/applications/utilities/postProcessing/wall/wallHeatFlux/Make/options
+++ b/applications/utilities/postProcessing/wall/wallHeatFlux/Make/options
@@ -1,6 +1,9 @@
 EXE_INC = \
     -I$(LIB_SRC)/TurbulenceModels/turbulenceModels/lnInclude \
+    -I$(LIB_SRC)/TurbulenceModels/incompressible/lnInclude \
     -I$(LIB_SRC)/TurbulenceModels/compressible/lnInclude \
+    -I$(LIB_SRC)/transportModels \
+    -I$(LIB_SRC)/transportModels/incompressible/singlePhaseTransportModel \
     -I$(LIB_SRC)/thermophysicalModels/specie/lnInclude \
     -I$(LIB_SRC)/thermophysicalModels/reactionThermo/lnInclude \
     -I$(LIB_SRC)/transportModels/compressible/lnInclude \
@@ -11,6 +14,8 @@ EXE_INC = \
 
 EXE_LIBS = \
     -lturbulenceModels \
+    -lincompressibleTransportModels \
+    -lincompressibleTurbulenceModels \
     -lcompressibleTurbulenceModels \
     -lreactionThermophysicalModels \
     -lgenericPatchFields \
diff --git a/applications/utilities/postProcessing/wall/wallHeatFlux/createFields.H b/applications/utilities/postProcessing/wall/wallHeatFlux/createFields.H
index 02043ef..7d8ee10 100644
--- a/applications/utilities/postProcessing/wall/wallHeatFlux/createFields.H
+++ b/applications/utilities/postProcessing/wall/wallHeatFlux/createFields.H
@@ -1,28 +1,14 @@
-autoPtr<basicThermo> thermo
-(
-    basicThermo::New(mesh)
-);
+autoPtr<surfaceScalarField> heatFlux;
 
-const volScalarField& h = thermo->he();
-
-// Register copy of thermo density
-volScalarField rho
+autoPtr<basicThermo> thermo
 (
-    IOobject
-    (
-        "rho",
-        runTime.timeName(),
-        mesh
-    ),
-    thermo->rho()
+    basicThermo::NewIfPresent(mesh,word::null,false)
 );
 
-// Construct turbulence model (if fluid)
 autoPtr<volVectorField> UPtr;
 autoPtr<surfaceScalarField> phiPtr;
-autoPtr<compressible::turbulenceModel> turbulence;
 
-if (isA<fluidThermo>(thermo()))
+if( (!thermo.valid()) || isA<fluidThermo>(thermo()) )
 {
     UPtr.reset
     (
@@ -39,22 +25,133 @@ if (isA<fluidThermo>(thermo()))
             mesh
         )
     );
-    const volVectorField& U = UPtr();
+}
+
+// Construct turbulence model (if fluid)
+autoPtr<volScalarField> rhoPtr;
+autoPtr<compressible::turbulenceModel> compressibleTurbulence;
+
+if (thermo.valid() && isA<fluidThermo>(thermo()))
+{
+    // Register copy of thermo density
+    rhoPtr.reset
+    (
+        new volScalarField
+        (
+            IOobject
+            (
+                "rho",
+                runTime.timeName(),
+                mesh
+            ),
+            thermo->rho()
+        )
+    );
+
+    const volScalarField rho = rhoPtr();
+    const volVectorField U = UPtr();
 
     #include "compressibleCreatePhi.H"
     // Copy phi to autoPtr. Rename to make sure copy is now registered as 'phi'.
     phi.rename("phiFluid");
     phiPtr.reset(new surfaceScalarField("phi", phi));
 
-    turbulence = compressible::turbulenceModel::New
+    compressibleTurbulence = compressible::turbulenceModel::New
     (
-        rho,
-        U,
+        rhoPtr(),
+        UPtr(),
         phiPtr(),
         refCast<const fluidThermo>(thermo())
     );
 }
 
+//If thermo package not found try Boussinesq
+autoPtr<volScalarField> T;
+autoPtr<volScalarField> alphat;
+autoPtr<volScalarField> rhok;
+
+autoPtr<incompressible::turbulenceModel> incompressibleTurbulence;
+autoPtr<singlePhaseTransportModel> laminarTransport;
+
+dimensionedScalar Pr("Pr", dimless, 1);
+dimensionedScalar rhoRef("rhoRef", dimDensity, 0);
+dimensionedScalar CpRef("CpRef", dimSpecificHeatCapacity, 0);
+
+if(!thermo.valid())
+{
+    T.reset
+    (
+        new volScalarField
+        (
+            IOobject
+            (
+                "T",
+                runTime.timeName(),
+                mesh,
+                IOobject::MUST_READ,
+                IOobject::NO_WRITE
+            ),
+            mesh
+        )
+    );
+
+    alphat.reset
+    (
+        new volScalarField
+        (
+            IOobject
+            (
+                "alphat",
+                runTime.timeName(),
+                mesh,
+                IOobject::MUST_READ,
+                IOobject::NO_WRITE
+            ),
+            mesh
+        )
+    );
+
+    rhok.reset
+    (
+        new volScalarField
+        (
+            IOobject
+            (
+                "rhok",
+                runTime.timeName(),
+                mesh,
+                IOobject::READ_IF_PRESENT,
+                IOobject::NO_WRITE
+            ),
+            mesh,
+            dimensionedScalar("rhok",dimless,1.0)
+        )
+    );
+
+    const volVectorField U = UPtr();
+
+    #include "createPhi.H"
+
+    phi.rename("phiFluid");
+    phiPtr.reset(new surfaceScalarField("phi", phi));
+
+    laminarTransport.reset
+    ( 
+        new singlePhaseTransportModel(UPtr(), phiPtr())
+    );
+
+    incompressibleTurbulence = incompressible::turbulenceModel::New
+    (
+        UPtr(), 
+        phiPtr(), 
+        laminarTransport()
+    );
+
+    laminarTransport().lookup("Pr") >> Pr;
+    laminarTransport().lookup("rhoRef") >> rhoRef;
+    laminarTransport().lookup("CpRef") >> CpRef;
+}
+
 //Read radiative heat flux if available
 volScalarField Qr
 (
diff --git a/applications/utilities/postProcessing/wall/wallHeatFlux/wallHeatFlux.C b/applications/utilities/postProcessing/wall/wallHeatFlux/wallHeatFlux.C
index b087331..096b5bf 100644
--- a/applications/utilities/postProcessing/wall/wallHeatFlux/wallHeatFlux.C
+++ b/applications/utilities/postProcessing/wall/wallHeatFlux/wallHeatFlux.C
@@ -33,6 +33,8 @@ Description
 
 #include "fvCFD.H"
 #include "turbulentFluidThermoModel.H"
+#include "turbulentTransportModel.H"
+#include "singlePhaseTransportModel.H"
 #include "solidThermo.H"
 #include "wallFvPatch.H"
 
@@ -55,20 +57,40 @@ int main(int argc, char *argv[])
 
         #include "createFields.H"
 
-        surfaceScalarField heatFlux
-        (
-            fvc::interpolate
+        if(thermo.valid())
+        {
+            heatFlux.reset
             (
+                new surfaceScalarField
                 (
-                    turbulence.valid()
-                  ? turbulence->alphaEff()()
-                  : thermo->alpha()
+                    fvc::interpolate
+                    (
+                        (
+                            compressibleTurbulence.valid()
+                          ? compressibleTurbulence->alphaEff()()
+                          : thermo->alpha()
+                        )
+                    )*fvc::snGrad(thermo->he())
                 )
-            )*fvc::snGrad(h)
-        );
+            );
+        }
+        else
+        {
+            heatFlux.reset
+            (
+                new surfaceScalarField
+                (
+                    CpRef*rhoRef*
+                    fvc::interpolate
+                    (
+                        rhok()*(incompressibleTurbulence->nu()/Pr + alphat())
+                    )*fvc::snGrad(T())
+                )
+            );
+        }
 
         const surfaceScalarField::GeometricBoundaryField& patchHeatFlux =
-            heatFlux.boundaryField();
+            heatFlux().boundaryField();
 
         const volScalarField::GeometricBoundaryField& patchRadHeatFlux =
             Qr.boundaryField();
@@ -102,7 +124,7 @@ int main(int argc, char *argv[])
                 mesh
             ),
             mesh,
-            dimensionedScalar("wallHeatFlux", heatFlux.dimensions(), 0.0)
+            dimensionedScalar("wallHeatFlux", heatFlux().dimensions(), 0.0)
         );
 
         forAll(wallHeatFlux.boundaryField(), patchi)
@@ -123,7 +145,7 @@ int main(int argc, char *argv[])
                     mesh
                 ),
                 mesh,
-                dimensionedScalar("totalWallHeatFlux", heatFlux.dimensions(), 0.0)
+                dimensionedScalar("totalWallHeatFlux", heatFlux().dimensions(), 0.0)
             );
 
             forAll(totalWallHeatFlux.boundaryField(), patchi)
-- 
1.9.1

From 1fbd03c8fb2aa11b164f4fb498f35bf1388019a0 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Daniel=20Jasi=C5=84ski?= <daniel.jasinski@gmail.com>
Date: Thu, 29 Oct 2015 14:00:12 +0100
Subject: [PATCH] Writing total heat flux

---
 .../wall/wallHeatFlux/wallHeatFlux.C               | 23 ++++++++++++++++++++++
 1 file changed, 23 insertions(+)

diff --git a/applications/utilities/postProcessing/wall/wallHeatFlux/wallHeatFlux.C b/applications/utilities/postProcessing/wall/wallHeatFlux/wallHeatFlux.C
index 876affa..b087331 100644
--- a/applications/utilities/postProcessing/wall/wallHeatFlux/wallHeatFlux.C
+++ b/applications/utilities/postProcessing/wall/wallHeatFlux/wallHeatFlux.C
@@ -111,6 +111,29 @@ int main(int argc, char *argv[])
         }
 
         wallHeatFlux.write();
+
+        if(Qr.headerOk())
+        {
+            volScalarField totalWallHeatFlux
+            (
+                IOobject
+                (
+                    "totalWallHeatFlux",
+                    runTime.timeName(),
+                    mesh
+                ),
+                mesh,
+                dimensionedScalar("totalWallHeatFlux", heatFlux.dimensions(), 0.0)
+            );
+
+            forAll(totalWallHeatFlux.boundaryField(), patchi)
+            {
+                totalWallHeatFlux.boundaryField()[patchi] = 
+                    patchHeatFlux[patchi] + patchRadHeatFlux[patchi];
+            }
+
+            totalWallHeatFlux.write();
+        }
     }
 
     Info<< "End" << endl;
-- 
1.9.1

From c4d42fbff748127b3166de4723e8a2c3f2a4a191 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Daniel=20Jasi=C5=84ski?= <daniel.jasinski@gmail.com>
Date: Thu, 29 Oct 2015 20:35:42 +0100
Subject: [PATCH] Final check based on pressure dimensions

---
 .../wall/wallHeatFlux/createFields.H               | 22 ++++++++++++++++++++++
 1 file changed, 22 insertions(+)

diff --git a/applications/utilities/postProcessing/wall/wallHeatFlux/createFields.H b/applications/utilities/postProcessing/wall/wallHeatFlux/createFields.H
index 7d8ee10..07370f1 100644
--- a/applications/utilities/postProcessing/wall/wallHeatFlux/createFields.H
+++ b/applications/utilities/postProcessing/wall/wallHeatFlux/createFields.H
@@ -25,6 +25,28 @@ if( (!thermo.valid()) || isA<fluidThermo>(thermo()) )
             mesh
         )
     );
+
+    //Final check if it is really compressible
+    volScalarField p
+    (
+        IOobject
+        (
+            "p",
+            runTime.timeName(),
+            mesh,
+            IOobject::MUST_READ,
+            IOobject::NO_WRITE
+        ),
+        mesh
+    );
+
+    if(p.dimensions() != dimPressure)
+    {
+        //- It is not a compressible pressure
+        //- probably an old thermodynamiProperties file
+        //- was found in the directory
+        thermo.clear();
+    }
 }
 
 // Construct turbulence model (if fluid)
-- 
1.9.1