矩形腔体流场模拟及噪声研究

用大涡模拟方法对低速湍流引起的矩形腔体内流动进行了模拟,并应用FW-H声学类比方程分析了由流动诱发的气动噪声.数值模拟观察到了涡结构的脱体及腔体内部的自激振荡过程,通过分析得出了由流动诱发噪声的声压-频率曲线.研究发现在流速30 m/s时,流动噪声声压级在60 dB以下,348.48 Hz及其高次谐波是噪声的主要来源,流场与声场表现出耦合关系,辐射声场具有明显的方向性.腔体噪声的风洞实验研究得到了与数值模拟吻合的结果.     

超声速三维空腔流气动噪声被动控制

抑制超声速武器舱空腔流噪声是航空领域中一项重要课题。大量研究表明在空腔前缘采用主/被动控制技术可以在一定程度上抑制腔内噪声水平。利用大涡模拟(large eddy simulation, LES)技术计算分析了Mach 1.4开式矩形方腔及波形、弧形两种前后壁几何修形后空腔的流动及噪声, 探索超声速来流条件下几何修形被动控制技术对开式方腔流噪声的抑制能力。计算结果表明波形和弧形空腔对腔内噪声均具有一定的抑制作用, 且波形空腔噪声控制效果更优。分析认为空腔几何修形能够改变空腔上方剪切层及腔内大尺度涡结构的发展演化, 进而实现对腔内噪声的控制。此外, 还应用LES方法计算分析了增厚的来流边界层条件下超声速方腔流, 发现来流边界层增厚可显著降低腔内噪声水平。      

开式轴流风扇气动噪声预测

本文采用LES/FW-H的匹配方法,研究了开式轴流风扇内部旋涡流动特征及其与叶片表面干涉引起的气动噪声之间的联系,同时进行了远场噪声预测,探讨了叶轮不同表面辐射噪声时的频谱分布特征.研究结果表明,开式轴流风扇吸力面附近形成的叶尖涡和前缘分离涡在吸力面叶片表面相应位置形成大压力波动,形成主要噪声源;叶片吸力面的辐射噪声可以通过改善吸力面附近的旋涡流动来降低;低速轴流叶轮由叶轮壁面辐射的噪声以宽频成分为主.     

贯流风机涡结构与噪声特性的数值研究

本文采用大涡模拟对空调室内机中贯流风机的内流进行了数值计算,结果显示了以偏心涡为代表的复杂非定常流动细节.计算结果表明,贯流风机气流两次进出叶轮,叶片尾缘、蜗舌处出现明显的脱落涡结构.叶轮周围监测点上出现了三个特征频率,分别对应叶片通过频率、叶片脱落涡频率以及蜗舌后缘的脱落涡频率,不同的特征点上表现出不同的频谱特性.另外,通过LEE方程中的声源项分布得到了贯流风机的主要噪声源区域,继而对蜗舌和叶轮等主要音源表面的远场辐射噪声频谱进行了分析,为后续的实验研究提供参考.     

喷管内壁多孔处理的圆口喷流大涡模拟

本文对一亚音速圆口喷管内壁面施加盲孔,使其壁面多孔化,采用大涡模拟结合FW-H方法研究了马赫数M=0.6时喷管外的非稳态流场与远场噪声。研究结果显示,这种内壁小孔结构对喷管外流场的湍流脉动有一定的削弱作用。时空相关系数的对比表明,内壁小孔结构同时减弱了脉动量各自在时间和空间上的相关性,削弱了大尺度的涡结构。通过对涡强、膨胀率和Q涡识别准则的可视化分析,研究了剪切层中涡结构在喷流轴向的产生、发展与演化过程,结果表明内壁小孔结构抑制了涡配对现象,加强了流动掺混,从而加速了剪切层的发展。     

自激脉冲喷嘴发生机理数值模拟

本文揭示了自激脉冲喷嘴产生自激振荡的根本原因。由于腔内产生了低于汽化压力的负压,引起低压区液体的汽化,加上射流不稳定剪切层的作用,最终形成大小周期性变化的轴对称空化气囊;对脉冲喷嘴腔内流场进行了二维非定常数值模拟,通过对腔内压力和流场的分析,验证了上述结论;试验研究通过对腔内压力的数据采集及对脉冲现象的同步观察,进一步说明了本文结论的正确性。     

分流流道内流动结构的测试与分析

