Anderson《空气动力学基础》读书笔记（一）：流体力学基础——空气动力学的语言

发表于2026-04-14|更新于2026-06-02|空气动力学基础

|浏览量:

本文是 John D. Anderson《Fundamentals of Aerodynamics》（第六版）读书笔记系列的第一篇，覆盖教材第一章（空气动力学概述）与第二章（流体力学基础）的核心内容。阅读本文无需流体力学基础，所有概念从零定义。

第一部分：空气动力学是什么——Anderson 的视角

一、空气动力学的研究对象

1.1 概念引入

你站在室外，一阵风吹来——风对你施加了力。一辆汽车在高速公路上开到 120 km/h，阻力会消耗近一半的燃油。一架波音 787 起飞时，机翼上下的压力差托起了 250 吨的机身。

这三个现象的背后是同一个科学问题：

物体在空气（或其他气体）中运动时，气体对物体施加了多大的力？这些力如何产生？如何预测和控制？

这就是空气动力学回答的核心问题。

Anderson 在教材开篇给出了一个简洁的定义：

空气动力学（Aerodynamics）是研究物体在气体中运动时气体对物体施加的力与力矩的科学。

这个定义听起来简单，但其内涵贯穿了整个航空航天工程：升力有多高？阻力有多大？翼型怎么设计才能让升力大、阻力小？超声速飞行时为什么会产生激波？

1.2 历史的五分钟

Anderson 讲述空气动力学史的方式很有特点——他不只是罗列时间点，而是通过**五大”空气动力学里程碑时刻”**来串联整个学科的发展：

时间	人物	贡献	意义
1687	牛顿（Isaac Newton）	提出阻力与速度平方成正比的”冲击理论”	最早的空气动力学理论模型
1738	伯努利（Daniel Bernoulli）	发表《流体动力学》，提出伯努利方程	从能量守恒理解流速与压力关系
1745	达朗贝尔（Jean le Rond d’Alembert）	发现”达朗贝尔佯谬”：理论计算的阻力为零	迫使学界承认黏性（摩擦）的关键作用
1752	皮托（Henri Pitot）	发明皮托管	有了测量气体速度的实用工具
1903	莱特兄弟	首次动力飞行	第一次在工程上验证了空气动力学理论

这五个里程碑构成了空气动力学的前理论时代。真正的理论突破要到 1904 年普朗特（Ludwig Prandtl）提出边界层理论，空气动力学才成为一门定量的工程科学。

1.3 空气动力学在无人机仿真中的出现

在已发表的博客文章中，空气动力学的核心概念频繁出现：

文章	使用的空气动力学概念
前进比完全解读：从螺旋桨空气动力学到高速四旋翼设计	推力系数 $C_T$ 、前进比 $J$ 、推力-阻力平衡、失速
四旋翼飞行力学基础	螺旋桨升力产生、反扭矩
四旋翼飞行动力学建模——AirSim源码深度解析	推力公式 $T=C_T\rho\omega^2R^4$ 、阻力、空气密度
截击无人机核心算法深度解析	空气动力学对制导与控制的约束
可微物理引擎完全解析	仿真器中的空气动力学模型（CFD 减少/代理模型）

二、流体力学的基本哲学——连续介质假设

2.1 概念引入

空气由分子构成——分子之间有间隙。但如果我们试图用每个分子的运动来描述空气流过一架飞机，那我们需要跟踪 $10^{44}$ 个分子（标准状态下 $1\text{ m}^3$ 空气约含 $2.5\times10^{25}$ 个分子）。这显然是不可能完成的任务。

空气动力学的核心假设——连续介质假设——解决了这个问题：

我们把空气看作一个”连续的介质”，忽略分子之间的间隙，认为流体在空间中的每一点都有定义的密度 $\rho$ 、压力 $p$ 、速度 $\mathbf{V}$ 和温度 $T$ ，而不需要考虑单个分子的行为。

这四个量 $\{\rho, p, \mathbf{V}, T\}$ 称为流场的基本变量。空气动力学的全部任务，就是求解这些变量在时空中的分布。

2.2 定义

连续介质假设的形式化表述：

