求OpenFOAM中MULES算法的解析



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


  • OpenFOAM讲师



  • @yhdthu 感谢老兄的文档,收获很多!帮助很大!因为没有计算流体力学的基础,看了三次才看完
    另外,有一个地方可能有问题:下图中lambda的上角标的+-号写反了。希望再确认一下:
    a6290ff3-8e64-4f1f-8b2b-c7e96e0cc3f6-image.png
    另外,这篇文章对于MULES算法通量修正算法写的很详细!与东岳写的interFoam解析搭配看,效果更佳!强推!
    interFoam解析
    为自己顶贴:duang:



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

    我这周刚把原始的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的最小值

    捕获.JPG


Log in to reply
 

CFD中文网 2016 - 2020 | 京ICP备15017992号-2