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

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

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

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

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

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

基于NPLS技术的可压缩湍流机理实验研究新进展

可压缩湍流机理的实验研究是一件难度很大的工作, 其主要的难度在于高时空分辨率的可压缩湍流结构非接触精细测试技术和低噪声的高速风洞设备技术. 近几年来, 由于在低噪声的超声速、高超声速风洞研究和可压缩流动精细结构测量技术研究方面取得的重要进展及其在可压缩湍流机理研究方面的应用, 超声速流动转捩与湍流的机理研究取得了较大的进展. 本文介绍了最近几年高速流动非接触精细测试技术, 尤其是基于纳米粒子的平面激光散射技术(nano-tracer planar laser scattering, NPLS)、背景导向纹影技术(background oriented schlieren, BOS) 和超声速流场的粒子图像测速技术(particle image velocimetry, PIV)的研究进展和发展前景, 以及基于这些技术, 在可压缩湍流机理实验研究方面的进展和发展前景, 其中包括在超声速混合层转捩、超声速绕流与尾流结构、超声速边界层转捩、激波边界层干扰等典型流场的机理研究方面, 以及气动光学机理研究方面的研究进展. 最后, 展望了目前湍流机理实验研究对湍流工程模型研究的可能贡献.      

可压缩燃烧反应转捩混合层直接数值模拟

针对三维时间发展可压缩氢/氧非预混燃烧反应平面自由剪切混合层,采用5阶迎风/6阶对称紧致混合差分格式以及3阶显式Runge-Kutta时间推进方法,直接数值模拟了伴随燃烧产物生成和反应能量释放, 流动受扰动激发失稳并转捩的演化过程. 在转捩初期, 获得了${\it\Lambda}$涡、马蹄涡等典型的大尺度拟序结构,观察到了流动失稳后发生双马蹄涡三维对并的现象, 大尺度结构呈较好的对称性.在流动演化后期, 大尺度结构逐次破碎形成小尺度结构, 混合层进入转捩末期,呈明显的不对称性.      

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

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

可压缩气固混合层中离散相与连续相的相互作用研究

尽管已有许多文献采用数值模拟方法研究两相流问题,但主要是集中不可压流动方面.本文采用Eul-er-Lagrange颗粒-轨道双向耦合模型对时间模式下含有固粒的二维可压缩混合层流动进行了研究.气相流场采用非定常全Navier-Stokes方程描述,并应用具有空间三阶精度的WNND(Weighted Non-Oscillatory, Contai-ning No Free Parameters and Dissipative)格式进行数值高散.固相方程采用二阶单边三点差分离散.在考虑流场对固粒作用的同时,也计及颗粒对流场的反作用.主要研究混合层大尺度涡对颗粒扩散特性的影响及颗粒对流场结构的影响问题.在对流马赫数为0.5时,研究不同Stokes数颗粒在连续流场中的扩散特性,而在对流马赫数为0.8时研究了不同Stokes数颗粒对流场小激波结构的影响.     

可压缩自由剪切流混合转捩大涡模拟

针对湍流气动光学效应与冲压发动机气体混合机理问题,开展了可压缩混合层流动空间模式大涡模拟和时间模式直接数值模拟研究.通过对流场(包含亚/亚混合、超/亚混合两种情况)失稳、转捩直至完全湍流的空间发展过程的研究表明,对流Mach数0.4状态下流场失稳以二维最不稳定扰动为主;非线性发展中,基频涡对并及展向涡撕裂主控流动转捩,流场发生混合转捩;转捩后脉动流场基本达到局部各向同性,此时,湍流Mach数低于0.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.     

PIV measurements of turbulent flow in planar mixing layer

