rhoCentralFoam 为什么做成半隐式？
-
大家好，
rhoCentralFoam 在处理inviscid flow时完全显式，对于粘性项做隐式处理。
好奇为什么不直接做成 fully implicit? 因为作者在文中也提到考虑到步长问题，将粘性项做了隐式处理。还是说 central-upwind 的 KT/KNP 在全隐式时不可行？
欢迎大家讨论，谢谢！
-
这里附上@东岳 前辈的 rhoCentralFoam解析
以及rhoCentralFoam作者的 文章
-
一方面对于存在激波这种双曲欧拉方程，显性离散虽然耦合性不强但对于间断的分辨率更好。对于有粘流动，本身激波要被抹掉，不需要捕获间断，隐形格式更加耦合
另一方面KT算法将压力项和对流项放在一起，是要隐性求解速度呢还是压力呢？
-
-
我见过有文献做成piso+central/upwind格式的sci
rhoCentralFoam里面隐式的求解应该对速度修正吧？加入粘性项对速度的影响。
没有修正，并没有采用迭代的方式处理压力速度耦合
-
@东岳 是啊我也看到了，改成压力基说是在low mach时更好。感觉各种拼凑。。
