Anderson《空气动力学基础》读书笔记（三）：势流理论（下）——环量与升力

本文是 John D. Anderson《Fundamentals of Aerodynamics》（第六版）读书笔记系列的第三篇，覆盖教材第3章后半部分与第4章前半部分。阅读本文前需熟悉速度势、流函数、点涡和偶极子的概念。

---

第三篇：环量与升力——空气动力学最核心的数学定理

势流理论的价值不是写一些优美的数学公式——它的终极成就是回答了空气动力学最基本的问题：

升力到底从哪里来？它的大小如何精确计算？

本篇将一步步推导出空气动力学中最重要的公式——库塔-儒科夫斯基定理（Kutta–Joukowski Theorem），以及它背后的物理和数学。

---

一、达朗贝尔佯谬——理论预测的”零升力零阻力”

1.1 问题

如果仅用均匀流 + 偶极子叠加出圆柱绕流，会得到什么结果？

均匀流 + 偶极子（偶极矩 $\mu = 2\pi V_\infty R^2$，其中 $R$ 为圆柱半径）：

速度势：

$$\phi = V_\infty\left(r + \frac{R^2}{r}\right)\cos\theta$$

流函数：

$$\psi = V_\infty\left(r - \frac{R^2}{r}\right)\sin\theta$$

表面速度（在 $r=R$ 处）：

$$V_\theta|_{r=R} = \left.\frac{1}{r}\frac{\partial \phi}{\partial \theta}\right|_{r=R} = -2V_\infty\sin\theta$$

由伯努利方程得到表面压力：

$$p(R,\theta) = p_\infty + \frac{1}{2}\rho V_\infty^2 - \frac{1}{2}\rho(2V_\infty\sin\theta)^2$$

积分压力得到总力——升力为 0，阻力也为 0。

这就是达朗贝尔佯谬（d’Alembert’s Paradox）：理想势流预测一个物体在流体中受到的阻力为零。

1.2 历史

时间	人物	进展
1745	达朗贝尔	在《流体阻力新论》中首次证明了”阻力为零”的矛盾结果
1752	达朗贝尔	进一步指出”我无法解释为什么流体会产生阻力”
1868	亥姆霍兹	提出间断流理论，尝试解决但未完全成功
1904	普朗特	边界层理论说明了阻力的真实来源（黏性）
1902	库塔	提出库塔条件——解决了”升力多大”的问题

达朗贝尔本人对这个问题非常苦恼。他在著作中写道：“这一悖论让流体动力学理论在面对工程应用时几乎完全失效。” 这个困境一直持续了 150 多年，直到普朗特的边界层理论和库塔的条件才真正解决。

1.3 解决问题的两条路径

1. 阻力问题（边界层理论）：黏性 → 壁面摩擦 + 压力阻力 → 总阻力 > 0。这在第12讲详细展开。
2. 升力问题（环量理论）：在圆柱绕流上叠加一个点涡，产生非对称的压力分布 → 升力 > 0。这就是本讲内容。

---

二、带环量的圆柱绕流——升力的诞生

2.1 概念引入

点涡是一个中心有环量 $\Gamma$ 的势流解。现在，如果我们把三个解叠加在一起：

均匀流 + 偶极子 + 点涡 = 带环量的圆柱绕流

2.2 速度势与流函数

速度势：

$$\phi = V_\infty\left(r + \frac{R^2}{r}\right)\cos\theta + \frac{\Gamma}{2\pi}\theta$$

流函数：

$$\psi = V_\infty\left(r - \frac{R^2}{r}\right)\sin\theta - \frac{\Gamma}{2\pi}\ln\frac{r}{R}$$

验证：圆柱表面（$r=R$）仍然是流线（$\psi = \text{const}$）：

$$\psi|_{r=R} = V_\infty(R-R)\sin\theta - \frac{\Gamma}{2\pi}\ln 1 = 0$$ $\psi=0$ 是流线，圆柱表面就是这条流线。壁面无穿透条件自动满足。

