噪声与振动CFD模拟.pdf

《噪声与振动CFD模拟.pdf》由会员分享，可在线阅读，更多相关《噪声与振动CFD模拟.pdf（14页珍藏版）》请在三一文库上搜索。

1、渐渐 307 版权所有，专做 CFD 电线气动噪声分析与研究电线气动噪声分析与研究 目录目录 一、引言一、引言 1.1 气动声学概况 1.2 气动声学的发展 二、二、噪声模型介绍噪声模型介绍 2.1 直接模拟模型（CAA 模型） 2.2 噪声比拟模型（FW-H 模型） 2.3 宽频噪声模型 三、电线气动噪声模型分析与计算三、电线气动噪声模型分析与计算 3.1 问题描述 3.2 几何建模与网格划分 3.3 模型计算 3.4 结果分析与讨论 四、总结与思考四、总结与思考 五、参考文献五、参考文献 渐渐 307 版权所有，专做 CFD 一、引言一、引言 1.1 气动声学概况气动声学概况 从史前以来，

2、 人们对由流体运动而产生的声音已十分熟悉。 在空气中的声音， 人们称之为气体动力声，在液体中的声音则称它为流体（水）动力声。风吹过树 枝所产生的啸声、管乐器和水壶的唿哨声以及人的口哨声等，都是普通的气体动 力声源的实例。尽管人们长期与这种不可见的声源相处在一起，可是直到第二次 世界大战结束时，人们对它的了解仍然很少。涡轮喷气式飞机和火箭，这两个世 界上最大的气体动力声源的出现促使一些优秀的科学家进行了一系列关于气动 声学的研究，并用了几年的时间初步建立起了气动声学领域的理论基础。气动声 学它研究的内容是流动与声的相互作用问题， 重点放在研究流动及其与物体作用 导致噪声产生的机理。所包含的基本物

3、理过程是流体中的波动、涡和声等运动之 间的相互作用。 具体问题的气动噪声是多种多样的， 如叶轮机械噪声、 喷气噪声、 螺旋桨噪声、管道噪声等。气动噪声不仅引起环境污染，还会造成结构的疲劳和 破坏。在军事上，预防并减少噪声显得尤为突出。所以，理解、预测并最终能控 制气动噪声具有非常重要的意义。由此可见，气动声学是介于气体力学和声学基 础这两门学科之间的一门交叉性和边缘性学科。 当今国内外从事气体动力学和声 学研究的工作者，都怀着很大的兴趣对气动声学进行研究，当然，研究的侧重点 各有不同。 气动声学在近二三十年才获得发展，其主要原因有两个：一是由于第二次世 界大战后，大型喷气式飞机问世所带来的巨大

4、噪声污染，以及近几十年来，随着 世界各国对环境保护工作的重视而不断要求降低诸如通风、 空调、 冷却器、 汽车、 高速列车等系统中出现的气动噪声污染所致。 二是气动声学发展必须建立在气体 力学和声学这两门学科本身发展的基础上， 两门学科都有一些尚未完全解决的一 些难题，如气体力学中湍流流动和非定常作用力的研究，以及声学中的湍流气流 中运动的固体壁面所产生的声辐射场的估计等。因此气动声学的研究，无论从理 论计算还是从实验研究上来说都是有相当难度的1。 1.2 气动声学的发展气动声学的发展 随着计算机能力的提高，数值模拟方法越来越多的被应用到柱体绕流噪声的 研究中。 对柱体噪声进行的研究工作根据噪声

5、预测的复杂和精确程度可分为三大 类： 第一类是利用不完全分析方法，它把动力声源项和声传播过程分开来进行计 算， 声传播过程的数值计算是基于源项包含不稳定流动变量的波动方程的积分公 式， 在各种各样的波动方程类型中， 应用最广的要属 Lighthill 的声类比波动方程。 随后 Curle 针对静止固体壁面边界进行了扩展，Ffowcs Williams 和 Hawkings 则 针对运动固体壁面边界进行了扩展。在这一类方法中，声源数据的获得来自于计 算流体动力学 （CFD） 技术的发展， Hardin和Lamkin以及Cox、 Brentlner和Rumsey 就利用此方法对圆柱绕流辐射噪声进行

6、了数值研究。与其它方法相比，利用这种 方法计算非常有效。然而，在利用此方法的大多数研究中，研究主要集中在预测 基于 Curle 声类比观察点的噪声大小而没有对噪声产生机制进行详细的研究2。 第二类方法是基于流动变量由“基流”分量和“残流”分量表示的混合数值方 法， 由此产生了分别代表流场和声场的控制方程。利用此方法的最大好处是其数 值运算方法最适合求解粘性流场和声场。Hardin 和 Pope 针对低速不可压流动情 渐渐 307 版权所有，专做 CFD 况运用此方法进行了扩展研究，他们利用从不可压 N-S 方程获得的整个不可压 粘性流场来预测声场。一旦水力密度常数被准确获得，声辐射噪声就能从因

