Improved nonlinear parabolized stability equations approach for hypersonic boundary layers

The nonlinear parabolized stability equations(NPSEs) approach is widely used to study the evolution of disturbances in hypersonic boundary layers owing to its high computational efficiency. However, divergence of the NPSEs will occur when disturbances imposed at the inlet no longer play a leading role or when the nonlinear effect becomes very strong.Two major improvements are proposed here to deal with the divergence of the NPSEs. First, all disturbances are divided into two types: dominant waves and non-dominant waves. Disturbances imposed at the inlet or playing a leading role are defined as dominant waves, with all others being defined as non-dominant waves. Second, the streamwise wavenumbers of the non-dominant waves are obtained using the phase-locked method, while those of the dominant waves are obtained using an iterative method. Two reference wavenumbers are introduced in the phase-locked method, and methods for calculating them for different numbers of dominant waves are discussed. Direct numerical simulation(DNS) is performed to verify and validate the predictions of the improved NPSEs in a hypersonic boundary layer on an isothermal swept blunt plate.The results from the improved NPSEs approach are in good agreement with those of DNS, whereas the traditional NPSEs approach is subject to divergence, indicating that the improved NPSEs approach exhibits greater robustness.     

多孔表面推迟高超声速边界层转捩的机理

利用线性稳定性理论分析（LST）结合直接数值模拟（DNS）研究高超声速多孔表面边界层流动的失稳特征,分析多孔表面推迟高超声速边界层转捩的机理.在Ma=6,Re=2.0×10⁴（参考长度为入口处边界层位移厚度）条件下获得平板边界层及不同孔隙排列情形下平板边界层的典型流动特征,并采用LST方法分析光滑平板及多孔平板扰动的增长率及累计放大率.研究表明三维顺排及错排多孔表面都可以抑制第二模扰动的发展,推迟高超声速边界层转捩,但顺排多孔表面推迟高超声速边界层转捩能力更强.     

高超声速椭圆锥边界层横流转捩特性大涡模拟

横流效应显著影响高超声速飞行器的三维边界层转捩过程, 深化对该流动机制的认识有助于提升和改善飞行器气动性能及热力学环境. 针对HIFiRE5椭圆锥绕流问题, 采用大涡模拟方法计算分析了超声速边界层横流转捩特性, 并揭示其中的流动机理. 参考HIFiRE5风洞模型试验条件, 数值模拟中椭圆锥来流入口处施加人工速度扰动以激发边界层内不稳定扰动波, 进而预测了高超声速边界层流动横流失稳、转捩过程等基本流动特征, 并基于转捩热流分布形态对比, 获得了与试验数据基本吻合的计算结果. 研究发现, 椭圆锥中心线流动汇聚形成的流向涡结构非常容易失稳, 另外在中心线及侧缘之间的中部区域存在较强的横流不稳定性, 两种机制共同作用影响边界层转捩过程. 此外, 分析了来流扰动幅值对边界层横流失稳转捩的影响, 并发现静来流条件下, 横流区域出现两组独立的定常横流涡结构, 而强噪声来流条件下, 中心线主涡和中部横流涡均发生失稳转捩, 且在椭圆锥表面形成多峰状的转捩阵面. 最后, 深入分析流场的压力脉动动力学特性, 揭示了三维边界层发生失稳转捩的非线性演化机制.      

Stability characteristic of hypersonic flow over a blunt wedge under freestream pulse wave

To investigate the stability characteristic of hypersonic flow under the action of a freestream pulse wave, a high-order finite difference method was employed to do direction numerical simulation (DNS) of hypersonic unsteady flow over an 8° half-wedge-angle blunt wedge with freestream slow acoustic wave. The evolution of disturbance wave modes in the boundary layer under a pulse wave and a continuous wave are compared, and the wall temperature effect on the hypersonic boundary layer stability for a pulse wave disturbance is discussed. Results show that, both for a pulse wave and a continuous wave in freestream, the disturbance waves inside the nose boundary layer are mainly a fundamental mode; the Fourier amplitude of pressure disturbance mode in the boundary layer for a pulse wave is far less than that for a continuous wave, and the band frequency of the former is wider than that of the latter. All disturbance modes decay rapidly along the streamwise in the nose boundary layer. In the non-nose boundary layer, the dominant mode is transferred from fundamental mode into second harmonic. The transformation of dominant mode for a pulse wave appears much earlier than that for a continuous wave. Different frequency disturbance modes present different changes along streamline in the boundary layer, and the frequency band narrows around the second harmonic mode along the streamwise. Keen competition and the transformation of energy exist among different modes in the boundary layer. Wall temperature modifies the stability characteristic of the hypersonic boundary layer, which presents little effect on the development of fundamental modes and cooling wall could accelerates the growth of the high frequency mode as well as the dominant mode transformation.     

