第五章 边界层理论
第五章 边界层理论

1Transport Phenomena, Xu Jian, 2009第五章边界层理论边界层概念 边界层方程 边界层分离2Transport Phenomena, Xu Jian, 20095.1 边界层概念在上述层流动量传递的若干实例的分析中,(1)形状简单;(2)引入了假设:管道无限长、忽略进口段影响。
实际问题要复杂得多。
边界层理论,粘滞力对动量传递影响的一般理论,是粘性流体力学的基础,也与热量传递过程和质量传递过程有着密切的关系。
3Transport Phenomena, Xu Jian, 20095.1 边界层概念Prandtl(1904)提出边界层概念,把统一的流场,划分成两个区域,边界层和外流区;其流体流动(沿流动方向和沿与流动方向垂直的方向)有不同的特点。
边界层:流体速度分布明显受到固体壁面影响的区域。
边界层的形成:¾壁面处流体的“不滑脱”no-slip ¾流体的“内摩擦”作用 边界层厚度δ¾U =0Æ0.99 U 04Transport Phenomena, Xu Jian, 20095Transport Phenomena, Xu Jian, 20095.1 边界层概念流过一物体壁面的流体分成两部分¾边界层,粘性流体,不能忽略粘滞力¾外流区,理想流体,可以忽略粘滞力6Transport Phenomena, Xu Jian, 2009边界层理论的要点边界层厚度δ的变化¾前缘处,δ=0¾x ↑, δ↑;沿壁面的法向将有更多的流体被阻滞¾δ<<x边界层内,δ<<x (距离很小);0Æ0.99 U 0(速度变化大)¾速度梯度很大,剪切力很大¾流体速度减慢Æ惯性力<<层外,惯性力与粘性力数量级相当7Transport Phenomena, Xu Jian, 2009边界层流动的转变x<x c (临界距离)层流边界层 过渡区 湍流边界层转变判据:¾临界值:5×105;¾特征长度:距前缘的距离;¾特征速度:来流速度0Re xU ρμ=8Transport Phenomena, Xu Jian, 2009圆管进口段效应靠近管壁部分:边界层,速度减慢;厚度不断增大,进口段长度之后,汇交在管中心处;充分发展段的流动状态取决于交汇处边界层的流动状态;进口段的中心部分:无粘性流动区,速度均匀,区域不断缩小,在边界层汇交时消失;沿程速度不断增大Î压降增大(附加压降);9Transport Phenomena, Xu Jian, 20095.2 边界层方程普兰德边界层方程:量级比较 边界层积分动量方程:动量衡算沿平壁层流边界层的计算:动量积分方程的应用10Transport Phenomena, Xu Jian, 20095.2.1 普兰德边界层方程2222222211x x x x xy y y y y x y u u u u P u u x y x x y u u u u P u u x y y x y μρρμρρ⎛⎞∂∂∂∂∂+=−++⎜⎟∂∂∂∂∂⎝⎠⎛⎞∂∂∂∂∂+=−++⎜⎟⎜⎟∂∂∂∂∂⎝⎠讨论不可压缩流体在平板壁面上的稳态二维层流2222221x x x x Du u u u PDt x x y z υρ⎛⎞∂∂∂∂=−+++⎜⎟∂∂∂∂⎝⎠2222221y y y yDu u u uPDtyx y z υρ⎛⎞∂∂∂∂=−+++⎜⎟⎜⎟∂∂∂∂⎝⎠不可压缩流体的Navier-Stocks 方程不可压缩流体在边界层中作稳态二维流动,方程简化为:y0x u u x y∂∂+=∂∂连续性方程:11Transport Phenomena, Xu Jian, 20095.2.1 普兰德边界层方程普兰德首先发现可以通过比较数量级简化方程:¾Re 较大时,边界层的厚度δ<<x¾边界层内的惯性力和粘性力数量级相当 标准数量级:¾x 为距离的标准数量级,记为x=O(1)¾u 0为速度的标准数量级,记为u 0=O(1)¾边界层厚度δ的数量级记为δ= O(δ),远远小于O(1) 其他物理量的数量级:¾u x 与u 0是一个数量级,记为u x =O(1)¾y 与u 0是一个数量级,记为u x =O(1)12Transport Phenomena, Xu Jian, 2009其他物理量的数量级(1)(1)(1)x x u u O O x x O ∂Δ≈==∂Δ()222(1)(1)(1)(1)x x u u O O x O O x ∂Δ≈==∂Δyx u u x y ∂∂+=∂∂(1)x u O x∂=∂+(1)y u O y∂=∂()y u O δ=(1)1()()x x u u O O y y O δδ∂Δ≈==∂Δ()22222(1)1()()x x u u O O y O y δδ∂Δ≈==∂Δ22221x x x x xy u u u u P u u x y x x y μρρ⎛⎞∂∂∂∂∂+=−++⎜⎟∂∂∂∂∂⎝⎠数量级(1)(1)×1()()δδ×(1)21()δ13Transport Phenomena, Xu Jian, 2009其他物理量的数量级22221x x x x xy u u u u P u u x y x x y μρρ⎛⎞∂∂∂∂∂+=−++⎜⎟∂∂∂∂∂⎝⎠(1)(1)×1()()δδ×(1)21()δInertial Force=Viscous Force:2()O μδρ=1(1)PO xρ∂≤∂22221y yy yx y u u u u P u u x y y xy μρρ⎛⎞∂∂∂∂∂+=−++⎜⎟⎜⎟∂∂∂∂∂⎝⎠(1)()δ×()(1)δ×()δ1()δ2()δ()δ≤()δ(1)14Transport Phenomena, Xu Jian, 2009普兰德边界层方程2210x x xxy yx u u u dP u u x ydx y u u x yμρρ∂∂∂+=−+∂∂∂∂∂+=∂∂000x y x y u u y u u ====∞=时,时,普兰德边界层方程B.C.通过数量级比较得到的简化方程:应用条件:不可压缩流体在边界层中作稳态二维流动,而且Re 比较大15Transport Phenomena, Xu Jian, 20095.2.2 边界层积分动量方程卡门避开使用N-S 方程,直接对边界层进行衡算x 方向质量衡算:¾左侧进入:¾右侧流出:¾上部外流区进入yxz dxdy 1个单位距离δlyu 0, ρμlx u dy ρ∫()00ll x xu dy u dy dxxρρ∂+∂∫∫()lx u dy dxxρ∂∂∫()()2220000000u (-u )ll l l x x x xlx x u dy u dy dx u dy u dy dxx xdx u u dyx ρρρρρ∂∂+−−∂∂∂=∂∫∫∫∫∫x 向净动量变化率:不可压缩流体沿平板壁面的稳态二维流动16Transport Phenomena, Xu Jian, 2009边界层积分动量方程作用于控制体的x 向外力¾壁面剪切力:¾作用在左右侧面的压力差:1s dx τ−⋅⋅1Pdx l x∂−⋅⋅∂00(u )l x x s Pu u dy l x xρτ∂∂−=+∂∂∫0[,]x y l u u δ∈=00(u )x x sP u u dy x xδρδτ∂∂−=+∂∂∫只考虑x 方向的流动00(u )x x s d dPu u dy dx dxδρδτ−=+∫边界层内外压力近似相等00(u )x x sd u u dy dx δρτ−=∫卡门边界层积分动量方程17Transport Phenomena, Xu Jian, 2009边界层积分动量方程可以求出边界层厚度、流体阻力、曳力系数等;方程有u x ,τw ,δ三个变量,需要补充u x =f 1(y),τw =f 2(δ)的关系;需要预先假定一个速度分布方程才能求解,故只能算是一种近似的方法。
第五章 边界层理论