A turbulent mixing layer consists of two different flow types, i.e. shear layer (shear-flow turbulence) and free stream regions (nearly homogeneous turbulence). The inherent non-uniform seeding tracer distributions observed around the interfaces between the shear layer and two free stream regions usually lead to a difficulty in particle image velocimetry (PIV) measurements. A parametric study on the application of PIV to the measurement of velocity field in a planar mixing layer is made by means of six factors, including interrogation window size, aspect ratio of interrogation window, interrogation window offset, threshold of data validation, sharpening spatial filters (Prewitt and Sobel masks), and smoothing spatial filter (median mask). The objective of this study is to obtain accurate turbulent measurements in both mean and fluctuating velocities using PIV under an appropriate parametric setting. The optimal levels, which are trade-off in between the accuracy and fine spatial resolution of velocity field measurements, are determined with the aid of the Taguchi method. It is shown that the PIV measurements made with this optimal set of parameters are in good agreement with the measurements made by a two-component hot-wire anemometer. Case independency of the proposed optimal set of parameters on the flow condition of the mixing layer is validated through the applications to two additional tests under the different experimental conditions in changing solely either velocity ratio of high-speed to low-speed free stream velocities or Reynolds number.     

Flow property and self-similarity in steady hydraulic jumps

The flow structure in a steady hydraulic jump in both the non-aerated and aerated regions was measured using the image-based particle image velocimetry and bubble image velocimetry techniques, respectively. Three highly aerated steady jumps with Froude numbers varying from 4.51 to 5.35 were tested, and a weak jump with a Froude number of 2.43 was generated for comparison. Mean velocities and turbulence statistics were obtained by ensemble averaging the repeated velocity measurements. Based on the mean velocities, the flow structure in the steady jumps was classified into four regions to distinguish their distinct flow behaviors; they are the potential core region, the boundary layer region, the mixing layer region, and the recirculation region. The flow structure in the weak jump features only three regions without the recirculation region. In addition, spatial variations of mean velocities, turbulence intensity, and Reynolds stresses were also presented. It was observed that the maximum horizontal bubble velocity and maximum horizontal water velocity occur at the same location in the overlapping regions of potential core and mixing layer. The ratio between the maximum horizontal bubble velocity and maximum horizontal water velocity is between 0.6 and 0.8, depending on the Froude number. Examining the mean horizontal bubble velocities in the mixing layer, a similarity profile was revealed with representative mixing layer thickness as the characteristic length scale and the difference between the maximum positive and maximum negative velocities as the characteristic velocity scale. It was also found that the mean horizontal water velocities in the near-wall region are self-similar and behave like a wall jet. Further analyzing autocorrelation functions and energy spectra of the water and bubble velocity fluctuations found that the energy spectra in the water region follow the ?5/3 slope, whereas the spectra in the bubble region follow a ?2/5 slope. In addition, the integral length scale of bubbles is one order of magnitude shorter than that of water.     

DEM simulation of particle mixing in a sheared granular flow

Mixing behaviors of particles are simulated in a sheared granular flow using differently colored but otherwise identical glass spheres, with five different bottom wall velocities. By DEM simulation, the solid fractions, velocities, velocity fluctuations and granular temperatures are measured.The mixing layer thicknesses are compared with the calculations from a simple diffusion equation using the data of apparent self-diffusion coefficients obtained from the current simulation measurements. The calculations and simulation results showed good agreements, demonstrating that the mixing process of granular materials occurred through the diffusion mechanism.     

DEM simulation of particle mixing in a sheared granular flow

Mixing behaviors of particles are simulated in a sheared granular flow using differently colored but otherwise identical glass spheres, with five different bottom wall velocities. By DEM simulation, the solid fractions, velocities, velocity fluctuations and granular temperatures are measured. The mixing layer thicknesses are compared with the calculations from a simple diffusion equation using the data of apparent self-diffusion coefficients obtained from the current simulation measurements. The calculations and simulation results showed good agreements, demonstrating that the mixing process of granular materials occurred through the diffusion mechanism.     

Experiments to study interactions between baroclinic lower flows and a stably stratified upper layer

Experiments were conducted on a rotating fluid annulus to study the basic interactions between baroclinic lower flows and a stably stratified upper layer. Sufficiently stable stratification is necessary for steady flows to emerge in the lower layer. Upward fluid motions make the baroclinic flows permeate into the upper layer. The stable stratification, however, suppresses upward motions so that zonal fluid velocities decrease with height. In fact, their maximum appears at the top level of the baroclinic lower layer and the sign of the radial temperature gradient changes there; namely, it is warmer on the inner side of the annulus in the upper layer. This temperature profile is reflected in a meridional fluid circulation mixing both layers. In the upper layer of the wave flow, there exists a critical level below and above which the zonal fluid velocities have opposite directions for the wave to have a phase shift of half a wavelength in appearance. The experimental results correspond to real atmospheric phenomena.     

