超声速边界层/混合层组合流动的稳定性分析简

利用可压缩线性稳定性理论研究了超声速混合层考虑壁面影响流动时的失稳特性. 基本流场选取了具有不同速度特征的2 股均匀来流,进入存在上下壁面的流道中. 混合层与边界层的距离为1~3 个边界层厚度,其中壁面取为绝热壁. 分析了该流动在超声速情况下的稳定性特征,同时还讨论了不同波角下的三维扰动波的演化特点,并与二维扰动波进行了比较和分析. 研究结果表明,在此流动情况下,边界层流动和混合层流动的稳定性特征同时存在,并互有影响,其流动稳定性特征既有别于单纯的平板边界层,也有别于单纯的平面混合层,呈现出了新的稳定性特征.     

超声速混合层中扰动增强混合实验

以基于纳米技术的平面激光散射(nano-based planar laser scattering, NPLS)流动显示技术定性研究了隔板扰动对超声速混合层($Mc=0.5$)的混合增强效果. 首先通过系列实验优化设计了扰动参数. 实验结果表明,超声速混合层对于从隔板引入的扰动非常敏感. 二维扰动的混合强化机制是提前混合层失稳位置,增厚混合层;而三维扰动的混合强化机制主要是通过诱导流向涡和展向运动,促进流动三维性质的发展. 总体而言,三维扰动的混合强化效果优于二维扰动. 由于是超声速混合层,隔板上的扰动片虽然很薄,但同样会引起激波的产生,是该方法中总压损失的主要原因.      

采用不同黏性处理方法的宽无叶扩压器不稳定流动研究

径向无叶扩压器的全局稳定性可能受到核心主流失稳,出口回流与壁面边界层分离等因素影响,对于宽无叶扩压器,无黏核心主流与壁面边界层流动对不稳定扰动诱发的作用机理是当前研究的重点.本文首先通过数值计算获得了大宽度比孤立无叶扩压器平均流动,然后基于小扰动理论和周向均质假设,分别对欧拉方程与 Navier-Stokes 方程进行线性化,建立了基于无黏核心流动的稳定性分析方法,以及基于涡黏性与分子黏性的混合稳定性分析方法;通过与实验结果的对比,验证了混合稳定性分析方法预测所得流动失稳频率和全局直接模态的准确性;最后基于伴随方法获得了特征值的结构敏感性,揭示了不同黏性处理条件下宽无叶扩压器内全局不失稳扰动的源发区域.在只考虑核心主流的无黏条件下,宽无叶扩压器内流动不稳定扰动来源于流场中部,为二维的离心失稳;在同时考虑核心主流与边界层的作用时,宽无叶扩压器不稳定扰动不仅来源于扩压器流场中部的核心主流,壁面回流对于不稳定扰动的产生了重要影响.      

采用不同黏性处理方法的宽无叶扩压器不稳定流动研究

径向无叶扩压器的全局稳定性可能受到核心主流失稳,出口回流与壁面边界层分离等因素影响,对于宽无叶扩压器,无黏核心主流与壁面边界层流动对不稳定扰动诱发的作用机理是当前研究的重点.本文首先通过数值计算获得了大宽度比孤立无叶扩压器平均流动,然后基于小扰动理论和周向均质假设,分别对欧拉方程与Navier-Stokes方程进行线性化,建立了基于无黏核心流动的稳定性分析方法,以及基于涡黏性与分子黏性的混合稳定性分析方法;通过与实验结果的对比,验证了混合稳定性分析方法预测所得流动失稳频率和全局直接模态的准确性;最后基于伴随方法获得了特征值的结构敏感性,揭示了不同黏性处理条件下宽无叶扩压器内全局不失稳扰动的源发区域.在只考虑核心主流的无黏条件下,宽无叶扩压器内流动不稳定扰动来源于流场中部,为二维的离心失稳;在同时考虑核心主流与边界层的作用时,宽无叶扩压器不稳定扰动不仅来源于扩压器流场中部的核心主流,壁面回流对于不稳定扰动的产生了重要影响.     

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

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

纳米示踪平面激光散射技术在激波复杂流场测量中的应用