对于任意一个宏观上足够小（远小于流动特征尺寸）、微观上足够大（远大于分子自由程）的微元体，该微元体内部的物理量（密度、压力、速度、温度）可以用连续函数来描述。

这个假设的适用范围有一个明确的判据——克努森数（Knudsen number）：

Kn = \frac{\lambda}{L}

其中 $\lambda$ 是气体分子的平均自由程（分子两次碰撞之间飞过的平均距离）， $L$ 是流动的特征尺寸（如翼弦长度、管径等）。

$Kn < 0.01$ ：连续介质假设成立（绝大多数航空航天问题满足）
$0.01 < Kn < 0.1$ ：滑移流区（高超声速稀薄气体）
$Kn > 0.1$ ：自由分子流区（需要分子动力学）

标准状态下：

海平面空气的 $\lambda \approx 6.5\times 10^{-8}\text{ m}$ （65 nm）
翼根弦长 $L \approx 0.1\text{ m}$ （以四旋翼桨叶为例）
$Kn \approx 6.5\times 10^{-7} \ll 0.01$ → 连续介质假设完美成立

2.3 性质

宏观可测性：连续介质假设下的 $\rho, p, \mathbf{V}, T$ 都可以用宏观仪器（皮托管、温度计、热风速仪）直接测量。
数学可微性：流场变量是空间坐标 $(x,y,z)$ 和时间 $t$ 的连续（通常还可微）函数，这使得我们可以用偏微分方程来描述流体运动。
适用范围明确：在导弹再入大气层（80 km 以上）、高海拔无人机（30 km 以上微小型）等场景下， $Kn$ 增大，连续介质假设可能失效。

2.4 表示方法

在空气动力学中，流场变量通常表示为：

变量	符号	单位（SI）	物理意义
密度	$\rho$	$\text{kg/m}^3$	单位体积的空气质量
压力	$p$	$\text{N/m}^2$ 或 $\text{Pa}$	单位面积的气体法向力
速度	$\mathbf{V}$	$\text{m/s}$	气体的宏观运动速度（矢量）
温度	$T$	$\text{K}$	气体分子热运动强度的度量（开尔文）

这四个变量并非独立——它们通过状态方程关联：

p = \rho R T

其中 $R$ 是气体常数。对于空气， $R = 287\ \text{J/(kg·K)}$ 。这个方程是空气动力学中连接热力学和流体动力学的一座桥梁。

2.5 手算例子：海平面空气状态

已知：海平面标准条件下， $p = 101325\ \text{Pa}$ ， $T = 288.15\ \text{K}$ （15°C）。

求：空气密度 $\rho$ 。

解：

\rho = \frac{p}{RT} = \frac{101325}{287 \times 288.15} \approx 1.225\ \text{kg/m}^3

这个值—— $1.225\ \text{kg/m}^3$ ——是整个空气动力学中最常见的常数之一。记住它：海平面标准空气密度就是 $1.225\ \text{kg/m}^3$ 。

物理直觉： $1\ \text{m}^3$ 空气大约重 $1.225\ \text{kg}$ ，差不多是一瓶大可乐的重量。这个”重量”就是推动机翼产生升力的物质基础——没有质量的空气是无法产生力的。

三、速度场与流线——描述流动的语言

3.1 概念引入

如果你对着水面撒一把细沙，沙粒的轨迹会”画出”水流的路径。同样地，在风洞中从壁面喷出烟雾，烟雾条纹会”画出”空气的流动轨迹。

我们需要一种方式来描述流体的运动。Anderson 引入了两种视角：

欧拉（Euler）视角：站在一个固定点，观察流体流过这个点。就像河边的人看水流过桥墩。
拉格朗日（Lagrange）视角：跟着一个特定的流体”微团”（parcel）走，观察它在运动中的变化。就像在河上放一个漂流瓶，跟着它漂。

3.2 定义

速度场（Velocity Field）：

在欧拉视角下，流场中每一点 $\mathbf{r} = (x, y, z)$ 在每一时刻 $t$ 都有一个速度矢量 $\mathbf{V}(x,y,z,t)$ 。这个矢量函数给出了整个流场的完整信息：

