请教关于基于压力的求解方法用于可压缩问题时，公式如何推导

[Jacobian]

像PISO算法，一般是用于求解不可压问题的。但经过一定的改造，是可以求解可压问题的，比如openfoam中的rhoPimpleFoam。
但我现在不是很明确速度预测方程，以及两次修正步的压力方程应该怎么写。大家能否为我解惑，或者提供一下这方面的权威文献？

[李东岳]

这个不是在icoFoam里面已经有了标准代码了么？
要看PISO，找找关键词R I Issa 1987 PISO

[Jacobian]

@李东岳 感谢东岳老师的回复
但我还是不太明白。我再把疑问写得具体一点，请您提出批评指正

按照Issa所述，可压时的速度预测方程为：

可以整理为：

openfoam中通常令：

而经过第一次修正后的速度，密度应当满足：

整理为：

右侧的相邻单元修正未知，u_n两星 用预测步得到的 u_n星 代替：

现在眼看着就能将它代入连续性方程：

但还有一个问题：

该项中还存在着一个未知的密度rho星，以至于和预测步公式的形式上产生了不统一，不能凑出HbyA了，

Issa 在 Solution of the Implicitly Discretised Fluid Flow Equations by Operator-Splitting中的公式（46）似乎直接用n时间步的rho代替rho星了？还是我上面的推导思路不对？

[李东岳]

目前我在出差，涉及到公式的以及翻文献的还讨论不了，很抱歉！等我出差回北京我看看。
希望论坛其他人有思路。

[Jacobian]

@李东岳 好的，谢谢！