第6章：势流理论的工程应用

当你观察河水绕过桥墩时的流动模式，或者看到旗帜在风中飘扬时在旗杆后方形成的规律涡旋，你正在目睹势流理论描述的现象。尽管真实流体都有粘性，但在许多工程问题中，势流理论提供的无粘流动近似却能给出惊人准确的结果。本章将探讨如何用简单的数学工具——叠加原理，来构建和理解复杂的流场，以及这些理论如何帮助我们设计更好的飞机、桥梁和建筑物。

6.1 引言与学习目标

势流理论的工程价值

势流理论假设流体无粘、无旋，这看似是对真实流体的过度简化，但它在以下工程场景中仍然极其有用：

  1. 流线型物体的设计：当雷诺数很高且没有流动分离时，势流理论能准确预测压力分布。这就是为什么NACA翼型系列的设计最初就基于势流理论。

  2. 初步设计阶段：在详细CFD分析之前，势流理论提供快速的定性理解。比如评估不同机翼平面形状的升力分布。

  3. 流动控制策略：理解基本流动元素（源、汇、涡）的特性，有助于设计涡流发生器、吹吸控制等流动控制装置。

  4. 非定常流动分析：拍动翼、振动结构周围的流场，在附加质量效应占主导时，势流理论提供了关键见解。

本章学习目标

完成本章学习后，你将能够：

  • 使用叠加原理构建复杂流场，并理解各基本流动元素的物理意义
  • 分析圆柱绕流问题，预测升力和涡脱落频率
  • 应用镜像法处理壁面附近的流动问题
  • 理解卡门涡街的形成机理及其工程影响
  • 掌握保角变换的基本概念，了解如何将圆柱绕流解推广到任意形状

6.2 叠加原理与复杂流场构造

线性叠加的数学基础

势流理论的核心优势在于拉普拉斯方程的线性特性：

$$\nabla^2 \phi = 0$$ 如果 $\phi_1$ 和 $\phi_2$ 都满足拉普拉斯方程，那么它们的任意线性组合 $\phi = a\phi_1 + b\phi_2$ 也满足该方程。这意味着我们可以将复杂流场分解为简单流动的叠加。

想象一下：就像用几个基本的乐高积木可以搭建复杂的结构，我们用几个基本流动元素就能构造出飞机周围的整个流场。

基本流动元素的组合

最常用的基本流动元素包括：

  1. 均匀流 ($\phi = V_\infty x$)：代表远场来流
  2. 点源/汇 ($\phi = \pm \frac{m}{2\pi} \ln r$)：模拟喷射或吸入
  3. 点涡 ($\phi = \frac{\Gamma}{2\pi} \theta$)：产生环量和升力
  4. 偶极子 ($\phi = \frac{\mu \cos\theta}{2\pi r}$)：源汇对的极限情况

通过组合这些元素，我们能构造出：

  • 兰金半体：均匀流 + 点源（模拟钝头体）
  • 兰金卵形体：均匀流 + 源汇对（模拟流线型物体）
  • 有升力圆柱：均匀流 + 偶极子 + 点涡（马格努斯效应）

工程应用实例

案例1：风力发电机叶片设计

风力发电机的叶片可以用一系列沿叶片分布的束缚涡来模拟。每个涡的强度决定了局部升力，通过调整涡强度分布，我们能优化叶片的扭转角和弦长分布，使功率输出最大化。

     来流 V∞
        ↓
    ┌───┼───┐  
    │   ↓   │  束缚涡 Γ(r)
    │  ~~~  │  
    │   ↓   │  
    └───┼───┘
        ↓
    尾流涡系

案例2：潜艇外形优化

潜艇的流线型外形可以用分布的源汇来表示。通过调整源汇强度分布，我们能：

  • 最小化形状阻力
  • 控制压力恢复，避免流动分离
  • 优化推进器入流条件

案例3：建筑物风载评估

高层建筑可简化为垂直线上的偶极子分布，地面用镜像法处理。这种方法能快速估算：

  • 不同高度的风压分布
  • 建筑群之间的干扰效应
  • 行人高度的风环境

实用叠加技巧

  1. 先远后近：先确定远场边界条件（均匀流），再添加局部修正（源、涡等）

  2. 对称利用：对于对称问题，只需构造一半流场，另一半由对称性确定

  3. 迭代修正：初始近似 → 检查边界条件 → 添加修正项 → 重复

  4. 物理直觉：每个流动元素都应有明确的物理意义，如模拟厚度效应用源汇，模拟升力用环量

