CFD 基 础(流体力学).doc

《CFD 基 础(流体力学).doc》由会员分享，可在线阅读，更多相关《CFD 基 础(流体力学).doc（53页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、-作者xxxx-日期xxxxCFD 基 础(流体力学)【精品文档】第1章 CFD 基 础计算流体动力学(computational fluid dynamics，CFD)是流体力学的一个分支，它通过计算机模拟获得某种流体在特定条件下的有关信息，实现了用计算机代替试验装置完成“计算试验”，为工程技术人员提供了实际工况模拟仿真的操作平台，已广泛应用于航空航天、热能动力、土木水利、汽车工程、铁道、船舶工业、化学工程、流体机械、环境工程等 领域。本章介绍CFD一些重要的基础知识，帮助读者熟悉CFD的基本理论和基本概念，为计算时设置边界条件、对计算结果进行分析与整理提供参考。1.1 流体力学的基本概念1

2、.1.1 流体的连续介质模型流体质点(fluid particle)：几何尺寸同流动空间相比是极小量，又含有大量分子的微 元体。连续介质(continuum/continuous medium)：质点连续地充满所占空间的流体或固体。连续介质模型(continuum/continuous medium model)：把流体视为没有间隙地充满它所占据的整个空间的一种连续介质，且其所有的物理量都是空间坐标和时间的连续函数的一种假设模型：u =u(t,x,y,z)。1.1.2 流体的性质1. 惯性惯性(fluid inertia)指流体不受外力作用时，保持其原有运动状态的属性。惯性与质量有关，质量越大

3、，惯性就越大。单位体积流体的质量称为密度(density)，以r表示，单位为kg/m3。对于均质流体，设其体积为V，质量为m，则其密度为 (1-1)对于非均质流体，密度随点而异。若取包含某点在内的体积，其中质量，则该点密度需要用极限方式表示，即 (1-2)2. 压缩性作用在流体上的压力变化可引起流体的体积变化或密度变化，这一现象称为流体的可压缩性。压缩性(compressibility)可用体积压缩率k来量度 (1-3)式中：p为外部压强。在研究流体流动过程中，若考虑到流体的压缩性，则称为可压缩流动，相应地称流体为可压缩流体，例如高速流动的气体。若不考虑流体的压缩性，则称为不可压缩流动，相应地

4、称流体为不可压缩流体，如水、油等。3. 粘性粘性(viscosity)指在运动的状态下，流体所产生的抵抗剪切变形的性质。粘性大小由粘度来量度。流体的粘度是由流动流体的内聚力和分子的动量交换所引起的。粘度有动力粘度和运动粘度之分。动力粘度由牛顿内摩擦定律导出： (1-4)式中：为切应力，Pa；为动力粘度，Pas；为流体的剪切变形速率。运动粘度与动力粘度的关系为 (1-5)式中：为运动粘度，m2/s。在研究流体流动过程中，考虑流体的粘性时，称为粘性流动，相应的流体称为粘性流体；当不考虑流体的粘性时，称为理想流体的流动，相应的流体称为理想流体。根据流体是否满足牛顿内摩擦定律，将流体分为牛顿流体和非牛

5、顿流体。牛顿流体严格满足牛顿内摩擦定律且保持为常数。非牛顿流体的切应力与速度梯度不成正比，一般又分为塑性流体、假塑性流体、胀塑性流体3种。塑性流体，如牙膏等，它们有一个保持不产生剪切变形的初始应力，只有克服了这个初始应力后，其切应力才与速度梯度成正比，即 (1-6)假塑性流体，如泥浆等，其切应力与速度梯度的关系是 (1-7)胀塑性流体，如乳化液等，其切应力与速度梯度的关系是 (1-8)1.1.3 流体力学中的力与压强1. 质量力与流体微团质量大小有关并且集中在微团质量中心的力称为质量力(body force)。在重力场中有重力mg；直线运动时，有惯性力ma。质量力是一个矢量，一般用单位质量所具

