跳到主要内容
Sandbox Physics

非定常空气动力学 · 扑翼风洞

Flapping Wing

一片三维翼段在透明风洞中做正弦升沉和俯仰。改变频率、振幅、相位和转轴位置，观察环量升力、附加质量、交替脱落的尾迹涡，以及整周期推力与功率代理怎样同步改变。

交互模型Flapping Wing
面板壁面残差00
已完成平均周期00
运动相位0∘0^\circ
约化频率 kk0.180.18
Strouhal 数 StSt0.080.08
瞬时升力 CLC_L0.0000.000
平均推力代理 CT∗‾\overline{C_T^{\ast}}0.0000.000
平均功率代理 CP∗‾\overline{C_P^{\ast}}0.0000.000
推进效率代理 η∗\eta^{\ast}0%0\%
Kelvin 环量账本漂移 ϵΓ\epsilon_\Gamma0.0×10−80.0\times10^{-8}

物理教程

扑翼风洞：相位怎样把振荡变成推进

背景默认模式联立求解薄翼面板的无穿透条件与 Kelvin 环量守恒，并让脱落涡反过来作用于翼段。可切换的 Jones–Wagner 谐波参考来自 Theodorsen 的 NACA 496 框架；两种方法不是同一个模型。

为什么值得理解升沉或俯仰本身只是往复运动；真正决定尾迹方向、升力相位和平均流向力的，是两种运动的幅值、频率和相位组合。

先抓重点

核心问题
为什么同一片翼用同样的振幅和频率，只改变俯仰相位，就能从制造阻力切换到产生推进力？
一句话直觉
升沉速度改变有效攻角，计算得到的尾迹速度改变面板环量，而 Kelvin 定理要求束缚环量的变化以相反环量进入尾迹；因此相位 ϕ\phi 同时重排力与尾迹。

核心数学模型

离散面板与尾迹反馈

(U∞+ub+uw−vwing)⋅ni=0,ΔpiΔsi=ρ[Ut,iΓi+Δsi∂∂t∑j≤iΓj](\mathbf U_\infty+\mathbf u_b+\mathbf u_w-\mathbf v_{\mathrm{wing}})\cdot\mathbf n_i=0,\qquad \Delta p_i\Delta s_i=\rho\left[U_{t,i}\Gamma_i+\Delta s_i\frac{\partial}{\partial t}\sum_{j\leq i}\Gamma_j\right]

四分之一面板放置束缚涡，四分之三面板执行无穿透。压强跳跃的环量变化率项是非定常项，并不等同于独立的纯附加质量。

升沉、俯仰与有效攻角

h(t)=h0sin⁡ωt,θ(t)=θ0sin⁡(ωt+ϕ),αeff=θ−tan⁡−1 ⁣(h˙U∞)h(t)=h_0\sin\omega t,\qquad \theta(t)=\theta_0\sin(\omega t+\phi),\qquad \alpha_{\mathrm{eff}}=\theta-\tan^{-1}\!\left(\frac{\dot h}{U_\infty}\right)

翼段向上运动时，它看到的相对来流向下偏转；俯仰可以增强或抵消这个运动攻角。

谐波参考：Theodorsen 尾迹记忆

C(k)≈1−0.165 ik0.0455+ik−0.335 ik0.3+ik,k=ωc2U∞C(k)\approx1-\frac{0.165\,ik}{0.0455+ik}-\frac{0.335\,ik}{0.3+ik},\qquad k=\frac{\omega c}{2U_\infty}

Jones–Wagner 有理近似把尾迹造成的幅值衰减和相位滞后压缩成一个可实时计算的复响应。

谐波参考：环量与附加质量升力

L′=πρb2 ⁣(−U∞θ˙−h¨+abθ¨)+2πρU∞b C(k) ⁣[U∞θ−h˙+b(12−a)θ˙]L'=\pi\rho b^2\!\left(-U_\infty\dot\theta-\ddot h+ab\ddot\theta\right)+2\pi\rho U_\infty b\,C(k)\!\left[U_\infty\theta-\dot h+b\left(\frac12-a\right)\dot\theta\right]

第一项是不需要稳态环量的加速度力，第二项是带尾迹记忆的环量力；这里半弦长为 b，转轴位置由 a 表示。

Kelvin 环量账本

Γb+∑jΓw,j=constant,ΔΓw=−ΔΓb\Gamma_b+\sum_j\Gamma_{w,j}=\mathrm{constant},\qquad \Delta\Gamma_w=-\Delta\Gamma_b

束缚环量每次变化，后缘就释放相反的增量；场景中的双色尾迹线就是这本离散账本。

常见难点

反向涡街不等于完整推进预测

容易误解尾迹看起来像推进涡街，因此效率读数就是现实飞行器的效率。

正确理解默认模式从薄涡片压强跳跃计算力，但没有前缘吸力修正、黏性与失速；参考模式使用升力投影和固定阻力代理。效率只作相位比较，既非现实飞行效率，也不保证落在零到一。

线性模型不会自动产生失速

容易误解把俯仰振幅推到最大，就能用同一个模型研究动态失速涡。

正确理解Theodorsen 理论假设小扰动和附着势流。峰值有效攻角超过界限时，Lab 会继续算以暴露外推趋势，但不会把它称为可信失速结果。

动手实验

  1. 01

    拆开两种运动

    先比较“纯升沉”和“纯俯仰”，观察尾迹环量、瞬时升力与平均流向力。

    应该观察到： 两种运动都制造交替环量。先等待完整周期，再比较实际计算的均值；开启谐波参考后，力与尾迹反馈的关系也会改变。
  2. 02

    组成推进相位

    选择“组合相位”，保持振幅不变，观察后缘包络、推力代理与功率代理。

    应该观察到： 俯仰可以抵消或增强运动攻角。比较整周期力和输入功率，不能仅凭尾迹形状判断推进。
  3. 03

    只反转相位

    切换到“反向相位”，不要改变频率和振幅。

    应该观察到： 尾迹组织和力的相位一起改变；记录完整周期均值，检查相位改变是否改变流向力的符号。

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