Numerical simulation of boundary-layer transition at Mach two

The linear and early nonlinear stages of boundary-layer transition at free-stream Mach numberM ==2.0 are investigated by direct numerical simulation of the compressible Navier-Stokes equations. Results from simulations with a large computational box and small-amplitude random initial conditions are compared with linear stability theory. The growth rates of oblique waves are reproduced correctly. Two-dimensional waves show a growth that is modulated in time, indicating the presence of an extra unstable mode which moves supersonically relative to the free stream. Further simulations are conducted to investigate the nonlinear development of two- and three-dimensional disturbances The transition due to oblique disturbance waves is the most likely cause of transition at this Mach number, and is found to lead to the development of strong streamwise vortices.     

基于直接数值模拟的可压缩湍流模型评估和改进

对来流Mach数2.25和6的平板边界层湍流进行了直接数值模拟, 并通过与理论、实验及他人计算结果的对比对数值结果进行了验证. 基于直接数值模拟得到的湍流数据库, 对常用的湍流模型进行了先验评估. 评估的湍流模型有k-εvarepsilon模型(包括标准k-εvarepsilon 模型、可实现的k-εvarepsilon模型及低Reynolds数k-εvarepsilon模型)、SA模型及BL模型. 结果显示, 对于Mach2.25的平板边界层, 可实现的k-εvarepsilon 模型及低Reynolds 数k-εvarepsilon模型具有较好的预测能力, 而标准k-εvarepsilon模型预测的湍流黏性系数偏高; SA模型在边界层内层预测准确度较高, 而在外层预测值偏高. 而对于Mach6的平板边界层, k-εvarepsilon模型及SA模型预测的湍流黏性系数均偏高, 尤其是标准k-εvarepsilon模型. 对于Mach6的平板边界层, BL模型低估了内-外层交界位置, 造成湍流黏性系数预测值严重偏低. 作者通过修改模型系数及内-外层交界位置对BL模型进行了修改, 修改后模型预测的湍流黏性系数与DNS给出的值吻合较好.     

Spatial direct numerical simulation of high-speed boundary-layer flows part I: Algorithmic considerations and validation

A highly accurate algorithm for the direct numerical simulation (DNS) of spatially evolving high-speed boundary-layer flows is described in detail and is carefully validated. To represent the evolution of instability waves faithfully, the fully explicit scheme relies on non-dissipative high-order compact-difference and spectral collocation methods. Several physical, mathematical, and practical issues relevant to the simulation of high-speed transitional flows are discussed. In particular, careful attention is paid to the implementation of inflow, outflow, and far-field boundary conditions. Four validation cases are presented, in which comparisons are made between DNS results and results obtained from either compressible linear stability theory or from the parabolized stability equation (PSE) method, the latter of which is valid for nonparallel flows and moderately nonlinear disturbance amplitudes. The first three test cases consider the propagation of two-dimensional second-mode disturbances in Mach 4.5 flat-plate boundary-layer flows. The final test case considers the evolution of a pair of oblique second-mode disturbances in a Mach 6.8 flow along a sharp cone. The agreement between the fundamentally different PSE and DNS approaches is remarkable for the test cases presented.     

超声速平面剪切层声辐射涡模态数值分析

对Mc = 1.2二维超声速空间发展平面自由剪切层, 进行了扰动模态及流动结构的数值分析. 采用时空三阶改进MacCormack格式, 差分求解可压缩扰动Navier-Stokes方程, 直接数值模拟入口不同基频谐波扰动的非线性演化特征. 采用空间线性稳定性理论证明, 计算所促发的扰动波是声辐射涡模态. 扰动参数及特征函数分析显示, 声辐射涡模态是弱色散的快／慢两种外部模态, 在扰动对流Mach数为超声速一侧呈膨胀／压缩状辐射. 单频受迫扰动可无相差地促发多模态混合扰动波, 而在自然扰动条件下, 剪切层的稳定性受慢模态主导.     

Direct numerical simulation of the transition to turbulence in a supersonic boundary layer on smooth and rough surfaces