6、有的质量力来表示，其形式如下： (1-9)式中：，为单位质量力在各轴上的投影。2. 表面力大小与表面面积有关而且分布作用在流体表面上的力称为表面力(surface force)。表面力按其作用方向可以分为两种：一是沿表面内法线方向的压力，称为正压力；另一种是沿表面切向的摩擦力，称为切向力。对于理想流体的流动，流体质点只受到正压力，没有切向力；对于粘性流体的流动，流体质点所受到的作用力既有正压力，也有切向力。作用在静止流体上的表面力只有沿表面内法线方向的正压力。单位面积上所受到的表面力称为这一点处的静压强。静压强具有两个特征：静压强的方向垂直指向作用面； 流场内一点处静压强的大小与方向无关。3.

7、 表面张力在液体表面，界面上液体间的相互作用力称为张力。在液体表面有自动收缩的趋势，收缩的液面存在相互作用的与该处液面相切的拉力，称为液体的表面张力(surface tension)。正是这种力的存在，引起弯曲液面内外出现压强差以及常见的毛细现象等。试验表明，表面张力大小与液面的截线长度L成正比，即 (1-10)式中：为表面张力系数，它表示液面上单位长度截线上的表面张力，其大小由物质种类决定，其单位为N/m。4. 绝对压强、相对压强及真空度标准大气压的压强是101325Pa(760mm汞柱)，通常用patm表示。若压强大于大气压，则以该压强为计算基准得到的压强称为相对压强(relative p

8、ressure)，也称为表压强，通常用pr表示。若压强小于大气压，则压强低于大气压的值就称为真空度(vacuum)，通常用pv表示。如以压强0Pa为计算的基准，则这个压强就称为绝对压强(absolute pressure)，通常用ps表示。这三者的关系如下： (1-11) (1-12)在流体力学中，压强都用符号p表示，但一般来说有一个约定：对于液体，压强用相对压强；对于气体，特别是马赫数大于0.1的流动，应视为可压缩流，压强用绝对压强。压强的单位较多，一般用Pa，也可用bar，还可以用汞柱、水柱，这些单位换算如下：1Pa=1N/m21bar=105Pa1patm2O=101325Pa5. 静压

9、、动压和总压对于静止状态下的流体，只有静压强。对于流动状态的流体，有静压强(static pressure)、动压强(dynamic pressure)、测压管压强(manometric tube pressure)和总压强(total pressure)之分。下面从伯努利(Bernoulli)方程(也有人称其为伯努里方程)中分析它们的意义。伯努利方程阐述一条流线上流体质点的机械能守恒，对于理想流体的不可压缩流动其表达式如下： (1-13)式中：称为压强水头，也是压能项，为静压强；称为速度水头，也是动能项；称为位置水头，也是重力势能项，这三项之和就是流体质点的总的机械能；H称为总的水头高。将式

10、(1-13)两边同时乘以，则有 (1-14)式中：称为静压强，简称静压；称为动压强，简称动压；称为总压强，简称总压。对于不考虑重力的流动，总压就是静压和动压之和。1.1.4 流体运动的描述1. 流体运动描述的方法描述流体物理量有两种方法，一种是拉格朗日描述；一种是欧拉描述。拉格朗日(Lagrange)描述也称随体描述，它着眼于流体质点，并将流体质点的物理量认为是随流体质点及时间变化的，即把流体质点的物理量表示为拉格朗日坐标及时间的函数。设拉格朗日坐标为(a,b,c)，以此坐标表示的流体质点的物理量，如矢径、速度、压强等等在任一时刻t的值，便可以写为a、b、c及t的函数。若以f表示流体质点的某一

11、物理量，其拉格朗日描述的数学表达式为 (1-15)例如，设时刻t流体质点的矢径即t时刻流体质点的位置以r表示，其拉格朗日描述为 (1-16)同样，质点的速度的拉格朗日描述是 (1-17)欧拉描述，也称空间描述，它着眼于空间点，认为流体的物理量随空间点及时间而变化，即把流体物理量表示为欧拉坐标及时间的函数。设欧拉坐标为(q1,q2,q3)，用欧拉坐标表示的各空间点上的流体物理量如速度、压强等，在任一时刻t的值，可写为q1、q2、q3及t的函数。从数学分析知道，当某时刻一个物理量在空间的分布一旦确定，该物理量在此空间形成一个场。因此，欧拉描述实际上描述了一个个物理量的场。若以f表示流体的一个物理量

