CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    求OpenFOAM中MULES算法的解析

    OpenFOAM
    3
    4
    1565
    正在加载更多帖子
    • 从旧到新
    • 从新到旧
    • 最多赞同
    回复
    • 在新帖中回复
    登录后回复
    此主题已被删除。只有拥有主题管理权限的用户可以查看。
    • 史
      史浩 讲师 最后由 编辑

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

      让我们随波逐流

      1 条回复 最后回复 回复 引用
      • yhdthu
        yhdthu 讲师 最后由 编辑

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

        长风破浪会有时,直挂云帆济沧海

        史 1 条回复 最后回复 回复 引用
        • 史
          史浩 讲师 @yhdthu 最后由 编辑

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

          让我们随波逐流

          1 条回复 最后回复 回复 引用
          • 李东岳
            李东岳 管理员 最后由 李东岳 编辑

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

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

            线上CFD课程开始报名:http://www.dyfluid.com/class.html

            CFD高性能服务器 http://dyfluid.com/servers.html

            1 条回复 最后回复 回复 引用
            • First post
              Last post