43 pEqn.solve(
mesh.solver(
p.select(
pimple.finalInnerIter())));
45 if (
pimple.finalNonOrthogonalIter())
51#include "continuityErrs.H"
57Urel.correctBoundaryConditions();
adjustPhi(phiHbyA, U, p_rgh)
tmp< volScalarField > rAtUrel(rAUrel)
volScalarField rAUrel(1.0/UrelEqn.A())
tmp< fvVectorMatrix > tUrelEqn(fvm::ddt(Urel)+fvm::div(phi, Urel)+turbulence->divDevReff(Urel)+SRF->Su()==fvOptions(Urel))