2.3 速度分布

径向和切向速度：

$$V_r(r,\theta) = \frac{\partial \phi}{\partial r} = V_\infty\left(1 - \frac{R^2}{r^2}\right)\cos\theta$$

在 $r=R$ 表面：

$$V_r|_{r=R} = 0 \quad \text{（符合壁面无穿透）}$$ $$\boxed{V_\theta|_{r=R} = -2V_\infty\sin\theta - \frac{\Gamma}{2\pi R}}$$

这个公式是本节最重要的结果。注意：环量 $\Gamma$ 的引入使得表面速度不再对称——上表面（$\theta \approx \pi/2$ 附近）速度加快，下表面（$\theta \approx -\pi/2$ 附近）速度减慢。

2.4 压力分布与升力

由伯努利方程，表面压力系数：

$$C_p = \frac{p - p_\infty}{\frac{1}{2}\rho V_\infty^2} = 1 - \left(\frac{V_\theta}{V_\infty}\right)^2$$

代入速度：

$$C_p(\theta) = 1 - \left(-2\sin\theta - \frac{\Gamma}{2\pi R V_\infty}\right)^2$$

垂直于来流方向的压力积分——就是升力：

$$L' = -\int_0^{2\pi} p(R,\theta) \sin\theta \, (R\,d\theta)$$

经过代数运算：

