UEqn.H 385 Bytes
Newer Older
Andrew Heather's avatar
Andrew Heather committed
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
    // Solve the momentum equation

    tmp<fvVectorMatrix> UEqn
    (
        fvm::ddt(U)
      + fvm::div(phi, U)
      + turbulence->divDevReff(U)
    );

    UEqn().relax();

    solve
    (
        UEqn()
      ==
       -fvc::reconstruct
        (
            (
                fvc::snGrad(pd)
              - betaghf*fvc::snGrad(T)
            ) * mesh.magSf()
        )
    );