6.3 偶极子、源、汇的物理意义

点源与点汇

点源代表流体从一点向外均匀喷射，就像：

  • 水下的气泡释放点
  • 喷泉的喷嘴
  • 爆炸的中心

其流函数和势函数为： $$\psi = \frac{m}{2\pi}\theta, \quad \phi = \frac{m}{2\pi}\ln r$$ 其中 $m$ 是源强度（体积流量）。径向速度 $v_r = m/(2\pi r)$，切向速度为零。

点汇则相反，代表流体向一点汇聚，如：

  • 排水口
  • 龙卷风的底部（简化模型）
  • 真空吸尘器的吸入口

点汇的数学表达式与点源相同，只是 $m$ 为负值。

工程意义：在飞艇设计中，我们用分布的源来模拟飞艇的体积效应。源强度分布决定了飞艇的外形。类似地，在进气道设计中，汇可以模拟发动机的吸入效应。

偶极子的形成与特性

偶极子是一对强度相等、符号相反的源汇无限接近时的极限情况。想象把一个小磁铁的南北极无限接近，但保持磁矩不变。 $$\phi = \frac{\mu \cos\theta}{2\pi r} = \frac{\mu x}{2\pi(x^2+y^2)}$$ 其中 $\mu$ 是偶极矩（$\mu = m \cdot d$，当 $d \to 0$, $m \to \infty$）。

偶极子的流线形态像一组相切的圆：

      ↑ y
      │
   ○  │  ○     上半平面：逆时针圆
  ○   │   ○    
──────┼──────→ x
  ○   │   ○    
   ○  │  ○     下半平面：顺时针圆
      │

实际流动中的对应关系

  1. 圆柱绕流 = 均匀流 + 偶极子

当均匀流遇到圆柱时，圆柱的存在效应可以用一个偶极子来表示： $$\phi = V_\infty r \cos\theta + \frac{V_\infty R^2 \cos\theta}{r}$$ 这给出了完美的无粘圆柱绕流解，前后驻点处速度为零，上下最宽处速度达到 $2V_\infty$。

  1. 三维球体 = 均匀流 + 三维偶极子

类似地，球体绕流可用三维偶极子表示。这解释了为什么雨滴在高速下落时会变形——前驻点高压使雨滴变扁。

  1. 机翼的厚度效应

薄翼型的厚度分布可以用沿弦向分布的源汇来模拟。NACA四位数字翼型族就是这样设计的：

  • 源分布决定厚度
  • 涡分布决定弯度和升力
  1. 地效飞行器

当飞行器接近地面时，其诱导的下洗流被地面"反射"，可用镜像偶极子系统来分析。这解释了为什么地效飞行器能以更小的功率维持飞行。

物理直觉培养

记住这个类比：

  • 源 = 充气（增加体积）
  • 汇 = 放气（减少体积）
  • 偶极子 = 推开流体（物体存在）
  • 涡 = 旋转（产生升力）

通过这些基本元素的组合，我们能够直观理解复杂流场。比如，潜望镜伸出水面的流动：

  • 垂直圆柱部分 → 偶极子（推开水流）
  • 顶部整流罩 → 源（使流线外扩）
  • 尾流区 → 汇（补充动量亏损）

6.4 圆柱绕流与卡门涡街

无环量圆柱绕流

理想的无粘、无环量圆柱绕流是势流理论的经典问题。通过叠加均匀流和偶极子，我们得到： $$\phi = V_\infty \left(r + \frac{R^2}{r}\right)\cos\theta$$ 这产生了著名的达朗贝尔佯谬：圆柱不受任何阻力！表面压力分布完全对称： $$C_p = 1 - 4\sin^2\theta$$ 在前后驻点 ($\theta = 0, \pi$)，$C_p = 1$；在最宽处 ($\theta = \pm\pi/2$)，$C_p = -3$。

为什么现实不同？ 真实流体的粘性导致边界层分离，破坏了前后对称性，产生压差阻力。这就像理论告诉我们"完美的球应该没有空气阻力"，但高尔夫球的凹坑恰恰是为了减少阻力而设计的。

有环量圆柱绕流（马格努斯效应）

