00001 fvVectorMatrix UEqn 00002 ( 00003 fvm::ddt(rho, U) 00004 + fvm::div(phi, U) 00005 + turbulence->divDevRhoReff(U) 00006 ); 00007 00008 if (momentumPredictor) 00009 { 00010 solve(UEqn == -fvc::grad(p)); 00011 }