Buoyant solvers: Add special handling for ghRef in the case g = (0 0 0)

applications/solvers/combustion/fireFoam/createFields.H

@@ -103,11 +103,7 @@
 
     #include "readGravitationalAcceleration.H"
     #include "readhRef.H"
-
+    #include "gh.H"
-    Info<< "Calculating field g.h\n" << endl;
-    dimensionedScalar ghRef(g & (cmptMag(g.value())/mag(g.value()))*hRef);
-    volScalarField gh("gh", (g & mesh.C()) - ghRef);
-    surfaceScalarField ghf("ghf", (g & mesh.Cf()) - ghRef);
 
 
     volScalarField p_rgh

applications/solvers/combustion/reactingFoam/rhoReactingBuoyantFoam/createFields.H

@@ -65,12 +65,7 @@ reaction->setTurbulence(turbulence());
 
 #include "readGravitationalAcceleration.H"
 #include "readhRef.H"
-
+#include "gh.H"
-Info<< "Calculating field g.h\n" << endl;
-dimensionedScalar ghRef(g & (cmptMag(g.value())/mag(g.value()))*hRef);
-volScalarField gh("gh", (g & mesh.C()) - ghRef);
-surfaceScalarField ghf("ghf", (g & mesh.Cf()) - ghRef);
-
 
 Info<< "Reading field p_rgh\n" << endl;
 volScalarField p_rgh

applications/solvers/heatTransfer/buoyantBoussinesqPimpleFoam/createFields.H

@@ -82,11 +82,7 @@
 
     #include "readGravitationalAcceleration.H"
     #include "readhRef.H"
-
+    #include "gh.H"
-    Info<< "Calculating field g.h\n" << endl;
-    dimensionedScalar ghRef(g & (cmptMag(g.value())/mag(g.value()))*hRef);
-    volScalarField gh("gh", (g & mesh.C()) - ghRef);
-    surfaceScalarField ghf("ghf", (g & mesh.Cf()) - ghRef);
 
 
     volScalarField p

applications/solvers/heatTransfer/buoyantBoussinesqSimpleFoam/createFields.H

@@ -82,11 +82,7 @@
 
     #include "readGravitationalAcceleration.H"
     #include "readhRef.H"
-
+    #include "gh.H"
-    Info<< "Calculating field g.h\n" << endl;
-    dimensionedScalar ghRef(g & (cmptMag(g.value())/mag(g.value()))*hRef);
-    volScalarField gh("gh", (g & mesh.C()) - ghRef);
-    surfaceScalarField ghf("ghf", (g & mesh.Cf()) - ghRef);
 
 
     volScalarField p

applications/solvers/lagrangian/reactingParcelFilmFoam/createFields.H

@@ -83,11 +83,7 @@
 
     #include "readGravitationalAcceleration.H"
     #include "readhRef.H"
-
+    #include "gh.H"
-    Info<< "Calculating field g.h\n" << endl;
-    dimensionedScalar ghRef(g & (cmptMag(g.value())/mag(g.value()))*hRef);
-    volScalarField gh("gh", (g & mesh.C()) - ghRef);
-    surfaceScalarField ghf("ghf", (g & mesh.Cf()) - ghRef);
 
 
     volScalarField p_rgh

applications/solvers/multiphase/compressibleInterFoam/createFields.H

@@ -67,11 +67,7 @@
 
     #include "readGravitationalAcceleration.H"
     #include "readhRef.H"
-
+    #include "gh.H"
-    Info<< "Calculating field g.h\n" << endl;
-    dimensionedScalar ghRef(g & (cmptMag(g.value())/mag(g.value()))*hRef);
-    volScalarField gh("gh", (g & mesh.C()) - ghRef);
-    surfaceScalarField ghf("ghf", (g & mesh.Cf()) - ghRef);
 
 
     // Mass flux

applications/solvers/multiphase/compressibleMultiphaseInterFoam/createFields.H

@@ -51,11 +51,7 @@
 
     #include "readGravitationalAcceleration.H"
     #include "readhRef.H"
-
+    #include "gh.H"
-    Info<< "Calculating field g.h\n" << endl;
-    dimensionedScalar ghRef(g & (cmptMag(g.value())/mag(g.value()))*hRef);
-    volScalarField gh("gh", (g & mesh.C()) - ghRef);
-    surfaceScalarField ghf("ghf", (g & mesh.Cf()) - ghRef);
 
 
     // Construct compressible turbulence model

applications/solvers/multiphase/driftFluxFoam/createFields.H

@@ -93,11 +93,7 @@
 
     #include "readGravitationalAcceleration.H"
     #include "readhRef.H"
