关于自由表面流动的压力

twophase

我在用vof,fvm,roe scheme,preconditioning and coupled solver开发有自由表面的水和空气的流动。
假设u v速度都是0，那么v 动量方程就只有
dp/dy=-rou*g.
假如我在压力初场给定的是静水压力，那么方程可以写成
-d(rou)/dy*y-rou*g=-rou*g.
这意味着只要有自由表面存在，在自由表面附近的网格里面在v方向有剩下的量。
我看fluent也有一些说明但是不是很清楚怎么处理，有没有高手指点一二，不剩感谢！！