[求救]水气两相流的压力出口的udf？

表面的和平

我是用VOF模拟二维水气两相流，出口上部是空气（主相），下部是水但是出口需要给压力设置udf：
水的压力：P＝Po+ρgh,
我的想法是利用y坐标与水深比较，高于水深的压力为0pa，低于水深的压力P＝Po+ρgh，
但是因为水深是随着时间，从0增加到某一未知水深，所以不知道用什么函数来求？
本来想用求副相体积分数的宏 CVOF(c,pt[1])来算边界上的水深，但是帮助里说它只能求出一个单元上水的体积分数。
恳求高人帮我指点～  先谢谢了！

pvxd

实在要算，可以试试计算和界相邻的单元cell的second phase体积分数。
不过觉得这个判断方法不是很好。
另外，不是很明白你的问题的物理过程，呵呵