飞行器大攻角复杂流动的POD和DMD对比分析

基于非结构/混合网格、耗散自适应2阶混合格式以及脱体涡模拟（detached eddy simulation,DES）方法开展了现代战斗机模型复杂分离流动的数值模拟,并与有限的平均气动力试验数据进行了对比,结果表明计算具有合理性,在此基础上进一步应用本征正交分解（proper orthogonal decomposition,POD）和动力学模态分解（dynamic mode decomposition,DMD）方法对数值模拟流场的非定常特性进行了对比分析.研究表明飞行器背风区流场由一对边条涡的螺旋运动主导,旋涡破裂前在横向空间截面上流场是中性稳定的,同时主涡核的运动是多频耦合的.POD和DMD的对比分析则表明:两者模态配对的方式不同,但主要模态之间具有一定相关性;POD模态中包含多种频率的运动,而且能量较集中于主模态,流场重构效率更高;DMD则将流场的主要特征运动提取为一些单频模态的组合,同时能够给出模态的稳定性.      

超燃火焰稳定凹腔自激振荡特性的数值研究

对超声速冷流条件下用于超燃冲压发动机的凹腔火焰稳定器的自激振荡特性进行研究.采用混合RANS/LES方法对非定常流场进行数值模拟,考虑了凹腔的长深比和后缘角度两个关键参数.混合RANS/LES方法很好的捕捉到流场非定常大尺度结构并揭示了凹腔自由剪切层的演化过程.对凹腔压力振荡历程进行幅频分析,所得到的频率和理论分析结果与文献计算结果符合的很好.结果表明,凹腔的长深比和后缘角度对凹腔自激振荡特性都有很大的影响.随着凹腔长深比的减小,振荡能量趋于集中到某些频率对应的振荡模式上.随着凹腔后缘倾角的减小,大部分频率对应的振荡很快的被削弱;相对于陡后缘凹腔,小角度后缘凹腔只有较高频率对应的振荡模式存在.     

二维圆柱绕流的动模态分析研究

运用Fluent对单/双圆柱绕流非定常流场进行数值模拟计算,分析了涡脱落的发展演变过程。采用DMD方法对圆柱绕流流场进行稳定性判断,利用分解得到的DMD模态数据对流场进行还原,并对还原流场的误差进行分析。分析结果表明:DMD方法得到的模态包含模态的频率信息;DMD方法能通过提取出来的主要结构完成对流场的还原及预测工作。对于复杂流动结构的DMD分解,所提取的数据量极大程度的影响了模态分解的结果,是能否精确还原实验流场的关键因素。     

低速开式空腔自激反馈流场结构与流致噪声的风洞试验研究

本文主要针对低速开式空腔流动自激振荡产生噪声问题,在0.55 m×0.4 m航空声学风洞开展了不同低马赫数(0.1/0.15/0.2/0.25)条件下长深比为2的空腔腔内流场结构和噪声特性风洞试验研究。通过利用高频粒子图像测速技术捕捉腔内流场结构,分析了腔内声波传递路径;完成空腔远场噪声和壁面压力测试,分析了噪声自激振荡模态和简正波模态,并对空腔壁面脉动压力和远场噪声进行压/声相关性研究。结果表明:空腔内部除主涡外,在腔口前缘处剪切涡与腔口后缘处碰撞涡明显存在;在875 Hz,1288 Hz,1875 Hz,2050 Hz四个频率附近出现了由声腔共振所致的单频噪声;壁面压力与远场噪声密切相关,在壁面压力主频位置有明显单频噪声出现。      

凹腔非定常特性的数值模拟

采用基于Menter SST两方程湍流模型的DES方法,数值模拟开式凹腔在跨声速条件下的非定常流动特性.计算凹腔底部和后壁面上的点的声压级频谱以及总声压级,证明在第二噪声模态上的声压级最大.     

基于大涡模拟下风力机翼型非定常转捩流动的动态模态分解分析

翼型绕翼流动对风力机整机性能产生重要影响.本文基于大涡模拟方法,得到风力机翼型非定常转捩流动.通过对压力流场的动态模态分解(DMD)分析,发现翼型层流分离泡的生成和发展为流动主要非定常特征,且该特征具有主要频率.预估得到高增长率T-S扰动波频率与DMD模态频率接近,发现高频变化的分离泡由边界层分离点下游不远处的T-S扰动波诱导K-H不稳定性而主导,以及中低频的DMD模态表现为T-S扰动波所引起的湍流结构.对比不同攻角下的DMD分解结果,发现上表面分离泡会逐渐向前移动,长度变短;而下表面分离泡略微后移,长度变化不大.     

基于凹腔驻涡的无焰燃烧室数值模拟

