某型空调轴流风扇的气动噪声仿真分析

1、某型空调轴流风扇的气动噪声仿真分析 詹福良 1 游 斌 2 (1 LMS 国际公司北京代表处 , 2 美的制冷家电集团技术研发中心 1、概述 研究对象是带短导管轴流风扇的气动噪声分析问题。这里主要介绍使用 LMS 公司著名声学软件 SYSNOISE 的流体声学模块生成气动噪声声源 ,然后使用 SYSNOISE 强大的声学边界元 (Acoustic BEM 功能进行整个声场谐波分析的过程和 结果。其中带短导管轴流风扇的流场分析使用 FLUENT 软件 ,分析模型和输入数据 由美的研发中心的游斌博士提供。 SYSNOISE 模型的网格在 FLUENT 模型网格的 基础上快速生成得到。 本文目的在于

2、 SYSNOISE 流体声学功能演示和抛砖引玉 ,并未对计算结果的准 确程度特别关注。实际上本题计算结果的准确性由 FLUENT 流场计算的精度和 SYSNOISE 声学计算精度共同决定 ,我们这里只选取了 FLUENT 计算初期的部分流 场结果 ,初始条件扰动较大 ,导致噪声计算的结果可能偏大。共计算了 62300Hz 之 间的噪声分布 ,这里只列出部分结果。 根据与游斌博士的交流 ,做了两种不同网格密度的 SYSNOISE 模型进行验证。 一种是严格按照 CFD 导出网格的密度 ,直接生成对应的 SYSNOISE 模型;另一种是 按照声学分析理论进行了网格稀疏化的 SYSNOISE 模型。

3、对两种模型的内声场计 算结果进行了对比 ,结果基本完全一致。但是第二种稀疏化模型 (有效频率已经达到 8000Hz 的计算速度大幅度增加。实 际应用中建议使用第二种模型。 2、分析流程 C F D c o m p u t a t i o n i n t i m e d o m a i n (e .g . L E S L a r g e E d d y S i m u l a t i o n ; D E S 图 1SYSNOISE 的流体声学分析流程图 具体的分析流程如下 : A在 FLUENT模型网格基础上快速生成各种密度的 SYSNOISE模型网格。 B 使用 FLUENT 软件对带短导管轴

4、流风扇的流场进行非稳态分析 , 并在时域内 输出流 场的分析结果。 (FLUENT 分析模型和输入数据由美的研发中心游斌博士提 C使用 SYSNOISE 的流体声学模块直接读入 FLUENT 的流场分析结果 ,并生成 相应的 气动噪声声源 ,这里主要是壁面流体压力脉动产生的二极子声源。 D 使用 SYSNOISE 强大的 ( 耦合声学边界元 (Acoustic BEM 功能进行整个声场的 谐波分析计算和结果后处理输出。本文使用内外 DBEM 模型在短导管的端部 进行耦合来模拟开放导管情况 ,内模型由导管壁面与风扇之间构成的空腔组成 (法线 方向指向空腔内侧 ,外模型由导管壁面组成 (法线方向指

5、向空腔外侧。 图 2 内外声场耦合的边界元模型及不同的网格密度对比 3、结果分析 3.1 二极子声源提取结果 下面列出的是 SYSNOISE 进行二极子声源提取的结果 ,为了对比说明 SYSNOISE 的计算流程中对 CFD 网格进行疏化处理 ,然后得到适当密度的声学网格 的适用性 ,这里包含了两种网格密度模型的对比 :（左列稀疏网格 ;右列CFD 网格, 因为网格太密影响云图显示 ,所以特别将网格隐藏不显示 : i_zim ENGINEERING INNOVATION SYS*：SC Ker b 6 Opt in BEI Vat aaf aoaalC B 区! FUe Gomrtry Mod

6、el Analysis Inqutrr Tools View Postprocess DisplayHelp FF 产冋亍上厂厂厂J aFFKPTT7 CKKMBC 001 6ZDC8 3 IMS-01 im ENGINEERING INNOVATION SYSB01SK Ker b. B Opt ion Bkl V.“ FHr Gromrlry Model Analysis Inquorv Tools Viuw Postprocess Display ?IMOOD IWIMKD i2：SC Ker b 6 Opt in BEI Vat aaf aoaalC B 区! FUe Gomrtry

7、Model Analysis Inqutrr Tools View Postprocess DisplayHelp FF 产冋亍上厂厂厂J aFFKPTT7 CKKMBC 001 6ZDC8 3 IMS-01 im ENGINEERING INNOVATION SYSB01SK Ker b. B Opt ion Bkl V.“ FHr Gromrlry Model Analysis Inquorv Tools Viuw Postprocess Display ?IMOOD IWIMKD i2IS trr A6 叶BU 7 fit Cewctry Uo4ri Amh frt xh Vlrr Pw

8、ttt Dtiplay 冋示wf rrrr卜卜i川s卜it ENGINEERING INNOVATION -丸 IrrBU Ttfaatamal fit GetMtry LMelbfr xh ViewDrtpliy X8 imte ISMWKv i_zim ENGINEERING INNOVATION * 5T5IS trr A6 叶BU 7 fit Cewctry Uo4ri Amh frt xh Vlrr Pwttt Dtiplay 冋示wf rrrr卜卜i川s卜it ENGINEERING INNOVATION -丸 IrrBU Ttfaatamal fit GetMtry LMelbf

9、r xh ViewDrtpliy X8 imte ISMWKv 图 3 封闭内声场模型的计算结果及不同的网格密度对比 从结果图形可以直观 看出 ,两种不同密度的网格所得到的内部声场特征也基本完全一致 ,无论从分布趋势 和数值大小 ,都表现出很好的一致性。但是稀疏化模型 （有效频率已经达到 8000Hz 的计算速度大幅度增加 ,所以实际应用中一般都要使用稀疏化的声学网格模型。 某型空调轴流风扇的气动噪声仿真分析 2006 LMS 首届用户大会论文集 3.3 真 实声场模型（内外耦合）的计算结果 最后， 利用稀疏化的网格模型建立了真实的 轴流风扇内外声场耦合模型， 进行完整的轴流 风扇噪声辐射模拟

10、分析。下面列出 的是 SYSNOISE 进行内外部声场分布计算的结果（左列 内声场；右列 外声场 6 某型空调轴流风扇的气动噪声仿真分析 2006 LMS 首届用户大会论文集 7 某型空调轴流风扇的气动噪声仿真分析 2006 LMS 首届用户大会论文集 图 4 真实内外声场耦合模型的计算结果 4、结论与建议 从计算结果看到， 气动噪声的 声源主要来自风扇迎风面的中上部、 以及对应的管路壁面部 位。内部声场的气动 噪声主要分布在出风侧，幅值较高。而外部声场的气动噪声主要分布在 风扇平面 内，而不在风扇的流场方向上。本结论与航空领域的螺旋桨平面噪声现象比较一 致。 根据上面的分析对比过程， SYSNOISE 可以非常方便地解决这类流体声学分 析问题，高 效准确地得到气动声学的内外声场分布。关键的是， SYSYNO