7、变量 被分成水力项和声学项的可压缩 Euler 方程预测到。但后来 Shen,Sorensen 发现 Hardin 和 Pope 的提出的扩展公式中存在矛盾并提议修正予以克服。该修正公式 被应用到层流绕流单圆柱和湍流绕流机翼中。另一种数值混合方法是基于 Ewert 和 Schroder 的声扰动方程，这些方程从可压缩和不可压缩流动预测产生的声源 项驱动。他们运用此方法对层流绕流单柱体进行了研究。然而，获得的结果只是 定性的，没有提供声场及其产生机制的详细论述。第三类方法是利用非稳态可压 N-S 方程直接计算声场的计算声学（CAA）技术来预测气动噪声。最近在空气动 力噪声产生机制上的研究上的进展

8、主要是基于计算声学技术。然而，解决低马赫 数气动噪声问题需求庞大的计算资源，对其进行直接数值计算还是非常困难的。 大多数的研究工作局限在层流范围内低马赫数气动噪声问题。最近，Inoue 和 Hatakeyam 利用直接数值模拟方法对低马赫数下层流绕流单圆柱产生的气动噪 声辐射问题进行了二维数值计算。其结果与用 Curle 声学类比方法获得的结果比 较发现：Curle 方法的数值解，包含多普勒效应，不仅精确地描述了声音的产生 机制， 而且还描述了其传播过程；但他们的结果不能充分描述出气动噪声产生的 机制3。这部分是因为 DNS 在评估噪声源的相对贡献上存在内在的困难，因为 它不能辨别出单个的噪声

9、源。 二、噪声模型介绍二、噪声模型介绍 2.1 直接模拟模型（直接模拟模型（CAA 模型）模型） 在这种方法中声音的产生和传播直接通过求解合适的流体动力学方程获得。 声波的预测要求控制方程的时间精确解。进一步讲，在大多数直接模型的实际应 用中，必须借助于能够模拟粘滞效应和湍流效应的控制方程，例如非稳态 N-S 方程，雷诺时均方程以及过 DES 和 LES 使用的过滤方程。 直接模型需要高精度的求解方法，非常细密的计算网格以及声音无反射边界条 件，所以计算代价大。当预测远场噪声（几百倍的机翼弦长处得噪声）计算代价 更大。当计算近场噪声，直接方法就变的可行，如舱室噪音。对于许多近场噪声 的计算中，

10、由于局部压力波动导致的噪声是可以通过 fluent 准确计算的4。 2.2 噪声比拟模型（噪声比拟模型（FW-H 模型）模型） 对于中场和近场噪声，fluent 采用基于 Ligthill 的声比拟方法，它是直接模拟 的一个很好的补充。在该方法中，近场流场从控制方程中获得，如非稳态的雷诺 平均方程，过滤的 DES 和 LES 方程，然后把求解结果作为噪声源，通过求解波 动方程得到解析解，这样就把流动求解过程从声学分析中分离出来。 Ansys Fluent 采用基于 FW-H 的方程，FW-H 方程采用最通用的 lighthill 的噪 声比拟方法，可以求解单极子、偶极子和四极子产生的噪声传播。

11、Ansys Fluent 采用时域积分的方法（声压、噪声信号与时间相关） ，通过面积分计算指定位置 的噪声。流场变量，如压力，速度时间精度的解的获得需要求解面积分。时间精 度的解可以从非稳态雷诺平均方程，大涡模拟，或分离涡模型求解获得，可以捕 捉精确的流动特征如涡脱落等现象。Ansys Fluent 中的噪声积分源面不仅可以放 在不可穿透壁面上，也能放在内部可穿透面上，这样就可以考虑源面包围的四极 子噪声的贡献。Ansys fluent 中先进的基于“可穿透面积分”的 FW-H 模型可以很 渐渐 307 版权所有，专做 CFD 好解决高、亚速气动噪声问题。 Ansys Fluent 中的 FW

12、-H 噪声模型可以选择多个源面和接收位置，也可以保 存噪声源数据，或在瞬态流动求解过程中同时执行噪声计算。声压信号可以快速 傅里叶变换得到，结合 Ansys Fluent 后处理获得全部声压标准（SPL）和能量谱 范围下的噪声数据图。 声比拟模型基于两步法：首先采用 CFD 方法在噪声源附近精确的计算瞬态 流场，其次是从声源处到接收处噪声传播通过求解波动方程获得4。 2.3 宽频噪声模型宽频噪声模型 在许多工程应用中的湍流，噪声没有明显的频段，声波能量连续分布在一个 宽频段范围内按频率连续分布，这就涉及到宽频噪声问题。在 Ansys Fluent 的宽 频噪声模型中， 湍流参数通过雷诺时均方程

