diff --git a/applications/solvers/heatTransfer/buoyantPimpleFoam/pEqn.H b/applications/solvers/heatTransfer/buoyantPimpleFoam/pEqn.H index a3cddc05d8..2ea2dfcb72 100644 --- a/applications/solvers/heatTransfer/buoyantPimpleFoam/pEqn.H +++ b/applications/solvers/heatTransfer/buoyantPimpleFoam/pEqn.H @@ -98,6 +98,9 @@ else phi = phiHbyA + p_rghEqn.flux(); +// Explicitly relax pressure for momentum corrector +p_rgh.relax(); + p = p_rgh + rho*gh + pRef; if (mesh.steady()) @@ -120,9 +123,6 @@ else #include "compressibleContinuityErrs.H" } -// Explicitly relax pressure for momentum corrector -p_rgh.relax(); - // Correct the momentum source with the pressure gradient flux // calculated from the relaxed pressure U = HbyA + rAU*fvc::reconstruct((phig + p_rghEqn.flux())/rhorAUf);