intuitive_cfd

第16章：非定常流动现象

当你站在风中，感受到阵阵吹来的风力时强时弱；当你看到旗帜在风中猎猎作响；当你听到电线在风中发出呜呜的哨音——这些都是非定常流动的表现。与我们之前讨论的定常流动不同，非定常流动中的流场参数随时间变化，这种时变特性带来了许多独特而重要的现象。从工程灾难（如塔科马海峡大桥的倒塌）到生物进化的奇迹（如昆虫的飞行），非定常流动塑造着我们的世界。本章将探讨这些动态现象背后的物理机制，以及如何在工程设计中考虑和利用它们。

16.1 涡激振动

16.1.1 涡脱落的物理机制

当流体绕过钝体（如圆柱）时，在一定的雷诺数范围内（约 $40 < Re < 10^5$），会发生周期性的涡脱落现象。这个过程可以这样理解：

流动方向 →
        ┌─────┐
        │     │  ← 圆柱
        │     │
        └─────┘
           ↓
    ╭──╮    ╭──╮
    │  │    │  │  ← 交替脱落的涡
    ╰──╯    ╰──╯
       ↘    ↙
        涡街

涡脱落的过程：

  1. 边界层分离：流体在圆柱背面分离，形成两个剪切层
  2. 剪切层卷起：上下剪切层交替卷起形成涡
  3. 涡的脱落：当涡增长到一定强度，它会从物体脱落
  4. 周期性重复：这个过程周期性地在上下两侧交替发生

16.1.2 斯特劳哈尔数

涡脱落频率由无量纲的斯特劳哈尔数（Strouhal number）描述：

\[St = \frac{fD}{V}\]

其中：

对于圆柱绕流，在很宽的雷诺数范围内，$St \approx 0.2$。这个经验值非常有用：

实例估算：一根直径10cm的烟囱在10m/s的风中： \(f = \frac{St \cdot V}{D} = \frac{0.2 \times 10}{0.1} = 20 \text{ Hz}\)

16.1.3 锁频现象

当结构的固有频率接近涡脱落频率时，会发生”锁频”（lock-in）现象：

  1. 频率捕获：涡脱落频率被”拉”到结构振动频率
  2. 振幅放大：结构振动与流动形成正反馈
  3. 共振区间：锁频发生在 $0.8f_n < f_v < 1.2f_n$ 范围内

锁频时的振幅可以用下式估算：

\[\frac{A}{D} \approx \frac{1}{2\zeta} \cdot \frac{\rho V^2 D^2}{\rho_s t^2 f_n^2 D^2}\]

其中 $\zeta$ 是阻尼比，$\rho_s$ 是结构密度，$t$ 是壁厚。

16.1.4 工程实例

海洋立管： 深海石油平台的立管在洋流作用下会发生涡激振动。设计时需要：

高层建筑和烟囱： 圆形截面的高层建筑特别容易发生涡激振动。常用对策：

经验法则：

16.2 颤振与驰振

16.2.1 颤振的物理机理

颤振（Flutter）是一种由气动力、弹性力和惯性力相互作用引起的自激振动。以机翼为例：

初始扰动 → 机翼扭转 → 攻角变化 → 升力变化
    ↑                              ↓
    └──── 能量输入 ←─── 相位关系合适 ←┘

颤振的关键在于运动与气动力之间的相位关系。当系统从流动中吸收的能量大于结构阻尼耗散的能量时，振动会不断增大。

16.2.2 临界颤振速度

经典的颤振速度估算（对于机翼）：

\[V_f \approx 2\pi f_\alpha \sqrt{\frac{I_\alpha}{\rho b^3 S_{C_L}}}\]

其中：

简化判据：

16.2.3 驰振的条件

驰振（Galloping）是另一种气动不稳定现象，主要发生在非流线型截面：

Glauert-Den Hartog判据： \(\frac{dC_L}{d\alpha} + C_D < 0\)

这意味着当升力系数随攻角的变化率（负值）超过阻力系数时，会发生驰振。

典型易发生驰振的截面：

16.2.4 工程实例

悬索桥设计： 现代悬索桥采用扁平箱梁截面，具有良好的气动稳定性：

输电线防舞动：

经验数据：

16.3 拍动翼空气动力学

16.3.1 昆虫飞行的独特机制

昆虫飞行打破了传统定常空气动力学的限制。一只蜜蜂按照定常理论计算，其翅膀面积完全不足以产生支撑体重的升力。关键在于非定常机制：

拍动周期：
     ↑ 上拍
    ╱ ╲     前缘涡
   ╱   ╲    附着
  ╱     ╲
 ╱_______╲  翻转
     ↓ 下拍

