UEqn.H

Go to the documentation of this file.

    MRF.correctBoundaryVelocity(U);
 
    tmp<fvVectorMatrix> UEqn
    (
        fvm::div(phi, U)
      + MRF.DDt(rho, U)
      + turbulence->divDevRhoReff(U)
     ==
        rho.dimensionedInternalField()*g
      + parcels.SU(U)
      + fvOptions(rho, U)
    );
 
    UEqn().relax();
 
    fvOptions.constrain(UEqn());
 
    solve(UEqn() == -fvc::grad(p));
 
    fvOptions.correct(U);