@李东岳 在 openFoam与Fluent模拟简单的Laval蒸汽喷管 中说:
上面这个图是OpenFOAM-9算的。使用OpenFOAM-10计算的话,直接发散掉了。经对比应该是pEqn.H的方程,OpenFOAM-10里面这一段代码有变动:
const fvScalarMatrix divPhidp(fvm::div(phid, p));
phiHbyA -= divPhidp.flux();
我把他改回去后,OpenFOAM-10可以收敛。但是结果异常。
针对这部分内容,今天碰到一个bug处理的帖子,确实是这部分原因导致的。在新版OpenFOAM-10中已经删掉
https://github.com/OpenFOAM/OpenFOAM-10/commit/910a10bf395da16de85ad12e88bc5cb9abb4f7e0