Skip to content
Snippets Groups Projects
UEqn.H 343 B
Newer Older
// Solve the Momentum equation

Henry's avatar
Henry committed
MRF.correctBoundaryVelocity(U);

tmp<fvVectorMatrix> UEqn
(
Henry's avatar
Henry committed
    fvm::ddt(U) + fvm::div(phi, U)
  + MRF.DDt(U)
  + turbulence->divDevReff(U)
    fvOptions(U)
fvOptions.constrain(UEqn());
if (pimple.momentumPredictor())
    solve(UEqn() == -fvc::grad(p));
    fvOptions.correct(U);