\mathbf{V} = \mathbf{V}(x,y,z,t)

流线（Streamline）：

流线是这样一个曲线，它在每一点的切线方向与该点的速度方向相同。

数学上，流线满足：

\left.\frac{dy}{dx}\right|_{z=\text{const}} = \frac{v}{u}, \quad \left.\frac{dz}{dx}\right|_{y=\text{const}} = \frac{w}{u}

其中 $(u,v,w)$ 是速度 $\mathbf{V}$ 的三个分量。

流线与速度场的区别：

瞬时的快照：流线不随时间变化（定常流）或反映某个瞬间的形状
路径线（Pathline）：单个流体微团的实际运动轨迹（拉格朗日视角）
条纹线（Streakline）：从同一点连续释放的所有流体微团在某个时刻的位置——风洞烟雾实验看到的实际就是条纹线

对于定常流动（流场不随时间变化），这三种线完全重合。

3.3 性质

流线不相交：除非在速度为 0 的”驻点”，否则两条流线不会相交。
流线密的地方流速快：流线间距反映了流速场的局部信息——间距缩小意味着流速增大（后续会解释为什么）。
固体壁面本身就是一条流线：气流不会穿过固体壁面，所以壁面处的切向速度方向与壁面平行。

3.4 应用：为什么四旋翼的桨叶是”越往外越宽”？

在《前进比完全解读》文章中我们讨论过：桨叶根部的叶素速度较低（ $\omega \cdot r$ 小），而叶尖速度最高。从流管的角度看，这相当于”流管截面积变化”导致速度分布变化——这是连续的流动结构。

作为反面例子：如果四旋翼的桨叶在高速前进时处于过大迎角，流线在桨叶上表面分离，推力急剧下降——这就是失速。失速的本质就是流线不再贴着翼型表面行走。

四、流体运动的基本方程——控制体积法

4.1 概念引入

Anderson 在第二章采用了控制体积方法来推导流体运动方程。这种方法的哲学是：

我们把空间中的一块固定区域（控制体积，Control Volume）作为考察对象，观察什么”流进”、什么”流出”，以及内部的积累或损失。

就像你数一个水池的进出水量：流入的水量减去流出的水量，等于水池中水量的变化。这就是守恒定律在流体力学中的表现。

空气动力学的基本方程就是三条守恒定律在流体上的具体表达：

质量守恒 → 连续方程
动量守恒（牛顿第二定律） → 动量方程（Euler/Navier-Stokes）
能量守恒（热力学第一定律） → 能量方程

4.2 定义：连续方程

物理表述：

控制体积内质量的变化率 = 通过控制体积表面净流入的质量流量

数学表达（积分形式）：

\frac{\partial}{\partial t}\iiint_V \rho\, dV + \iint_S \rho\mathbf{V}\cdot d\mathbf{S} = 0

这个方程的意思是：

第一项：控制体积内部质量的变化率
第二项：通过控制体积边界净流出的质量流量
两者之和等于零——质量不生不灭

对定常流动的简化：

如果流动不随时间变化（定常流动）， $\frac{\partial}{\partial t}=0$ ，则：

\iint_S \rho\mathbf{V}\cdot d\mathbf{S} = 0

即流入控制体积的质量流量 = 流出控制体积的质量流量。

对于一维管流，这个方程简化为著名的流管方程：

\rho_1 V_1 A_1 = \rho_2 V_2 A_2

其中 $A$ 是流管的横截面积。

4.3 定义：动量方程

物理表述（牛顿第二定律）：

控制体积内流体动量的变化率 = 作用在流体上的合力

如果忽略黏性（无黏流动，即 Euler 方程），则只考虑两种力：

体积力（Body Force）：如重力
表面力（Surface Force）：压力

一维定常、忽略重力的欧拉方程是最简单也最重要的形式：

p + \frac{1}{2}\rho V^2 = \text{constant} \quad \text{(Bernoulli 方程)}

这就是伯努利方程——空气动力学中最出名的一个方程。注意它的三个限制条件：

无黏（不考虑摩擦）
定常（不随时间变化）
不可压缩（密度为常数）

4.4 伯努利方程——用手算出升力

