00001 fvVectorMatrix UEqn 00002 ( 00003 fvm::ddt(rho, U) 00004 + fvm::div(phi, U) 00005 + turbulence->divDevRhoReff(U) 00006 ); 00007 00008 UEqn.relax(); 00009 00010 solve 00011 ( 00012 UEqn 00013 == 00014 fvc::reconstruct 00015 ( 00016 ( 00017 - ghf*fvc::snGrad(rho) 00018 - fvc::snGrad(p_rgh) 00019 )*mesh.magSf() 00020 ), 00021 mesh.solver(U.select(oCorr == nOuterCorr - 1)) 00022 );