12、，其欧拉描述的数学表达式是(设空间坐标取用直角坐标) (1-18)如流体速度的欧拉描述是 (1-19)2. 拉格朗日描述与欧拉描述之间的关系拉格朗日描述着眼于流体质点，将物理量视为流体坐标与时间的函数；欧拉描述着眼于空间点，将物理量视为空间坐标与时间的函数。它们可以描述同一物理量，必定互相相关。设表达式表示流体质点(a,b,c)在t时刻的物理量；表达式表示空间点(x,y,z)在时刻t的同一物理量。如果流体质点(a,b,c)在t时刻恰好运动到空间点(x,y,z)上，则应有 (1-20) (1-21)事实上，将式(1-16)代入式(1-21)左端，即有 (1-22)或者反解式(1-16)，得到 (

13、1-23)将式(1-23)代入式(1-21)的右端，也应有 (1-24)由此，可以通过拉格朗日描述推出欧拉描述，同样也可以由欧拉描述推出拉格朗日 描述。3. 随体导数流体质点物理量随时间的变化率称为随体导数(substantial derivative)，或物质导数、质点导数。按拉格朗日描述，物理量f表示为，f的随体导数就是跟随质点(a,b,c)的物理量f对时间t的导数。例如，速度是矢径对时间的偏导数， (1-25)即随体导数就是偏导数。按欧拉描述，物理量f表示为，但并不表示随体导数，它只表示物理量在空间点上的时间变化率。而随体导数必须跟随t时刻位于空间点上的那个流体质点，其物理量f的时间变化

14、率。由于该流体质点是运动的，即x、y、z是变的，若以a、b、c表示该流体质点的拉格朗日坐标，则x、y、z将依式(1-16)变化，从而f=F(x,y,z,t)的变化依连锁法则处理。因此，物理量f=F(x,y,z,t)的随体导数是 (1-26)式中：表示随体导数。从中可以看出，对于质点物理量的随体导数，欧拉描述与拉格朗日描述大不相同。前者是两者之和，而后者是直接的偏导数。4. 定常流动与非定常流动根据流体流动过程以及流动过程中的流体的物理参数是否与时间相关，可将流动分为定常流动(steady flow)与非定常流动(unsteady flow)。定常流动：流体流动过程中各物理量均与时间无关，这种流

15、动称为定常流动。非定常流动：流体流动过程中某个或某些物理量与时间有关，则这种流动称为非定常流动。5. 流线与迹线常用流线和迹线来描述流体的流动。迹线(track)：随着时间的变化，空间某一点处的流体质点在流动过程中所留下的痕迹称为迹线。在t=0时刻，位于空间坐标(a,b,c)处的流体质点，其迹线方程为 (1-27)式中：u、v、w分别为流体质点速度的三个分量；x、y、z为在t时刻此流体质点的空间 位置。流线(streamline)：在同一个时刻，由不同的无数多个流体质点组成的一条曲线，曲线上每一点处的切线与该质点处流体质点的运动方向平行。流场在某一时刻t的流线方程为 (1-28)对于定常流动，

16、流线的形状不随时间变化，而且流体质点的迹线与流线重合。在实际流场中除驻点或奇点外，流线不能相交，不能突然转折。6. 流量与净通量流量(flux)：单位时间内流过某一控制面的流体体积称为该控制面的流量Q，其单位为m3/s。若单位时间内流过的流体是以质量计算，则称为质量流量Qm；不加说明时“流量”一词概指体积流量。在曲面控制面上有 (1-29)净通量(net flux)：在流场中取整个封闭曲面作为控制面A，封闭曲面内的空间称为控制体。流体经一部分控制面流入控制体，同时也有流体经另一部分控制面从控制体中流出，此时流出的流体减去流入的流体，所得出的流量称为流过全部封闭控制面A的净流量(或净通量)，通过

