vof法求解自由表面中遇到的一个问题，求高手指点！

jcn

我用的是有限元和vof求解自由表面问题。应该说采用什么方法求解流场和采用什么方法求解自由表面没有什么关系。毕竟现在的vof法都是基于几何的方法求解的，也就是通过几何的方法求解下一时刻的单元的体积分数。
但是，从我现在的情况来看，用vof求解自由表面容易产生“气泡”现象，即单元过空或单元过满。从误差分析来看，这种现象主要是由于体积通亮不守恒引起的。但是，同样的程序，在没有自由表面的情况下（比如求解圆柱绕流问题），结果非常好，涡街的形状和脱落周期都还算完美。但是一同vof一起求解就会出现以上问题。我的vof程序也是经过验证的（在流畅速度已知，且严格守恒情况下），精度较高。现在还没有解决这个问题，希望高手能够指点一二。