A2 0.332
x
v0
是平板流动边界层微 分方程解的最终结论。
5.0
5 .0
5.3. 边界层内积分方程
1.边界层积分方程的建立
M x ux dy
0 2 Wx uxux dy ux dy 0 0 l l
l
M x x
d l x u x dy u dy 0 0 x dx
速度的0.99处到固体壁面间的距离定义为边界层的厚度。
层流 底层
5.1. 边界层理论的基本概念
2 边界层的形成与特点:
① 形成:
流体流过平板,与平板紧临的流体受平板阻力而与平
板相对静止,边界层其余内各层流体自上而下依次受到 下层流体的粘性力作用而速度逐渐减小,这样就产生了 速度梯度较大的边界层。
5.1. 边界层理论的基本概念
d u0 u x u x dy 0 0 dx 3 ux 3 y 1 y u 2 2 0 u x y 0 0 a0 u x y u 0 2 3 u u a by cy dy y u0 b 3 u0 x y x y 0 2 y u x 2 b 2cy 3dy 2 0 u0 y y ux d y 0 3 2 y 0 2 y 2u x 2c 6dy y 0 2c 0 2 c0 y y 0
长度L,宽度B的平板总阻力
积分方程的解
4.64 x
v0 4.64 1 Re x
3
S
B
0
L
0
y 0
dxdz
3 0.646 v0 LB
空气动力学基础第五章边界层理论及其近似解读