17、式(1-30)计算： (1-30)对于不可压缩流体来说，流过任意封闭控制面的净通量等于0。7. 有旋流动与有势流动由速度分解定理，流体质点的运动可以分解为：(1)随同其他质点的平动；(2)自身的旋转运动；(3)自身的变形运动(拉伸变形和剪切变形)。在流动过程中，若流体质点自身做无旋运动(irrotational flow)，则称流动是无旋的，也就是有势的，否则就称流动是有旋流动(rotational flow)。流体质点的旋度是一个矢量，通常用表示，其大小为 (1-31)若=0，则称流动为无旋流动，否则就是有旋流动。与流体的流线或迹线形状无关；粘性流动一般为有旋流动；对于无旋流动，伯努利方程适

18、用于流场中任意两点之间；无旋流动也称为有势流动(potential flow)，即存在一个势函数，满足： (1-32)即 (1-33)8. 层流与湍流流体的流动分为层流流动(laminar flow)和湍流流动(turbulent flow)。从试验的角度来看，层流流动就是流体层与层之间相互没有任何干扰，层与层之间既没有质量的传递也没有动量的传递；而湍流流动中层与层之间相互有干扰，而且干扰的力度还会随着流动而加大，层与层之间既有质量的传递又有动量的传递。判断流动是层流还是湍流，是看其雷诺数是否超过临界雷诺数。雷诺数的定义如下： (1-34)式中：V为截面的平均速度；L为特征长度；为流体的运动粘

19、度。对于圆形管内流动，特征长度L取圆管的直径d。一般认为临界雷诺数为2320，即 (1-35)当Re2320时，管中是湍流。对于异型管道内的流动，特征长度取水力直径dH，则雷诺数的表达式为 (1-36)异型管道水力直径的定义如下： (1-37)式中：A为过流断面的面积；S为过流断面上流体与固体接触的周长。临界雷诺数根据形状的不同而有所差别。根据试验几种异型管道的临界雷诺数如 表1-1所示。表1-1 几种异型管道的临界雷诺数管道截面形状正方形正三角形偏心缝隙207019301000对于平板的外部绕流，特征长度取沿流动方向的长度，其临界雷诺数为51053106。1.2 CFD基本模型流体流动所遵循

20、的物理定律，是建立流体运动基本方程组的依据。这些定律主要包括质量守恒、动量守恒、动量矩守恒、能量守恒、热力学第二定律，加上状态方程、本构方程。在实际计算时，还要考虑不同的流态，如层流与湍流。1.2.1 基本控制方程1. 系统与控制体在流体力学中，系统是指某一确定流体质点集合的总体。系统以外的环境称为外界。分隔系统与外界的界面，称为系统的边界。系统通常是研究的对象，外界则用来区别于系统。系统将随系统内质点一起运动，系统内的质点始终包含在系统内，系统边界的形状和所围空间的大小可随运动而变化。系统与外界无质量交换，但可以有力的相互作用，及能量(热和功)交换。控制体是指在流体所在的空间中，以假想或真实

21、流体边界包围，固定不动形状任意的空间体积。包围这个空间体积的边界面，称为控制面。控制体的形状与大小不变，并相对于某坐标系固定不动。控制体内的流体质点组成并非不变的。控制体既可通过控制面与外界有质量和能量交换，也可与控制体外的环境有力的相互作用。2. 质量守恒方程(连续性方程)在流场中，流体通过控制面A1流入控制体，同时也会通过另一部分控制面A2流出控制体，在这期间控制体内部的流体质量也会发生变化。按照质量守恒定律，流入的质量与流出的质量之差，应该等于控制体内部流体质量的增量，由此可导出流体流动连续性方程的积分形式为 (1-38)式中：V表示控制体，A表示控制面。等式左边第一项表示控制体V内部质

