为了能够简单高效地处理两相流动问题,Hirt和Nichols于1981年开创性地提出了VOF(volume-of-fluid)方法,并引入了流体体积函数C的概念,流体体积函数C等于一个单元内目标流体体积与该单元体积之比。
C=1表示单元内充满目标流体,C=0表示单元内不含目标流体,而0
连续性方程和动量方程中的密度和粘度则通过对两相所占的比例进行加权平均得到。
显然,与单相流动相比,两相流动中多了一个变量C,那么根据多一个变量就需要多一个方程的原则,还需要列出一个关于C的方程。不过还好,这个方程是比较容易找的。
因为根据质量守恒的原则,网格内C的增加量就等于流入网格的C(如果存在相变,还需要考虑源项,这里先不予讨论),即