关于李东岳老师icoFoam求解器解析细节问题

[金石为开]

关于李东岳老师icoFoam求解器解析细节问题，我在网站上看到他的解析，看到细节有点不明白，想跟网友确认一下是笔误还是自己理解的错误呢？
还有后面的：加入压力梯度，我看一开始公式是有压力除以密度rho的，后面离散加的时候没有rho，这是咋回事？是漏写了吗还是就是这样子的？谢谢

[李东岳]

非常感谢！
已经将公式(13)和公式(14)更新。

[李东岳]

另外，在你发帖的时候，我应该把公式的速度重新换了一个符号（粗体），看截图还是之前的？是另存的pdf？

在前天更新的icoFoam解析我把速度矢量$\mathbf{U}$和其他文档的速度分量$u$进行了区分。并且强调了PISO循环。

[金石为开]

@李东岳 您的意思是变成了粗体的U来表示矢量吗？我的图就是截图您之前的内容。另外非常感谢你做的工作，让我们这些小研究生学习openfoam过程中少走弯路，真的很感谢。另外您以后会讲讲动网格，重叠网格相关东西可以吗？比如interDyMFoam的东西，我会用到动模型，动网格等。谢谢~~

[李东岳]

@金石为开
可以参考这里，和这里
interDyMFoam就是interFoam和DyMFoam的结合。

重叠网格OpenFOAM目前没有，我也没做过植入，不过extend版本的浸没边界法很有意思。

[金石为开]

@李东岳 我也第一次看这个piso循环。之前问题差不多明白了，想问下关于循环的。for (int corr=0; corr<nCorr; corr++)，for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)，if (nonOrth == nNonOrthCorr)等关于些循环步分别是什么意思呢？PISO循环的截止条件是什么？nCorr,nNonOrthCorr指什么呢？这个循环和非正交网格有关系？问题有点多哈，不过确实不理解。。。

[李东岳]

@金石为开
for (int corr=0; corr<nCorr; corr++)为PISO的内迭代，for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)为非正交修正。PISO循环无截止条件。不过Issa在1985年的论文里面已经指出2-3次即可。在我们使用中，需要强制定义2次或者3次。

有关非正交修正，需要码公式。在我处理完LES文档后说一下。也欢迎交流。

[金石为开]

@李东岳 恩恩，非常感谢，也就是这里nCorr的值也就是2-3的意思吧？

[李东岳]

@金石为开

是的，为2或者3或者其他。

[金石为开]

@李东岳 东岳老师我想说一下自己的理解，不知道对不对：最后速度更新的时候，U = HbyA - rAUfvc::grad(p)中HbyA是通过HbyA = rAUUEqn.H();//公式(15)得出的吧，也就是Ur,预测的速度而不是n+1步的速度对不？我还没看SIMPLE，不过按照我的理解是不是HbyA要随着速度变化不断更新？那内迭代的最终截止条件是什么呢？非常感谢！