分叉管道常见于工程上的流量分配装置、飞行器内外涵道结构,其中分流流道的流动结构影响着主流道的通流量,且回流涡的产生与扩大会引起分流流道的阻塞,使之失去分流的作用。本文采用高速摄影对分流流道内部的流动结构进行识别,发现在其入口处存在明显的回流涡,进而对回流涡处壁面进行了压力动态测量。压力动态测试与高速摄影的结果显示出分流通道中回流涡流场变化的频率信息。在不同的管道入口雷诺数(Re_(in)=80249到179414)下,回流涡处的壁面压力变化具有混沌特性,随着入口雷诺数的增加,压力脉动的幅值增大,而随机性却减小,确定性和稳定性增强,且在回流涡产生位置尤为显著,影响整个分流流道的通流能力。     

开孔结构流致噪声的数值模拟和机理分析

飞机机体表面的开孔设计会形成空腔结构,产生空腔流致噪声。空腔噪声的控制需要彻底认识其流动和噪声机理。以飞机的功能性开孔为例,通过半经验公式分析了其空腔噪声频率随速度的变化规律,预测了出现流声共振的工况。空腔发生流声共振时,特定频率的纯音噪声会被放大。为此,采用脱体涡模拟方法开展了开孔结构流声共振的三维非定常数值计算,分析了其流场和声场特性。其中,数值方法的准确性通过圆形空腔标模计算进行验证。结果表明,在一定速度下剪切层内的扰动将诱发空腔深度方向声模态,出现流声共振现象。此时,剪切层表现为强烈的周期性上下拍动,空腔底部和后缘区域的局部压力脉动幅值较大,声波主要由空腔后缘向上游方向辐射,上游噪声大于下游。     

槽道流POD重构及湍流动能耗散率分析

采用二维粒子图像测速仪（2DPIV）对槽道内涡波流场进行实验研究,用POD技术对2DPIV瞬态速度矢量场进行主导模态重构,得到槽道内的平均流速和湍流动能分布;采用大涡PIV方法对湍流动能耗散率分布进行计算.结果表明：重构流场表征了原始流场的主导结构,剔除了噪声等干扰信息;大涡PIV方法能有效地估算动能耗散率的分布;湍流动能在壁面附近较小,在接近槽道中心区域湍流动能越来越大,呈现出射流的特征;动能耗散率的峰值出现在壁面附近和槽道中心区域,动能耗散率随着远离壁面程度的增加先降低后逐渐增加直至达到峰值.     

襟翼侧缘噪声机理及修型降噪设计

襟翼侧缘噪声是飞机起降阶段机体噪声的重要噪声源。采用极大涡模拟对襟翼侧缘非定常流场进行数值模拟,分析其噪声产生机理.基于此,提出了两种襟翼侧缘修型方式,应用虚拟渗透面的Ffowcs Williams and Hawkings(FW-H)声比拟方法将修型构型的远场噪声频谱特性和指向性与基准构型对比分析,研究其降噪效果。通过流场和声场的数值模拟表明,襟翼侧缘噪声属于宽频噪声。不同的襟翼侧缘形状改变了流场形态、侧缘涡结构以及涡系的发展过程,进而对声源分布和远场噪声特性产生影响。结果表明:在给定的5°计算迎角下,两种襟翼侧缘修型方式在保证增升装置的原有升阻气动特性的前提下,能达到减小全场总声压级1~2 dB的降噪效果。      

利用质量流注入的流激孔腔噪声共振峰抑制

通过数值仿真揭示了开口前缘垂直注入质量流和前壁面平行注入质量流抑制流激孔腔噪声的机制,研究了多参数影响下脉动压力峰值降噪量和总降噪量随质量流注入速度的变化规律。开口前缘垂直注入质量流通过抬升剪切层,避免漩涡冲击开口后缘,抑制流激孔腔噪声脉动压力峰值;在一定范围内质量流注入速度越大,脉动压力峰值降噪量越大,但是低频部分引起的抬升也会越高,导致总降噪量先增大后减小;经优化后的峰值降噪量和总降噪量分别可以达到15 dB和9.5 dB。开口前壁面平行注入质量流则是通过加强开口处剪切层的稳定性,避免发生漩涡脱落,达到抑制流激孔腔噪声的目的;当质量流入口面积大于孔腔开口前壁面积2/3时,不仅可以显著降低流激孔腔噪声脉动压力的峰值,并且可以很好地抑制其它频段噪声的抬升;质量流注入速度为来流速度的0.5倍时,脉动压力峰值降噪量和总降噪量分别可以达到18 dB和15.4 dB。     

对称槽道涡波流场流动特征的POD分析

