本文是 John D. Anderson《Fundamentals of Aerodynamics》(第六版)读书笔记系列的第二篇,覆盖教材第3章《势流理论》的前半部分。阅读本文前,建议先了解流体力学基础中的连续方程和伯努利方程的基本概念。
第二部分:势流理论——理解升力的数学语言
在第一篇中我们提到:如果忽略黏性,流体的动量方程简化为欧拉方程;进一步,如果流动无旋(irrotational),我们就可以引入一个叫做速度势的函数 ϕ,使得流场的所有信息浓缩在一个标量拉普拉斯方程中。这就是势流理论的起点。
势流理论在 Anderson 的书中占据了第3-4章,是全书数学最密集的部分之一,也是理解升力产生的数学本质的必经之路。
一、涡量与无旋流动——为什么能”简化”?
1.1 概念引入
想象一个河面上的小水车。如果水车在流动中旋转,就说明这个地方的水有”旋涡”;如果水车平移而不旋转,就说明这里的水是”无旋”的。
这个区别非常重要。有旋流动和无旋流动,在数学上完全是两个世界:
- 有旋流动:需要用完整的 Navier-Stokes 方程,四个方程解四个未知数,数学上极其复杂
- 无旋流动:可以引入速度势,化简为一个标量拉普拉斯方程,数学上大幅简化
而幸运的是,许多工程上重要的流动——包括机翼周围的附着流动——在很大的程度上可以视为无旋流动。
1.2 定义
一个流动是否”有旋”,由**涡量 ω(Vorticity)**来度量:
ω=∇×V=i∂x∂uj∂y∂vk∂z∂w
三维展开:
ω=(∂y∂w−∂z∂v)i+(∂z∂u−∂x∂w)j+(∂x∂v−∂y∂u)k
无旋流动的定义:ω=0(流场每一点的涡量为零)。
在二维流动中,这个条件简化为:
∂x∂v−∂y∂u=0
1.3 性质
- 无旋 ≠ 无涡:无旋流动并不意味着流线是直的。绕圆柱的流动可以是无旋的,但流线是弯曲的——“无旋”指的是流体微团不绕自身旋转,而不是流线不弯曲。
- 无旋具有”记忆”:如果均匀来流是无旋的,它流过光滑物体后仍保持无旋——这个性质称为开尔文环量定理(Kelvin’s circulation theorem)。
- 边界层内部是有旋的,外部是无旋的:黏性只在边界层内产生旋涡,边界层以外的势流区域可以视为无旋流动。这就是为什么势流理论仍然实用的根本原因。
1.4 历史
| 时间 |
人物 |
贡献 |
| 1781 |
拉普拉斯(Pierre-Simon Laplace) |
建立了拉普拉斯方程 ∇2ϕ=0 |
| 1815 |
拉格朗日(Joseph-Louis Lagrange) |
证明了无旋流动存在速度势 |
| 1858 |
亥姆霍兹(Hermann von Helmholtz) |
建立了涡量动力学,证明”涡管强度守恒” |
| 1867 |
开尔文(Lord Kelvin) |
证明环量定理——无旋流动保持无旋 |
| 1902 |
兰金(William John Macquorn Rankine) |
提出了半无限体绕流的组合解 |
1.5 应用
在博客文章《四旋翼飞行动力学建模——AirSim源码深度解析》中,旋翼的诱导速度模型假设旋翼远场流动是无旋的,因此在动量理论中使用了伯努利方程来描述旋翼盘前后的压力差。实际上这个假设在中等精度下成立,但在旋翼盘附近(近场)流动是有旋的——这就是为什么动量理论需要加上”涡流因子”来修正。
二、速度势——把向量问题变成标量问题
2.1 概念引入
一个三维流场 V(x,y,z) 是一个矢量函数——有三个分量 (u,v,w),需要三个方程才能解。但如果我们能找到一个标量函数 ϕ(x,y,z),使得 V=∇ϕ,那么我们就可以用一个标量函数代替三个分量函数。
这是空气动力学中最重要的”降维”技巧。
2.2 定义
速度势(Velocity Potential)ϕ 定义为满足以下关系的函数:
V=∇ϕ
在笛卡尔坐标系下展开:
u=∂x∂ϕ,v=∂y∂ϕ,w=∂z∂ϕ
在极坐标系下(二维):
Vr=∂r∂ϕ,Vθ=r1∂θ∂ϕ
存在性定理:流动无旋(∇×V=0)⟺ 存在速度势 ϕ 使得 V=∇ϕ。
这个定理来自向量微积分的基本定理:旋度为零的向量场一定是某个标量函数的梯度。
2.3 从速度势到拉普拉斯方程
将速度势代入连续方程(不可压缩):
∇⋅V=0
∇⋅(∇ϕ)=0
这就得到了势流理论中最核心的方程——拉普拉斯方程:
∇2ϕ=0
二维展开:
∂x2∂2ϕ+∂y2∂2ϕ=0
这是一个二阶级线性偏微分方程。它的重要意义在于:
- 线性:解可以叠加——复杂流动 = 基本流动解的线性组合
- 成熟的理论:拉普拉斯方程的数学理论在 19 世纪就已经相当完善
- 一个约束明确边值问题:只需要在物体表面施加一个边界条件(∂ϕ/∂n=0,即壁面无穿透条件),就可以唯一确定流场
2.4 性质
速度势的物理含义:
- ϕ=const 的等值面称为**等势面**。速度方向总是垂直于等势面。
- ϕ 的大小本身没有物理意义,有意义的是**它的梯度**(速度)。
- 在不可压缩无旋流动中,ϕ 满足拉普拉斯方程——这是整个势流理论的数学基础。
对比速度势与流函数(流函数将在下节介绍):
| 特性 |
速度势 ϕ |
流函数 ψ |
| 定义 |
V=∇ϕ |
u=∂y∂ψ,v=−∂x∂ψ |
| 存在条件 |
无旋(∇×V=0) |
不可压缩(∇⋅V=0) |
| 基本方程 |
∇2ϕ=0(拉普拉斯方程) |
∇2ψ=0(拉普拉斯方程) |
| 等值线意义 |
等势线 |
流线 |
| 物理意义 |
速度的”势” |
流量 |
2.5 手算例子:一维均匀流的速度势
最简单的势流——均匀流 u=V∞,v=0。
验证无旋:
∂x∂v−∂y∂u=0−0=0⇒无旋条件满足
求速度势:
u=∂x∂ϕ=V∞⇒ϕ=V∞x+f(y)
v=∂y∂ϕ=f′(y)=0⇒f(y)=const
取常数项为零,则:
ϕ=V∞x
验证拉普拉斯方程:
∂x2∂2ϕ+∂y2∂2ϕ=0+0=0✓
物理意义:均匀流的速度势是沿流向的线性函数。ϕ=const 的等势线是垂直于 x 轴的平面。
三、流函数——另一种描述流动的语言
3.1 概念引入
速度势要求流动无旋。但如果我们对可压缩性更敏感——要求流动不可压缩——那么我们可以引入另一个函数 ψ,称为流函数(Stream Function)。它自动满足连续方程(质量守恒),因此在不可压缩流动中始终可以定义。
3.2 定义
对于二维不可压缩流动,流函数 ψ 定义为:
u=∂y∂ψ,v=−∂x∂ψ
在极坐标下:
Vr=r1∂θ∂ψ,Vθ=−∂r∂ψ
自动满足连续方程:
∇⋅V=∂x∂u+∂y∂v=∂x∂(∂y∂ψ)+∂y∂(−∂x∂ψ)=∂x∂y∂2ψ−∂y∂x∂2ψ=0✓
3.3 流函数的意义
流函数 ψ 最有魅力的性质是:它的等值线就是流线。
证明:
沿着一条 ψ=const 的曲线:
dψ=∂x∂ψdx+∂y∂ψdy=−vdx+udy=0
整理:
dxdy=uv
这正是流线的微分方程。所以 ψ=const 就是流线。
另一个重要性质:两条流线之间的体积流量 Q˙ 等于流函数之差:
Q˙=ψ2−ψ1
这个性质非常实用——你只需要读出两条流线的 ψ 值,就能知道它们之间的流量。
3.4 手算例子:均匀流的流函数
对于均匀流 u=V∞,v=0:
u=∂y∂ψ=V∞⇒ψ=V∞y+g(x)
v=−∂x∂ψ=−g′(x)=0⇒g(x)=const
取常数项为零:
ψ=V∞y
验证:ψ=const 的等值线是 y=const 的水平直线——这正是均匀流的流线。
四、势流理论的”乐高积木”——四种基本流动解
势流理论最优雅的地方在于:拉普拉斯方程是线性的,所以我们可以把多个基本解叠加起来,构造出更复杂的流动!
Anderson 教材中介绍了四种基本流动解,它们是势流理论的”乐高积木”:
| 基本解 |
速度势 ϕ |
流函数 ψ |
| 均匀流 |
V∞x |
V∞y |
| 源/汇 |
2πΛlnr |
2πΛθ |
| 偶极子 |
2πμx2+y2x |
−2πμx2+y2y |
| 涡(点涡) |
2πΓθ |
−2πΓlnr |
其中:
- Λ 是源强(m2/s),Λ>0 为源,Λ<0 为汇
- μ 是偶极矩强度
- Γ 是环量(m2/s),Γ>0 为逆时针涡
下面我们逐一分析每一种基本解。
定义:
均匀流是最简单的基本流动——全流场速度恒定,方向一致。
速度势:ϕ=V∞x=V∞rcosθ
流函数:ψ=V∞y=V∞rsinθ
速度场:
u=V∞,v=0
物理意义:这代表了远离物体时”未受扰动”的自由来流。在风洞中,风机产生的就是近似均匀流。
4.2 源流与汇流(Source and Sink)
4.2.1 概念引入
想象在平静的水面上有一个不断冒水的孔——水从孔中均匀地向四面八方流出。这就是源流。如果孔在吸水(水从四面八方流入孔中),这就是汇流。
4.2.2 定义
源(Source):流体从中心点径向向外流出,流线是从中心点发出的射线。
速度势(极坐标):
ϕ=2πΛlnr
流函数:
ψ=2πΛθ
速度场:
Vr=∂r∂ϕ=2πrΛ,Vθ=r1∂θ∂ϕ=0
源强 Λ 的物理含义:单位长度(三维中是单位深度)上从源流出的体积流量。单位是 m2/s(二维)。
汇(Sink):Λ<0。流体径向流入中心点。
4.2.3 性质
- 速度随半径增大而减小:Vr∝1/r。离源越远,强度越弱。
- 在中心处(r=0),Vr→∞。这称为奇点——理论上的缺陷,但远离奇点时解非常有用。
- 流线是径向射线(θ=const),等势线是同心圆(r=const)。
4.2.4 手算例子
一个源强 Λ=2 m2/s 的源,求距中心 r=0.5 m 处的速度大小。
解:
Vr=2πrΛ=2π×0.52=π2≈0.637 m/s
这个速度大约相当于微风。如果离源 1 m 远,速度降为 0.318 m/s——确实符合 1/r 衰减规律。
4.3 偶极子(Doublet)
4.3.1 概念引入
偶极子是一个源和一个等强度的汇无限接近的极限情况。它产生的流场形状就像”∞”形:
- 流体从一边”吐出”,再从另一边”吸入”
- 流线是从一点出发、回到另一点的闭合曲线
4.3.2 定义
将强度为 Λ 的源放在 (−ϵ,0),强度为 −Λ 的汇放在 (ϵ,0)。让 ϵ→0 的同时保持 μ=2ϵΛ 为常数,得到偶极子。
偶极矩强度 μ(m3/s,二维时)。
速度势:
ϕ=2πμx2+y2x=2πμrcosθ
流函数:
ψ=−2πμx2+y2y=−2πμrsinθ
速度场(极坐标):
Vr=∂r∂ϕ=−2πμr2cosθ
Vθ=r1∂θ∂ϕ=−2πμr2sinθ
4.3.3 性质
- 偶极子的速度也按 1/r2 衰减——比源/汇(1/r)衰减更快。
- 在原点处的速度无穷大(另一个奇点)。
- 偶极子的方向:这里的 μ 沿 x 轴正方向(从汇指向源),如果 μ 取负值,则方向相反。
- 偶极子最重要的应用:均匀流 + 偶极子 = 圆柱绕流(将在超声速理论中进一步讨论)。
4.4 点涡(Vortex)
4.4.1 概念引入
不同于前面三种基本解(它们的速度场是辐射状或偶极状),点涡引入了一个旋转的流动——流体绕着中心点做圆周运动。
关键区别:
- 源/汇/偶极子:不可压缩,无旋
- 点涡:不可压缩,除中心点外无旋
是的——除了中心那个奇点,点涡流场中每一个微团的涡量都为零!这是一个”整体有旋、局部无旋”的微妙情况。
4.4.2 定义
速度势:
ϕ=2πΓθ
流函数:
ψ=−2πΓlnr
速度场:
Vr=0,Vθ=−r1∂θ∂ϕ=−2πrΓ
其中 Γ 是环量(Circulation),单位为 m2/s。
4.4.3 环量——最核心的概念
环量的定义:
Γ=∮CV⋅dl
即沿闭合曲线 C 对速度的线积分。单位是 m2/s。
对于点涡,取包围中心的任意闭合路径计算环量:
Γ=∮CVθ⋅rdθ=∫02π2πrΓ⋅rdθ=Γ
环量 = 常数!与路径的半径无关。这就是环量最关键的性质——它不随路径衰减。
4.4.4 性质
- 除奇点外无旋:在 r>0 的任何一点,∇×V=0。这很好验证——速度只有 θ 分量,且 Vθ=Γ/(2πr)。
- 速度按 1/r 衰减:远离涡中心速度逐渐减小。
- 环量守恒(开尔文定理):对于无黏正压流动,沿闭合流体线(跟随流体运动的闭合回路)的环量不随时间变化。
- 流线是同心圆(r=const),等势线是径向射线(θ=const)——与源/汇正好相反。
4.4.5 手算例子
一个环量 Γ=50 m2/s 的点涡,求 r=0.1 m 处的速度。
解:
Vθ=2πrΓ=2π×0.150=0.62850≈79.6 m/s
这是一个非常大的速度(约 286 km/h)。如果距离增大到 r=1 m,速度降为约 7.96 m/s。
这个例子说明:点涡的诱导速度在近处非常强、远处迅速衰减。这就是为什么翼尖涡对附近飞行的飞机会产生很大的影响,但随着距离增大很快减弱。
五、叠加原理——用”乐高积木”搭出真实流动
5.1 概念引入
拉普拉斯方程的线性特性意味着:如果 ϕ1 和 ϕ2 都是拉普拉斯方程的解,那么它们的和 ϕ1+ϕ2 也是解。
这就给了我们叠加的威力:用简单的基本解来构造复杂的流动!
5.2 定义
若 ∇2ϕ1=0,∇2ϕ2=0,则:
∇2(ϕ1+ϕ2)=∇2ϕ1+∇2ϕ2=0
对于速度场:
V=∇(ϕ1+ϕ2)=∇ϕ1+∇ϕ2=V1+V2
流场中任意一点的速度 = 各基本解在该点速度的矢量叠加。
5.3 手算例子:源 + 均匀流 = 半无限体
这是叠加法最简单的应用之一。
叠加:均匀流(V∞ 沿 x 正方向)+ 源强 Λ 的源(位于原点)
组合速度势:
ϕ=V∞x+2πΛlnr=V∞rcosθ+2πΛlnr
组合流函数:
ψ=V∞y+2πΛθ=V∞rsinθ+2πΛθ
速度场:
Vr=V∞cosθ+2πrΛ
Vθ=−V∞sinθ
物理意义:这个叠加描述了一个半无限体(Rankine half-body)外部流动——类似于一个鼻锥形物体放在均匀来流中。源将均匀流”推开”,形成一条 ψ=0 的流线,代表物体的表面。
驻点位置:
在流场中寻找 V=0 的点。令 Vr=0,Vθ=0:
从 Vθ=0 得 θ=0(x 正半轴方向)。
代入 Vr:
V∞+2πrΛ=0⇒r=−2πV∞Λ
由于 Λ>0,r 为正值意味着驻点在 θ=π(x 负半轴方向)。驻点位于源的上游。
这种叠加的魅力:用两个简单解的线性组合,我们就能模拟一个”航天器头锥”外形的绕流场!
六、完整概念地图
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21
| 黏性 → 极薄边界层(Re→∞) → 边界层外无旋 ↓ 无旋流动: ∇×V = 0 / \ ↓ ↓ ∇²φ = 0 ∇²ψ = 0 (拉普拉斯方程) (拉普拉斯方程) ↓ ↓ 基本解(乐高积木) 叠加原理 → 复杂流动 / | \ \ ↓ ↓ ↓ ↓ 均匀流 源/汇 偶极子 点涡 | | | | ↓ ↓ ↓ ↓ V∞x (Λ/2π)lnr μcosθ/(2πr) Γθ/(2π) | 叠加 + 均匀流 ↓ [势流理论的纵深] 均匀流+偶极子→圆柱绕流 圆柱绕流+点涡→产生升力! ← 库塔-儒科夫斯基定理
|
七、各概念在已发表文章中的出现
| 概念 |
出现文章 |
使用方式 |
| 速度势 |
— |
本文首次系统性引入 |
| 流函数 |
— |
本文首次系统性引入 |
| 源/汇 |
《可微物理引擎完全解析》 |
用于仿真代理模型中的诱导速度场参数化 |
| 环量 |
《前进比完全解读》 |
间接涉及——桨叶环量与推力关系 |
| 点涡 |
《四旋翼飞行力学基础》 |
翼尖涡与诱导阻力(简化讨论) |
| 叠加原理 |
— |
本文首次系统性引入 |
八、核心公式速查卡
| 公式 |
含义 |
使用场景 |
| ∇2ϕ=0 |
拉普拉斯方程 |
势流的控制方程 |
| ω=∇×V |
涡量定义 |
判断有无旋 |
| u=∂ϕ/∂x |
速度势提取速度 |
已知ϕ求u,v |
| u=∂ψ/∂y, v=−∂ψ/∂x |
流函数定义 |
已知ψ求u,v |
| Γ=∮V⋅dl |
环量定义 |
升力理论核心量 |
| ϕuniform=V∞x |
均匀流速度势 |
来流叠加 |
| ϕsource=(Λ/2π)lnr |
源流速度势 |
叠加构造物体外形 |
| ϕdoublet=μcosθ/(2πr) |
偶极子速度势 |
圆柱绕流 |
| ϕvortex=Γθ/(2π) |
点涡速度势 |
产生环量→升力 |
九、给自学者的建议
势流理论是空气动力学中数学密度最高的部分之一。自学时请注意以下几点:
- 先理解”为什么需要这个”:势流理论不是纯数学游戏——它的目标是构造机翼绕流场。始终记住:我们要的是机翼上的速度分布,进而用伯努利方程得到压力分布、计算升力。
- 做手算:拿一张纸,把四种基本解的速度势、流函数、速度场自己推一遍。拉普拉斯方程的解法其实简单——关键是熟悉极坐标下的算符。
- 关注叠加法的精髓:均匀流 + 偶极子 = 圆柱绕流,再加一个点涡 = 带环量的圆柱(即产生升力的圆柱),这是整个升力理论的数学核心。
- 记住环量的意义:环量 Γ 不是”旋涡强度”——它是绕物体积分出来的一个量,形式上定义简单,物理上承载了全部升力信息。
参考文献
- Anderson, J. D. (2010). Fundamentals of Aerodynamics (6th ed., Chapter 3). McGraw-Hill. ISBN: 978-0-07-339810-5.
- Batchelor, G. K. (2000). An Introduction to Fluid Dynamics (Chapters 2, 6). Cambridge University Press. ISBN: 978-0-521-66396-0.
- Lamb, H. (1932). Hydrodynamics (6th ed., Chapter 3). Cambridge University Press.
- Milne-Thomson, L. M. (1996). Theoretical Hydrodynamics (5th ed., Chapters 4-6). Dover Publications. ISBN: 978-0-486-68970-8.
- Prandtl, L., & Tietjens, O. G. (1934). Fundamentals of Hydro- and Aerodynamics (Chapters 1-4). Dover Publications.
下一节内容:环量与升力
我们将用均匀流 + 偶极子构造出圆柱绕流的精确解(无需准备——直接给出),然后精雕细琢最关键的部分:在圆柱绕流上加一个点涡→产生非对称压力分布→产生升力。这将引出空气动力学最重要的定理之一——库塔-儒科夫斯基定理:L′=ρV∞Γ。然后再讨论:为什么机翼周围存在环量?库塔条件如何指定了环量的大小?最后推广到保角变换和儒科夫斯基翼型。