求OpenFOAM中MULES算法的解析

[史浩]

哪里有OpenFOAM中MULES算法的详细解析？现在看源代码快看吐了，Henrik Rusche论文中关于这一部分解释和代码对应不起来，哪位大佬可以详细的讲解一下？

[yhdthu]

https://zhuanlan.zhihu.com/p/25025837

[史浩]

@yhdthu 感谢老兄的文档，收获很多！帮助很大！因为没有计算流体力学的基础，看了三次才看完
另外，有一个地方可能有问题：下图中lambda的上角标的+-号写反了。希望再确认一下：

另外，这篇文章对于MULES算法通量修正算法写的很详细！与东岳写的interFoam解析搭配看，效果更佳！强推！
interFoam解析
为自己顶贴

[东岳]

@史浩 他那个应该是个笔误。写反了。

我这周刚把原始的FCT算法植入进去了，OpenFOAM的算法跟原始的不太一样，主要在于

1. 最大值和最小值的确定，这对应不同的方法，参考Zalesak1979 342页末

2. 反扩散系数的确定（你贴的按个公式）

这也是他之所以叫做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的最小值