From 9c3a404533e4454eb29d489dc92028ec8b5d82d8 Mon Sep 17 00:00:00 2001
From: Henry Weller <http://cfd.direct>
Date: Fri, 15 Sep 2023 08:13:13 +0100
Subject: [PATCH] multiphaseEuler::cellPressureCorrector: Added p_rgh update
 before PISO loop

to ensure consistency between p_rgh, p and the current phase densities updated
after the phase energy solution.
---
 .../multiphaseEuler/cellPressureCorrector.C              | 9 +++++----
 1 file changed, 5 insertions(+), 4 deletions(-)

diff --git a/applications/modules/multiphaseEuler/multiphaseEuler/cellPressureCorrector.C b/applications/modules/multiphaseEuler/multiphaseEuler/cellPressureCorrector.C
index ef49235fa9..b7d26cfb27 100644
--- a/applications/modules/multiphaseEuler/multiphaseEuler/cellPressureCorrector.C
+++ b/applications/modules/multiphaseEuler/multiphaseEuler/cellPressureCorrector.C
@@ -45,6 +45,11 @@ void Foam::solvers::multiphaseEuler::cellPressureCorrector()
 {
     volScalarField& p(p_);
 
+    volScalarField rho("rho", fluid.rho());
+
+    // Correct p_rgh for consistency with the current density
+    p_rgh = p - rho*buoyancy.gh - buoyancy.pRef;
+
     // Face volume fractions
     PtrList<surfaceScalarField> alphafs(phases.size());
     forAll(phases, phasei)
@@ -108,8 +113,6 @@ void Foam::solvers::multiphaseEuler::cellPressureCorrector()
         fluid.invADVs(As, HVms, invADVs, invADVfs);
     }
 
-    volScalarField rho("rho", fluid.rho());
-
     // Explicit force fluxes
     PtrList<surfaceScalarField> alphaByADfs;
     PtrList<surfaceScalarField> FgByADfs;
@@ -219,7 +222,6 @@ void Foam::solvers::multiphaseEuler::cellPressureCorrector()
         }
     }
 
-
     // --- Pressure corrector loop
     while (pimple.correct())
     {
@@ -511,7 +513,6 @@ void Foam::solvers::multiphaseEuler::cellPressureCorrector()
         // Correct p_rgh for consistency with p and the updated densities
         rho = fluid.rho();
         p_rgh = p - rho*buoyancy.gh - buoyancy.pRef;
-        p_rgh.correctBoundaryConditions();
     }
 
     UEqns.clear();