无焰燃烧的关键在于高温烟气与新鲜空气燃料的充分掺混。本文设计了一种基于凹腔驻涡的无焰燃烧室,并对其进行了零维和三维数值计算。零维数值计算运用简单的化学反应网络模型预测了燃烧室NO_x排放,凹腔当量比为0.5时,NO_x排放低于5 mg/m³（@15%O₂）。凹腔产生的高温烟气,可以使燃烧室后部温度分布均匀,消除了温度峰值,燃烧室实现了无焰燃烧。双凹腔结构有利于改善燃烧室出口温度分布。     

旋翼尾流时空发展二维速度场的实验研究

旋翼作为直升机的主要升力和操纵部件, 具有复杂的流场结构, 如非定常性, 桨-涡干扰和桨尖涡等, 导致旋翼流场研究十分困难. 针对这一问题, 结合锁相技术和粒子图像测速(particle image velocimetry, PIV)技术开展了悬停状态下旋翼流场的实验研究, 并通过本征正交分解(proper orthogonal decomposition, POD)提取主要含能模态, 刻画流场时空演化. 结果显示, 旋翼尾流发展过程中向旋转轴靠近, 二维流场结构呈现倒三角结构, 即扩展到三维流动中会呈现倒锥型结构的特性; 通过POD进行含能模态分析, 旋翼尾流中对湍动能贡献最大的为桨叶涡结构, 其次是桨尖涡结构.      

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

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

风力机翼型在失速工况下非定常流场的本征正交分解分析

风力机气动性能受静态失速与动态失速影响很大,对风力机翼型的失速问题研究具有重要意义。本文通过计算流体力学方法得到的风力机翼型在固定大攻角工况,以及大攻角震荡工况下的非定常流场,来研究翼型静态失速与动态失速。采用本征正交分解方法(POD),对非定常流场降阶,得到流场的POD模态以及对应的系数。POD模态结果表明在静态失速下,主要非定常流动结构是尾迹区域交替脱落的涡结构;在动态失速下,除了尾迹区域,前缘和整个吸力面都存在流动分离结构。     

Self-sustained oscillations of turbulent pipe flow terminated by an axisymmetric cavity

Turbulent flow through a long pipe terminated by an axisymmetric cavity can give rise to self-sustained oscillations exhibiting a very strong coherence, as evidenced by the narrow-band character of corresponding amplitude spectra. These oscillations, associated with the turbulent axisymmetric jet passing through the cavity, are strongly influenced by the acoustic modes of the pipe. The frequencies of oscillation lie within or near the range of most “unstable” frequencies of the turbulent jet previously predicted by using concepts of inviscid hydrodynamic stability theory; consequently, these experiments show truly self-excited and strongly coherent “instability” of a fully turbulent, low Mach number (～10^?2), axisymmetric flow undergoing separation, corroborating previous experiments involving the external forcing of free turbulent jets. As flow velocity or cavity length is varied, both upward and downward jumps in oscillation frequency are observed; the sign (up or down) of these jumps tends to systematically alternate with increase of velocity or length. The role of these frequency jumps is, in effect, to allow the oscillation of the flow to remain “locked-on” to a pipe mode over a wide range of impingement length or flow velocity. Moreover, these jumps exhibit two types of behavior: for the first kind, the predominant frequency makes a relatively continuous transition between stages and the frequency of the neighboring stage appears as a secondary component; for the second kind, there is a dead zone (where no oscillation occurs) between stages. The consequence of externally exciting the system is strongly dependent on whether the self-sustaining oscillation is relatively near, or well away from, a frequency jump. During excitation, the amplitudes of pressure fluctuations in the cavity substantially exceed the corresponding no-flow values only in regions away from the frequency jumps; at locations of jumps, there can be significant attenuation of the no-flow excitation amplitude. For the type of frequency jump involving a “dead zone”, enhancement of a given mode of oscillation can be achieved by externally exciting not only the given mode, but also neighboring modes. For the other type of jump, involving a relatively continuous transition from one stage to the next, the predominant mode of oscillation following the jump is that mode giving maximum amplitude response to excitation before the jump.     

ANALYSIS OF THE PERIODIC PRESSURE FLUCTUATIONS INDUCED BY FLOW OVER A CAVITY

Subsonic flows over Helmholtz resonators often cause strong periodic pressure fluctuations inside the resonators over a range of outer flow velocities. The flow-excitation mechanism is known to be governed by both the shedding of discrete vortices within the shear layer over the orifice and the acoustic response of the cavity. This self-sustained oscillation phenomenon is often analyzed by using a feedback loop model where the flow excitation and the acoustic response of the resonator are approximately modelled as a forward gain function and as a backward gain function respectively. In the present work, a similar approach was followed and a new forward gain function was derived based on the concept of “vortex sound” to model the flow excitation. The formulation combined this forward gain function with a backward gain function from previous work, within the framework of the feedback loop analysis. The approximate method allowed the frequency and the relative amplitude of the cavity pressure fluctuations to be predicted for a range of flow velocities. In addition, the extended Nyquist stability criterion was used to estimate the onset and the termination velocities of the first two modes of the shear layer flow oscillations. Experimental data were obtained using a rigid-walled cavity in a low-speed wind tunnel. The results showed that the model predictions were in reasonably good agreement with the experimental data.     

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.     