22、量的增量；第二项表示通过控制表面流入控制体的净通量。根据数学中的奥-高公式，在直角坐标系下可将其化为微分形式： (1-39)对于不可压缩均质流体，密度为常数，则有 (1-40)对于圆柱坐标系，其形式为 (1-41)对于不可压缩均质流体，密度为常数，则有 (1-42)3. 动量守恒方程(运动方程)动量守恒是流体运动时应遵循的另一个普遍定律，描述为：在一给定的流体系统，其动量的时间变化率等于作用于其上的外力总和，其数学表达式即为动量守恒方程，也称为运动方程，或N-S方程，其微分形式表达如下： (1-43)式中：、分别是单位质量流体上的质量力在三个方向上的分量；是流体内应力张量的分量。动量守恒方程在

23、实际应用中有许多表达形式，其中比较常见的有如下几种。(1)可压缩粘性流体的动量守恒方程 (1-44)(2)常粘性流体的动量守恒方程 (1-45)(3)常密度常粘性流体的动量守恒方程 (1-46)(4)无粘性流体的动量守恒方程(欧拉方程) (1-47)(5)静力学方程 (1-48)(6)相对运动方程在非惯性参考系中的相对运动方程是研究像大气、海洋及旋转系统中流体运动的所必须考虑的。由理论力学得知，绝对速度为相对速度及牵连速度之和，即 (1-49)其中，为运动系中的平动速度，是其转动角速度，为质点矢径。而绝对加速度为相对加速度、牵连加速度及科氏加速度之和，即 (1-50)其中，。将绝对加速度代入运

24、动方程，即得到流体的相对运动方程 (1-51)4. 能量守恒方程将热力学第一定律应用于流体运动，把式(1-51)各项用有关的流体物理量表示出来，即是能量方程。如式(1-52)所示。 (1-52)式中：；是有效热传导系数，其中是湍流热传导系数，根据所使用的湍流模型来定义；是组分j的扩散流量；包括了化学反应热以及其他用户定义的体积热源项；方程右边的前3项分别描述了热传导、组分扩散和粘性耗散带来的能量输运。1.2.2 湍流模型 湍流是自然界广泛存在的流动现象。大气、海洋环境的流动，飞行器和船舰的绕流，叶轮机械、化学反应器、核反应器中的流体运动都是湍流。湍流流动的核心特征是其在物理上近乎于无穷多的尺度

25、和数学上强烈的非线性，这使得人们无论是通过理论分析、实验研究还是计算机模拟来彻底认识湍流都非常困难。回顾计算流体力学的发展，特别是活跃的20世纪80年代，不仅提出和发展了一大批高精度、高分辨率的计算格式，从主控方程看相当成功地解决了欧拉方程的数值模拟，可以说欧拉方程数值模拟方法的精度已接近于它有效使用范围的极限；同时还发展了一大批有效的网格生成技术及相应的软件，具体实现了工程计算所需要的复杂外形的计算网格；且随着计算机的发展，无论从计算时间还是从计算费用考虑，欧拉方程都已能适用于各种实践所需。在此基础上，20世纪80年代还进行了求解可压缩雷诺平均方程及其三维定态粘流流动的模拟。20世纪90年代

26、又开始一个非定常粘流流场模拟的新局面，这里所说的粘流流场具有高雷诺数、非定常、不稳定、剧烈分离流动的特点，显然需要继续探求更高精度的计算方法和更实用可靠的网格生成技术。但更为重要的关键性的决策将是，研究湍流机理，建立相应的模式，并进行适当的模拟仍是解决湍流问题的重要途径。1. 湍流模型分类湍流流动模型很多，但大致可以归纳为以下3类。第一类是湍流输运系数模型，即将速度脉动的二阶关联量表示成平均速度梯度与湍流粘性系数的乘积，用笛卡儿张量表示为 (1-53)模型的任务就是给出计算湍流粘性系数的方法。根据建立模型所需要的微分方程的数目，可以分为零方程模型(代数方程模型)、单方程模型和双方程模型。第二类