给圆柱添加环量 $\Gamma$（比如让圆柱旋转），势函数变为： $$\phi = V_\infty \left(r + \frac{R^2}{r}\right)\cos\theta + \frac{\Gamma}{2\pi}\theta$$ 这产生了垂直于来流的升力（马格努斯力）： $$L = \rho V_\infty \Gamma$$ 日常例子：

  • 乒乓球的弧圈球：上旋使球下沉，下旋使球上飘
  • 足球的香蕉球：侧旋产生横向力，使球轨迹弯曲
  • 棒球的曲球：投手通过控制球的旋转来迷惑打者
  • Flettner转子船：利用旋转圆柱代替帆，节能环保

驻点位置随环量变化：

  • 无环量：前后驻点在 $\theta = 0, \pi$
  • 小环量：驻点向下移动
  • 临界环量 ($\Gamma = 4\pi V_\infty R$)：两驻点在底部汇合
  • 大环量：驻点离开圆柱表面

卡门涡街的形成机理

当 $Re > 40$ 时，圆柱后方开始周期性地脱落反向旋转的涡，形成卡门涡街。这种现象普遍存在：

    流向 →

    ┌─┐     ↻     ↺     ↻     ↺
    │ │  ───────────────────────→
    └─┘     ↺     ↻     ↺     ↻
    圆柱   
           交替脱落的涡（卡门涡街）

形成过程：

  1. 低雷诺数 ($Re < 5$)：流动附着，无分离
  2. $5 < Re < 40$：形成稳定的分离泡
  3. $40 < Re < 200$：层流涡街，规律脱落
  4. $200 < Re < 3×10^5$：湍流涡街，脱落变得不规则
  5. $Re > 3×10^5$：边界层转捩，涡街重新规律化

斯特劳哈尔数与涡脱落频率

涡脱落频率由斯特劳哈尔数表征： $$St = \frac{fD}{V_\infty}$$ 其中 $f$ 是脱落频率，$D$ 是圆柱直径。对于圆柱：

  • $St ≈ 0.2$ 在很宽的雷诺数范围内（$200 < Re < 10^5$）

工程应用：

  1. 烟囱和冷却塔设计 - 涡脱落频率：$f = 0.2 V/D$ - 如果接近结构固有频率 → 共振 → 破坏 - 解决方案：螺旋条板、调谐质量阻尼器

  2. 潜望镜和天线 - 潜艇潜望镜的"唱歌"现象 - 汽车天线的振动噪音 - 解决：螺旋缠绕、变截面设计

  3. 电力线的舞动 - 风致振动导致短路 - 防振锤、间隔棒的设计原理

  4. 斜拉桥缆索 - 雨振现象（水线改变截面形状） - 螺旋线、凹坑表面处理

经验法则：

  • 估算涡脱落频率：$f ≈ 0.2V/D$（圆柱）
  • 避免共振：结构固有频率应远离涡脱落频率（至少相差30%）
  • 抑制方法：
  • 破坏展向相关性（螺旋条板）
  • 改变截面形状（椭圆、方形带倒角）
  • 主动控制（吹吸、振动）

涡激振动的控制

当涡脱落频率接近结构固有频率时，会发生"锁频"现象，振幅急剧增大。1940年塔科马海峡大桥的倒塌就是典型案例。

被动控制方法：

  • 螺旋列板：破坏涡的展向相关性，减少升力脉动
  • 整流罩：改变钝体为流线型
  • 分流板：阻止剪切层相互作用
  • 多孔罩：减弱涡强度

主动控制方法：

  • 合成射流：在分离点注入动量
  • 等离子体激励：改变边界层特性
  • 自适应外形：根据流动状态改变几何

6.5 镜像法与壁面效应

镜像法原理

镜像法是处理壁面边界条件的优雅技巧。核心思想：在壁面另一侧放置"镜像"流动元素，使壁面成为流线（满足无穿透条件）。

基本规则：

  • 源 → 镜像源（强度相同）
  • 汇 → 镜像汇（强度相同）
  • 顺时针涡 → 逆时针镜像涡（强度相反）

例如，壁面附近的点涡：

    真实涡 Γ ↻        
         │            
    ─────┼───────────  壁面（y=0）
         │            
    镜像涡 -Γ ↺

合成流场在壁面处只有切向速度，垂直速度为零，完美满足边界条件。

地面效应

当飞行器接近地面时，地面的存在显著改变流场，这就是地面效应。

机翼的地面效应：

高度 $h$ 处的机翼可用镜像法分析：

  • 真实机翼：束缚涡 $\Gamma$，位于 $(0, h)$
  • 镜像机翼：束缚涡 $-\Gamma$，位于 $(0, -h)$

