计算流体力学4.pptVIP

4.3三维湍流模型及其在CFD中的应用湍流流动是自然界常见的流动现象，在多数工程问题中流体的流动往往处于湍流状态，湍动特性在工程中占有重要的地位．因此，湍流研究一直被研究者高度重视。但由于湍流本身的复杂性，直到现在仍有一些基本问题尚未解决。本章不深入涉及湍流的结构及发生的机理，主要从工程实际应用的观点介绍不可压流体湍流流动与换热的常用数值模拟方法。4.4湍流流动的特征4.4.1湍流及其数学描述流体试验表明，当Reynolds数小于某一临界值时，相邻的流体层彼此有序地流动，这种流动称作层流。当Re数大于临界值时，会出现一系列复杂的变化，最终导致流动特征的本质变化，流动呈无序的混乱状态。这时．即使是边界条件保持不变，流动也是不稳定的，速度等流动特性都随机变化这种状态称为湍流。图4.1是在湍流状态下在某一点测得的速度随时间的变化情况。可以看出，速度值的脉动性很强。湍流是一种高度复杂的三维非稳态、带旋转的不规则流动。既有大尺度旋涡，又有小尺度旋涡。观测表明，湍流带有旋转流动结构，这就是所谓的湍流涡(turbulenteddies)，简称涡(eddy)。从物理结构上看，可以把湍流看成是由各种不同尺度的涡叠合而成的流动，这些涡的大小及旋转轴的方向分布是随机的。大尺度的涡主要由流动的边界条件所决定，其尺寸可以与流场的大小相比拟，它主要受惯性影响而存在，是引起低频脉动的原因；小尺度的涡主要是由粘性力所决定，其尺寸可能只有流场尺度的干分之一的量级，是引起高频脉动的原因。大尺度的涡破裂后形成小尺度的涡，较小尺度的涡破裂后形成更小尺度的涡。在充分发展的湍流区域内，流体涡的尺寸可在相当宽的范围内连续变化。大尺度的涡不断地从主流获得能量，通过涡间的相互作用，能量逐渐向小尺寸的涡传递。最后由于流体粘性的作用，小尺度的涡不断消失，机械能就转化(或称耗散)为流体的热能。同时由于边界的作用、扰动及速度梯度的作用，新的涡旋又不断产生，这就构成了湍流运动。4.4.2湍流的基本方程一般认为，无论湍流运动多么复杂，非稳态的连续方程和N-S方程对于湍流的瞬时运动仍然是适用的。在此，考虑不可压流动，使用笛卡尔坐标系，速度矢量u在x、y和z方向的分量为u、v和w，写出湍流瞬时控制方程如下：为了考察脉动的影响，目前广泛采用的方法是时间平均法，即把湍流运动看作由两个流动叠加而成，一是时间平均流动，二是瞬时脉动流动。时均值：脉动值：瞬时值：现引入Re平均法，任一变量φ的时间平均值定义为：到目前为止，我们一直假定流体密度为常数，但实际流动中，密度可能是变化的。Bradshaw等指出，细微的密度变动并不对流动造成明显影响，在此，忽略密度脉动的影响，但考虑平均密度的变化，写出可压湍流平均流动的控制方程如下（注意，为方便起见，除脉动值的时均值外，下式中去掉了表示时均值的上划线符号“一”)：方程(4．9)是时均形式的连续方程，方程(4．10)是时均形式的N－S方程。由于在式(4.3）中采用的是Re平均法，因此，方程（4.10）被称为Re时均N-S方程(简称RANS)，常直接称为Re方程。方程(4.11)是标量φ的时均输运方程。上面三式就是用张量的指标形式表示的时均连续方程、Re方程和标量的时均输运方程。这里的i和j指标取值范围是(1，2，3)。可以看到，时均流动的方程里多出与有关的项，我们定义该项为Re应力连续性方程动量方程（张量形式）其他变量方程这里，τij实际对应6个不向的Re应力项，即3个正应力和3个切应力。由式(4．12)、(4．13)和(4．14)构成的方程组共有5个方程(Re方程实际是3个)，现在新增了6个Re应力，再加上原来的5个时均未知量(u,v,w,p和φ)，总共有10个未知量，因此，方程组不封闭，必须引入新的湍流模型(方程)才能使方程组(4.12)、(4.13)和(4.14)封闭。脉动值乘积的时均值：湍流应力Reynolds应力湍流模型：把湍流的脉动值附加项与时均值联系起来的特定关系式Reynolds应力方程法湍流粘性系数法4.5湍流的数值模拟方法简介4.5.1三维湍流数值模拟方法的分类目前的湍流数值模拟方法可以分为直接数值模拟方法和非直接数值模拟方法。所谓直接数值模拟方法是指直接求解瞬时湍流控制方程(4.1)和(4.2)。而非直接数值模拟方法就是不直接计算湍流的脉动特性，而是设法对湍流作某种程度的近似和简化处理，例如，采用4．1节给出的时均性质的Re方程就是其中一种典型作法。依赖所果采用的近似和简化方法不同，非直接数值模拟方法分为大涡模