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