在激波以及激波边界层相互作用这类含激波的复杂流场中,流场结构具有明显的三维特征.研究这类流场,采用纹影、阴影和干涉等传统流动显示技术空间分辨率较低,难以分辨流场的三维特性.基于纳米示踪的平面激光散射技术(nano-tracer planar laser scattering,NPLS),是作者近年来开发的一种新的研究超声速流场的测试与显示技术,可对超声速复杂三维流场进行高时空分辨率流动显示与测量.NPLS技术的特点使其成为测量激波复杂流场的有力手段.近年来,作者以NPLS技术为主要手段,对航空航天领域典型的激波复杂流场进行了试验研究,包括超声速弹头绕流、超声速混合层、超声速边界层,以及激波边界层相互作用流场,显示出NPLS技术在激波复杂流场精细测试与流动显示中优势.本文简要介绍NPLS技术在激波复杂流场测量中应用的研究进展.      

在超声速流动中激波与边界层的干扰特性

激波与物面边界层的干扰涉及可压缩流动的稳定性、转捩、分离等问题,直接影响到飞行器的阻力、表面热防护和飞行性能等工程技术问题。首先总结了前人对于激波与边界层的干扰所做的工作,之后重点研究和对比分析了超声速与跨声速流动中,正激波、斜激波以及头部激波对于飞行器层流和湍流边界层的干扰影响。激波强度的不同对边界层干扰作用不同,在强干扰情况下将会引起边界层分离和翼型失速。     

INTERFERENCE OF SHOCK AND BOUNDARY LAYER IN SUPERSONIC FLOW

激波与物面边界层的干扰涉及可压缩流动的稳定性、转捩、分离等问题,直接影响到飞行器的阻力、表面热防护和飞行性能等工程技术问题。首先总结了前人对于激波与边界层的干扰所做的工作,之后重点研究和对比分析了超声速与跨声速流动中,正激波、斜激波以及头部激波对于飞行器层流和湍流边界层的干扰影响。激波强度的不同对边界层干扰作用不同,在强干扰情况下将会引起边界层分离和翼型失速。     

高超声速粗糙元诱导转捩的数值模拟及机理分析

采用直接数值模拟方法细致刻画了钻石型粗糙元诱导的高超声速边界层从层流到湍流的转捩过程,从拓扑结构稳定性和边界层流动稳定性两个角度分析了钻石型粗糙元诱导转捩的机理. 流动结构的拓扑分析表明,钻石型粗糙元头部区域和底部区域分别存在不稳定的鞍点-鞍点(SS) 型轨线和鞍点-结点-鞍点(SNS) 型轨线,在扰动的作用下其会形成非定常、非对称的振荡结构. 边界层流动失稳过程计算分析表明,钻石型粗糙元会产生高波数扰动,并发现在扰动发展过程中大尺度结构会破碎. 两种不同类型的流动失稳效应同时存在. 此外,通过不同类型粗糙元(圆柱、斜坡及钻石型) 的对比,揭示了不同类型粗糙元诱导转捩机理的差异,为高超声速人工转捩装置设计提供了基础理论支撑.      

微槽-吸气组合控制平板边界层多模态稳定性研究

边界层转捩会使高超声速飞行器壁面摩阻和热流显著增加,因此在高超声速飞行器设计过程中往往占据重要地位.针对高超声速飞行器多模态转捩控制问题,提出了微槽道(1 mm)与边界层吸气的组合控制方法,并通过直接数值模拟和线性稳定性理论研究了Ma=4.5平板边界层的稳定性及组合控制效果.边界层在无控状态时,同时存在失稳的第一、二模态波,且二维第二模态波最不稳定;单纯施加微槽道控制时,边界层第二模态波会被抑制但第一模态波会被略微激发.对比而言,采用“微槽-吸气”组合控制后,不仅增强了对第二模态波的抑制效果,而且减弱了第一模态波的激发程度;同时随着吸气强度的增加,第二模态波不稳定区域明显收缩、频率显著增高,而第一模态波则变化不明显.相较于单纯的微槽道,吸气增强了“微槽吸收”与“声波散射”作用,因此中等吸气强度下该组合控制方法对第一和第二模态波的增长率分别实现了12.63%和28.02%的抑制效果.以上结果表明“微槽-吸气”组合控制手段具有适用宽频、布置区域灵活的优点,展现出了一定的多模态控制效果.     

On the inviscid acoustic-mode instability of supersonic shear flows

