本文是 John D. Anderson《Fundamentals of Aerodynamics》(第六版)读书笔记系列的第二篇,覆盖教材第3章《势流理论》的前半部分。阅读本文前,建议先了解流体力学基础中的连续方程和伯努利方程的基本概念。


第二部分:势流理论——理解升力的数学语言

在第一篇中我们提到:如果忽略黏性,流体的动量方程简化为欧拉方程;进一步,如果流动无旋(irrotational),我们就可以引入一个叫做速度势的函数 ϕ\phi,使得流场的所有信息浓缩在一个标量拉普拉斯方程中。这就是势流理论的起点。

势流理论在 Anderson 的书中占据了第3-4章,是全书数学最密集的部分之一,也是理解升力产生的数学本质的必经之路。


一、涡量与无旋流动——为什么能”简化”?

1.1 概念引入

想象一个河面上的小水车。如果水车在流动中旋转,就说明这个地方的水有”旋涡”;如果水车平移而不旋转,就说明这里的水是”无旋”的。

这个区别非常重要。有旋流动和无旋流动,在数学上完全是两个世界

  • 有旋流动:需要用完整的 Navier-Stokes 方程,四个方程解四个未知数,数学上极其复杂
  • 无旋流动:可以引入速度势,化简为一个标量拉普拉斯方程,数学上大幅简化

而幸运的是,许多工程上重要的流动——包括机翼周围的附着流动——在很大的程度上可以视为无旋流动。

1.2 定义

一个流动是否”有旋”,由**涡量 ω\boldsymbol{\omega}(Vorticity)**来度量:

ω=×V=ijkxyzuvw \boldsymbol{\omega} = \nabla \times \mathbf{V} = \begin{vmatrix} \mathbf{i} & \mathbf{j} & \mathbf{k} \\ \frac{\partial}{\partial x} & \frac{\partial}{\partial y} & \frac{\partial}{\partial z} \\ u & v & w \end{vmatrix}

三维展开:

ω=(wyvz)i+(uzwx)j+(vxuy)k \boldsymbol{\omega} = \left(\frac{\partial w}{\partial y} - \frac{\partial v}{\partial z}\right)\mathbf{i} + \left(\frac{\partial u}{\partial z} - \frac{\partial w}{\partial x}\right)\mathbf{j} + \left(\frac{\partial v}{\partial x} - \frac{\partial u}{\partial y}\right)\mathbf{k}

无旋流动的定义:ω=0\boxed{\boldsymbol{\omega} = \mathbf{0}}(流场每一点的涡量为零)。

在二维流动中,这个条件简化为:

vxuy=0 \frac{\partial v}{\partial x} - \frac{\partial u}{\partial y} = 0

1.3 性质

  1. 无旋 ≠ 无涡:无旋流动并不意味着流线是直的。绕圆柱的流动可以是无旋的,但流线是弯曲的——“无旋”指的是流体微团不绕自身旋转,而不是流线不弯曲。
  2. 无旋具有”记忆”:如果均匀来流是无旋的,它流过光滑物体后仍保持无旋——这个性质称为开尔文环量定理(Kelvin’s circulation theorem)。
  3. 边界层内部是有旋的,外部是无旋的:黏性只在边界层内产生旋涡,边界层以外的势流区域可以视为无旋流动。这就是为什么势流理论仍然实用的根本原因。

1.4 历史

时间 人物 贡献
1781 拉普拉斯(Pierre-Simon Laplace) 建立了拉普拉斯方程 2ϕ=0\nabla^2\phi=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)\mathbf{V}(x,y,z) 是一个矢量函数——有三个分量 (u,v,w)(u,v,w),需要三个方程才能解。但如果我们能找到一个标量函数 ϕ(x,y,z)\phi(x,y,z),使得 V=ϕ\mathbf{V} = \nabla \phi,那么我们就可以用一个标量函数代替三个分量函数。

这是空气动力学中最重要的”降维”技巧。

2.2 定义

速度势(Velocity Potential)ϕ\phi 定义为满足以下关系的函数:

V=ϕ \mathbf{V} = \nabla \phi

在笛卡尔坐标系下展开:

u=ϕx,v=ϕy,w=ϕz u = \frac{\partial \phi}{\partial x}, \quad v = \frac{\partial \phi}{\partial y}, \quad w = \frac{\partial \phi}{\partial z}

在极坐标系下(二维):

Vr=ϕr,Vθ=1rϕθ V_r = \frac{\partial \phi}{\partial r}, \quad V_\theta = \frac{1}{r}\frac{\partial \phi}{\partial \theta}

