VOF方法中相函数的控制方程的求解方法..doc

VOF方法中相函数的控制方程的求解方法 一、VOF方法的基本原理 VOF方法的基本原理是通过研究网格单元中流体和网格体积比函数F来构造和追踪自由界面的，而非追踪自由液面上质点的运动。VOF方法是在整个流场中定义一个函数F，在每个网格中，这个函数定义为一种流体（目标流体）的体积与网格体积的比值。只要求得每个网格上的F值，就可以构造出运动界面，然后求解物理方程时，可以在界面附近作特殊的精细处理，提高分辨率和精度。由于每个网格只需保存一个F值，从而可以大大减少储存量。 设计算区域为，流体A所在区域记为，流体B所在的区域记为，首先定义一个函数， (1) 对于两种不相容流体组成的流场，满足拉格朗日流体体积的特性： (2) 其中是流体的速度场，为单个网格的体积，在每个网格上定义为在网格上的积分， (3) 称之为VOF函数，同样它满足： (4) 称之为VOF(Volume of Fluid)方程。 由上可以看出每个流体上的体积函数实质上是 显然，F=1的网格充满流体A，我们称之为流体网格(fluid)；当F=0时为空网格(void)；而0F1的网格，则是含有流体界面的网格，称为界面网格(surface)。 VOF方法是通过求解流体体积函数，实现对运动界面的追踪。由于流体体积函数是单个网格中一种流体体积与网格体积之比，是具有物理意义的。但在求解VOF方程时，要保持它的这一物理意义是困难的。首先要求求解格式是守恒的，而且每一时间步得到的函数值在跨过界面时，即在含有界面的网格上，VOF函数要0F1，在其它网格上必须是1或0。 二、流体体积方程的求解方法 Hirt和Nichols针对界面追踪的一般方法以及通过装配格子内界面的斜率来确定运动界面的斜率和位置，提出了VOF方法。在这以后，AShgriz和Poo, Kim等人通过对格子内界面斜率更加精确的构造得到了其它不同的重构方法，尽管他们的精确性和形式是不同的。Rider和Kothe, Gueyffier和Lijie, Youngs等人利用了运动界面的法向矢量近似值来构造界面。Rudman, Vincent和Caltagirone, Ubbink和Issa等人采用了更精细的方法求解流体体积函数方程。对于流体体积方程的求解归纳起来主要有两种方法：一种是基于几何学的求解方法，这种方法又可以成为重构技术；另一种则是基于数学的求解方法。 1. Hirt Nichols法： 图1 Hirt-Nichols的重构模板 这是一种基于几何学的求解方法，该方法用水平和竖直两种直线近似一个网格内的流体界面。采用施主--受主(Donor-Acceptor)的差分逼近格式计算运动界面的流体体积函数方程。 界面的重构如下：将流体的自由表面看作局部的单值函Y(x)和X(x)i, i+1网格的列Yl值和j-1, j, j+1网格列的Xl值，估算出每个网格上界面的斜率值dY/dx和dX/dy，然后根据流体体积函数和斜率的大小确定网格(i, j)上的自由面的位置和方向。 (5) (6) (7) (8) 如果：； 则定义网格内的界面是水平的；否则，网格内的界面是竖直的。 2. FLAIR法： 这种求解流体体积方程的方法也是基于几何学理论。1991年，Ashgriz和Poo的FLAIR(Flux Line-segment model for Advection and Interface Reconstruction)方法通过构造跨过任意网格边界的两个相邻网格的带有倾角的直线段来近似界面。由于界面构造涉及到两个网格，它分成多种情况分别进行计算。首先利用网格边界的流体速度，确定施主单元和受主单元，分别用FD, FA表示其流体体积函数，然后按以下几种情况计算： （1）施主网格是满网格，即FD=1.0; （2）施主网格是空网格，即FD=0.0; （3）施主网格和受主网格都是半网格，即0.0FD1.0，0.0FA1.0; （4）施主网格是半网格，受主网格为空或者满网格，即0.0FD1.0，FA=0.0，或FA=1.0; 对于前两种情况，单位时间内，流体体积流通量分别为，0.0。第三种情况又可分为十六中情形，经过对称和反转可以分为四类(如图2所示)。此时，且目标流体为图中阴影部分。近似界面直线的斜率可以通过两边的流体体积函数值来计算。从而可以计算出网格边界的流量运输，即的长度上阴影部分的面积(左边为施主网格)。 图2 FLAIR的界面重构 3. Youngs法： Youngs的PLIC (Piecewise Liner Interface Calculation) 方法，这种方法也是基于几何学原理的，