16.3.2 前缘涡机制

拍动翼在大攻角下产生强烈的前缘涡（Leading Edge Vortex, LEV）：

  1. 涡的形成：急剧的攻角变化在前缘产生分离涡
  2. 涡的稳定：翼尖流动和科氏力稳定前缘涡
  3. 升力增强：前缘涡产生的低压区大幅增加升力

升力系数可达： \(C_L \approx 2-4\) 远超定常流动的失速值（$C_L \approx 1.2$）

16.3.3 克拉普机制

“拍手-投掷”（Clap and Fling）机制，常见于小型昆虫：

  1. 拍手阶段：两翅在背部接触，挤出中间的流体
  2. 投掷阶段：两翅快速分开，形成强烈的环量

这种机制可以产生瞬时的高升力，特别适合悬停和起飞。

16.3.4 尺度效应

雷诺数对拍动翼性能影响巨大：

飞行者 翼展(mm) 频率(Hz) Re数 主导机制
果蝇 3 200 100 粘性主导
蜜蜂 15 230 1000 前缘涡
蜂鸟 100 25 5000 混合机制
鸽子 300 5 50000 准定常

16.3.5 仿生应用

微型飞行器（MAV）：

设计准则：

16.4 脉动流的特征

16.4.1 周期性流动参数

脉动流中，速度可表示为： \(V(t) = V_0 + V_1 \sin(\omega t)\)

关键无量纲参数：

沃默斯利数（Womersley number）： \(Wo = R\sqrt{\frac{\omega}{\nu}}\) 表征非定常惯性力与粘性力之比。

速度振幅比： \(\beta = \frac{V_1}{V_0}\)

16.4.2 准定常假设的适用性

准定常假设认为瞬时流场等同于相同瞬时条件下的定常流场。

适用条件：

当 $St > 1$ 时，流动完全由非定常效应主导。

16.4.3 相位滞后

在脉动流中，压力梯度与流量之间存在相位差：

\[\Delta \phi = \arctan\left(\frac{Wo^2}{8}\right)\]

这种相位差在：

16.4.4 生理流动实例

心血管系统：

呼吸系统：

设计人工器官时必须匹配这些特征时间尺度。

16.5 历史人物：塔科马海峡大桥倒塌事件（1940）

16.5.1 “舞动格蒂”的故事

1940年11月7日，华盛顿州塔科马海峡大桥在仅19m/s（约68km/h）的风速下发生剧烈振动并最终倒塌。这座桥因其建成后就经常在风中起伏而被昵称为”舞动格蒂”（Galloping Gertie）。

大桥的关键参数：

16.5.2 倒塌原因分析

长期以来，教科书将其归因于共振，但实际机理更复杂：

初期（建成-1940.11）：

倒塌当天：

关键因素：

  1. 截面形状：H型截面气动性能差
  2. 刚度不足：扭转刚度严重不足
  3. 气动耦合：扭转与竖弯模态耦合

16.5.3 工程教训

这次事故彻底改变了桥梁设计理念：

设计改进：

现代标准：

检验风速等级：
- 施工阶段：10年一遇
- 运营阶段：100年一遇
- 颤振检验：1.2×100年一遇

16.5.4 科学遗产

塔科马大桥事故推动了多个学科的发展：

  1. 气动弹性理论：Theodorsen和Scanlan的开创性工作
  2. 风工程学科：成为独立的工程分支
  3. 计算方法：CFD-CSD耦合分析
  4. 实验技术：节段模型和全桥气弹模型试验

这个事件提醒我们：大自然的力量常常超出我们的直觉，谦逊和严谨的科学态度是工程安全的基石。

16.6 高级话题：流固耦合与气动弹性

16.6.1 流固耦合的分类

根据耦合强度和方式，流固耦合（FSI）可分为：

单向耦合：

双向耦合：

强耦合vs弱耦合：

16.6.2 气动弹性三角

气动弹性现象由三种力的相互作用决定：

        空气动力
           /\
          /  \
         /    \
        /  AE  \
       /________\
    弹性力    惯性力

不同区域的现象：

16.6.3 降阶模型方法

对于复杂的FSI问题，降阶模型（ROM）提供了高效的解决方案：

本征正交分解（POD）： 将流场分解为基本模态： \(\mathbf{u}(\mathbf{x},t) = \bar{\mathbf{u}}(\mathbf{x}) + \sum_{i=1}^{N} a_i(t)\boldsymbol{\phi}_i(\mathbf{x})\)

系统识别方法：

16.6.4 工程应用实例

飞机机翼设计：

风力发电机叶片：

高速列车受电弓：

16.6.5 数值方法进展

分区耦合策略：