存在性定理:流动无旋(×V=0\nabla \times \mathbf{V} = 0\Longleftrightarrow 存在速度势 ϕ\phi 使得 V=ϕ\mathbf{V} = \nabla \phi

这个定理来自向量微积分的基本定理:旋度为零的向量场一定是某个标量函数的梯度。

2.3 从速度势到拉普拉斯方程

将速度势代入连续方程(不可压缩):

V=0 \nabla \cdot \mathbf{V} = 0 (ϕ)=0 \nabla \cdot (\nabla \phi) = 0

这就得到了势流理论中最核心的方程——拉普拉斯方程

2ϕ=0 \boxed{\nabla^2 \phi = 0}

二维展开:

2ϕx2+2ϕy2=0 \frac{\partial^2 \phi}{\partial x^2} + \frac{\partial^2 \phi}{\partial y^2} = 0

这是一个二阶级线性偏微分方程。它的重要意义在于:

  1. 线性:解可以叠加——复杂流动 = 基本流动解的线性组合
  2. 成熟的理论:拉普拉斯方程的数学理论在 19 世纪就已经相当完善
  3. 一个约束明确边值问题:只需要在物体表面施加一个边界条件(ϕ/n=0\partial\phi/\partial n = 0,即壁面无穿透条件),就可以唯一确定流场

2.4 性质

速度势的物理含义:

  1. ϕ=const\phi = \text{const} 的等值面称为**等势面**。速度方向总是垂直于等势面。
  2. ϕ\phi 的大小本身没有物理意义,有意义的是**它的梯度**(速度)。
  3. 在不可压缩无旋流动中,ϕ\phi 满足拉普拉斯方程——这是整个势流理论的数学基础。

对比速度势与流函数(流函数将在下节介绍):

特性 速度势 ϕ\phi 流函数 ψ\psi
定义 V=ϕ\mathbf{V} = \nabla \phi u=ψy,v=ψxu = \frac{\partial\psi}{\partial y}, v = -\frac{\partial\psi}{\partial x}
存在条件 无旋(×V=0\nabla \times \mathbf{V}=0 不可压缩(V=0\nabla \cdot \mathbf{V}=0
基本方程 2ϕ=0\nabla^2 \phi = 0(拉普拉斯方程) 2ψ=0\nabla^2 \psi = 0(拉普拉斯方程)
等值线意义 等势线 流线
物理意义 速度的”势” 流量

2.5 手算例子:一维均匀流的速度势

最简单的势流——均匀流 u=Vu = V_\inftyv=0v=0

验证无旋

vxuy=00=0无旋条件满足 \frac{\partial v}{\partial x} - \frac{\partial u}{\partial y} = 0 - 0 = 0 \quad \Rightarrow \quad \text{无旋条件满足}

求速度势

u=ϕx=Vϕ=Vx+f(y) u = \frac{\partial \phi}{\partial x} = V_\infty \quad \Rightarrow \quad \phi = V_\infty x + f(y) v=ϕy=f(y)=0f(y)=const v = \frac{\partial \phi}{\partial y} = f'(y) = 0 \quad \Rightarrow \quad f(y) = \text{const}

取常数项为零,则:

ϕ=Vx \boxed{\phi = V_\infty x}

验证拉普拉斯方程

2ϕx2+2ϕy2=0+0=0 \frac{\partial^2 \phi}{\partial x^2} + \frac{\partial^2 \phi}{\partial y^2} = 0 + 0 = 0 \quad \checkmark

物理意义:均匀流的速度势是沿流向的线性函数。ϕ=const\phi = \text{const} 的等势线是垂直于 xx 轴的平面。


三、流函数——另一种描述流动的语言

3.1 概念引入

速度势要求流动无旋。但如果我们对可压缩性更敏感——要求流动不可压缩——那么我们可以引入另一个函数 ψ\psi,称为流函数(Stream Function)。它自动满足连续方程(质量守恒),因此在不可压缩流动中始终可以定义。

3.2 定义

对于二维不可压缩流动,流函数 ψ\psi 定义为:

u=ψy,v=ψx u = \frac{\partial \psi}{\partial y}, \quad v = -\frac{\partial \psi}{\partial x}

在极坐标下:

Vr=1rψθ,Vθ=ψr V_r = \frac{1}{r}\frac{\partial \psi}{\partial \theta}, \quad V_\theta = -\frac{\partial \psi}{\partial r}

自动满足连续方程

V=ux+vy=x(ψy)+y(ψx)=2ψxy2ψyx=0 \nabla \cdot \mathbf{V} = \frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} = \frac{\partial}{\partial x}\left(\frac{\partial \psi}{\partial y}\right) + \frac{\partial}{\partial y}\left(-\frac{\partial \psi}{\partial x}\right) = \frac{\partial^2\psi}{\partial x\partial y} - \frac{\partial^2\psi}{\partial y\partial x} = 0 \quad \checkmark

3.3 流函数的意义

流函数 ψ\psi 最有魅力的性质是:它的等值线就是流线

证明:

沿着一条 ψ=const\psi = \text{const} 的曲线:

dψ=ψxdx+ψydy=vdx+udy=0 d\psi = \frac{\partial \psi}{\partial x} dx + \frac{\partial \psi}{\partial y} dy = -v\,dx + u\,dy = 0

整理:

dydx=vu \frac{dy}{dx} = \frac{v}{u}

这正是流线的微分方程。所以 ψ=const\psi = \text{const} 就是流线。

另一个重要性质:两条流线之间的体积流量 Q˙\dot{Q} 等于流函数之差

Q˙=ψ2ψ1 \dot{Q} = \psi_2 - \psi_1

这个性质非常实用——你只需要读出两条流线的 ψ\psi 值,就能知道它们之间的流量。

3.4 手算例子:均匀流的流函数

对于均匀流 u=Vu = V_\inftyv=0v=0

u=ψy=Vψ=Vy+g(x) u = \frac{\partial \psi}{\partial y} = V_\infty \quad \Rightarrow \quad \psi = V_\infty y + g(x) v=ψx=g(x)=0g(x)=const v = -\frac{\partial \psi}{\partial x} = -g'(x) = 0 \quad \Rightarrow \quad g(x) = \text{const}

取常数项为零:

ψ=Vy \boxed{\psi = V_\infty y}

验证:ψ=const\psi = \text{const} 的等值线是 y=consty = \text{const} 的水平直线——这正是均匀流的流线。


四、势流理论的”乐高积木”——四种基本流动解

势流理论最优雅的地方在于:拉普拉斯方程是线性的,所以我们可以把多个基本解叠加起来,构造出更复杂的流动!

Anderson 教材中介绍了四种基本流动解,它们是势流理论的”乐高积木”:

基本解 速度势 ϕ\phi 流函数 ψ\psi
均匀流 VxV_\infty x VyV_\infty y
源/汇 Λ2πlnr\frac{\Lambda}{2\pi}\ln r Λ2πθ\frac{\Lambda}{2\pi}\theta
偶极子 μ2πxx2+y2\frac{\mu}{2\pi}\frac{x}{x^2+y^2} μ2πyx2+y2-\frac{\mu}{2\pi}\frac{y}{x^2+y^2}
涡(点涡) Γ2πθ\frac{\Gamma}{2\pi}\theta Γ2πlnr-\frac{\Gamma}{2\pi}\ln r

其中:

  • Λ\Lambda 是源强(m2/s\text{m}^2/\text{s}),Λ>0\Lambda>0 为源,Λ<0\Lambda<0 为汇
  • μ\mu 是偶极矩强度
  • Γ\Gamma 是环量(m2/s\text{m}^2/\text{s}),Γ>0\Gamma>0 为逆时针涡

下面我们逐一分析每一种基本解。

4.1 均匀流(Uniform Flow)

定义

均匀流是最简单的基本流动——全流场速度恒定,方向一致。

速度势ϕ=Vx=Vrcosθ\phi = V_\infty x = V_\infty r\cos\theta

流函数ψ=Vy=Vrsinθ\psi = V_\infty y = V_\infty r\sin\theta

速度场

u=V,v=0 u = V_\infty, \quad v = 0

物理意义:这代表了远离物体时”未受扰动”的自由来流。在风洞中,风机产生的就是近似均匀流。

4.2 源流与汇流(Source and Sink)

4.2.1 概念引入

想象在平静的水面上有一个不断冒水的孔——水从孔中均匀地向四面八方流出。这就是源流。如果孔在吸水(水从四面八方流入孔中),这就是汇流

4.2.2 定义

源(Source):流体从中心点径向向外流出,流线是从中心点发出的射线。

速度势(极坐标):

ϕ=Λ2πlnr \phi = \frac{\Lambda}{2\pi}\ln r

流函数

ψ=Λ2πθ \psi = \frac{\Lambda}{2\pi}\theta

速度场

Vr=ϕr=Λ2πr,Vθ=1rϕθ=0 V_r = \frac{\partial \phi}{\partial r} = \frac{\Lambda}{2\pi r}, \quad V_\theta = \frac{1}{r}\frac{\partial \phi}{\partial \theta} = 0

源强 Λ\Lambda 的物理含义:单位长度(三维中是单位深度)上从源流出的体积流量。单位是 m2/s\text{m}^2/\text{s}(二维)。

汇(Sink)Λ<0\Lambda < 0。流体径向流入中心点。

4.2.3 性质

  1. 速度随半径增大而减小:Vr1/rV_r \propto 1/r。离源越远,强度越弱。
  2. 在中心处(r=0r=0),VrV_r \to \infty。这称为奇点——理论上的缺陷,但远离奇点时解非常有用。
  3. 流线是径向射线(θ=const\theta = \text{const}),等势线是同心圆(r=constr = \text{const})。

4.2.4 手算例子

一个源强 Λ=2 m2/s\Lambda = 2\ \text{m}^2/\text{s} 的源,求距中心 r=0.5 mr=0.5\ \text{m} 处的速度大小。

Vr=Λ2πr=22π×0.5=2π0.637 m/s V_r = \frac{\Lambda}{2\pi r} = \frac{2}{2\pi\times 0.5} = \frac{2}{\pi} \approx 0.637\ \text{m/s}

这个速度大约相当于微风。如果离源 1 m 远,速度降为 0.318 m/s——确实符合 1/r1/r 衰减规律。

4.3 偶极子(Doublet)

4.3.1 概念引入

偶极子是一个源和一个等强度的汇无限接近的极限情况。它产生的流场形状就像”∞”形:

  • 流体从一边”吐出”,再从另一边”吸入”
  • 流线是从一点出发、回到另一点的闭合曲线

4.3.2 定义

将强度为 Λ\Lambda 的源放在 (ϵ,0)(-\epsilon,0),强度为 Λ-\Lambda 的汇放在 (ϵ,0)(\epsilon,0)。让 ϵ0\epsilon \to 0 的同时保持 μ=2ϵΛ\mu = 2\epsilon\Lambda 为常数,得到偶极子。

偶极矩强度 μ\mum3/s\text{m}^3/\text{s},二维时)。

速度势

ϕ=μ2πxx2+y2=μ2πcosθr \phi = \frac{\mu}{2\pi}\frac{x}{x^2+y^2} = \frac{\mu}{2\pi}\frac{\cos\theta}{r}

流函数

ψ=μ2πyx2+y2=μ2πsinθr \psi = -\frac{\mu}{2\pi}\frac{y}{x^2+y^2} = -\frac{\mu}{2\pi}\frac{\sin\theta}{r}

速度场(极坐标):

Vr=ϕr=μ2πcosθr2 V_r = \frac{\partial \phi}{\partial r} = -\frac{\mu}{2\pi}\frac{\cos\theta}{r^2} Vθ=1rϕθ=μ2πsinθr2 V_\theta = \frac{1}{r}\frac{\partial \phi}{\partial \theta} = -\frac{\mu}{2\pi}\frac{\sin\theta}{r^2}

4.3.3 性质

  1. 偶极子的速度也按 1/r21/r^2 衰减——比源/汇(1/r1/r)衰减更快。
  2. 在原点处的速度无穷大(另一个奇点)。
  3. 偶极子的方向:这里的 μ\mu 沿 xx 轴正方向(从汇指向源),如果 μ\mu 取负值,则方向相反。
  4. 偶极子最重要的应用:均匀流 + 偶极子 = 圆柱绕流(将在超声速理论中进一步讨论)。

4.4 点涡(Vortex)

4.4.1 概念引入

不同于前面三种基本解(它们的速度场是辐射状或偶极状),点涡引入了一个旋转的流动——流体绕着中心点做圆周运动。

关键区别:

  • 源/汇/偶极子:不可压缩,无旋
  • 点涡:不可压缩,除中心点外无旋

是的——除了中心那个奇点,点涡流场中每一个微团的涡量都为零!这是一个”整体有旋、局部无旋”的微妙情况。

4.4.2 定义

速度势

ϕ=Γ2πθ \phi = \frac{\Gamma}{2\pi}\theta

流函数

ψ=Γ2πlnr \psi = -\frac{\Gamma}{2\pi}\ln r

速度场

Vr=0,Vθ=1rϕθ=Γ2πr V_r = 0, \quad V_\theta = -\frac{1}{r}\frac{\partial\phi}{\partial\theta} = -\frac{\Gamma}{2\pi r}

其中 Γ\Gamma环量(Circulation),单位为 m2/s\text{m}^2/\text{s}

4.4.3 环量——最核心的概念

环量的定义:

Γ=CVdl \boxed{\Gamma = \oint_C \mathbf{V} \cdot d\mathbf{l}}

即沿闭合曲线 CC 对速度的线积分。单位是 m2/s\text{m}^2/\text{s}

对于点涡,取包围中心的任意闭合路径计算环量:

Γ=CVθrdθ=02πΓ2πrrdθ=Γ \Gamma = \oint_C V_\theta \cdot r\,d\theta = \int_0^{2\pi} \frac{\Gamma}{2\pi r} \cdot r\,d\theta = \Gamma

环量 = 常数!与路径的半径无关。这就是环量最关键的性质——它不随路径衰减。

4.4.4 性质

  1. 除奇点外无旋:在 r>0r>0 的任何一点,×V=0\nabla \times \mathbf{V} = 0。这很好验证——速度只有 θ\theta 分量,且 Vθ=Γ/(2πr)V_\theta = \Gamma/(2\pi r)
  2. 速度按 1/r1/r 衰减:远离涡中心速度逐渐减小。
  3. 环量守恒(开尔文定理):对于无黏正压流动,沿闭合流体线(跟随流体运动的闭合回路)的环量不随时间变化。
  4. 流线是同心圆r=constr = \text{const}),等势线是径向射线(θ=const\theta = \text{const})——与源/汇正好相反。

