高等工程流体力学

1、高等工程流体力学粘性流动康 顺华北电力大学能源与动力工程系学院Kangs精选文库内容提纲? 边界层及其方程? 层流边界层流动转捩? 湍流边界层结构? 流动分离、二次流动与旋涡能源动力领域流动问题的主要特征? 全三维? 非定常? 粘性高雷诺数，边界层边界层：层流、转捩、湍流（紊流），分离流动，旋涡运动叶轮机械（透平和压气机等）大多由单个或多个级组成。每个级含有一排静子叶片列和一排转子叶片列。在级内的气流场中，一般至少有以下几种流动现象发生： 1、前缘马蹄涡； 2、通道涡； 3、顶部间隙涡； 4、边界层转捩； 5、叶片尾迹； 6、旋涡、尾迹等与叶片列周期性非定常相互作用。激波、激波与边界层相互作用

2、-2精选文库边界层流动边界层边界层概念：粘性很小的流体以大雷诺数运动时，在大部分流场上可以略去粘性的作用；但在物面附近的很薄的一层流体内必须考虑粘性作用。这一薄层流体称为边界层。平板边界层示意图有边界的流动图谱如右上图所示：流动分为三个区：边界层，尾迹区，位流区（外部势流区）二维平板的边界层微分方程设直匀流以零迎角平行流过一块长度为l 的平板 ,如左下图所示，人为规定，当某个y 处的速度达到层外自由流的99时，这一点到物体表面的距离（即y）称为边界层在改点的厚度，记为。显然，边界层的厚度是与X 有关的，所以可以写成（ x） 。平板边界层边界层的厚度（x）很小，满足此关系式：（x）<<

