heatTransfer solvers: Added optional pressure reference

Added optional pressure reference pRef to p_rgh in buoyantPimpleFoam,
buoyantSimpleFoam and chtMultiRegionFoam which handles cases in which the
pressure variation is small compared to the pressure level more accurately.

The pRef value is provided in the optional constant/pRef file.

All tutorials and templates have been updated to use pRef as appropriate.

applications/solvers/heatTransfer/buoyantPimpleFoam/createFields.H

@@ -62,6 +62,7 @@ autoPtr<fluidThermophysicalTransportModel> thermophysicalTransport
 #include "readGravitationalAcceleration.H"
 #include "readhRef.H"
 #include "gh.H"
+#include "readpRef.H"
 
 
 Info<< "Reading field p_rgh\n" << endl;
@@ -79,7 +80,7 @@ volScalarField p_rgh
 );
 
 // Force p_rgh to be consistent with p
-p_rgh = p - rho*gh;
+p_rgh = p - rho*gh - pRef;
 
 pressureControl pressureControl
 (

applications/solvers/heatTransfer/buoyantPimpleFoam/pEqn.H

@@ -103,7 +103,7 @@ else
 
 phi = phiHbyA + p_rghEqn.flux();
 
-p = p_rgh + rho*gh;
+p = p_rgh + rho*gh + pRef;
 
 if (mesh.steady())
 {
@@ -146,7 +146,7 @@ if (adjustMass && !thermo.incompressible())
 {
     p += (initialMass - fvc::domainIntegrate(thermo.rho()))
         /fvc::domainIntegrate(psi);
-    p_rgh = p - rho*gh;
+    p_rgh = p - rho*gh - pRef;
     p.correctBoundaryConditions();
 }
 

applications/solvers/heatTransfer/buoyantSimpleFoam/createFields.H

@@ -57,6 +57,7 @@ autoPtr<fluidThermophysicalTransportModel> thermophysicalTransport
 #include "readGravitationalAcceleration.H"
 #include "readhRef.H"
 #include "gh.H"
+#include "readpRef.H"
 
 
 Info<< "Reading field p_rgh\n" << endl;
@@ -74,7 +75,7 @@ volScalarField p_rgh
 );
 
 // Force p_rgh to be consistent with p
-p_rgh = p - rho*gh;
+p_rgh = p - rho*gh - pRef;
 
 pressureControl pressureControl
 (

applications/solvers/heatTransfer/buoyantSimpleFoam/pEqn.H

@@ -77,7 +77,7 @@ else
 
 phi = phiHbyA + p_rghEqn.flux();
 
-p = p_rgh + rho*gh;
+p = p_rgh + rho*gh + pRef;
 
 #include "incompressible/continuityErrs.H"
 
@@ -98,7 +98,7 @@ if (closedVolume && !thermo.incompressible())
 {
     p += (initialMass - fvc::domainIntegrate(thermo.rho()))
         /fvc::domainIntegrate(psi);
-    p_rgh = p - rho*gh;
+    p_rgh = p - rho*gh - pRef;
     p.correctBoundaryConditions();
 }
 

applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/createFluidFields.H

@@ -5,6 +5,7 @@ PtrList<volVectorField> UFluid(fluidRegions.size());
 PtrList<surfaceScalarField> phiFluid(fluidRegions.size());
 PtrList<uniformDimensionedVectorField> gFluid(fluidRegions.size());
 PtrList<uniformDimensionedScalarField> hRefFluid(fluidRegions.size());
+PtrList<uniformDimensionedScalarField> pRefFluid(fluidRegions.size());
 PtrList<volScalarField> ghFluid(fluidRegions.size());
 PtrList<surfaceScalarField> ghfFluid(fluidRegions.size());
 PtrList<compressible::momentumTransportModel>
@@ -123,6 +124,24 @@ forAll(fluidRegions, i)
         )
     );
 
