From 486a82fad6547977ebb5a10069faceba36e5cccd Mon Sep 17 00:00:00 2001
From: sergio <sergio>
Date: Wed, 7 Jun 2017 16:26:43 +0100
Subject: [PATCH] Updating rho with thermo rho at the end of the pEq.H,
independently of simpleRho or limiting p.
---
.../solvers/compressible/rhoPimpleFoam/pEqn.H | 14 +++-----------
1 file changed, 3 insertions(+), 11 deletions(-)
diff --git a/applications/solvers/compressible/rhoPimpleFoam/pEqn.H b/applications/solvers/compressible/rhoPimpleFoam/pEqn.H
index ab89fcd8115..80d5d9b7946 100644
--- a/applications/solvers/compressible/rhoPimpleFoam/pEqn.H
+++ b/applications/solvers/compressible/rhoPimpleFoam/pEqn.H
@@ -98,19 +98,11 @@ K = 0.5*magSqr(U);
if (pressureControl.limit(p))
{
p.correctBoundaryConditions();
- thermo.correctRho(psi*p - psip0, rhoMin, rhoMax);
- rho = thermo.rho();
-}
-else if (pimple.SIMPLErho())
-{
- thermo.correctRho(psi*p - psip0, rhoMin, rhoMax);
- rho = thermo.rho();
-}
-else
-{
- thermo.correctRho(psi*p - psip0, rhoMin, rhoMax) ;
}
+thermo.correctRho(psi*p - psip0, rhoMin, rhoMax) ;
+rho = thermo.rho();
+
if (thermo.dpdt())
{
dpdt = fvc::ddt(p);
--
GitLab