4.4.5 手算例子

一个环量 Γ=50 m2/s\Gamma = 50\ \text{m}^2/\text{s} 的点涡,求 r=0.1 mr=0.1\ \text{m} 处的速度。

Vθ=Γ2πr=502π×0.1=500.62879.6 m/s V_\theta = \frac{\Gamma}{2\pi r} = \frac{50}{2\pi \times 0.1} = \frac{50}{0.628} \approx 79.6\ \text{m/s}

这是一个非常大的速度(约 286 km/h)。如果距离增大到 r=1 mr=1\ \text{m},速度降为约 7.96 m/s。

这个例子说明:点涡的诱导速度在近处非常强、远处迅速衰减。这就是为什么翼尖涡对附近飞行的飞机会产生很大的影响,但随着距离增大很快减弱。


五、叠加原理——用”乐高积木”搭出真实流动

5.1 概念引入

拉普拉斯方程的线性特性意味着:如果 ϕ1\phi_1ϕ2\phi_2 都是拉普拉斯方程的解,那么它们的 ϕ1+ϕ2\phi_1 + \phi_2 也是解。

这就给了我们叠加的威力:用简单的基本解来构造复杂的流动!

5.2 定义

2ϕ1=0\nabla^2 \phi_1 = 02ϕ2=0\nabla^2 \phi_2 = 0,则:

