非紧致气动噪声半空间传播基本解的边界积分方法

为了获得低马赫数流动诱发的非紧致气动噪声在半空间内传播的基本解,结合复等效源方法和边界元方法建立了半空间精确格林函数的边界积分方程,当半空间边界为刚度型阻抗边界时可避免奇异性积分。同时基于等效源方法提出一种半空间二维非紧致圆柱声散射模型,推导了静止介质中声散射基本解的理论表达式。对静止介质中的二维圆柱声散射,数值解在研究的频率与观察点处与理论解一致。采用数值方法计算二维圆柱绕流诱发的半空间声场基本解,结果显示半空间边界强化了声源的声辐射能力,小于马赫数0.2的介质运动对声传播的影响可以忽略。     

基于声扰动方程的气动噪声传播积分计算方法

传统声比拟思想假设声传播介质均匀,所获得的远场噪声积分方法无法考虑近场流动非均匀性引起的声折射效应,声场预测结果往往存在较大的数值误差。为解决上述问题,本文采用声扰动方程描述声波运动,并利用声学互易原理,在频域内构建关联近场声源与远场声压的伴随格林函数,从而将非均匀介质折射效应纳入伴随格林函数。基于伴随格林函数和声扰动方程源项的数值结果,利用积分方程计算气动噪声的空间传播。结合上述方法对二维圆柱和NACA0012翼型绕流噪声开展了数值预测研究,远场声压预测结果与CFD直接计算值吻合。数值研究表明,本文提出的基于声扰动方程的气动噪声传播积分计算方法可以考虑非均匀介质对声传播的影响,能够提高气动噪声的预测精度。     

方腔流动气动噪声边界散射数值预测方法研究北大核心CSCD

方腔流动噪声问题因其流场变化复杂且剧烈而倍受关注,本文主要开展高亚声速方腔流动气动噪声数值预测方法研究。基于任意边界条件的格林函数解和Lighthill声模拟理论,提出可以考虑空间边界影响的气动噪声积分计算方法。数值模拟包含流动和噪声计算两部分,通过二阶精度的DDES模型进行流动数值模拟,边界积分方法计算散射声场分布。数值结果显示声场分布随时间呈现周期性变化,与流场的脉动及其脉动周期一致。观察点的声压级随频率逐渐下降且在谐波频率突然增大。本文计算结果与高精度计算气动声学方法计算结果相符,表明该方法合理、可靠,并且具有较高的计算效率。     

方腔流动气动噪声边界散射数值预测方法研究

方腔流动噪声问题因其流场变化复杂且剧烈而倍受关注,本文主要开展高亚声速方腔流动气动噪声数值预测方法研究。基于任意边界条件的格林函数解和Lighthill声模拟理论,提出可以考虑空间边界影响的气动噪声积分计算方法。数值模拟包含流动和噪声计算两部分,通过二阶精度的DDES模型进行流动数值模拟,边界积分方法计算散射声场分布。数值结果显示声场分布随时间呈现周期性变化,与流场的脉动及其脉动周期一致。观察点的声压级随频率逐渐下降且在谐波频率突然增大。本文计算结果与高精度计算气动声学方法计算结果相符,表明该方法合理、可靠,并且具有较高的计算效率。     

复杂几何条件下伴随格林函数的数值求解

为简化气动噪声预测中复杂边界的网格处理,采用计算气动声学(CAA)方法与浸入式边界方法(IBM)相结合数值求解气动噪声预测所需的伴随格林函数.根据伴随格林函数的基本形式,文中设计了特定的圆柱声散射算例,并将数值结果与解析解对比,分析该方法的计算精度,验证了该方法在复杂几何边界条件下伴随格林函数求解中的适用性.最后应用该方法求解了考虑吊挂安装效应的锯齿型喷管喷流剪切层内的伴随格林函数.结果表明:由于喷流剪切层的散射效应,伴随格林函数在喷流内分布不均匀,最大值比最小值高出3倍以上,反映了喷流内不同区域声源对远场噪声贡献的差别,降噪设计可以重点考虑降低格林函数较大区域的声源强度以达到降低噪声的目的.      

非紧致气动噪声传播的研究进展与分析

非均匀运动介质中非紧致气动噪声的传播是气动声学的一个研究热点与难点。本文围绕理论模型和数值方法,对非紧致气动噪声传播的研究进展进行了综述与分析。首先回顾了声传播理论模型的发展历程,对各模型的优缺点进行了总结分析;然后着重讨论了近20年来采用混合计算气动声学方法研究非紧致气动噪声传播的数值预测方法,包括有限元法、边界元法和精确格林函数法;最后对非紧致气动噪声传播数值方法的发展进行了展望。     

一种基于平均流场的翼型尾缘宽频噪声预测方法