效应：

  1. 诱导阻力减小：下洗减弱，诱导阻力可减少40-50%
  2. 升力增加：地面阻挡下洗，等效迎角增大
  3. 俯仰力矩变化：压力中心前移

经验法则：

  • 显著地效区：$h/b < 0.25$（$h$=高度，$b$=翼展）
  • 强地效区：$h/b < 0.1$
  • 地效增升比：$\Delta C_L/C_L ≈ 0.3(b/h)$（$h/b < 0.3$时）

工程应用：

  1. 地效飞行器（WIG） - 里海怪物（苏联ekranoplan） - 利用地效可载重500吨，巡航高度2-5米 - 升阻比可达25-30（普通飞机15-20）

  2. 赛车的地面效应 - F1赛车底部文丘里通道 - 产生向下的"负升力"增加抓地力 - 离地间隙的微小变化会导致下压力剧变

  3. 直升机的地面效应 - 悬停功率在地效内可减少15-25% - IGE（地效内）vs OGE（地效外）性能图表 - 起降时的"地面共振"问题

壁面干扰修正

风洞测试中，洞壁的存在会影响测量结果，需要修正。

闭口风洞修正：

实体壁面的堵塞效应使有效速度增加： $$V_{eff} = V_\infty(1 + \epsilon)$$ 其中堵塞系数： $$\epsilon = \frac{K \cdot S_{model}}{S_{tunnel}}$$

  • 二维翼型：$K ≈ 0.25$
  • 三维机翼：$K ≈ 0.15$
  • 钝体：$K ≈ 0.4-0.6$

开口风洞修正：

射流边界的影响相反，需要不同的修正公式。

角区效应：

风洞角区的二次流会影响：

  • 边界层发展
  • 转捩位置
  • 分离特性

修正方法：

  • 角区吸气
  • 涡流发生器
  • 填角处理

管道和腔体中的镜像

管道流动：

圆管中的点涡需要无穷多个镜像才能满足壁面条件。实践中常用：

  • 前3-5个镜像（误差<1%）
  • 格林函数方法
  • 数值保角变换

腔体共振：

方腔中的振荡流会激发特定模态：

┌─────────┐
│ n=1,1   │  基频
│  ╱╲╱╲   │  
└─────────┘

┌─────────┐
│ n=2,1   │  二次谐波
│ ╱╲╱╲╱╲  │  
└─────────┘

应用：

  • 汽车天窗啸叫
  • 飞机弹舱气动噪声
  • 管风琴和长笛的发声原理

多体干扰的镜像处理

双圆柱干扰：

两个平行圆柱的相互作用可用镜像涡系分析：

  • 串列：尾流干扰，阻力变化
  • 并列：间隙流加速，可能吸引
  • 交错：复杂的三维效应

建筑群风场：

城市建筑群的风场评估：

  1. 每栋建筑用偶极子/源模拟
  2. 地面用镜像法处理
  3. 叠加得到总流场
  4. 识别危险区域（下冲流、角区加速）

实用设计准则：

  • 建筑间距 > 2倍特征尺寸，干扰可忽略
  • 通风廊道宽度 > 建筑高度的0.3倍
  • 裙楼可减缓下冲流影响

6.6 历史人物：西奥多·冯·卡门与涡街的发现

卡门的生平与贡献

西奥多·冯·卡门（Theodore von Kármán, 1881-1963）是20世纪最伟大的流体力学家之一。生于匈牙利布达佩斯的犹太家庭，他的一生横跨两次世界大战，见证并推动了航空航天工业从莱特兄弟的木制双翼机到超声速飞行的巨大飞跃。

主要成就：

  • 1911年：发现并解释卡门涡街现象
  • 1921年：提出边界层动量积分方程（卡门-波尔豪森方法）
  • 1930年：创立湍流的统计理论基础
  • 1946年：预见超声速飞行时代，提出"卡门线"（100km高度，太空边界）

卡门不仅是理论家，更是工程实践的推动者。他创立了加州理工学院的古根海姆航空实验室（GALCIT），后来发展成NASA喷气推进实验室（JPL）。他的学生包括钱学森、郭永怀等中国航天事业的奠基人。

涡街现象的发现过程

1908年，年轻的卡门在哥廷根大学普朗特的实验室工作。一个看似简单的问题困扰着他：为什么圆柱后的涡会交替脱落？

关键洞察：