Direct numerical simulations of instability development and transition to turbulence in a supersonic boundary layer on a flat plate are performed. The computations are carried out for moderate supersonic (free-stream Mach number M = 2) and hypersonic (M = 6) velocities. The boundary layer development is simulated, which includes the stages of linear growth of disturbances, their nonlinear interaction, stochastization, and turbulent flow formation. A laminar–turbulent transition initiated by distributed roughness of the plate surface at the Mach number M = 2 is also considered.     

Experimental study and direct numerical simulation of the evolution of disturbances in a viscous shock layer on a flat plate

The evolution of disturbances in a hypersonic viscous shock layer on a flat plate excited by slow-mode acoustic waves is considered numerically and experimentally. The parameters measured in the experiments performed with a free-stream Mach number M _∞ = 21 and Reynolds number Re _L = 1.44 · 10⁵ are the transverse profiles of the mean density and Mach number, the spectra of density fluctuations, and growth rates of natural disturbances. Direct numerical simulation of propagation of disturbances is performed by solving the Navier-Stokes equations with a high-order shock-capturing scheme. The numerical and experimental data characterizing the mean flow field, intensity of density fluctuations, and their growth rates are found to be in good agreement. Possible mechanisms of disturbance generation and evolution in the shock layer at hypersonic velocities are discussed. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 47, No. 5, pp. 3–15, September–October, 2006.     

Film cooling on a convex surface: influence of external pressure gradient and Mach number on film cooling performance

 The film cooling performance on a convex surface subjected to zero and favourable pressure gradient free-stream flow was investigated. Adiabatic film cooling effectiveness values were obtained for five different injection geometries, three with cylindrical holes and two with shaped holes. Heat transfer coefficients were derived for selected injection configurations. CO₂ was used as coolant to simulate density ratios between coolant and free-stream close to gas turbine engine conditions. The film cooling effectiveness results indicate a strong dependency on the free-stream Mach number level. Results obtained at the higher free-stream Mach number show for cylindrical holes generally and for shaped holes at moderate blowing rates significant higher film cooling effectiveness values compared to the lower free-stream Mach number data. Free-stream acceleration generally reduced adiabatic film cooling effectiveness relative to constant free-stream flow conditions. The different free-stream conditions investigated indicate no significant effects on the corresponding heat transfer increase due to film injection. The determined heat flux ratios or film cooling performance indicated that coolant injection with shaped film cooling holes is much more efficient than with cylindrical holes especially at higher blowing rates. Heat flux penalties can occur at high blowing rates when using cylindrical holes. Received on 29 May 2000     

Self-consistent parabolized stability equation (PSE) method for compressible boundary layer

The parabolized stability equation (PSE) method has been proven to be a useful and convenient tool for the investigation of the stability and transition problems of boundary layers. However, in its original formulation, for nonlinear problems, the complex wave number of each Fourier mode is determined by the so-called phase-locked rule, which results in non-self-consistency in the wave numbers. In this paper, a modification is proposed to make it self-consistent. The main idea is that, instead of allowing wave numbers to be complex, all wave numbers are kept real, and the growth or decay of each mode is simply manifested in the growth or decay of the modulus of its shape function. The validity of the new formulation is illustrated by comparing the results with those from the corresponding direct numerical simulation (DNS) as applied to a problem of compressible boundary layer with Mach number 6.     

Linear spatial instability analysis in 3D boundary layers using plane-marching 3D-LPSE

It is widely accepted that a robust and efficient method to compute the linear spatial amplified rate ought to be developed in three-dimensional (3D) boundary layers to predict the transition with the e ^N method, especially when the boundary layer varies significantly in the spanwise direction. The 3D-linear parabolized stability equation (3D-LPSE) approach, a 3D extension of the two-dimensional LPSE (2D-LPSE), is developed with a plane-marching procedure for investigating the instability of a 3D boundary layer with a significant spanwise variation. The method is suitable for a full Mach number region, and is validated by computing the unstable modes in 2D and 3D boundary layers, in both global and local instability problems. The predictions are in better agreement with the ones of the direct numerical simulation (DNS) rather than a 2D-eigenvalue problem (EVP) procedure. These results suggest that the plane-marching 3D-LPSE approach is a robust, efficient, and accurate choice for the local and global instability analysis in 2D and 3D boundary layers for all free-stream Mach numbers.     

