第二章 流体力学基础回顾：连续性方程、伯努利方程、动量定理与动量矩定理
各位好，欢迎来到叶片空气动力学核心解析的第二讲。
说实话，很多搞风电或航空的同行，一提到流体力学就头疼。我当年刚入行时也一样，觉得那些方程离实际工程太远了。但干久了你会发现——叶片设计的灵魂，就藏在这几个基础方程里。今天咱们就把它们掰开揉碎了讲清楚。
2.1 连续性方程：质量守恒的直观表达
先问个问题：为什么水管捏扁了，水流反而变急？
答案就在连续性方程里。说白了就是一句话：流进多少，就得流出多少。
对于一维定常流动，数学形式很简单：
ρ₁ · A₁ · V₁ = ρ₂ · A₂ · V₂其中 ρ 是密度，A 是截面积，V 是速度。
如果流体不可压缩（比如低速下的空气或水），密度不变，公式就简化为：
A₁ · V₁ = A₂ · V₂这意味着：截面积越小，流速越大。
核心理解：连续性方程不是算出来的，是“守恒”出来的。质量既不会凭空产生，也不会凭空消失。
我的经验：我在做风力机叶片设计时，经常用这个方程估算叶片不同截面处的来流速度变化。尤其是在叶片根部到尖端的过渡区域，面积变化剧烈，速度分布必须用连续性方程先过一遍，否则后面算载荷全是错的。
2.2 伯努利方程：能量守恒的工程利器
伯努利方程，可能是流体力学里最出名的一个方程了。飞机为什么能飞起来？叶片为什么能产生升力？都跟它有关。
对于不可压缩、无粘性、定常流动，沿流线有：
p + ½ρV² + ρgz = 常数三项分别代表：压力能、动能、势能。它们之和沿流线保持不变。
实际工程中，忽略重力项（气体），我们常用的是：
p₁ + ½ρV₁² = p₂ + ½ρV₂²流速快的地方，压力低；流速慢的地方，压力高。这就是升力产生的根本原因。
注意：伯努利方程只适用于无粘、不可压缩、定常流动。在叶片表面边界层内部，粘性力占主导，伯努利方程就不成立了。我曾经见过有新手直接用伯努利方程算叶片表面压力分布，结果跟实验数据差了30%以上——就是因为没考虑粘性。
避坑指南：我曾经在调试一个小型风力机时，发现叶片表面压力分布总是不对。后来仔细排查，原来是叶片前缘有微小污染，导致边界层提前转捩。伯努利方程在层流区还好用，一到湍流区就失效了。所以用之前，先确认流动状态。
2.3 动量定理：力从何来
动量定理讲的是力与动量变化率的关系。对于控制体，有：
ΣF = ∂/∂t ∫∫∫ ρV dV + ∫∫ ρV (V·n) dA看着复杂？其实核心就一句话：流体动量变化，等于作用在它上面的合力。
在叶片分析中，我们常用的是简化形式：
F = ṁ · (V₂ - V₁)其中 ṁ 是质量流量，V₁ 和 V₂ 分别是进出口速度。
这个公式在风力机致动盘理论中极其重要。通过动量定理，我们可以直接算出风轮从风中提取了多少能量。
工程应用：用动量定理计算风力机推力时，不需要知道叶片具体形状，只需要知道来流速度和尾流速度。这就是为什么Betz极限（59.3%）可以用纯动量理论推导出来——跟叶片形状无关。
2.4 动量矩定理：扭矩的来源
动量矩定理，说白了就是动量定理的“旋转版本”。它描述的是：角动量的变化率等于外力矩。
数学形式：
ΣM = ∂/∂t ∫∫∫ (r × ρV) dV + ∫∫ (r × ρV) (V·n) dA对于叶轮机械，简化后常用：
M = ṁ · (r₂Vθ₂ - r₁Vθ₁)其中 Vθ 是切向速度分量。
这个方程直接给出了叶片扭矩的计算方法。你想想看，风力机叶片为什么是扭转的？就是因为不同半径处需要不同的切向速度变化，才能产生均匀的扭矩输出。
我的习惯：在做叶片气动设计时，我通常先用动量矩定理估算出需要的扭矩分布，然后再反推叶片的弦长和扭角分布。这样设计出来的叶片，气动效率往往更高。记住：扭矩是设计目标，几何是设计手段。
2.5 四个方程的内在联系
这四个方程不是孤立的。它们分别从不同角度描述了同一个物理过程：
| 方程 | 守恒量 | 工程应用 |
|---|---|---|
| 连续性方程 | 质量 | 速度-面积关系 |
| 伯努利方程 | 机械能 | 压力-速度关系 |
| 动量定理 | 动量 | 推力/升力计算 |
| 动量矩定理 | 角动量 | 扭矩计算 |
在叶片分析中，我们通常需要联合使用这些方程。比如用连续性方程确定速度，用伯努利方程求压力，再用动量定理算力，最后用动量矩定理算扭矩。一环扣一环，缺一不可。
总结一下：这四个方程是叶片空气动力学的“四根柱子”。你把这四个方程吃透了，后面学叶素动量理论、涡流理论都会轻松很多。我建议你花时间把每个方程的物理意义和适用条件搞清楚，而不是死记公式。
好了，这一章的内容就到这里。记住：基础不牢，地动山摇。这四个方程你掌握得越扎实，后面学叶片设计就越轻松。有什么问题，欢迎随时交流。
公众号:蓝海资料掘金营，微信deep3321