Receptivity to free-stream disturbance waves for hypersonic flow over a blunt cone

A high-order shock-fitting finite difference scheme is studied and used to do direction numerical simulation (DNS) of hypersonic unsteady flow over a blunt cone with fast acoustic waves in the free stream, and the receptivity problem in the blunt cone hypersonic boundary layers is studied. The results show that the acoustic waves are the strongest disturbance in the blunt cone hypersonic boundary layers. The wave modes of disturbance in the blunt cone boundary layers are first, second, and third modes which are generated and propagated downstream along the wall. The results also show that as the frequency decreases, the amplitudes of wave modes of disturbance increase, but there is a critical value. When frequency is over the critial value, the amplitudes decrease. Because of the discontinuity of curvature along the blunt cone body, the maximum amplitudes as a function of frequencies are not monotone. Supported by the National Natural Science Foundation of China (Grant Nos. 10632050 and 10502052)     

Transition prediction of a hypersonic boundary layer over a cone at small angle of attack—with the improvement of e<Superscript><Emphasis Type="Italic">N</Emphasis></Superscript> method

The problem of transition prediction for hypersonic boundary layers over a sharp cone has been studied in this work. The Mach number of the oncoming flow is 6, the cone half-angle is 5°, and the angle of attack is 1°. The conventional e ^N method is used, but the transition location so obtained is obviously incorrect. The reason is that in the conventional method, only the amplifying waves are taken into account, while in fact, for different meridians the decay processes of the disturbances before they begin to grow are different. Based on our own previous work, new interpretation and essential improvement for the e ^N method are proposed. Not only the amplification process but also the decay process is considered. The location, where by linear stability theory, the amplitude of disturbance wave is amplified from its initial small value to 1%, is considered to be the transition location. The new result for transition prediction thus obtained is found to be fairly satisfactory. It is also indicated that for the calculation of base flow, boundary layer equations can be used for a small angle of attack. Its computational cost is much smaller than those for DNS. Supported by the National Natural Science Foundation of China (Grant Nos. 10632050 and 90716007), the Special Foundation for the Authors of National Excellent Doctoral Dissertations (Grant No. 200328), and the Foundation of Liu-Hui Center of Applied Mathematics of Nankai University and Tianjin University     

A high-order cut-cell method for numerical simulation of hypersonic boundary-layer instability with surface roughness

Laminar-turbulent transition of hypersonic boundary layers can be affected significantly by the existence of surface roughness. Currently many important mechanisms of roughness-induced transition are not well understood. In recent years, direct numerical simulation (DNS) has been extensively applied for investigating instability and transition mechanisms of hypersonic boundary layers. Most of the past DNS studies, however, have been based on body-fitted grids for smooth surfaces without roughness. Due to complex geometry of embedded roughness, the use of body-fitted grids can be very difficult for flow with arbitrary surface roughness. In this paper, we present a new high-order cut-cell method to overcome the natural complexities in grid generation around arbitrary surface of roughness. The new method combines a non-uniform-grid finite-difference method for discrete grid points near the solid boundary and a shock-fitting method for the treatment of the bow shock. The non-uniform-grid finite-difference formulas are expressed in a general explicit form so that they can be applied to different multi-dimensional problems without any modification. The computational accuracy of new algorithms of up to O(h⁴) are tested on several one- and two-dimensional elliptic equations in irregular domains. In addition, the new method is applied to the simulation of the receptivity process of a Mach 5.92 flow over a flat plate under the combined effect of an isolated surface roughness element and surface blow and suction. A good agreement is found between the unsteady flow results and those obtained by a Linear Stability Theory (LST).     

Asymptotic stabilization of a nonlinear axially moving string by adaptive boundary control