发展了一种基于平均流场的翼型尾缘宽频噪声预测方法.可压缩Navier-Stokes方程被重组为能描述气动噪声产生和传播的形式,包括准确描述声传播的线性算子和雷诺应力构成的声源项.为减少计算时间和避免格林函数空间奇异性,采用伴随格林函数方法求解格林函数;同时采用基于平均流场的声源模化方法生成声源项.为了验证该方法,预测了NACA0012翼型尾缘宽频噪声,其中平均流场使用OpenFOAM计算,格林函数采用均匀流动下半无限大平板近似.预测结果表明:对于不同来流速度和攻角,该方法不仅能准确预测噪声幅值及峰值频率,而且噪声频谱也都在实验结果的误差范围之内(±3dB),证明了该预测方法的精度和适用性.      

飞机螺桨声场的数值模拟方法

本文给出了预测飞机螺桨在自由空间内幅射的气动声场的一种数值方法。即时域内进行延迟时间积分的一类方法。在该方法中使用片条假设,并考虑厚度声源和负荷声源项作为非紧致声源。计算与实验结果对比表明,该方法和相应软件可作为一种工程手段来使用。      

离心风扇非定常气动力和气动噪声特性数值研究

将Lighthill方程转变为频域Helmholtz弱积分形式并采用Galerkin方法离散.基于声学有限元方法考虑声波在复杂固壁(叶轮和蜗壳)内的散射和反射等作用,利用Ffowcs Williams-Hawking(FW-H)方程耦合非定常流场计算结果数值预测了某离心风扇的噪声辐射,流场计算结果和蜗壳壁面动态压力测量结果在基频上吻合较好.结果表明:基频压力脉动分量在噪声源特性中占据主导地位且靠近叶轮前盖板对应位置的蜗舌区域(叶轮出口宽度范围内)是最主要的噪声源区域;声学有限元方法和实验吻合较好,复杂固壁对声传播影响不容忽略.叶轮出口不稳定气流对蜗壳周期性冲击引起的转/静干涉噪声远大于叶片偶极子源噪声是离心风扇最主要的噪声辐射分量且噪声主要从风扇管道出口方向传播.      

远场假设对喷流噪声预测中格林函数求解的影响

目前基于雷诺平均Navier-Stokes（RANS）的喷流噪声预测方法在格林函数求解时, 为简化求解过程,通常对喷流流动做平行流假设, 对观测点做远场假设。随着格林函数求解方法发展,近年来的研究表明平行流假设对下游观测点格林函数的计算会引起较大偏差,而目前远场假设对格林函数求解的影响仍不清楚。为研究远场假设对喷流格林函数求解的影响,以二维喷流为例, 采用计算气动声学方法（CAA）分别数值求解了观测点远场假设条件与实际条件下90°~150°方向喷流内伴随格林函数,进而分析远场假设对格林函数求解的影响。研究结果表明,对于不同方向的观测点,由观测远场假设导致的伴随格林函数求解偏差不尽相同,且对于越靠近喷流中心线方向的观测点,远场假设导致的偏差越大,其中150°方向观测点,采用远场假设后,格林函数计算结果最大偏差达到-15 dB以上。因此,对于靠近喷流中心线方向的噪声观测点而言,为避免预测偏差,应采用实际观测条件求解喷流格林函数。     

Aeroacoustic integral methods,formulation and efficient numerical implementation

The paper makes a contribution to the clearer understanding of the physical meaning and domain of applicability of the Ffowcs Williams–Hawkings (FW-H) and Kirchhoff aeroacoustic integral methods. The analytical relationship between the two approaches, and the physical implications involved in the use of the Kirchhoff approach are discussed. In particular, where the Kirchhoff surface cuts into a domain where non-negligible volume sources of Lighthill's Acoustic Analogy are present (domain crossed by a shock, non-uniform flow or vortices for instance), the Kirchhoff integral must be complemented by additional surface integrals to recover the FW-H surface integral. The paper goes on to describe a fast, robust integration technique for computing FW-H or Kirchhoff surface integrals. This integration technique starts from emission time, and treats each integration surface element as non-compact. It obviates the need to build retarded acoustic surfaces and avoids some of the drawbacks of current methods, such as the Doppler singularity when using supersonically rotating grids. The paper then recalls formulas for calculations in the aircraft reference frame, and the computational efficiency of the proposed integration method is demonstrated by its application to the prediction of helicopter rotor noise directivity contours.     

Numerical analysis on noise of rotor with unconventional blade tips based on CFD/Kirchhoff method

A solver is developed aiming at efficiently predicting rotor noise in hover and forward flight. In this solver, the nonlinear near-field solutions are calculated by a hybrid approach which includes the Navier-Stokes and Euler equations based on a moving-embedded grid system and adaptive grid methodology. A combination of the third-order upwind scheme and flux-difference splitting scheme, instead of the second-order center-difference scheme which may cause larger wake dissipation, has been employed in the present computational fluid dynamics (CFD) method. The sound pressure data in the near field can be calculated directly by solving the Navier-Stokes equations, and the sound propagation can be predicted by the Kirchhoff method. A harmonic expansion approach is presented for rotor far-field noise prediction, which gives an analytical expression for the integral function in the Kirchhoff formula. As a result, the interpolation process is simplified and the efficiency and accuracy of the interpolation are improved. Then, the high-speed impulsive (HIS) noise of a helicopter rotor at different tip Mach numbers and on different observers is calculated and analyzed in hover and forward flight, which shows a highly directional characteristic of the rotor HIS noise with a maximum value in the rotor plane, and the HSI noise weakens rapidly with the increasing of the directivity angle. In order to investigate the effects of the rotor blade-tip shape on its aeroacoustic characteristics, four kinds of blade tips are designed and their noise characteristics have been simulated. At last, a new unconventional CLOR-II blade tip has been designed, and the noise characteristics of the presented CLOR-II model rotor have been simulated and measured compared to the reference rotors with a rectangular or swept-back platform blade tip. The results demonstrate that the unconventional CLOR-II blade tip can significantly reduce the HSI noise of a rotor.     