多孔介质壁面剪切湍流速度时空关联的研究

作为一个基础统计量,时空关联函数在湍流问题的研究中有着广泛的应用,是研究湍流噪声、湍流中物质扩散和大涡模拟亚格子模型等问题的重要参考.本文通过建立三维多孔结构壁面剪切湍流模型,采用含Darcy-Brinkman-Forchheimer作用力项的格子Boltzmann方程对无穷大多孔介质平行板之间壁湍流进行了数值模拟,进而研究其速度脉动时空关联函数的统计特性.一方面,根据计算得到的流场数据,对比分析了常规槽道湍流与多孔介质壁面槽道湍流的时间关联函数.另一方面,计算并讨论了不同孔隙率和渗透率的多孔介质壁面对速度脉动时空关联性的影响.通过研究表明:多孔结构壁面剪切湍流的时空关联函数等值线与椭圆理论相符;在研究参数范围内,多孔介质壁面的速度时空关联系数随着孔隙率增大而增大,随着渗透率增大而减小.同时发现在槽道壁面的近壁区、过渡区、对数律区和中心区等不同位置处,速度时空关联呈现较大差异性:越远离壁面位置(对数律区和中心区),其时空关联函数所呈现的关联等值线椭圆越细长,高值相关等值线越集中.多孔介质主要改变速度时空关联椭圆图像的椭圆率,说明多孔介质壁面主要影响湍流横扫速度.     

An LES Turbulent Inflow Generator using A Recycling and Rescaling Method

The present paper describes a recycling and rescaling method for generating turbulent inflow conditions for Large Eddy Simulation. The method is first validated by simulating a turbulent boundary layer and a turbulent mixing layer. It is demonstrated that, with input specification of mean velocities and turbulence rms levels (normal stresses) only, it can produce realistic and self-consistent turbulence structures. Comparison of shear stress and integral length scale indicates the success of the method in generating turbulent 1-point and 2-point correlations not specified in the input data. With the turbulent inlet conditions generated by this method, the growth rate of the turbulent boundary/mixing layer is properly predicted. Furthermore, the method can be used for the more complex inlet boundary flow types commonly found in industrial applications, which is demonstrated by generating non-equilibrium turbulent inflow and spanwise inhomogeneous inflow. As a final illustration of the benefits brought by this approach, a droplet-laden mixing layer is simulated. The dispersion of droplets in the near-field immediately downstream of the splitter plate trailing edge where the turbulent mixing layer begins is accurately reproduced due to the realistic turbulent structures captured by the recycling/rescaling method.     

稳定分层流动中湍流特性的研究综述

本文总结了近60 年来分层流动中湍流特性研究的成果. 主要从两个方面进行了综述:(1) 分层流动中湍流场的演变和混合. 在这方面主要分析稳定分层对湍流混合和湍流结构的影响, 以及混合层内湍流结构的特性和混合层的演化规律. (2) 分层流动中湍流的扩散和输运. 动量和标量的逆梯度输运特性是分层湍流研究的一个重要方向;分析分层对湍流扩散的影响. 并指出了一些值得今后进一步研究的方向.     

混合层中的流向涡与Rayleigh离心不稳定

用有限差分方法求解三维Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程和连续性方程，对时间发展的平面混合层中流向涡的产生原因进行了分析．将Ｒａｙｌｅｉｇｈ的轴对称无粘离心不稳定理论推广应用于分析混合层的二维基本流，并导出一无量纲量Ｒａｙ＝－（ｒ／νθ）νθ／ｒ，其中νθ为一流体质点相对于平均速度的速度，ｒ为该质点迹线的曲率半径．当Ｒａｙ＞１．０时就会发生离心不稳定．采用这一判别式后发现混合层中展向涡的周围，尤其是在辫带区，的确存在离心不稳定区域，而过去的实验结果也表明三维不稳定产生于展向涡之间的辫带区．因此有理由认为，除非雷诺数特别小，离心不稳定是流向涡产生和发展的主要原因