diff --git a/src/finiteVolume/fields/fvPatchFields/derived/prghPressure/prghPressureFvPatchScalarField.C b/src/finiteVolume/fields/fvPatchFields/derived/prghPressure/prghPressureFvPatchScalarField.C
index 2338247e7883896265b82c85b5f768e7c666b18c..c86986f711090937b0ad9552b6dfde32f32df278 100644
--- a/src/finiteVolume/fields/fvPatchFields/derived/prghPressure/prghPressureFvPatchScalarField.C
+++ b/src/finiteVolume/fields/fvPatchFields/derived/prghPressure/prghPressureFvPatchScalarField.C
@@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
- \\ / A nd | Copyright (C) 2013 OpenFOAM Foundation
+ \\ / A nd | Copyright (C) 2013-2015 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@@ -152,7 +152,17 @@ void Foam::prghPressureFvPatchScalarField::updateCoeffs()
const uniformDimensionedVectorField& g =
db().lookupObject<uniformDimensionedVectorField>("g");
- operator==(p_ - rhop*((g.value() & patch().Cf())));
+ const uniformDimensionedScalarField& hRef =
+ db().lookupObject<uniformDimensionedScalarField>("hRef");
+
+ dimensionedScalar ghRef
+ (
+ mag(g.value()) > SMALL
+ ? g & (cmptMag(g.value())/mag(g.value()))*hRef
+ : dimensionedScalar("ghRef", g.dimensions()*dimLength, 0)
+ );
+
+ operator==(p_ - rhop*((g.value() & patch().Cf()) - ghRef.value()));
fixedValueFvPatchScalarField::updateCoeffs();
}
diff --git a/src/finiteVolume/fields/fvPatchFields/derived/prghPressure/prghPressureFvPatchScalarField.H b/src/finiteVolume/fields/fvPatchFields/derived/prghPressure/prghPressureFvPatchScalarField.H
index 9431ae325e26e1cd2f9bc3fd211cf35b345a79c5..bad53659f0059939a2515f2de75c86b988a77553 100644
--- a/src/finiteVolume/fields/fvPatchFields/derived/prghPressure/prghPressureFvPatchScalarField.H
+++ b/src/finiteVolume/fields/fvPatchFields/derived/prghPressure/prghPressureFvPatchScalarField.H
@@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
- \\ / A nd | Copyright (C) 2013 OpenFOAM Foundation
+ \\ / A nd | Copyright (C) 2013-2015 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@@ -32,7 +32,7 @@ Description
calculated as:
\f[
- p_rgh = p - \rho g h
+ p_rgh = p - \rho g (h - hRef)
\f]
where
@@ -40,8 +40,9 @@ Description
p_rgh | Pseudo hydrostatic pressure [Pa]
p | Static pressure [Pa]
h | Height in the opposite direction to gravity
+ hRef | Reference height in the opposite direction to gravity
\rho | density
- g | acceleration due to gravity [m/s2]
+ g | acceleration due to gravity [m/s^2]
\endtable
\heading Patch usage