为深入分析层流状态下对称槽道内涡波流场的流动特性及其变化规律,对流场进行了二维粒子图像测速(2DPIV)测量获取瞬态速度矢量数据,利用本征正交分解(POD)技术进行模态分解以及涡波流场的重构,然后根据重构的流场对对称槽道内涡波流场进行了平均速度剖面、流场脉动强度以及特征点的速度和频谱分布等方面的分析。结果表明:POD的前15阶模态能够表征涡波流场的主导结构,第1,3阶模态主要表现为一对旋向相反的涡对特征,第2阶模态具有涡旋和波状主流的特征;提取了5个涡旋涡核的位置作为流场流动特性的特征点;根据POD重构流场分析发现流向平均速度呈抛物线形状分布,法向平均速度呈对称分布特征;流向脉动强度受壁面的影响较大,法向脉动强度呈现抛物线形状分布;距离中心主流较近的1#,4#,5#特征点的速度脉动程度受主流的脉动强度影响较大,速度的脉动主频0.15 Hz与次频、流场的自然频率0.35 Hz共同影响特征点的速度分布;2#,3#特征点的流向速度呈衰减趋势,法向速度在初期幅度变化较大。     

DMD和POD对超燃冲压发动机凹腔流动的稳定性分析

开式凹腔作为超燃冲压发动机中增加掺混和稳焰的装置, 其流动稳定性的研究对深入理解凹腔增加掺混和稳焰机理以及凹腔的设计有着重要的学术意义和工程应用价值.基于大涡模拟方法对超燃冲压发动机开式凹腔流动进行数值模拟, 分别采用动力学模态分解(dynamic mode decomposition, DMD)和本征正交分解方法(proper orthogonal decomposition, POD)对自激振荡流动进行稳定性分析. DMD方法可准确提取凹腔的振荡频率, 与Rossiter模型以及压力脉动FFT分析得到的频率吻合较好, 且DMD中对应Rossiter前3阶频率的模态在流动中的主导作用顺序也与FFT分析结果一致, 自激振荡中RossiterⅢ模态占据主导作用, 同时DMD方法对Rossiter 3阶以上模态频率的预测能力明显强于FFT分析方法.在对低频的提取方面, DMD方法比Rossiter模型更具有优势.与前6阶Rossiter模态对应DMD模态均缓慢收敛, 主要表现为剪切层中的分离涡结构和中部及下游区域中的涡结构.前3阶不稳定模态中的分离涡结构主要集中在中部剪切层以及后缘附近区域. POD方法中较少的模态包含流场绝大部分的能量.但是, 通过POD方法提取的模态频率在分辨率上效果不佳, 提取到最低频率为Rossiter 3阶模态对应的频率, 且模态中均存在次频, 次频与主频之间的耦合导致模态的形态相差较大.另外, 与DMD方法相比POD方法无法判断所提取的模态的稳定性.      

Attenuation of sound due to vortex shedding from a splitter plate in a mean flow duct

A theoretical analysis is given of an experiment being performed at the University of Southampton [1] as part of a programme to quantify the effectiveness of perforated screens in dissipating sound in the presence of tangential mean flow. In the experiment vorticity is generated at the trailing edge of a splitter plate in a mean flow duct by a plane sound wave incident from upstream, acoustic energy being ceded to the kinetic energy of the vortex field. An expression is derived for the dissipated sound power at arbitrary subsonic mean flow Mach number and frequency. The calculation is performed both by a consideration of the net flux of acoustic energy into the trailing edge region of the splitter plate, and by evaluating the rate of working of the vortex lift forces in the field of the acoustic particle velocity. In particular, it is shown that the absorption is independent of frequency, provided the frequency does not exceed the minimum cut-on frequency of transverse acoustic modes within the duct.     

An aeroacoustically driven thermoacoustic heat pump

The mean flow of gas in a pipe past a cavity can excite the resonant acoustic modes of the cavity--much like blowing across the top of a bottle. The periodic shedding of vortices from the leading edge of the mouth of the cavity feeds energy into the acoustic modes which, in turn, affect the shedding of the next vortex. This so-called aeroacoustic whistle can excite very high amplitude acoustic standing waves within a cavity defined by coaxial side branches closed at their ends. The amplitude of these standing waves can easily be 20% of the ambient pressure at optimal gas flow rates and ambient pressures within the main pipe. A standing wave thermoacoustic heat pump is a device which utilizes the in-phase pressure and displacement oscillations to pump heat across a porous medium thereby establishing, or maintaining, a temperature gradient. Experimental results of a combined system of aeroacoustic sound source and a simple thermoacoustic stack will be presented.     

The effect of a cavity on airfoil tones

The presence of a cavity in the pressure surface of an airfoil has been found via experiment to play a role in the production of airfoil tones, which was attributed to the presence of an acoustic feedback loop. The cavity length was sufficiently small that cavity oscillation modes did not occur for most of the investigated chord-based Reynolds number range of 70,000–320,000. The airfoil tonal noise frequencies varied as the position of the cavity was moved along a parallel section at the airfoil's maximum thickness: specifically, for a given velocity, the frequency spacing of the tones was inversely proportional to the geometric distance between the cavity and the trailing edge. The boundary layer instability waves considered responsible for the airfoil tones were only detected downstream of the cavity. This may be the first experimental verification of these aspects of the feedback loop model for airfoil tonal noise.     

