v2015两数值模拟第9讲vof方法及其应用.pptx

第九讲：VOF ( Volume-Of-Fluid )方法及其应用;　;气(液)－液两相流的相界面描述方法可分为两大体系：即 　　（１）界面捕捉法 　　　　　(Front-Capturing Methods) 　　（２）界面跟踪法 　　　　　(Front-Tracking Methods)。;　　界面追踪方法，直观的理解，属于Lagarange类方法。 　　如果有足够多双眼睛盯住／“锁定”界面上的各个点，通过对各个点的观察、记录就可以确定整个界面的运动过程和轨迹。这就是界面追踪方法产生的最初的出发点和本质。 　　 这类方法中最具代表性的是： 　　　　PIC (Particle-in-Cell) 　　　　MAC (Mark-and-Cell);界面捕捉法中最具有代表性的是： １ Level Set ２　VOF ３　高度函数法　　 　;１．1　VOF方法的基本思想： 　　 1981年，C. W. Hirt 和B. D. Nichols在 International Journal of Computational Physics 杂志上首先正式发表了著名的VOF论文，开创性地提出了用VOF方法进??运动相界面追踪的思想，并首先用VOF方法对溃坝和浪涌（Broken Dam, Breaking bore）及Rayleigh Taylor不稳定性现象进行了成功的数值模拟模拟。 ;VOF方法的基本思想 VOF方法在整个流场定义了一个相函数f，在每个网格中，这个函数的定义为一种流体（目标流体）的体积与网格体积的比值。相函数满足一定的输运方程，在得出相函数输运方程之前先定义一个染色函数 从而在控制体σ上，目标流体的体积函数表达式为 由拉格朗日流体体积的特性可知 ;将上述全导数展开可得输运方程为 相函数f的定义式为 因此可得相函数满足的输运方程为 当f=0时，网格内没有目标流体，称为空网格；当f=1时，网格被目标流体充满，称为满网格；当0f1时，网格含有目标流体，但未被充满，称之为界面网格。;（１）VOF方法用相函数(Phase Function)F取代了MAC方法中的虚拟无质量彩色粒子。 　　　从这个意义上说，VOF方法可看作是MAC方法的一个变种。 （２）在一种流体相（比如说，液体）中，相函数F取值为1，而在另一种流体（比如，气体或另一种液体）中取值为0；在相函数取0到1之间的数值的地方即为相界面位置。 （３）相界面的取向可由界面附近各点上的F值来确定。; 　　 （4） VOF方法避开了采用工程浩大的Marker点的方案，比MAC方法对计算机内存的要求低，更容易实施。 （5） Hirt和Nichols在最初的VOF模型中，设计了类似于MAC方法和PIC方法的Staggered型差分格式——即将压力和相函数（体积分数）定义在格子的中心处，而将X方向的速度定义在格子的左、右格边中点，将Y方向的速度定义在格子的上、下格边中点。;１．２　关于相函数的概念： （１）相函数类似于气－液两相流中的截面含气率（容积含气率），表示某一相介质占据网格面积（二维）或体积（三维）的分数。 （２）相函数是以一个网格为单元来定义的，与含气率的概念不同。;１．２　关于相函数的概念： （３）相函数就是一个介质指针，指示着占据某一望个的介质种类。但这种指示不仅是定性的，而且是定量的。 （４）对应于同一个相函数值，气－液界面在网格内的形状和方位是多值的。如图所示。 （５）F必须而且只能在0和1之间。;VOF方法中的基本问题： 1.　气体的动量／运动控制方程 2.　液体的动量／运动控制方程 3.　相函数F的控制方程及其求解 4.　如何由相函数F的分布获得气－液相界面——相界面的构造问题。;三、气－液两相流的动量方程　　——两相流运动控制方程;三、气－液两相流的动量方程　　——两相流运动控制方程;三、气－液两相流的动量方程　　——两相流运动控制方程;三、气－液两相流的动量方程　　——两相流运动控制方程;四．气－液相界面的控制方程;四．气－液相界面的控制方程;五．气－液相界面控制方程的求解方法;五．气－液相界面控制方程的求解方法;五．气－液相界面控制方程的求解方法;六．气－液相界面的构造方法;六．气－液相界面的构造方法;六．气－液相界面的构造方法;六．气－液相界面的构造方法;六．气－液相界面的构造方法;六．气－液相界面的构造方法;六．气－液相界面的构造方法;六．气－液相界面的构造方法;六．气－液相界面的构造方法;六．气－液相界面的构造方法;六．气－液相界面的构造方法;六．气－液相界面的构造方法;六．气－液相界面的构造方法;方程（４）中的最小函数“Min｛s｝”主要是为了防止从“施主”网格中过多的流出流体;——在一个时间单元内，流过（流出／流入）