27、是抛弃了湍流输运系数的概念，直接建立湍流应力和其他二阶关联量的输运 方程。第三类是大涡模拟。前两类是以湍流的统计结构为基础，对所有涡旋进行统计平均。大涡模拟把湍流分成大尺度湍流和小尺度湍流，通过求解三维经过修正的Navier-Stokes方程(纳维-斯托克斯方程，简称N-S方程)，得到大涡旋的运动特性，而对小涡旋运动还采用上述的模型。实际求解中，选用什么模型要根据具体问题的特点来决定。选择的一般原则是精度要高，应用简单，节省计算时间，同时也具有通用性。Fluent 提供的湍流模型包括：单方程(Spalart-Allmaras)模型、双方程模型(标准模型、重整化群模型、可实现模型)及雷诺应力模型

28、和大涡模拟，如图1-1所示。图1-1 湍流模型详解2. 平均量输运方程雷诺平均就是把Navier-Stokes方程中的瞬时变量分解成平均量和脉动量两部分。对于速度，有 (1-54)式中：和分别是平均速度和脉动速度()。类似地，对于压力等其他标量，也有 (1-55)式中：表示标量，如压力、能量、组分浓度等。把上面的表达式代入瞬时的连续与动量方程，并取平均(去掉平均速度上的横线)，可以把连续与动量方程写成如下的笛卡儿坐标系下的张量形式： (1-56) (1-57)上面两个方程称为雷诺平均的Navier-Stokes(RANS)方程。它们和瞬时Navier-Stokes方程有相同的形式，只是速度或其

29、他求解变量变成了时间平均量。额外多出来的项是雷诺应力，表示湍流的影响。对于密度变化的流动过程，如燃烧问题，需要采用法夫雷(Favre)平均才可以求解。法夫雷平均就是除了压力和密度本身以外，所有变量都用密度加权平均。变量的密度加权平均定义如下： (1-58)式中：符号表示密度加权平均，对应于密度加权平均值的脉动值用表示，有。显然，这种脉动值的简单平均值不为零，但它的密度加权平均值等于零，即。为了求解方程(1-57)，必须模拟雷诺应力项以使方程封闭。通常的方法是应用Boussinesq假设，认为雷诺应力与平均速度梯度成正比，表达式如下： (1-59)Boussinesq假设被用于单方程模型和双方程

30、模型。这种近似方法好处是与求解湍流粘性系数有关的计算时间比较少。例如，在Spalart-Allmaras单方程模型中只多求解一个表示湍流粘性的输运方程；在双方程模型中只需多求解湍动能k和耗散率两个方程，湍流粘性系数用湍动能k和耗散率的函数来描述。Boussinesq假设的不足之处是假设是个各向同性标量，对于一些复杂流动，该条件并不是严格成立，所以具有其应用局限性。另外的近似方法是求解雷诺应力各分量的输运方程。这也需要额外再求解一个标量方程，通常是耗散率方程。这就意味着对于二维湍流流动问题，需要多求解4个输运方程，而三维湍流问题需要多求解7个方程，需要较多的计算时间，要求更高的计算机内存。在很多

31、情况下基于Boussinesq假设的模型很好用，而且计算量并不是很大。但是，如果湍流场各向异性很明显，如强旋流动以及应力取得的二次流等流动中，求解RSM模型可以得到更好的结果。3. 常用湍流模型简介1)单方程(Spalart-Allmaras)模型单方程模型求解变量是，表征出了近壁(粘性影响)区域以外的湍流运动粘性系数。的输运方程为 (1-60)式中：是湍流粘性产生项；是由于壁面阻挡与粘性阻尼引起的湍流粘性的减少；和是常数；是分子运动粘性系数。湍流粘性系数，其中，是粘性阻尼函数，定义为，。而湍流粘性产生项模拟为，其中，和是常数，是计算点到壁面的距离；，。在Fluent软件中，考虑到平均应变率对

