企业官网建设流程全解析

简介：面向风机与轴流风机气动噪声仿真入门及进阶人员的视频教程，系统讲解基于大涡模拟（LES）与FW-H方程的风扇旋转机械气动噪声模拟完整流程。教程围绕FLUENT旋转机械模拟通用流程、气动噪声计算设置、CFD后处理以及噪声指向图展开，覆盖从网格处理、求解设置到结果解读的关键环节，适合具备一定CFD基础、希望掌握噪声评估方法的工程师和研究人员学习。压缩包共17个文件，包含7张演示截图、4个文本说明、4个HTML配套页面和2个Word文档，整体大小约4.31MB，内容紧凑、便于按需查阅。目前已有162人学习使用，资源以视频讲解配合图文笔记的形式呈现，可帮助读者快速建立气动噪声仿真总体思路并复现计算流程，节省自主摸索时间。 做风扇和轴流风机噪声仿真，很多朋友一上来就问“能不能直接用CFD把噪声算出来”。这话问得没毛病，但真在工程里跑起来，绝大多数人都绕不开大涡模拟加FW-H方程这条路线。我自己从用RANS稳态算流场，到被客户追着要频谱，再到现在能把一套LES+FWH流程跑完并出报告，中间踩过的坑比教程里写明白的要多得多。所以这篇文章我不打算复述软件帮助文档里的菜单操作，而是把一套可复用的通用流程拆开讲清楚——从为什么必须这么搭配，到网格怎么画、时间步怎么定、声源面和接收点怎么放，再到后处理里那些容易出错的细节。如果你正准备做风扇类旋转机械的气动噪声模拟，或者已经跑出过一份“看起来不太对”的频谱，这篇文章应该能帮你省点时间。

1. 为什么风扇噪声仿真绕不开“LES+FWH”这个组合

先说结论：直接气动声学计算在风扇这类工业设备上，工程成本高到不现实；而单纯用稳态RANS流场，又根本拿不到非定常压力脉动，也就没有噪声频谱可言。LES负责把流场里的非定常涡结构解析出来，FW-H则把这些涡产生的压力脉动积分成远场声压。两者各干各的活，拼起来就是当前旋转机械噪声仿真的主流工程方案。

1.1 直接算声波为什么在工程上不划算

理论上最“直接”的噪声模拟是让计算域同时包含流动声源和声传播路径，网格细化到能解析声波波长，然后用可压缩N-S方程一口气把声波也算出来。听着很完美，但实际算一下就知道问题在哪：流场里的涡尺度很小，而声波波长远大于涡尺度；为了同时解析两者，网格量会高到离谱，时间步又要被声波CFL条件限制得极小，最后计算资源根本跟不上。对轴流风机这种叶尖马赫数不高的低速设备，声源和传播路径分开处理更划算——流动计算只覆盖近场非定常区域，声传播靠积分公式在几乎没有额外网格成本的情况下完成。这就是声类比方法的核心思路。

1.2 FW-H声类比到底在算什么

FW-H方程是Ffowcs Williams和Hawkings在1969年从Lighthill方程推导出来的，本质上是把运动固体边界对声场的作用等价成分布在边界上的单极子源和偶极子源。对低速旋转机械来说，厚度噪声（叶片的体积排开产生的单极子项）和载荷噪声（叶片表面压力脉动产生的偶极子项）占主导，而湍流涡产生的四极子噪声在低马赫数下很弱，工程上经常忽略。所以你在后处理里看到的接收点声压，其实是通过一个包围声源的控制面或叶片表面上的瞬态压力信号积分出来的，不是直接在流场里“测”到的。这也是为什么FW-H设置里声源面的位置和完整性直接影响结果精度。

1.3 什么情况下这个组合会翻车

LES+FWH不是万能的。如果风扇叶尖马赫数很高，四极子噪声占比明显上升，忽略四极子项会低估总声级；如果壳体和管道壁面的反射、声模态效应显著，FW-H积分本身不天然包含壁面反射，需要额外处理；如果流场里有明显的强逆压梯度分离，LES对网格和数值格式的敏感性会放大误差。做之前先判断流动工况属于哪一类，别把方法用错场景。

2. 从几何到网格：旋转域划分与声源区网格策略

2.1 计算域搭建和旋转界面怎么处理

