@wsxfyy 谢谢，我试一下

@李东岳 好的，谢谢李老师，我试一下

@李东岳 谢谢李老师，壁面设置一种润湿角我是知道怎么设置的，在一个壁面上设置不同的润湿角，我能想到的是把一个壁面分为好几个壁面，然后再分别设置润湿角，但是这样比较麻烦。我想问的是有没有方法就是只有一个壁面，然后把润湿角设置成位置的函数

使用InterFoam模拟液滴在壁面的润湿性，如果在一个壁面上设置不同的润湿性呢，比如一个表面0到50 cm润湿角是60度， 50到100cm润湿角是120度，100到150cm润湿角是60度，谢谢。

这个结果很正常，毕竟很难严格保界，但是这里的-1.1e-11不算越界，毕竟这个数基本都是0

@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