伯努利方程之所以如此重要，是因为它直接给出了压力与速度的定量关系：

流速快的地方压力低，流速慢的地方压力高。

这正是产生”升力”的核心机制。来看一个典型应用——估算机翼升力。

手算例子：

一架四旋翼无人机在悬停时，旋翼平均下洗速度 $V = 10\ \text{m/s}$ ，旋翼上方自由流速度近似为 0，压力为当地大气压 $p_\infty$ 。

求：旋翼上方的压力变化。

解：

在自由流处： $p_\infty + \frac{1}{2}\rho (0)^2 = \text{constant}$

在旋翼处： $p + \frac{1}{2}\rho V^2 = p_\infty$

p = p_\infty - \frac{1}{2}\rho V^2 = p_\infty - \frac{1}{2} \times 1.225 \times 10^2

p - p_\infty = -61.25\ \text{Pa}

即旋翼上方形成了一个 61.25 Pa 的低压区。如果旋翼盘面积为 $0.1\ \text{m}^2$ ，则升力为 $61.25 \times 0.1 \times 2 \approx 12.25\ \text{N}$ （乘以 2 是因为上下表面都有压差）。这个力恰好支撑 1.25 kg 左右的四旋翼——这正是很多小型四旋翼的典型重量。

4.5 性质对比：三种方程

方程	物理原理	数学形式（一维定常）	工程应用
连续方程	质量守恒	$\rho_1 V_1 A_1 = \rho_2 V_2 A_2$	喷管设计、流管分析
动量方程	$F=ma$	$p + \frac12\rho V^2 = \text{const}$	升力/阻力计算、管道设计
能量方程	能量守恒	$h + \frac12 V^2 = \text{const}$	可压缩流、激波计算

4.6 历史：从 Euler 到 Bernoulli

经典错误纠正：很多人以为伯努利方程是由 Daniel Bernoulli 独立完成的。事实是：

1738年：Daniel Bernoulli 出版了《流体动力学》（Hydrodynamica），提出了能量形式的方程
1749年：Leonhard Euler 给出了伯努利方程更严格的数学推导，并建立了更一般的 Euler 方程

实际上，我们今天使用的 $p + \frac12\rho V^2 = \text{const}$ 形式主要来自 Euler 的数学工作，Bernoulli 的原著使用的表达方式不同。Anderson 在书中的注释明确指出了这一点。

五、可压缩与不可压缩——什么时候密度不再是常数

5.1 概念引入

如果我用针筒压缩空气，活塞被推向里——空气体积缩小了，密度增加了。但如果我用水做同样的事，水几乎不会被压缩。

空气本来就是可压缩的。但在低速流动中，压力变化很小，密度变化可以忽略。这引出了空气动力学的一个基本分类：

流动类型	密度变化	典型场景	马赫数范围
不可压缩	可忽略（ $\rho\approx\text{const}$ ）	低速飞机、四旋翼	$Ma < 0.3$
可压缩	显著	客机巡航、超声速飞行	$Ma > 0.3$

5.2 定义

马赫数（Mach Number） 是判断可压缩性是否起作用的标尺：

Ma = \frac{V}{a}

其中 $a$ 是声速。在海平面标准大气下， $a \approx 340\ \text{m/s}$ （约 1224 km/h）。

为什么马赫数能标记可压缩性？因为密度变化的百分比大致正比于 $Ma^2$ ：

\frac{\Delta\rho}{\rho} \propto Ma^2

所以：

$Ma = 0.3$ 时，密度变化约 4.5% → 不可压缩假设误差可接受
$Ma = 0.8$ 时，密度变化约 32% → 必须考虑可压缩
$Ma = 2.0$ 时，密度变化超过 200% → 完全进入可压缩领域

5.3 四旋翼的马赫数：一个被低估的问题

很多四旋翼设计者可能没有意识到——四旋翼桨尖的马赫数经常超过 0.6。

以常见的 7 英寸（177.8 mm）桨为例：

手算例子：

若转速 $\Omega = 10000\ \text{RPM}$ ，桨尖速度：

V_{\text{tip}} = \Omega \times R = \frac{10000}{60} \times 2\pi \times 0.0889 \approx 93\ \text{m/s}

