下载文档原格式
湍流模型概述湍流是一种复杂的非稳态三维流动,通常把瑞流定义为具有随机性、扩散性、高雷诺数、三维祸量脉动性、耗散性及连续性特征的复杂流动。
虽然瑞流具有多种特性,但瑞流不是流体本身具有的某些特征而是流体流动的特征,仍是一种连续流动,仍然同层流一样满足流动的基本方程。
从数学的观点看,瑞流是N-S方程的通解,求解端流与求解层流无本质区别,目前己具有足以求解瑞流问题的有关方程式。
端流还可以看作是由多种大尺度祸流和小尺度祸流组成的特殊流动。
大尺度的祸流主要由流动的边界条件和流动区域的几何形状所决定,是引起流场中低频脉动的主要原因;小尺度的祸流主要是點性力所决定,是引起流场中高频脉动的主要原因。
瑞流的物理量的脉动特点就是由于流体内各种不同尺度祸流的随机运动造成。
用数值方法直接计算瑞流单元运动规律时,计算网格尺寸要小于瑞流单元尺度,并在瑞流单元尺度内计算N-S方程的通解。
但是在实际工程中具有重要意义的不是端流的精细结构,而是瑞流对于时间的平均(时均)效应。
因此,雷诺首先提出了将N-S方程对某一时间比例尺取平均,得到时均N-S方程。
虽然瑞流的N-S方程经过时均化处理后方程式的形式可以保持不变,但是出现了脉动应力项(雷诺应力),因此需要提出相应的端流模型(一个或一组数学方程)使时均方程得到封闭。
这种方法按雷诺应力方程模型化方法的不同可分为两类:一类是直接就雷诺应力建立模型化方程的雷诺应力方程模型;另一类是在雷诺应力与局部时均速度梯度成比例的Boussinesq假设下引入的瑞流黏度系数模型。
另一种瑞流数值计算方法是亚网格尺度模拟,即大祸模拟(LES),由N-S方程出发直接模拟大尺度祸流,小尺度祸流的影响可以通过近似模型来考虑。
但是由于大祸模拟计算量仍很大,也只能模拟一些简单的情况。
工程上通常需要深入了解的是温度场、时均速度场、瑞流脉动时均特性等,并不需要了解瑞流产生和发展的详细过程。
因此,利用雷诺提出的时均值的概念来研究瑞流运动的方法是一种有效的简化,从N-S方程导出瑞流平均运动方程和雷诺方程,还导出了连续性方程和能量方程等基本方程。
湍流模型目前计算流体力学常用的湍流的数值模拟方法主要有以下三种：直接模拟（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 其中，φ表示标量，如压力、能量、组分浓度等。
湍流模型推导对纳维斯托克斯方程做时间平均处理，即采用雷诺平均法（RANS ：Reynolds-Averaged Navier-Stokes ），可以得到湍流基本方程。
对于任意变量φ，按照雷诺时间平均法，可以拆分为如下格式：φφφ'+=“-” 表示对时间的平均，上标“’”代表脉动量。
按照dt TTt tφφ⎰∆+∆=1计算平均值，将流动变量i u 和p 转换成时间平均和脉动值之和u u u i '+=，p p p '+=为了使方程组更具有封闭性，必须模化雷诺应力，引入模型使方程组封闭。
其方法之一是湍流粘性系数法。
按照基于Boussinesq 的涡粘假设湍流粘性系数法有ij i i t i jj i t j i x u k x u xu u u δμρμρ⎪⎪⎭⎫ ⎝⎛∂∂+-⎪⎪⎭⎫ ⎝⎛∂∂+∂∂=''-32 上述方程式中t μ为涡粘系数，i u 为时均速度，ij δ是Kronecker 符号，k 为湍流动能（当j i =时，1=ij δ；当j i ≠时，0=ij δ）。
2i i u u k ''=确定涡粘性系数t μ就是整个湍流模型的目标关键，确定湍流粘性系数法具体可以分为零方程模型、一方程模型、和二方程模型等等。
一 零方程模型零方程模型也可称作代数模型，直接建立雷诺应力和时均值的代数关系，从而把涡粘系数和时均值联系到一起的模型。
1 混合长度模式混合长度模式是基于分子运动的比拟，在二维剪切层中导出的。
混合长度l 类比分子运动自由程，在经历混合长度的横向距离上，脉动速度正比于混合长度及流向平均速度梯度，即：yUlu ∂∂∝' （1.1-1） 而粘性系数应当正比于脉动速度和混合长度之积（分子粘性系数正比于自由程和分子热运动速度之积），从而涡粘系数有如下的估计式：yUl l u v t ∂∂∝'∝2（1.1-2） 在湍流输运中，涡粘系数和沃扩散系数之比定义为普朗特数t Pr ，即：t t t v κ=Pr （1.1-3）工程计算中通常采用0.1~8.0Pr =t 。
湍流模型介绍
湍流模型是数学模型的一种，用于描述液体或气体中的湍流运动。
湍流是一种不规律的、难以预测的流体运动，通常是由于速度、密度或温度的不规则分布引起的。
湍流模型通过使用一系列方程，描述流体的速度、压力和密度等参数之间的相互作用，以预测和模拟流体的复杂运动行为。
湍流模型主要分为两类：基于雷诺平均的模型（如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）方法是通过滤波处理计算湍流的，其主要思想是大涡结构（又称拟序结构）受流场影响较大，小涡则可以认为是各向同性的，因而可以将大涡计算与小涡计算分开处理，并用统一的模型计算小涡。
湍流模型公式湍流模型公式是研究复杂流场的理论模型，由美国科学家亨利莱茨史密斯（Henry L. Smith）在1941年提出，它可以有效地用于模拟和研究湍流性质。
湍流模型公式也被称为史密斯-史密斯-贝尔（SSB）湍流模型，它是要全面描述湍流的最常用的模型之一。
湍流模型公式的核心是以下一维湍流方程：u_t+(uU)_x=νu_xx+ f(x,t) (1)其中u(x,t)为流体速度，U(x,t)为水流平均速度，ν为流的粘性系数，f(x,t)为外力。
另外，有必要提到的是，湍流模型公式还包括一维流体剪切应力方程和一维湍流能量方程：τ_t+(uτ)_x =ντ_xx+συw (2)E_t+(uE)_x =νE_xx+σu^2w (3) 其中τ(x,t)为流体剪切应力，E(x,t)为湍流能量，σ为流的粘性系数，u(x,t)为流体速度，U(x,t)为水流平均速度，υw为对流湍流损失（eddy viscosity）。
湍流模型公式最引人注目的特点是它准确地描述了湍流动力学的复杂性，因此它可以用来模拟流动场中的湍流运动，以提出解决实际技术问题的方案。
因此，湍流模型公式已广泛应用于工程和科学领域，如水资源开发、船舶控制、海洋工程、航空设计、气动学等。
湍流模型公式的求解用到了数值方法，即有限差分法。
它可以将湍流动力学的三维湍流问题简化为一维湍流问题，然后使用有限求差分法（FD）或有限元法（FEM）对参数进行有效求解。
其主要过程如下：1. 使用有限差分法或有限元法对定义的流动模型进行离散；2.用叙述方法或迭代方法求解湍流模型公式；3. 使用反演方法或最优控制方法优化参数；4.析和模拟湍流的空间和时间变化及其特性。
湍流模型公式的应用范围越来越广泛，它对于实际工程外面的许多流体动力学问题具有重要的理论意义和工程意义，这些问题包括流体运动的抗阻，水动力学的结构优化设计，湍流动态特性，空气动力学的优化控制，结构调节系数的估计，气动学的分析等等。