基于T-S模糊模型的非线性分布参数系统极点配置P-sD控制CSCD

采用Lyapunov直接法讨论分布参数系统的稳定性,建立分布参数系统的Lyapunov函数。为保证分布参数系统稳定,应使Lyapunov函数对时间的微分小于0。由于分布参数系统中存在空间一次微分项与二次微分项,Lyapunov函数对时间的微分中将出现常数项与积分项,针对常数项,引入空间一次微分项来抵消;针对积分项,引入对应的状态反馈来使系统稳定。利用边界条件量化Lyapunov函数对时间的微分中的各项,从而设计控制器,这是一种新的设计P-sD状态控制的方式。其中状态反馈的部分采用极点配置的方法来设计,当分布参数系统中出现状态变量的非线性函数时,采用T-S模糊模型表达,可以通过状态变量的线性组合精确描述非线性项,以便极点配置设计状态反馈。针对两个非线性分布参数数学模型进行仿真,结果证明设计的P-sD状态控制器可以使分布参数系统稳定,并达到期望的效果。     

基于边缘射流的起落架气动噪声控制研究

:起落架噪声可以看成是一系列结构件单独引起的气体扰流噪声以及这些不同结构件相对位置引起的干扰噪声的耦合.为了降低起落架气动噪声,提出一种基于边缘射流的主动控制技术.以某型飞机前起落架为研究对象,在其扭力臂背风面施加射流,利用分离涡脱模拟方法对其支柱和扭力臂结构简化模型的周围流场进行非定常计算,获取声源分布,采用FW-H积分获得远场噪声特性.结果表明:边缘射流能够有效抑制干扰噪声和支柱噪声,起落架的中频噪声得到一定幅度的下降,宽频噪声强度也有所减弱;射流改变扭力臂尾涡的脱落状况,可以减轻甚至消除涡脱落对支柱的冲击,从而减弱了支柱表面由于撞击而产生的脉动压力,达到降低声源强度的目的.     

基于基尔霍夫面的四维声学频域公式

流动诱发噪声问题是实际工程领域极为普遍的问题之一,经典的声比拟模型仅以声学压力为参考来评估声场特征分布还远远不够。从声压和声学速度矢量为变量的四维线性波动方程出发,选择包围非线性声源的基尔霍夫面为积分面,并结合对流格林函数,给出均匀运动介质的四维声学频域积分公式;针对静止、旋转单极子源和偶极子源开展数值预测研究。结果表明：本文获得的声压和声学速度分布与解析解吻合,均匀来流情形下静止点源的声场分布表现出典型的对流效应;受均匀来流、点源的自激频率、谐波阶次和旋转频率的共同影响,旋转点源的声场分布则表现出明显的多普勒效应和对流效应。本文对均匀流矢量气动噪声开展的精细化研究能够为声能量的评估和传输路径预测提供技术支撑,进一步为降噪研究提供理论依据。     

螺桨飞机缩比座舱模型噪声传播的试验研究

本文使用Ｙ１２飞机的１／３缩比模型桨进行了螺桨地面噪声试验，使用圆柱型缩比座舱模型测量了舱内声压。试验中发现时间平均增强法是去除地面湍流影响的有效方法，利用舱内传声器阵可以测量舱内声场分布，经改进，使用户外螺桨模型噪声试验来作为螺桨噪声预测结果的验证手段是可行的。     

Transfer element method with application to acoustic design of aeroengine nacelle

In the present survey, various methods for the acoustic design of aeroengine nacelle are first briefly introduced along with the comments on their advantages and disadvantages for practical application, and then detailed analysis and discussion focus on a kind of new method which is called ‘‘transfer element method'(TEM) with emphasis on its application in the following three problems: turbomachinery noise generations, sound transmission in ducts and radiation from the inlet and outlet of ducts, as well as the interaction between them. In the theoretical frame of the TEM, the solution of acoustic field in an infinite duct with stator sound source or liner is extended to that in a finite domain with all knows and unknowns on the interface plane, and the relevant acoustic field is solved by setting up matching equation. In addition, based on combining the TEM with the boundary element method(BEM) by establishing the pressure and its derivative continuum conditions on the inlet and outlet surface, the sound radiation from the inlet and outlet of ducts can also be investigated. Finally, the effects of various interactions between the sound source and acoustic treatment have been discussed in this survey. The numerical examples indicate that it is quite important to consider the effect of such interactions on sound attenuation during the acoustic design of aeroengine nacelle.