u v 0 x y
ue u u u ue 2u u v ue 2 t x y t x y
在定常流动情况下,有
u v 0 x y
ue u u 2u u v ue 2 x y x y
5.2、平面不可压缩流体层流边界层方程
0 1 p y
这说明,在高Re数情况下,在边界层内压力沿法向是不变的。
5.2、平面不可压缩流体层流边界层方程
边界层内的压力分布与边界层外边界线上的压力分布相等。也就是, p与y无关,仅是x和t的函数。即
p pe ( x, t )
忽略质量力,Prandtl边界层方程变为
u v 0 x y
0
1 1
0
u dy e ue
这部分主流区增加的流体厚度是由边界层流体排挤入主流区造成的。因 此,称其为排移厚度。 (b)边界层动量损失厚度 在边界层内,在质量流量不变的条件下,理想流体通过的动量为
K i u e udy
0
由于粘性的存在,实际流体通过的动量为
v v v 1 p 2v 2v u v fy 2 2 t x y y x y
选取长度特征L,速度尺度ue,时间尺度t=L/ue,边界层近似假定:
5.2、平面不可压缩流体层流边界层方程
(1)根据边界层定义,纵向偏导数远远小于横向偏导数。
5.1、边界层近似及其特征
Prandtl的边界层概念,为人们如何计入粘性的作用开辟了划时代的途 径,因此称其为粘性流体力学之父。对整个流场提出的基本分区是: (1)整个流动区域可分成理想流体的流动区域(势流区)和粘性流体的 流动区域(粘流区)。 (2)在远离物体的理想流体流动区域,可忽略粘性的影响,按势流理论 处理。 (3)粘性流动区域仅限于物面近区的薄层内,称为边界层。既然是粘流 区,粘性力的作用不能忽略,与惯性力同量级,流体质点作有旋运 动。 2、边界层的特征 (1)边界层定义 严格而言,边界层区与主流区之间无明显界线,通常以速度达到主 流区速度的0.99U作为边界层的外缘。由边界层外缘到物面的垂直距离称 为边界层名义厚度。
第五章边界层理论

Y方向
按边界层概念: 边界层以外势流区的速度u∞不变,所以也不存在压力梯度 进一步简化:
H.布拉修斯对上述方程组进行了解析,引入流函数ψ(x,y),将 偏微分方程组化为可以解的常微分方程:
通常规定:u=0.99 u∞的位置为边界层的外边界线
5.2 平面层流边界层微分方程
以不可压稳态层流边界层为例: 1.微分方程建立与简化:
控制方程(二维,不可压,稳态,层流,不考虑质量力)
v x vy 0 x y
连续性方程
N-S方程
2v x 2v x v x v x 1 p vx vy 2 2 X方向 x y x x y
1.328 C f 1.328 0 L Re L
x 4.64 Re x
其中:Re L
不可压层流平板绕流摩擦阻力系数:
0 L
v
其总阻力:S
Cf 2
2 0 LB 其中L为平板长度,B为平
板宽度。
1. 平板紊(湍)流中速度分布与边界层厚度关系:
x y 17 ( ) 0
将流函数带入上面的方程组 并认为层流边界层内沿x轴各截面的速度分布图象相似 vx y F( ) v 又依
x
1 Re
则
x Re
y
y x Re
5.3 不同条件下边界层厚度与摩擦阻力系数
1. 平板层流中速度分布与边界层厚度关系:
x 3 y 1 y 3 ( ) ( ) 0 2 2
第五章 边界层理论
王连登 liandeng@ 13506970553
边界层理论