CFD求解器 ←→ 接口 ←→ CSD求解器
    ↓         ↓         ↓
  流场网格   数据传递   结构网格

时间推进格式：

网格变形技术：

16.7 本章小结

非定常流动现象展现了流体力学的动态美感和复杂性。本章的关键要点：

核心概念

  1. 涡激振动：
    • 斯特劳哈尔数 $St \approx 0.2$ 对圆柱
    • 锁频发生在 $0.8f_n < f_v < 1.2f_n$
    • 螺旋列板可有效抑制
  2. 颤振与驰振：
    • 颤振：自激振动，需要气动-弹性耦合
    • 驰振：Glauert-Den Hartog判据
    • 临界风速是设计关键
  3. 拍动翼机制：
    • 前缘涡提供高升力（$C_L = 2-4$）
    • 尺度效应决定主导机制
    • 为MAV设计提供灵感
  4. 脉动流特征：
    • 沃默斯利数判断流态
    • 准定常假设：$St < 0.1$
    • 相位关系影响能量传递

关键无量纲数

参数 定义 物理意义 临界值
$St$ $fD/V$ 非定常性 0.1-1
$Wo$ $R\sqrt{\omega/\nu}$ 脉动流特征 1, 10
$k$ $\omega L/V$ 减缩频率 0.1
$V_r$ $V/(f_n D)$ 约化速度 5

工程启示

16.8 练习题

基础题

16.1 涡脱落频率计算 一根直径20cm的圆形烟囱在风速15m/s的环境中，计算涡脱落频率。如果烟囱的固有频率为18Hz，是否会发生锁频现象？

提示 使用斯特劳哈尔数 St = 0.2 计算涡脱落频率，然后检查是否在锁频范围内。
答案 涡脱落频率： $$f = \frac{St \cdot V}{D} = \frac{0.2 \times 15}{0.2} = 15 \text{ Hz}$$ 锁频范围：$0.8 \times 18 = 14.4$ Hz 到 $1.2 \times 18 = 21.6$ Hz 由于15 Hz在[14.4, 21.6]范围内，会发生锁频现象，振动会被放大。

16.2 昆虫翅膀雷诺数 一只蜜蜂翅膀长度15mm，拍动频率230Hz，飞行速度3m/s。估算其雷诺数并判断主导的流动机制。（运动粘度 $\nu = 1.5 \times 10^{-5}$ m²/s）

提示 特征速度应该考虑拍动速度而非飞行速度。拍动速度 = 2πfA，其中A是拍动幅度。
答案 拍动幅度约为翼长：A ≈ 15mm = 0.015m 拍动速度：$V_{flap} = 2\pi f A = 2\pi \times 230 \times 0.015 \approx 21.7$ m/s 这远大于飞行速度，所以用拍动速度计算： $$Re = \frac{V_{flap} \cdot L}{\nu} = \frac{21.7 \times 0.015}{1.5 \times 10^{-5}} \approx 2170$$ 在Re ≈ 1000-5000范围，前缘涡机制主导。

16.3 脉动流相位差 在直径2cm的血管中，血流脉动频率为1.2Hz（心率72次/分），血液运动粘度 $\nu = 4 \times 10^{-6}$ m²/s。计算沃默斯利数并估算压力梯度与流速的相位差。

提示 先计算角频率ω，然后用Wo数公式，最后估算相位差。
答案 角频率：$\omega = 2\pi f = 2\pi \times 1.2 = 7.54$ rad/s 半径：R = 0.01 m 沃默斯利数： $$Wo = R\sqrt{\frac{\omega}{\nu}} = 0.01 \sqrt{\frac{7.54}{4 \times 10^{-6}}} = 13.7$$ 相位差： $$\Delta \phi = \arctan\left(\frac{Wo^2}{8}\right) = \arctan\left(\frac{13.7^2}{8}\right) = \arctan(23.5) \approx 87°$$ 压力梯度领先流速约87°，接近90°，说明惯性效应主导。

挑战题

16.4 悬索桥颤振分析 某悬索桥主梁宽30m，高3m，质量线密度8000 kg/m，扭转惯量 $I_\alpha = 1.2 \times 10^7$ kg·m²/m，扭转频率0.3Hz。估算颤振临界风速。（空气密度1.2 kg/m³）

提示 使用简化的颤振速度公式，考虑质量比和频率的影响。对于扁平箱梁，升力线斜率约为2π。
答案 半宽度 b = 15m 质量比：$\mu = \frac{m}{\rho b^2} = \frac{8000}{1.2 \times 15^2} = 29.6$ 颤振速度估算： $$V_f \approx 2\pi f_\alpha b \sqrt{\frac{\mu}{2\pi}} = 2\pi \times 0.3 \times 15 \times \sqrt{\frac{29.6}{2\pi}}$$ $$V_f \approx 28.3 \times 2.17 = 61.4 \text{ m/s}$$ 这相当于221 km/h的风速，属于超强台风级别。

