5 - fvm::Sp(fvc::div(phi), h)
6 - fvm::laplacian(turbulence->alphaEff(), h)
8 fvc::div(phi/fvc::interpolate(rho), p, "div(U,p)")
9 - p*fvc::div(phi/fvc::interpolate(rho))
14 eqnResidual = hEqn.solve().initialResidual();
15 maxResidual = max(eqnResidual, maxResidual);