化工传递过程基础2课件.ppt

1、化工传递过程基础化工传递过程基础 第五章第五章 边界层流动边界层流动 NS方程式反映了流体流动规律，但其解只在某些特殊情况下才能获得，方程式反映了流体流动规律，但其解只在某些特殊情况下才能获得，对很小对很小Re的爬流结果正确，而对的爬流结果正确，而对Re很大的势流导致错误的结果，对此很大的势流导致错误的结果，对此1904年年Prandtl提出边界层学说后才得以解释。提出边界层学说后才得以解释。y u0 第一节第一节 边界层的概念边界层的概念 1、流动现象流动现象 当流体遇到壁面时，由于流体当流体遇到壁面时，由于流体内部粘性力的作用，流速将从壁面处的内部粘性力的作用，流速将从壁面处的0逐渐逐渐

2、ux 增加到增加到u0。即在整个流层中，沿垂直于流动方。即在整个流层中，沿垂直于流动方向产生了速度梯度。向产生了速度梯度。2、提出论点提出论点 Prandtl提出的论点是：假定提出的论点是：假定 ux速度梯度全部集中在紧靠壁面的一薄层流体中，速度梯度全部集中在紧靠壁面的一薄层流体中，该薄层称为边界层，在边界层以外流速不再变化。为此将流动划分为两个区域：该薄层称为边界层，在边界层以外流速不再变化。为此将流动划分为两个区域：边界层（粘性效应起作用，存在明显速度梯度的区域）和主流区。边界层（粘性效应起作用，存在明显速度梯度的区域）和主流区。3、应用应用 边界层理论为许多试验所证实，一些复杂的传递现象

3、可获得解决。边界层理论为许多试验所证实，一些复杂的传递现象可获得解决。4、边界层的形成和发展边界层的形成和发展 形成形成：壁面的粘附作用壁面的粘附作用；流体具有粘性流体具有粘性。发展发展：边界层在一定距离内变化，然后趋于稳定。：边界层在一定距离内变化，然后趋于稳定。99.00uux 在发展过程，边界层内的流动可能由层流转化为湍流，即由层流边界层转在发展过程，边界层内的流动可能由层流转化为湍流，即由层流边界层转为湍流边界层，但在靠近壁面处仍然存在一层层流内层。开始转变的距离称为湍流边界层，但在靠近壁面处仍然存在一层层流内层。开始转变的距离称为为临界距离临界距离xc，转变点取决于临界，转变点取决于

4、临界Rec=5105 。u0 y u0 xc u0 ux 层流边界层层流边界层 过渡区过渡区 湍流边界层湍流边界层 x 在管内流动时，管内壁面形成边界层，而且逐渐加厚，在离进口某一段距在管内流动时，管内壁面形成边界层，而且逐渐加厚，在离进口某一段距离离Le处边界层在管中心汇合，此后的流动称为充分发展了的流动。从管入口处边界层在管中心汇合，此后的流动称为充分发展了的流动。从管入口到汇合处的距离称为到汇合处的距离称为进口段长度进口段长度，以，以Le表示，用于流体物理量的测量时，要表示，用于流体物理量的测量时，要求测点超过求测点超过Le才结果准确。层流时才结果准确。层流时Le=0.05dRe；湍流时

5、；湍流时Le50d。u0 umax 湍流核心湍流核心 Le Le 5、边界层厚度的定义、边界层厚度的定义 一般取流速达到一般取流速达到u0 的的99%处距离壁面的垂直距离（处距离壁面的垂直距离（y方向）为边界层厚度方向）为边界层厚度，即：，即：虽然很小，但对流体的流动阻力，传热、虽然很小，但对流体的流动阻力，传热、传质过程的速率有重要影响，其大小与传质过程的速率有重要影响，其大小与 流体流动时的湍动程度有关。流体流动时的湍动程度有关。第二节第二节 Prandtl边界层方程式边界层方程式 不可压缩流体沿壁面作稳态（层流边界层）流动时，可看作二维流动过程，不可压缩流体沿壁面作稳态（层流边界层）流动