Spatial direct numerical simulation of high-speed boundary-layer flows Part II: Transition on a cone in Mach 8 flow

The laminar breakdown of the boundary-layer flow of an axisymmetric sharp cone in a Mach 8 flow is simulated by a synergistic approach that combines the parabolized stability equation (PSE) method and spatial direct numerical simulation (DNS). The transitional state is triggered by a symmetric pair of oblique second-mode disturbances whose nonlinear interactions generate strong streamwise vorticity, which leads in turn to severe spanwise variations in the flow and eventual laminar breakdown. The PSE method is used to compute the weakly and moderately nonlinear initial stages of the transition process and, thereby, to derive a harmonically rich inflow condition for the DNS. The strongly nonlinear and laminar-breakdown stages of transition are then computed by well-resolved DNS, with a highly accurate algorithm that exploits spectral collocation and high-order compact-difference methods. Evolution of the flow is presented in terms of modal energies, mean quantities (e.g., skin friction), Reynolds stresses, turbulent kinetic energy, and flow visualization. The numerical test case is an approximate computational analog of one of the few stability experiments performed for hypersonic boundary-layer flows. Comparisons and contrasts are drawn between the experimental and the computational results. Rope-like waves similar to those observed in schlieren images of high-speed transitional flows are also observed in the numerical experiment and are shown to be visual manifestations of second-mode instability waves.This research was supported under NASA Contracts NAS1-19831 and NAS1-20059 for the first and second authors, respectively.     

Studies on stability and dynamics of a swirling jet

The temporal instability and nonlinear evolution of the swirling jet near a nozzle exit are studied by both normal-mode method and three dimensional direct numerical simulation (3D DNS). It is found that the swirl enhances the maximum linear growth rates for negative helical modes, while decreases the growth rate for axisymmetric mode. Numerical simulations show that the evolution in early stage is compared well with the linear stability theory. In nonlinear stage, the swirl promotes the breakup of 3D large scale organized structures in the flow into small eddies. The project supported by the National Natural Science Foundation of China (19772052)     

Effect of Heat Supply on the Stability of the Supersonic Swept Atiachment-Line Boundary Layer

The stability of a boundary layer with volume heat supply on the attachment line of a swept wing is investigated within the framework of the linear theory at supersonic inviscid-free-stream Mach numbers. The results of numerical calculations of the flow stability and neutral curves are presented for the flow on the leading edge of a swept wing with a swept angle χ=60° at various free-stream Mach numbers. The effect of volume heat supply on the characteristics of boundary layer stability on the attachment line is studied at a surface temperature equal to the temperature of the external inviscid flow. It is shown that in the case of a supersonic external inviscid flow volume heat supply may result in an increase in the critical Reynolds number and stabilization of disturbances corresponding to large wave numbers. For certain energy supply parameters the situation is reversed, the unstable disturbances corresponding to the main flow-instability zone are stabilized but another zone of flow-instability with small wave numbers and a significantly lower critical Reynolds number appears.     

A numerical study for small amplitude T-S waves in a supersonic boundary layer

ＩｎｔｒｏｄｕｃｔｉｏｎＲｅｃｅｎｔｌｙ ,ｔｈｅｐｒｏｂｌｅｍｏｆｔｒａｎｓｉｔｉｏｎｉｎｓｕｐｅｒｓｏｎｉｃｂｏｕｎｄａｒｙｌａｙｅｒｓｈａｓａｔｔｒａｃｔｅｄｍｏｒｅａｎｄｍｏｒｅａｔｔｅｎｔｉｏｎｄｕｅｔｏｉｔｓｔｅｃｈｎｉｃａｌｉｍｐｏｒｔａｎｃｅ .Ｕｐｔｏｎｏｗ ,ｍｏｓｔａｎａｌｙｓｅｓｗｅｒｅｌｉｎｅａｒａｎａｌｙｓｉｓ,ａｎｄｄｅｔａｉｌｅｄｅｘｐｅｒｉｍｅｎｔｓａｒｅｖｅｒｙｒａｒｅ .Ｓｏｍｅｐｅｏｐｌｅｈａｖｅｔｒｉｅｄｔｏｅｘｔｅｎｄｔｈｅｎｏｎｌｉｎｅａｒｔｈｅｏｒｙｏｆｈ…     