马赫数：

Ma_{\text{tip}} = \frac{93}{340} \approx 0.27

这个速度还处于不可压缩范围。但如果是高速四旋翼截击机，情况完全不同。

高速场景（以 Thunder STD100 为例，巡航 350 km/h ≈ 97 m/s）：

前向飞行时，前进桨叶的桨尖速度是旋转速度和前飞速度的矢量叠加：

V_{\text{tip,adv}} = \omega R + V_\infty = 93 + 97 = 190\ \text{m/s}

Ma_{\text{tip,adv}} = \frac{190}{340} \approx 0.56

这个马赫数已经进入了可压缩区域（ $Ma > 0.3$ ），桨尖处会出现可压缩性效应——推力系数下降、阻力系数上升。这就是为什么很多高速无人机在加速到某个速度后，推力提升变得非常低效。

在《前进比完全解读》文章中，我们讨论了前进比 $J = V_\infty/(nD)$ 对推力和效率的影响。实际上， $J$ 和 $Ma$ 是等价的参数——在给定桨径和转速下，两者都正比于飞行速度。Anderson 教材进一步告诉我们：马赫数效应将改变整个推力-空气动力学结构。

六、黏性——空气动力学中”不完美”的美

6.1 概念引入

达朗贝尔在 1745 年发现了一个令人困惑的结果：如果用无黏理论计算一个物体的阻力，结果是——零。

但现实世界中，任何物体在空气中运动都有阻力。为什么？

答案是黏性（Viscosity）。空气是有”粘稠度”的——虽然很小，但不可忽略。黏性使得紧贴着物体表面的空气”粘”在表面不动（无滑移条件），形成了一个速度从表面零到外界自由流速度的过渡区域——这就是边界层（Boundary Layer）。

6.2 定义

黏性系数（Dynamic Viscosity） $\mu$ 衡量流体的”粘稠程度”。空气的 $\mu$ 在标准状态下约为 $1.79\times 10^{-5}\ \text{Pa·s}$ ——比水的黏性低约 50 倍。

Sutherland 公式（Anderson 书中给出的经验关系）：

\frac{\mu}{\mu_0} = \left(\frac{T}{T_0}\right)^{3/2} \frac{T_0 + S}{T + S}

其中对空气： $\mu_0 = 1.7894\times 10^{-5}\ \text{Pa·s}$ ， $T_0 = 288.15\ \text{K}$ ， $S = 110.4\ \text{K}$ 。

这个公式告诉我们：温度升高，空气的黏性增大——与直觉相反，热空气比冷空气”更粘稠”。

运动黏性系数（Kinematic Viscosity） $\nu$ ：

\nu = \frac{\mu}{\rho}

标准状态下， $\nu \approx 1.46\times 10^{-5}\ \text{m}^2/\text{s}$ 。

6.3 区分流动性质的关键：雷诺数

雷诺数（Reynolds Number） 是空气动力学中最重要的无量纲数（没有之一）：

Re = \frac{\rho V L}{\mu} = \frac{V L}{\nu}

其中 $L$ 是特征长度（翼弦长度、球体直径等）。

雷诺数的物理意义：惯性力与黏性力的比值。

$Re$ 很大（ $> 10^5$ ）：惯性力主导，流动大多是湍流
$Re$ 很小（ $< 1$ ，如微型昆虫飞行）：黏性力主导，流动是层流

6.4 手算例子：两种完全不同尺度的流动

例1：波音 787 机翼

V = 250\ \text{m/s}

（巡航），

L = 6\ \text{m}

（平均气动弦长）

Re = \frac{1.225 \times 250 \times 6}{1.79\times 10^{-5}} \approx 1.03 \times 10^8

Re = 10^8

→ 几乎完全湍流边界层。

例2：四旋翼桨叶

V = 10\ \text{m/s}

（h 工况），

L = 0.02\ \text{m}

（桨叶弦长）

Re = \frac{1.225 \times 10 \times 0.02}{1.79\times 10^{-5}} \approx 1.37 \times 10^4

Re = 1.4 \times 10^4

→ 层流边界层，非常容易失速。

这两个例子解释了为什么：