13、求出， 再用一定的半经验修正模型 （如 Proudman 方程模型、边界层噪声源模型、线性 Euler 方程源项模型、Lilley 方程 源项模型）计算表面单元或体积单元的噪声功率。 三、电线气动噪声模型分析与计算三、电线气动噪声模型分析与计算 3.1 问题描述问题描述 如图所示为电线气动噪声的几何模型的示意图，电线地的直径为 1.9cm，计 算流体区域宽 9cm，长 38cm。气流以 69.2m/s 的速度从左端入口流入，由右端 出口流出。圆柱为刚性材料。不考虑变形，本次模拟采用 LES 模型进行模拟计 算分析。同时启用噪声模型模拟气动噪声。 图 3-1 电线气动噪声模型的示意图 3.2 几

14、何建模与网格划分几何建模与网格划分 首先在 Gambit 中对模拟模型进行建模，采用 Tri 类型网格，划分网格如下： D=1.9cm Inlet Outlet 渐渐 307 版权所有，专做 CFD 图 3-2 模型网格划分图 3.3 模型计算模型计算 边界条件：计算模型中采用速度入口 Velocity inlet，压力出口 Pressure outlet 边界采用对称边界 symetry。 Fluent 的求解设置，打开 Fluent16.0 选 2D 求解器，其余保持默认，然后导入 划分好的模型网格格式为 msh 的文件。在 General 面板中单击 Scale，对模型的 尺寸进行缩放。

15、 在 General 面板中将 Solver 中的 time 设置为 Transient（非定常） ，其余保持默认。 选择求解模型 Large Eddy Simulation（LES）模型： 渐渐 307 版权所有，专做 CFD 流体介质选择 air 空气。操作条件 Opreation Conditions 中设置（Operating Pressure）为 101325Pa。 在边界条件 Boundary Conditions 中速度进口设置气流速度为 62.9m/s,出口为压力 出口，压力为 0。具体设置如下图所示： 在求解方法中 Solution Methods 选项中，Pressure-

16、Velocity Coupling 的 Scheme 下拉列表 中选择 PISO 选项，在 Spatial Discretization 的 Pressure 中选择 PRESTO！其余保持默认。在 渐渐 307 版权所有，专做 CFD 求解控制中 Solution Controls 设置 Pressure 为 0.75 其余保持不变。具体如下： 在残差监控窗口中设置收敛标准。选择 Drag 和 Lift 两个监视数值。 运行计算，在 Run Calculation 选项中，选择迭代次数 4000 步，进行计算。 3.4 结果分析与讨论结果分析与讨论 流场计算完成后，选择 Results 中

17、Plots 选项，双击 Histogram 选项，选择 Velocity 和 Cell Courant Number 选项。得到如下图所示： 图 3-2 Histogram of Cell Courant Number 曲线 从图中可以看出Cell Courant Number 集中在 1 以下。 流场计算完成之后，进行噪声计算，选择模型 Models 中 Acoustics-off，选中 渐渐 307 版权所有，专做 CFD Ffowcs-Williams&Hawkings，在 Export Option 中选择 Export Acoustic Source Data in ASD Form

18、at选项。 然后继续在Acoutics Model对话框中， 将Far field Sound Speed 修改为 340，Reference Acoustic Pressure 修改为 2e-05，Source Correlation Length 修改为 0.095，如下图所示： 计算之后，分别对电线附近的速度云图、压力云图、以及湍流图进行显示。 图 3-3 速度云图 渐渐 307 版权所有，专做 CFD 图 3-4 压力云图 图 3-5 湍流程度云图 渐渐 307 版权所有，专做 CFD 图 3-6 Cell Courant Number 从图中的直方图可以看出，Cell Courant

19、 Number 集中在 1 以下。 分别对两个不同的采集点进行观察，得到声压曲线如下： 图 3-7 在位置 1 的声压曲线图 渐渐 307 版权所有，专做 CFD 图 3-8 在位置 1 考察局部的声压曲线图 图 3-9 在位置 2 的声压曲线图 渐渐 307 版权所有，专做 CFD 图 3-10 在位置 2 考察局部的声压曲线图 图 3-11 在位置 1 的 Struhal Number 曲线 渐渐 307 版权所有，专做 CFD 图 3-12 在位置 2 的 Struhal Number 曲线 四、总结与思考四、总结与思考 通过对于电线的气动噪声进行了模拟，熟悉了噪声模型的种类，以及噪声模

20、 型的使用范围，相关的噪声计算参数的设置，对于现实的实际问题有了更加深刻 的理解。通过对电线绕流的气动噪声数值模拟，验证了该方法可以准确地模拟层 流和湍流近场的声场分布， 且该方法可以和计算流场精度较高的亚格子湍流模型 相结合，更准确地计算近场气动噪声的产生机理和传播过程。 对于后期的模型可以进一步的优化，以确保能够更加准确的接近真实模拟体 系，更加有助于解决实际的问题。 渐渐 307 版权所有，专做 CFD 五五、参考文献、参考文献 1马大猷.现代声学理论基础.科学出版社，2005. 2倪大明，张文平，明平剑. 低速湍流中气动噪声预测方法J,2013. 3杨南奇,串列双圆柱绕流噪声数值模拟与实验研究,学位论文，2009. 4朱红均，Fluent 流场分析与实战指南M,人民邮电出版社，北京 2014.