y 0
h
λ:流体导热系数
t
t y
y 0
∂t/∂y: 贴壁流体层的温度梯度 注意与导热问题第三类边界条件的区别
SJTU-OYH
对流传热问题的数学描述 假设: 流体为连续介质,流动为二维; 流体为不可压缩牛顿流体; 常物性、无内热源; 忽略粘性耗散热;
忽略辐射换热。 四个未知量:u, v, p, t。
x<xc, Re<Rec 层流 x>xc, Re>Rec 湍流
层流底层(粘性底层):紧靠壁面处,粘性力占主导地位,使 粘附于壁的一极薄层仍然会保持层流特征。层流底层内具有最 大的速度梯度。 SJTU-OYH
4-2边界层型对流传热问题的数学描述 波尔豪森热边界层的概念: 实验发现:流体对流换热时温 度梯度主要存在于近壁面的薄 层,主流区温度梯度几乎为零。
4-2边界层型对流传热问题的数学描述 波尔豪森热边界层的概念: 与边界层内速度分布一样,热边界 层内的温度分布也与流动形态密切 相关。 层流:温度呈抛物线分布 湍流:温度呈幂函数分布
y y w ,t w ,l
hw,t hw,l
u ' v' + =0 x' y '
u ' u ' p ' 1 2u ' 2u ' u ' +v' =- Eu 2 2 x' y ' x' Re x' y '
欧拉数 雷诺数 普朗特数 贝克莱数 努谢尔数
2 Eu p /( u )
v v v p 2v 2v ( u v ) Fy ( 2 2 ) x y y x y
17589第五章边界层理论及层流边界层中的传递现象