$$\boxed{L' = \rho V_\infty \Gamma}$$

这就是库塔-儒科夫斯基定理（Kutta–Joukowski Theorem）。

2.5 手算例子

均匀流速度 $V_\infty = 10\ \text{m/s}$，圆柱半径 $R = 0.1\ \text{m}$，空气密度 $\rho = 1.225\ \text{kg/m}^3$。

添加环量 $\Gamma = 5\ \text{m}^2/\text{s}$。

求：单位展长的升力。

解：

$$L' = \rho V_\infty \Gamma = 1.225 \times 10 \times 5 = 61.25\ \text{N/m}$$

即每米展长产生 61.25 N 的升力，约等于 6.25 kgf/m。

验证环量对速度分布的影响：

在圆柱顶部（$\theta = 90^\circ$）：$V_\theta = -2V_\infty(1) - \Gamma/(2\pi R)$

$$V_\theta|_{\text{top}} = -2\times 10 - \frac{5}{2\pi\times 0.1} \approx -20 - 7.96 = -27.96\ \text{m/s}$$

在圆柱底部（$\theta = -90^\circ$）：$V_\theta = -2V_\infty(-1) - \Gamma/(2\pi R)$

$$V_\theta|_{\text{bottom}} = 20 - 7.96 = 12.04\ \text{m/s}$$

顶部速度 27.96 m/s（远快于来流），底部速度 12.04 m/s（略快于来流但比顶部慢）。顶部的低压区与底部的高压区之间的压差，产生了向上的净力——升力。

---

三、库塔条件——环量的大小由谁决定

3.1 问题

库塔-儒科夫斯基定理的形式之美在于：$L' = \rho V_\infty \Gamma$。

但它立刻引出一个问题：环量 $\Gamma$ 是多少？

对于圆柱，$\Gamma$ 可以由我们任意指定（这是数学上的自由）。但对于一个尖后缘的翼型，大自然给出了一个确定的答案。

3.2 库塔条件的表述

库塔条件（Kutta Condition）：

对于一个有尖后缘的翼型，在给定迎角下，流场会自动调整环量，使得后缘处的速度是有限的。更具体地，上表面和下表面的气流平滑地在后缘汇合——后缘不再是奇点。

数学表达：

$$\boxed{V_{\text{upper, TE}} = V_{\text{lower, TE}}}$$

即在翼型后缘处，上表面和下表面的速度大小相等、方向一致。

3.3 物理直觉

为什么库塔条件成立？

想象一个翼型从静止开始加速到运动速度 $V_\infty$。在加速的瞬间，流体在翼型后缘从下表面绕到上表面——这本质上是黏性在起作用，会形成一个起动涡（starting vortex）。

根据开尔文环量定理，总环量必须守恒（初始为0），所以翼型周围必然产生一个与起动涡大小相等、方向相反的附着涡（bound vortex）——这就是翼型的环量。

这个过程 Anderson 在书中用实验照片生动地展示了——翼型在起动时，后缘真的会”释放”出一个旋涡。

3.4 应用不可压流动

对于薄翼型，库塔条件加上势流理论，可以形成完整的升力预测：

1. 翼型后缘处，满足库塔条件
2. 由此确定附着涡的强度 $\Gamma$
3. 用库塔-儒科夫斯基定理：$L' = \rho V_\infty \Gamma$

这个链条是第5讲（薄翼理论）的核心基础。

---

四、保角变换——从圆柱到翼型的数学桥梁

4.1 概念引入

如果我们能够精确求解圆柱绕流，但是想要研究翼型绕流——怎么办？

答案是一种精巧的数学工具：保角变换（Conformal Mapping）。它的基本思想是：

通过一个解析变换 $w = f(z)$，把物理平面（$z$ 平面）上复杂的翼型形状映射到计算平面（$\zeta$ 平面）上一个简单的圆。

在圆上我们已经知道了全部解（上节内容），然后反变换回去，就得到了翼型绕流的精确解。

4.2 儒科夫斯基变换

最著名的保角变换是儒科夫斯基变换（Zhukovsky Transform，以俄罗斯空气动力学先驱命名）：

$$\zeta = z + \frac{c^2}{z}$$

其中 $c$ 是变换常数。

这个变换有这样的特性：

• $z$ 平面上的圆 $\rightarrow$ $\zeta$ 平面上的**儒科夫斯基翼型**
• 圆上 $\theta = 0$（最后缘）的点被映射到翼型后缘
• 圆上 $\theta = \pi$（最前缘）的点被映射到翼型前缘

借助这个变换，空气动力学理论可以精确地计算出任意儒科夫斯基翼型的升力系数！

4.3 升力系数的理论值

通过保角变换得到的理论升力系数公式：

$$\boxed{C_l = 2\pi\alpha}$$

其中 $\alpha$ 是迎角（单位：弧度）。

对二维薄翼型，这个公式预测升力系数与迎角成正比，斜率为 $2\pi$（每弧度约 0.1096/度）。

这是一个极其重要的结果——它意味着：

• 迎角增加 1° → $C_l$ 增加约 0.11
• 迎角 5° → $C_l \approx 0.55$
• 迎角 10° → $C_l \approx 1.1$

当然，这个理论值只适用于未失速的小迎角范围。在第5讲薄翼理论中，我们会看到这个 $2\pi$ 斜率的更严格推导。

4.4 手算例子：一个薄翼的升力

已知：巡航迎角 $\alpha = 4^\circ$，弦长 $c = 0.3\ \text{m}$，来流速度 $V_\infty = 15\ \text{m/s}$。

求：单位展长升力。

解：

$$\alpha = 4^\circ = 4 \times \frac{\pi}{180} = 0.0698\ \text{rad}$$ $$C_l = 2\pi\alpha = 2\pi \times 0.0698 \approx 0.438$$

升力：

$$L' = C_l \cdot \frac{1}{2}\rho V_\infty^2 \cdot c = 0.438 \times \frac12 \times 1.225 \times 15^2 \times 0.3$$ $$L' = 0.438 \times 0.5 \times 1.225 \times 225 \times 0.3 \approx 0.438 \times 41.34 \approx 18.1\ \text{N/m}$$

物理直觉：对于一个展长 2 m 的机翼，总升力为 36.2 N，约 3.7 kgf——刚好支撑一个小型无人机的重量。

---

五、库塔-儒科夫斯基定理的推广

5.1 一般形式

库塔-儒科夫斯基定理不仅仅适用于圆柱或翼型，它是二维势流中最普遍的升力定理：

在不可压缩、无黏、定常的势流中，任何柱体所受的升力都与绕柱体的环量成正比：

$$\boxed{L' = \rho V_\infty \Gamma}$$

其中 $\Gamma$ 是绕物体的环量（沿包围物体的闭合路径积分），与路径的选取无关。

5.2 升力与阻力的方向

库塔-儒科夫斯基定理只给出了升力，没有阻力。这符合达朗贝尔佯谬——势流中确实没有阻力。

一个重要推论：在势流中，阻力为零；升力等于环量乘以密度和流速。

5.3 力和力矩

除了升力，翼型还有力矩（pitching moment）：

$$M'_{\text{LE}} = -\frac{1}{2}\rho V_\infty^2 c^2 \pi \alpha$$

其中负号表示低头力矩——这正是为什么需要水平尾翼（或前置鸭翼）来配平。

---

六、完整概念地图

达朗贝尔佯谬
    ↓
圆柱绕流: 均匀流 + 偶极子
    ↓ (对称, 零升力 + 零阻力)
    ↓
需要非对称 → 加一个点涡(环量Γ)
    ↓
带环量的圆柱绕流
    ↓
压力积分 → ρV∞Γ = 升力
    ↓
库塔-儒科夫斯基定理: L' = ρV∞Γ
    ↓
问题: Γ 等于多少?
    ↓
库塔条件: 尖后缘 → 流速有限 → 确定Γ
    ↓
保角变换: 圆 → 翼型
    ↓
理论结果: Cl = 2πα
    ↓
[真实翼型的Cl与Cd曲线]

---

七、各概念在已发表文章中的出现

概念	出现文章	使用方式
环量	《前进比完全解读》	间接涉及——桨叶截面环量与推力
升力线斜率 $2\pi$	《四旋翼飞行动力学建模》	简易升力模型参数
库塔条件	—	本文首次出现
保角变换	—	本文首次出现

---

八、核心公式速查卡

公式	含义	备注
$V_\theta|_{r=R} = -2V_\infty\sin\theta - \frac{\Gamma}{2\pi R}$	圆柱表面速度分布	环量使分布不对称
$L' = \rho V_\infty \Gamma$	库塔-儒科夫斯基定理	二维升力公式
$C_l = 2\pi\alpha$	薄翼升力系数理论值	$\alpha$ 用弧度制
$V_{\text{upper,TE}} = V_{\text{lower,TE}}$	库塔条件	后缘光滑流动
$\zeta = z + c^2/z$	儒科夫斯基变换	圆→翼型的映射

---

参考文献

1. Anderson, J. D. (2010). Fundamentals of Aerodynamics (6th ed., Chapters 3-4). McGraw-Hill.
2. Kutta, W. M. (1902). “Auftriebskräfte in strömenden Flüssigkeiten”. Illustrierte Aeronautische Mitteilungen, 6, 133-135.
3. Joukowski, N. E. (1906). “On the Contours of the Supporting Surfaces of Aeroplanes”. Bulletin de l’Institut Aérodynamique de Koutchino, 1.
4. Milne-Thomson, L. M. (1996). Theoretical Hydrodynamics (5th ed., Chapter 9). Dover.
5. Prandtl, L., & Tietjens, O. G. (1934). Applied Hydro- and Aeromechanics (Chapter 4). Dover.

---

下一节：翼型空气动力学与薄翼理论

从真实翼型的几何定义（中弧线、弦线、厚度分布、弯度、迎角）出发，深入推导薄翼理论：用涡面法将翼型离散为一系列附着涡，通过求解积分方程得到环量分布，进而得到升力系数、力矩系数和压力中心。最后讨论翼型的失速特性和边界层分离机制。