32、湍流产生也起到很大作用，其中，=2.0，平均应变率。在涡量超过应变率的计算区域计算出来的涡旋粘性系数变小。这适合涡流靠近涡旋中心的区域，那里只有“单纯”的旋转，湍流受到抑止。包含应变张量的影响更能体现旋转对湍流的影响。忽略了平均应变，估计的涡旋粘性系数产生项偏高。湍流粘性系数减少项为，其中，、是常数，在计算r时用到的受平均应变率的影响。上面的模型常数在Fluent软件中默认值为，。2)标准模型标准模型需要求解湍动能及其耗散率方程。湍动能输运方程是通过精确的方程推导得到的，但耗散率方程是通过物理推理，数学上模拟相似原形方程得到的。该模型假设流动为完全湍流，分子粘性的影响可以忽略。因此，标准模型只

33、适合完全湍流的流动过程模拟。标准模型的湍动能k和耗散率方程为如下形式： (1-61) (1-62)式中：表示由于平均速度梯度引起的湍动能产生，表示由于浮力影响引起的湍动能产生；表示可压缩湍流脉动膨胀对总的耗散率的影响。湍流粘性系数。在Fluent中，作为默认值常数，=， =， = ，湍动能k 与耗散率的湍流普朗特数分别为=，=。3)重整化群模型重整化群模型是对瞬时的Navier-Stokes方程用重整化群的数学方法推导出来的模型。模型中的常数与标准模型不同，而且方程中也出现了新的函数或者项。其湍动能与耗散率方程与标准模型有相似的形式： (1-63) (1-64)式中：表示由于平均速度梯度引起的

34、湍动能产生，表示由于浮力影响引起的湍动能产生；表示可压缩湍流脉动膨胀对总的耗散率的影响，这些参数与标准模型中相同。和分别是湍动能k 和耗散率的有效湍流普朗特数的倒数。湍流粘性系数计算公式为，其中，。对于前面方程的积分，可以精确到有效雷诺数(涡旋尺度)对湍流输运的影响，这有助于处理低雷诺数和近壁流动问题的模拟。对于高雷诺数，上面方程可以给出：，。这个结果非常有意思，和标准模型的半经验推导给出的常数非常近似。在Fluent中，如果是默认设置，用重整化群模型时是针对的高雷诺数流动问题。如果对低雷诺数问题进行数值模拟，必须进行相应的设置。4)可实现模型可实现模型的湍动能及其耗散率输运方程为 (1-65

35、) (1-66)式中：，。在上述方程中，表示由于平均速度梯度引起的湍动能产生，表示由于浮力影响引起的湍动能产生；表示可压缩湍流脉动膨胀对总的耗散率的影响；和是常数；和分别是湍动能及其耗散率的湍流普朗特数。在Fluent中，作为默认值常数，=，=，=，=。该模型的湍流粘性系数与标准模型相同。不同的是，粘性系数中的不是常数，而是通过公式计算得到，其中， ，表示在角速度旋转参考系下的平均旋转张量率。模型常数，式中，。从这些式子中发现，是平均应变率与旋度的函数。在平衡边界层惯性底层，可以得到，与标准模型中采用的常数一样。该模型适合的流动类型比较广泛，包括有旋均匀剪切流、自由流(射流和混合层)、腔道流动

36、和边界层流动。对以上流动过程模拟结果都比标准模型的结果好，特别是可实现模型对圆口射流和平板射流模拟中，能给出较好的射流扩张角。双方程模型中，无论是标准模型、重整化群模型还是可实现模型，三个模型有类似的形式，即都有k和的输运方程，它们的区别在于：计算湍流粘性的方法不同；控制湍流扩散的湍流普朗特数不同；方程中的产生项和Gk关系不同。但都包含了相同的表示由于平均速度梯度引起的湍动能产生，表示由于浮力影响引起的湍动能产生；表示可压缩湍流脉动膨胀对总的耗散率的影响。湍动能产生项 (1-67) (1-68)式中：是能量的湍流普特朗数，对于可实现模型，默认设置值为；对于重整化群模型，。热膨胀系数，对于理想气

