00001 // Solve the Momentum equation 00002 00003 tmp<fvVectorMatrix> UEqn 00004 ( 00005 fvm::div(phi, U) 00006 + turbulence->divDevReff(U) 00007 ); 00008 00009 UEqn().relax(); 00010 00011 eqnResidual = solve 00012 ( 00013 UEqn() == -fvc::grad(p) 00014 ).initialResidual(); 00015 00016 maxResidual = max(eqnResidual, maxResidual); 00017 00018 00019 // ************************ vim: set sw=4 sts=4 et: ************************ //