+    Info<< "    Adding to pRefFluid\n" << endl;
+    pRefFluid.set
+    (
+        i,
+        new uniformDimensionedScalarField
+        (
+            IOobject
+            (
+                "pRef",
+                runTime.constant(),
+                fluidRegions[i],
+                IOobject::READ_IF_PRESENT,
+                IOobject::NO_WRITE
+            ),
+            dimensionedScalar(dimPressure, 0)
+        )
+    );
+
     dimensionedScalar ghRef(- mag(gFluid[i])*hRefFluid[i]);
 
     Info<< "    Adding to ghFluid\n" << endl;
@@ -200,7 +219,7 @@ forAll(fluidRegions, i)
     );
 
     // Force p_rgh to be consistent with p
-    p_rghFluid[i] = thermoFluid[i].p() - rhoFluid[i]*ghFluid[i];
+    p_rghFluid[i] = thermoFluid[i].p() - rhoFluid[i]*ghFluid[i] - pRefFluid[i];
 
     fluidRegions[i].setFluxRequired(p_rghFluid[i].name());
 

applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/setRegionFluidFields.H

@@ -40,6 +40,7 @@
     volScalarField& p_rgh = p_rghFluid[i];
 
     const dimensionedVector& g = gFluid[i];
+    const dimensionedScalar& pRef = pRefFluid[i];
     const volScalarField& gh = ghFluid[i];
     const surfaceScalarField& ghf = ghfFluid[i];
 

etc/templates/singleFluidCHT/templates/0/fluid/p

@@ -25,13 +25,13 @@ boundaryField
     outlet
     {
         type            calculated;
-        value           uniform $pInitial;
+        value           $internalField;
     }
 
     inlet
     {
         type            calculated;
-        value           uniform $pInitial;
+        value           $internalField;
     }
 
     wall

etc/templates/singleFluidCHT/templates/0/fluid/p_rgh

@@ -16,22 +16,22 @@ FoamFile
 
 #include "$FOAM_CASE/constant/initialConditions"
 
-dimensions      [ 1 -1 -2 0 0 0 0 ];
+dimensions      [1 -1 -2 0 0 0 0];
 