卡门没有满足于观察，而是问了一个深刻的问题："什么样的涡排列是稳定的？"通过理论分析，他发现：

  1. 单排涡列不稳定：会迅速扩散消失
  2. 对称双排也不稳定：微小扰动就会破坏
  3. 交错排列可以稳定：但只在特定间距比 $h/l = 0.281$ 时

其中 $h$ 是两排涡的横向间距，$l$ 是同排涡的纵向间距。这个神奇的比值来自动量守恒和能量考虑。

实验验证：

卡门用一个优雅的实验证实理论：

     水槽
    ┌────────────────┐
    │  →→→→→→→→→→→   │
    │     ┃          │  圆柱
    │  →→→┃→→→→→→   │
    │     ┃  ↻ ↺ ↻  │  涡街
    └────────────────┘
      铝粉显示流线

他发现自然形成的涡街恰好满足 $h/l ≈ 0.28$！这种理论预测与实验的完美吻合，成为流体力学史上的经典案例。

对现代工程的影响

卡门涡街的发现不仅是学术成就，更带来了深远的工程影响：

  1. 结构设计革命

1940年塔科马海峡大桥倒塌后，工程界才真正认识到涡激振动的危险。今天，任何暴露在风中的细长结构都必须考虑卡门涡街：

  • 摩天大楼的阻尼器设计
  • 海洋平台的立管疲劳分析
  • 输电线路的防振设计
  1. 流量测量技术

涡街流量计基于斯特劳哈尔数的恒定性：

  • 精度高（±1%）
  • 无活动部件，可靠性好
  • 广泛用于石化、电力行业
  1. 仿生学应用

鱼类游动时主动利用涡街：

  • 鱼群中的个体利用前鱼的涡街节省能量
  • 仿生推进器设计借鉴这一原理
  • 风力发电场的优化布置

卡门的哲学：

卡门常说："科学家研究已存在的世界，工程师创造从未存在的世界。"他强调：

  1. 理论必须为实践服务："没有应用的理论是空洞的，没有理论的应用是盲目的"

  2. 跨学科思维：流体力学问题常需要数学、物理、工程的综合

  3. 国际合作：他创立了国际理论与应用力学联合会（IUTAM），促进全球科学交流

轶事：

  • 卡门是出了名的"夜猫子"，常在深夜2-3点工作，早上11点才起床
  • 他的办公室总是烟雾缭绕（大烟斗不离手），学生们戏称要在他办公室讨论问题需要"防毒面具"
  • 尽管是理论大师，他坚持亲自做实验，认为"手感"对理解流动至关重要

卡门的遗产不仅是涡街理论，更是一种将严谨理论与工程直觉完美结合的研究方法。正如他的墓志铭所言："他破解了湍流的密码，让人类飞得更高更远。"

6.7 高级话题：保角变换与复势理论

复变函数方法

二维势流的优雅之处在于可以用复变函数统一处理。定义复势： $$W(z) = \phi + i\psi$$ 其中 $z = x + iy$，$\phi$ 是速度势，$\psi$ 是流函数。复速度： $$\frac{dW}{dz} = u - iv$$ 基本流动的复势：

  • 均匀流：$W = V_\infty z e^{-i\alpha}$（$\alpha$ 是迎角）
  • 源/汇：$W = \frac{m}{2\pi}\ln z$
  • 点涡：$W = -\frac{i\Gamma}{2\pi}\ln z$
  • 偶极子：$W = \frac{\mu}{2\pi z}$

复变函数的威力在于：

  1. 自动满足拉普拉斯方程（解析函数的性质）
  2. 保角变换保持势流特性
  3. 奇点理论提供物理洞察

儒科夫斯基变换

儒科夫斯基变换是航空史上最重要的数学工具之一： $$\zeta = z + \frac{c^2}{z}$$ 它将圆变换成翼型！

变换过程：

   z平面（圆）           ζ平面（翼型）
      ○                    ⟵===⟶
   偏心圆               对称翼型

   过(-c,0)的圆          平板

   不过(-c,0)的圆        钝尾缘翼型

关键参数：

  • 圆心位置 $(x_0, y_0)$ 决定厚度和弯度
  • 圆半径 $R$ 决定翼型尺寸
  • 后缘条件：圆必须过点 $z = -c$

库塔条件的实现：

环量通过库塔条件确定——后缘速度有限： $$\Gamma = 4\pi V_\infty R \sin(\alpha + \beta)$$ 其中 $\beta = \arcsin(y_0/R)$。

