湍流模型概述
- 格式:doc
- 大小:653.50 KB
- 文档页数:17
第三章湍流模型湍流模型第一节前言湍流流动模型很多,但大致可以归纳为以下三类:第一类是湍流输运系数模型,是Boussinesq于1877年针对二维流动提出的,将速度脉动的二阶关联量表示成平均速度梯度与湍流粘性系数的乘积。
即:3-1 推广到三维问题,若用笛卡儿张量(笛卡尔坐标系)表示,即有:3-2 为DELT函数,一般i=j时为1,否则为0. 模型的任务就是给出计算湍流粘性系数的方法。
根据建立模型所需要的微分方程的数目,可以分为零方程模型(代数方程模型),单方程模型和双方程模型。
(模拟大空间建筑空气流动)μt=0.038 74ρvl (模拟通风空调室内的空气流动)比例系数由直接数值模拟的结果拟合而得,其中:v为当地时均速度,l为当地距壁面最近的距离。
第二类是抛弃了湍流输运系数的概念,直接建立湍流应力和其它二阶关联量的输运方程。
第三类是大涡模拟。
前两类是以湍流的统计结构为基础,对所有涡旋进行统计平均。
大涡模拟把湍流分成大尺度湍流和小尺度湍流,通过求解三维经过修正的Navier-Stokes方程,得到大涡旋的运动特性,而对小涡旋运动还采用上述的模型。
实际求解中,选用什么模型要根据具体问题的特点来决定。
选择的一般原则是精度要高,应用简单,节省计算时间,同时也具有通用性。
参见:湍流模型的选择资料。
FLUENT提供的湍流模型包括:单方程(Spalart-Allmaras)模型、双方程模型(标准κ-ε模型、重整化群κ-ε模型、可实现(Realizable)κ-ε模型)及雷诺应力模型和大涡模拟。
Zero-Equation Models One-Equation Models Spalart-Allmaras Two-Equation Models Standard k-e RNG k-e Realizable k-e Reynolds-Stress Model Large-Eddy Simulation Direct Numerical Simulation 包含更多物理机理每次迭代计算量增加FLUENT提供的模型选择RANS-based models 湍流模型种类示意图大涡模拟启动需要用命令:(rpsetvar 'les-2d? #t) 第二节平均量输运方程输运过程的粘滞系数、扩散系数和热传导率,故称为输运方程雷诺平均就是把Navier-Stokes方程中的瞬时变量分解成平均量和脉动量两部分。
9个湍流模型介绍
好的,为你介绍9个湍流模型:
1. Reynolds平均的NS方程(Reynolds-Averaged Navier-Stokes,RANS):Reynolds 提出了平均法,将“瞬时值=平均值+脉动值”带入不可压缩流体控制方程中,得到了一个更复杂的方程。
对于可压缩流体,假设瞬时密度的变化对流动影响不大,忽略其影响。
2. Reynolds应力模型(RSM):模仿控制方程的样子,搞出一个针对Reynolds应力的输运方程。
3. 代数应力模型(ASM):简化Reynolds应力方程的对流项和扩散项。
此外,还有一些其他湍流模型,如Spalart-Allmaras模型、k-双方程模型等。
这些模型都有各自的特点和适用范围,可根据具体问题选择合适的湍流模型进行计算。
湍流模型目前计算流体力学常用的湍流的数值模拟方法主要有以下三种:直接模拟(direct numerical simulation, DNS)直接数值模拟(DNS)特点在湍流尺度下的网格尺寸内不引入任何封闭模型的前提下对Navier-Stokes方程直接求解。
这种方法能对湍流流动中最小尺度涡进行求解,要对高度复杂的湍流运动进行直接的数值计算,必须采用很小的时间与空间步长,才能分辨出湍流中详细的空间结构及变化剧烈的时间特性。
基于这个原因,DNS目前仅限于相对低的雷诺数中湍流流动模型。
另外,利用DNS模型对湍流运动进行直接的数值模拟对计算工具有很高的要求,计算机的内存及计算速度要非常的高,目前DNS模型还无法应用于工程数值计算,还不能解决工程实际问题。
大涡模拟(large eddy simulation, LES)大涡模拟(LES)是基于网格尺度封闭模型及对大尺度涡进行直接求解N-S方程,其网格尺度比湍流尺度大,可以模拟湍流发展过程的一些细节,但其计算量仍很大,也仅用于比较简单的剪切流运动及管流。
大涡模拟的基础是:湍流的脉动与混合主要是由大尺度的涡造成的,大尺度涡是高度的非各向同性,而且随流动的情形而异。
大尺度的涡通过相互作用把能量传递给小尺度的涡,而小尺度的涡旋主要起到耗散能量的作用,几乎是各向同性的。
这些对涡旋的认识基础就导致了大涡模拟方法的产生。
Les大涡模拟采用非稳态的N-S方程直接模拟大尺度涡,但不计算小尺度涡,小涡对大涡的影响通过近似的模拟来考虑,这种影响称为亚格子Reynolds应力模型。
大多数亚格子Reynolds模型都是将湍流脉动所造成的影响用一个湍流粘性系数,既粘涡性来描述。
LES对计算机的容量和CPU的要求虽然仍然很高,但是远远低于DNS方法对计算机的要求,因而近年来的研究与应用日趋广泛。
应用Reynolds时均方程(Reynolds-averaging equations)的模拟方法许多流体力学的研究和数值模拟的结果表明,可用于工程上现实可行的湍流模拟方法仍然是基于求解Reynolds时均方程及关联量输运方程的湍流模拟方法,即湍流的统观模拟方法。
第三章 湍流模型第一节 前言湍流流动模型很多,但大致可以归纳为以下三类:第一类是湍流输运系数模型,是Boussinesq 于1877年针对二维流动提出的,将速度脉动的二阶关联量表示成平均速度梯度与湍流粘性系数的乘积。
即:2121x u u u t ∂∂=''-μρ 3-1 推广到三维问题,若用笛卡儿张量表示,即有:ij ijj i t j i k x u xu u u δρμρ32-⎪⎪⎭⎫ ⎝⎛∂∂+∂∂=''- 3-2 模型的任务就是给出计算湍流粘性系数t μ的方法。
根据建立模型所需要的微分方程的数目,可以分为零方程模型(代数方程模型),单方程模型和双方程模型。
第二类是抛弃了湍流输运系数的概念,直接建立湍流应力和其它二阶关联量的输运方程。
第三类是大涡模拟。
前两类是以湍流的统计结构为基础,对所有涡旋进行统计平均。
大涡模拟把湍流分成大尺度湍流和小尺度湍流,通过求解三维经过修正的Navier-Stokes 方程,得到大涡旋的运动特性,而对小涡旋运动还采用上述的模型。
实际求解中,选用什么模型要根据具体问题的特点来决定。
选择的一般原则是精度要高,应用简单,节省计算时间,同时也具有通用性。
FLUENT 提供的湍流模型包括:单方程(Spalart-Allmaras )模型、双方程模型(标准κ-ε模型、重整化群κ-ε模型、可实现(Realizable)κ-ε模型)及雷诺应力模型和大涡模拟。
湍流模型种类示意图Direct Numerical Simulation包含更多 物理机理每次迭代 计算量增加提的模型选RANS-based models第二节 平均量输运方程雷诺平均就是把Navier-Stokes 方程中的瞬时变量分解成平均量和脉动量两部分。
对于速度,有:i i i u u u '+= 3-3其中,i u 和i u '分别是平均速度和脉动速度(i=1,2,3)类似地,对于压力等其它标量,我们也有:φφφ'+= 3-4 其中,φ表示标量,如压力、能量、组分浓度等。
湍流模型介绍
湍流模型是数学模型的一种,用于描述液体或气体中的湍流运动。
湍流是一种不规律的、难以预测的流体运动,通常是由于速度、密度或温度的不规则分布引起的。
湍流模型通过使用一系列方程,描述流体的速度、压力和密度等参数之间的相互作用,以预测和模拟流体的复杂运动行为。
湍流模型主要分为两类:基于雷诺平均的模型(如k-ε模型、k-ω模型)和直接数值模拟(DNS)。
每种模型都有其适用的范围和局限性,需要根据具体问题的特性选择合适的模型。
湍流模型在气象、水文、工程、航空航天等领域中得到了广泛应用。
湍流模型介绍因为湍流现象是高度复杂的,所以至今还没有一种方法能够全面、准确地对所有流动问题中的湍流现象进行模拟。
在涉及湍流的计算中,都要对湍流模型的模拟能力以及计算所需系统资源进行综合考虑后,再选择合适的湍流模型进行模拟。
FLUENT 中采用的湍流模拟方法包括Spalart-Allmaras模型、standard(标准)k −ε模型、RNG(重整化群)k −ε模型、Realizable(现实)k −ε模型、v2 −f 模型、RSM(Reynolds Stress Model,雷诺应力模型)模型和LES(Large Eddy Simulation,大涡模拟)方法。
7.2.1 雷诺平均与大涡模拟的对比因为直接求解NS 方程非常困难,所以通常用两种办法对湍流进行模拟,即对NS 方程进行雷诺平均和滤波处理。
这两种方法都会增加新的未知量,因此需要相应增加控制方程的数量,以便保证未知数的数量与方程数量相同,达到封闭方程组的目的。
雷诺平均NS 方程是流场平均变量的控制方程,其相关的模拟理论被称为湍流模式理论。
湍流模式理论假定湍流中的流场变量由一个时均量和一个脉动量组成,以此观点处理NS 方程可以得出雷诺平均NS 方程(简称RNS 方程)。
在引入Boussinesq 假设,即认为湍流雷诺应力与应变成正比之后,湍流计算就归结为对雷诺应力与应变之间的比例系数(即湍流粘性系数)的计算。
根据计算中使用的变量数目和方程数目的不同,湍流模式理论中所包含的湍流模型又被分为二方程模型、一方程模型和零方程模型(代数模型)等大类。
FLUENT 中使用的三种k −ε模型、Spalart-Allmaras 模型、k −ω模型及雷诺应力模型RSM)等都属于湍流模式理论。
大涡模拟(LES)方法是通过滤波处理计算湍流的,其主要思想是大涡结构(又称拟序结构)受流场影响较大,小涡则可以认为是各向同性的,因而可以将大涡计算与小涡计算分开处理,并用统一的模型计算小涡。
大多数飞行器都是在高Re数下飞行,表面的流态是湍流。
为了准确地确定湍流流态下的摩阻、热流,湍流成为一个重要而困难的研究课题。
(一)DNS目前处理湍流数值计算问题有三种方法,第一种方法即所谓直接数值模拟方法(DNS方法),直接求解湍流运动的N-S方程,得到湍流的瞬时流场,即各种尺度的随机运动,可以获得湍流的全部信息。
随着现代计算机的发展和先进的数值方法的研究,DNS方法已经成为解决湍流的一种实际的方法。
但由于计算机条件的约束,目前只能限于一些低Re数的简单流动,不能用于工程应用。
目前国际上正在做的湍流直接数值模拟还只限于较低的需诺数(Re~200)和非常简单的流动外形,如平板边界层、完全发展的槽道流,以及后台阶流动等。
用直接数值模拟方法处理工程中的复杂流动问题,即使是当前最先进的计算机也还差三个量级。
(二)LES另一种方法称做大涡模拟方法(LES方法)。
这是一种折衷的方法,即对湍流脉动部分直接地模拟,将N-S方程在一个小空间域内进行平均(或称之为滤波),以使从流场中去掉小尺度涡,导出大涡所满足的方程。
小涡对大涡的影响会出现在大涡方程中,再通过建立模型(亚格子尺度模型)来模拟小涡的影响。
由于湍流的大涡结构强烈地依赖于流场的边界形状和边界条件,难以找出普遍的湍流模型来描述具有不同的边界特征的大涡结构,宜做直接模拟。
相反地,小尺度涡对边界条件不存在直接依赖关系,而且一般具有各向同性性质。
所以亚格子模型具有更大的普适性,比较容易构造,这是它比雷诺平均方法要优越的地方。
自从1970年Deardorff第一次给出具有工程意义的LES计算以来,LES方法已经成为计算湍流的最强有力的工具之一,应用的方向也在逐步扩展,但是仍然受计算机条件等的限制,使之成为解决大量工程问题的成熟方法仍有很长的路要走。
(三)RANS目前能够用于工程计算的方法就是模式理论。
所谓湍流模式理论,就是依据湍流的理论知识、实验数据或直接数值模拟结果,对Reynolds应力做出各种假设,即假设各种经验的和半经验的本构关系,从而使湍流的平均Reynolds方程封闭。
由于航发燃烧室中的流动特性极其复杂,要想提高数值计算的预测能力,必须要慎重选择湍流模型。
用四种不同的湍流模型对带双径向旋流杯的下游流场进行数值模拟,将计算结果与实验结果作对比,比较各湍流模型的原理和物理基础,优劣,并分析流场速度分布和回流区特性。
涉及的湍流模型:标准k-ε湍流模型(SKE)1标准k-ε湍流模型有较高的稳定性,经济性和计算精度,应用广泛,适合高雷诺数湍流,但不适合旋流等各向异性较强的流动。
2简单的湍流模型是两个方程的模型,需要解两个变量,即速度和长度。
在fluent中,标准k-ε湍流模型自从被Launder and Spalding 提出之后,就变成流场计算中的主要工具。
其在工业上被普遍应用,其计算收敛性和准确性都非常符合工程计算的要求。
3但其也有某些限制,如ε方程包含不能在壁面计算的项,因此必须使用壁面函数。
另外,其预测强分离流,包含大曲率的流动和强压力梯度流动的结果较弱。
它是个半经验的公式,是从实验现象中总结出来的。
动能输运方程是通过精确的方程推导得到,耗散率方程是通过物理推理,数学上模拟相似原型方程得到的。
应用范围:该模型假设流动为完全湍流,分子粘性的影响可以忽略,此标准κ-ε模型只适合完全湍流的流动过程模拟。
可实现的k-ε模型是才出现的,比起标准k-ε模型来有两个主要的不同点:·可实现的k-ε模型为湍流粘性增加了一个公式。
·为耗散率增加了新的传输方程,这个方程来源于一个为层流速度波动而作的精确方程。
术语“realizable”,意味着模型要确保在雷诺压力中要有数学约束,湍流的连续性。
应用范围:可实现的k-ε模型直接的好处是对于平板和圆柱射流的发散比率的更精确的预测。
而且它对于旋转流动、强逆压梯度的边界层流动、流动分离和二次流有很好的表现。
可实现的k-ε模型和RNG k-ε模型都显现出比标准k-ε模型在强流线弯曲、漩涡和旋转有更好的表现。
由于带旋流修正的k-ε模型是新出现的模型,所以还没有确凿的证据表明它比RNG k-ε模型有更好的表现。
流体的湍流模型和湍流模拟流体力学是研究流体的运动规律和性质的学科,其中湍流模型和湍流模拟是其中非常重要的研究方向。
湍流是流体力学中一种复杂而普遍存在的现象,它具有不规则、无序和随机性等特点。
湍流模型和湍流模拟的发展,对于理解和预测真实世界中的湍流现象,以及涉及湍流的工程设计和应用具有重要意义。
一、湍流模型湍流模型是描述湍流现象的数学模型,在流体力学中起着扮演着非常重要的作用。
根据流体力学理论,湍流是由于流体中微小尺度的速度涡旋突然出现和消失所导致的现象。
由于湍流涡旋的尺度范围很广,从而难以直接模拟和计算。
因此,使用湍流模型来近似描述湍流现象,成为了一种常用的方法。
常见的湍流模型包括雷诺平均湍流模型(Reynolds-averaged Navier-Stokes equations, RANS)和大涡模拟(large eddy simulation, LES)等。
雷诺平均湍流模型是基于平均流场的统计性质,通过求解雷诺平均速度和湍流应力来评估湍流效应。
而大涡模拟是将湍流现象分解为不同尺度的涡旋,并通过直接模拟大涡旋来研究湍流运动。
二、湍流模拟湍流模拟是利用计算机来模拟湍流现象的方法,通常基于数值方法对流体力学方程进行求解。
湍流模拟分为直接数值模拟(direct numerical simulation, DNS)、雷诺平均湍流模拟和大涡模拟等。
直接数值模拟是将流场划分为网格,并通过离散化流体力学方程和湍流模型来求解湍流流场的详细信息。
由于该方法需要计算微小尺度的细节,计算量非常大,限制了其在实际工程中的应用。
因此,直接数值模拟主要用于湍流现象的基础研究和理论验证。
相比之下,雷诺平均湍流模拟和大涡模拟能够更有效地模拟湍流现象。
雷诺平均湍流模拟通过对湍流参数进行求解,来描述平均的湍流效应。
而大涡模拟则将湍流现象分为大涡旋和小涡旋,通过模拟大涡旋来捕获湍流流场的主要特征。
三、湍流模型与湍流模拟的应用湍流模型和湍流模拟在工程设计和应用中有着广泛的应用。
流体的湍流模型和雷诺平均法在流体力学中,湍流是一种复杂的流动形式,通常包括随机性和混沌性。
湍流模型是用来描述和计算湍流过程的数学模型。
本文将介绍两种被广泛使用的湍流模型,分别是湍流模型和雷诺平均法。
一、湍流模型湍流模型是用来描述湍流流动的数学方程。
最常用的湍流模型是雷诺应力传输方程,其中流体中的湍流应力可以分解为三个部分:紊流应力、剪切应力和正压力梯度。
湍流模型通过模拟这些应力的传输过程来描述湍流的产生和发展。
湍流模型可以分为两大类:一是基于经验参数的现象学模型,二是基于湍流方程的直接数值模型。
现象学模型是基于实验数据和经验参数的统计方法,可以用于工程实践中。
直接数值模型则是基于湍流方程的求解,可以提供更精确的湍流计算结果。
在实际工程应用中,湍流模型的选择要根据具体情况和需求进行。
一些常用的湍流模型包括:k-ε模型、k-ω模型、Reynolds stress model (RSM)等。
这些模型基于不同的假设和方程,适用于不同的流动条件和复杂性。
二、雷诺平均法雷诺平均法是湍流模拟中的一种重要方法,它通过将流体流动分解为平均分量和脉动分量来描述湍流流动。
在雷诺平均法中,流体的物理量(如速度、压力)被分解为时均分量和涡旋分量。
雷诺平均法通过时间平均和空间平均的操作,将湍流流动中的涡旋分量消除,从而得到平均流动的描述。
利用统计学的方法,雷诺平均法可以获得平均流速、湍流能量和相关涡旋结构等湍流统计信息。
雷诺平均法的优点是可以较好地描述流体流动的平均特性,并具有较低的计算成本。
然而,雷诺平均法忽略了湍流中的空间和时间涨落,对于某些复杂的湍流流动问题,其精度可能不够高。
三、小结湍流是流体力学中常见的复杂流动形式,湍流模型和雷诺平均法是描述和计算湍流流动的重要工具。
湍流模型通过数学方程模拟湍流的产生和发展,可以精确描述湍流流动的特性。
而雷诺平均法则通过将流动分解为平均分量和脉动分量来描述湍流流动的平均特性。
在实际应用中,需要根据具体问题和要求选择合适的湍流模型和计算方法。
湍流模型构建一、湍流模型概述湍流是指流体在运动过程中出现的不规则、无序的运动状态。
由于湍流的不稳定性和复杂性,使得研究湍流问题成为流体力学中的难点之一。
为了描述湍流运动,需要建立适当的数学模型,即湍流模型。
目前常用的湍流模型主要有直接数值模拟(DNS)、大涡模拟(LES)和雷诺平均Navier-Stokes方程(RANS)三种。
二、雷诺平均Navier-Stokes方程1.基本原理雷诺平均Navier-Stokes方程是一种基于统计平均方法来描述湍流运动的数学模型。
该模型假设了在一个足够长时间内,湍流中各个位置上的速度和压力都会发生变化,并且这些变化都是随机性的。
因此,可以通过对时间进行平均来消除这种随机性,并得到一个稳定的平均场。
2.方程形式雷诺平均Navier-Stokes方程包含了连续性方程、动量守恒方程和能量守恒方程三个部分。
其中,连续性方程描述了质量守恒;动量守恒方程描述了动量守恒;能量守恒方程描述了能量守恒。
这三个方程的具体形式如下:连续性方程:$$\frac{\partial \rho}{\partial t}+\nabla \cdot (\rho u)=0$$动量守恒方程:$$\rho \frac{\partial u}{\partial t}+\rho u \cdot \nabla u=-\nabla p+\mu\nabla^2u+\rho g$$能量守恒方程:$$\rho c_p(\frac{\partial T}{\partial t}+u \cdot \nablaT)=\nabla\cdot(k\nabla T)+Q$$其中,$\rho$为流体密度,$u$为流速,$p$为压力,$\mu$为粘性系数,$g$为重力加速度,$c_p$为比热容,$T$为温度,$k$为热导率,$Q$为单位时间内的热源或热汇。
3.湍流模型雷诺平均Navier-Stokes方程中包含了湍流运动的统计平均过程。
大多数飞行器都是在高Re数下飞行,表面的流态是湍流.为了准确地确定湍流流态下的摩阻、热流,湍流成为一个重要而困难的研究课题。
(一)DNS目前处理湍流数值计算问题有三种方法,第一种方法即所谓直接数值模拟方法(DNS方法),直接求解湍流运动的N-S方程,得到湍流的瞬时流场,即各种尺度的随机运动,可以获得湍流的全部信息。
随着现代计算机的发展和先进的数值方法的研究,DNS方法已经成为解决湍流的一种实际的方法。
但由于计算机条件的约束,目前只能限于一些低Re数的简单流动,不能用于工程应用.目前国际上正在做的湍流直接数值模拟还只限于较低的需诺数(Re~200)和非常简单的流动外形,如平板边界层、完全发展的槽道流,以及后台阶流动等。
用直接数值模拟方法处理工程中的复杂流动问题,即使是当前最先进的计算机也还差三个量级.(二)LES另一种方法称做大涡模拟方法(LES方法).这是一种折衷的方法,即对湍流脉动部分直接地模拟,将N—S方程在一个小空间域内进行平均(或称之为滤波),以使从流场中去掉小尺度涡,导出大涡所满足的方程.小涡对大涡的影响会出现在大涡方程中,再通过建立模型(亚格子尺度模型)来模拟小涡的影响。
由于湍流的大涡结构强烈地依赖于流场的边界形状和边界条件,难以找出普遍的湍流模型来描述具有不同的边界特征的大涡结构,宜做直接模拟。
相反地,小尺度涡对边界条件不存在直接依赖关系,而且一般具有各向同性性质。
所以亚格子模型具有更大的普适性,比较容易构造,这是它比雷诺平均方法要优越的地方。
自从1970年Deardorff第一次给出具有工程意义的LES计算以来,LES方法已经成为计算湍流的最强有力的工具之一,应用的方向也在逐步扩展,但是仍然受计算机条件等的限制,使之成为解决大量工程问题的成熟方法仍有很长的路要走。
(三)RANS目前能够用于工程计算的方法就是模式理论。
所谓湍流模式理论,就是依据湍流的理论知识、实验数据或直接数值模拟结果,对Reynolds应力做出各种假设,即假设各种经验的和半经验的本构关系,从而使湍流的平均Reynolds方程封闭.随着计算流体力学的发展,湍流模式理论也有了很大的进步,有了非常丰硕的成果。
从对模式处理的出发点不同,可以将湍流模式理论分类成两大类:一类称为二阶矩封闭模式,另一类称涡粘性封闭模式。
(1)雷诺应力模式所谓二阶矩封闭模式,是从Reynolds应力满足的方程出发,将方程右端未知的项(生成项,扩散项,耗散项等)用平均流动的物理量和湍流的特征尺度表示出来。
典型的平均流动的变量是平均速度和平均温度的空间导数.这种模式理论,由于保留了Reynolds应力所满足的方程,如果模拟的好,可以较好地反映Reynolds应力随空间和时间的变化规律,因而可以较好地反映湍流运动规律。
因此,二阶矩模式是一种较高级的模式,但是,由于保留了Reynolds应力的方程,加上平均运动的方程整个方程组总计15个方程,是一个庞大的方程组,应用这样一个庞大的方程组来解决实际工程问题,计算量很大,这就极大地限制了二阶矩模式在工程问题中的应用。
(2)涡粘性模式在工程湍流问题中得到广泛应用的模式是涡粘性模式。
这是由Boussinesq 仿照分子粘性的思路提出的,即设Reynolds 应力为,ij ij k k i j j i T j i k U U U u u δδν32)32(,,,+++-= () 这里j i u u k 21=是湍动能,T ν称为涡粘性系数,这是最早提出的基准涡粘性模式,即假设雷诺应力与平均速度应变率成线性关系,当平均速度应变率确定后,六个雷诺应力只需要通过确定一个涡粘性系数T ν就可完全确定,且涡粘性系数各向同性,可以通过附加的湍流量来模化,比如湍动能k ,耗散率ε,比耗散率以及其它湍流量ετ/k =,ε/2/3k l =,k q =,根据引入的湍流量的不同,可以得到不同的涡粘性模式,比如常见的ε-k ,k —w 模式,以及后来不断得到发展的τ-k ,q -w,k —l 等模式,涡粘性系数可以分别表示为 ενμ/2k C T =,ωνμkC T =,τνμk C T =,ωνμ2q C T =,.l k C T μν=为了使控制方程封闭,引入多少个附加的湍流量,就要同时求解多少个附加的微分方程,根据求解的附加的微分方程的数目,一般可将涡粘性模式划分为三类:零方程模式,半方程模型,一方程模式,两方程模式。
1) 零方程模式所谓零方程模式是试图直接用平均流动物理量模化T ν,而不引入任何湍流量(如ε,k 等).例如,Prandttl 的混合长理论就是一种零方程模式:yUl T ∂∂∝2ν (5.7) 式中l 称为混合长。
在零方程模式的框架下,得到最为广泛应用的是Baldwin-Lomax 模式[22]。
该模式是对湍流边界层的内层和外层采用不同的混合长假设.这是因为靠近壁面处,湍流脉动受到很大的抑制,含能涡的尺度减小很多,因此长度尺度减小很多;另一方面,在边界层外缘,湍流呈间歇状,质量、动量和能量的输运能力大大下降,即湍流的扩散能力减小。
这样,应用混合长理论来确定涡粘性系数在这两个不同的区域应该有不同的形式。
Baldwin —Lomax 模式的具体数学描述如下。
⎩⎨⎧>≤=contT c inn T T y y y y )()(ννν (5。
8)这里c y 是ont inn )()(T T νν=的离壁面最小距离y 值. 对于内层,即c y y ≤,有Ω=2T )(l inn ν (5.9)Ω是涡量,l U j k ijk ,,ε=Ω是长度尺度))/(1(++--=A y exy ky l (5.10) 其中k=0.4是Karman 常数,A +是模化常数,+y 是无量纲法向距离:w y U y ντ/=+而τu 是摩擦速度,其含义为,wyU u ∂∂=ντ此处下标w 表示壁面. 对于外层,即c y y >,有)(()(T y F F kleb wake out =ν (5.11) 其中)/,min(max 2max max max F U y C F y F dif wk wake =m ax F 是下列函数的最大值: ))/ex p(1()(++--Ω=A y y y F而m ax y 是)(y F 达到最大值的位置。
kleb F 是所谓的Klebanoff 间歇函数:16max )(5.51)(-⎪⎪⎭⎫⎝⎛⋅+=y y C y F kleb kleb dif U 是平均速度分布中最大值和最小值之差。
几个模化常数的值如下:.4.0;,0.1;3.0;02668.0;0.26=====+K C C C A wk kleb由上述模化关系中可以看出,Reynolds 应力完全地由当时当地的平均流参数用代数关系式所决定。
平均流场的任何变化立刻为当地的湍流所感知,这表明零方程模式是一个平衡态模式,假定湍流运动永远处于和平均运动的平衡之中。
实际上对大多数湍流运动而言,并非如此,特别是对平均流空间和时间有剧烈变化的情形,再有因为坐标y 显式地出现在湍流模式中,零方程模式不具有张量不变性,当将它应用到复杂几何外形的流动的数值模拟会带来困难。
当流动发生分离时,Baldwin —Lomax 模式会遇到困难,这是因为在分离点和再附点附近,摩擦速度τu 为零,此时要引入一些人为的干涉来消除这些困难.计算实践表明,只要流动是附体的,零方程模式一般都可以较好地确定压强分布,但是摩阻和传热率的估算不够准确,特别是当流动有分离和再附时。
这是因为附体流压强分布对湍流应力不敏感。
总之,对附体流动,如果只关心压强分布,应用零方程模式通常可以给出满意的结果,而且模式应用起来十分简便。
但是对于我们计算摩阻的需求,零方程模式是不能满足要求。
对于有分离、再附等复杂流动,零方程模式是不适用的。
2)半方程模式为了能计算具有较强压强梯度,特别是较强逆压梯度的非平衡湍流边界层,Johnson-King于1985年提出了一个非平衡代数模型,该模型仍采用涡粘性假设,把涡粘性的分布与最大剪切应力联系在一起,内层涡粘性与外层涡粘性分布用一个指数函数作光滑拟合,外层涡粘性系数作为一个自由参数,由描述最大剪切应力沿流向变化的常微分方程来确定,此常微分方程是由湍流动能方程导出的,故此模型又称为半方程模型。
JK模型虽然仍采用涡粘性假设,却包含有雷诺应力模型的特点.由于求解常微分方程比一方程,二方程模型中求解偏微分方程要简单,省时的多,故用JK模型的工作量只略高于通常平衡状态的零方程代数模型的工作量JK模型后又经不断修正,发展了JK1990A,JK1990J以及JK1992等改进型3)一方程模式Baldwin—Barth(BB) 模型是在二方程模型中,将某一导出的应变量作为基本物理量而得到的,应用此一方程模型可避免求解两方程时会遇到的某些数值困难。
BB一方程模型所选择的导出应变量为“湍流雷诺数”Rt。
BB模型对计算网格的要求低,壁面的网格可以与采用BL代数模型的相当,而不象两方程k-e模型那样要求壁面网格很细,这样就避免了在k-e模型中流场求解的刚性问题。
Spalart-Allmaras(SA)模型与BB模型不同,不是直接利用k—e模型两方程模型加于简化而得,而是从经验和量纲分析出发,由针对简单流动在逐渐补充发展而适用于带有层ν相关的量ν~,流流动的固壁湍流流动的一方程模型,模型中选用的应变量是与涡粘性T除在粘性次层外,ν~与Tν是相等的。
上述两种一方程模型具有相似的特点,它们不象代数模型那样需要分为内层模型,外层模型或壁面模型,尾流模型,同时亦不需要沿法向网格寻找最大值,因此易于用到非结构网格中去;但由于在每个时间步长内,需要对整个流场求解一组偏微分方程,故比BL和JK模型更费机时4)两方程模式2.1 k-两方程模式2。
1。
1 标准k—两方程模式k-模式是最为人所知和应用最广泛的两方程涡粘性模式,为积分到壁面的不可压缩/可压缩湍流的两方程涡粘性模式,各种不同版本的k—模式常见于各种文献中,选择Jones-Launder模式作为一般性介绍。
k-模式最初的发展是为了改善混合长(mixing—length)模式和避免复杂流动中湍流长度尺度(turbulent length scale)的代数表示(algebraic prescription)。
它求解两个湍流标量k和的输运方程.k方程表示湍动能输运方程,方程表示湍动能的耗散率。
该模式对较小压力梯度(relatively small pressure gradients)下的自由剪切流(free —shear-layer flows )具有较好的结果。
对于壁面流动(wall bounded flows ),在零或者小平均压力梯度下,模式结果和实验结果符合得较为一致,但是对大的逆压梯度(adverse pressure gradients ),其结果就不太正确了。