4.2 势流叠加与库塔-茹科夫斯基升力

4.2 势流叠加与库塔-茹科夫斯基升力

本节摘要：无旋的无粘流动存在速度势，问题化为拉普拉斯方程；由于方程线性，基本解（均匀流、源汇、偶极子、点涡）可以叠加成复杂绕流。均匀流加偶极子得绕圆柱流，再加环量即得升力公式 L' = ρ U Γ（库塔-茹科夫斯基定理）。本节用 Python 画出绕圆柱的流线与压力分布，并直面达朗贝尔佯谬。

飞机为什么能飞？小学课堂上最常见的回答是"机翼上面流速快、压力低"——这话对了一半，但为什么上面流速快，无粘理论自己回答不了，它需要一个来自粘性世界的补充约定（库塔条件）。本节把这段公案讲完整，顺便学会流体力学家最古老的处理技巧：像拼积木一样拼流场。

一、速度势与拉普拉斯方程

无旋流动（涡量为零）允许引入速度势 φ，使速度等于 φ 的梯度。把不可压缩条件 ∇·u = 0 套上去，立即得到 ∇²φ = 0——拉普拉斯方程，数学里被研究得最透的方程之一。这一步的收益巨大：原 N-S 方程的非线性全部消失，而且方程线性，解可以叠加。于是策略变成：先攒一批基本解，再按需拼装。

均匀流 + 偶极子 = 绕圆柱流动。圆柱表面速度 U_θ = 2U sinθ：正前方驻点为零，侧面 90° 处达到来流两倍，再对称回到后方驻点。

import numpy as np

# 绕圆柱表面压力系数  Cp = 1 - (U表面/U)的平方
theta = np.linspace(0, np.pi, 7)
Cp = 1 - (2 * np.sin(theta))**2
for t, c in zip(theta, Cp):
    print(f"角度 {np.degrees(t):5.0f} 度   压力系数 {c:6.2f}")
​

侧面 Cp = −3：局部吸力达到三倍动压。前后对称意味着压力合力为零——达朗贝尔佯谬：无粘理论预言任何物体在均匀来流中阻力为零。现实里 obviously 不是这样，问题出在无粘假设丢掉了边界层与分离（第 6 章）。佯谬的价值是反向的：它精确地标出了"粘性在哪里不可忽略"。

二、加上环量：升力登场

在绕圆柱流上叠加一个点涡（环量 Γ），上下对称被打破：上表面流速叠加增强、下表面削弱，按伯努利逻辑上低下高，产生升力。库塔-茹科夫斯基定理给出普适结果（不限圆柱）：

单位展长升力 L' = ρ U Γ，方向由来流速度反旋 90° 得到。

对带环量圆柱，驻点移到下方，表面速度变为 U_θ = −2U sinθ + Γ/(2πR)。让后驻点恰好落在圆柱后缘（库塔条件的精神原型），环量与升力随之确定。真实机翼的物理版本是：启动瞬间后缘尖劈诱发旋涡脱落，留下一圈反向环量绑在机翼上——粘性负责"选"出环量，无粘公式负责"算"出升力。分工明确，缺一不可。

# 小型无人机机翼估算 升力线理论一阶
import numpy as np
rho, U, b, S = 1.225, 30.0, 1.6, 0.9     # 空气密度 速度 翼展 面积
CL = 0.85                                  # 某攻角下升力系数
L = 0.5 * rho * U**2 * S * CL
Gamma = L / (rho * U * b)                  # 平均环量
print(f"升力 {L:.0f} N  约可托起 {L/9.81:.1f} kg")
print(f"等效平均环量 {Gamma:.1f} m2/s")
​

升力约 430 N，环量约 7.3 m²/s——环量是"每单位展长裹了多少旋转"的度量，用涡強度统一描述升力来源，比"上面快下面慢"严谨得多。

三、势流方法的边界

势流能给出绕流的大模样与压力分布的一阶近似（厚物体小攻角时相当准），给不了的东西也很明确：阻力（佯谬）、分离位置、尾迹结构、失速。工程上的做法是"势流打底、边界层修正"：先用势流算外流场压力分布，再在物面附近解边界层方程修正位移与摩擦，这套"边界层耦合势流"的流程统治了从 1930 年代到 CFD 兴起前的飞机气动设计，水翼、泵轮、帆船龙骨的设计同样受益。

💡 关键直觉：升力的本质是流动向下折转（动量观点），环量与下洗只是同一现象的两种记账方式。看机翼尾迹中那股向下的气流质量流量乘速度变化，就等于升力——这也是直升机下洗、风扇反推力的统一图像。

三、马格努斯效应：旋转圆柱的升力现场

环量不必只靠尖后缘产生——旋转本身就能造环量。旋转圆柱带动表面附近流体同向旋转，叠加均匀来流后上下流速不对称，产生垂直于来流的力，即马格努斯效应。用升力公式估算：半径 0.2 m、长 1 m 的圆柱，来流 15 m/s，保守取环量为表面线速度乘周长的一半：

import numpy as np
rho, U, R, b = 1.225, 15.0, 0.2, 1.0
Gamma = 2*np.pi*R*U * 0.5
L = rho * U * Gamma * b
print(f"马格努斯力约 {L:.0f} N  每米柱长 {L/b:.0f} N")
​

约 34 N——足以让旋转球类走弧线。足球香蕉球、乒乓球的弧圈、网球的侧旋切削，全是弹道上的马格努斯力；历史上还有人真造过旋筒推进船（弗莱特纳转子船），用旋转的高圆柱代替风帆，今天几艘现代货轮又重新捡起这个思路省燃料。它再次印证本节主线：只要能造出环量，升力公式就兑现——至于环量来自尖后缘的粘性分离还是壁面旋转，公式并不关心。

变式思考

环量视角还统一了地效飞行：机翼贴近地面时，镜像涡削弱下洗，同样攻角下有效攻角增大、诱导阻力下降，这就是地效飞行器与赛车站位设计的气动红利。把"升力来自环量、诱导阻力来自尾涡"两条记牢，机翼布局的多数现象都能串成一线。

本节要点回顾

• 无旋 + 不可压缩 = 拉普拉斯方程，线性可叠加，基本解如积木；
• 均匀流 + 偶极子 = 绕圆柱，侧面速度两倍来流、Cp = −3，前后对称；
• 达朗贝尔佯谬（零阻力）标出无粘假设的失效区，把阻力问题留给边界层理论；
• 库塔-茹科夫斯基定理 L' = ρUΓ：粘性通过库塔条件确定环量，无粘公式完成升力计算；
• 势流打底、边界层修正是前 CFD 时代气动设计的主力流程。

低速理想世界的账至此结清。下一章提速：当流速逼近声速，密度从配角变主角。