37、体，浮力引起的湍动能产生项变为 (1-69)5)雷诺应力模型雷诺应力模型(RSM)是求解雷诺应力张量的各个分量的输运方程。具体形式为 (1-70)式中：左边的第二项是对流项，右边第一项是湍流扩散项，第二项是分子扩散项，第三项是应力产生项，第四项是浮力产生项，第五项是压力应变项，第六项是耗散项，第七项系统旋转产生项。在式(1-69)中，、不需要模拟，而、需要模拟以封闭方程。下面简单对几个需要模拟项进行模拟。可以用Delay和Harlow的梯度扩散模型来模拟，但这个模型会导致数值不稳定，在Fluent中是采用标量湍流扩散模型： (1-71)式中：湍流粘性系数用来计算，根据Lien和Leschzin

38、er，这和标准模型中选取1.0有所不同。压力应变项可以分解为三项，即 (1-72)式中：、和分别是慢速项、快速项和壁面反射项，具体表述可以参见文献2。浮力引起的产生项模拟为 (1-73)耗散张量模拟为 (1-74)式中：，是马赫数；标量耗散率用标准模型中采用的耗散率输运方程求解。6)大涡模拟湍流中包含了不同时间与长度尺度的涡旋。最大长度尺度通常为平均流动的特征长度尺度。最小尺度为Komogrov尺度。LES的基本假设是：动量、能量、质量及其他标量主要由大涡输运；流动的几何和边界条件决定了大涡的特性，而流动特性主要在大涡中体现；小尺度涡旋受几何和边界条件影响较小，并且各向同性，大涡模拟(LES)

39、过程中，直接求解大涡，小尺度涡旋模拟，从而使得网格要求比DNS低。LES的控制方程是对Navier-Stokes方程在波数空间或者物理空间进行过滤得到的。过滤的过程是去掉比过滤宽度或者给定物理宽度小的涡旋，从而得到大涡旋的控制方程： (1-75) (1-76)式中：为亚网格应力，。很明显，上述方程与雷诺平均方程很相似，只不过大涡模拟中的变量是过滤过的量，而非时间平均量，并且湍流应力也不同。1.2.3 初始条件和边界条件计算流体动力学(CFD)分析中，初始条件和边界条件的正确设置是关键的一步。现有的CFD软件都提供了现成的各种类型的边界条件，这里对有关的初始条件和边界条件作一般讨论。1. 初始条

40、件顾名思义，初始条件就是计算初始给定的参数，即时给出各未知量的函数分布，如 (1-77)很明显，当流体运动定常时，无初始条件问题。2. 边界条件所谓边界条件就是流体力学方程组在求解域的边界上，流体物理量应满足的条件。例如，流体被固壁所限，流体将不应有穿过固壁的速度分量；在水面这个边界上，大气压强认为是常数(一般在距离不大的范围内可如此)；在流体与外界无热传导的边界上，流体与边界之间无温差，如此等。由于各种具体问题不同，边界条件提法千差万别，一般要保持恰当：保持在物理上是正确的；要在数学上不多不少，刚好能用来确定积分微分方程中的积分常数，而不是矛盾的或有随意性。通常流体边界分为流固交界面和流流(

41、液液、液气)交界面，下面分别讨论。1)流固分界面边界条件飞机、船舶在空气及水中运动时的流固分界面，水在岸边及底部的流固分界面，均属这一类。一般而言，流体在固体边界上的速度依流体有无粘性而定。对于粘性流体，流体将粘附于固体表面(无滑移)，即 (1-78)式中：是流体速度；是固壁面相应点的速度。式(1-78)表明，在流固边界面上，流体在一点的速度等于固体在该点的速度。对于无粘性流体，流体可沿界面滑移，即有速度的切向分量，但不能离开界面，也就是流体的法向速度分量等于固体的法向速度分量，即 (1-79)另外，也可视所给条件，给出无温差条件： (1-80)式中：是流体温度，是固壁面相应点的温度。2)液液分界面边界条件密度不同的两种液体的分界面就属于这一类。一般而言，对分界面两侧的液体情况经常给出的条件是 (1-81)对应力及传导热情况给出的条件是 (1-82) (1-83)3)液气分界面边界条件液气分界面最典型的是水与大气的分界面，即自由面。由于自由面本身是运动和变形的，而且其形状常常也是一个需要求解的未