任意翼型的势流解

薄翼理论（线化理论）：

对于小弯度、小厚度翼型，可将翼型表示为：

  • 弦线上的涡分布 $\gamma(x)$（产生升力）
  • 源汇分布 $q(x)$（产生厚度）

升力系数： $$C_L = 2\pi\alpha + 2\pi\int_0^c \frac{dy_c}{dx}\frac{dx}{c}$$ 其中 $y_c(x)$ 是中弧线。

面元法：

现代方法将翼型表面离散成小面元，每个面元上放置：

  • 源（满足厚度条件）
  • 涡（满足环量条件）

通过满足表面无穿透条件，得到线性方程组，求解得到压力分布。

工程应用实例

  1. 多段翼型设计

高升力系统（襟翼、缝翼）可用多个儒科夫斯基翼型叠加：

  • 主翼：大环量，提供主要升力
  • 襟翼：后置小翼，增加弯度
  • 缝翼：前缘小翼，延迟失速
  1. 螺旋桨叶片设计

叶素理论结合保角变换：

  • 每个半径站位用不同翼型
  • 考虑诱导速度的影响
  • 优化扭转分布
  1. 风洞壁面修正

保角变换将有限风洞映射到无限流场： $$z_{corrected} = z_{tunnel} \cdot f(S_{model}/S_{tunnel})$$

数值实现技巧

快速多极子方法（FMM）：

  • 将远场影响分组计算
  • 复杂度从 $O(N^2)$ 降到 $O(N\log N)$
  • 可处理10万+面元

自适应网格：

  • 前缘、后缘加密
  • 分离区自动识别
  • 误差估计与网格调整

复势理论的局限性

  1. 粘性效应：无法预测分离、失速
  2. 三维效应：只适用于二维或准二维流动
  3. 可压缩效应：需要Prandtl-Glauert修正
  4. 非定常效应：需要扩展到非定常势流理论

现代发展

共形映射与机器学习：

  • 神经网络学习最优变换
  • 自动化翼型优化
  • 逆设计问题

与CFD的结合：

  • 势流提供初始解
  • 边界层修正
  • 粘性-无粘耦合方法

记住：保角变换不仅是数学技巧，更是理解流动拓扑结构的强大工具。通过将复杂几何变换到简单区域，我们能洞察流动的本质特征。

6.8 本章小结

势流理论虽然忽略了粘性，但它为我们提供了理解流体运动的强大框架。通过本章学习，我们掌握了：

核心概念：

  • 叠加原理：复杂流场 = 简单流动元素的线性组合
  • 基本元素：源（体积生成）、汇（体积消失）、偶极子（物体存在）、涡（环量和升力）
  • 镜像法：优雅处理壁面边界条件的数学技巧
  • 卡门涡街：周期性涡脱落现象，$St ≈ 0.2$ 的普适性

关键公式：

  • 复势：$W = \phi + i\psi$
  • 马格努斯力：$L = \rho V_\infty \Gamma$
  • 儒科夫斯基变换：$\zeta = z + c^2/z$
  • 涡脱落频率：$f = St \cdot V/D$

工程智慧：

  1. 势流理论在高雷诺数、无分离流动中惊人准确
  2. 地面效应可显著改善飞行器性能（诱导阻力减少40-50%）
  3. 涡激振动是细长结构的主要威胁，需要仔细设计避免共振
  4. 保角变换不仅是数学工具，更是洞察流动本质的窗口

实用经验法则：

  • 圆柱涡脱落：$St ≈ 0.2$（$200 < Re < 10^5$）
  • 地效显著区：$h/b < 0.25$
  • 风洞堵塞修正：$\epsilon ≈ 0.25 \times S_{model}/S_{tunnel}$（二维）
  • 建筑间距 > 2倍特征尺寸可忽略干扰

势流理论如同流体力学的"理想国"——虽然现实更加复杂，但这个理想化的世界让我们看清了流动的骨架和脉络。

6.9 练习题

基础题（1-4题）

题1：一个半径 $R = 1m$ 的圆柱在 $V_\infty = 10 m/s$ 的均匀流中。计算： a) 无环量时，最高点的压力系数 b) 如果圆柱以 500 rpm 顺时针旋转，估算升力（空气密度 $\rho = 1.2 kg/m^3$）

提示：使用 $C_p = 1 - 4\sin^2\theta$ 和马格努斯力公式

答案