2(ϕ1+ϕ2)=2ϕ1+2ϕ2=0 \nabla^2 (\phi_1 + \phi_2) = \nabla^2 \phi_1 + \nabla^2 \phi_2 = 0

对于速度场:

V=(ϕ1+ϕ2)=ϕ1+ϕ2=V1+V2 \mathbf{V} = \nabla (\phi_1 + \phi_2) = \nabla \phi_1 + \nabla \phi_2 = \mathbf{V}_1 + \mathbf{V}_2

流场中任意一点的速度 = 各基本解在该点速度的矢量叠加

5.3 手算例子:源 + 均匀流 = 半无限体

这是叠加法最简单的应用之一。

叠加:均匀流(VV_\infty 沿 xx 正方向)+ 源强 Λ\Lambda 的源(位于原点)

组合速度势

ϕ=Vx+Λ2πlnr=Vrcosθ+Λ2πlnr \phi = V_\infty x + \frac{\Lambda}{2\pi}\ln r = V_\infty r\cos\theta + \frac{\Lambda}{2\pi}\ln r

组合流函数

ψ=Vy+Λ2πθ=Vrsinθ+Λ2πθ \psi = V_\infty y + \frac{\Lambda}{2\pi}\theta = V_\infty r\sin\theta + \frac{\Lambda}{2\pi}\theta

