激波和剪切层相互作用下的超音速射流

欠膨胀超音速射流处于螺旋模式下的中度欠膨胀时,其入射剪切层的激波具有很高强度,激波和剪切层发生了强烈的相互作用,远场辐射的拢动波出现了大间隔、交错的上下行类似螺旋锥面波形图像,该扰动波具有很强的向上游传播的指向性,导致上游噪声高于垂直喷嘴方向的声压级.而在相对压比较低的低度欠膨胀情况下,或高压比下的高度欠膨胀的情形,入射剪切层激波强度相对较弱,远场辐射没有大间隔、交错的上下行远场辐射 关键词： 超音速射流 啸叫 扰动波 激波 剪切层     

非定常欠膨胀射流的正格式数值模拟

 采用二阶正格式方法对非定常欠膨胀射流进行了数值模拟。将二维守恒方程的正格式方法推广到轴对称Euler方程组的求解,并对不同马赫数下的燃气射流进行了数值计算。计算结果表明,该方法能够较好地捕捉到包含膨胀波、入射激波、反射激波、马赫盘、射流边界以及三波点等复杂射流流场的波系结构,与实验照片反映的流动特征以及已有的数值结果相吻合。表明该方法对间断解具有较强的捕捉能力,在激波阵面上不会出现数值振荡。     

超声速后台阶流动/射流相互作用的数值模拟

采用高精度格式求解二维Navier-Stokes方程研究超声速射流与同向超声速后台阶流动相互作用的流场基本结构及规律,分别应用5阶WENO格式、6阶中心差分格式离散对流项和黏性项,时间推进采用3阶Runge-Kutta格式,并应用消息传递接口（message passing interface,MPI）非阻塞式通信实现并行化.分别研究了超声速后台阶流动、超声速射流的基本结构特征,以此讨论和分析超声速后台阶流动/射流相互作用的特征,以及不同来流条件对波系结构、涡结构、剪切层、膨胀扇等的影响,尤其是来流剪切层和射流剪切层的相互作用,形成复杂的波系结构及相互干扰的流动现象.      

欠膨胀氢气射流激波结构数值模拟研究

本文利用数值模拟的方法研究了高压氢气泄漏所产生的欠膨胀氢气射流的激波结构。采用二维轴对称几何模型,分别模拟了不同压力比下,出口直径为1 mm的氢气射流。模拟得到了不同压力比下激波结构的形态和特征尺寸,包括马赫盘的位置、直径和边界层的厚度等。通过计算结果与实验测量的比较表明,当泄漏源压力小于3 MPa时,数值模拟不能够得到实验中观测到的典型激波结构;而当泄漏源压力达到3 MPa及以上时,数值模拟就可以成功地计算欠膨胀射流的激波结构.此外,模拟还得到了沿射流中心线的气流马赫数、压力和温度等一系列气流参数。     

合成冷/热射流控制超声速边界层流动稳定性

为了探究超声速边界层流动稳定性及其转捩控制机理,提出基于合成冷/热射流的边界层速度-温度耦合控制方法,并通过数值模拟研究了Ma=4.5超声速平板边界层不稳定波的传播,采用线性稳定性理论中的时间模式分析了壁面吹吸、射流温度、扰动频率、扰动振幅等对不稳定波控制效果的影响.结果表明:无射流控制时,边界层内同时存在不稳定的第一模态扰动波和第二模态扰动波,且二维波形式的第二模态占主导地位;壁面吹吸作用下,仅出现更加不稳定的第二模态,第一模态被抑制;速度-温度耦合控制下,射流温度对扰动模态的不稳定区域大小及扰动增长率影响显著,射流温度与来流温度不同时,温度的脉动使得流动转捩为湍流的速度加快,边界层速度型更加饱满,抗干扰能力增强,流动稳定性提高;高频的吹吸扰动对流场的控制效果优于低频扰动,扰动频率超过400 Hz时,第二模态扰动波时间增长率降低,扰动分量对边界层速度剖面和温度剖面的修正加快,第二模态更加稳定;扰动振幅减小为主流速度的1%时,仅出现时间增长率较小的第二模态,控制效果较好,进一步减小时,第一模态重新出现,并且波数范围与第二模态先重合后分离,对应的时间增长率先增加后减小.研究结果为边界层转捩控制技术提供了新的思路.     

波后滑移线引起的界面变形和R-M不稳定性

