求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的最小值
-
