pimple算法的流程图
-
dear forumers:
写论文想用一下pimple算法的流程图,这样直观一些,大家有没有借我参考一下,我自己对这个算法感觉也是有点懵。。 -
@cccrrryyy 谢谢回答。我的理解是PIMPLE有一个外循环和内循环。内循环是压力修正,外循环是动量预测,但是松弛系数在哪里添加我不清楚。。
-
苏老师培训讲义里有这个流程图,如果你能搞到 话。
-
@Samuel-Tu 基本就是这样。参考rhoPimpleFoam的源代码。外循环外面的Compute rho(Move mesh下一个)那个地方严格来说有一个判据的,但好像不同版本之间这个地方也是改来改去。Solve energy之后是物性的更新,包括黏度、比热之类的,然后在压力方程的最开始更新密度(rho = thermo.rho())。
-
@cccrrryyy 谢谢了,这个回答对我很有帮助。因为我是对不可压缩流体进行了pimple算法的开发,所以没有rho和energy,也没有move mesh。我主要是关注u和p的两个循环对不对。这个图跟我理解的差不多。。我改一改嘿嘿
-
有个地方不太对劲,内循环里的New p guess,似乎应该是new u guess,这个地方的公式应该是重算了面通量。。
-
@Samuel-Tu 搞错了,最后一次迭代才重建面通量,这里我感觉就是new u guess
-
@Samuel-Tu 大佬,期待你搞一个附带公式的流程图
-
我毕业论文做的气泡流, 写了一个simple, piso,和pimple的框图,也不知道对不对,发在这里给大家参考一下。
https://de.overleaf.com/read/dnknwzrgmdkw!
SIMPLE
diagramm_of_simple.png
PISO
PIMPLE
-
@qingdong_wang 感谢分享!
-
@winsway_zero 😓,我可不是大佬,等我论文发了,就来补充
-
@qingdong_wang 看起来很对,如果有松弛的描述就更好了
-
@samuel-tu 为什么感觉PIMPLE才是真正的PISO算法,而PISO算法只是进行了多次压力和速度修正,在OpenFOAM中并没有判断一个时间步内收敛与否。从流程图上看也没有区别,https://marinecfd.xyz/post/piso-algorithm/ 这个帖子还对PISO的细节进行了分析。
-
@qingdong_wang 您好。我是小白,最近在恶补这些基础知识。对于您这个piso算法这个循环这里,我注意您这里到达了循环次数corrector number之后还进行了一个convergence判断。相当于两个要求都要满足。但是我看有些博主说法是到达循环次数或者到达收敛判断就能输出了(这个部分似乎和pimple图是意思一样)所以应该是那种符合piso算法呢? corrector number和convergence符合一个即可输出还是两个都需要符合
-
@qingdong_wang 还有这个pimple算法这,按我理解不应该是速度预测,压力修正,速度修正,迭代不足回到压力修正这种?为什么您这个是速度直接求出来了,然后直接就压力值了,迭代不成功回压力修正
7/18