如果不使用staggered grid，速度存储在网格中心，怎么模拟不可压流体

[xiangzi]

问题

模拟不可压缩的流体，如果速度存储在网格中心， 那么project以后，怎么修改网格的速度呢？

[houkensjtu]

相对于交差网格，将速度和压力都保存在网格中心的形式称为同位网格（collocated grid）。为了保证在同位网格下不出现阶梯状的波状压力场，需要在动量方程的压力梯度项里引入相邻两点的压力差（如果按照交错网格的方法去近似压力梯度，压力项会缺少当前网格点本身的压力，只有周围压力的值，从而导致波状压力出现）。
具体的做法可以参考一些文献，比如提出动量插值法的Rhie & Chow的论文"A numerical study of the turbulent flow past an isolated airfoil with trailing edge separation"，知乎上的这篇回答也很详细。