求OpenFOAM中MULES算法的解析
-
哪里有OpenFOAM中MULES算法的详细解析?现在看源代码快看吐了,Henrik Rusche论文中关于这一部分解释和代码对应不起来,哪位大佬可以详细的讲解一下?
-
-
@yhdthu 感谢老兄的文档,收获很多!帮助很大!因为没有计算流体力学的基础,看了三次才看完
另外,有一个地方可能有问题:下图中lambda的上角标的+-号写反了。希望再确认一下:
另外,这篇文章对于MULES算法通量修正算法写的很详细!与东岳写的interFoam解析搭配看,效果更佳!强推!
interFoam解析
为自己顶贴 -
@史浩 他那个应该是个笔误。写反了。
我这周刚把原始的FCT算法植入进去了,OpenFOAM的算法跟原始的不太一样,主要在于
-
最大值和最小值的确定,这对应不同的方法,参考Zalesak1979 342页末
-
反扩散系数的确定(你贴的按个公式)
这也是他之所以叫做MULES的原因吧,不过思路都是一样的,Weller为了加快计算速度也优化的非常好。原始FCT已经完全可以对变量进行有界了,我对比原始的FCT和MULES差异非常小,原始FCT也不需要对$\lambda$进行迭代。我要用原始FCT验证几个算例分析一下。
在测试可压缩传输的时候,我发现可压缩FCT和不可压缩FCT的密度处理需要特殊处理。也就是说 http://www.cfd-china.com/topic/1029/mules算法中的const-rhotype-rho/3 帖子中的可压缩不能直接这么处理,要把rho和alpha分开。不能看做对rhoAlpha的传输。具体的我还在研究
下图FCT1 和FCT2对应原始FCT和MULES,差别微小。抽取的damBreak的alpha的最小值
-