风扇噪声仿真的计算域基本沿用旋转机械CFD的经典结构：一个包含叶片轮毂和叶轮旋转区域的网格子域，一个包含进出口管道或自由空间的静止子域，两者之间通过旋转界面（sliding mesh，也叫瞬态转子-静子界面）衔接。滑移网格在旋转机械噪声仿真里几乎是必选项，因为只有它能真实传递非定常压力脉动；而混合平面或冻结转子界面只适合稳态流场，用在瞬态噪声计算上会直接切掉频谱里的关键成分。

计算域尺寸一定要给足。入口和出口如果离旋转域太近，远场反射的压力脉动会污染FW-H声源面信号。我自己的习惯是：入口至少延伸5倍叶轮直径，出口至少延伸8到10倍，侧向空间也要留够，否则计算域边界对声传播路径的“假约束”会让你拿到的频谱明显偏高或偏低。尤其是做自由场环境模拟时，出口段的长度不足会让尾迹涡在出口边界附近产生非物理反射。

2.2 LES网格的硬指标：y+、增长率与叶顶间隙加密

LES对近壁面网格的依赖比RANS苛刻得多。要捕捉叶片表面的压力脉动，网格的y+需要控制在1量级甚至更低，否则近壁区湍流结构会被壁面模型抹掉，噪声频谱的宽频基底会掉一截。同时，网格法向增长率建议不超过1.2，过渡太陡会人为耗散湍流涡，导致高频段能量损失。

叶片表面流向和展向网格尺寸按叶片弦长和展长的比例来规划，至少保证弦向有200到300个节点，展向在主要负荷区有50个以上节点。叶顶间隙区域特别关键——轴流风机宽频噪声很大一部分来自叶尖泄漏涡，间隙里至少要保证5到8层网格，否则叶尖涡直接被数值耗散掉，频谱上1kHz以上那一段会明显“塌陷”，跟实测相差好几个分贝。

2.3 入口边界条件和湍流量的隐藏影响

很多人做风扇噪声只关注旋转域和叶片近壁面网格，忽略入口和出口边界条件的影响。入口如果直接给一个均匀速度，没有设置合理湍流强度，自由来流湍流对叶片产生的来流湍流噪声就会被漏掉；反过来，给了一个明显偏高的湍流强度，又会叠加一层非物理的宽频噪声。比较妥当的做法是：入口设置速度入口或质量流量入口，并指定一个实际工况合理的湍流强度（轴流风机管路里通常1%到5%之间），湍流尺度也不要乱填，按入口水力直径的比例估算；出口用压力出口，回流的湍流参数也要给，否则计算后期出口回流会产生虚假的宽频噪声源。

3. 求解器设置：时间步长、亚格子模型与统计稳态判断

3.1 时间步长怎么定：旋转角度法与CFL法各取所需

瞬态LES的时间步长决定了两件事：一是每个时间步内流动物理过程能否被解析，二是最终噪声频谱能达到多高的有效频率。工程上常用“旋转角度法”来约束：一个时间步内转子的旋转角度建议控制在0.5°到1°之间，叶片通道内的流场变化才能被充分捕捉。角度定下来，时间步长就是转速的简单除法。

举一个示意算例：某轴流风机叶片数Z=7，转速n=3000rpm，换算成每秒50转，一个时间步旋转0.5°对应的时间步Δt = 1/(50×720) ≈ 2.78e-5秒。这个时间步对应的采样率是36kHz，按Nyquist频率算，能保证18kHz以内的高频信号不被混叠。如果目标只想分析到10kHz，这个步长依然合适；但如果要分析到20kHz以上，就把每步旋转角降到0.2°或更小。

另一方面还要检查CFL条件。LES对时间推进的稳定性要求比RANS高，尤其是近壁面网格很细的区域，CFL数如果超过1，数值振荡会被误当成噪声源进入FW-H信号。所以时间步长选择要把旋转角度法和CFL两者都过一遍，取更严的那个。

3.2 WALE还是Smagorinsky：旋转机械场景的取舍

亚格子模型选型直接影响LES流场质量。Fluent里常见的Smagorinsky-Lilly模型实现简单、鲁棒性好，但在近壁区需要人为阻尼函数修正，对旋转机械这种壁面边界层和叶顶涡同时存在的复杂流场，容易在壁面附近给出偏高的湍流粘性，把压力脉动抹平。相比之下，WALE模型不需要阻尼函数，对壁面附近湍流行为的还原更准，在旋转机械LES里用得更多。如果你用的求解器还支持结构化LES或SBES这类混合模型，对网格量的需求会友好一些，但纯LES+FWH的分析路径下，WALE是足够稳妥的默认选型。