实验研究复杂波形结构引起平面界面变形和反射激波冲击下的R-M不稳定性的问题.在竖直激波管中生成稳定的N₂/SF₆平面界面,激波在圆柱绕射后,冲击平面界面,由此研究复杂激波引起的界面变形.平面激波在圆柱绕射后的流场,演化成具有初始入射波、三波点、弯曲反射波、Mach波和Mach反射产生的滑移线等复杂结构.研究复杂结构激波对界面的作用,对认识界面扰动的生成具有较大帮助.绕柱激波冲击后,平面界面仅在两对滑移线内部发生变形.绕柱激波冲击界面后,两对滑移线将界面分成"内界面"和"外界面",界面变形形态同滑移线和界面相交位置相关.反射激波二次冲击下,界面扰动的增长与Jacobs-Sheeley涡量模型较吻合.      

超音速射流啸叫模式切换不稳定性

超音速射流的噪声辐射有3个主要成分,即湍流的混合噪声、与激波结构有关的宽带噪卢以及离散的单频啸叫。其中啸叫离散基频特性具有多种模式,模式之间的切换是一个不稳定过渡过程,与流动纹影显示的不稳定相对应。模式切换过程中在频率特性上反映出一种模式的饱和、衰减,而另一种模式的强化,最终成为主导。这种切换过程中的频率变化很可能和流动失稳涡结构的破碎、配对相关联,利用线性不稳定性理论得到了与实验相符的不稳定波的增长率,分析了不稳定波最大增长率与啸叫频率切换的对应关系。     

两电极等离子体合成射流激励器工作特性研究

采用放电测量和高速阴影技术对两电极等离子体合成射流激励器工作特性进行了系统实验研究.实验表明:激励器工作击穿电压和放电峰值电流随激励器所处环境压强的降低和放电频率的增大而减小,激励器腔体内的放电过程为火花电弧放电.典型的等离子体合成射流流场包含有一道前驱激波和一股呈蘑菇状的高速射流.在整个射流发展过程中,前驱激波以当地声速恒速传播,不随激励器条件的改变而变化,波的强度则随着激励器出口直径的减小、腔体体积的增大、环境压强的降低和放电频率的升高而减小.激励器腔体体积和放电频率的增加会降低腔内气体的加热效果,并减小射流速度.激励器出口直径和环境压强对射流速度的影响按规律变化且存在最佳值.本文实各验条件下激励器都产生了明显的前驱激波和高速射流,具有实现高速流场主动流动控制的应用潜能.     

矩形喷口欠膨胀超声速射流对撞的实验研究

在不同喷口间距和射流压力下开展了矩形喷口欠膨胀超声速射流对撞实验并与自由射流进行了对比. 实验表明:超声速射流对撞的辐射噪声中存在四种不同的啸音模式, 且随喷口距离和射流压力的变化在不同模式间切换. 在射流压力大于0.5 MPa且喷口间距小于50 mm时, 射流对撞面在两个喷口外形成两道正激波之间, 啸音基频维持在3 kHz左右. 随喷口间距的增大或射流压力的降低, 射流对撞面在一侧喷口外的弓形激波与另一侧喷口外的正激波之间. 对撞面也有可能出现在两个弓形激波之间, 对应的啸音基频约为9 kHz, 但容易受扰动而回到喷口一侧或是在喷口之间大幅度振荡. 当射流压力小于0.36 MPa且喷口间距大于70 mm后, 对撞面在两个喷口之间大幅度振荡, 产生基频在1 kHz左右并随射流压力的降低和喷口间距的增大而降低的啸音. 关键词： 超声速射流 啸音 射流对撞 激波     

流经矩形喷嘴的超音速射流啸叫模式切换的实验研究

通过实验分析比较了对于相同高度不同宽度的四种矩形喷嘴, 当压力在0.2 MPa到0.8 MPa 之间变动时, 欠膨胀超音速自由射流的啸叫特性和对应的流场纹影结构.结果表明: 不同宽高比喷嘴的超音速自由射流辐射噪声中的单频离散啸叫存在两种不同的啸叫模式, 且随着射流压力的变化会出现模式间的切换.所谓模式切换是指不同模式的轮流占优和消失的现象.啸叫模式间的切换及占优区间的宽度随着喷嘴宽高比的减小而缩短.其中, 宽高比为2的射流啸叫模式中的一种模式所占的射流压降区间异常小, 此现象未在相关文献中提及; 喷嘴宽高比为4的射流啸叫占优区间内, 啸叫基频-射流压力曲线在0.49 MPa时出现了间断、跳跃现象.随着压力的降低激波纹影的轴线出现了抖动, 不同宽高比下流场结构的稳定性随压力变化的规律各异.射流压力在0.70 MPa到0.45 MPa区间内, 随着宽高比减小, 第一波节格栅的激波致密度减弱, 且出现轴向脉动, 第二波节后方的流场变得紊乱; 当射流压力低于0.45 MPa 时, 激波串结构随着宽高比的增大而趋于稳定, 在此压力区间内周期性激波格栅结构较射流压力在0.45 MPa以上时有所减弱.结合啸叫频谱及纹影图分析, 可初步认为, 第二、三波节也会对啸叫频率的声压幅值起到反馈增强作用.     