大型飞机可以在大迎角下安全飞行（高雷诺数→湍流→延迟分离）
四旋翼在突然加速或风中倾斜时容易推力崩溃（低雷诺数→层流分离→失速）

七、完整概念地图

连续介质假设
     ↓
流场变量 {ρ, p, V, T}
     ↓
    ┌────────────────────────────────┐
    │    描述流动的两种方法            │
    ├──────────────┬─────────────────┤
    │ 欧拉视角      │ 拉格朗日视角      │
    │ (速度场 流线)  │ (路径线)          │
    └──────────────┴─────────────────┘
     ↓
    ┌────────────────────────────────┐
    │    三大守恒方程                  │
    ├──────────────┬─────────────────┤
    │ 质量守恒      │ 动量守恒          │
    │ (连续方程)    │ (欧拉/N-S方程)    │
    └──────────────┴─────────────────┘
     ↓                      ↓
    伯努利方程        黏性 → 边界层
    (压力-速度关系)          ↓
                        雷诺数
     ↓
    ┌────────────────────────────────┐
    │    可压缩性的标尺                │
    │    马赫数 Ma                    │
    │ Ma<0.3: 不可压缩                │
    │ 0.3<Ma<0.8: 亚声速可压缩         │
    │ 0.8<Ma<1.2: 跨声速               │
    │ 1.2<Ma<5.0: 超声速               │
    │ Ma>5.0: 高超声速                  │
    └────────────────────────────────┘

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

概念	出现文章	使用方式
伯努利方程	《前进比完全解读》	解释桨叶截面压力分布
连续方程	《四旋翼飞行动力学建模》	推导旋翼诱导速度
马赫数	《前进比完全解读》	讨论桨尖可压缩效应
雷诺数	《空气动力学在仿真中的验证》	解释边界层转捩
黏性	《可微物理引擎完全解析》	代理模型中的阻力项

九、核心公式速查卡

公式	含义	适用条件
$p = \rho RT$	气体状态方程	热力学平衡
$\frac{\partial}{\partial t}\iiint_V \rho dV + \iint_S \rho\mathbf{V}\cdot d\mathbf{S} = 0$	连续方程（一般形式）	一切流动
$\rho_1 V_1 A_1 = \rho_2 V_2 A_2$	连续方程（一维定常）	管流
$p + \frac12\rho V^2 = \text{const}$	伯努利方程	无黏、定常、不可压缩
$Ma = V/a$	马赫数	一切流动
$a = \sqrt{\gamma R T}$	声速	理想气体
$Re = \rho V L / \mu$	雷诺数	一切流动
$\mu/\mu_0 = (T/T_0)^{3/2}(T_0+S)/(T+S)$	Sutherland 公式	空气黏性温度依赖

十、给自学者的建议

如果你想继续深入学习空气动力学，建议按以下顺序阅读 Anderson 教材：

第1-2章（本篇覆盖）：打好基础，建立空气动力学的物理直觉
第3-4章（势流理论）：真正理解升力产生的数学本质
第5-6章（翼型理论）：从理论到实用的桥梁
第7-8章（有限翼理论）：真正三维机翼的分析
第9章（可压缩流）：超声速世界的入口
第12章（边界层）：理解阻力来源

每读一章，对照本章中的公式速查卡做 2-3 个手算练习，理解每个参数的量级。

参考文献

Anderson, J. D. (2010). Fundamentals of Aerodynamics (6th ed.). McGraw-Hill. ISBN: 978-0-07-339810-5.
Anderson, J. D. (2016). Introduction to Flight (8th ed.). McGraw-Hill. ISBN: 978-0-07-802767-3.
Prandtl, L. (1904). “Über Flüssigkeitsbewegung bei sehr kleiner Reibung”. Verhandlungen des III. Internationalen Mathematiker-Kongresses, Heidelberg.
Newton, I. (1687). Philosophiæ Naturalis Principia Mathematica. London. Book II, Section VII.
Euler, L. (1757). “Principles généraux du mouvement des fluides”. Mémoires de l’Académie des Sciences de Berlin, 11, 274-315.

下一节：势流理论——从拉普拉斯方程到基本流动解