diff --git a/applications/solvers/compressible/rhoPimpleFoam/pEqn.H b/applications/solvers/compressible/rhoPimpleFoam/pEqn.H index d4c4f442b1..9c6785109f 100644 --- a/applications/solvers/compressible/rhoPimpleFoam/pEqn.H +++ b/applications/solvers/compressible/rhoPimpleFoam/pEqn.H @@ -85,6 +85,11 @@ else } } +if (pressureControl.limit(p)) +{ + p.correctBoundaryConditions(); +} + // Thermodynamic density update thermo.correctRho(psi*p - psip0); @@ -99,12 +104,7 @@ U.correctBoundaryConditions(); fvOptions.correct(U); K = 0.5*magSqr(U); -if (pressureControl.limit(p)) -{ - p.correctBoundaryConditions(); - rho = thermo.rho(); -} -else if (pimple.simpleRho()) +if (pimple.simpleRho()) { rho = thermo.rho(); } diff --git a/applications/solvers/compressible/rhoPimpleFoam/pcEqn.H b/applications/solvers/compressible/rhoPimpleFoam/pcEqn.H index c7c740eff8..62033c1ae8 100644 --- a/applications/solvers/compressible/rhoPimpleFoam/pcEqn.H +++ b/applications/solvers/compressible/rhoPimpleFoam/pcEqn.H @@ -99,6 +99,11 @@ else } } +if (pressureControl.limit(p)) +{ + p.correctBoundaryConditions(); +} + // Thermodynamic density update thermo.correctRho(psi*p - psip0); @@ -113,12 +118,7 @@ U.correctBoundaryConditions(); fvOptions.correct(U); K = 0.5*magSqr(U); -if (pressureControl.limit(p)) -{ - p.correctBoundaryConditions(); - rho = thermo.rho(); -} -else if (pimple.simpleRho()) +if (pimple.simpleRho()) { rho = thermo.rho(); }