4、叶片噪声数值预测方法：计算气动声学（CWA）基础、声类比方法（Lighthill、FW-H方程）、基于CFD的噪声预测流程、大涡模拟（LES）与分离涡模拟（DES）在噪声预测中的应用
叶片噪声的数值预测，说白了就是「用电脑算声音」。你可能会问：声音不是靠耳朵听吗？怎么算？
嗯，这里的关键在于——我们不是直接算声音，而是先算流场，再从流场里「提取」出声源，最后传播到远场。这套思路，就是计算气动声学（CWA）的核心。
核心观点：噪声预测不是玄学，是物理。你算得准不准，取决于你对流场的理解有多深。
4.1 计算气动声学（CWA）基础
计算气动声学，英文叫 Computational Aeroacoustics，简称 CWA。它和普通 CFD 最大的区别是什么？
我个人习惯这样区分：CFD 算的是「力」和「速度」，CWA 算的是「压力脉动」——也就是声音。声音的本质，就是压力在介质中的波动。
CWA 有几个基本要求，我得提醒你注意：
- 数值耗散要低：声音的能量很小，如果算法本身耗散大，声音还没传出去就被「算没了」。
- 网格分辨率要高：一个声波周期至少需要 10-20 个网格点来分辨。我在项目中遇到过，有人用粗网格算噪声，结果出来的频谱全是数值噪声——白忙一场。
- 边界条件要处理好：声波遇到边界会反射，处理不好就是「假回声」。
小技巧：如果你刚开始做 CWA，建议先用简单的二维算例验证你的算法。我曾经在三维模型上折腾了两周，结果发现是边界条件设错了——教训啊。
4.2 声类比方法：Lighthill 与 FW-H 方程
声类比方法，是气动声学里最经典、最实用的工具。它的思路很巧妙：
我们不直接求解声波方程，而是把流场中的「声源项」提取出来，然后用一个波动方程去算声音的传播。你想想看，这相当于把「发声」和「传声」分开了，大大降低了计算难度。
4.2.1 Lighthill 声类比
Lighthill 在 1952 年提出了这个方程，可以说是气动声学的开山之作。它的形式是这样的：
∂²ρ'/∂t² - c₀² ∇²ρ' = ∂²T_ij/∂x_i∂x_j其中 T_ij 是 Lighthill 应力张量，包含了雷诺应力、粘性应力和非线性项。说白了，右边就是「声源」。
Lighthill 方程适用于自由射流噪声，比如喷气发动机的尾流噪声。但它的局限性也很明显——它不考虑固体边界的影响。你想想，叶片表面就在那里，声音会被反射、散射，Lighthill 方程处理不了。
4.2.2 FW-H 方程
Ffowcs Williams 和 Hawkings 在 1969 年把 Lighthill 方程推广到了包含运动固体边界的情况。这就是 FW-H 方程，也是目前叶片噪声预测最常用的方法。
FW-H 方程包含三项声源：
- 单极子源（厚度噪声）：由叶片运动排开流体产生。说白了，就是叶片「推」空气的声音。
- 偶极子源（载荷噪声）：由叶片表面的非定常力产生。这是叶片噪声的主要来源，尤其是旋转机械。
- 四极子源（非线性噪声）：由湍流中的非线性效应产生。高速流动时不可忽略。
重要提醒：在低速风扇或风机中，偶极子源占主导。但在高速涡轮中，四极子源可能成为主要噪声源。我在做航空发动机噪声预测时，就吃过这个亏——一开始只算了偶极子，结果和实验差了好几个分贝。
4.3 基于 CFD 的噪声预测流程
好了，理论讲完了，咱们来看看实际怎么做。基于 CFD 的噪声预测，我一般按以下流程走：
- 流场计算：先用 CFD 算出叶片周围的非定常流场。这一步最关键，流场不准，后面全是白搭。
- 声源提取：从流场中提取 FW-H 方程需要的声源数据——比如叶片表面的压力脉动、体积内的湍流应力等。
- 声传播计算：用 FW-H 方程或其它声传播模型，把声源传播到远场观测点。
- 频谱分析：对远场声压信号做 FFT，得到声压级频谱。
- 结果验证：和实验数据对比，看看算得准不准。
这个流程看起来简单，但每一步都有坑。我建议你重点关注第一步——流场计算。你想想看，如果流场本身就不准，那声源项能对吗？
避坑指南：我曾经在某个项目中，为了节省计算时间，用了稳态 RANS 算流场，然后用经验公式估算声源。结果呢？算出来的噪声频谱和实验完全对不上。后来改用非定常计算，才勉强吻合。所以——非定常流场是噪声预测的基础，别偷懒。
4.4 大涡模拟（LES）与分离涡模拟（DES）在噪声预测中的应用
说到非定常流场计算，就不得不提 LES 和 DES。这两种方法，是目前叶片噪声预测的主流工具。
4.4.1 大涡模拟（LES）
LES 的思路是：大尺度的湍流结构直接计算，小尺度的用模型近似。为什么这样分？因为大尺度结构携带了主要的能量，也是噪声的主要来源。小尺度结构虽然多，但能量低，对噪声贡献小。
LES 的优点很明显：
- 能捕捉到主要的湍流结构，声源计算更准确
- 适用于宽频噪声预测
但缺点也很突出：
- 计算量巨大，尤其是高雷诺数流动
- 近壁面网格要求极高，y+ 通常需要小于 1
我的经验：LES 适合做研究，但不适合做工程。我在做风机噪声优化时，用 LES 算一个工况就要一周，根本没法做参数化研究。后来改用 DES，效率提升了好几倍。
4.4.2 分离涡模拟（DES）
DES 是 RANS 和 LES 的混合方法。在近壁面区域用 RANS，在远离壁面的分离区用 LES。这样既保证了壁面附近的计算效率，又保留了分离区的湍流结构。
DES 在叶片噪声预测中的应用：
- 适合高雷诺数、有分离流动的叶片
- 计算量比 LES 小一个量级
- 宽频噪声预测精度可接受
但 DES 也有问题——RANS 和 LES 的过渡区域处理不好，会产生「灰色区域」，导致湍流粘度不连续，影响声源计算。
总结一下：如果你追求精度、不差计算资源，用 LES。如果你要做工程优化、需要快速迭代，用 DES。我个人习惯是：先用 DES 做参数扫描，找到最优方案后，再用 LES 做最终验证。
4.5 本章知识体系
下面这张图，是我整理的本章知识体系。你可以把它当作一个「导航图」，看看各个知识点之间的关系。
这张图把本章的内容串起来了。你可以看到，CWA 是基础，声类比方法和 CFD 是两大工具，LES 和 DES 是具体实现手段。它们之间的关系，不是替代，而是互补。
最后说一句：噪声预测这个领域，理论和实践差距很大。书上写的公式再漂亮，到了实际工程中，网格、边界条件、时间步长……每一个细节都可能让你翻车。我的建议是：多动手，多对比实验数据，慢慢积累自己的「感觉」。
公众号:蓝海资料掘金营，微信deep3321