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