-internalField   uniform $pInitial;
+internalField   uniform 0;
 
 boundaryField
 {
     outlet
     {
         type            fixedValue;
-        value           uniform $pInitial;
+        value           $internalField;
     }
 
     inlet
     {
         type            fixedFluxPressure;
-        value           uniform $pInitial;
+        value           $internalField;
     }
 
     wall

etc/templates/singleFluidCHT/templates/constant/fluid/pRef

@@ -0,0 +1,23 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  dev
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       uniformDimensionedScalarField;
+    location    "constant";
+    object      pRef;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#include "$FOAM_CASE/constant/initialConditions"
+
+dimensions      [1 -1 -2 0 0 0 0];
+value           $pInitial;
+
+// ************************************************************************* //

tutorials/heatTransfer/buoyantPimpleFoam/hotRoom/0/p_rgh

@@ -16,26 +16,26 @@ FoamFile
 
 dimensions      [1 -1 -2 0 0 0 0];
 
-internalField   uniform 1e5;
+internalField   uniform 0;
 
 boundaryField
 {
     floor
     {
         type            fixedFluxPressure;
-        value           uniform 1e5;
+        value           $internalField;
     }
 
     ceiling
     {
         type            fixedFluxPressure;
-        value           uniform 1e5;
+        value           $internalField;
     }
 
     fixedWalls
     {
         type            fixedFluxPressure;
-        value           uniform 1e5;
+        value           $internalField;
     }
 }
 

tutorials/heatTransfer/buoyantPimpleFoam/hotRoom/constant/pRef

@@ -0,0 +1,21 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  dev
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       uniformDimensionedScalarField;
+    location    "constant";
+    object      pRef;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [1 -1 -2 0 0 0 0];
+value           1e5;
+
+// ************************************************************************* //

tutorials/heatTransfer/buoyantSimpleFoam/buoyantCavity/0/p_rgh

@@ -17,32 +17,32 @@ FoamFile
 
 dimensions      [1 -1 -2 0 0 0 0];
 
-internalField   uniform 1e5;
+internalField   uniform 0;
 
 boundaryField
 {
     frontAndBack
     {
         type            fixedFluxPressure;
-        value           uniform 1e5;
+        value           $internalField;
     }
 
     topAndBottom
     {
         type            fixedFluxPressure;
-        value           uniform 1e5;
+        value           $internalField;
     }
 
     hot
     {
         type            fixedFluxPressure;
-        value           uniform 1e5;
+        value           $internalField;
     }
 
     cold
     {
         type            fixedFluxPressure;
-        value           uniform 1e5;
+        value           $internalField;
     }
 }
 

tutorials/heatTransfer/buoyantSimpleFoam/buoyantCavity/constant/pRef

@@ -0,0 +1,21 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  dev
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       uniformDimensionedScalarField;
+    location    "constant";
+    object      pRef;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [1 -1 -2 0 0 0 0];
+value           1e5;
+
+// ************************************************************************* //

tutorials/heatTransfer/buoyantSimpleFoam/circuitBoardCooling/0.orig/p

@@ -17,7 +17,7 @@ FoamFile
 
 dimensions      [1 -1 -2 0 0 0 0];
 
-internalField   uniform 101325;
+internalField   uniform 1e5;
 
 boundaryField
 {
@@ -27,28 +27,31 @@ boundaryField
     floor
     {
         type            calculated;
-        value           uniform 101325;
+        value           $internalField;
     }
+
     ceiling
     {
         type            calculated;
-        value           uniform 101325;
+        value           $internalField;
     }
+
     inlet
     {
         type            calculated;
-        value           uniform 101325;
+        value           $internalField;
     }
+
     outlet
     {
         type            calculated;
-        value           uniform 101325;
+        value           $internalField;
     }
+
     fixedWalls
     {
         type            empty;
     }
-
 }
 
 

tutorials/heatTransfer/buoyantSimpleFoam/circuitBoardCooling/0.orig/p_rgh

@@ -17,7 +17,7 @@ FoamFile
 
 dimensions      [1 -1 -2 0 0 0 0];
 
-internalField   uniform 101325;
+internalField   uniform 0;
 
 boundaryField
 {
@@ -27,31 +27,31 @@ boundaryField
     floor
     {
         type            fixedFluxPressure;
-        gradient        uniform 0;
+        value           $internalField;
-        value           uniform 101325;
     }
+
     ceiling
     {
         type            fixedFluxPressure;
-        gradient        uniform 0;
+        value           $internalField;
-        value           uniform 101325;
     }
+
     inlet
     {
         type            fixedFluxPressure;
-        gradient        uniform 0;
+        value           $internalField;
-        value           uniform 101325;
     }
+
     outlet
     {
         type            fixedValue;
-        value           uniform 101325;
+        value           $internalField;
     }
+
     fixedWalls
     {
         type            empty;
     }
-
 }
 
 

tutorials/heatTransfer/buoyantSimpleFoam/circuitBoardCooling/constant/pRef

@@ -0,0 +1,21 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  dev
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       uniformDimensionedScalarField;
+    location    "constant";
+    object      pRef;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [1 -1 -2 0 0 0 0];
+value           1e5;
+
+// ************************************************************************* //

tutorials/heatTransfer/buoyantSimpleFoam/comfortHotRoom/0/p_rgh

@@ -16,7 +16,7 @@ FoamFile
 
 dimensions      [1 -1 -2 0 0 0 0];
 
-internalField   uniform 1e5;
+internalField   uniform 0;
 
 boundaryField
 {

tutorials/heatTransfer/buoyantSimpleFoam/comfortHotRoom/constant/pRef

@@ -0,0 +1,21 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  dev
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       uniformDimensionedScalarField;
+    location    "constant";
+    object      pRef;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [1 -1 -2 0 0 0 0];
+value           1e5;
+
+// ************************************************************************* //

tutorials/heatTransfer/buoyantSimpleFoam/externalCoupledCavity/0/p_rgh

@@ -17,32 +17,32 @@ FoamFile
 
 dimensions      [1 -1 -2 0 0 0 0];
 
-internalField   uniform 1e5;
+internalField   uniform 0;
 
 boundaryField
 {
     frontAndBack
     {
         type            fixedFluxPressure;
-        value           uniform 1e5;
+        value           $internalField;
     }
 
     topAndBottom
     {
         type            fixedFluxPressure;
-        value           uniform 1e5;
+        value           $internalField;
     }
 
     hot
     {
         type            fixedFluxPressure;
-        value           uniform 1e5;
+        value           $internalField;
     }
 
     cold
     {
         type            fixedFluxPressure;
-        value           uniform 1e5;
+        value           $internalField;
     }
 }
 

tutorials/heatTransfer/buoyantSimpleFoam/externalCoupledCavity/constant/pRef

@@ -0,0 +1,21 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  dev
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       uniformDimensionedScalarField;
+    location    "constant";
+    object      pRef;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [1 -1 -2 0 0 0 0];
+value           1e5;
+
+// ************************************************************************* //

tutorials/heatTransfer/buoyantSimpleFoam/hotRadiationRoom/0/p

@@ -16,7 +16,7 @@ FoamFile
 
 dimensions      [1 -1 -2 0 0 0 0];
 
-internalField   uniform 100000;
+internalField   uniform 1e5;
 
 boundaryField
 {

tutorials/heatTransfer/buoyantSimpleFoam/hotRadiationRoom/0/p_rgh

@@ -16,32 +16,32 @@ FoamFile
 
 dimensions      [1 -1 -2 0 0 0 0];
 
-internalField   uniform 100000;
+internalField   uniform 0;
 
 boundaryField
 {
     floor
     {
         type            fixedFluxPressure;
-        value           uniform 100000;
+        value           $internalField;
     }
 
     ceiling
     {
         type            fixedFluxPressure;
-        value           uniform 100000;
+        value           $internalField;
     }
 
     fixedWalls
     {
         type            fixedFluxPressure;
-        value           uniform 100000;
+        value           $internalField;
     }
 
     box
     {
         type            fixedFluxPressure;
-        value           uniform 100000;
+        value           $internalField;
     }
 }
 

tutorials/heatTransfer/buoyantSimpleFoam/hotRadiationRoom/constant/pRef

@@ -0,0 +1,21 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  dev
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       uniformDimensionedScalarField;
+    location    "constant";
+    object      pRef;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [1 -1 -2 0 0 0 0];
+value           1e5;
+
+// ************************************************************************* //

tutorials/heatTransfer/buoyantSimpleFoam/hotRadiationRoomFvDOM/0/p

@@ -16,32 +16,32 @@ FoamFile
 
 dimensions      [1 -1 -2 0 0 0 0];
 
-internalField   uniform 100000;
+internalField   uniform 1e5;
 
 boundaryField
 {
     floor
     {
         type            calculated;
-        value           uniform 100000;
+        value           $internalField;
     }
 
     ceiling
     {
         type            calculated;
-        value           uniform 100000;
+        value           $internalField;
     }
 
     fixedWalls
     {
         type            calculated;
-        value           uniform 100000;
+        value           $internalField;
     }
 
     box
     {
         type            calculated;
-        value           uniform 100000;
+        value           $internalField;
     }
 }
 

tutorials/heatTransfer/buoyantSimpleFoam/hotRadiationRoomFvDOM/0/p_rgh

@@ -16,32 +16,32 @@ FoamFile
 
 dimensions      [1 -1 -2 0 0 0 0];
 
-internalField   uniform 100000;
+internalField   uniform 0;
 
 boundaryField
 {
     floor
     {
         type            fixedFluxPressure;
-        value           uniform 100000;
+        value           $internalField;
     }
 
     ceiling
     {
         type            fixedFluxPressure;
-        value           uniform 100000;
+        value           $internalField;
     }
 
     fixedWalls
     {
         type            fixedFluxPressure;
-        value           uniform 100000;
+        value           $internalField;
     }
 
     box
     {
         type            fixedFluxPressure;
-        value           uniform 100000;
+        value           $internalField;
     }
 }
 

tutorials/heatTransfer/buoyantSimpleFoam/hotRadiationRoomFvDOM/constant/pRef

@@ -0,0 +1,21 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  dev
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       uniformDimensionedScalarField;
+    location    "constant";
+    object      pRef;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [1 -1 -2 0 0 0 0];
+value           1e5;
+
+// ************************************************************************* //

tutorials/heatTransfer/buoyantSimpleFoam/hotRoomBoussinesq/0/p_rgh

@@ -23,19 +23,19 @@ boundaryField
     floor
     {
         type            fixedFluxPressure;
-        value           uniform 0;
+        value           $internalField;
     }
 
     ceiling
     {
         type            fixedFluxPressure;
-        value           uniform 0;
+        value           $internalField;
     }
 
     fixedWalls
     {
         type            fixedFluxPressure;
-        value           uniform 0;
+        value           $internalField;
     }
 }
 

tutorials/heatTransfer/buoyantSimpleFoam/iglooWithFridges/0/p_rgh

@@ -26,7 +26,7 @@ boundaryField
     wall
     {
         type            fixedFluxPressure;
-        value           uniform 0;
+        value           $internalField;
     }
 }
 

tutorials/heatTransfer/chtMultiRegionFoam/coolingSphere/templates/0/fluid/p_rgh

@@ -18,20 +18,20 @@ FoamFile
 
 dimensions      [ 1 -1 -2 0 0 0 0 ];
 
-internalField   uniform $pInitial;
+internalField   uniform 0;
 
 boundaryField
 {
     outlet
     {
         type            fixedValue;
-        value           uniform $pInitial;
+        value           uniform 0;
     }
 
     inlet
     {
         type            fixedFluxPressure;
-        value           uniform $pInitial;
+        value           uniform 0;
     }
 
     wall

tutorials/heatTransfer/chtMultiRegionFoam/coolingSphere/templates/constant/fluid/pRef

@@ -0,0 +1,23 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  dev
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       uniformDimensionedScalarField;
+    location    "constant";
+    object      pRef;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#include "$FOAM_CASE/constant/initialConditions"
+
+dimensions      [1 -1 -2 0 0 0 0];
+value           $pInitial;
+
+// ************************************************************************* //

tutorials/heatTransfer/chtMultiRegionFoam/heatExchanger/0/air/p_rgh

@@ -17,7 +17,7 @@ FoamFile
 
 dimensions      [1 -1 -2 0 0 0 0];
 
-internalField   uniform 100000;
+internalField   uniform 0;
 
 boundaryField
 {

tutorials/heatTransfer/chtMultiRegionFoam/heatExchanger/constant/air/pRef

@@ -0,0 +1,21 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  dev
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       uniformDimensionedScalarField;
+    location    "constant";
+    object      pRef;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [1 -1 -2 0 0 0 0];
+value           1e5;
+
+// ************************************************************************* //

tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/p_rgh

@@ -17,7 +17,7 @@ FoamFile
 
 dimensions      [ 1 -1 -2 0 0 0 0 ];
 
-internalField   uniform 1e6;
+internalField   uniform 0;
 
 boundaryField
 {

tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/pRef

@@ -0,0 +1,21 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  dev
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       uniformDimensionedScalarField;
+    location    "constant";
+    object      pRef;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [1 -1 -2 0 0 0 0];
+value           1e6;
+
+// ************************************************************************* //