速度场

Vr=Vcosθ+Λ2πr V_r = V_\infty\cos\theta + \frac{\Lambda}{2\pi r} Vθ=Vsinθ V_\theta = -V_\infty\sin\theta

物理意义:这个叠加描述了一个半无限体(Rankine half-body)外部流动——类似于一个鼻锥形物体放在均匀来流中。源将均匀流”推开”,形成一条 ψ=0\psi=0 的流线,代表物体的表面。

驻点位置

在流场中寻找 V=0\mathbf{V} = 0 的点。令 Vr=0V_r = 0Vθ=0V_\theta = 0

Vθ=0V_\theta = 0θ=0\theta = 0xx 正半轴方向)。

代入 VrV_r

V+Λ2πr=0r=Λ2πV V_\infty + \frac{\Lambda}{2\pi r} = 0 \quad \Rightarrow \quad r = -\frac{\Lambda}{2\pi V_\infty}

由于 Λ>0\Lambda > 0rr 为正值意味着驻点在 θ=π\theta = \pixx 负半轴方向)。驻点位于源的上游。

这种叠加的魅力:用两个简单解的线性组合,我们就能模拟一个”航天器头锥”外形的绕流场!


六、完整概念地图

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\nabla^2 \phi = 0 拉普拉斯方程 势流的控制方程
ω=×V\boldsymbol{\omega}=\nabla\times\mathbf{V} 涡量定义 判断有无旋
u=ϕ/xu=\partial\phi/\partial x 速度势提取速度 已知ϕ\phiu,vu,v
u=ψ/yu=\partial\psi/\partial y, v=ψ/xv=-\partial\psi/\partial x 流函数定义 已知ψ\psiu,vu,v
Γ=Vdl\Gamma = \oint \mathbf{V}\cdot d\mathbf{l} 环量定义 升力理论核心量
ϕuniform=Vx\phi_{\text{uniform}} = V_\infty x 均匀流速度势 来流叠加
ϕsource=(Λ/2π)lnr\phi_{\text{source}} = (\Lambda/2\pi)\ln r 源流速度势 叠加构造物体外形
ϕdoublet=μcosθ/(2πr)\phi_{\text{doublet}} = \mu\cos\theta/(2\pi r) 偶极子速度势 圆柱绕流
ϕvortex=Γθ/(2π)\phi_{\text{vortex}} = \Gamma\theta/(2\pi) 点涡速度势 产生环量→升力

