rhoCentralFoam 为什么做成半隐式？

[sibo]

大家好，

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

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

欢迎大家讨论，谢谢！

[sibo]

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

[东岳]

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

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

[sibo]

@东岳 感谢前辈回复！

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

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

[东岳]

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

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

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

[sibo]

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