8第八章-边界层理论基础和绕流运动
- 格式:doc
- 大小:425.00 KB
- 文档页数:7
第八章 边界层理论§8—1 边界层的基本概念实际流体和理想流体的本质区别就是前者具有粘性。
对层流而言,单位面积摩擦力的大小yud d μτ=,可以看出,对于确定的流体的等温流场,摩擦力的大小与速度梯度有关,其比例函数即动力粘度。
速度梯度yud d 大,粘性力也大,此时的流场称为粘性流场。
若速度梯度yud d 很小,则粘性力可以忽略,称为非粘性流场。
对于非粘性流场,则可按理想流体来处理。
则N-S 方程可由欧拉方程代替,从而使问题大为简化。
Vlv l lV v A y u V l tVl t u mρρμρρ======2223d d d d 粘性力惯性力当空气、蒸汽,水等小粘度的流体与其它物体作高速相对运动时,一般雷诺数很大。
由vVl==粘性力惯性力Re ,则在这些流动中,惯性力〉〉粘性力,所以可略去粘性力。
但在紧靠物体壁面存在一流体薄层,粘性力却与惯性力为同一数量级。
所以,在这一薄层中,两者均不能略去。
这一薄层就叫边界层,或叫速度边界层,由普朗特在1904年发现.a .流体流过固体壁面,紧贴壁面处速度从零迅速增至主流速度,这一流体薄层,就叫边界层或速度边界层。
b .整个流场分为两部分 层外,0=∂∂yu,粘性忽略,无旋流动。
层内,粘性流,主要速度降在此,有旋流动.c .由边界层外边界上∞=V u %99,来定义δ,δ为边界层厚度。
d .按流动状态,边界层又分为层流边界层和紊流边界层。
由于在边界层内,流体在物体表面法线方向(即yu∂∂)速度梯度很大,所以,边界层内的流体具有相当大的旋涡强度;而在层外,由于速度梯度很小。
所以,即使对于粘度很大的流体,粘性力也很小,故可忽略不计,所以可认为,边图8-2空气沿平板边界层速度分布外部区域边界层界层外的流动是无旋的势流.边界层的基本特征有: (1)1<<Lδ⇒薄层性质,其中L 为物体的长度;沿流方向↑↑→δx 。
(2) 层内yu∂∂很大, 边界层内存在层流和紊流两种流态。
第八章 边界层理论§8-1 边界层的基本概念实际流体和理想流体的本质区别就是前者具有粘性。
对层流而言,单位面积摩擦力的大小yud d μτ=,可以看出,对于确定的流体的等温流场,摩擦力的大小与速度梯度有关,其比例函数即动力粘度。
速度梯度yud d 大,粘性力也大,此时的流场称为粘性流场。
若速度梯度yud d 很小,则粘性力可以忽略,称为非粘性流场。
对于非粘性流场,则可按理想流体来处理。
则N-S 方程可由欧拉方程代替,从而使问题大为简化。
Vlv l lV v A y u V l tVl t u mρρμρρ======2223d d d d 粘性力惯性力当空气、蒸汽,水等小粘度的流体与其它物体作高速相对运动时,一般雷诺数很大。
由vVl==粘性力惯性力Re ,则在这些流动中,惯性力>>粘性力,所以可略去粘性力。
但在紧靠物体壁面存在一流体薄层,粘性力却与惯性力为同一数量级。
所以,在这一薄层中,两者均不能略去。
这一薄层就叫边界层,或叫速度边界层,由普朗特在1904年发现。
a .流体流过固体壁面,紧贴壁面处速度从零迅速增至主流速度,这一流体薄层,就叫边界层或速度边界层。
b .整个流场分为两部分 层外,0=∂∂yu,粘性忽略,无旋流动。
层内,粘性流,主要速度降在此,有旋流动。
c .由边界层外边界上∞=V u %99,来定义δ,δ为边界层厚度。
d .按流动状态,边界层又分为层流边界层和紊流边界层。
由于在边界层内,流体在物体表面法线方向(即yu∂∂)速度梯度很大,所以,边界层内的流体具有相当大的旋涡强度;而在层外,由于速度梯度很小。
所以,即使对于粘度很大的流体,粘性力也很小,故可忽略不计,所以可认为,图8-2空气沿平板边界层速度分布外部区域边界层边界层外的流动是无旋的势流。
边界层的基本特征有: (1)1<<Lδ⇒薄层性质,其中L 为物体的长度;沿流方向↑↑→δx 。
(2) 层内yu∂∂很大, 边界层内存在层流和紊流两种流态。
杭州电子科技大学全国硕士研究生招生考试业务课考试大纲考试科目名称:工程流体力学科目代码:882第一章绪论1-1 工程流体力学的学科任务1-2 连续介质假设,流体的主要物理性质1-3 作用在流体上的力1-4 工程流体力学的研究方法第二章流体静力学2-1 流体静压强特性2-2 流体的平衡微分方程及积分式、等压面方程2-3 流体静力学基本方程及物理意义和几何意义,压强的计算单位和表示方法,静压强的分布图、测压计原理2-4 液体的相对平衡2-5 作用在平面上的液体总压力表示方法2-6 作用在曲面上的液体总压力计算,虚、实压力体区别2-7 阿基米德原理,浮力和潜体及浮体的稳定性第三章流体运动学3-1 描述流体运动的两种方法及其特点,迹线、流线、脉线的表示3-2 描述流体运动的一些基本概念3-3 流体运动的类型3-4 流体运动的连续性方程的表示3-5 流体微元运动的基本形式及与速度变化的关系3-6 无涡流和有涡流,速度势和速度环量第四章理想流体动力学和平面势流4-1 理想流体的运动微分方程—欧拉运动微分方程,伯努利方程及其条件4-2 理想流体元流的伯努利方程及其物理、几何意义,皮托管原理4-3 恒定平面势流,速度势和流函数的性质及其两者的关系第五章实际流体动力学基础5-l实际流体的运动微分方程——纳维一斯托克斯方程,流体质点的应力状态及压应力的特性5-2 实际流体元流的伯努利方程及其物理、几何意义5-3 实际流体总流的伯努利方程及应用条件,文丘里管工作原理,有能量输入和输出的伯努利方程5-5 总流的动量方程及其应用条件和方法第六章量纲分析和相似原理6-l 量纲分析,量纲和单位,量纲和谐原理种类和区别6-2 流动相似原理6-3 相似准则6-5 模型试验第七章流动阻力和能量损失7-1 流体的两种流动形态——层流和湍流,流态的判别准则7-2 恒定均匀流基本方程,沿程损失的普遍表示式7-3 层流沿程损失的分析和计算,圆管层流的沿程损失系数7-4 湍流理论基础,湍流的脉动和时均法,湍流附面层分区的判别标准7-5 湍流沿程损失的分析和计算7-6 局部损失的分析和计算第八章边界层理论基础和绕流运动8-1 边界层的基本概念8-3 边界层的动量积分方程8-4 平板上的边界层8-5 边界层的分离现象和卡门涡街8-6 绕流运动参考书目:工程流体力学(水力学)(第2版)(上册),闻德荪,高等学校教材,第三版,2010年。
流体⼒学第⼋章绕流运动第⼋章绕流运动⼀、应⽤背景1、问题的⼴泛存在性:在⾃然界和⼯程实际中,存在着⼤量的流体绕物体的流动问题(绕流问题),如:飞机在空⽓中的飞⾏、河⽔流过桥墩、⼤型建筑物周围的空⽓流动、植物护岸(消浪,船⾏波),粉尘颗粒在空⽓中的飞扬和沉降,⽔处理中固体颗粒污染物在⽔中的运动。
(⼀种:流体运动;另外⼀种:物体运动),我们研究,将坐标系固结于物体上,将物体看成静⽌的,讨论流体相对于物体的运动。
2、问题的复杂性上⼀章的内容中可以看出,流体⼒学的问题可以归结为求解在⼀定边界条件和初始条件下偏微分⽅程组的求解。
但描述液体运动的⽅程式⾮常复杂的:⼀⽅⾯,是⽅程的⾮线性性质,造成⽅程求解的困难;另⼀⽅⾯,复杂的边界条件和初始条件都给求解流体⼒学造成了很多⿇烦。
迄今为⽌,只有很少数的问题得到了解决。
平⾯泊萧叶流动,圆管coutte流动等等。
⽽我们所要解决的绕流问题正是有着⾮常复杂的边界条件。
3、问题的简化及其合理性流体⼒学对此的简化则是,简化原⽅程,建⽴研究理想液体的势流理论。
实际液体满⾜势流运动的条件:粘性不占主导地位,或者粘性还没有开始起作⽤。
正例:远离边界层的流体绕流运动、地下⽔运动、波浪运动、物体落⼊静⽌⽔体中,⽔的运动规律研究。
反例:研究阻⼒规律、能量损失、内能转换等等。
圆柱绕流(经典之⼀)半⽆限长平板绕流(经典之⼆)分成两个区域:⼀个区域是远离边界的地⽅,此区域剪切作⽤不明显,⽽且流体惯性⼒的影响远远⼤于粘性⼒的影响(理想液体)(引导n-s⽅程);另⼀个是靠近边界的地⽅(附⾯层,粘性底层),此区域有很强烈的剪切作⽤,粘性⼒的影响超强,据现代流体⼒学的研究表明,此区域是产⽣湍流的重要区域,有强烈的剪切涡结构,但此区域只有⾮常薄的厚度。
此区域对绕流物体的阻⼒、能量耗损、扩散、传热传质都产⽣重要影响。
4、本章的主要研究内容(1)外部:理想液体,(简化⽅法,求解⽅式)、(2)内部:附⾯层理论,(简化⽅法,求解⽅式,求解内容,现象描述)(3)两者的衔接。
8第⼋章-边界层理论基础和绕流运动第⼋章边界层理论基础和绕流运动8—1 设有⼀静⽌光滑平板宽b =1m ,长L =1m ,顺流放置在均匀流u =1m/s 的⽔流中,如图所⽰,平板长边与⽔流⽅向⼀致,⽔温t =20℃。
试按层流边界层求边界层厚度的最⼤值δmax 和平板两侧所受的总摩擦阻⼒F f 。
解:20℃⽔的运动粘度ν=1.003?10-6 m 2/s 密度3998.2/kg m ρ=6119970091.00310ν-?===?L uLRe 因为 56310997009310?<=按层流边界层计算。
max 1/25.4470.0055m Re L L δ===3f 1/21.46 1.4610-===?L C Re 223998.2122 1.461011N 1.46N 22f ff u F C A ρ-?=== 8—2 设有极薄的静⽌正⽅形光滑平板,边长为a ,顺流按⽔平和铅垂⽅向分别置放于⼆维恒定均速u 的⽔流中,试问:按层流边界层计算,平板两种置放分别所受的总摩擦阻⼒是否相等,为什么?解:因为两种置放情况的物理模型和数学模型及其分析、推导所得计算公式是相同的,所以两种情况平板所受的总摩擦阻⼒相等。
8—3 设有⼀静⽌光滑平板,如图所⽰,边长1m,上宽0.88m,下宽0.38m,顺流铅垂放置在均匀流速u =0.6m/s 的⽔流中,⽔温t=15℃。
试求作⽤在平板两侧的总摩擦阻⼒F f 。
注:若为层流边界层,C f 按式(8—24)计算。
解:由表1—1查得,15℃时⽔的密度ρ=999.13/kg m ,运动粘度ν=1.139×10-6m 2/s 。
⾸先判别流态,计算平板上宽雷诺数560.60.884635655101.13910-?===Re ,按层流边界层计算。
设z 轴铅垂向上,平板宽度x 为0.38+0.5z ,阻⼒系数C f 按式(8-24)计算,即12f 60.6(0.380.5)1.328 1.13910--?+??==z C1521.328 5.2677810(0.380.5)z -轾=创?犏臌总摩擦阻⼒F f 按式(8—20)计算,f f12012(0.380.5)d 2F C u z z r =?ò11522 1.328 5.2677810(0.380.5)z -轾=创创+犏臌ò题8-1图21999.10.6(0.380.5)d 2z z 创创+ 110.658(0.380.5)d z z =ò。
边界层及绕流由于流体粘滞性的存在,紧靠平板的一层流体质点将附着于平板表面上,与平板表面无U,相对运动,流速为0,而在距平板法线方向一定距离处流速仍为未受扰动的原有流速因此从平板表面到未扰动的流体之间存在着一个流速分布不均匀的区域,这个区域就是水流受平板影响的范围叫边界层。
边界层厚度常用符号δ表示。
边界层的厚度是沿平板而变化的。
因为粘滞流体流经平板时有内摩擦阻力发生,克服阻力必耗损一部分能量,以致平板附近部分水流的流速变缓,流经平板距离越长,耗损能量越多,水流受平板影响范围也越大,所以边界层的厚度总是沿板端的距离x而增加的。
边界层内的流体形态可能是层流,也可能是紊流。
在板端附近边界层极薄,流速自0U,因此流速剃度极大,以致产生很大的内摩擦阻力,所以板端附近边界层内的迅速增至流体往往是层流。
沿板端距离越远,边界层厚度越厚。
流速剃度随边界层厚度增加而变小,内摩擦阻力也相应减小,边界层内的流体可自层流逐渐过渡到紊流。
但在紊流边界层中靠近固体表面仍有一层极薄的粘性存在,如图所示若雷诺数用下列形式表示:0Re x U xγ=则距板端距离越远,雷诺数也越大。
当雷诺数达到某一临界值时,流体即自层流转变为紊流。
据实验结果临界雷诺数约在5*510~610之间,如流体非常平静,最高的临界雷诺数也可超过610。
根据边界层的概念,可把粘滞流体分成两个区域:在边界层外,流速剃度为0,无内摩擦力发生,因而也可视为理想流体的流动,符合势流的运动规律;在边界层以内,流速自0增至0U ,流速剃度很大,内摩擦力十分显著。
因此,分析边界层内的运动规律时,必须以粘滞流体所服从的定律(纳为-斯托克斯方程式)为依据。
边界层的分离现象及绕流阻力流体压强在驻点N 处最大,在较高压强作用下,流体由此分道向圆柱体两侧流动。
由于圆柱面的阻滞作用便形成了边界层。
边界层内的特点是流体流动时有能量损失,从N 点起向下游达到A 或B 以前,由于圆柱表面的弯曲,使流体挤压,流速沿程增加,故沿边界层的外边界上0U x ∂∂=正值,p x∂∂=负值,即在外边界上压强是沿程下降的,由此可知在NA 或NB 一段边界层内的流体是处于加速减压状态的,也就是说,在该段边界层内用压强下降来补偿能量损失外,尚有一部分压能转变为动能。
第八章 边界层理论基础和绕流运动8—1 设有一静止光滑平板宽b =1m ,长L =1m ,顺流放置在均匀流u =1m/s 的水流中,如图所示,平板长边与水流方向一致,水温t =20℃。
试按层流边界层求边界层厚度的最大值δmax 和平板两侧所受的总摩擦阻力F f 。
解:20℃水的运动粘度ν=1.003⨯10-6 m 2/s 密度3998.2/kg m ρ=6119970091.00310ν-⨯===⨯L uLRe 因为 56310997009310⨯<=<⨯L Re按层流边界层计算。
max 1/25.4470.0055m Re L L δ===3f 1/21.46 1.4610-===⨯L C Re 223998.2122 1.461011N 1.46N 22f ff u F C A ρ-⨯==⨯⨯⨯⨯⨯= 8—2 设有极薄的静止正方形光滑平板,边长为a ,顺流按水平和铅垂方向分别置放于二维恒定均速u 的水流中,试问:按层流边界层计算,平板两种置放分别所受的总摩擦阻力是否相等,为什么?解:因为两种置放情况的物理模型和数学模型及其分析、推导所得计算公式是相同的,所以两种情况平板所受的总摩擦阻力相等。
8—3 设有一静止光滑平板,如图所示,边长1m,上宽0.88m,下宽0.38m,顺流铅垂放置在均匀流速u =0.6m/s 的水流中,水温t =15℃。
试求作用在平板两侧的总摩擦阻力F f 。
注:若为层流边界层,C f 按式(8—24)计算。
解:由表1—1查得,15℃时水的密度ρ=999.13/kg m ,运动粘度ν=1.139×10-6m 2/s 。
首先判别流态,计算平板上宽雷诺数560.60.884635655101.13910ν-⨯===<⨯⨯L uLRe ,按层流边界层计算。
设z 轴铅垂向上,平板宽度x 为0.38+0.5z ,阻力系数C f 按式(8-24)计算,即12f 60.6(0.380.5)1.328 1.13910--⨯+⎡⎤==⨯⎢⎥⨯⎣⎦z C1521.328 5.2677810(0.380.5)z -轾=创?犏臌总摩擦阻力F f 按式(8—20)计算,f f12012(0.380.5)d 2F C u z z r =?ò11522 1.328 5.2677810(0.380.5)z -轾=创创+犏臌ò题8-1图21999.10.6(0.380.5)d 2z z 创创+ 1120.658(0.380.5)d z z =?ò。
因0.380.5x z =+,所以d 0.5d x z = ,或d 2d =z x 。
代入上式得0.88130.8822f 0.380.3820.6582 1.3163=⨯⨯=⨯⎰F x dx x0.88(0.830.23)N 0.528N =?=8—4 油的动力粘度μ=50×10-3Pa ·s ,密度r =990kg/m 3,流速u =0.3m/s ,流过一水平放置的静止光滑平板。
试求距离平板始端150mm 处的边界层厚度δ以及边界层厚度为50mm 处距离平板始端的距离L 。
解:(1)30.39900.158915010x u x Re r m -创===´,为层流边界层。
5.4770.028m d ==?(2)0.05m d =时,假设仍为层流边界层0.05===0.495m L =30.39900.49529405010L Re -创==´,为层流边界层。
0.05m d ==8—5 试按光滑平板上的湍流边界层计算习题8—1中平板上边界层厚度的最大值maxd 和平板两侧所受的总摩擦阻力F f 。
解:max 150.3810.024m Re L Ld ===3f 150.074= 4.6710LC Re -==? ()2f f f =22u F C A r 两侧23998.212 4.671011N4.66N 2-´=创创?max d 、F f 值均大于习题8—1按层流边界层计算所得的值。
8—6 空气的温度t =0℃,流速u =30m/s ,在一个标准大气压下,流过一水平放置的(静止)光滑平板。
已知距平板始端4m 处的某点流速u x =27m/s ,试求该点距平板的垂直距离y 。
解:t =0℃时,空气的动力粘度μ=1.71×10-5Pa ·s ,密度31.293kg/m r =。
51.29330490736841.7110x ux Re r m -创===´,在5731010´范围内。
按湍流边界层计算150.3810.3810.062m xx Re d ==?17x y u u d骣÷ç=÷ç÷ç桫 77270.062m 0.03m 30x u y ud 骣骣÷÷çç=??÷÷çç÷÷çç桫桫8—7 有一宽b =2.5m ,长L =30m 的光滑平板潜没在静水中,以5m/s 的速度等速拖曳,平板长边与运动方向一致,水温为20℃,试求光滑平板的总摩擦阻力F f 。
解:t =20℃时,水的运动粘度n =1.003×10-6m 2/s ,密度998.2r =kg/m 3。
65301.00310o L U L Re n -´===´149551346>107,按湍流边界层计算。
()fm 2.580.4550.002lg L C Re == 220fm f 998.25220.002 2.530N 3743.25N 22U F C A r ´==创创=总8—8 空气的温度t =40℃,流速U 0=60m/s ,流过一长L =6m ,宽b =2m 的光滑平板, 平板长边与流速方向一致。
设平板边界层由层流转变为湍流的条件为60crcr 10x U x Re n==。
试求平板两侧所受的总摩擦阻力F f (注:按混合边界层计算)。
解:t =40℃时,空气的运动粘度521.6810m /s n -=?,密度31.128kg/m r =。
60560621428571101.6810L U L Re n -´===>´,按混合边界层计算。
Re x cr =106,由表8-1可查得A=3300()fm 11550.0740.07433000.002352142857121428571L LA C Re Re =-=-= 220f fm 1.12860220.0023526N 114.5N 22U F C bLr ´=?创创=8—9 空气的温度为293K ,流速u =30m/s,在一个标准大气压下,流过一水平放置的光滑平板。
层流边界层转变为湍流边界层的临界雷诺数cr x Re 5510=?,试求(1)边界层流态转变处离平板始端距离x cr 和该处离平板垂直距离y =1mm 处的流速u x ;(2)离平板始端1m 处的边界层厚度和每米宽平板所需的总拖曳力F f 。
(按混合边界层计算)解:(1)t =293K 时,空气的动力粘度51.8110Pa s m -=醋,密度31.205kg/m r =。
cr cr =x Re ux rm55cr cr 510 1.8110m 0.25m 1.20530x Re x u m r -创?===´cr 1/2cr 5.4775.4770.00194m x x Re d ==? 2222300.001()(0.001)m/s 22.96m/s 20.0019420.00194x u y u y d d ´=-=-=´(2) 551.205301199********.8110L uL Re r m -创===>?´为湍流边界层。
1/5=0.3810.3810.0209m L x Re d =?f m 1/50.074170017000.003221997238L L C Re Re =-==22f fm 1.20530220.0032211N 3.49N 22u F C bL r ´==创创=8—10 设有一宽b =2.5m ,长L =30m 的粗糙平板潜没在静水中,以5m/s 的速度等速拖曳,平板宽边b 与运动方向一致,水温为20℃,平板当量粗糙度∆=0.3mm 。
试求粗糙平板的总摩擦阻力F f 。
解:由表1—1查得,水温t =20℃时,水的密度ρ= 998.2 kg/m 3,运动粘度ν=1.003×10-6m 2/s 。
5065 2.5124626125101.00310ν-⨯===>⨯⨯b U bRe ,为湍流边界层。
允许粗糙度∆'650100100 1.00310m 2.00610m 0.02mm 0.3mm /5U ν--⨯⨯≤==⨯≈<∆=>14×0.02=0.28mm ,粗糙平板,且可认为属于湍流边界层粗糙区。
层流边界层长度 56cr cr 0510 1.00310m 0.1m 5-⨯⨯⨯===x Re v x U与平板宽边b =2.5m 相比,可略去不计。
按湍流边界层粗糙区计算摩阻系数C f ,即2.5 2.5f 2.5(1.62lg1.89)(1.62lg 1.89)0.0003--=+=+∆b C 31013.5-⨯= 2320f f 2 5.1310998.25 2.5302N 22U F C A r -创创创== 9601.44N =8—11 球形尘粒密度s ρ= 2.5×103kg/m 3,在20℃的大气中等速自由沉降。
若空气阻力可按斯托克斯阻力公式计算,试求尘粒最大直径d max 和自由沉降速度u f 。
解:由表1-2查得空气的运动粘度ν=1.5×10-5m 2/s ,密度ρ=1.205 kg/m 3 。
50515.8410m58.4μm d -===?max 58.4μm d =2f 1()18s u d g r r nr=- 2510.000058418 1.510 1.205-=?创? (2.51000 1.205)9.8m/s 创-? f u 0.257m/s = 8—12 球形水滴在20℃的大气中等速自由沉降,若空气阻力可按斯托克斯阻力公式计算,试求水滴最大直径d max 和自由沉降速度u f 。
解:(1)由表1-2查得空气的运动粘度ν=1.5×10-5m 2/s ,密度ρ=1.205 kg/m 3 ,水的密度F ρ= 998.2kg/m 3。