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