diff --git a/applications/solvers/compressible/rhoPimpleFoam/pcEqn.H b/applications/solvers/compressible/rhoPimpleFoam/pcEqn.H index 55be311749..867fb9f641 100644 --- a/applications/solvers/compressible/rhoPimpleFoam/pcEqn.H +++ b/applications/solvers/compressible/rhoPimpleFoam/pcEqn.H @@ -109,19 +109,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); diff --git a/applications/solvers/compressible/rhoPimpleFoam/rhoPimpleDyMFoam/pEqn.H b/applications/solvers/compressible/rhoPimpleFoam/rhoPimpleDyMFoam/pEqn.H index d1dcc87b41..df14de0e76 100644 --- a/applications/solvers/compressible/rhoPimpleFoam/rhoPimpleDyMFoam/pEqn.H +++ b/applications/solvers/compressible/rhoPimpleFoam/rhoPimpleDyMFoam/pEqn.H @@ -101,7 +101,7 @@ if (pressureControl.limit(p)) p.correctBoundaryConditions(); } -thermo.correctRho(psi*p - psip0, rhoMin, rhoMax) ; +thermo.correctRho(psi*p - psip0, rhoMin, rhoMax); rho = thermo.rho(); {