In this paper, a robust adaptive boundary control for an axially moving string that shows nonlinear behavior resulting from spatially varying tension is investigated. A hydraulic actuator equipped with a damper is used as the control actuator at the right boundary of the string. The Lyapunov redesign method is employed to derive a robust control algorithm employing adaptation laws that estimate three unknown system parameters (mass per unit length of string, lumped mass of hydraulic actuator, and damping coefficient of damper) and an unknown boundary disturbance. The uniform asymptotic stability (when the three parameters are all unknown), the exponential stability (when they are known), and the uniform ultimate boundedness (with a bounded boundary disturbance) of the closed loop system are investigated. The convergence of the parameter estimates to the true values is shown. Numerical simulations are performed to demonstrate the effectiveness of the proposed robust adaptive boundary control.     

激光聚焦扰动作用下高超声速边界层稳定性实验研究

在马赫数6、单位雷诺数3.1×10⁶/m的条件下对半锥角7°直圆锥边界层稳定性开展了实验研究.以激光聚焦于流场中局部空间而产生的膨胀冲击波作为人工添加的小扰动,分析了该扰动对高超声速圆锥边界层流动稳定性的影响.实验中利用响应频率达到兆赫兹量级的高频压力传感器对圆锥壁面脉动压力进行测量,通过对压力数据进行短时傅里叶分析和功率谱分析发现,相比于不添加激光聚焦扰动的结果,添加激光聚焦扰动使边界层中第二模态波的出现位置提前,且扰动波的幅值大幅度地增加,在相同的流向范围内,激光聚焦扰动将边界层中的扰动波从线性发展阶段推进到非线性发展阶段,其对边界层中扰动波发展的促进效果明显.同时,激光聚焦位置的不同对边界层中扰动波的发展也具有不同的影响.当激光直接聚焦于圆锥壁面X=100 mm位置时,边界层中频率为90 kHz的扰动波幅值增长最快,在X=500 mm的位置处其幅值放大倍数为3.81,相比而言当激光聚焦位置位于圆锥前方自由来流中时,边界层幅值增长最快的扰动波频率大幅减小为73 kHz,相同范围内,其幅值放大倍数为4.51倍.由此可见,当激光聚焦位置位于圆锥上游的自由来流中时,其对边界层中扰动波的影响更为显著.     

Nonlinear interaction mechanisms of disturbances in supersonic flat-plate boundary layers

Due to the complexity of compressible flows, nonlinear hydrodynamic stability theories in supersonic boundary layers are not sufficient. In order to reveal the nonlinear interaction mechanisms of the rapidly amplified 3-D disturbances in supersonic boundary layers at high Mach numbers, the nonlinear evolutions of different disturbances in flat-plate boundary layers at Mach number 4.5, 6 and 8 are analyzed by numerical simulations. It can be concluded that the 3-D disturbances are amplified rapidly when the amplitude of the 2-D disturbance reaches a certain level. The most rapidly amplified 3-D disturbances are Klebanoff type (K-type) disturbances which have the same frequency as the 2-D disturbance. Among these K-type 3-D disturbances, the disturbances located at the junction of upper branch and lower branch of the neutral curve are amplified higher. Through analyzing the relationship between the amplification rate and the spanwise wavenumber of the 3-D disturbances at different evolution stages, the mechanism of the spanwise wavenumber selectivity of K-type 3-D disturbances in the presence of a finite amplitude 2-D disturbance is explained.     

Wavelet Cross-Spectrum Analysis of Multi-Scale Disturbance Instability and Transition on Sharp Cone Hypersonic Boundary Layer

Experimental measurement of hypersonic boundary layer stability and transition on a sharp cone with a half angle of 5° is carried out at free-coming stream Mach number 6 in a hypersonic wind tunnel. Mean andfluctuation surface-thermal-flux characteristics of the hypersonic boundary layer flow are measured by Pt-thin-film thermocouple temperature sensors installed at 28 stations on the cone surface along longitudinal direction. At hypersonic speeds, the dominant flow instabilities demonstrate that the growth rate of the second mode tends to exceed that of the low-frequency mode. Wavelet-based cross-spectrum technique is introduced to obtain the multi-scale cross-spectral characteristics of the fluctuating signals in thefrequency range of the second mode. Nonlinear interactions both of the second mode disturbance and the first mode disturbance are demonstrated to be dominant instabilities in the initial stage of laminar-turbulence transition for hypersonic shear flow.     

High-order non-uniform grid schemes for numerical simulation of hypersonic boundary-layer stability and transition

