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