VOF模型介绍要点

2数值计算模型和方法2.1基本理论实验测量、理论分析和数值模拟是研究流体常用的三种方法。三种方法都有各自的优缺点。实验测量的结果具有直观性和正确性；理论研究具有普遍性，能够为其他的研究方法提供理论支持；而数值模拟是最先进的不受设备条件要求的研究方法。目前数值模拟的方法因其方便且受限制小而被大家广泛应用，但是数值结果也需要与实验数据进行对比，通过实验来验证模型的正确可行性。正确的模拟方法可以为我们研究的对象提供实验所观察不到的有价值的研究和预测。三种方法在研究中是相辅相成的。目前数值计算的研究方法已经有了长足的发展，得到了广泛的应用。此种方法将需要研究的问题在虚拟的环境中进行计算，排除了时间、空间等各种因素的限制，我们可以对不同参数和工况进行取值，得到以下在实验中无法观察得到的数据，因此具有实验无法达到的优越性。在模拟计算中，只要建立合理的数学模型和物理模型，就能借助计算机就可以得到理想的结果。对于流体换热或流动问题，都可以采用CFD方法解决，CFD方法计算过程如下：选择合理的数学模型来描述简化之后的实际问题，只有建立其合理的数学模型，才能对问题进行分析。选择合理的计算方法，建立描述数学模型的方程，选择合理的方法对计算区域进行合理的处理。同时选择合理的求解方法、设置边界条件、微分方程的离散化方法以及坐标系的建立等等也是至关重要的。直接关系着结果的正确与否。利用Gambit、ICEM等软件对计算区域进行网格划分。在FLUENT中对区域进行定义和参数的设定，这是模型改变最多的部分，也是CFD中最重要的部分。得到计算结果。可以利用FLUENT自带的后处理软件对数据进行处理，结果可以通过图表或者散点图等形式表现出来，通过给出的结果分析各项参数的变化情况。2.2VOF模型（VolumeofFluentModel）VOF模型，是建立在固定的欧拉网格下的表面跟踪办法。建立在两种或者多种流体（或相）不相互混合的前提下。当需要得到一种或者多种互不相融的流体交界面时，可以采用这种模型。在VOF模型中，不同的流体组分共用、着一套动量方程，通过引进相体积分数这一变量，实现对每一个计算单元相界面的追踪。在每个控制容积内，所有相体积分数额总和为1，。所有变量及其属性正在控制容积内各相共享，并且代表了容积平均值。这样，在任何给定控制容积内的变量及其属性纯粹的代表了一相或者相的混合，并且由相体积分数决定。换句话说，在单元中，若第q相流体体积分数为，那么可能存在以下三种情况：=0：单元里不存在第q相流体。=1：单元里充满了第q相流体。0<<1:单元里包含了第q相流体和一相或者其他多相流体的界面。基于的局部值，适当的属性和变量在一定范围内分配给每一个控制单元。2.1体积分数方程（连续性方程）跟踪相之间的界面是他通过求解一相或者多相的容积比率的连续方程来完成的。对第q相，有（2-1）其中：为质量源项，在默认情形下方程2-1右端源项为零，但当你给每一相指定常数或用户定义的质量源，则右端不为零。主相的体积分数的计算基于如下约束：（2-2）2.2属性计算（2-10）其中：和是两种流体界面两侧的压力。在Fluent中，使用CSF公式时，这里的表面曲率hi从垂直于界面的表面局部梯度计算的。n为表面法线，定义为第q相的体积分数的梯度，则有：（2-11）表面曲率是为了区别单位法向量而定义的：（2-12）其中：（2-13）表面张力也可以根据越过表面的压力变化得出。表面力使用散度定力可以表示为体积力，它有如下形式：（2-14）若一个单元中只有两相，那么有；。因此方程2-14可简化为：（2-15）其中ρ是使用方程2-14计算的平均密度。方程2-15说明一个单元表面张力源项是和单元平均密度成比例。2）壁面粘附假定流体与壁面产生的接触角常用于调节壁面附近单位表面的法向而不是加强壁面本身的边界条件，这个所谓的动力壁面的边界条件熬制了壁面附近表面曲率的改变。若是壁面的接触角，那么挨着壁面的实际单元表面法向为：（2-18）这里和分别是壁面额单位法向量和切向量这个接触角与一个单元正常计算的表面法向量联合决定了表面的局部曲率，这个曲率常用于调整表面张力计算中的体积力项。接触角壁面和壁面上界面切线的夹角如图2-2所示。图2-2示意图2.3湍流模型图2-3湍流模拟分类湍流模型在数值模拟计算中至关重要，模型的选择直接影响计算结果的正确性。具体模型分类详见图2-3。直接模拟就是直接对N-S方程进行求解，理论上能得到精确值解，但计算难度高不容易实现。非直接模拟是对某部分进行简化处理，从而得到近似解。实际应用中，比较常用的有标准k-ε模型和RNGk-ε模型。本文采用了标准k-ε模型，下面将对标准k-ε进行具体阐述：标准k-ε模型由湍动能k和湍流耗散率ε两个方程组成，其中湍动粘度为k和ε的函数 （2-19）标准的k-ε模型的方程可以表示为 （2-20） （2-21）式中，是速度梯度引起的湍能项，是由浮力引起的湍能项。 （2-22）对可压缩流体，有式中，Pr是湍动普朗特数。是膨胀系数，。为脉动扩张项。对可压缩流体，。其中为湍动马赫数，。当流体不可压缩时，，。当忽略源项时，，，上述方程简化为： （2-23） （2-24）