Anderson《空气动力学基础》读书笔记（十一）：黏流与边界层——普朗特的开创

本文是 John D. Anderson《Fundamentals of Aerodynamics》（第六版）读书笔记系列的第十一篇，覆盖教材第12章《黏流与边界层》的核心内容。这是理解真实流动物理的关键——解释了为什么势流理论预测的”零阻力”在现实中并不存在。

---

第十一篇：黏流与边界层——普朗特的开创

一、边界层——空气动力学史上最重要的概念

1.1 概念引入

1904 年，在海德堡举行的第三届国际数学大会上，一位 29 岁的德国学者提出了一个看似简单的想法，却永远改变了空气动力学的面貌。

Ludwig Prandtl 指出：

对于黏性很小的流体（如空气），黏性效应局限在物体表面附近一个很薄的区域内。在这个区域以外，流动可以看作无黏的势流。

这个”薄区域”就是边界层。

1.2 边界层的定义

边界层是紧贴在物体壁面上的流体层，其中：

1. 速度梯度很大——从壁面零速度（无滑移条件）变化到外部势流速度
2. 黏性剪切力不可忽略
3. 厚度 $\delta(x)$ 沿流向增大

边界层厚度 $\delta$ 的定义通常为：从壁面到速度达到 $0.99V_e$ 的距离，其中 $V_e$ 是边界层外缘速度。

1.3 无滑移条件

边界层理论的基石是无滑移条件（No-slip Condition）：

$$V_{\text{wall}} = 0$$

在壁面上，流体分子与壁面之间的黏附力使得紧贴壁面的空气”粘”在表面上。这个条件使得壁面附近形成巨大的速度梯度，从而产生了黏性剪切力——壁面摩擦阻力。

---

二、边界层的基本方程——Prandtl 边界层方程

2.1 从 Navier-Stokes 到边界层方程

完整的 Navier-Stokes 方程极其复杂。Prandtl 的关键洞察是：在边界层内，流向尺度远大于法向尺度，因此可以简化：

边界层假设（Prandtl, 1904）：

1. 边界层厚度 $\delta \ll L$（特征长度）
2. 法向速度 $v \ll u$（流向速度）
3. 法向压力梯度 $\partial p/\partial y \approx 0$（压力穿透边界层）

在这些假设下，不可压缩边界层方程为：

流向动量方程：

$$u\frac{\partial u}{\partial x} + v\frac{\partial u}{\partial y} = -\frac{1}{\rho}\frac{dp}{dx} + \nu\frac{\partial^2 u}{\partial y^2}$$

连续方程：

$$\frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} = 0$$

边界条件：

$$y=0: u=0, v=0$$ $$y\to\infty: u \to V_e(x)$$

这些方程虽然在形式上是偏微分方程，但可以通过相似性变换（Blasius 解）简化为常微分方程。

---

三、层流边界层——Blasius 解

3.1 Blasius 的贡献

1908 年，Prandtl 的学生 Heinrich Blasius 在他的博士论文中，成功求解了平板层流边界层方程。

Blasius 引入了一个无量纲变量：

$$\eta = y\sqrt{\frac{V_\infty}{\nu x}}$$

并假设流函数形式为：

$$\psi = \sqrt{\nu V_\infty x}\, f(\eta)$$

3.2 Blasius 方程的解

Blasius 将偏微分方程转化为一个三阶常微分方程：

