@aiweimo 你说的对，这里应该是一个typo，其实t和n都表示前一个时间步的值（已知时间步的值）

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

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 = phiHbyA - pEqn.flux()修正 ，为什么不直接在下面这条语句的后面直接用phi = fvc::flux(U)来更新呢，这样直接顺带把边界也更新了？

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

@dzw05 @李东岳 @wwzhao

同求alphaEqn.H的解析

怎么用OpenFOAM自带的SnappyHexMesh对边界网格进行局部加密，或者其他的网格剖分工具也可以