流体力学边界层理论

于是
τ 0 = 0.332
μρU 2 x
上式可看出平板层流边界层局部摩擦切应力与ｘ坐标的平方根成反比的规
律随着ｘ的增加而减小。
现计算整个平板上总摩擦阻力。设板长为Ｌ，板宽为ｂ，则平板单面总摩擦
阻力是：
∫ ∫ Rf =
Lτ
0
0bdx
=b
L
0.332
0
μρU 3 dx = 0.664 x
μρ LU 3
总摩擦阻力系数 C f 由下式确定：
2
则：
vx
(
x,
y)
=
U
⋅
1 2
ϕ ′(η )
设 U=25 km/h，ν=0.15cm2/s, x=3m,ｙ=5mm，
求：Vx＝？
解：U=25×1000/3600=6.95m/s, ν=0.0015m2/s,
x=3m, y=0.005m,
代入η中得：
η = 1 × 5×10−3 × 2
6.95 0.15 ×10−4
（11-14）式应采取如下形式：
ϕ(x, y) = xϕ( y ) x
（11-16）
返回为有量纲解时，不出现Ｌ，即 ：
ϕ = ν U x ϕ (η )
η=1y U 2 νx
（11-18）
通过以上分析，来求解下列形式的ψ。
⎡y⎤
ϕ=
νUL
x
⎢ ⎢
L⎢
⎢ ⎣
νL ⎥
U ⎥=
x⎥
L
⎥ ⎦
⎡ νUxϕ ⎢ y
Ｕ（起参数作用），ν和Ｕ不同时，同一空间点上ψ的值不同。
现设法将方程和边界条件中各个物理量无量纲化，不再出现ν和Ｕ。
选特征量：
Ｌ：ｘ的比例尺
ν L ：y 的比例尺， U
Ψ: ψ的比例尺，Ψ为常数
若用ψ , x, y 表示ψ，ｘ及ｙ的无量纲值，则有
x= x L
y= y νL U
ϕ=ϕ ψ
(11-8)
x = Lx
Cf
=
Rf
1 2
ρU
2bL
= 1.328 Re
（11-21）
式中为按平板板长计算的雷诺数。算出摩擦阻力系数后，可确定平板层流 边界层情况下的摩擦阻力为：
Rf
= Cf
⋅ 1 ρU 2bL
2
（11-22）
虽然边界层基本微分方程比 N-S 方程要简单得多，但求解问题仍有很大困难 尚且如此之大，因此,发展求解边界层问题的近似方法便具有很大的理论与实际 意义。
讨论：
Prandtl 边界层方程中第二个方程： ∂p = 0 。说明了什么？ ∂y
说明了： P1 = P2 = P3 = 0
p0 p1
p2
Prandtl 边界层方程的求解:
Blasius 解----顺流放置无限长平板上的层流边界层流动。 均匀来流平行于平板，ｘ轴平行于板面，原点在平板前缘，平板极薄且无曲 度，边界层外缘处速度为来流速度Ｕ。沿边界层外缘上各点上压力相同，即 dp = 0 。 dx
1
δ′
∂v′x + ∂v′y = 0 ∂x′ ∂y′
1
1
因为 δ ~ 1 ，所以 Re ~ δ ′2 L Re
因为 0 ≤ x ≤ L ，所以 x′ = x ~ 1 L
因为 y′ = y ~ 1， 0 ≤ y ≤ δ 所以 y′ ~ δ = δ ′
L
L
因为 0 ≤ Vx
≤U
，所以Vx′ =
Vx U
~1
Rekp
=
(Ux ν
)
kp
= Uxkp ν
= 5×105
层流边界层转为湍流边界层转捩点的位置坐标
§8-2
xkp
=
5 ×105
ν U
边界层基本微分方程
（11-1）
粘性不可压缩流体,不计质量力，定常流过小曲率物体，物体表 面可近似当作平面。
取物面法线为ｙ轴。在大 Re 数情况下的边界层流动有下面两个主要性质：
所以
∂v′x = ∂v′x , ∂y′ ∂x′
∂v′x ~ 1, ∂x′
所以
∂v′y ~ 1, ∂y′
v′y ~ δ ′
∂ 2 v′x ∂x′2
~ 1，
∂v′y ~ δ ′,
∂x′
∂v′y ~ 1 ,
∂y′ δ ′ ∂2v′y ~ δ ′
∂x′2
∂v′x ~ 1
∂y′ δ ′′
(a )
(b )
(c)
平板上u=0u=0
边界层内粘 性力不可忽
一薄层内速度这
∂vx
梯度 ∂y 很大
边界层外边界
U99％
边界层名义厚度 :外边界上流速达到 U99％的点到物面的法向距离 边界层厚度：
根据速度分布的特点，可将流场分为两个区域：
一、边界层：
∂vx 1.这一薄层内速度梯度 ∂y 很大。
2.边界层内的流动是有旋流动 ωz
=
1 (∂vy 2 ∂x
−
∂vx ) ∂y
=
− ∂vx ∂y
。
二、边界层外部区域
边界层外部粘性影响很小，μ可以忽略不计，可认为边界层外部的流动是 理想流体无旋势流。
重要推论：
（1）边界层内各截面上压力等于同一截面上边界层外边界上的压力
即： P1 = P2 = P3
P P1
P2
x
（2）势流的近似计算中，可略去边界层的厚度，解出沿物体表面的流速和压力 分布，并认为就是边界层边界上的速度和压力分布，据此来计算边界层。
（3）根据边界层厚度极薄的基本假设，可将 N-S 方程化简，获得边界层的基本 微分方程。
边界层内的流动状态：
层流边界层，湍流边界层均存在粘性底层（层流底层） ，其厚度与 Re 有 关。
层流边界层转变为湍流边界层的判别准则：
雷诺数：
Re
=
Ux ν
（ｘ为离平板前缘点的距离）
对于平板，层流转变为湍流的临界雷诺数为：
∂v′x ∂x′
+ v′y
∂v′x ∂y′
=
−
∂p′ ∂x′
+
1 Re
(
∂ 2 v′x ∂x′2
+
∂ 2 v′x ∂y′2
)
1⋅1
δ2⋅1
δ
(δ 2 )
1
1
δ2
v
′
x
∂
v
′
y
∂x′
+
v
′
y
∂
v
′
y
∂y′
=
−
∂p′ ∂y′
+
1 Re
(
∂
2
v
′
y
∂x′2
+
∂
2
v
′
y
∂y′2
)
1⋅δ δ ⋅ 1
(δ 2 ) δ ′
上述边界层方程简化为：
vx
∂vx ∂x
+
vy
∂vx ∂y
=ν
∂ 2vx ∂x2
∂vx + ∂vy = 0
∂x
∂y
(11-5)
边界条件： ｙ＝０， Vx = 0,Vy = 0 ； ｙ→∞，Vx = U 。
严格上，速度从零增至Ｕ须经过无限远距离，近似认为ｙ＝δ，Vx = U 。
引入流函数ψ，与速度的关系为：
ux
=
∂ψ ∂y
ux
=
−
∂ψ ∂x
(11-6)
将其代入简化后的边界层方程第一式有：
∂ ψ ∂ 2ψ − ∂ ψ ∂ 2ψ = ν ∂ 3ψ
∂y ∂x∂y ∂x ∂y 2
∂y3
(11-7)
边界条件： y = 0
∂ψ = 0 ∂y
∂ψ = 0 x>0 ∂x
y→∞
∂ψ = U ∂y
若求出了流函数ψ，便可求出速度，ψ应是ｘ，ｙ的函数，且ψ中包含ν和
边界层：在固体壁面附近，显著地受到粘性影响的这一薄层。从边界层厚
度很小这个前提出发，Prandtl 率先建立了边界层内粘性流体运动的简化方程，
开创了近代流体力学的一个分支—边界层理论。
均匀来流绕一薄平板流动，微型批托管测得沿平板垂直方向的速度分布如 下图：
与来流速度相同的量级,U99％ 均 匀 来 流 速 度
Karman 动量积分方程方程，就是一种近似求解边界层问题的方法。
§8-3 边界层动量积分方程
应用动量定理来研究边界层内单位时间内沿ｘ方向的动量变化 和外力之间的关系。
设流动定
控制体边界ABCD
单位时间内经过ＡＢ面流入的质量和带入的动量分别为：
∫ ∫ mAB =
δ 0
ρuxdy
K AB =
δ 0
1) 边界层厚度较物体特征长度小得多，即：
δ′= δ
1
L
2)边界层内粘性力和惯性力具有相同的数量级以此作为基本假定，将
N-S 方程（二维）化简：
vx
∂vx ∂x
+
vy
∂vx ∂y
=−
1 ρ
∂p ∂x
+ν
(
∂ 2v x ∂x 2
+
∂ 2v x ∂y 2
)
vx
∂vy ∂x
+
vy
∂vy ∂y
=−
1 ρ
∂p ∂y
表 11-1 给出问题的数值解，其中 1 ϕ′(η) = vx 就是边界层内无量纲的速
2
U
度分布。
例 11.1 本例说明上表 11-1 的用法。
(1) 欲求边界层内点（x,y）的速度 Vx（x,y）可将ｘ及ｙ的值代入η = 1 y U ， 2 νx