From 2d41f24639a3cae82b31218965d3f80074c81c0e Mon Sep 17 00:00:00 2001 From: Henry Weller Date: Tue, 15 Aug 2017 10:01:57 +0100 Subject: [PATCH] chtMultiRegionFoam: Corrected pressure and density update order Resolves bug-report https://bugs.openfoam.org/view.php?id=2663 --- .../solvers/heatTransfer/chtMultiRegionFoam/fluid/pEqn.H | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/pEqn.H b/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/pEqn.H index 4dec94266..8773c152b 100644 --- a/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/pEqn.H +++ b/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/pEqn.H @@ -70,12 +70,12 @@ constrainPressure(p_rgh, rho, U, phiHbyA, rhorAUf, MRF); } } + p = p_rgh + rho*gh; + // Thermodynamic density update thermo.correctRho(psi*p - psip0); } -p = p_rgh + rho*gh; - // Update pressure time derivative if needed if (thermo.dpdt()) {