UEqn.H 349 Bytes
Newer Older
henry's avatar
henry committed
1 2 3 4 5 6 7 8 9
fvVectorMatrix UEqn
(
    fvm::ddt(rho, U)
  + fvm::div(phi, U)
  + turbulence->divDevRhoReff(U)
);

UEqn.relax();

sergio's avatar
sergio committed
10 11 12
if (pimple.momentumPredictor())
{
    solve
henry's avatar
henry committed
13
    (
sergio's avatar
sergio committed
14 15 16
        UEqn
    ==
        fvc::reconstruct
17
        (
sergio's avatar
sergio committed
18 19 20 21 22 23 24
            (
            - ghf*fvc::snGrad(rho)
            - fvc::snGrad(p_rgh)
            )*mesh.magSf()
        )
    );
}