跳到主要内容
Sandbox Physics

流体力学 · 翼型测压工作台

Airfoil Circulation Studio

保留可旋转的三维风洞，把上下表面压力、可移动测压点、表面压力积分与环流升力接到同一套 Joukowski 解。改变设计和采样精度，检查升力与零阻力的闭环误差；关闭 Kutta 条件，观察后缘奇性，而不是被截断后的“正常”压力。

交互模型Airfoil Circulation Studio
壁面／后缘残差（相对来流）analytic\text{analytic}
实际采样／面板数128128
净源通量残差（归一化）analytic\text{analytic}
环流 Γ\Gamma0 m2 s−10\,\mathrm{m^2\,s^{-1}}
理论升力 CLC_L00
压力积分升力 CL,pC_{L,p}00
力积分闭环绝对误差00
压力积分阻力 CD,pC_{D,p}00
后缘极限／后缘面板估计 VTE/U∞V_{\mathrm{TE}}/U_\infty00
测点上下压差 ΔCp\Delta C_p00
单位翼展升力 L′L'0 N m−10\,\mathrm{N\,m^{-1}}
模型边界二维势流

物理教程

翼型环流：后缘条件怎样选出升力

背景不可压缩、无黏、无旋流动可以先在圆柱外部求解，再通过保角映射搬到翼型周围。偏心圆的位置同时控制翼型弯度和厚度，而尖锐后缘对应映射的一个临界点。

为什么值得理解这个模型不会预测失速，却能把几何、环流、表面速度、压力差与升力连成一条完整因果链。

先抓重点

核心问题
同一个翼型外部为什么存在许多环流解，而真实附着流动只选择其中一个？
一句话直觉
Kutta 条件要求速度有限且流动平滑离开后缘，因此选出唯一的环流 Γ\Gamma；随后压力差给出升力。

核心数学模型

NACA 几何与真实源／涡面板

yt=5t(0.2969x−0.1260x−0.3516x2+0.2843x3−0.1036x4),V⋅n=0,Vt,1+Vt,N=0y_t=5t\left(0.2969\sqrt{x}-0.1260x-0.3516x^2+0.2843x^3-0.1036x^4\right),\quad \mathbf V\cdot\mathbf n=0,\quad V_{t,1}+V_{t,N}=0

使用 NASA TM 4741 的四位数中弧线和厚度多项式，把末项修改为闭合后缘。Hess–Smith 用各面板源强与一个均匀涡层强度联立无穿透和 Kutta 条件；表面压力、流线与环流都来自同一个数值解。面板残差、源通量误差及半分辨率力积分与解析 Joukowski 模式分开显示；小残差不能代替网格收敛。

有限后缘速度不是零速度

VTEU∞=bR∣cos⁡(α−β)∣,β=arg⁡(b−ζc)\frac{V_{\mathrm{TE}}}{U_\infty}=\frac{b}{R}\left|\cos(\alpha-\beta)\right|,\qquad \beta=\arg(b-\zeta_c)

在 Kutta 解中，复速度分子与映射导数有共同的后缘零点。解析消去后得到有限速度，通常并不为零。关闭 Kutta 后若共同零点不再存在，就没有有限极限；诊断会停用力积分，而不是截断奇性。

独立的表面压力力积分

Cx=−∮Cp d(y/c),Cy=∮Cp d(x/c),CL,p=−Cxsin⁡α+Cycos⁡αC_x=-\oint C_p\,\mathrm{d}(y/c),\quad C_y=\oint C_p\,\mathrm{d}(x/c),\quad C_{L,p}=-C_x\sin\alpha+C_y\cos\alpha

沿逆时针翼型轮廓积分未截断的压力，再把力投影到垂直来流的升力方向。诊断显示与环流公式的差异和数值阻力；半采样数误差是实际重算，不是预设精度估计。参考 MIT 18.354J 第 20 讲的压力轮廓积分。

Joukowski 保角映射

z=ζ+b2ζz=\zeta+\frac{b^2}{\zeta}

原像平面中的偏心圆被映成具有圆钝前缘和尖锐后缘的翼型截面。

圆柱外部的复势

W(ζ)=U∞(e−iαη+eiαR2η)+iΓ2πlog⁡η,η=ζ−ζcW(\zeta)=U_\infty\left(e^{-i\alpha}\eta+e^{i\alpha}\frac{R^2}{\eta}\right)+\frac{i\Gamma}{2\pi}\log\eta,\qquad \eta=\zeta-\zeta_c

均匀来流、圆柱镜像项和点涡环流共同构成圆柱平面中的解析解。

压力与升力

Cp=1−∣VU∞∣2,L′=ρU∞Γ,CD=0C_p=1-\left\lvert\frac{V}{U_\infty}\right\rvert^2,\qquad L'=\rho U_\infty\Gamma,\qquad C_D=0

Bernoulli 关系把表面速度变成压力；在后缘有界的解中，Kutta–Joukowski 给出单位翼展升力，闭合压力积分给出零阻力。该结果暴露了模型缺失黏性的边界；奇性解不适用这里的闭合力积分对照。

常见难点

Kutta 条件不是额外的推力

容易误解打开 Kutta 条件等于给翼型施加一个隐藏的升力源。

正确理解它不是外力，而是从许多数学允许的环流解中，挑出后缘速度有限且离开方向平滑的那一个。

零阻力不是现实结论

容易误解压力图显示升力却没有阻力，所以真实翼型也能无损飞行。

正确理解无黏势流忽略边界层、尾迹与表面摩擦，因此必然遗漏阻力与失速；这个矛盾就是 D’Alembert 悖论。

动手实验

  1. 01

    建立零升力基线

    选择“对称 · 零升力”，同时保留 Kutta 条件。

    应该观察到： 上下表面流线和压力近似对称，环流与升力读数接近零。
  2. 02

    只增加攻角

    选择“巡航翼型”，观察上表面的低压测点与升力矢量。

    应该观察到： 上表面流动加速、压力系数降低，选出的环流和单位翼展升力同步增大。
  3. 03

    关闭后缘条件

    选择“关闭 Kutta”，比较几何不变时的流线和压力。

    应该观察到： 零环流使理论升力为零，但后缘出现速度奇性；压力曲线可能超出画幅，原始数据仍保留。压力积分显示不适用，不能把截断后的有限数值当作正常解。
  4. 04

    用两条计算路径检查精度

    回到巡航翼型，把积分采样从 N=32N=32 加密至 N=128N=128；再移动测压位置。

    应该观察到： 压力积分升力趋近环流公式，数值阻力趋近零。解析 Joukowski 模式直接评估翼面，面板模式对计算出的面板中点压强插值，因此局部测压值可能随加密改变。

概念灵感：FluidMotion。模拟与图形均为原创，模型适用边界见上文。