关于icoFOAM中通量phi的更新

[dabiu]

icoFOAM中更新通量phi采用的如下语句（注：我理解这个语句）

if (piso.finalNonOrthogonalIter())
{
    phi = phiHbyA - pEqn.flux();
}

为什么不在下面这条语句的后面，直接用phi = fvc::flux(U)来更新通量，这样更新跟上面那条语句有什么区别吗？用phi=fvc::flux(U)会有什么问题吗？

U = HbyA - rAU*fvc::grad(p);
U.correctBoundaryConditions();

初学icoFOAM，求各位高手解答，OpenFOAM用的版本为OpenFOAM 7.

[李东岳]

这里面phi是严格守恒的，U不存在守恒的概念。所以你用phi = fvc::flux(U)，也可以，但是没有原来那样consistent

[dabiu]

@李东岳 感谢李老师的回复，谢谢

[AppleKiller]

@dabiu 在 关于icoFOAM中通量phi的更新 中说：

icoFOAM中更新通量phi采用的如下语句（注：我理解这个语句）

if (piso.finalNonOrthogonalIter())
{
    phi = phiHbyA - pEqn.flux();
}

为什么不在下面这条语句的后面，直接用phi = fvc::flux(U)来更新通量，这样更新跟上面那条语句有什么区别吗？用phi=fvc::flux(U)会有什么问题吗？

U = HbyA - rAU*fvc::grad(p);
U.correctBoundaryConditions();

初学icoFOAM，求各位高手解答，OpenFOAM用的版本为OpenFOAM 7.

我觉得这里的phi是经过Rhie-Chow插值修正过的面通量。