[求助]有关simple中重力方程的处理

lyzhang98

本人正在用simple方法模拟一个三维流动，遇到困难，不知道如何处理，请高手指点。
在三维方程中，有一维的速度是重力方向上的，这的速度的微分方程中必须有重力项（体积力），或者单独的列一个重力方程。我采用了用重力代替压力的影响，就是将压力离散项：dv(i,j,k)*(p(i,j,k-1)-p(i,j,k))用：densit*9.8*vol
这样的话，压力方程的重力方向的压力相当于没有用，而且程序也开始出错.
请您帮忙分析一下，在下感激不尽！！