Direct Numerical Simulation of Hypersonic Boundary-Layer Flow on a Flared Cone

The forced transition of the boundary layer on an axisymmetric flared cone in Mach 6 flow is simulated by the method of spatial direct numerical simulation (DNS). The full effects of the flared afterbody are incorporated into the governing equations and boundary conditions; these effects include nonzero streamwise surface curvature, adverse streamwise pressure gradient, and decreasing boundary-layer edge Mach number. Transition is precipitated by periodic forcing at the computational inflow boundary with perturbations derived from parabolized stability equation (PSE) methodology and based, in part, on frequency spectra available from physical experiments. Significant qualitative differences are shown to exist between the present results and those obtained previously for a cone without afterbody flare. In both cases, the primary instability is of second-mode type; however, frequencies are much higher for the flared cone because of the decrease in boundary-layer thickness in the flared region. Moreover, Goertler modes, which are linearly stable for the straight cone, are unstable in regions of concave body flare. Reynolds stresses, which peak near the critical layer for the straight cone, exhibit peaks close to the wall for the flared cone. The cumulative effect appears to be that transition onset is shifted upstream for the flared cone. However, the length of the transition zone may possibly be greater because of the seemingly more gradual nature of the transition process on the flared cone. Received 20 March 1997 and accepted 21 May 1997     

Linear instability of the supersonic wake behind a flat plate aligned with a uniform stream

In this paper we present a theoretical and numerical study of the growth of linear disturbances in the high Reynolds number laminar compressible wake behind a flat plate which is aligned with a uniform stream. No ad hoc assumptions are made as to the nature of the undisturbed flow (in contrast to previous investigations) but instead the theory is developed rationally by use of proper wake profiles which satisfy the steady equations of motion. The initial growth of near-wake perturbations is governed by the compressible Rayleigh equation which is studied analytically for long and short waves. These solutions emphasize the asymptotic structures involved and provide a rational basis for a nonlinear development. The phenomenon of enhanced stability with increasing Mach number observed in compressible free shear-layers is demonstrated analytically for short- and long-wavelength disturbances. The evolution of arbitrary wavelength perturbations is addressed numerically and spatial stability solutions are presented that account for the relative importance of the different physical mechanisms present, such as three-dimensionality, increasing Mach numbers, and the nonparallel nature of the mean flow. Our findings indicate that for low enough (subsonic) Mach numbers, there exists a region of absolute instability very close to the trailing edge with the majority of the wake being convectively unstable. At higher Mach numbers (but still not large—hypersonic) the absolute instability region seems to disappear and the maximum available growth rates decrease considerably. Three-dimensional perturbations provide the highest spatial growth rates.This work was carried out while the author was a summer visitor at the Institute for Computer Applications in Science and Engineering, NASA Langley Research Center under NASA Contract No. NAS1-18605.     

Simulation and Modelling of Turbulent Trailing-Edge Flow

Computations of turbulent trailing-edge flow have been carried out at a Reynolds number of 1000 (based on the free-stream quantities and the trailing-edge thickness) using an unsteady 3D Reynolds-Averaged Navier–Stokes (URANS) code, in which two-equation (k–ε) turbulence models with various low-Re near wall treatments were implemented. Results from a direct numerical simulation (DNS) of the same flow are available for comparison and assessment of the turbulence models used in the URANS code. Two-dimensional URANS calculations are carried out with turbulence mean properties from the DNS used at the inlet; the inflow boundary-layer thickness is 6.42 times the trailing-edge thickness, close to typical turbine blade flow applications. Many of the key flow features observed in DNS are also predicted by the modelling; the flow oscillates in a similar way to that found in bluff-body flow with a von Kármán vortex street produced downstream. The recirculation bubble predicted by unsteady RANS has a similar shape to DNS, but with a length only half that of the DNS. It is found that the unsteadiness plays an important role in the near wake, comparable to the modelled turbulence, but that far downstream the modelled turbulence dominates. A spectral analysis applied to the force coefficient in the wall normal direction shows that a Strouhal number based on the trailing-edge thickness is 0.23, approximately twice that observed in DNS. To assess the modelling approximations, an a priori analysis has been applied using DNS data for the key individual terms in the turbulence model equations. A possible refinement to account for pressure transport is discussed. This revised version was published online in July 2006 with corrections to the Cover Date.     