The same methods used previously to study acoustic-mode instability in supersonic boundary layers are applied to free shear layers, and new calculations are made for boundary layers with cooling and suction. The objective is to obtain additional information about acoustic-mode instability, and to find what features of the instability are common to boundary layers and free shear flows. Acoustic modes exist whenever there is an embedded region of locally supersonic flow relative to the phase speed of the instability wave. Consequently, they can be found in boundary layers, wakes, and jets, but not in mixing layers unless the flow is confined. In this first part of a two-part paper, attention is directed principally to two-dimensional waves. The linear, inviscid stability theory is used to calculate spatial amplification rates at Mach number 3 for the sinuous and varicose modes of a single wake flow and a single jet flow, each made up of the same mixing-layer profile plus a central region of uniform flow. Along with sequences of sinuous and varicose unstable modes clearly identifiable as acoustic modes, both of these flows, unlike the boundary layer, have a lowest sinuous mode that is the most unstable. The unstable modes include both subsonic and radiating disturbances with large amplification rates. The latter phenomenon is also found for highly cooled boundary layers with suction. In these boundary layers, suction is generally stabilizing for nonradiating acoustic disturbances, but destabilizing for radiating disturbances.The work described in this paper was carried out at the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration (NASA). Support from the Aerodynamics Division of the Office of Aeronautics and Exploration Technology is gratefully acknowledged. A preliminary version of this paper was presented at the Fourth Symposium on Numerical and Physical Aspects of Aerodynamic Flows, California State University, Long Beach, CA, 16–19 January 1989.     

The effect of inflow conditions on the transition to turbulence in large eddy simulations of spatially developing mixing layers

Large Eddy Simulations (LES) of spatially developing turbulent mixing layers have been performed for flows of uniform density and Reynolds numbers of up to 50,000 based on the visual thickness of the layer and the velocity difference across it. On a fine LES grid, a validation simulation performed with a hyperbolic tangent inflow profile produces flow statistics that compare extremely well with reference Direct Numerical Simulation (DNS) data. An inflow profile derived from laminar Blasius profiles produces a flow that is significantly different to the reference DNS, particularly with respect to the initial development of the flow. When compared with experimental data, however, it is the boundary layer-type inflow simulation produces the better prediction of the flow statistics, including the mean transition location. It is found that the boundary layer inflow condition is more unstable than the hyperbolic tangent inlet profile. A suitably designed coarse LES grid produces good predictions of the mean transition location with boundary layer inflow conditions at a low computational cost. The results suggest that hyperbolic tangent functions may produce unreliable DNS data when used as the initial condition for studies of the transition in the mixing layer flow.     

The effect of boundary layer fluctuations on the streamwise vortex structure in simulated plane turbulent mixing layers

This paper details the influence of the magnitude of imposed inflow fluctuations on Large Eddy Simulations of a spatially developing turbulent mixing layer originating from laminar boundary layers. The fluctuations are physically-correlated, and produced by an inflow generation technique. The imposed high-speed side boundary layer fluctuation magnitude is varied from a low-level, up to a magnitude sufficiently high that the boundary layer can be considered, in a mean sense, as nominally laminar. Cross-plane flow visualisation shows that each simulation contains streamwise vortices in the laminar and turbulent regions of the mixing layer. Statistical analysis of the secondary shear stress reveals that mixing layers originating from boundary layers with low-level fluctuations contain a spatially stationary streamwise structure. Increasing the high-speed side boundary layer fluctuation magnitude leads to a weakening of this stationary streamwise structure, or its removal from the flow entirely. The mixing layer growth rate reduces with increasing initial fluctuation level. These findings are discussed in terms of the available experimental data on mixing layers, and recommendations for both future experimental and numerical research into the mixing layer are made.     

Direct numerical simulation of stable and unstable turbulent thermal boundary layers

This paper presents direct numerical simulations (DNS) of stable and unstable turbulent thermal boundary layers. Since a buoyancy-affected boundary layer is often encountered in an urban environmental space where stable and unstable stratifications exist, exploring a buoyancy-affected boundary layer is very important to know the transport phenomena of the flow in an urban space. Although actual observation may qualitatively provide the characteristics of these flows, the relevant quantitative turbulent quantities are very difficult to measure. Thus, in order to quantitatively investigate a buoyancy-affected boundary layer in detail, we have here carried out for the first time time- and space-developing DNS of slightly stable and unstable turbulent thermal boundary layers. The DNS results show the quantitative turbulent statistics and structures of stable and unstable thermal boundary layers, in which the characteristic transport phenomena of thermally stratified boundary layers are demonstrated by indicating the budgets of turbulent shear stress and turbulent heat flux. Even though the input of buoyant force is not large, the influence of buoyancy is clearly revealed in both stable and unstable turbulent boundary layers. In particular, it is found that both stable and unstable thermal stratifications caused by the weak buoyant force remarkably alter the structure of near-wall turbulence.     