第五章边界层理论及层流边界层中的传递现象5.1 边界层理论的要点5.1.1 问题的提出前述,Re∝惯性力/粘性力当Re<1时,惯性力<<粘性力,可用“爬流”模型,略去惯性力项,N-S方程==>爬流方程(stokes近似),解决一些实际问题(沉降、润滑、渗流等),获得比较满意的结果。
但工程流动问题,绝大多数的Re很大。
这时,是否可以完全略去粘性力,使Navier-Stokes方程==>Euler方程(理想流体)。
但是,这样的结果与实际情况相差很大。
突出的一例即“达朗倍尔佯谬(paradox)——在流体中作等速运动的物体不受阻力”。
究竟应当怎样才能正确地处理大Re数的流动呢?这个矛盾一直到1904年,德国流体学家普兰德(Prandtl)提出了著名的边界层理论(大Re数的流动中,大部分区域的惯性力>>粘性力,但在紧靠固——流边界的极薄流层中,惯性力≈粘性力),才令人满意地解决了大Re数的流动的阻力问题。
后人把Prandtl 提出的流动边界层概念,推广到流动系统的传热边界层和传质边界层,从而确定传热、传质的速率以及了解有关的影响因素。
还有人研究了边界层中的化学反应,解决了一些实际问题。
因此,边界层理论被认为是近代流体力学的奠基石。
5.1.2 流动边界层(速度边界层)以平板流动为例,x方向一维稳态流动,在垂直壁面的y方向上,流动可划分为性质不同的两个区域:(1)y<δ(边界层):受壁面影响,法向速度变化急剧,du/dy很大,粘性力大(与惯性同阶),不能忽略。
(2)y>δ(层外主流层):壁面影响很弱,法向速度基本不变,du/dy≈0。
所以可忽略粘性力(即忽略法向动量传递)。
按理想流体处理,Euler方程适用。
这两个区域在边界层的外缘衔接起来,由于层内的流动趋近于外流是渐进的,不是突变的,因此,通常约定:在流动边界层的外缘处(即y=δ处),u x=0.99u∞,δ——流动边界层厚度,δ=δ(x)。
第5章-边界层理论基础PPT课件
虽然对Re很小的流动,惯性力可以忽略, 但对于Re很大的流动,粘性力却不能忽略, 否则会带来很大的误差,这是何故?
如水和空气,其粘度都很小,在处理其高
速流动时,如果忽略粘性力的影响,就会
导致与实际不符的错误结果。这个矛盾在
普兰德(Plandt)提出边界层学说之后,才获
得令人满意的解答。 -
-
20
卡门边界层方程即适用于层流,也适用 于湍流。
例:流体沿平板壁面流动时层流边界层 的计算,主要目标是边界层厚度和曳力 子数的计算
大量观察和测量得知ux与y的关系与抛 物线近似,因此可假设:
uxabycy2dy3 a,b,c,d 待定
边界条件:
-
21
y 0处ux 0 a 0
dux dy
-
5
随着边界层的厚度逐渐增加,边界层内
部也会发生变化,在边界层厚度较小处,
其内部流动为层流,该区域称为层流边
界层,当其厚度达到其临界厚度δc或临
界距离xc时,其内的流动逐渐经过一过
渡区转变为湍流,此后的边界层称为湍
流边界层,即使在这区域靠近壁面极薄
的一层流体内,仍然维持层流,称为层
流内层。
-
6
临界距离xc的长度与壁面前缘的形状、粗 糙度、流体性质和流速大小有关。壁面愈 粗糙xc愈短。
-
10
但实际中流速ux接近u0到一定程度时,便 可赋予其有应用价值的边界层厚度定义:
(1)
取ux达到u0的99%时的y值,即
ux u0
0 .9 9
处,y的值即为边界层厚度。
(2)可假设一个表示边界层内速度分布的
公式,如抛物线方程,计算当ux达到
u0时的y值,即为边界层厚度。
空气动力学:第5章 边界层理论及其近似
物体的特征长度。
EXIT
5.1、边界层近似及其特征
EXIT
5.1、边界层近似及其特征
EXIT
5.1、边界层近似及其特征
EXIT
5.1、边界层近似及其特征 (4)边界层各种厚度定义 (a)边界层位移厚度
假设某点P处的边界层厚度是 ,则在 的范围内 以速度 ue 流动的质量流量是:
EXIT
5.1、边界层近似及其特征
y 2
0 ,因而
z
t
0
。
在
z
的极大值点,2
x2
z
2
y 2
z
0 ,因而
z
t
0
。
这就说明了,在粘性流体中,不均匀的涡量场 是不断变化的,涡较强的部分要变弱,而涡较弱的 部分要变强。总的说来,趋向于涡量场强度“拉平 ”,就好像旋涡在扩散一样。
EXIT
5.1、边界层近似及其特征
(3)边界层厚度的量级估计 根据边界层内粘性力与惯性力同量级的条件,可估算边
(2)法向速度远远小于纵向速度。
v~
t
~
L / ue
L
ue
,
v ~
ue
1, Re
u
~
L t
ue ,
v u
(3)边界层内的压强量级与外流速度的平方成正比。
p ~ ue2
将这些量级关系式代入到N-S方程中,得到
EXIT
5.2、平面不可压缩流体层流边界层方程
N-S方程组与各项量级比较:
u v 0 x y
EXIT
5.1、边界层近似及其特征 (b)边界层各种厚度的定义式,既适用于层流,
也适用于湍流。 (c)边界层各种厚度的大小与边界层内流速分布
边界层理论Ch05-Boundary-Layer-Theory
4
Boundary layer equations
Applying our analysis to the momentum equations gives (x direction) :
μ ⎛ ∂ 2u x ∂ 2u x ⎞ 1 ∂p ∂u ∂u + = ux x + u y x ⎜ ⎟− ρ ⎝ ∂x 2 ∂y 2 ⎠ ρ ∂x ∂x ∂y
y ≈ dy ≈ u y ≈ du y ≈ O(δ h )
BUCT
BSc 2011
19
BUCT
BSc 2011
20
Boundary layer equations
2. The velocity in the x direction (along the plate) is much greater than the velocity in the y direction (perpendicular to the plate)
N-S Equations at x direction:
μ⎜
⎛ ∂ u x ∂ u x ∂ u x ⎞ ∂pd + 2 + 2 ⎟− 2 ∂y ∂z ⎠ ∂x ⎝ ∂x
2 2 2
Boundary layer equations x direction
Stationary flow:
∂u x =0 ∂t
u x ≈ du x ≈ x ≈ dx ≈ p ≈ O (1)
BUCT
BSc 2011
21
BUCT
BSc 2011
22
Boundary layer equations
第五章边界层理论解读
第三节
以二维绕平面 流动为例来导出边 界层积分方程,如 固5-2所示。 首先对控制体 (单元体)做动量平 衡计算(在计算过程 中取垂直于纸面 z 方 向为单位长度):
边界层内积分方程
1)流体从AB面单位时间流入的动量记为 Mx 。由 图5-2知,从 AB 面单位时间流入的质量为
(5-10) 2)流体从 CD 面单位时间流出的动量记为 Mx+∆x: 从 CD 面单位时间流出的质量为
(3)湍流区:随着进流尺寸的进一步增加,使得Rex > 3×106,这时边界层内流动形态已进入湍流状态,边界 层的厚度随进流长度的增加而迅速增加。
应当注意,无论是对过渡区还是湍流区,边界层 最靠近壁面的一层始终做层流流动,这一层称为层流 底层,这主要是因为在最靠近壁面处壁面的作用使该 层流体所受的粘性力永远大于惯性力所致。这里要特 别说明的是,边界层与层流底层是两个不同的概念。 层流底层是根据有无脉动现象来划分,而边界层则是 根据有无速度梯度来划分的。因此,边界层内的流动 既可以为层流,也可以为湍流。
(5-14)
由动量守恒可得结
本章重点叙述了边界层的概念、特点,建立了边界 层的微分方程、积分方程,并介绍了求解方法。对平板 绕流摩擦阻力的计算也进行了简要介绍。实际上,边界 层理论是在数值模拟技术没有发展起来之前,人们为了 运用流体流动的控制方程去解决工程实际问题的一部分 重要工作。尽管现在数值模拟技术已经能够处理某些真 三维实际流体的运动规律,但是通过边界层理论的学习, 仍然可以领略前人在处理实际流体流动问题上的输妙简 化与抽象思考,这是科学方法最突出的特征,这是精确 的数值模拟所不能替代的。
(1)层流区:流体统流进入平板后,当进流长度不是 很长,x<xc(xc为对应Rex=2×105的进流深度),这时 Rex < 2×105,边界层内部为层流流动,这一个区域称 为层流区。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
1Transport Phenomena, Xu Jian, 2009第五章边界层理论边界层概念 边界层方程 边界层分离2Transport Phenomena, Xu Jian, 20095.1 边界层概念在上述层流动量传递的若干实例的分析中,(1)形状简单;(2)引入了假设:管道无限长、忽略进口段影响。
实际问题要复杂得多。
边界层理论,粘滞力对动量传递影响的一般理论,是粘性流体力学的基础,也与热量传递过程和质量传递过程有着密切的关系。
3Transport Phenomena, Xu Jian, 20095.1 边界层概念Prandtl(1904)提出边界层概念,把统一的流场,划分成两个区域,边界层和外流区;其流体流动(沿流动方向和沿与流动方向垂直的方向)有不同的特点。
边界层:流体速度分布明显受到固体壁面影响的区域。
边界层的形成:¾壁面处流体的“不滑脱”no-slip ¾流体的“内摩擦”作用 边界层厚度δ¾U =0Æ0.99 U 04Transport Phenomena, Xu Jian, 20095Transport Phenomena, Xu Jian, 20095.1 边界层概念流过一物体壁面的流体分成两部分¾边界层,粘性流体,不能忽略粘滞力¾外流区,理想流体,可以忽略粘滞力6Transport Phenomena, Xu Jian, 2009边界层理论的要点边界层厚度δ的变化¾前缘处,δ=0¾x ↑, δ↑;沿壁面的法向将有更多的流体被阻滞¾δ<<x边界层内,δ<<x (距离很小);0Æ0.99 U 0(速度变化大)¾速度梯度很大,剪切力很大¾流体速度减慢Æ惯性力<<层外,惯性力与粘性力数量级相当7Transport Phenomena, Xu Jian, 2009边界层流动的转变x<x c (临界距离)层流边界层 过渡区 湍流边界层转变判据:¾临界值:5×105;¾特征长度:距前缘的距离;¾特征速度:来流速度0Re xU ρμ=8Transport Phenomena, Xu Jian, 2009圆管进口段效应靠近管壁部分:边界层,速度减慢;厚度不断增大,进口段长度之后,汇交在管中心处;充分发展段的流动状态取决于交汇处边界层的流动状态;进口段的中心部分:无粘性流动区,速度均匀,区域不断缩小,在边界层汇交时消失;沿程速度不断增大Î压降增大(附加压降);9Transport Phenomena, Xu Jian, 20095.2 边界层方程普兰德边界层方程:量级比较 边界层积分动量方程:动量衡算沿平壁层流边界层的计算:动量积分方程的应用10Transport Phenomena, Xu Jian, 20095.2.1 普兰德边界层方程2222222211x x x x xy y y y y x y u u u u P u u x y x x y u u u u P u u x y y x y μρρμρρ⎛⎞∂∂∂∂∂+=−++⎜⎟∂∂∂∂∂⎝⎠⎛⎞∂∂∂∂∂+=−++⎜⎟⎜⎟∂∂∂∂∂⎝⎠讨论不可压缩流体在平板壁面上的稳态二维层流2222221x x x x Du u u u PDt x x y z υρ⎛⎞∂∂∂∂=−+++⎜⎟∂∂∂∂⎝⎠2222221y y y yDu u u uPDtyx y z υρ⎛⎞∂∂∂∂=−+++⎜⎟⎜⎟∂∂∂∂⎝⎠不可压缩流体的Navier-Stocks 方程不可压缩流体在边界层中作稳态二维流动,方程简化为:y0x u u x y∂∂+=∂∂连续性方程:11Transport Phenomena, Xu Jian, 20095.2.1 普兰德边界层方程普兰德首先发现可以通过比较数量级简化方程:¾Re 较大时,边界层的厚度δ<<x¾边界层内的惯性力和粘性力数量级相当 标准数量级:¾x 为距离的标准数量级,记为x=O(1)¾u 0为速度的标准数量级,记为u 0=O(1)¾边界层厚度δ的数量级记为δ= O(δ),远远小于O(1) 其他物理量的数量级:¾u x 与u 0是一个数量级,记为u x =O(1)¾y 与u 0是一个数量级,记为u x =O(1)12Transport Phenomena, Xu Jian, 2009其他物理量的数量级(1)(1)(1)x x u u O O x x O ∂Δ≈==∂Δ()222(1)(1)(1)(1)x x u u O O x O O x ∂Δ≈==∂Δyx u u x y ∂∂+=∂∂(1)x u O x∂=∂+(1)y u O y∂=∂()y u O δ=(1)1()()x x u u O O y y O δδ∂Δ≈==∂Δ()22222(1)1()()x x u u O O y O y δδ∂Δ≈==∂Δ22221x x x x xy u u u u P u u x y x x y μρρ⎛⎞∂∂∂∂∂+=−++⎜⎟∂∂∂∂∂⎝⎠数量级(1)(1)×1()()δδ×(1)21()δ13Transport Phenomena, Xu Jian, 2009其他物理量的数量级22221x x x x xy u u u u P u u x y x x y μρρ⎛⎞∂∂∂∂∂+=−++⎜⎟∂∂∂∂∂⎝⎠(1)(1)×1()()δδ×(1)21()δInertial Force=Viscous Force:2()O μδρ=1(1)PO xρ∂≤∂22221y yy yx y u u u u P u u x y y xy μρρ⎛⎞∂∂∂∂∂+=−++⎜⎟⎜⎟∂∂∂∂∂⎝⎠(1)()δ×()(1)δ×()δ1()δ2()δ()δ≤()δ(1)14Transport Phenomena, Xu Jian, 2009普兰德边界层方程2210x x xxy yx u u u dP u u x ydx y u u x yμρρ∂∂∂+=−+∂∂∂∂∂+=∂∂000x y x y u u y u u ====∞=时,时,普兰德边界层方程B.C.通过数量级比较得到的简化方程:应用条件:不可压缩流体在边界层中作稳态二维流动,而且Re 比较大15Transport Phenomena, Xu Jian, 20095.2.2 边界层积分动量方程卡门避开使用N-S 方程,直接对边界层进行衡算x 方向质量衡算:¾左侧进入:¾右侧流出:¾上部外流区进入yxz dxdy 1个单位距离δlyu 0, ρμlx u dy ρ∫()00ll x xu dy u dy dxxρρ∂+∂∫∫()lx u dy dxxρ∂∂∫()()2220000000u (-u )ll l l x x x xlx x u dy u dy dx u dy u dy dxx xdx u u dyx ρρρρρ∂∂+−−∂∂∂=∂∫∫∫∫∫x 向净动量变化率:不可压缩流体沿平板壁面的稳态二维流动16Transport Phenomena, Xu Jian, 2009边界层积分动量方程作用于控制体的x 向外力¾壁面剪切力:¾作用在左右侧面的压力差:1s dx τ−⋅⋅1Pdx l x∂−⋅⋅∂00(u )l x x s Pu u dy l x xρτ∂∂−=+∂∂∫0[,]x y l u u δ∈=00(u )x x sP u u dy x xδρδτ∂∂−=+∂∂∫只考虑x 方向的流动00(u )x x s d dPu u dy dx dxδρδτ−=+∫边界层内外压力近似相等00(u )x x sd u u dy dx δρτ−=∫卡门边界层积分动量方程17Transport Phenomena, Xu Jian, 2009边界层积分动量方程可以求出边界层厚度、流体阻力、曳力系数等;方程有u x ,τw ,δ三个变量,需要补充u x =f 1(y),τw =f 2(δ)的关系;需要预先假定一个速度分布方程才能求解,故只能算是一种近似的方法。
00(u )x x sd u u dy dx δρτ−=∫18Transport Phenomena, Xu Jian, 20095.2.3 沿平壁层流边界层的计算待定系数法:0nix i i u a y ==∑需要n +1个边界条件以确定n +1个待定系数23x u a by cy dy =+++边界条件:i. y=0, u x =0 (不滑脱)ii. y=0,220d 0x y u dy =⎛⎞=⎜⎟⎝⎠iii. y =δ, u x =u 0iv. y =δ, d 0xudy =3031u 22x u y y δδ⎛⎞=−⎜⎟⎝⎠19Transport Phenomena, Xu Jian, 2009沿平壁层流边界层的计算()300022000000011220031(u )223939280280314012134.641,0,0,10 4.64 4.64Re x x x x s x x sx s y xu d y y u u dy dx u d u u u dy u u dx du u d dxdy u x c x c u x u x δδρτδδδδρτμτμμδδδρμδμδδρρ=−⎛⎞⎛⎞−==−⎜⎟⎜⎟⎝⎠⎝⎠⇒−=⇒=⎛⎞==⇒=⎜⎟⎝⎠⎛⎞⇒=+===⇒==⎜⎟⎝⎠∫∫++处的剪切应力:1220300200.323Re 0.646,,1.2921Re 2sx xL sx D Lu b dx b Lu b L C u bL τρτμρρ−=====∫d d 总曳力F 宽度长度为F 20Transport Phenomena, Xu Jian, 2009【例】沿平壁层流边界层的计算温度为20℃的空气在常压下以5m/s 的速度流过一块宽1 m 的平板壁面。
试计算距平板前缘0.5m 处的边界层厚度及进入边界层内的质量流率,并计算这一段平板壁面的曳力系数与承受的摩擦曳力。
假设临界雷诺数Re xc =5×105。
()53550.5120.530.50000120.51.8110. 1.205/Re 1.664105100.52 4.64Re 0.00569315(3)0.0214/228(4) 1.292.Re 0.x D C Pa s kg m x my y u bdy u bdy u b kg sC δδμρδωρρρδδδ−−−=×==×<×==⎡⎤⎛⎞==−==⎢⎥⎜⎟⎝⎠⎢⎥⎣⎦==∫∫o 解:(1)判断边界层流型:20空气,处的边界层为层流边界层20003170.02382Du C bL Nρ==d F 21Transport Phenomena, Xu Jian, 20095.3 边界层分离边界层分离的必要条件¾流体具有粘性¾流动过程中存在逆压梯度。