16.5 涡激振动抑制设计 某海洋立管直径0.5m，长度100m，在1m/s的洋流中工作。设计螺旋列板参数来抑制涡激振动。分析列板对阻力的影响。

提示 螺旋列板的标准设计：螺距5D，高度0.1D，三股螺旋。阻力会增加20-30%。
答案 螺旋列板设计： - 螺距：5 × 0.5 = 2.5 m - 列板高度：0.1 × 0.5 = 0.05 m - 螺旋数：3股，相位差120° - 螺旋角：$\theta = \arctan(2.5/(\pi \times 0.5)) \approx 58°$ 效果分析： - 破坏涡的相关长度，从100m降至约2.5m - 振幅减少80-90% - 阻力系数从1.2增至约1.5（增加25%） - 年疲劳损伤降低100倍以上 权衡：阻力增加的代价远小于疲劳寿命延长的收益。

16.6 拍动翼MAV设计 设计一个翼展10cm的微型飞行器，重量5g。基于昆虫飞行原理，确定拍动频率、拍动幅度和所需功率。

提示 使用频率-尺度关系 f ∝ L^(-1)，参考蜜蜂的数据进行缩放。功率密度约为30-50 W/kg。
答案 基于尺度律（参考蜜蜂）： - 蜜蜂：翼展30mm，频率230Hz - MAV缩放：$f = 230 \times (30/100) = 69$ Hz 设计参数： - 拍动频率：70 Hz - 拍动幅度：±60°（约±50mm行程） - 翼面积：约25 cm²（双翼） - 展弦比：4 升力估算： - 需要升力：5g × 9.8 = 0.049 N - 升力系数（时均）：$\bar{C_L} \approx 1.5$ - 验证速度：$V = \sqrt{\frac{2mg}{\rho S C_L}} \approx 2.5$ m/s 功率需求： - 比功率：40 W/kg - 总功率：40 × 0.005 = 0.2 W - 考虑效率（30%）：需要0.7 W输入功率 这在现有电池和电机技术范围内可实现。

16.7 开放思考题：流固耦合的未来 随着计算能力的提升和新材料的发展，讨论以下问题：

  1. 智能材料（如形状记忆合金、压电材料）如何改变流动控制策略？
  2. 机器学习在预测和控制非定常流动中的潜力？
  3. 仿生设计中，哪些生物的流固耦合机制最值得深入研究？
思考方向 1. 智能材料应用： - 自适应翼型：根据流动状态改变形状 - 主动涡发生器：按需激活 - 智能蒙皮：分布式传感和致动 2. 机器学习潜力： - 降阶模型的自动构建 - 非线性系统识别 - 实时流动状态预测和控制 3. 仿生研究重点： - 蝙蝠翼膜的被动适应性 - 鱼类侧线系统的流动感知 - 鸟类羽毛的多尺度流动控制

16.9 常见陷阱与错误

概念误区

  1. 共振 ≠ 颤振
    • 错误：塔科马大桥是共振导致的倒塌
    • 正确：颤振是自激振动，能量来自流动
  2. 涡脱落频率不是常数
    • 错误：固定St数适用于所有情况
    • 正确：St数随Re变化，特别在临界Re附近
  3. 准定常假设的滥用
    • 错误：拍动翼可以用准定常理论分析
    • 正确：St > 0.1时必须考虑非定常效应

计算错误

  1. 特征速度选择
    • 错误：总是使用来流速度
    • 正确：拍动问题用拍动速度，可能远大于飞行速度
  2. 相位关系忽略
    • 错误：脉动流中压力和流速同相
    • 正确：高Wo数时相位差接近90°
  3. 尺度效应忽视
    • 错误：大尺度经验直接用于微尺度
    • 正确：雷诺数相似才能保证机理相同

设计陷阱

  1. 单一模态考虑
    • 错误：只检查基频共振
    • 正确：高阶模态和模态耦合同样危险
  2. 阻尼的过度依赖
    • 错误：增加阻尼总是有益的
    • 正确：过大阻尼可能激发其他不稳定性
  3. 风洞试验的局限
    • 错误：缩尺模型完全代表实际
    • 正确：需要考虑雷诺数效应和弹性相似

16.10 最佳实践检查清单

设计阶段

分析验证

控制措施

文档记录

记住：非定常流动现象往往是灾难性失效的根源，但也可能是创新设计的灵感。理解其物理本质，合理利用和控制，是流体工程师的核心能力。