6、时，可看作二维流动过程，若流动方向若流动方向x，与壁面垂直方向，与壁面垂直方向y，则，则NaverStokes方程式及连续性方程式为：方程式及连续性方程式为：%990uuyx)(12222yuxuxpyuuxuuxxxyxx)(12222yuxuypyuuxuuyyyyyx0yuxuyx1、Prandtl边界层方程式的推导边界层方程式的推导 采用采用数量级分析法数量级分析法：当流体流动的：当流体流动的Re很大时，很大时，x，甚至可以忽略不计。，甚至可以忽略不计。因此对式中各项进行数量级分析，使方程式简化。因此对式中各项进行数量级分析，使方程式简化。（采用（采用O代表数量级）代表数量级）（1）取

7、）取x为距离的标准数量级，用为距离的标准数量级，用O（1）表示，记）表示，记 x=O（1）；）；（2）取）取u0为速度的标准数量级，用为速度的标准数量级，用O（1）表示，记）表示，记 u0=O（1）及）及 ux=O（1）；（3）取）取的数量级为的数量级为O（），记），记 =O（）及）及y=O（）；（4）由二维连续性方程式）由二维连续性方程式 知：知：（5）其余数量级：）其余数量级：0yuxuyx)1(Oxux)1(Oyuy)(Ouy)1(Oyux)1(222Oyux)(22Oxuy)1(22Oyuy)1(22Oxux)(Oxuy根据以上讨论，对根据以上讨论，对NaverStokes方程式中各项

8、数量级之间的关系标注为：方程式中各项数量级之间的关系标注为：（1）（）（1）（）（）（）（1/）（1）（2）（1）（1/2）由于：由于：因此方程式简化为：因此方程式简化为：同理：同理：（1）（）（）（）（）（1）（）（2）（）（）（1/）由此数量级分析可得到的结论是：由此数量级分析可得到的结论是：第二个方程式与第一个方程式相比，可以略去；第二个方程式与第一个方程式相比，可以略去；)(12222yuxuxpyuuxuuxxxyxx2222yuxuxx221yuxpyuuxuuxxyxx)(12222yuxuypyuuxuuyyyyyx0)(Oypdxdpxp0dxdp 因此根据数量级分析得出的因

9、此根据数量级分析得出的 Prandtl边界层方程式为：边界层方程式为：以及连续性方程式以及连续性方程式：满足的边界条件：满足的边界条件：y=0，ux=0，uy=0；y=（），），ux=u02、Prandtl边界层方程式的数学解边界层方程式的数学解 将将 代入到边界层方程式得：代入到边界层方程式得：Blasuis采用相似变换法将其转变为常微分方程，进行积分求解。采用相似变换法将其转变为常微分方程，进行积分求解。（1）寻找变量）寻找变量 通过相似变换通过相似变换 用无因次变量代替用无因次变量代替x、y：22yuyuuxuuxxyxx0yuxuyxyuxxuy33222yyxyxyxuyyx0),(

10、xuf0)(过程：过程：通过因次分析，引入变量通过因次分析，引入变量 经分析经分析以质量以质量M、时间、时间及及 x、y、z方向上的长度方向上的长度Lx，Ly，Lz为基本因次，代入：为基本因次，代入：根据因次一致性原则，解得：根据因次一致性原则，解得：),(yxedcbaxyxBuu0ezyxdzxycybxaxxLLLMLLMLLLLL)()()(ca211cb21cd21ce21cccccxyxBuu2121212110即：即：式中：式中：引入流函数引入流函数，找出，找出与与 的关系：的关系：)()()(2102100FxuyFxuyBuucxxxyxuyyxRe),(0),(yx)(00

11、FuxuyyuxxuF0)()()(00fxudFxuxuf0)(（2）引入变量）引入变量 和和，对各项进行变换：，对各项进行变换：),(yxfuxufxuyyux000)(2100ffxufxuxxuyyxfxuxfufuxxux 2022002)(22002200)(yfxuuyfufuyyux 332033000022)(yfxuyfxuufxuuyyux （3）代入到）代入到得：得：（4）解方程式：）解方程式：Blasuis应用级数衔接法，在应用级数衔接法，在=0附近按附近按Taler级数将级数将f（）展）展开，方程的边界条件为：开，方程的边界条件为：22yuyuuxuuxxyxx02

12、 ff f0,0,000uufxuyyx0,0,000fxuufxuyyy0,000uuuufyx在在=0附近按附近按Taler级数将级数将f（）展开：）展开：由边界条件由边界条件：y=0，=0，f（0）=0，c0=0 由边界条件由边界条件：y=0，=0，f（0）=0，c1=0 代入并且整理：代入并且整理：55443322105432ccccccf45342321432cccccf 46352432432cccccf 58475625435432ccccccf041)42(31)24(21)2(224234273632252243ccccccccccc为使上式成立，各项系数等于零，即：为使上式

13、成立，各项系数等于零，即：c3=0 ，c4=0 ，c6=0 ，c7=0 ，式中：式中：A0=1，A1=1，A2=11，c2 由由时的边界条件确定，其求时的边界条件确定，其求解结果为：解结果为：实际计算时可通过查取表实际计算时可通过查取表4-1进行。进行。22521cc231283252222)23()21(814115122nnnnncAcccf542105943.416603.0f11886104277.1104972.23、Prandtl边界层方程式的应用边界层方程式的应用（1）边界层中的速度分布边界层中的速度分布ux，uy：（2）边界层厚度边界层厚度：（3）曳力系数曳力系数CD：设平壁宽

14、度：设平壁宽度b，长度，长度L，流体受到的总阻力为：，流体受到的总阻力为：xyxyuuffxuuRe083.0)(212200000000Re332.0332.0332.0uxyuxuxyufuuxxxuuux005,5,99.05.00Re55xxuxbdxyubdxAuCFLyxLsxDd0002021其中：其中：第三节第三节 Karman边界层积分动量方程式边界层积分动量方程式1、Karman边界层积分动量方程式的推导边界层积分动量方程式的推导 方法：对方法：对Prandtl边界层方程从边界层方程从y=0到到y=进行积分，然后根据速度分布求解。进行积分，然后根据速度分布求解。xuufxu

15、uyuyxsx020000332.0)0(5.0200020Re664.0332.0LLdbLuxdxbuuF5.0205.020Re328.1Re664.02LLDbLubLuCPrandtl边界层方程式左侧积分：边界层方程式左侧积分：其中：其中：000000)(yxyxxxxyxxxyxxduuuudyxuuduudyxuudyyuuxuudyyuudyxuuyxxx00dyxuudyyuuuuuuxyyxyx0000000dyxuudyyuuduuxxyxyx000dyuuuxdyxuudyxuudyyuuxuuxxxxxxyxx)(2)(000000Prandtl边界层方程式右侧积分：

16、边界层方程式右侧积分：因此因此 Karman边界层积分动量方程式：边界层积分动量方程式：若已知若已知uxy的关系，通过对的关系，通过对Karman边界层动量方程式积分，可得速度分布等。边界层动量方程式积分，可得速度分布等。2、流体沿平版壁面流动时层流边界层的近似解、流体沿平版壁面流动时层流边界层的近似解（1）速度分布速度分布：不可压缩流体作稳态二维流动时，根据实验测定层流边界层内：不可压缩流体作稳态二维流动时，根据实验测定层流边界层内速度分布与抛物线形状相似，即：速度分布与抛物线形状相似，即：其中系数其中系数ai由相应的边界条件确定，见由相应的边界条件确定，见87-89页。页。ssyxyxxy

17、uyudyyu)0()()(0022sxxdyuuudxd)(00iniixyau0设速度分布方程式为：设速度分布方程式为：根据边界条件：根据边界条件：得层流边界层内速度分布方程式：得层流边界层内速度分布方程式：32dycybyaux0,0,0auyx00,0,022cyuuuyxyx30,dbuuyx0,30dbyuuuyxx320ub 302ud30)(21)(23yyuux（2）边界层厚度边界层厚度：将边界层内速度分布方程代入：将边界层内速度分布方程代入Karman边界层动量方程式中边界层动量方程式中当当 x=0 时，时，=0，故，故 c1=0sxxdxdudyyyyyudxddyuuu

18、dxd28039)2231)(223()(2033330200023)2323(0023000uyuudyduyyxs2328039020udxdudxud0131401064.4cux5.00Re64.464.4xxux（3）曳力系数曳力系数CD：设平壁宽度：设平壁宽度b，长度，长度L，流体受到的总阻力为：，流体受到的总阻力为：其它情况下的速度分布、边界层厚度、曳力系数见表其它情况下的速度分布、边界层厚度、曳力系数见表4-2中。中。bdxyubdxAuCFLyxLsxDd00020213000646.023LubbdxuFLd5.002030Re292.1292.1646.02LDLubLu

19、LubC 第四节第四节 边界层分离边界层分离 当流体绕过圆柱或球体等流动时，当流体绕过圆柱或球体等流动时，Re很小时阻力由粘性力引起；很小时阻力由粘性力引起；Re较大时摩擦阻力和形体阻力都有影响，而形体阻力取决于边界层分离。较大时摩擦阻力和形体阻力都有影响，而形体阻力取决于边界层分离。1、现象分析现象分析 流体流过平行置于流场中的薄平板时，沿流动方向边界层外流体流过平行置于流场中的薄平板时，沿流动方向边界层外的速度、压力保持不变，即的速度、压力保持不变，即dp/dx=0；但当流过曲面时，边界层外的流速、；但当流过曲面时，边界层外的流速、压力沿流动方向发生不断变化，由压力沿流动方向发生不断变化，

20、由Benulii方程式：方程式：2、结论结论 对边界层外的对边界层外的加速加速过程，边界层内外为过程，边界层内外为减压减压过程，压力梯度为负；过程，压力梯度为负；而对边界层外的而对边界层外的减速减速过程，边界层内外均为过程，边界层内外均为加压加压过程，压力梯度为正。过程，压力梯度为正。3、影响影响 流体流过曲面时，夹在主流和固体表面间的边界层，在加速减压流体流过曲面时，夹在主流和固体表面间的边界层，在加速减压阶段，虽受到粘性力的作用而减小，但仍能向下游流动；而在减速加压阶阶段，虽受到粘性力的作用而减小，但仍能向下游流动；而在减速加压阶段，同时受到粘性力和逆向压力的作用，紧贴壁面的流体速度迅速下

21、降，段，同时受到粘性力和逆向压力的作用，紧贴壁面的流体速度迅速下降，当到达当到达S点时所有的动能耗尽，出现点时所有的动能耗尽，出现停滞停滞。但后面的流体继续流动，在惯。但后面的流体继续流动，在惯性力的作用下，使性力的作用下，使边界层流体脱离了固体壁面边界层流体脱离了固体壁面，该现象，该现象称为边界层分离称为边界层分离。Kup2021xuuxp00边界层开始与固体表面分离的点边界层开始与固体表面分离的点S称为称为分离点分离点，其上，其上4、边界层分离的结果边界层分离的结果 产生倒流和大量旋涡，形成极不规则的湍流区，产生倒流和大量旋涡，形成极不规则的湍流区，使得使得能量损失急剧加大能量损失急剧加大

22、。5、形成边界层分离的必要条件、形成边界层分离的必要条件 流体具有粘性流体具有粘性；存在逆向压力梯度存在逆向压力梯度。边界层。边界层分离是形成旋涡的重要来源，旋涡导致形体阻力，为产生局部阻力的主要原分离是形成旋涡的重要来源，旋涡导致形体阻力，为产生局部阻力的主要原因。因。6、应用、应用 用于计算局部阻力，工程上为减小阻力采取相应措施。用于计算局部阻力，工程上为减小阻力采取相应措施。u0 加速减压加速减压 减速减速 加压加压 00yxyus 第六章第六章 湍流湍流 湍流是指湍流是指Re4000（圆形直管内）的流动，质点间碰撞混合程度剧烈，（圆形直管内）的流动，质点间碰撞混合程度剧烈，阻力要大于层

23、流。研究湍流的内容是：导致发生原因，特征，流动规律。阻力要大于层流。研究湍流的内容是：导致发生原因，特征，流动规律。第一节第一节 湍流的特点、形成、表征湍流的特点、形成、表征一、湍流的特点一、湍流的特点 湍流是在高湍流是在高Re数下发生的流动过程，数下发生的流动过程，特点特点流体向前流动时流体向前流动时伴随不规则的脉动，混合剧烈，流动参数随时发生变化。伴随不规则的脉动，混合剧烈，流动参数随时发生变化。其其基本特征是质点的脉动基本特征是质点的脉动。脉动的结果导致：脉动的结果导致：流动阻力加大流动阻力加大；速度分布均匀速度分布均匀 （但在近壁处存在层流内层）。（但在近壁处存在层流内层）。二、湍流的

24、形成二、湍流的形成 形成湍流具备的条件：形成湍流具备的条件：旋涡的形成旋涡的形成；旋涡的运动旋涡的运动。1、旋涡的形成、旋涡的形成 （1）流体具有粘性流体具有粘性，相邻流层间构成力偶，是产生旋涡的基，相邻流层间构成力偶，是产生旋涡的基本因素；本因素；（2）流层的波动流层的波动（或产生边界层分离），在横向压力和剪应（或产生边界层分离），在横向压力和剪应力的双重作用下导致了旋涡的形成。力的双重作用下导致了旋涡的形成。2、旋涡的运动、旋涡的运动 由于旋涡的形成，使附近流层的速度分布改变，由于旋涡的形成，使附近流层的速度分布改变，产生了压力差，促使旋涡脱离原来的流层进入邻近产生了压力差，促使旋涡脱离原

25、来的流层进入邻近的流层，各流层间旋涡的不断交换形成了旋涡。的流层，各流层间旋涡的不断交换形成了旋涡。三、湍流的表征三、湍流的表征1、时均量、脉动量和瞬时量、时均量、脉动量和瞬时量 ux 湍流中质点的运动极不规则，为非稳定流动，采湍流中质点的运动极不规则，为非稳定流动，采用统计方法或取平均值的方法进行处理。用测速仪用统计方法或取平均值的方法进行处理。用测速仪测出某段时间内流体测出某段时间内流体瞬时速度瞬时速度ux随时间变化关系如图，随时间变化关系如图，ux随时间虽变化频繁，但总是围绕随时间虽变化频繁，但总是围绕“平均值平均值”在波动。在波动。0 （1）时均量）时均量 取取0内内ux的时间平均值的

26、时间平均值，称为，称为时均速度时均速度：01duuxx（2）脉动量）脉动量 实际速度和时均速度之差称为实际速度和时均速度之差称为脉动速度脉动速度 （其值可正可负）：（其值可正可负）：且：且：（3）瞬时量）瞬时量 瞬时速度等于时均速度与脉动速度之和。瞬时速度等于时均速度与脉动速度之和。区别：区别：瞬时量瞬时量指某时刻运动参数的真实值；指某时刻运动参数的真实值；时均量时均量指某时段内瞬时量的平均值；指某时段内瞬时量的平均值；脉动量脉动量指某时刻运动参数的真实值与时均值的差值（可正可负）。指某时刻运动参数的真实值与时均值的差值（可正可负）。2、湍动强度（湍流的激烈程度）、湍动强度（湍流的激烈程度）湍

27、动强度湍动强度I=脉动速度脉动速度/时均速度时均速度 用用 代替代替 则：则：xuxxxuuuxxxuuu_2xu _xu0_xuxzyxuuuuI)(31_2_2_2 第二节第二节 流体湍流时的运动方程式流体湍流时的运动方程式 引入瞬时速度等于时均速度与脉动速度之和，且各脉动速度的时均值引入瞬时速度等于时均速度与脉动速度之和，且各脉动速度的时均值为零，可将流体的湍流流动理解为按时均速度在流动，使得问题简化。但为零，可将流体的湍流流动理解为按时均速度在流动，使得问题简化。但因湍流的本质是因湍流的本质是质点的脉动，因此必须考虑脉动。质点的脉动，因此必须考虑脉动。Reynold将将瞬时速度等瞬时速

28、度等于时均速度与脉动速度的方程代入到以应力表示的运动微分方程式中，然于时均速度与脉动速度的方程代入到以应力表示的运动微分方程式中，然后取时均值，导出相应的湍流运动方程式，过程称为雷诺转换。后取时均值，导出相应的湍流运动方程式，过程称为雷诺转换。一、一、Reynold方程式方程式1、时均值的有关运算法则：、时均值的有关运算法则：设设 f1 和和 f2 代表湍流运动代表湍流运动 时的两个物理量，而且：时的两个物理量，而且：则有：则有：（1）（2）（3）（4）111fff222fff11ff 21_21ffff21_21ffff0_1f（5）（6）2、对不可压缩流体的连续性方程式进行雷诺转换、对不可

29、压缩流体的连续性方程式进行雷诺转换(2、6）：）：即湍流时的时均速度仍然满足连续性方程式。即湍流时的时均速度仍然满足连续性方程式。_2121_21ffffffxfxf1_1yfyf1_1zfzf1_10zuyuxuzyx0_zuyuxuzyx0zuyuxuzyx3、对以应力表示的运动微分方程式进行雷诺转换（、对以应力表示的运动微分方程式进行雷诺转换（x方向）：方向）：由于湍流时包括脉动量，对两侧各项时均化，运用法则（由于湍流时包括脉动量，对两侧各项时均化，运用法则（2）、（）、（6）、（）、（5）得：）得：zyxXzuuyuuxuuzxyxxxxzxyxx)(0)(zuyuxuuzyxxzyx

30、Xzuuyuuxuzxyxxxxzxyx)()(2zyxXzuuyuuxuzxyxxxxzxyx)(_2zyxXzuuyuuxuzuuyuuxuzxyxxxxzxyxxzxyx)()()(_22将含脉动量的各项移到右侧，展开左侧第一项，得：将含脉动量的各项移到右侧，展开左侧第一项，得：左侧第一项左侧第一项 而且：而且：即为不可压缩流体稳态湍流时的时均运动方程式（即为不可压缩流体稳态湍流时的时均运动方程式（x方向）方向），称为雷诺方，称为雷诺方程式。在程式。在y、z方向可得到类似的方程式。方向可得到类似的方程式。)()(zuuyuuxuuzuuyuuxuuzxyxxxxzxyxx)(zuuyuu

31、xuuDuDxzxyxxx)()()(_2zxzxyxyxxxxxuuzuuyuxXDuD0)(zuuyuuxuuzxyxxx二、雷诺应力二、雷诺应力 上述方程式多出上述方程式多出3项，因此可推知；湍流时所产生的应力除和项，因此可推知；湍流时所产生的应力除和层流相同的部分外，还存在一部分附加应力。即一个法向附加应层流相同的部分外，还存在一部分附加应力。即一个法向附加应力力 和两个切向附加应力和两个切向附加应力 ，称为雷诺应力或表观应力。湍流时雷诺应力较粘性应力大得多。称为雷诺应力或表观应力。湍流时雷诺应力较粘性应力大得多。在在x方向方向的雷诺应力，的雷诺应力，总的时均应力总的时均应力可表示为：

32、可表示为：rzxzxtzxxzrzxryxyxtyxyxryxrxxxxtxxxrxxuuuuu_2_2xu_zxyxuuuu而在三维流动时，诸雷诺应力的应力矩阵表示为：而在三维流动时，诸雷诺应力的应力矩阵表示为：由上面看出，一般由上面看出，一般雷诺应力前均加一个负号雷诺应力前均加一个负号，为什么？，为什么？分析获得分析获得：在层流内层，仅粘性应力起作用，雷诺应力不存在；在层流内层，仅粘性应力起作用，雷诺应力不存在；在湍流区，主要雷诺应力起作用，粘性应力很小；在湍流区，主要雷诺应力起作用，粘性应力很小；在过渡层，粘性应力和雷诺应力同时起作用。在过渡层，粘性应力和雷诺应力同时起作用。rzzryz

33、rxzrzyryyrxyrzxryxrxx 第三节第三节 涡流粘度与涡流粘度与Plandtl混合长混合长一、湍流应力一、湍流应力 1877年年Boussinesq提出假设，类似于粘性应力，雷诺应力可表示为：提出假设，类似于粘性应力，雷诺应力可表示为：二、二、Plandtl混合长混合长 1925年，年，Plandtl据层流和湍流之间动量传递机理据层流和湍流之间动量传递机理的类似性，将分子动量传递过程中平均自由程的概念的类似性，将分子动量传递过程中平均自由程的概念用于湍流，提出了混合长的假设。即：用于湍流，提出了混合长的假设。即：脉动过程流体脉动过程流体微团保持原微团保持原x方向时均流速（动量）不

34、变时的脉动垂直方向时均流速（动量）不变时的脉动垂直距离距离，称，称混合长。混合长。假定假定混合长足够小，则：混合长足够小，则：若由下向上脉动：若由下向上脉动：，若由上向下脉动：若由上向下脉动：dyudxryx)(_yxuu xuxu2xu1xudyudluuxxx12dyudluxxdyudluxx根据质量守恒定律，根据质量守恒定律，y方向的脉动必引起方向的脉动必引起x方向的脉动，假定：方向的脉动，假定：则：则：故雷诺应力故雷诺应力比较可得：比较可得：第四节第四节 圆管中的稳态湍流流动圆管中的稳态湍流流动（注：在以后的讨论中将上下标略去，表示一维流动，速度均指时均速度）（注：在以后的讨论中将上

35、下标略去，表示一维流动，速度均指时均速度）一、通用速度分布方程式一、通用速度分布方程式（x方向）方向）xyxyuuuu,22)(dyudluuxxyryx_yxuu 22)(dyudlxdyudxdyudlx21、层流内层、层流内层 令：令：常数，在常数，在 0y范围内积分：范围内积分：采用无因次形式表达时，采用无因次形式表达时，令令 称为摩擦速度（称为摩擦速度（m/s）分别称为无因次速度、无因次距离。分别称为无因次速度、无因次距离。因而：因而：即为即为层流内层层流内层通通 用速度分布方程用速度分布方程 式式：dydusyusyuusyuyuuu,yuyyuuusss yu2、湍流中心、湍流中

36、心Plandtl根据混合长的学说，假设：根据混合长的学说，假设：常数；在管径范围内常数；在管径范围内因此：因此：积分：积分：即为即为湍流中心湍流中心通用速度分布方程式通用速度分布方程式：22)(dydulSKyl dyduKyudyduKydyduyKss,)(,)(2222CnyuKu12111111CnyKCunKyunKCnyKuu21CnyKu Nikurade采用实验方法在半对数坐标上对上述关系进行了描绘，得到：采用实验方法在半对数坐标上对上述关系进行了描绘，得到：所以：所以：由实验结果（图由实验结果（图5-9）可看到，当）可看到，当 与实验结果相吻合；与实验结果相吻合；当当 也与实

37、验结果相吻合；但当也与实验结果相吻合；但当 时两式均不适时两式均不适用，其原因是存在一个过渡区。用，其原因是存在一个过渡区。Karmen等建议在等建议在Re=40003200000范围内，将范围内，将光滑管中的流动划分为三个区域：光滑管中的流动划分为三个区域：（1）层流区：）层流区：（2）过渡区：）过渡区：（3）湍流区：）湍流区：5.5,4.02CK5.55.2nyuyuy55.55.230nyuy305yyuy505.55.230nyuy05.30.5305nyuy由此可以得到：由此可以得到：（1）层流内层厚度：）层流内层厚度：过渡区厚度：过渡区厚度：湍流中心厚度：湍流中心厚度：（2）圆管中

38、心流速）圆管中心流速 或者：或者：uuyubb5,5uuuyuym25530,30)(2mbcR)5.55.2(maxRunuuyRnuuu5.2maxnRyuu1max)(二、光滑管中的平均速度和流动阻力二、光滑管中的平均速度和流动阻力根据根据Fanning摩擦系数的定义：摩擦系数的定义：代入：代入：即为光滑管中即为光滑管中Fanning摩擦系数摩擦系数 f 的半经验公式，其它的有的半经验公式，其它的有Blasius方程式方程式等，注意此处：等，注意此处：)75.15.2()5.55.2(211022RunurdryunuRddrurRudAAuRAAbfufuufsbbs2222)75.1

39、5.2(2Runufu601.0)(Re768.11fnf41f5/1Re046.0fRe16f4/1Re079.0f三、粗糙管中的速度分布和流动阻力三、粗糙管中的速度分布和流动阻力Nikurade对粗糙管（对粗糙管（e/d0）的流动阻力进行了研究，得出的结论是：）的流动阻力进行了研究，得出的结论是：1、层流区：、层流区：e/d对流动阻力无影响；对流动阻力无影响；2、过渡区：、过渡区：e/d对流动阻力基本无关；对流动阻力基本无关；3、湍流区：从某一、湍流区：从某一Re开始，开始，e/d对流动阻力产生影响；甚至当对流动阻力产生影响；甚至当Re达到一定程度达到一定程度后，流动阻力完全取决于后，流动阻力完全取决于e/d而与而与 Re无关，因而将湍流区分为三种状态。无关，因而将湍流区分为三种状态。（1）水力光滑状态：）水力光滑状态：（2）过渡状态：）过渡状态：（3）完全粗糙状态：）完全粗糙状态：(Re),50ffeu)(Re,705deffeu)Re51.27.3(21deg3116)Re1020000(10055.0de)(,70deffdeu66.3277.11ednf