3.3 “什么时候开始算噪声”：统计稳态判断不能省

这是最容易导致噪声频谱整体“虚高”的坑之一。启动阶段流场从初始条件向真实流场过渡，会产生剧烈的非物理压力脉动，如果一开始就打开FW-H信号记录，这段过渡都会被算进噪声里，频谱的低频段会被明显抬高。正确做法是：先让流场跑足够长的物理时间，等监测点压力、转矩或流量系数进入统计稳态后，再清零FW-H信号开始记录。

判断统计稳态没有绝对标准，我个人的经验是至少让转子转满20到30圈，再用监控变量做滑动平均，确认均值不再漂移。另外要注意：统计平均的时长不够会让FW-H信号里混入周期性成分，频谱上出现假的窄带尖峰，这个现象在低速风扇里很容易跟BPF的谐波搞混。

4. FW-H噪声计算设置：声源面、接收点与有效频带之间的配合

4.1 声源面放哪里才不会漏掉主要噪声源

FW-H声源面的选择是设置里最影响结果的一步。常用有两种策略：一是直接把叶片表面作为声源面，此时偶极子载荷项和厚度项被精确积分；二是用一个包住旋转域的闭合多孔面（porous surface），此时可以额外捕获叶顶涡、尾迹等叶片表面以外区域的噪声贡献。对轴流风机来说，叶顶间隙泄漏涡和叶片尾迹往往是宽频噪声的重要来源，单靠叶片表面会漏掉这部分，所以更推荐“叶片表面+闭合多孔面”组合，或者在旋转域外边界建一个包裹完整的积分面。

多孔面不能跟叶片表面重叠，也不要离旋转域太近，以免滑移界面上的数值非物理脉动被当成声源；但也不能放得太远，否则声源区外的湍流伪声会被积分进来。工程上我习惯把闭合积分面放在旋转域外边界再往外偏一两个网格层的位置，保证能包住主要湍流区，同时又不把下游远尾迹的大涡结构全部包进来。

4.2 接收点布置的工程惯例

接收点（observer point）代表实际关心的噪声测量位置。很多标准里测风机噪声会把麦克风放在离风扇中心1米处的多个方向，比如标准中常见的前方1米、后方1米、侧向0.5到1米。在FW-H设置里，你可以直接定义这些远场点，比如：

• 入口轴向上游1m处（P1）
• 出口轴向下游1m处（P2）
• 径向45°方向1m处（P3）
• 径向90°方向0.5m处（P4）

这些位置都要落在计算域网格之外，FW-H的远场积分公式才能应用——如果接收点放在网格内部，那算出来的是流场里的伪声压，不是传播声压。

4.3 有效频率范围：采样率与记录时长的联合约束

噪声频谱的两个基本约束都来自信号处理：频率上限由时间步长对应的采样率决定（Nyquist频率），频率分辨率由总记录时长决定（Δf=1/T）。还是拿上面示例算：Δt≈2.78e-5秒，采样率36kHz，可分析频率上限18kHz；如果记录0.5秒物理时间，频率分辨率2Hz，足够把BPF=350Hz及其谐波尖峰看清楚。如果要求1Hz分辨率，就得记录满1秒物理时间，这对LES来说计算量直接翻倍。所以设置记录时常一定要跟你的频率分辨率需求对着算，不要盲目拉长。

5. 后处理：从压力脉动到SPL频谱的完整链路

5.1 别把流场内压力脉动直接当声压

新手最容易犯的错误：直接在流场里取一个监测点，把它上的压力脉动做FFT，然后当噪声频谱来用。这个做法混淆了“伪声”和“辐射声”——流场内监测点感受到的压力脉动，很大一部分是随流动速度一起运动的流体动力学扰动，不是以声速向外传播的声波。只有在FW-H积分得到接收点声压之后，那才是可以跟麦克风测量对比的远场声压。所以后处理的第一步永远是确认数据来源：你手里是一段接收点的时域声压信号，而不是流场监测点的静压。

5.2 FFT参数、窗函数与频谱呈现方式

拿到接收点时域声压信号后，做FFT之前先处理两件事：去均值、加窗函数。去均值是去掉信号里的直流分量，否则零频处会有一个巨大的伪峰；窗函数则为了抑制频谱泄漏，汉宁窗是气动声学分析中最常用的选择。轴流风机噪声频谱里BPF（叶片通过频率）的尖峰特别尖锐，如果不加窗，这个尖峰的能量会泄漏到相邻频段，把谐波和附近宽频基底全部抬高。