-
+    #include "gh.H"
-    Info<< "Calculating field g.h\n" << endl;
-    dimensionedScalar ghRef(g & (cmptMag(g.value())/mag(g.value()))*hRef);
-    volScalarField gh("gh", (g & mesh.C()) - ghRef);
-    surfaceScalarField ghf("ghf", (g & mesh.Cf()) - ghRef);
 
 
     volScalarField p

applications/solvers/multiphase/interFoam/createFields.H

@@ -76,13 +76,11 @@
         incompressible::turbulenceModel::New(U, phi, mixture)
     );
 
+
     #include "readGravitationalAcceleration.H"
     #include "readhRef.H"
+    #include "gh.H"
 
-    Info<< "Calculating field g.h\n" << endl;
-    dimensionedScalar ghRef(g & (cmptMag(g.value())/mag(g.value()))*hRef);
-    volScalarField gh("gh", (g & mesh.C()) - ghRef);
-    surfaceScalarField ghf("ghf", (g & mesh.Cf()) - ghRef);
 
     volScalarField p
     (

applications/solvers/multiphase/interFoam/interMixingFoam/createFields.H

@@ -81,11 +81,7 @@
 
     #include "readGravitationalAcceleration.H"
     #include "readhRef.H"
-
+    #include "gh.H"
-    Info<< "Calculating field g.h\n" << endl;
-    dimensionedScalar ghRef(g & (cmptMag(g.value())/mag(g.value()))*hRef);
-    volScalarField gh("gh", (g & mesh.C()) - ghRef);
-    surfaceScalarField ghf("ghf", (g & mesh.Cf()) - ghRef);
 
     volScalarField p
     (

applications/solvers/multiphase/interPhaseChangeFoam/createFields.H

@@ -69,11 +69,7 @@
 
     #include "readGravitationalAcceleration.H"
     #include "readhRef.H"
-
+    #include "gh.H"
-    Info<< "Calculating field g.h\n" << endl;
-    dimensionedScalar ghRef(g & (cmptMag(g.value())/mag(g.value()))*hRef);
-    volScalarField gh("gh", (g & mesh.C()) - ghRef);
-    surfaceScalarField ghf("ghf", (g & mesh.Cf()) - ghRef);
 
 
     volScalarField p

applications/solvers/multiphase/multiphaseInterFoam/createFields.H

@@ -55,11 +55,7 @@
 
     #include "readGravitationalAcceleration.H"
     #include "readhRef.H"
-
+    #include "gh.H"
-    Info<< "Calculating field g.h\n" << endl;
-    dimensionedScalar ghRef(g & (cmptMag(g.value())/mag(g.value()))*hRef);
-    volScalarField gh("gh", (g & mesh.C()) - ghRef);
-    surfaceScalarField ghf("ghf", (g & mesh.Cf()) - ghRef);
 
 
     volScalarField p

applications/solvers/multiphase/twoLiquidMixingFoam/createFields.H

@@ -72,11 +72,7 @@
 
     #include "readGravitationalAcceleration.H"
     #include "readhRef.H"
-
+    #include "gh.H"
-    Info<< "Calculating field g.h\n" << endl;
-    dimensionedScalar ghRef(g & (cmptMag(g.value())/mag(g.value()))*hRef);
-    volScalarField gh("gh", (g & mesh.C()) - ghRef);
-    surfaceScalarField ghf("ghf", (g & mesh.Cf()) - ghRef);
 
 
     volScalarField p

applications/solvers/multiphase/twoPhaseEulerFoam/createFields.H

@@ -30,12 +30,7 @@
         fluid.lookup("pMin")
     );
 
-
+    #include "gh.H"
-    Info<< "Calculating field g.h\n" << endl;
-    dimensionedScalar ghRef(g & (cmptMag(g.value())/mag(g.value()))*hRef);
-    volScalarField gh("gh", (g & mesh.C()) - ghRef);
-    surfaceScalarField ghf("ghf", (g & mesh.Cf()) - ghRef);
-
 
     rhoThermo& thermo1 = phase1.thermo();
     rhoThermo& thermo2 = phase2.thermo();

src/finiteVolume/cfdTools/general/include/gh.H

@@ -0,0 +1,9 @@
+    Info<< "Calculating field g.h\n" << endl;
+    dimensionedScalar ghRef
+    (
+        mag(g.value()) > SMALL
+      ? g & (cmptMag(g.value())/mag(g.value()))*hRef
+      : dimensionedScalar("ghRef", g.dimensions()*dimLength, 0)
+    );
+    volScalarField gh("gh", (g & mesh.C()) - ghRef);
+    surfaceScalarField ghf("ghf", (g & mesh.Cf()) - ghRef);