非均匀进流条件下螺旋桨无空化噪声的预报与校验.docVIP

———————————————————————————————————————————————非均匀进流条件下螺旋桨无空化噪声的预报与校验 杨琼方 王永生 海军工程大学船舶与动力学院 武汉 均匀来流下的螺旋桨无空化噪声对应为单极厚度噪声和偶极Gutin噪声。对于船舶螺旋桨而言，Gutin噪声可忽略不计。非均匀进流时，螺旋桨无空化辐射声场主要对应为桨叶非定常负载产生的偶极噪声（负载噪声）。 对于低马赫数流动的螺旋桨噪声而言，其辐射声场对源流场的反馈影响作用可以忽略不计，同时，湍流四极声源发声对声远场的贡献量相比于负载偶极声源也可以忽略不计，因此，可以采用计算流体力学（Computational fluid dynamics, CFD）和计算声学(Computational acoustics， CA)分步耦合方法来预报螺旋桨远场辐射噪声。流场CFD计算可采用不可压流动假设，计算得出桨叶表面的时域脉动压力，但声场计算时必须考虑可压延迟的影响。 为满足源流场的计算精度，CFD计算一般采用非定常雷诺时均（URANS）、尺度适应模拟（Scale-Adaptive Simulation, SAS）、分离涡模拟（Detached Eddy Simulation, DES）或大涡模拟（Large Eddy Simulation, LES）方法，度依次提高。现有的船舶工程应用中辐射声场计算多采用基于Ffowcs Williams-Hawkings (FW-H)方程（包括porous FWH方程和Kirchhoff-FWH方程）的面积且至少对应于螺旋桨转一圈进行分析。而频域方法是直接预测螺旋桨各谐次频率处噪声强度，将桨叶上的旋转源分布用桨盘面上的静态分布进行代替，这样处理的优势在于在计算中消除了源运动的旋转分量，使得作为时间函数的运动源分布被其傅利叶序列所代替，从而直接采用模态分量参与计算，易于分析。对于无限定边界声场计算而言，直接边界元法（Direct boundary element method, DBEM）频域内计算可以较好地解决工程问题。 Seol等(2002, 2005)利用面元法求解流场和FWH积分方法求解声传播的分步耦合方法计算了非均匀进流条件下DTRC 4119无侧斜标准桨无空化噪声和空化噪声的总声级中厚度噪声和负载噪声各自所占的份额以及声指向性，其计算结果证明了非空化条件下厚度噪声可以忽略不计的结论。本案例同样选用该桨，但采用LES模拟结合边界元数值声学的混合方法，计算软件采用LMS Virtual lab中acoustics模块。 1 1 流噪声预报的理论模型 近几年国内外预报流体噪声(包括湍流噪声和螺 旋桨噪声)的基本思路及可用方法如图1所示。Ω代表计算域，Γ代表其边界。下标S、F和A分别表示声源、流体和声。积分法预报主要对应为该图左边的时域计算，而以BEM数值声学方法为代表的混合法预报则对应为右边的频域计算。本案例采用在频域内分步耦合的方法进行分析。 ΩA ΓSΓF ΩF ΓA 图1 流体噪声耦合方法预报思路及工程中可用方法 瞬态声模拟直接耦合 谐次声模拟分步耦合 图2 流-声耦合计算时分析流程 因该桨几何相对比较简单，容易同时满足大网格密度、小网格尺度和好网格正则度的要求，所以在流场粘性模拟时选择LES模拟，且采用动力学Smagorinsky-Lilly（DSM）模型来模拟亚格子应力： τij=δijτkk?2μtSij （1） 其中，δij为克罗内克符号，τkk为各向同性的亚格子应力，μt为亚格子湍流粘度，Sij为所求解尺度涡的应变率张量。 μt=(CΔ)213 采用边界元数值声学方法频域内计算时，fs对应为桨叶及桨毂表面的脉动压力。当桨毂两端面取为封闭面时，螺旋桨噪声对应为外部辐射声问题，可采用直接边界元积分公式，其一般形式为： rrC(r)?p(r)= Ωa 1???j ，Sij=?i+ 2???xj?xi? ?，=2SijSij ?? ∫ rr rrr?r?r?G(r,ra) +jρωG(r,ra)?vν(ra)?dΩ(ra)（4）?p(ra)? ?ν?? 其中，V为积分式所取的有界流域体积，p(r)为点r处 ?jkr?r rrrre 声压，p(ra)为点ra处压力，G(r,ra)=为格林 4πr?ra a rr rr C为Smagorinsky系数，本文计算中，C=0.1；Δ为 局部网格尺度，取为单元体积的1/3。 在得到声源节点