buoyantBoussinesqSimpleFoam的问题

[Vortex]

最近在阅读buoyantBoussinesqSimpleFoam代码遇到困惑，为什么公式（11）右端前两项对应的代码不使用

-fvc::grad(p_rgh))-gh*fvc::grad(rhok))

而使用的是面法向梯度snGrad的数据

fvc::reconstruct
            (
                (
                  - ghf*fvc::snGrad(rhok)
                  - fvc::snGrad(p_rgh)
                )*mesh.magSf()
            )

请各位大佬赐教~

[东岳]

http://dyfluid.com/docs/RhieChow.html 最后一段解释
主要是抹掉震荡

[Vortex]

@东岳 谢谢李老师~