icoFoam源代码解读

[风云5091]

各位老师，大佬晚上好！openfoam小白最近在阅读icoFoam源代码时遇到了一点问题，根据李老师对icoFoam求解器的理论推导，它的压力泊松方程如下：

在源代码中首先定义了phiHbyA,应该是HbyA的面通量:

但是,它的压力泊松方程里面的用的是phiHbyA而不是HbyA,这里我不是很能理解：

另外按理说piso算法每算一次U、P，按理来说应该要更新HbyA一次：

但是源代码里面更新的好像是phiHbyA（速度变了，且是个通量）
在我的理解里，这段代码每次进行piso算法时，HbyA似乎没变，毕竟它是从之前的速度方程那里得到的？（对UEqn.A这个函数不太了解）

但是piso算完以后的U用的又是和理论推导的结果一致，用的是HbyA，而不是phiHbyA:

恳请各位老师，大佬解惑，谢谢各位！

[李东岳]

我才看到你的帖子。

它的压力泊松方程里面的用的是phiHbyA而不是HbyA,这里我不是很能理解：

如果HbyA到phiHbyA里面没有那个ddtcorr操作，那么你可以给HbyA。但更OpenFOAM的写法是，div()里面加面通量，这可以很好的理解为满足守恒关系。

你要注意HbyA是体矢量场，phiHbyA是面标量场。两个区别很大

[风云5091]

@李东岳 哦哦原来是这样，了解了，谢谢李老师！