The direct numerical simulation of receptivity, instability and transition of hypersonic boundary layers requires high-order accurate schemes because lower-order schemes do not have an adequate accuracy level to compute the large range of time and length scales in such flow fields. The main limiting factor in the application of high-order schemes to practical boundary-layer flow problems is the numerical instability of high-order boundary closure schemes on the wall. This paper presents a family of high-order non-uniform grid finite difference schemes with stable boundary closures for the direct numerical simulation of hypersonic boundary-layer transition. By using an appropriate grid stretching, and clustering grid points near the boundary, high-order schemes with stable boundary closures can be obtained. The order of the schemes ranges from first-order at the lowest, to the global spectral collocation method at the highest. The accuracy and stability of the new high-order numerical schemes is tested by numerical simulations of the linear wave equation and two-dimensional incompressible flat plate boundary layer flows. The high-order non-uniform-grid schemes (up to the 11th-order) are subsequently applied for the simulation of the receptivity of a hypersonic boundary layer to free stream disturbances over a blunt leading edge. The steady and unsteady results show that the new high-order schemes are stable and are able to produce high accuracy for computations of the nonlinear two-dimensional Navier–Stokes equations for the wall bounded supersonic flow.     

基于温敏漆技术的圆锥高超声速大攻角绕流背风面流动结构实验研究

高超声速流动中, 大攻角下圆锥背风面边界层会存在流动分离与再附、边界层转捩等多种流动现象, 进而对圆锥表面温度分布产生显著的影响。为了对这一复杂流动规律及其对表面温升分布的影响进行讨论, 研究基于温敏漆技术, 得到了在Mach数为6的低湍流度来流条件下, 攻角为10°的圆锥背风面温升分布结果。通过对不同位置、不同方位角处温升分布曲线的分析, 对大攻角下圆锥背风面边界层流动发展过程及不同发展阶段的流动特征进行了讨论。同时, 通过对来流总压的调节, 得到了不同Reynolds数下的圆锥背风面温升分布结果, 总结了Reynolds数对流动的影响规律。研究发现, 高超声速大攻角圆锥背风面边界层流动发展过程中会依次出现层流分离、定常横流涡影响、转捩以及湍流分离与再附等流动特征, 而在不同的Reynolds数下, 各个流动特征产生影响的范围不同, 随着Reynolds数的降低, 层流范围和定常横流涡影响范围均有所增加, 而从观察到横流影响到转捩开始发生的范围基本相同。      

高精度非定常激波装配法

为正确模拟高超声速绕流中,来流小扰动与弓形激波之间的干扰对流动特征的影响,将弓形激波作为动边界,利用非定常特征关系处理激波处的边界条件.应用五阶精度迎风紧致格式和六阶精度的对称格式与三阶精度的R-K方法相结合,建立高精度非定常激波装配方法.采用该方法数值模拟钝锥高超声速定常流场和二维抛物外形高超声速边界层流动的感受性问题,数值模拟来流小扰动与弓形激波干扰激波后非定常扰动流场,研究扰动波进入边界层产生边界层不稳定波的特征.     

高超声速边界层中由粗糙元引起强制转捩的机理

由大粗糙元引起的高超声速边界层强制转捩在航天技术中有实际应用, 因而近年来受到人们的广泛关注.虽然目前导致该转捩过程的内在机理尚不完全清楚, 但有一点是明确的, 即粗糙元的尾迹流场中存在强对流不稳定性.文章的出发点是研究这种对流不稳定模态是如何触发转捩的.首先通过CFD方法, 计算出高超声速边界层中粗糙元的尾迹流场, 并对其进行二维稳定性分析.结果发现, 在传统不稳定Tollmien-Schlichting(T-S)模态出现的临界Reynolds数之前, 存在高增长率的无黏不稳定模态, 表现为对称的余弦模态和反对称的正弦模态.然后对该不稳定模态在粗糙元尾迹流中的演化进行了模拟, 验证了二维稳定性分析的结果, 并考察了非平行性效应的影响.最后通过直接数值模拟, 研究由这些不稳定模态触发转捩的全过程.结果表明, 对流不稳定模态确实是导致边界层转捩的关键机制.该转捩过程的特点是, 局部湍斑首先在不稳定模态特征函数的峰值附近出现, 然后向全流场扩散.就文章研究的工况而言, 余弦和正弦模态的相互作用对转捩的影响并不明显, 且后者在转捩过程中起主导作用.      