The structure and dynamics of reacting plane mixing layers

The reacting two-dimensional plane mixing layer has been studied in two configurations: a rearward facing step and a two-stream mixing layer. Observations have been made of the steady state behavior, and the unsteady behavior when the flow was forced by a specific acoustic frequency. The steady behavior of the mean properties of the reacting flows is similar to that of non-reacting free shear flows except for the global effects of thermodynamic property changes. The structure of these flows is qualitatively similar to that of non-reacting flows. Vortices form by the two-dimensional Kelvin-Helmholtz instability and grow by subharmonic combination until the mixing layer interacts with the walls. Entrainment is dominated by the two-dimensional vortex motion. Three-dimensional instabilities give rise to secondary vortices which are coherent over several Kelvin-Helmholtz structures and dominate the fine scale mixing process. The mixing transition corresponds to a loss of coherence in the layer. Unsteady behavior occurs when there are resonant interactions with the Kelvin-Helmholtz instability or the instability associated with the recirculation vortex in the rearward facing step flow. Modeling efforts are reported which show promise of simulating the essential features of plane mixing layers.A version of this paper was presented at the ASME Winter Annual Meeting of 1984 and printed in AMD-Vol. 66     

可压缩边界层与混合层失稳结构的研究进展及其工程应用

本文回顾了可压缩边界层与混合层中失稳结构及其应用的研究进展. 这些工作包括人们对高超声速平板边界层失稳特性、高超声速圆锥边界层转捩攻角效应的产生机制和可压缩混合层失稳特性的研究认识,以及这些相关认识的3 个应用方向. 这些相关工作丰富了人们对高速流动转捩和湍流拟序结构的认识.     

Near-field turbulence properties of single- and two-stream plane mixing layers

Detailed mean flow and turbulence measurements have been made in the near-field of two plane mixing layers in air with a maximum velocity of 21 m/s. The experimental rig enabled mixing layers of velocity ratios 0 and 0.46 to be generated simultaneously. Cases with both tripped and untripped initial boundary layers were studied. In all cases, it was found that the two-stream layer developed to the self-preserving state in a distance much shorter than the single-stream layer, which followed accepted criteria for the development distance. The asymptotic levels of the turbulence quantities in the two-stream layer and the development of the single-stream layer showed agreement with existing data. The results suggest that the two-stream mixing layer should provide a better test case for the development of turbulence models and calculation methods than the single-stream mixing layer.     

Numerical solution of the problem of the mixing of the boundary layers shed from the trailing edge of a wing

During the mixing of viscous incompressible flows with different velocities, in the vicinity of a trailing edge an interaction region with a three-layer structure is formed, similar to that in the case of symmetric shedding with equal velocities. The boundary layers developing on the upper and lower sides of the airfoil form a viscous mixing layer, or vortex sheet, which separates the flows downstream of the trailing edge. The boundary value problem corresponding to the flow in the viscous sublayer in the vicinity of the trailing edge of a flat plate is solved for high Reynolds numbers using an efficient numerical method for solving the equations of asymptotic interaction theory.     

Numerical Simulation of Particle Dispersion in a Spatially Developing Mixing Layer

Although there have been several numerical studies on particle dispersion in mixing layers, most of them have been conducted for temporally evolving mixing layers. In this study, numerical simulations of a spatially developing mixing layer are performed to investigate particle dispersion under various conditions. The full compressible Navier--Stokes equations are solved with a high-order compact finite difference scheme, along with high-order time-integration. Accurate non-reflecting boundary conditions for the fluid flow are used, and several methods for introducing particles into the computational domain are tested. The particles are traced using a Lagrangian approach assuming one-way coupling between the continuous and the dispersed phases. The study focuses on the roles of the large-scale vortex structures in particle dispersion at low, medium and high Stokes numbers, which highlights the important effects of interacting vortex structures in nearby regions in the spatially developing mixing layer. The effects of particles with randomly distributed sizes (or Stokes numbers) are also investigated. Both instantaneous flow fields and statistical quantities are analyzed, which reveals essential features of particle dispersion in spatially developing free shear flows, which are different from those observed in temporally developing flows. The inclusion of the gravity not only modifies the overall dispersion patterns, but also enhances stream-crossing by particles. Received 7 June 2001 and accepted 19 February 2002