频谱呈现方式要看你关心的指标：

• 窄带谱：直接按FFT频率分辨率显示SPL，每个频点一个值。适合观察BPF尖峰和主要谐波。
• 1/3倍频程谱：把每个倍频带内的能量积分，显示成每个频带的声压级。适合对照产品噪声限值标准，以及观察宽频噪声的整体分布。
• A计权或Z计权：如果需要评估人耳主观感受，对频谱做A计权后再积分得到总声压级，否则直接用线性计权。

5.3 用Python手动复算SPL的示例

虽然商用软件自带后处理绘图，但自己用脚本复算一遍，能帮你确认软件输出没被处理得“看不出问题”。下面这段逻辑可以放进Jupyter里直接用：

import numpy as np # 从FW-H输出加载接收点时域声压信号，单位Pa p = np.loadtxt("receiver_p.txt")[:, 1] fs = 36000.0 # 采样率，来自时间步长 # 去均值、加汉宁窗 p = p - np.mean(p) w = np.hanning(len(p)) p_win = p * w # FFT并转换到单边频谱，得到各频率分量的RMS值 N = len(p_win) Y = np.fft.rfft(p_win) freq = np.fft.rfftfreq(N, d=1/fs) p_rms = np.abs(Y) / (N/2) / np.sqrt(2) * np.sqrt(2 / np.sum(w**2) / N * N) # 考虑到汉宁窗的能量修正 p_rms = np.abs(Y) / (N/2) / np.sqrt(2) * np.sqrt(8/3) # 转为声压级SPL，参考声压20微帕 p_ref = 2e-5 SPL = 20 * np.log10(p_rms / p_ref + 1e-12) # 总声压级OASPL = 各频带能量求和 OASPL = 10 * np.log10(np.sum(10**(SPL/10))) print(f"OASPL: {OASPL:.1f} dB")

实际计算里要特别注意加窗后的能量修正系数，不同窗函数对应的能量修正不同，汉宁窗的修正因子是sqrt(8/3)。如果不做这个修正，谱级会被低估，总声级也偏小。上面代码里我写了两行计算方式，正确的取后面那行带窗函数修正的，前面那行是为了说明不加修正会得到什么结果。

6. 算完之后怎么判断结果可信：三个常见误差源

6.1 接收点与转子的距离够不够远

FW-H公式本身是一个远场近似，接收点距离声源不能太近。虽然软件不会限制你定义近距离接收点，但距离太近时，近场效应会显著，积分结果会偏离真实的辐射声场。做风扇噪声仿真时，接收点至少放在距离转子中心2到3倍叶轮直径以外，最好跟实际测试标准里的麦克风位置一致。对比仿真和测试时，如果两者位置对不上，光这个差异就能带来2到3dB的误差。

6.2 宽频噪声偏低：网格和声源面的典型症状

如果算出来的频谱在宽频段整体比实测低，先检查叶顶间隙网格。叶尖泄漏涡对宽频噪声贡献大，间隙处网格一旦不够，涡强度被数值耗散，高频宽频段会系统性低估。另一个常见问题是多孔声源面没有把叶顶间隙区域完整包住，导致这部分声源贡献没有被积分到。

6.3 从频谱形状读问题

对比仿真和实测时，频谱形状比总声压级更能说明问题。如果BPF尖峰位置对但幅值偏差大，往往要看载荷噪声的预测——叶片表面压力脉动幅值、相位和叶轮域的离散误差；如果尖峰频率明显偏移，先检查转速和叶片数是不是设置错了；如果宽频基底形状对了但整体平移，大概率是声源面选取或湍流耗散水平的问题。我给的排查建议是：先看频谱，用形状差异定位问题环节，而不是只盯着总声级这一个数。

气动噪声仿真的结果验证永远是个反复迭代的过程。我第一次跑完一套轴流风机的LES+FWH流程时，OASPL倒是跟实测很接近，但打开窄带谱一看，BPF附近的形状对不上，后来发现是统计稳态开启太早，启动阶段的非物理脉动把低频抬高了。后来把FW-H记录起点往后推了半圈，频谱形状一下就正常了。这类问题在教程里很少细讲，但实际项目里比调网格更磨人。希望这篇把流程和坑都在同一条线里讲清楚了，你操作的时候能少走点弯路。

本文还有配套的精品资源，点击获取