3、;l在忽略质量力的前提下，粘性平面不可压流的运动方程加上连续方程是:uuu1p2u2uuxy(x2y2 )tx1p22u2 )xy(x2ytyu0xy用边界层条件式（x）<<l上式， y 的数值限制在边界层之内，即0 y-3精选文库经过数量级分析，上面方程组化为：uu uu1 p2utxyxy2p0yu0xyp0 的物理意义： 在边界层内， 沿物体表面的发法线方向压强是不变的，亦即等于边y界层处自由流的压强。? 卡门动量积分关系解采用动量积分法得出控制面ABCD 的动量变化：dtdx( d2 dyddy)dx 0dx 0其中：为边界层边界上的流速。dp作用在 AB,BC,CD,AD

4、四个面上的力在 x 方向上投影的合力的冲量是：(w )dt dxdx根据动量定理得： d2 dyddy (dpw )dx 0dx 0dx即定常流动的边界层动量积分关系式，也叫卡门波尔豪森（Karman-Pohlhausen）动量积分关系式。该式不仅适用于层流边界层，也适用于准定常紊流边界层；不仅适用于平板， 也适用于微弯曲面；既适用于不可压流，也适用于可压流。引入两个新的概念：一个叫位移厚度，另一个叫动量损失厚度。分别定义如下：0(1-) dy0平板边界层示意图(1-) dy的物理意义：相当于流线在边界层边界处被推出边界层外的距离，为了把粘性考虑进去，就要把按理想流算出的流型各点的y 坐标都增

5、加当地位移厚度那么大的尺寸。的物理意义：2(2 )dy0-4精选文库等号右侧第一项是实际流量乘以层外流速这样一个假想动量，而第二项则是实际流量乘以实际流速，这是实际动量。二者之差就是层内那部分流量在没有粘性力作用时应有的动量与粘性力作用时的实际动量之差，也就是由于有粘性力作用而损失的动量，这些损失的动量折合成以流动的，厚度为的一层流体所具有的动量，动量损失厚度即为：202dy(1-) dy0然后，引用符号（形状因子）H/，边界层动量积分关系式可以改写为：d''d(2 H)wdx2其中dx采用控制体方法建立积分关系式在定常流动的流体中，沿边界层划出一个单位宽度的微小控制体，它的投

6、影面ABDC 。由作为 x 轴的物体壁面上的一微元距离 BD 、边界层的外边界 AC 和彼此相距 dx 的两直线 AB和 CD 所围成。应用动量方程来研究该控制体内的流体在单位时间内沿x 方向的动量变化和外力之间的关系。边界层厚度于雷诺数成反比，与流程的距离成正比：三种边界层厚度的比较*对于平板层流边界层：边界层厚度正比于 x1/25.0x*1.72xRexRex对于平板湍流边界层：边界层厚度正比于 x4/50.4x1Rex5湍流边界层厚度大于层流边界层厚度摩擦阻力系数与雷诺数成反比对于平板层流边界层：壁面摩察应力Cf 反比于 L1/2对于平板湍流边界层：壁面摩察应力Cf 反比于 L1/5边界

7、层的动量积分0.644xRexC f1.328ReLC f0.0741Re5L-5精选文库湍流壁面摩察阻力大于层流壁面摩察阻力流体在翼型的上表面形成较大的扩压区，以致引起边界层的分离。随着冲角的增大，分离点向前移动，在翼上表面的大部分及其后面形成很大的尾涡区，使翼型上、下表面的压差减小，因此升力和升力系数都急剧下降。最大升力系数的相应点称为失速点，因为超过该点飞机和涡轮机的性能均将恶化。再看翼型的阻力系数Cd 与冲角的关系。当冲角较小时，翼型前后的压差很小，总的阻力中摩擦阻力是主要的，阻力边界层转捩湍流局部换热系数-6精选文库稳定层流 T/S 波 展向旋涡三维涡破碎湍流斑充分发展湍流影响边界层

8、转捩的因素转捩是个复杂的物理过程，受到很多因素的影响：雷诺数自由流湍流强度压力梯度、表面曲率、加速率表面粗糙度、拌线壁面温度及热流噪声波激波与边界层相互干涉等等。为了准确预估转捩，必须综合考虑上述诸多因素，但至今还没有完整的理论。转捩的分类自然转捩（ Nature Transition ）当来流湍流度较低 (小于 0.1%) ，层流中的 Tollmien-Schlichting波或横流不稳定波的非线性指数增长将导致湍斑的出现，之后流动迅速发展为完全湍流状态，这种过程被称为自然转捩或横流转捩旁路（跨越）转捩（Bypass Transition ）当来流湍流度较高时，转捩会直接绕过T-S 波阶段，

9、在来流湍流度的影响下，湍流斑直接生成，而后迅速发展成为完全湍流状态，由于这种转捩跨越了自然转捩的一个或两个阶段，故被称为跨越转逆向转捩（ Reverse Transition ）当加速率足够大时发生。KdUU 2 dx-7精选文库分离泡转捩（ Separation Bubble Transition ）物体扰流的前部分边界层流动常处于层流状态，其抵抗逆压梯度的能力弱，容易产生分离，层流分离后，剪切层离开壁面， 由于其速度剖面不稳定， 很快导致转捩，转捩后的湍流裹入能量使边界层再附，形成层流分离泡，这类转捩被称为分离泡转捩。周期性非定常转捩（Perodic-Unsteady Transition

10、 ）-8精选文库粘性流动一切流体都具有粘性在高雷诺数条件下，粘性影响主要在固体壁面附近边界层内，以及尾迹区内。其它区域的流动可近似认为无粘。边界层和尾迹区内流动的主要特征：大的速度和总压梯度；产生摩擦阻力和能量损失边界层壁面速度为零、总压最小；尾迹中心线上速度和总压最小，损失最大粘性流动层流流动uy层流边界层内速度分布呈线性：如AU壁面摩擦力为：duwdy y 0m 为流体的粘性系数，一般随温度变化,例如：萨特兰（Sutherland）公式1.5T T 1100T0T0 110U描述层流流动的基本方程为NS（ Navier Stocks）方程tFIFVQU 为矢量守恒；UWxWyWzEFI 和

11、 FV 分别代表无粘矢通量和粘性矢通量，可表达为：FIFIx i FIy j FIz kFVFVx i FVy j FVz kWi0P xiWiWxixFIiP yiWiWyFViiyP ziWiWzizE P WiqiWz iz Wy iy Wz iz1,ijWiWj2TijWijqi0,ijijX i3X iX jNS 方程中的源项Q其中包含哥氏力和离心力-9精选文库002 X x2WyQ2 Wr2 X y2Wx12r20W22X xWxX yWy未知量为：密度r、三个速度分量和温度。压力可由状态方程求出。理想气体的状态方程为pR T六个未知数、六个方程粘性流动无量纲参数雷诺数：ReLUR

12、 为流体的密度， L 为特征长度（一般取为转子叶顶半径） ， U 为特征速度（一般取为转子叶顶的线速度）普朗特数：PrC p对于空气，普朗特数为0.708粘性流动湍流流动层流流动间歇湍流完全发展湍流热线测得的当地速度随时间的变化湍流流动时非定常流动湍流由不同尺度的旋涡形成。湍流的尺度由涡旋的尺度所确定。能量从大尺度涡旋向小尺度涡旋传递大尺度区雷诺数大，小尺度区雷诺数小-10精选文库湍流的时间尺度：k /湍流的速度尺度：uk 1/ 2湍流的长度尺度：lk3 / 2 /湍流雷诺数：Re,tk1 / 2lk2驱动湍流的能量栅湍流能量由平均流动产生，由能量栅传递，并由微尺度级的粘性影响而耗散雷诺平均N

13、S 方程RANS方程假设任何流动参数其时均值其脉动值时均值脉动值i, jl i , jt i, j可压缩形式湍动能：在湍流条件下，静压与总能量的表达式将为：在考虑湍动能时，雷诺应力为：特征长度和速度形湍流边界层的长度尺度：对 NS 方程进行时均处理后RANS 方程新的未知量:雷诺应力张量ui' u'juiu jt i , jtxix j粘性应力层流应力雷诺应力ij( luiu j2t )xiuijx j3粘性底层过渡层对数层yyu摩擦速度；uwuu近似速度分布：uy1u U ( )7y-11精选文库湍流相关量1 (u湍流度：Tu3对于各向同性湍流，u22vUv2 w 2 )2w

14、 2所以Tuu 2U湍动能： k1 (u 2v2w 2 )3湍动能耗散率2(uiu j ) 2x jxi雷诺应力特性平板湍流边界层中雷诺应力分量的变化：各雷诺应力在壁面处均为零，但在离壁面很近的区域内迅速增加，达到最大值后再下降。在边界层外沿处接近于零-12精选文库粘性底层过渡层对数层uuu湍流模型y高雷诺数模型低雷诺数模型不计及粘性底层不应用壁面函数应用壁面函数模拟对数率区所有区域都需要求解不需要精确捕捉对数率层的流动需要捕捉对数率层和粘性底层的流动y 30 50y 1本课所介绍的湍流模型均为低雷诺数模型注意：这里的雷诺数为湍流雷诺数问题：如何求解 t 湍流模型tC U*L*U *湍流参考速

15、度L* 湍流参考长度现有多种湍流模型, 若安附加的方程数分类，可有：代数模型或0方程模型例如 : Cebeci-Smith, Baldwin-Lomax,计算时均短1-方程模型例如 : Spalart-Allmaras,增加一个变量，一个方程。计算时均略长2-方程模型例如 : k-Models,增加两个变量，两个两个方程。计算时间长? 代数模型（ 0 方程模型）Baldwin Lomax 模型（简称 BL模型）两层模型：内层t 由普朗特混合长度模型确定(t )i , nnc外层 t 由平均流和特征长度确定t(t )o , nncn 离开壁面第一次 ( ) ( )的法向坐标ct it o( t )il 2, lkn(1e y / A ), yw w n,iu jijkxkw-13精选文库toKC cp Fwake FKleb ( n)特点：简单易用，可用于附着流动和小的分离流动。往往过度模拟分离流动? 1 方程模型Spalart-Allmaras模型（简称 SA 模型）引入湍流运动变量tt f v1湍流运动变量遵循如下的输运方程：1 cb2 QW l1 cb 2t其它变量的定义见Fine/Turbo 手册特点：加入一个控制方程，与 BL 模型比较，求解速度