Vortex shedding noise of low tip speed,axial flow fans

Noise and performance tests were conducted on three low tip speed, half-stage, axial flow fans to determine the nature of the vortex shedding noise mechanism. Each fan was 356 mm in diameter and had eight equally spaced, variable pitch blades. The noise measurements were made in a free field environment and the fan back pressure and speed were varied during the tests. An acenaphthene coating on the blades was used to determine the regions of laminar and turbulent flow.Vortex shedding can be a significant source of noise when the fan is operated in a lightly loaded condition. Essentially it is due to instabilities in the laminar boundary layer on the suction side of the blade where these instabilities are in the form of Tollmien-Schlichting (T-S) waves. These instabilities interact with the trailing edge of the blade and generate acoustic waves which radiate from the trailing edge and form a feedback loop with the source of the instabilities. Vortex shedding noise can contribute as much as 5 dB in overall noise level and up to 22 dB at higher frequencies (8–14 kHz).Serrations located at the leading edge, at the mid-chord, or near the trailing edge on the suction side were found to reduce the vortex shedding noise significantly. The mid-chord location was found to be the most satisfactory because, as well as eliminating the noise, the serrations provided a 3% improvement in peak efficiency. This improvement occurred because separation of the laminar boundary layer was prevented on the suction side. On the other hand, serrations placed at the other two locations tended to degrade fan performance.     

Azimuthal behaviour of flow-excited diametral modes of internal shallow cavities

The aerodynamic excitation of ducted cavity diametral modes gives rise to complex flow-sound interaction mechanisms, in which the axisymmetric free shear layer interacts with the asymmetric acoustic modes. This results in various azimuthal patterns and behaviours depending on different flow and geometrical parameters. The azimuthal behaviour of this self-excitation mechanism is investigated experimentally. Axisymmetric shallow cavities in a duct have been tested over the range of cavity length to depth ratio from 1 to 6 and at Mach numbers up to 0.4. A set of pressure transducers flush mounted to the cavity floor is used to determine the acoustic mode amplitude and orientation. The excited acoustic modes are classified into spinning, partially spinning, and stationary diametral modes. An analytical representation based on the duct acoustics theory is used to analyse the measurements and provides a physical explanation of the observed behaviour of the diametral modes. Splitter plates are installed inside the cavity to form a geometrical preference. The acoustic response of this geometrically altered case show that pressure oscillations at different azimuthal angles along the cavity circumference can be uncorrelated, or even oscillate at different frequencies, while the diametral modes are still strongly excited. Two hot-wire probes are also used in a separate set of measurements to investigate the azimuthal behaviour of the shear layer oscillation. The results show that the shear layer oscillation has the same azimuthal distribution as that of the excited acoustic modes, indicating that the shear layer oscillation at different azimuthal angles can be not only uncorrelated but also occur at different frequencies.     

Experimental investigation of ethylene/air combustion instability in a model scramjet combustor using image-based methods

The combustion instabilities of supersonic combustion were investigated experimentally in a laboratory-scale scramjet combustor with a cavity flame holder. Ethylene was injected transversely from an orifice to the supersonic flow of Mach 2 with a stagnation temperature of 1900 K and a total pressure of 0.37 MPa. The dynamic pressure, CH* chemiluminescence and shadowgraph images were measured with a pressure sensor and a high-speed video camera. Dynamic pressure was analyzed by fast Fourier transform, and time-resolved CH* chemiluminescence images were modally decomposed by the sparsity-promoting dynamic mode decomposition (SP-DMD). The results indicated that two combustion instabilities were observed for cavity shear-layer stabilized combustion and the oscillation between jet-wake stabilized and cavity shear-layer ram combustions for the power spectral density (PSD) of pressure. In the case of the combustion instability of cavity shear-layer stabilized combustion, a dominant peak of approximately 128 Hz was observed for the PSD of pressure. This instability corresponded to an entire flame oscillation of the cavity shear-layer stabilized combustion, which was validated by the SP-DMD and a low rank reproduction with 10 modes. This was driven by a fuel injection oscillation in the injection orifice. In the case of oscillation between the jet-wake stabilized and the cavity shear-layer ram combustions, peaks around 1600 Hz were observed for the PSD of pressure. This mechanism was also explained by the SP-DMD modes and a low rank reproduction using within 10 modes. The DMD and shadowgraph images indicated that the vortex formed by a separation of the boundary layer induced a strong jet-wake flame, resulting in the temporal thermal choke followed by cavity shear-layer stabilized ram combustion. The data-driven approach with SP-DMD clarified the combustion instability mechanisms of the supersonic combustion in detail.