CFD中文网

openfoam里面的rhoPimpleFOAM采用的PISO算法能适用于可压缩问题，但一直受诟病。最近在看一下投影法的文章，可以很方便实现RK4格式，但投影法是基于不可压缩的，能否推广可压缩问题，显式RK格式有没有植入到可压缩PISO的文章。

@李东岳 东岳老师，您有什么好的解决思路吗？

方程如图，等号最右端那一项是介质中气泡膨胀收缩的体积分数变化源项，会给方程引入比较大的压力变化，这种源项有什么好的处理方法。

@李东岳 感谢东岳老师解答

@李东岳 感谢东岳老师，这个审稿人应该是做高精度算法的，给的拒稿意见，对分离式求解，一阶时间精度意见比较大，比较愁人，另外他还觉得二阶TVD格式计算对流项更好，我使用的NVD gamma格式他也表示不理解，当时算的时候没太管这个，我理解的带限制的TVD/NVD格式不是差不多吗，精度介于一阶迎风和中心差分之间，不知道怎么回复了。

@李东岳 东岳老师，能不能帮忙解答一下，比较困惑，引的那两篇也是OpenFOAM做爆炸气泡比较经典的论文。

审稿人提出意见对OpenFOAM中常用的一阶时间精度与PIMPLE算法模拟水下爆炸质疑，具体如下：

1. Line 52 reads "Since the calculation of pressure wave propagation in a single medium is fully developed (Koch et al., 2016; Li et al., 2019; Miller et al., 2013; Wu et al., 2020)". Modelling pressure wave propagation remains a significant challenge for the semi-implicit methods cited here. Especially the work of Koch et al. and Miller et al. are very questionable works in this respect, as they use first-order time differencing which is arguably a poor choice if pressure waves are to be approximated accurately.
   2.Generally, the PIMPLE algorithm is widely used for pressure-velocity decoupling in incompressible flows. However, the subject of this study appears to be compressible multiphase flow. What is the purpose of using PIMPLE for pressure iteration correction in this context?
   如何回复为好。