[alundilong]

This might be silly question, but it confused me a lot..
these snippet code comes from buoyantSimpleFoam of OpenFOAM-2.1.1
3 tmp<fvVectorMatrix> UEqn
4 (
5 fvm::div(phi, U)
6 + turbulence->divDevRhoReff(U)
7 );
8
9 UEqn().relax();
10
11 if (simple.momentumPredictor())
12 {
13 solve
14 (
15 UEqn()
16 ==
17 fvc::reconstruct
18 (
19 (
20 - ghf*fvc::snGrad(rho)
21 - fvc::snGrad(p_rgh)
22 )*mesh.magSf()
23 )
24 );
25 }
shouldn't be this?
3 tmp<fvVectorMatrix> UEqn
4 (
5 fvm::div(phi, U)
6 + turbulence->divDevRhoReff(U)
7 );
8
9 UEqn().relax();
10
11 if (simple.momentumPredictor())
12 {
13 solve
14 (
15 UEqn()
16 ==
17 fvc::reconstruct
18 (
19 (
20 - ghf*fvc::snGrad(rho)
21 - fvc::snGrad(p_rgh)
22 )*mesh.magSf()
23 )
24 + rho*g
25 );
25 }