8阶群速度控制格式及其应用

构造了8阶精度的群速度控制型差分格式.根据激波捕捉能力及对小尺度分辨能力对格式系数进行了优化.采用该格式对可压缩均匀各向同性湍流进行了直接数值模拟,所计算的最大湍流马赫数达到0.95.相同计算条件下的结果与他人结果吻合较好,说明了格式在捕捉激波的同时对湍流小尺度有较好的分辨能力.     

Optimized group velocity control scheme and DNS of decaying compressible turbulence of relative high turbulent Mach number

To overcome the difficulty in the DNS of compressible turbulence at high turbulent Mach number, a new difference scheme called GVC8 is developed. We have succeeded in the direct numerical simulation of decaying compressible turbulence up to turbulent Mach number 0.95. The statistical quantities thus obtained at lower turbulent Mach number agree well with those from previous authors with the same initial conditions, but they are limited to simulate at lower turbulent Mach numbers due to the so‐called start‐up problem. The energy spectrum and coherent structure of compressible turbulent flow are analysed. The scaling law of compressible turbulence is studied. The computed results indicate that the extended self‐similarity holds in decaying compressible turbulence despite the occurrence of shocklets, and compressibility has little effects on relative scaling exponents when turbulent Mach number is not very high. Copyright © 2005 John Wiley & Sons, Ltd.     

Spanwise effects on instabilities of compressible flow over a long rectangular cavity

The stability properties of two-dimensional (2D) and three-dimensional (3D) compressible flows over a rectangular cavity with length-to-depth ratio of \(L/D=6\) are analyzed at a free-stream Mach number of \(M_\infty =0.6\) and depth-based Reynolds number of \(Re_D=502\). In this study, we closely examine the influence of three-dimensionality on the wake mode that has been reported to exhibit high-amplitude fluctuations from the formation and ejection of large-scale spanwise vortices. Direct numerical simulation (DNS) and bi-global stability analysis are utilized to study the stability characteristics of the wake mode. Using the bi-global stability analysis with the time-averaged flow as the base state, we capture the global stability properties of the wake mode at a spanwise wavenumber of \(\beta =0\). To uncover spanwise effects on the 2D wake mode, 3D DNS are performed with cavity width-to-depth ratio of \(W/D=1\) and 2. We find that the 2D wake mode is not present in the 3D cavity flow with \(W/D=2\), in which spanwise structures are observed near the rear region of the cavity. These 3D instabilities are further investigated via bi-global stability analysis for spanwise wavelengths of \(\lambda /D=0.5{-}2.0\) to reveal the eigenspectra of the 3D eigenmodes. Based on the findings of 2D and 3D global stability analysis, we conclude that the absence of the wake mode in 3D rectangular cavity flows is due to the release of kinetic energy from the spanwise vortices to the streamwise vortical structures that develops from the spanwise instabilities.     

超声速平板边界层斜波失稳转捩过程研究

以5阶迎风和6阶对称紧致格式混合差分求解三维可压缩滤波Navier-Stokes方程,对Mach 数为4.5, Reynolds数为10000的空间发展平板边界层湍流进行了大涡模拟. 时间推进采用 紧致存储3阶Runge-Kutta方法,亚格子尺度模型为修正Smagorinsky涡黏性模型. 通过在 入口边界叠加一对线性最不稳定第一模态斜波扰动,数值模拟得到了平板层流边界层失稳转 捩直至湍流的演化过程. 对流场转捩过程中瞬时量及统计平均量的分析表明,数值模拟结果 与理论吻合,得到的Y型剪切层、交替\Lambda涡结构以及转捩后期的发卡涡结构的发展 变化与相关文献结果一致,湍流流谱定性合理.