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

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

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

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

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

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

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

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

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

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

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

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

INTERFERENCE OF SHOCK AND BOUNDARY LAYER IN SUPERSONIC FLOW

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

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

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

进气道等离子体/磁流体流动控制研究进展

为实现高速飞行器的宽速域飞行,如何保证进气道在非设计状态下的性能至关重要。相比于传统被动控制方式,等离子体/磁流体流动控制技术作为新概念主动流动控制技术,由于其具有结构简单,快速响应,并可根据实际飞行条件进行反馈控制等优势,在国内外上得到了广泛关注。本文介绍了等离子体/磁流体在高超/超声速进气道的主要应用方式与等离子体/磁流体建模方法。当进气道处于超临界状态时,等离子体/磁流体流动控制主要通过热阻塞效应产生虚拟型面,从而将激波系推回至唇口,该技术有望在需要短时间流动控制的高马赫数导弹上走向工程应用;由于等离子体/磁流体激励器与壁面平齐安装,对于高超声速飞行条件,相比于粗糙元其对热防护的要求较低,并且通过超声速风洞实验初步证明了通过高频激励对边界层施加扰动的可行性,需要从稳定性理论的角度对其物理机制进行研究。在后续发展中需要进一步创新等离子体产生技术及激励方式,发展等离子体与流的全耦合计算模型等离子体与流的全耦合计算模型与高效算法,为指导工程应用提供依据.     

进气道等离子体/磁流体流动控制研究进展

为实现高速飞行器的宽速域飞行,如何保证进气道在非设计状态下的性能至关重要。相比于传统被动控制方式,等离子体/磁流体流动控制技术作为新概念主动流动控制技术,由于其具有结构简单,快速响应,并可根据实际飞行条件进行反馈控制等优势,在国内外上得到了广泛关注。本文介绍了等离子体/磁流体在高超/超声速进气道的主要应用方式与等离子体/磁流体建模方法。当进气道处于超临界状态时,等离子体/磁流体流动控制主要通过热阻塞效应产生虚拟型面,从而将激波系推回至唇口,该技术有望在需要短时间流动控制的高马赫数导弹上走向工程应用;由于等离子体/磁流体激励器与壁面平齐安装,对于高超声速飞行条件,相比于粗糙元其对热防护的要求较低,并且通过超声速风洞实验初步证明了通过高频激励对边界层施加扰动的可行性,需要从稳定性理论的角度对其物理机制进行研究。在后续发展中需要进一步创新等离子体产生技术及激励方式,发展等离子体与流的全耦合计算模型等离子体与流的全耦合计算模型与高效算法 ,为指导工程应用提供依据.      

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.     

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.     

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

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

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.     

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 Reynolds number on boundary layer turbulence

A new facility for studying high Reynolds number incompressible turbulent boundary layer flows has been constructed. It consists of a moderately sized wind tunnel, completely enclosed by a pressure vessel, which can raise the ambient air pressure in and around the wind tunnel to 8 atmospheres. This results in a Reynolds number range of about 20:1, while maintaining incompressible flow. Results are presented for the zero pressure gradient flat plate boundary layer over a momentum thickness Reynolds number range 1500–15?000. Scaling issues for high Reynolds number non-equilibrium boundary layers are discussed, with data comparing the three-dimensional turbulent boundary layer flow over a swept bump at Reynolds numbers of 3800 and 8600. It is found that successful prediction of these types of flows must include length scales which do not scale on Reynolds number, but are inherent to the geometry of the 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.     

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