A PIV study of a supersonic impinging jet

The characteristics of supersonic impinging jets are investigated using Particle Image Velocimetry (PIV). The purpose of the experiments is to understand the jet induced forces on STOVL aircraft while hovering close to the ground. For this purpose, a large diameter circular plate was attached at the nozzle exit. The oscillations of the impinging jet generated due to a feedback loop are captured in the PIV images. The instantaneous velocity field measurements are used to describe flow characteristics of the impinging jet. The important flow features such as oscillating shock waves, slipstream shear layers and large scale structures are captured clearly by the PIV. The presence of large scale structures in the impinging jet induced high entrainment velocity in the near hydrodynamic field, which resulted in lift plate suction pressures. A passive control device is used to interfere with the acoustic waves travelling in the ambient medium to suppress the feedback loop. As a consequence, the large scale vortical structures disappeared completely leading to a corresponding reduction in the entrainment.     

Shock oscillation in highly underexpanded jets

The oscillatory motions of shocks in highly underexpanded jets with nozzle pressure ratios of 5.60, 7.47, 9.34, and11.21 are quantitatively studied by using large eddy simulation. Two types of shock oscillations are observed: one is the Mach disk oscillation in the streamwise direction and the other is the shock oscillation in the radial direction. It is found that the Mach disk moves quickly in the middle of the oscillatory region but slowly at the top or bottom boundaries. The oscillation cycles of Mach disk are the same for different cases, and are all dominated by an axisymmetric mode of 5.298 k Hz. For the oscillation in the radial direction, the shocks oscillate more toward the jet centerline but less in the jet shear layer, and the oscillation magnitude is an increasing function of screech amplitude. The cycles of the radial shock oscillation switch randomly between the two screech frequencies for the first two cases. However, the oscillation periodicity is more complex for the jets with high nozzle pressure ratios of 9.34 and 11.21 than for the jets with the low nozzle pressure ratios of 5.6 and 7.47. In addition, the shock oscillation characteristics are also captured by coarse mesh and Smagorinsky model,but the coarse mesh tends to predict a slower and weaker shock oscillation.     

稠密颗粒射流倾斜撞击颗粒膜特征

采用高速摄像仪对稠密颗粒射流倾斜撞击形成的类液体颗粒膜特征进行实验研究,考察了颗粒粒径、射流速度以及射流含固率等因素对颗粒膜形态及动态特征的影响.结果表明:随着颗粒粒径增大,稠密颗粒倾斜撞击流由颗粒膜向散射模式转变;随着射流速度增加,气固不稳定增强,射流流量出现脉动,正面与侧面分别表现为颗粒膜的非轴对称振荡和表面波纹结构;颗粒膜非轴对称振荡的振荡频率和振荡幅度随射流速度的增大而增大;表面波纹速度和波长沿传播方向增大,波纹间存在叠加现象.颗粒膜出现非轴对称振荡主要是因为喷嘴出口由气固不稳定性引起的射流流量脉动,射流流量脉动频率与撞击面振荡频率基本相当.     

Application of frequency-domain linearized Euler solutions to the prediction of aft fan tones and comparison with experimental measurements on model scale turbofan exhaust nozzles