Low-frequency components and modulation processes in compressible cavity flows

The modulated character of cavity-flow oscillations is investigated through a so-called modulation analysis of spectral distributions. The approach is based on a recently proposed viewpoint on the Rossiter formula and concerns mid- to high-subsonic flows in shallow cavities. For this type of configuration, the spectra are mainly characterized by the presence of several dominant peaks (Rossiter modes) that are not in a harmonic relation but uniformly spaced at a distance equal to the fundamental frequency of the oscillation mechanism (aero-acoustic feedback loop). This feature is interpreted as the result of an amplitude modulation process and related to variations in the vortex-corner interaction in the downstream part of the cavity (γ modulation). A lower frequency modulation is identified through the secondary peak distribution. A detailed analysis of the spectral structure confirms the presence of a component at the corresponding frequency value (Δ_f mode). The assumption of a specific coupling between the two modulation processes is investigated. It leads to a new approximate form for the γ-modulation ratio that allows an explicit expression of the Rossiter constant γ related to aspect ratio of the cavity (L/D).     

Noisy inflows cause a shedding-mode switching in flow past an oscillating cylinder

Vortex streets formed behind oscillating bluff bodies consist of arrays of groups of two, three, or four vortices classified as 2S, P+S, and 2P shedding modes, respectively. The prevailing dominant mode depends primarily on the amplitude and the frequency of the oscillation and on the Reynolds number. We investigate the effect of noise at the inflow on the stability of these vortex modes in laminar flow past a circular cylinder. We employ stochastic simulations based on a new polynomial chaos method to study the shedding-mode switching from a P+S pattern to a 2S mode in the presence of noise.     

Instability and mode transition analysis of a hydrogen-rich combustion in a model afterburner

Combustion instabilities were investigated experimentally for a hydrogen-rich combustion in a model afterburner installed at the end of a high-enthalpy wind tunnel. Air was supplied at 0.3 MPa and 950 K. The combustion instabilities were studied with the time-resolved measurements of a near-infrared (NIR) emission from water molecules over 780 nm using a high-speed video camera. Pressure was also measured in the combustor. The pressure and the NIR images were analyzed by data-driven approach, which include the fast Fourier transform (FFT), the wavelet transform, the dynamic mode decomposition (DMD) and the Gaussian process latent variable methods (GP-LVM). Thermoacoustic instability was observed under a rich condition, and the amplitude of the pressure oscillation was the maximum at the overall equivalence ratio of approximately 2.4 or 2.7 as a result of the FFT. The combustion dynamics were investigated in detail for an experimental run at the equivalence ratio of 2.4. A pressure spectrogram indicated a flame–vortex interaction with a Strouhal number of 0.5 (2300 Hz), thermoacoustic instability (560 Hz), and their transitions with the wavelet transform. For NIR images, the same tendency was also observed in the spectrogram of the modes obtained by the Gabor-filtered DMD, which could clearly resolve the high-order harmonic modes of the flame–vortex interaction and the thermoacoustic instability. Furthermore, NIR images were analyzed with GP-LVM to study the evolution of the combustion dynamics in a three-dimensional latent space. Recurrence plots with the Euclidean distance function were used to visualize the evolutions of the combustion dynamics. A limit cycle behavior of the flame–vortex interaction was clearly observed, whereas the limit cycle of the thermoacoustic instability showed more complicated behaviors. The transition behaviors of the instabilities were observed in the recurrence plots in detail, indicating that the flame–vortex interaction excited the fourth harmonic mode of the thermoacoustic instability, followed by the basic mode.     

Experimental investigation on laminar separation control for flow over a two-dimensional bump

The laminar boundary layer separation flow over a two-dimensional bump controlled by synthetic jets is experimentally investigated in a water channel with hydrogen-bubble visualisation and particle image velocimetry (PIV) techniques. The two-dimensional synthetic jet is applied near the separation point. Two Reynolds numbers (Re = 700 and 1120) based on the bump height and free-stream velocity are adopted in this experiment, and seven different excitation frequencies at each Reynolds number are considered, focusing on the separation control as well as the vortex dynamics. The experimental results show that the optimal control can only be achieved within some excitation frequencies at both Reynolds numbers. However, beyond this range, further increasing the excitation frequency leads to an increase in the separation region. The proper orthogonal decomposition (POD) technique and vortex identification by swirling strength (Λ_ci) are applied for the deeper analysis of the separated flow. The reconstructed Λ_ci field by the first four POD modes is used and vortex lock-on phenomenon is observed. It is found that the negative synthetic jet vortex with clockwise rotation draws the separated wake shear layer as it is convected downstream, and then they syncretise together. Thus, the new vortex is induced and shedding downstream periodically.