ENH: Adding pRef handling in pEqn.H for incompressible eqOfState in buoyantSimpleFoam

2025-11-28 03:28:01 +00:00 · 2015-12-17 15:02:01 -08:00
parent fb504d8975
commit 23b99b5e3b
1 changed files with 15 additions and 3 deletions
--- a/applications/solvers/heatTransfer/buoyantSimpleFoam/pEqn.H
+++ b/applications/solvers/heatTransfer/buoyantSimpleFoam/pEqn.H
@ -70,10 +70,22 @@

    // For closed-volume cases adjust the pressure level
    // to obey overall mass continuity
-    if (closedVolume && compressible)
+    if (closedVolume)
    {
-        p += (initialMass - fvc::domainIntegrate(psi*p))
-            /fvc::domainIntegrate(psi);
+        if(!compressible)
+        {
+            p += dimensionedScalar
+            (
+                "p",
+                p.dimensions(),
+                pRefValue - getRefCellValue(p, pRefCell)
+            );
+        }
+        else
+        {
+            p += (initialMass - fvc::domainIntegrate(psi*p))
+                /fvc::domainIntegrate(psi);
+        }
        p_rgh = p - rho*gh;
    }