基于Ludwieg管的高超声速边界层转捩实验

高超声速边界层层/湍流转捩是高超声速飞行器气动力和气动热设计中的难点和热点问题.为了降低开展高超声速边界层不稳定性与转捩实验研究的门槛,研究基于Ludwieg管原理设计并建造了一座Mach 6高超声速管风洞,重点对Ludwieg管风洞的启动和运行过程开展了数值模拟,分析了储气段弯管布局对试验段流场的影响;之后,对该高超声速风洞的自由来流品质进行了静态和动态的标定,验证了风洞的设计Mach数,并给出了流场的动态扰动特征;最后,基于7°半张角尖锥标模开展了高超声速边界层转捩实验,通过表面齐平式安装的高频PCB传感器获得边界层不稳定波,分析了高超声速边界层不稳定波的演化特征.以上工作表明,Ludwieg管相对常规高超声速风洞具有建设和运行成本低、运行效率高、流场品质好等优点,适合开展高超声速边界层转捩等基础实验研究.      

Reduced aliasing formulations of the convective terms within the Navier–Stokes equations for a compressible fluid

The effect on aliasing errors of different formulations describing the cubically nonlinear convective terms within the discretized Navier–Stokes equations is examined in the presence of a non-trivial density spectrum. Fourier analysis shows that the existing skew-symmetric forms of the convective term result in reduced aliasing errors relative to the conservation form. Several formulations of the convective term, including a new formulation proposed for cubically nonlinear terms, are tested in direct numerical simulation (DNS) of decaying compressible isotropic turbulence both in chemically inert (small density fluctuations) and reactive cases (large density fluctuations) and for different degrees of resolution. In the DNS of reactive turbulent flow, the new cubic skew-symmetric form gives the most accurate results, consistent with the spectral error analysis, and at the lowest cost. In marginally resolved DNS and LES (poorly resolved by definition) the new cubic skew-symmetric form represents a robust convective formulation which minimizes both aliasing and computational cost while also allowing a reduction in the use of computationally expensive high-order dissipative filters.     

Stability of Large Amplitude Ekman-Hartmann Boundary Layers in MHD : The Case of Ill-Prepared Data

In this paper, we study an incompressible highly rotating fluid submitted to a high magnetic field between two planes with a Dirichlet boundary condition. We investigate the nonlinear stability of Ekman-Hartmann boundary layers under a spectral assumption for general initial data; this means that the data can be chosen as an arbitrary (but smooth enough) three-dimensional divergence free vector field independent of the small parameter.     

A numerical investigation of the evolution of 2-D disturbances in hypersonic boundary layers and the effect on the flow structure due to the existence of shocklets

The evolution of 2-D disturbances in hypersonic boundary layer with Mach number 6,8, and 10 was investigated numerically by three different numerical schemes. At the entrance, second mode T-S waves with different amplitudes were introduced, and the relation between the Mach number and the amplitude of the disturbance when shocklets started to appear was investigated. By comparing the disturbance velocity profiles with those provided by linear stability theory, the effects of shocklets on flow structures were also investigated.     

A DNS study on the flame structures and flame stabilization mechanism of a laboratory-scale lean premixed jet flame in crossflow

In the present work, direct numerical simulation (DNS) of a laboratory-scale lean premixed reacting jet flame in crossflow was performed to understand the flame structures and the flame stabilization mechanism. In the DNS, an ethylene-air jet with an equivalence ratio of 0.6 was injected into a hot vitiated crossflow. The jet Reynolds number reaches 6161. The DNS results were compared with those of the experiment with a good agreement. It was found that the windward and leeward branches of the flame show significantly different behaviors. The windward flame branch, appearing lifted and discontinuous, is located in the shear layer regions with high temperature, low vorticity and low scalar dissipation rate. The location of the peak heat release rate shifts to a higher mixture fraction with increasing distance from the jet exit. The leeward branch of the flame anchors in the shear layer near the jet exit. The recirculation zone in the wake of the jet facilitates the stabilization of the leeward flame. The chemical explosive mode analysis (CEMA) and species budget analysis were employed to characterize the local combustion mode. Auto-ignition plays a key role in the stabilization of the windward flame where a large range of extinction is also found due to the high strain rate. In contrast, premixed flame propagation is dominant on the leeward side.