Although it is widely accepted that aircraft noise needs to be further reduced, there is an equally important, on-going requirement to accurately predict the strengths of all the different aircraft noise sources, not only to ensure that a new aircraft is certifiable and can meet the ever more stringent local airport noise rules but also to prioritize and apply appropriate noise source reduction technologies at the design stage. As the bypass ratio of aircraft engines is increased - in order to reduce fuel consumption, emissions and jet mixing noise - the fan noise that radiates from the bypass exhaust nozzle is becoming one of the loudest engine sources, despite the large areas of acoustically absorptive treatment in the bypass duct. This paper addresses this ‘aft fan’ noise source, in particular the prediction of the propagation of fan noise through the bypass exhaust nozzle/jet exhaust flow and radiation out to the far-field observer. The proposed prediction method is equally applicable to fan tone and fan broadband noise (and also turbine and core noise) but here the method is validated with measured test data using simulated fan tones. The measured data had been previously acquired on two model scale turbofan engine exhausts with bypass and heated core flows typical of those found in a modern high bypass engine, but under static conditions (i.e. no flight simulation). The prediction method is based on frequency-domain solutions of the linearized Euler equations in conjunction with perfectly matched layer equations at the inlet and far-field boundaries using high-order finite differences. The discrete system of equations is inverted by the parallel sparse solver MUMPS. Far-field predictions are carried out by integrating Kirchhoff's formula in frequency domain. In addition to the acoustic modes excited and radiated, some non-acoustic waves within the cold stream-ambient shear layer are also captured by the computations at some flow and excitation frequencies. By extracting phase speed information from the near-field pressure solution, these non-acoustic waves are shown to be convective Kelvin-Helmholtz instability waves. Strouhal numbers computed along the shear layer, based on the local momentum thickness also confirm this in accordance with Michalke's instability criterion for incompressible round jets with a similar shear layer profile. Comparisons of the computed far-field results with the measured acoustic data reveal that, in general, the solver predicts the peak sound levels well when the farfield is dominated by the in-duct target mode (the target mode being the one specified to the in-duct mode generator). Calculations also show that the agreement can be considerably improved when the non-target modes are also included, despite their low in-duct levels. This is due to the fact that each duct mode has its own distinct directionality and a non-target low level mode may become dominant at angles where the higher-level target mode is directionally weak. The overall agreement between the computations and experiment strongly suggests that, at least for the range of mean flows and acoustic conditions considered, the physical aeroacoustic radiation processes are fully captured through the frequency-domain solutions to the linearized Euler equations and hence this could form the basis of a reliable aircraft noise prediction method.     

The connection between sound production and jet structure of the supersonic impinging jet

An experimental investigation into the sound-producing characteristics of moderately and highly underexpanded supersonic impinging jets exhausting from a round convergent nozzle is presented. The production of large plate tones by impingement on a square plate with a side dimension equal to 12 nozzle exit diameters is studied using random and phase-locked shadowgraph photography. Discrete frequency sound is produced in the near-wall region of the jet when a Mach disk occurs upstream of the standoff shock wave. Tones cease when the plate distance is approximately 2.2 free-jet cell lengths and the first and second shock waves are located in the free-jet positions. The production of impulsive sound appears to be associated with the collapse of the standoff shock wave during a portion of the oscillation cycle. Results from unsteady plate-pressure measurements indicate that plane-wave motion occurs in the impingement region and a secondary pressure maximum is observed on the plate adjacent to the flow region where sound appears to originate.     

Large-eddy simulation of impinging jets on smooth and rough surfaces

We performed large-eddy simulations (LES) of forced impinging jets over smooth and rough surfaces, containing large-scale, azimuthal vortices generated by the enhanced primary instability in the jet shear layer. The interaction between these vortices and the turbulence in the wall jet that is formed downstream of the impingement region determines their rate of decay. To explore the surface-roughness effects on the evolution of the vortices, sand-grain-like surfaces are generated using uniformly distributed but randomly oriented ellipsoids. The flow is compared to our previous LES of jets impinging on a smooth surface. In spite of the severe modification caused by the roughness on the near-wall flow, the vortex development is not significantly altered. Slightly faster decay of the primary vortices is observed in the rough-wall case compared to the smooth-wall one; the secondary vortex that detaches from the wall and is lifted up has larger vorticity. The highly disturbed near-wall flow is advected outward and affects the evolution of the primary vortex for a longer period during the vortex interaction. The robust turbulent generation mechanism in the outer shear layer, however, mitigates the changes in vortex behaviour. The momentum deficit and the enhancement of turbulence due to the surface roughness play a key role during this process.     

Experimental study of propagation of instability waves in a submerged jet under transverse acoustic excitation

An experimental study was conducted on the specific features of instability wave propagation in the mixing layer of a turbulent jet when the jet is excited by an external acoustic wave. We used the technique of conditional phase averaging of data obtained by particle image velocimetry using the reference signal of a microphone placed near the jet. The influence of the excitation frequency on the characteristics of large-scale structures in the mixing layer was investigated. It is shown that the propagation patterns of the instability waves agree well with previously obtained data on the localization of acoustic sources in turbulent jets.     

射流抛光多相紊流流场的数值模拟

 理论分析了射流抛光的紊动冲击射流特点,构建了射流抛光的垂直冲击射流模型和斜冲击射流模型。根据射流抛光冲击射流的特点,比较各种流体模型后,采用RNG k-e 模型应用于射流抛光模型的计算。利用计算流体力学理论的二阶迎风格式对抛光模型方程离散,用SIMPLEC数值计算方法对射流抛光过程的紊动冲击射流和离散相磨粒分布进行数值模拟,得到了射流抛光过程的连续流场和离散相磨粒与水溶液的耦合流场,同时计算出了抛光液射流在工件壁面上的压力、速度、紊动强度、剪切力分布和磨粒体积质量分布,分析了垂直射流抛光模型和斜冲击射流抛光模型紊流流场的特点。