沃纳·海森堡曾说:“当我见到上帝后,我一定要问他两个问题——什么是相对论,什么是湍流。但我相信他应该只对第一个问题有了答案。”
流体力学中为便于研究,对流体特性引入假设,而将流动分成两大类:无粘流动(有势流动,简称“势流”)和粘性流动(简称“粘流”,包括层流和湍流)。目前势流计算相对成熟,它具有简易、计算量小等特点而成为切实的工程方法。就CFD技术本身功能发展而言,目前主要是处理粘性流动问题,特别是湍流,理论上来说,对于任何空间点,湍流中流体运动都满足非线性N-S方程,所以至少可以用数值计算机做数值运算,然而这一工作量是十分巨大的。数值求解这一方程的方法有三大类:直接数值模拟方法(DNS),大尺度涡动模拟方法(LES),雷诺平均方法(RANS),由于计算机硬件资源的限制,前二种方法暂还没有得到广泛的应用,而雷诺平均方法在工程上得到了大量的应用。而雷诺平均方法(RANS)是一个湍流模型大家族,包含了很多的子模型,因而本文将主要对RANS模型下的一些子模型进行简要介绍。
控制方程
在介绍常用的湍流模型之前,先把描述流体流动的控制方程进行罗列。
连续方程
雷诺平均方法(RANS)
在雷诺平均法中,把瞬时速度u分解为时均速度 和脉动速度 ,则对于速度分量:
同样,对于压力和能量、相分数或其他标量也采用这种方式进行分解。
设及是两个瞬时值,及是相应的脉动值,则有以下关系成立:
对于连续性方程,将三个坐标方向的瞬时速度表示为时均速度和脉动速度之和并代入,再对该式作时均运算,得到:
这表明湍流速度的时均值和脉动值仍然满足连续性方程。
对于动量方程,以x方向为例,作类似于上面的处理可得
利用上面给出的关系式可得
将左端脉动分量乘积的时均值放到右端粘性项中,即
对其他两个方向也可以做类似的推导,对不可压缩流体,常写作张量符号形式,可得到下列时均形式的N-S方程,即Reynolds方程
上式的雷诺应力通常记为。写作分量形式如下
是对角对称的,包含六个未知量:三个正应力分量(对角)和三个切应力分量(非对角)。
要求解 RANS 方程,直观的做法是求解各应力分量的对流输运方程。对应力分量的对流输运方程进行雷诺平均,会得到更高阶的未知量 。实际上这样下去会得到无限多的未知量,整个方程组永远无法封闭。这就需要采用一些近似的模型来封闭方程组,比如用某些平均量来近似表示雷诺应力。这些不同的近似方法则被统称为湍流模型(turbulence models)。所谓湍流模型就是把湍流的脉动值附加项与时均值联系起来的一些特定关系式。
从数学形式上看,雷诺应力很像二阶应力张量。因此 Boussinesq 仿照分子粘性应力与速度变形率的关系,提出了关于涡粘性(eddy viscosity)的 Boussinesq 假设,认为雷诺应力与平均速度变形率成线性关系:
基于以上假设的模型被称为涡粘模型(Eddy Viscosity Model, EVM)。Boussinesq 假设被用于Spalart-Allmaras单方程模型和双方程模型,其优点是与求解湍流粘性系数有关的计算花费时间较少,但缺点是认为湍流粘性系数是各向同性标量,对于一些复杂的流动条件下并不严格成立,所以具有其应用限制性。工业界常用的涡粘模型有一方程 Spalart-Allmaras 模型,两方程 k−ε 及 k−ω 模型等。
如果湍流场各向异性很明显,如强旋流动以及应力驱动的二次流等流动中,这时可采用雷诺应力模型(Reynolds Stress Model, RSM)在更高阶的应力上封闭 RANS 方程。具体做法通常是求解六个未知应力分量的方程及其对应的耗散率方程。雷诺应力模型可以较好的模拟具有强烈非线性特征的复杂流动,然而由于它需要求解许多额外的方程,因此计算量相比涡粘模型大得多,对计算机内存也有更高的要求,因此在工程中应用不多。常见的雷诺应力模型有 LRR (Launder-Reece-Rodi)和 SSG(Speziale-Sarkar-Gatski) 模型。
零方程模型
所谓零方程模型是指确定湍流黏性系数不需要微分方程的模型,最简单的零方程模型是常系数模型,对自由剪切层流动,Prandtl提出在同一截面上为常数。
在二维坐标系中,湍流切应力表示成为
混合长度理论适用于一些比较简单的流动,如边界层类型流动;平直通道内的流动;回流较弱接近于边界层类型流动等。但是混合长度理论在物理概念上存在一些不足如,在管道中心线处速度梯度为零,但实际不为零;不考虑湍流的历史(上游情况)的影响;不考虑湍流强度的影响。因此需要引入偏微分方程来确定湍流黏性系数来适应这方面的需要。
一方程模型
以Spalart-Allmaras单方程模型为例,它具有良好的鲁棒性和数值收敛性,可以很好的模拟绝大部分的附着流动和薄层自由剪切流动。Spalart-Allmaras单方程模型基于线性涡粘假设,湍流粘度按以下公式计算
0.66666 | 0.41 | 0.1355 | 0.622 | 0.3 | 2.0 | 7.1 | 0.3 |
二方程模型
k−ε模型
在一方程模型中,湍流长度的标尺是由经验公式给出的。文献中广泛采用形如的公式来选择与湍流脉动的长度标尺有关的量。己经采用的变量的主要形式列于表中。Rodi指出,所有这些变量的微分方程形式均类似,但对靠近壁面地区来说方程计算最为方便。因而,在湍流的工程计算中,两方程模型使用最宽广。
变量 | ||||
---|---|---|---|---|
提出者 | Kolmogorov | Chou(周培源) | Rodi, Spalding | Spalding |
符号 | ||||
物理意义 | 涡旋频率 | 能量的耗散 | 能量与标尺之积 | 涡量脉动的时均方值 |
耗散率的严格定义:
耗散率的模拟定义:
标准模型需要求解湍动能及其耗散率方程。湍动能输运方程是通过精确的方程推导得到,但耗散率方程是通过物理推理,数学上模拟相似原形方程得到的。该模型假设流动为完全湍流,分子粘性的影响可以忽略。因此,标准模型只适合完全湍流的流动过程模拟。
标准模型的湍动能及其耗散率方程为如下形式
在上述方程中,表示由于平均速度梯度引起的湍动能产生,是用于浮力影响引 起的湍动能产生,可压速湍流脉动膨胀对总的耗散率的影响。湍流粘性系数
在FLUENT中,作为默认值常数,,,,湍动能与耗散率的湍流普朗特数分别为,可以通过调节“粘性模型”面板来调节这些常数值。
k−ω模型
湍动能和比耗散率由以下输运方程得出:
其中表示平均速度梯度产生的湍流动能;表示w的生成;和 分别代表和的有效扩散系数;和表示和在湍流作用下的耗散;和是用户定义的源项。
模型的有效扩散系数由下式给出:
和分别代表和的湍流普朗特数,结合和计算湍流粘度:
表示平均速度梯度产生的湍流动能。由的输运方程可知,这一项可以定义为:
的产生由下式给出:
的耗散为:
的耗散由下式给出:
可压缩函数由下式给出:
压缩效应在非常有限的自由剪切流实验中进行了校准,不推荐普遍使用。默认情况下是禁用的。
模型常数如下:
参考链接
Evaluation of Two-Equation Turbulence Models for Predicting Transitional Flows