rhoCentralFoam 为什么做成半隐式?



  • 大家好,

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

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

    欢迎大家讨论,谢谢!



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



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

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



  • @东岳 感谢前辈回复!

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

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



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

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

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



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


Log in to reply