九、给自学者的建议

势流理论是空气动力学中数学密度最高的部分之一。自学时请注意以下几点:

  1. 先理解”为什么需要这个”:势流理论不是纯数学游戏——它的目标是构造机翼绕流场。始终记住:我们要的是机翼上的速度分布,进而用伯努利方程得到压力分布、计算升力。
  2. 做手算:拿一张纸,把四种基本解的速度势、流函数、速度场自己推一遍。拉普拉斯方程的解法其实简单——关键是熟悉极坐标下的算符。
  3. 关注叠加法的精髓:均匀流 + 偶极子 = 圆柱绕流,再加一个点涡 = 带环量的圆柱(即产生升力的圆柱),这是整个升力理论的数学核心。
  4. 记住环量的意义:环量 Γ\Gamma 不是”旋涡强度”——它是绕物体积分出来的一个量,形式上定义简单,物理上承载了全部升力信息。

参考文献

  1. Anderson, J. D. (2010). Fundamentals of Aerodynamics (6th ed., Chapter 3). McGraw-Hill. ISBN: 978-0-07-339810-5.
  2. Batchelor, G. K. (2000). An Introduction to Fluid Dynamics (Chapters 2, 6). Cambridge University Press. ISBN: 978-0-521-66396-0.
  3. Lamb, H. (1932). Hydrodynamics (6th ed., Chapter 3). Cambridge University Press.
  4. Milne-Thomson, L. M. (1996). Theoretical Hydrodynamics (5th ed., Chapters 4-6). Dover Publications. ISBN: 978-0-486-68970-8.
  5. Prandtl, L., & Tietjens, O. G. (1934). Fundamentals of Hydro- and Aerodynamics (Chapters 1-4). Dover Publications.

下一节内容:环量与升力

我们将用均匀流 + 偶极子构造出圆柱绕流的精确解(无需准备——直接给出),然后精雕细琢最关键的部分:在圆柱绕流上加一个点涡→产生非对称压力分布→产生升力。这将引出空气动力学最重要的定理之一——库塔-儒科夫斯基定理:L=ρVΓL' = \rho V_\infty \Gamma。然后再讨论:为什么机翼周围存在环量?库塔条件如何指定了环量的大小?最后推广到保角变换和儒科夫斯基翼型。