a) 最高点 $\theta = \pi/2$：$C_p = 1 - 4 = -3$

b) 旋转产生环量：$\Gamma = 2\pi R \cdot \omega R = 2\pi \times 1 \times (500 \times 2\pi/60) \times 1 = 52.4 m^2/s$

升力：$L = \rho V_\infty \Gamma = 1.2 \times 10 \times 52.4 = 629 N/m$

题2：烟囱直径 $D = 2m$，风速 $V = 20 m/s$。估算： a) 涡脱落频率 b) 如果烟囱固有频率为 2.5 Hz，是否会发生共振？

提示：使用斯特劳哈尔数 $St = 0.2$

答案

a) $f = St \cdot V/D = 0.2 \times 20/2 = 2 Hz$

b) 涡脱落频率 2 Hz 与固有频率 2.5 Hz 相差 25%，接近共振区域，可能发生涡激振动，建议安装螺旋板或调谐质量阻尼器。

题3：机翼翼展 $b = 10m$，在离地高度 $h = 2m$ 飞行。估算地面效应导致的诱导阻力减少百分比。

提示：地效强度与 $h/b$ 相关

答案

$h/b = 2/10 = 0.2$，处于显著地效区。

诱导阻力减少约：$(1 - h/b) \times 50\% = 0.8 \times 50\% = 40\%$

这解释了为什么地效飞行器如此高效。

题4：点源强度 $m = 10 m^2/s$ 位于壁面上方 $h = 1m$ 处。求壁面上 $x = 2m$ 处的水平速度。

提示：使用镜像法，壁面下方有镜像源

答案

真实源在 $(0, 1)$，镜像源在 $(0, -1)$

点 $(2, 0)$ 到两源的距离：$r = \sqrt{4+1} = \sqrt{5}$

两源在该点产生的水平速度相加： $u = 2 \times \frac{m}{2\pi r^2} \times \frac{x}{r} = 2 \times \frac{10}{2\pi \times 5} \times \frac{2}{\sqrt{5}} = \frac{4}{\pi\sqrt{5}} ≈ 0.57 m/s$

挑战题（5-8题）

题5：设计一个兰金卵形体，长度 $4m$，最大直径 $1m$。确定源汇的位置和强度，使流线型外形最优。讨论如何评估"最优"。

提示：考虑压力恢复和分离点位置

答案

源在 $x = -1.5m$，汇在 $x = 1.5m$，强度 $m = 2\pi V_\infty$

"最优"标准：

  1. 压力梯度平缓，避免强逆压梯度
  2. 最大厚度在 $x/L ≈ 0.3$ 处
  3. 后部压力恢复渐进
  4. 形状阻力最小

可通过调整源汇位置微调外形。

题6：两个直径 $D = 0.5m$ 的平行圆柱，中心距 $S = 1.5m$，来流速度 $V = 15 m/s$。分析可能的流动模式和涡脱落特征。如果要避免共振，两圆柱的固有频率应如何设计？

提示：考虑间距比 $S/D$ 的影响

答案

$S/D = 1.5/0.5 = 3$，属于过渡区域。

可能的流动模式：

  • 同步脱落（同相或反相）
  • 准周期脱落
  • 间隙流偏转

涡脱落频率：$f ≈ 0.2 \times 15/0.5 = 6 Hz$

设计建议：

  • 两圆柱固有频率应避开 5-7 Hz 范围
  • 可设计为 $f_1 < 4 Hz$，$f_2 > 9 Hz$
  • 或使用不同直径打破对称性

题7：利用儒科夫斯基变换，设计一个升力系数 $C_L = 1.2$ 的对称翼型，弦长 $c = 2m$，来流速度 $V = 50 m/s$，迎角 $\alpha = 8°$。计算所需的环量和圆的参数。

提示：$C_L = 2\pi\sin(\alpha + \beta)$ 对于儒科夫斯基翼型

答案

对称翼型：$\beta = 0$

从 $C_L = 1.2 = 2\pi\sin(8°)$，验证：$2\pi \times 0.139 = 0.87$（需要更大迎角）

实际需要：$\sin\alpha = 1.2/(2\pi) = 0.191$，$\alpha ≈ 11°$

环量：$\Gamma = C_L \times \frac{1}{2}\rho V^2 c / (\rho V) = 0.5 \times 1.2 \times 50 \times 2 = 60 m^2/s$

圆参数：中心在原点，半径 $R = c/4 = 0.5m$

