1. 为什么旋转机械噪声仿真要选重叠网格做了这么多年CFD被问得最多的一个问题是风扇、风机这类旋转机械的噪声仿真到底怎么做才能又准又不折磨人说实话旋转机械气动噪声仿真说难不难说简单也不简单难点不在声学模型上而是卡在网格和运动处理上。今天这篇就把我踩过的坑和验证过的流程完整写出来重点讲清楚重叠网格这条技术路线。先明确一个概念气动噪声的本质是流场脉动压力引起的声辐射。旋转机械的噪声主要来自叶片周期性扫过时的压力脉动离散单音噪声加上湍流脉动产生的宽频噪声。要捕捉这些物理量流场计算必须足够精细这直接决定了网格和时间步长的下限。传统的做法是滑移网格Sliding Mesh把计算域分成旋转域和静止域交界面两侧网格节点不需要对齐通过插值传递流动信息。滑移网格在旋转机械稳态和瞬态流场计算里非常成熟但做噪声仿真时有个突出短板当几何外型比较复杂比如轴流风扇带导流罩、前缘锯齿、多级叶片或者需要处理叶片大幅度摆动、变形这类非定常运动时滑移网格的旋转域划分会非常别扭网格质量很难保证交界面插值误差也会直接污染声源的精度。这时候重叠网格Overset Mesh的优势就体现出来了。重叠网格的核心思路是把计算域拆成互相重叠但各自独立生成的网格块背景网格负责捕捉主流场物体表面的网格嵌套在背景网格里通过挖洞Hole Cutting和插值Interpolation建立网格间的信息传递。物体怎么动只需要修改嵌套网格的运动状态不用强制交界面节点匹配也不会因为网格畸变导致计算发散。我在多个轴流风机和离心风机的噪声仿真实例里对比过滑移网格和重叠网格结论是重叠网格在捕捉叶尖涡和叶片表面压力脉动方面更干净尤其是前缘和叶尖附近网格可以加密而不连累背景网格算噪声这种对压力脉动精度极其敏感的工况这个优势是决定性的。重叠网格并不只适用于气动噪声凡是涉及多体相对运动、大位移运动、运动边界穿越计算域的场景比如直升机旋翼、风力机尾流、船舶螺旋桨、汽车雨刮器绕流基本都是重叠网格的看家领域。对于旋转机械风扇风机的噪声仿真它解决问题的方式更直接旋转域和静止域独立画网格交界处重叠旋转运动由求解器自动处理。我把整个技术路线的核心逻辑总结成一句话噪声仿真的精度上限取决于声源区的网格与时间分辨率而重叠网格让你在声源区投入更多网格资源而不付出额外的网格生成代价。这是它被工业界越来越看重的原因。2. 整体方案设计从稳态流场到声学结果2.1 声学理论模型怎么选气动噪声仿真的数值方法大体分两类直接声学模拟CAA和混合声学方法。直接声学模拟要求流场和声场同时解析网格量和计算量超出大多数工程场景的承受范围这里不做讨论。工程上用得最多的是基于Lighthill声类比衍生出来的FW-H方程。FW-H方程的思路很直接把流场中固体壁面的脉动压力、湍流脉动等声源信息先算出来然后通过一个远场积分公式把声传播到监测点位置。它不需要在计算域里解析声波传播因此允许流场网格和声学网格采用不同尺度计算资源可以集中在声源区。旋转机械噪声仿真用FW-H是主流选择。Fluent里与FW-H配套的还有宽带噪声源模型Broadband Noise Source Model。这个模型不计算声传播只基于湍流量估算声源分布好处是便宜能在设计初期快速判断主要噪声源位置。但它给不了监测点声压级只能做定性参考。我的建议是先用宽带噪声模型快速定位叶片表面哪一段噪声贡献大再调整网格加密重点最后启动FW-H瞬态计算拿准确定量结果。这两者配合使用比一上来直接跑FW-H稳妥得多。2.2 三步走稳态流场到瞬态声源再到声学计算完整的旋转机械气动噪声仿真流程我习惯拆成三个阶段第一步稳态流场计算。用稳态RANS求解器湍流模型选k-omega SST或者标准k-epsilon把整个计算域的基本流场先算收敛。这一步的目的是给瞬态计算提供一个合理的初始流场避免LES或DES从零起步时经历漫长的瞬态过渡段。稳态收敛的判据我一般看两点残差降到1e-4以下同时进出口流量误差小于1%。第二步瞬态流场计算。把求解器切到LES或DES时间步长按叶片通过频率的需求设定跑足够多的物理时间让流场进入统计稳态。这一步产出的是随时间变化的流场数据也是声源信息的主要来源。瞬态段的计算时间通常占整个项目的70%以上急不得。第三步声学计算。在瞬态流场计算过程中同时打开FW-H声学模块设置监测点和声源面随时间步推进累积声源数据。瞬态计算结束后对FW-H积累的压力-时间历程做FFT变换得到监测点的声压频谱再按A计权或线性计权处理成需要的声学指标。这里特别注意FW-H的声源面必须覆盖所有主要噪声源通常包括叶片表面、轮毂表面和旋转域的外边界。我见过有人只在叶片表面设声源面结果是低频段声压级明显偏低原因就是旋转域运动引起的体积声源没被捕捉到。重叠网格方案下嵌套网格的外表面本身就是天然的运动声源面这又是一个加分项。2.3 一个可复制的参考对象轴流风扇噪声仿真拿一个典型的轴流风扇项目来举例。风扇直径D约0.2 m叶片数Z7设计转速3000 rpm叶片额定工况流量约0.05 m³/s。这个尺度下的主要噪声来源是叶片通过频率BPF Z×转频 7×50 350 Hz及其谐波。整个流程的物理时间和网格规模参考如下稳态计算域约500万网格瞬态阶段在叶片周围加密后到约1200万到1500万网格物理时间跑0.2秒风扇转10圈时间步长取对应1度转角的时间即约5.6e-5秒。这套配置在32核的机器上大约需要跑2到3天。这个量级和周期对大多数工程验证和论文课题来说是可接受的。3. 网格策略与重叠网格实操细节3.1 几何前处理和计算域划分做重叠网格之前先要把计算域划分清楚。我的做法是建立两个独立的计算域一个大的背景域静止域包含进出口流道和外围环境一个嵌套域旋转域包裹风扇叶片和轮毂。嵌套域的外边界距离叶片表面至少要留出0.3倍风扇直径的间距确保重叠区域的插值有足够的质量。几何处理有几个细节容易出错第一嵌套域外表面必须是封闭的不能有开口或重叠面。Fluent里检测嵌套域外表面是否有效主要看能不能自动识别为切割边界Overset Boundary识别不了十有八九是几何缝隙。第二背景域和嵌套域的网格在重叠区尺寸要尽量接近最好控制在两倍以内。如果背景网格太粗、嵌套网格太细插值误差会明显增加压力脉动会被抹平。这个原则直接关系到声源精度。第三嵌套域和背景域的单位要统一。听起来是废话但我真的见过有人模型一个用mm一个用m最后溢出插值找了一周的问题。3.2 背景网格和嵌套网格的尺寸设计背景网格的尺寸设计遵循一个原则声学关心的频段决定最大网格尺度。LES/DES要求网格能分辨涡的尺度而涡的尺度又对应着声源频率。工程上有个经验公式可以估算最大允许网格尺寸最大网格尺寸Δ U_min / (2 × f_max)其中U_min是关心区域内最小流速f_max是要分辨的最高频率。举个例子叶尖附近流速U30 m/s想分辨到3500 Hz十倍BPF以上那么Δ最大约为30/(2×3500)0.0043 m也就是4.3毫米。但要说明这只是估算起点实际需要在叶片表面附近再加密到上述值的1/5到1/10才能让叶尖涡和边界层分离被解析出来。嵌套网格的体网格我习惯用多面体网格Polyhedral。Fluent里多面体网格单元数平均比四面体少一半还多插值精度略高收敛性也更好。叶片表面第一层网格高度根据y来确定。对LES最好y接近1这样壁面模型不用额外修正对DES采用k-omega SSTy在1到5之间都可以接受。第一层网格高度的估算公式是Δy y × μ / (ρ × u_τ)其中u_τ是壁面摩擦速度。如果算不出精确的u_τ可以用经验值估算对叶尖速度30 m/s的工况第一层网格高度约在0.02到0.05毫米之间。这个数值不要盲信先算一遍看y分布再调。边界层内网格增长率控制在1.1到1.2之间沿叶高方向加密不少于40层。叶片尾缘附近的网格要单独加密因为尾缘脱落涡是旋转机械宽频噪声的重要来源。3.3 Fluent里重叠网格的创建与设置在Fluent里做重叠网格关键步骤我按操作顺序列一下第一步导入并排列网格。把背景网格和嵌套网格分别导入Fluent注意在Mesh Interfaces里把两个网格区分配成不同的Cell Zones。嵌套网格的旋转运动通过设定Cell Zone的运动类型来实现。第二步创建重叠网格界面。在Boundary Conditions面板里找到Overset Mesh选择嵌套域的所有外表面作为Overset Boundary。Fluent会自动执行切割Hole Cutting识别出哪些背景单元被覆盖哪些需要保留。这里有个小技巧切割结果可以用Overset Diagnostics来验证关注两个指标——被切割的单元数不能超过背景域单元总量的5%插值边缘Fringe上不能有未配对单元。第三步设定运动。旋转运动的定义有两种方式一是给定恒定转速刚性旋转适用于风扇风机这种匀速旋转二是通过UDF或者6DOF耦合计算运动轨迹适用于需要考虑气动反作用影响叶片转速的瞬态工况。风扇仿真绝大多数用恒定转速就够了。转速单位注意选对Fluent默认rad/s3000 rpm要填314.16。第四步设置插值格式。Fluent提供多种插值方法对旋转机械默认的Least Squares插值精度足够。如果发现重叠区出现压力振荡可以尝试更平滑的插值格式但计算量会涨不少。插值类型选三线性Trilinear还是反距离权重Inverse Distance Weighting我的实测结论是对同心旋转运动的旋转域Trilinear的守恒性更好推荐直接用。第五步监测重叠网格的稳定性。瞬态计算前100步是观察窗口如果残差震荡加剧、压力场出现非物理振荡优先检查插值边缘是否穿越了网格尺寸突变区域。常见的补救办法是把重叠区的背景网格局部加密而不是把嵌套域缩小。3.4 网格数量怎么估算才不浪费算力网格数量不是越多越好。工程上我通常按如下步骤反推先确定最高目标频率f_max。对风扇来说一般取十倍BPF也就是3500 Hz。再估算来流和叶尖相对速度的最大值U_max用上述公式算出最大允许背景网格尺寸。接着计算声源区包含的体积除以目标网格体积得到背景域大概需要的单元总量。最后加上嵌套域的加密网格数乘以一个1.2到1.5的余量系数就是总网格数上限。拿前面那个轴流风扇例子算背景域尺寸0.4 m × 0.4 m × 0.4 m目标背景网格4.3毫米理论单元量约80万但考虑到需要局部加密和过渡层实际背景网格做到约400万。嵌套域加密到叶片表面0.05毫米叶片附近单元数约800万。合计下来1200万左右和前面说的规模吻合。网格量的控制直接影响后续计算周期。500万网格转1200万网格计算时间通常是翻倍不止的因为时间步长也可能需要同步缩小。所以网格加密一定要锚定声学目标频率不要无脑加密。4. 求解器设置与参数选择4.1 稳态流场计算的设置要点稳态阶段采用密度基还是压力基取决于叶尖马赫数。风扇风机叶尖速度一般在30到80 m/s之间马赫数小于0.3属于典型不可压缩流动压力基求解器完全够用。只有转速特别高、叶尖进入跨声速范围才需要密度基。湍流模型稳态阶段用k-omega SST是稳妥选择。它对旋转机械近壁面的分离预测比k-epsilon准确而且和后续DES的切换是无缝的。如果建模对象主要是远场流道k-epsilon标准型也能胜任但一旦涉及叶尖涡捕捉SST明显更好。稳态计算收敛判据建议同时满足两条连续性残差和动量残差降低到初始值的1e-4以下风扇进出口质量流量差小于0.5%如果不满足第二条就强行进入瞬态很可能在瞬态启动阶段出现压力波传播的伪振荡噪声频谱开头全是杂波处理起来很麻烦。4.2 瞬态阶段LES还是DES瞬态流场质量直接决定FW-H声源的精度。在Fluent中旋转机械气动噪声的主流方案是近壁面采用RANS模式远离壁面的大分离区采用LES模式这就是DESDetached Eddy Simulation的核心思想。相比纯LESDES在同样的网格下能更好地处理高雷诺数壁面边界层而且不容易发散。如果计算资源充足、工况雷诺数不高也可以直接上LES。LES的亚格子模型在Fluent里最常用的是Smagorinsky-Lilly和WALE。我的实际经验是WALE更适合旋转机械它在近壁区域能自动恢复正确壁面渐近行为对壁面附近涡结构的捕捉更准在高剪切区域不会过度耗散。瞬态时间步长的确定是最容易出问题的一个环节。时间步长不够小高频噪声直接被滤掉太大残差和压力脉动不稳定。经验法则是时间步长对应旋转角度不大于1度。以3000 rpm为例每秒50转1度所需时间 1/(50×360) ≈ 5.56e-5秒即55.6微秒。如果要分辨到更高频率再按两倍甚至四倍加密时间步长。CFL条件也要检查Fluent里用压力基求解器时一般要求库朗数Courant Number小于1才能保证LES/DES的稳定性库朗数定义是U×Δt/Δx代入局部最大流速和最小网格尺寸核算看是否满足。内迭代次数我一般设20到30每时间步内保证残差至少降低一个数量级。如果达不到优先检查网格质量而不是盲目增加迭代次数。4.3 FW-H声学模块的配置FW-H在Fluent里是瞬态计算开始前就要打开的。配置有几个关键点监测点的位置要和实验测点相互呼应。风扇噪声实验标准里麦克风位置常在出风方向45度、1米处仿真里同样设置监测点频率分辨率才可比。监测点离声源太近会引入近场效应声压级不是真实的远场辐射。声源面选择。FW-H源面可以选固体表面也可以选渗透面Porosurface。对旋转机械我建议用叶片表面加轮毂表面作为主要源面同时把旋转域外表面作为渗透面一并纳入。渗透面能捕获四极子声源湍流脉动直接辐射的声在高风速高湍流工况下这类声源占比不低。FW-H的数据输出时间要和瞬态计算同步。Fluent在计算过程中按时间步累积声源的时域信号计算结束后可以对时域信号做FFT。需要提前设定好FFT的参数包括采样频率由时间步决定和采样时长由总物理时间决定。频谱分辨率等于采样时长的倒数如果想分辨1 Hz间隔至少需要1秒的物理时间。对旋转机械通常频率分辨率取BPF的1/10左右就够了比如目标分辨10 Hz需要0.1秒物理时间。4.4 压力-速度耦合和离散格式压力基求解器下压力-速度耦合算法选Coupled Scheme比SIMPLE稳得多尤其对重叠网格这种插值较多的工况。Coupled在每步内同步求解动量和压力修正对非定常计算的稳定性帮助明显。缺点是内存占用高一些但在现代工作站上不成问题。空间离散格式上压力项用Second Order动量项和湍流量项用Second Order Upwind。LES运行阶段动量项尽量用Bounded Central Differencing它能在一阶迎风和二阶中心格式间自适应切换保精度同时抑制非物理振荡。瞬态时间推进用Bounded Second Order Implicit比First Order Implicit在相同时间步下精度更高但要注意库朗数限制。梯度插值选Least Squares Cell Based对多面体网格兼容好。千万别用Green-Gauss Node Based它在重叠插值边缘会产生不必要的数值扩散。5. 常见问题与排查技巧实录5.1 重叠网格切割失败或插值不收敛现象Fluent在整个计算域初始化时报错提示某些单元没有接受插值数据或者残差在迭代初期就剧烈震荡。排查思路先在Overset Diagnostics里检查切割结果。常见原因有三个嵌套域外表面没有完全封闭背景网格在重叠区域太疏嵌套域网格嵌入后切割边缘出现孤立单元嵌套域初始位置和背景域重叠度太小低于10%会导致插值边缘质量差我的处置办法如果是封闭性问题回到网格生成软件里重构外表面如果是尺寸失配就把重叠区的背景网格局部加密。这事别拖切削失败越早发现代价越小。5.2 计算发散压力场崩溃或流速溢出现象算到某一步连续性残差突然飙升压力场出现正负交替的方波流速出现异常大值。最常见的原因一个是时间步长过大导致库朗数超过1尤其在叶尖高速区。另一个是稳态初始流场没收敛好带有大尺度分离涡瞬态启动后压力波剧烈响应。处理办法把时间步长缩小一半同时把动量方程和压力方程的松弛因子降到0.5左右先缓算500步等流场稳定后再恢复。如果还是发散回到稳态阶段把残差多降一个数量级或者直接换一个更高质量的初始流场比如先用k-epsilon算一遍再切到SST。5.3 噪声频谱出现大量伪峰现象FW-H的声压频谱在高频段出现密集的窄带峰值这些峰的位置不随网格变化物理上不合理。这种伪峰通常是声源面选择不善导致数值噪声被当成声源累积。排查时先把叶片表面源面单独保留去掉渗透面再算一遍对比频谱变化。如果伪峰消失说明是渗透面穿过了强湍流区域把伪压力脉动当成了声源。解决方法是把渗透面移动到一个湍流较弱的位置或者只保留固体表面源面。另一种常见情况是监测点离声源太近速度势的瞬态脉动直接污染了压力信号。把监测点移到更远位置伪峰一般会明显减小。5.4 声压级整体偏低或偏高这个现象多数时候和声源面设置及时间采样窗口有关。整体偏高通常是渗透面离叶片太近把湍流脉动当成声辐射整体偏低则常见于采样时长不够低频段分辨率不足导致FFT能量泄漏。给FFT加Hann窗可以减少泄漏但采样时长不够时根本解法还是延长物理时间确保至少覆盖10个目标最低周期。若只是低频段偏低重点检查是否遗漏了旋转域运动声源。旋转部件作为整体运动时会产生厚度声源和载荷声源忽略它们低频段会明显缺声音。5.5 常见问题速查表问题现象可能原因排查顺序初始残差不降卡在平台期稳态未收敛或初始流场含大分离检查稳态残差必要时重算稳态瞬态残差周期性跳动时间步长过大或内迭代不足缩小时间步增内迭代次数压力场出现棋盘分布网格质量差或插值格式不适合检查网格偏斜率替换插值格式叶尖区域出现负体积提示嵌套网格运动导致网格交叉检查嵌套域几何和运动速度缩小时间步FW-H报告零声压级声源面没有选中或未开启FW-H检查声源面的Zone归属和声学模块状态监测点时域信号缓慢漂移瞬态尚未进入统计稳态延长物理时间丢弃前段信号再做FFT高频段声压级明显过高网格太粗或湍流模型耗散不足加密叶尾和叶尖区域网格改用WALE计算内存溢出面插值和梯度存储过大降低梯度格式内存占用改用单元数据传递或换机器5.6 重叠网格监测和处理的经验总结几个实操中稳定提升成功率的小技巧第一瞬态开始前先冻结流场检查重叠网格插值边缘的数值通量是否平衡。Fluent里有Field Overlay功能可以可视化背景网格和嵌套网格的场变量如果插值边缘上压力等值线有明显错位说明网格尺寸匹配有问题先修这个再启动计算。第二输出中间文件。每500步保存一个data文件同时输出若干截面的瞬时压力云图。后面做噪声诊断时这些文件能帮你确认伪声源到底是来自叶尖涡破裂还是来自交界面插值。第三FFT前的数据截断窗口要统一。FW-H的时域信号在计算结束时会截断截断位置如果选在压力脉动幅值大的时刻频谱会出现明显旁瓣。用前一个周期的最后时刻对齐再做FFT旁瓣会小很多。写在实际操作之后重叠网格Fluent做旋转机械气动噪声仿真这条路我自己从最初的反复发散到现在基本一次跑通最大的体会是这个方法确实比滑移网格多了一组需要维护的网格拓扑关系前期准备成本高一些但换来的声源区网格自由度和瞬态计算稳定性非常值。尤其是风扇叶片带前缘锯齿、叶顶开槽这类需要精细捕捉局部流场细节的模型重叠网格几乎是唯一能同时保证网格质量、运动精度和计算效率的选择。最后分享一个小技巧当你刚开始搭这套流程时不要一上来就模拟全尺寸高转速工况。先用低速、叶片数少的小模型走通一遍完整流程把重叠网格设置、FW-H参数、监测点位置这些环节调试稳定再上真实转速和真实尺寸。这样排查问题时变量少定位快。我每次换新结构都先这么做虽然看起来多花一两天整体时间反而省不少。
旋转机械气动噪声仿真新思路：重叠网格在Fluent中的完整实践与参数指南
1. 为什么旋转机械噪声仿真要选重叠网格做了这么多年CFD被问得最多的一个问题是风扇、风机这类旋转机械的噪声仿真到底怎么做才能又准又不折磨人说实话旋转机械气动噪声仿真说难不难说简单也不简单难点不在声学模型上而是卡在网格和运动处理上。今天这篇就把我踩过的坑和验证过的流程完整写出来重点讲清楚重叠网格这条技术路线。先明确一个概念气动噪声的本质是流场脉动压力引起的声辐射。旋转机械的噪声主要来自叶片周期性扫过时的压力脉动离散单音噪声加上湍流脉动产生的宽频噪声。要捕捉这些物理量流场计算必须足够精细这直接决定了网格和时间步长的下限。传统的做法是滑移网格Sliding Mesh把计算域分成旋转域和静止域交界面两侧网格节点不需要对齐通过插值传递流动信息。滑移网格在旋转机械稳态和瞬态流场计算里非常成熟但做噪声仿真时有个突出短板当几何外型比较复杂比如轴流风扇带导流罩、前缘锯齿、多级叶片或者需要处理叶片大幅度摆动、变形这类非定常运动时滑移网格的旋转域划分会非常别扭网格质量很难保证交界面插值误差也会直接污染声源的精度。这时候重叠网格Overset Mesh的优势就体现出来了。重叠网格的核心思路是把计算域拆成互相重叠但各自独立生成的网格块背景网格负责捕捉主流场物体表面的网格嵌套在背景网格里通过挖洞Hole Cutting和插值Interpolation建立网格间的信息传递。物体怎么动只需要修改嵌套网格的运动状态不用强制交界面节点匹配也不会因为网格畸变导致计算发散。我在多个轴流风机和离心风机的噪声仿真实例里对比过滑移网格和重叠网格结论是重叠网格在捕捉叶尖涡和叶片表面压力脉动方面更干净尤其是前缘和叶尖附近网格可以加密而不连累背景网格算噪声这种对压力脉动精度极其敏感的工况这个优势是决定性的。重叠网格并不只适用于气动噪声凡是涉及多体相对运动、大位移运动、运动边界穿越计算域的场景比如直升机旋翼、风力机尾流、船舶螺旋桨、汽车雨刮器绕流基本都是重叠网格的看家领域。对于旋转机械风扇风机的噪声仿真它解决问题的方式更直接旋转域和静止域独立画网格交界处重叠旋转运动由求解器自动处理。我把整个技术路线的核心逻辑总结成一句话噪声仿真的精度上限取决于声源区的网格与时间分辨率而重叠网格让你在声源区投入更多网格资源而不付出额外的网格生成代价。这是它被工业界越来越看重的原因。2. 整体方案设计从稳态流场到声学结果2.1 声学理论模型怎么选气动噪声仿真的数值方法大体分两类直接声学模拟CAA和混合声学方法。直接声学模拟要求流场和声场同时解析网格量和计算量超出大多数工程场景的承受范围这里不做讨论。工程上用得最多的是基于Lighthill声类比衍生出来的FW-H方程。FW-H方程的思路很直接把流场中固体壁面的脉动压力、湍流脉动等声源信息先算出来然后通过一个远场积分公式把声传播到监测点位置。它不需要在计算域里解析声波传播因此允许流场网格和声学网格采用不同尺度计算资源可以集中在声源区。旋转机械噪声仿真用FW-H是主流选择。Fluent里与FW-H配套的还有宽带噪声源模型Broadband Noise Source Model。这个模型不计算声传播只基于湍流量估算声源分布好处是便宜能在设计初期快速判断主要噪声源位置。但它给不了监测点声压级只能做定性参考。我的建议是先用宽带噪声模型快速定位叶片表面哪一段噪声贡献大再调整网格加密重点最后启动FW-H瞬态计算拿准确定量结果。这两者配合使用比一上来直接跑FW-H稳妥得多。2.2 三步走稳态流场到瞬态声源再到声学计算完整的旋转机械气动噪声仿真流程我习惯拆成三个阶段第一步稳态流场计算。用稳态RANS求解器湍流模型选k-omega SST或者标准k-epsilon把整个计算域的基本流场先算收敛。这一步的目的是给瞬态计算提供一个合理的初始流场避免LES或DES从零起步时经历漫长的瞬态过渡段。稳态收敛的判据我一般看两点残差降到1e-4以下同时进出口流量误差小于1%。第二步瞬态流场计算。把求解器切到LES或DES时间步长按叶片通过频率的需求设定跑足够多的物理时间让流场进入统计稳态。这一步产出的是随时间变化的流场数据也是声源信息的主要来源。瞬态段的计算时间通常占整个项目的70%以上急不得。第三步声学计算。在瞬态流场计算过程中同时打开FW-H声学模块设置监测点和声源面随时间步推进累积声源数据。瞬态计算结束后对FW-H积累的压力-时间历程做FFT变换得到监测点的声压频谱再按A计权或线性计权处理成需要的声学指标。这里特别注意FW-H的声源面必须覆盖所有主要噪声源通常包括叶片表面、轮毂表面和旋转域的外边界。我见过有人只在叶片表面设声源面结果是低频段声压级明显偏低原因就是旋转域运动引起的体积声源没被捕捉到。重叠网格方案下嵌套网格的外表面本身就是天然的运动声源面这又是一个加分项。2.3 一个可复制的参考对象轴流风扇噪声仿真拿一个典型的轴流风扇项目来举例。风扇直径D约0.2 m叶片数Z7设计转速3000 rpm叶片额定工况流量约0.05 m³/s。这个尺度下的主要噪声来源是叶片通过频率BPF Z×转频 7×50 350 Hz及其谐波。整个流程的物理时间和网格规模参考如下稳态计算域约500万网格瞬态阶段在叶片周围加密后到约1200万到1500万网格物理时间跑0.2秒风扇转10圈时间步长取对应1度转角的时间即约5.6e-5秒。这套配置在32核的机器上大约需要跑2到3天。这个量级和周期对大多数工程验证和论文课题来说是可接受的。3. 网格策略与重叠网格实操细节3.1 几何前处理和计算域划分做重叠网格之前先要把计算域划分清楚。我的做法是建立两个独立的计算域一个大的背景域静止域包含进出口流道和外围环境一个嵌套域旋转域包裹风扇叶片和轮毂。嵌套域的外边界距离叶片表面至少要留出0.3倍风扇直径的间距确保重叠区域的插值有足够的质量。几何处理有几个细节容易出错第一嵌套域外表面必须是封闭的不能有开口或重叠面。Fluent里检测嵌套域外表面是否有效主要看能不能自动识别为切割边界Overset Boundary识别不了十有八九是几何缝隙。第二背景域和嵌套域的网格在重叠区尺寸要尽量接近最好控制在两倍以内。如果背景网格太粗、嵌套网格太细插值误差会明显增加压力脉动会被抹平。这个原则直接关系到声源精度。第三嵌套域和背景域的单位要统一。听起来是废话但我真的见过有人模型一个用mm一个用m最后溢出插值找了一周的问题。3.2 背景网格和嵌套网格的尺寸设计背景网格的尺寸设计遵循一个原则声学关心的频段决定最大网格尺度。LES/DES要求网格能分辨涡的尺度而涡的尺度又对应着声源频率。工程上有个经验公式可以估算最大允许网格尺寸最大网格尺寸Δ U_min / (2 × f_max)其中U_min是关心区域内最小流速f_max是要分辨的最高频率。举个例子叶尖附近流速U30 m/s想分辨到3500 Hz十倍BPF以上那么Δ最大约为30/(2×3500)0.0043 m也就是4.3毫米。但要说明这只是估算起点实际需要在叶片表面附近再加密到上述值的1/5到1/10才能让叶尖涡和边界层分离被解析出来。嵌套网格的体网格我习惯用多面体网格Polyhedral。Fluent里多面体网格单元数平均比四面体少一半还多插值精度略高收敛性也更好。叶片表面第一层网格高度根据y来确定。对LES最好y接近1这样壁面模型不用额外修正对DES采用k-omega SSTy在1到5之间都可以接受。第一层网格高度的估算公式是Δy y × μ / (ρ × u_τ)其中u_τ是壁面摩擦速度。如果算不出精确的u_τ可以用经验值估算对叶尖速度30 m/s的工况第一层网格高度约在0.02到0.05毫米之间。这个数值不要盲信先算一遍看y分布再调。边界层内网格增长率控制在1.1到1.2之间沿叶高方向加密不少于40层。叶片尾缘附近的网格要单独加密因为尾缘脱落涡是旋转机械宽频噪声的重要来源。3.3 Fluent里重叠网格的创建与设置在Fluent里做重叠网格关键步骤我按操作顺序列一下第一步导入并排列网格。把背景网格和嵌套网格分别导入Fluent注意在Mesh Interfaces里把两个网格区分配成不同的Cell Zones。嵌套网格的旋转运动通过设定Cell Zone的运动类型来实现。第二步创建重叠网格界面。在Boundary Conditions面板里找到Overset Mesh选择嵌套域的所有外表面作为Overset Boundary。Fluent会自动执行切割Hole Cutting识别出哪些背景单元被覆盖哪些需要保留。这里有个小技巧切割结果可以用Overset Diagnostics来验证关注两个指标——被切割的单元数不能超过背景域单元总量的5%插值边缘Fringe上不能有未配对单元。第三步设定运动。旋转运动的定义有两种方式一是给定恒定转速刚性旋转适用于风扇风机这种匀速旋转二是通过UDF或者6DOF耦合计算运动轨迹适用于需要考虑气动反作用影响叶片转速的瞬态工况。风扇仿真绝大多数用恒定转速就够了。转速单位注意选对Fluent默认rad/s3000 rpm要填314.16。第四步设置插值格式。Fluent提供多种插值方法对旋转机械默认的Least Squares插值精度足够。如果发现重叠区出现压力振荡可以尝试更平滑的插值格式但计算量会涨不少。插值类型选三线性Trilinear还是反距离权重Inverse Distance Weighting我的实测结论是对同心旋转运动的旋转域Trilinear的守恒性更好推荐直接用。第五步监测重叠网格的稳定性。瞬态计算前100步是观察窗口如果残差震荡加剧、压力场出现非物理振荡优先检查插值边缘是否穿越了网格尺寸突变区域。常见的补救办法是把重叠区的背景网格局部加密而不是把嵌套域缩小。3.4 网格数量怎么估算才不浪费算力网格数量不是越多越好。工程上我通常按如下步骤反推先确定最高目标频率f_max。对风扇来说一般取十倍BPF也就是3500 Hz。再估算来流和叶尖相对速度的最大值U_max用上述公式算出最大允许背景网格尺寸。接着计算声源区包含的体积除以目标网格体积得到背景域大概需要的单元总量。最后加上嵌套域的加密网格数乘以一个1.2到1.5的余量系数就是总网格数上限。拿前面那个轴流风扇例子算背景域尺寸0.4 m × 0.4 m × 0.4 m目标背景网格4.3毫米理论单元量约80万但考虑到需要局部加密和过渡层实际背景网格做到约400万。嵌套域加密到叶片表面0.05毫米叶片附近单元数约800万。合计下来1200万左右和前面说的规模吻合。网格量的控制直接影响后续计算周期。500万网格转1200万网格计算时间通常是翻倍不止的因为时间步长也可能需要同步缩小。所以网格加密一定要锚定声学目标频率不要无脑加密。4. 求解器设置与参数选择4.1 稳态流场计算的设置要点稳态阶段采用密度基还是压力基取决于叶尖马赫数。风扇风机叶尖速度一般在30到80 m/s之间马赫数小于0.3属于典型不可压缩流动压力基求解器完全够用。只有转速特别高、叶尖进入跨声速范围才需要密度基。湍流模型稳态阶段用k-omega SST是稳妥选择。它对旋转机械近壁面的分离预测比k-epsilon准确而且和后续DES的切换是无缝的。如果建模对象主要是远场流道k-epsilon标准型也能胜任但一旦涉及叶尖涡捕捉SST明显更好。稳态计算收敛判据建议同时满足两条连续性残差和动量残差降低到初始值的1e-4以下风扇进出口质量流量差小于0.5%如果不满足第二条就强行进入瞬态很可能在瞬态启动阶段出现压力波传播的伪振荡噪声频谱开头全是杂波处理起来很麻烦。4.2 瞬态阶段LES还是DES瞬态流场质量直接决定FW-H声源的精度。在Fluent中旋转机械气动噪声的主流方案是近壁面采用RANS模式远离壁面的大分离区采用LES模式这就是DESDetached Eddy Simulation的核心思想。相比纯LESDES在同样的网格下能更好地处理高雷诺数壁面边界层而且不容易发散。如果计算资源充足、工况雷诺数不高也可以直接上LES。LES的亚格子模型在Fluent里最常用的是Smagorinsky-Lilly和WALE。我的实际经验是WALE更适合旋转机械它在近壁区域能自动恢复正确壁面渐近行为对壁面附近涡结构的捕捉更准在高剪切区域不会过度耗散。瞬态时间步长的确定是最容易出问题的一个环节。时间步长不够小高频噪声直接被滤掉太大残差和压力脉动不稳定。经验法则是时间步长对应旋转角度不大于1度。以3000 rpm为例每秒50转1度所需时间 1/(50×360) ≈ 5.56e-5秒即55.6微秒。如果要分辨到更高频率再按两倍甚至四倍加密时间步长。CFL条件也要检查Fluent里用压力基求解器时一般要求库朗数Courant Number小于1才能保证LES/DES的稳定性库朗数定义是U×Δt/Δx代入局部最大流速和最小网格尺寸核算看是否满足。内迭代次数我一般设20到30每时间步内保证残差至少降低一个数量级。如果达不到优先检查网格质量而不是盲目增加迭代次数。4.3 FW-H声学模块的配置FW-H在Fluent里是瞬态计算开始前就要打开的。配置有几个关键点监测点的位置要和实验测点相互呼应。风扇噪声实验标准里麦克风位置常在出风方向45度、1米处仿真里同样设置监测点频率分辨率才可比。监测点离声源太近会引入近场效应声压级不是真实的远场辐射。声源面选择。FW-H源面可以选固体表面也可以选渗透面Porosurface。对旋转机械我建议用叶片表面加轮毂表面作为主要源面同时把旋转域外表面作为渗透面一并纳入。渗透面能捕获四极子声源湍流脉动直接辐射的声在高风速高湍流工况下这类声源占比不低。FW-H的数据输出时间要和瞬态计算同步。Fluent在计算过程中按时间步累积声源的时域信号计算结束后可以对时域信号做FFT。需要提前设定好FFT的参数包括采样频率由时间步决定和采样时长由总物理时间决定。频谱分辨率等于采样时长的倒数如果想分辨1 Hz间隔至少需要1秒的物理时间。对旋转机械通常频率分辨率取BPF的1/10左右就够了比如目标分辨10 Hz需要0.1秒物理时间。4.4 压力-速度耦合和离散格式压力基求解器下压力-速度耦合算法选Coupled Scheme比SIMPLE稳得多尤其对重叠网格这种插值较多的工况。Coupled在每步内同步求解动量和压力修正对非定常计算的稳定性帮助明显。缺点是内存占用高一些但在现代工作站上不成问题。空间离散格式上压力项用Second Order动量项和湍流量项用Second Order Upwind。LES运行阶段动量项尽量用Bounded Central Differencing它能在一阶迎风和二阶中心格式间自适应切换保精度同时抑制非物理振荡。瞬态时间推进用Bounded Second Order Implicit比First Order Implicit在相同时间步下精度更高但要注意库朗数限制。梯度插值选Least Squares Cell Based对多面体网格兼容好。千万别用Green-Gauss Node Based它在重叠插值边缘会产生不必要的数值扩散。5. 常见问题与排查技巧实录5.1 重叠网格切割失败或插值不收敛现象Fluent在整个计算域初始化时报错提示某些单元没有接受插值数据或者残差在迭代初期就剧烈震荡。排查思路先在Overset Diagnostics里检查切割结果。常见原因有三个嵌套域外表面没有完全封闭背景网格在重叠区域太疏嵌套域网格嵌入后切割边缘出现孤立单元嵌套域初始位置和背景域重叠度太小低于10%会导致插值边缘质量差我的处置办法如果是封闭性问题回到网格生成软件里重构外表面如果是尺寸失配就把重叠区的背景网格局部加密。这事别拖切削失败越早发现代价越小。5.2 计算发散压力场崩溃或流速溢出现象算到某一步连续性残差突然飙升压力场出现正负交替的方波流速出现异常大值。最常见的原因一个是时间步长过大导致库朗数超过1尤其在叶尖高速区。另一个是稳态初始流场没收敛好带有大尺度分离涡瞬态启动后压力波剧烈响应。处理办法把时间步长缩小一半同时把动量方程和压力方程的松弛因子降到0.5左右先缓算500步等流场稳定后再恢复。如果还是发散回到稳态阶段把残差多降一个数量级或者直接换一个更高质量的初始流场比如先用k-epsilon算一遍再切到SST。5.3 噪声频谱出现大量伪峰现象FW-H的声压频谱在高频段出现密集的窄带峰值这些峰的位置不随网格变化物理上不合理。这种伪峰通常是声源面选择不善导致数值噪声被当成声源累积。排查时先把叶片表面源面单独保留去掉渗透面再算一遍对比频谱变化。如果伪峰消失说明是渗透面穿过了强湍流区域把伪压力脉动当成了声源。解决方法是把渗透面移动到一个湍流较弱的位置或者只保留固体表面源面。另一种常见情况是监测点离声源太近速度势的瞬态脉动直接污染了压力信号。把监测点移到更远位置伪峰一般会明显减小。5.4 声压级整体偏低或偏高这个现象多数时候和声源面设置及时间采样窗口有关。整体偏高通常是渗透面离叶片太近把湍流脉动当成声辐射整体偏低则常见于采样时长不够低频段分辨率不足导致FFT能量泄漏。给FFT加Hann窗可以减少泄漏但采样时长不够时根本解法还是延长物理时间确保至少覆盖10个目标最低周期。若只是低频段偏低重点检查是否遗漏了旋转域运动声源。旋转部件作为整体运动时会产生厚度声源和载荷声源忽略它们低频段会明显缺声音。5.5 常见问题速查表问题现象可能原因排查顺序初始残差不降卡在平台期稳态未收敛或初始流场含大分离检查稳态残差必要时重算稳态瞬态残差周期性跳动时间步长过大或内迭代不足缩小时间步增内迭代次数压力场出现棋盘分布网格质量差或插值格式不适合检查网格偏斜率替换插值格式叶尖区域出现负体积提示嵌套网格运动导致网格交叉检查嵌套域几何和运动速度缩小时间步FW-H报告零声压级声源面没有选中或未开启FW-H检查声源面的Zone归属和声学模块状态监测点时域信号缓慢漂移瞬态尚未进入统计稳态延长物理时间丢弃前段信号再做FFT高频段声压级明显过高网格太粗或湍流模型耗散不足加密叶尾和叶尖区域网格改用WALE计算内存溢出面插值和梯度存储过大降低梯度格式内存占用改用单元数据传递或换机器5.6 重叠网格监测和处理的经验总结几个实操中稳定提升成功率的小技巧第一瞬态开始前先冻结流场检查重叠网格插值边缘的数值通量是否平衡。Fluent里有Field Overlay功能可以可视化背景网格和嵌套网格的场变量如果插值边缘上压力等值线有明显错位说明网格尺寸匹配有问题先修这个再启动计算。第二输出中间文件。每500步保存一个data文件同时输出若干截面的瞬时压力云图。后面做噪声诊断时这些文件能帮你确认伪声源到底是来自叶尖涡破裂还是来自交界面插值。第三FFT前的数据截断窗口要统一。FW-H的时域信号在计算结束时会截断截断位置如果选在压力脉动幅值大的时刻频谱会出现明显旁瓣。用前一个周期的最后时刻对齐再做FFT旁瓣会小很多。写在实际操作之后重叠网格Fluent做旋转机械气动噪声仿真这条路我自己从最初的反复发散到现在基本一次跑通最大的体会是这个方法确实比滑移网格多了一组需要维护的网格拓扑关系前期准备成本高一些但换来的声源区网格自由度和瞬态计算稳定性非常值。尤其是风扇叶片带前缘锯齿、叶顶开槽这类需要精细捕捉局部流场细节的模型重叠网格几乎是唯一能同时保证网格质量、运动精度和计算效率的选择。最后分享一个小技巧当你刚开始搭这套流程时不要一上来就模拟全尺寸高转速工况。先用低速、叶片数少的小模型走通一遍完整流程把重叠网格设置、FW-H参数、监测点位置这些环节调试稳定再上真实转速和真实尺寸。这样排查问题时变量少定位快。我每次换新结构都先这么做虽然看起来多花一两天整体时间反而省不少。
