UEqn.H

// Solve the Momentum equation
 
MRF.correctBoundaryVelocity(U);
 
tmp<fvVectorMatrix> tUEqn
(
    fvm::div(phi, U)
    + MRF.DDt(rho, U)
    + turbulence->divDevRhoReff(U)
    ==
    fvOptions(rho, U)
);
fvVectorMatrix& UEqn = tUEqn.ref();
 
UEqn.relax();
 
fvOptions.constrain(UEqn);
 
if (simple.momentumPredictor())
{
    uresidual_v = solve(UEqn == -fvc::grad(p)).initialResidual();
    fvOptions.correct(U);
}