题8：开放性问题：某体育场屋顶设计成流线型，跨度 200m，高度 40m。风速 30m/s 时，如何用势流理论评估屋顶的升力？讨论该方法的局限性和改进方案。

提示：考虑三维效应、分离、真实边界层

答案

势流评估方法：

  1. 将屋顶简化为儒科夫斯基翼型或椭圆
  2. 用面元法计算压力分布
  3. 积分得到升力

局限性：

  • 忽略粘性分离（实际可能在 $x/c > 0.7$ 分离）
  • 无法预测湍流影响
  • 三维效应（端部效应、横向流动）
  • 非定常脱落（动态载荷）

改进方案：

  1. 势流 + 边界层修正
  2. 考虑分离泡模型
  3. 风洞试验验证
  4. 非定常势流 + 涡方法
  5. 最终需要全CFD验证

工程估算：$C_L ≈ 0.3-0.5$（考虑分离），升力 $≈ 500 kN$（需要仔细的结构设计）

6.10 常见陷阱与错误

理论应用的误区

  1. 过度相信势流预测 - ❌ 错误："势流说没有阻力，所以流线型物体阻力为零" - ✅ 正确：势流只能预测压差阻力的理想情况，真实阻力包括摩擦和分离

  2. 忽视雷诺数范围 - ❌ 错误：在 $Re = 100$ 时用势流理论 - ✅ 正确：势流适用于 $Re > 10^4$ 且无分离的情况

  3. 环量的任意性 - ❌ 错误："可以任意选择环量大小" - ✅ 正确：环量由库塔条件唯一确定

数值计算陷阱

  1. 奇点太近边界 - ❌ 错误：源点放在物体表面上 - ✅ 正确：源点应在物体内部，距表面至少 0.1 倍特征长度

  2. 镜像法的误用 - ❌ 错误：弯曲壁面直接用平面镜像 - ✅ 正确：只有平面和圆形边界有简单镜像解

  3. 保角变换的数值问题 - ❌ 错误：在变换奇点附近计算 - ✅ 正确：识别奇点位置，特殊处理或避开

工程应用错误

  1. 共振评估不当 - ❌ 错误："涡脱落频率与固有频率相差 10% 就安全" - ✅ 正确：±30% 范围内都可能发生锁频，需要仔细分析

  2. 地面效应的误判 - ❌ 错误："只要离地飞行就有显著地效" - ✅ 正确：$h/b > 0.5$ 时地效可忽略

  3. 忽视三维效应 - ❌ 错误：用二维理论直接计算有限展长机翼 - ✅ 正确：需要考虑翼尖涡和诱导阻力修正

6.11 最佳实践检查清单

概念设计阶段

  • [ ] 确认雷诺数范围适合势流分析（$Re > 10^4$）
  • [ ] 识别可能的分离区域
  • [ ] 选择合适的基本流动元素组合
  • [ ] 考虑三维效应的影响程度
  • [ ] 评估非定常效应的重要性

建模与计算

  • [ ] 叠加顺序：先远场（均匀流），后近场（源、涡等）
  • [ ] 边界条件检查：法向速度为零
  • [ ] 库塔条件：后缘速度有限
  • [ ] 网格独立性：面元数量足够（前后缘加密）
  • [ ] 奇点处理：避免数值奇异性

结果验证

  • [ ] 质量守恒：源汇强度平衡
  • [ ] 环量守恒：总环量检查
  • [ ] 压力系数范围：$C_p$ 在合理范围（通常 $-3 < C_p < 1$）
  • [ ] 远场行为：扰动衰减符合理论（$\sim 1/r$ 或 $1/r^2$）
  • [ ] 对称性检查：对称问题的解应对称

工程评估

  • [ ] 涡脱落频率计算（$St \approx 0.2$）
  • [ ] 共振风险评估（避开固有频率 ±30%）
  • [ ] 地面/壁面效应修正
  • [ ] 堵塞效应评估（风洞或管道）
  • [ ] 安全系数：考虑理论简化带来的不确定性

报告与交流

  • [ ] 明确说明假设条件和适用范围
  • [ ] 提供误差估计或不确定性范围
  • [ ] 建议需要的试验验证
  • [ ] 指出需要更详细分析的区域
  • [ ] 给出清晰的工程建议

记住：势流理论是工具，不是目的。合理使用它来获得物理洞察，但始终要意识到其局限性。最好的工程师知道何时使用简单理论，何时需要更复杂的方法。