$$2f''' + f f'' = 0$$

边界条件：

$$\eta=0: f=0, f'=0$$ $$\eta\to\infty: f' \to 1$$

这个方程没有解析解，但 Blasius 用级数展开和数值积分求出了精确数值解。

3.3 层流边界层的关键结果

边界层厚度：

$$\delta(x) = \frac{4.91 x}{\sqrt{Re_x}} \approx \frac{5.0 x}{\sqrt{Re_x}}$$

其中 $Re_x = V_\infty x / \nu$。

壁面摩擦应力：

$$\tau_w(x) = \frac{0.332 \rho V_\infty^2}{\sqrt{Re_x}}$$

局部摩擦系数：

$$C_f(x) = \frac{\tau_w}{\frac12 \rho V_\infty^2} = \frac{0.664}{\sqrt{Re_x}}$$

平板总摩擦阻力系数：

$$C_D = \frac{1}{L}\int_0^L C_f(x)\,dx = \frac{1.328}{\sqrt{Re_L}}$$

3.4 手算例子

一块平板长 $L = 0.5\ \text{m}$，空气 $V_\infty = 20\ \text{m/s}$，$\nu = 1.46\times10^{-5}\ \text{m}^2/\text{s}$。

求层流边界层在板后缘的厚度和总摩擦阻力。

解：

雷诺数：

$$Re_L = \frac{20 \times 0.5}{1.46\times 10^{-5}} = 6.85 \times 10^5$$

在 $x=L$ 处：

$$\delta(L) = \frac{5.0 \times 0.5}{\sqrt{6.85\times 10^5}} = \frac{2.5}{827.7} = 0.00302\ \text{m} \approx 3\ \text{mm}$$

总摩擦阻力系数：

$$C_D = \frac{1.328}{\sqrt{Re_L}} = \frac{1.328}{827.7} = 0.00160$$

平板总摩擦阻力（单位展长）：

$$D_f = C_D \cdot \frac12 \rho V_\infty^2 \cdot L = 0.00160 \times 0.5 \times 1.225 \times 400 \times 0.5 \approx 0.196\ \text{N/m}$$

这个例子说明：在 $Re \approx 10^6$ 时，层流摩擦阻力相当小。但在实际工程中，$Re > 5\times 10^5$ 时边界层通常会转捩为湍流，阻力会显著增大。

---

四、湍流边界层——更厚的层、更高的阻力

4.1 转捩

层流到湍流的转捩（Transition）发生在临界雷诺数 $Re_{\text{crit}}$ 附近：

• 平板：$Re_{\text{crit}} \approx 5\times 10^5$（取决于来流湍流度）
• 机翼：$Re_{\text{crit}} \approx 10^5 - 10^6$
• 光滑球体：$Re_{\text{crit}} \approx 3\times 10^5$

4.2 湍流边界层特性

湍流边界层比层流厚得多，壁面摩擦大得多，但有一个重要的工程优势——更难分离。

湍流边界层的速度分布（1/7 次幂分布近似）：

$$\frac{u}{V_\infty} = \left(\frac{y}{\delta}\right)^{1/7}$$

厚度（1/7 次幂近似）：

$$\delta(x) = \frac{0.382 x}{Re_x^{1/5}}$$

局部摩擦系数：

$$C_f(x) = \frac{0.0592}{Re_x^{1/5}}$$

总摩擦阻力系数：

$$C_D = \frac{0.074}{Re_L^{1/5}}$$

4.3 手算例子对比

对同一块平板（$L=0.5\ \text{m}$，$V_\infty = 20\ \text{m/s}$）：

参数	层流	湍流	比值（湍/层）
$\delta$ 后缘	3.0 mm	7.7 mm	2.6 倍
$C_D$	0.00160	0.00635	4.0 倍
$D_f$	0.196 N/m	0.779 N/m	4.0 倍

湍流边界层的摩擦阻力是层流的约 4 倍。这就是为什么层流机身和机翼的设计可以显著降低摩擦阻力——但维持层流本身非常困难。

---

五、边界层分离——失速的物理根源

5.1 分离的物理机制

当外部势流中的逆压梯度（$dp/dx > 0$）足够强时，壁面附近的低速流体无法克服逆压而被”推回”，导致流线与壁面脱离。

分离判据：

$$\left.\frac{\partial u}{\partial y}\right|_{y=0} = 0 \quad \rightarrow \quad \text{分离点}$$

5.2 层流 vs 湍流分离

层流边界层：更容易分离（$C_f$ 小，壁面附近动量低）
湍流边界层：更难分离（$C_f$ 大，湍流混合将高能流体带入壁面附近）

这就是机翼的实际升力特性：

• 层流翼型：失速较突然（薄翼失速）
• 湍流翼型：失速较平缓（后缘失速）

5.3 分离控制

工程中控制边界层分离的常见方法：

方法	原理	代价
涡流发生器	增加壁面附近湍流度	增加基准阻力
前缘缝翼	向后缘注入高能气流	增加重量和阻力
边界层吹气	主动向壁面附近输送高能气体	系统重量和能耗
边界层吸气	移除低动量流体	系统重量

---

六、完整概念地图

Navier-Stokes 方程
        ↓
Prandtl (1904): 边界层假设
δ << L, v << u, ∂p/∂y ≈ 0
        ↓
Prandtl 边界层方程
        ↓
    ┌───────────┴───────────┐
    ↓                       ↓
层流边界层              湍流边界层
(solution by Blasius)    (empirical/semi-empirical)
    │                       │
    ├─ δ ∝ x/√Re            ├─ δ ∝ x/Re^(1/5)
    ├─ Cf ∝ 1/√Re           ├─ Cf ∝ 1/Re^(1/5)
    ├─ 阻力小 (1.328/√ReL)   ├─ 阻力大 (0.074/ReL^(1/5))
    └─ 易分离               └─ 难分离 (工程优势!)
                 ↓
          边界层分离
           ────┬────
               ↓
         失速 (Stall)
               ↓
    最大升力系数 Cl,max 的物理极限

---

七、核心公式速查卡

公式	含义	流动状态
$\delta = 5.0x/\sqrt{Re_x}$	层流边界层厚度	层流
$\delta = 0.382x/Re_x^{1/5}$	湍流边界层厚度	湍流
$C_f = 0.664/\sqrt{Re_x}$	层流局部摩擦系数	层流
$C_f = 0.0592/Re_x^{1/5}$	湍流局部摩擦系数	湍流
$C_D = 1.328/\sqrt{Re_L}$	层流平板总阻力系数	全程层流
$C_D = 0.074/Re_L^{1/5}$	湍流平板总阻力系数	全程湍流
$(\partial u/\partial y)|_{y=0} = 0$	分离判据	分离点
$Re_{\text{crit}} \approx 5\times 10^5$	平板转捩临界雷诺数	转捩

---

参考文献

1. Anderson, J. D. (2010). Fundamentals of Aerodynamics (6th ed., Chapter 12). McGraw-Hill.
2. Prandtl, L. (1904). “Über Flüssigkeitsbewegung bei sehr kleiner Reibung”. Verhandlungen des III. Internationalen Mathematiker-Kongresses, Heidelberg.
3. Blasius, H. (1908). “Grenzschichten in Flüssigkeiten mit kleiner Reibung”. Zeitschrift für Mathematik und Physik, 56, 1-37.
4. Schlichting, H., & Gersten, K. (2017). Boundary-Layer Theory (10th ed.). Springer.
5. White, F. M. (2006). Viscous Fluid Flow (3rd ed.). McGraw-Hill.

---

下一节：湍流模型与现代数值方法

从工程湍流模型（混合长度、k-ε、k-ω SST）到现代 CFD（有限体积法、网格技术、高阶格式），再到高保真仿真（LES、DNS、DES）的技术前沿。本系列的收官之作。