CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    rhoCentralFoam 为什么做成半隐式?

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

      大家好,

      rhoCentralFoam 在处理inviscid flow时完全显式,对于粘性项做隐式处理。

      好奇为什么不直接做成 fully implicit? 因为作者在文中也提到考虑到步长问题,将粘性项做了隐式处理。还是说 central-upwind 的 KT/KNP 在全隐式时不可行?

      欢迎大家讨论,谢谢!

      1 条回复 最后回复 回复 引用
      • sibo
        sibo 最后由 编辑

        这里附上@东岳 前辈的 rhoCentralFoam解析
        以及rhoCentralFoam作者的 文章

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

          一方面对于存在激波这种双曲欧拉方程,显性离散虽然耦合性不强但对于间断的分辨率更好。对于有粘流动,本身激波要被抹掉,不需要捕获间断,隐形格式更加耦合

          另一方面KT算法将压力项和对流项放在一起,是要隐性求解速度呢还是压力呢?:xinxin3:

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

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

          sibo 1 条回复 最后回复 回复 引用
          • sibo
            sibo @李东岳 最后由 编辑

            @东岳 感谢前辈回复!

            看来如果要做全隐式,还需要一些其他方法捕获间断。比如改成压力基然后用pimple?

            rhoCentralFoam里面隐式的求解应该对速度修正吧?加入粘性项对速度的影响。

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

              我见过有文献做成piso+central/upwind格式的sci

              rhoCentralFoam里面隐式的求解应该对速度修正吧?加入粘性项对速度的影响。

              没有修正,并没有采用迭代的方式处理压力速度耦合

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

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

              sibo 1 条回复 最后回复 回复 引用
              • sibo
                sibo @李东岳 最后由 编辑

                @东岳 是啊我也看到了,改成压力基说是在low mach时更好。感觉各种拼凑。。

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