热湍流研究的新十年:从聚焦传统到延伸拓展

热湍流(浮力驱动湍流)作为一种典型的湍流现象,广泛存在于自然界和工程应用中. Rayleigh-Bénard (RB)湍流是从众多自然现象中抽象出来研究热湍流的经典模型, RB湍流的典型特征是系统中存在大尺度环流和羽流等不同尺度的湍流结构,这些结构通过作用于边界层,影响RB湍流的输运效率.因此,明确不同尺度湍流结构的生成、演化和作用机理,对理解RB湍流的输运特性至关重要,也是通过控制湍流结构调控输运效率的科学基础.本文重点从湍流结构的时空演化规律、输运特性、湍流调控和热湍流在其他领域的拓展四个方面评述近十年来RB湍流研究所取得的新进展,并对今后的研究方向做出展望.     

稳定分层湍流中逆梯度输运的关联分析研究

本文采用关联分析方法研究了稳定温度分层湍流中的结构特性、输运特性，以及热量、动量逆梯度输运现象的尺度效应及其参数演化．首先采用大涡模拟方法对稳定分层湍流中的结构特性和输运特性进行了分析，将逆梯度输运发生的时间尺度作为已知条件，结合关联量分析方法在波数空间中的解析解，对逆梯度输运现象的尺度效应进行了分析研究．结果发现，稳定分层强度较大的流动中发生垂向热量及动量逆梯度输运现象，发生的结构尺度与关联分析所发现垂向热量、动量逆梯度输运的波数形成了呼应．随着分层强度增加，热量、动量的输运强度均受抑制，与逆梯度输运关联的流场结构尺度减小，同样的效应也发生在流场结构向下游演化的过程中．     

湍流逆梯度输运的理论与实验研究

湍流中存在的动量、标量的逆梯度输运是当前湍流研究中一个十分重要的方向,它具有深刻的理论意义和重要的实用背景．本文从实验、理论两方面综述了湍流逆梯度输运的研究成果,并指出了一些值得今后进一步研究的方向．     

大涡模拟研究展向旋转槽道湍流流动

本文用大涡模拟预测了以不同转速做展向旋转的槽道湍流流动,统计平均的流向速度型在壁面附近与已有实验数据符合很好,在通道中部的预测差异也能给出合理解释,对比不同转速的计算结果,表明展向旋转通道的湍流应力和壁面摩擦力在压力面附近提高、在吸力面附近降低,这些高阶湍流统计量的变化规律可以结合湍流应力输运方程加以解释,漩涡识别技术显示了近壁条带结构,其形态和猝发率受旋转附加力的影响发生改变,进而影响壁面摩擦速度的数值和分布,进一步考察垂直流动方向的截面内速度分布,发现旋转引起了垂直壁面方向的流动,形成正负相间排列的流向涡对,并随着转速的增加向压力面靠近。      

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

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

可压缩混合层流动近十年研究进展

从稳定性分析、实验研究和数值模拟3个方面回顾了过去近10年间关于可压缩混合层流动的研究.这些工作研究了可压缩性对混合层其它性质的影响,这些性质包括混合层增长率、不稳定模式、流场结构、平均速度和脉动速度、以及标志物的混合等.3种研究方法相互补充, 共同提供了认识可压缩混合层流动机制的大量信息.这些研究是湍流机理研究的重要组成部分,也是混合层流动在应用技术领域成功应用的前提.      

圆柱绕流尾迹中相干结构对湍流特性的影响

实验研究了圆柱尾迹中相干结构对湍流平均量的影响.用一个X热线在距离圆柱体y/d=4测量参考信号,用X热线和冷线相结合的三线探头测量从x/d=10～40的圆柱绕流尾迹中的速度和温度脉动,用条件平均的方法研充圆柱绕流近尾迹中的相干结构,并对相干结构对动量和热量的湍流输运过程的影响进行初步分析.结果表明:相干结构对尾迹中速度的横向脉动影响最大;相干结构对湍流参数的影响随x/d的变化明显.该工作还对相干结构引起的湍动能产生率变成负和逆温度梯度输运现象做了定性的解释.     

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

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

对流边界层湍流特性的数值研究

采用大涡模拟方法研究了存在逆温层的情况下大气对流边界层的湍流特性。实际大气边界层中出现逆温层是较常见的，逆温层会导致大气边界层湍流结构的变化，从而影响大气的湍流扩散和输运特性。本文比较了不同逆温梯度的工况，着重分析了逆温层对边界层中热量逆梯度输运（counter gradient heat transportation，CGHT）的影响。计算结果表明：逆温梯度越大，对流边界层的发展越受到抑制；逆温层高度降低会影响整个对流边界层的温度抬升；逆温梯度越大，垂直速度方差越小；在逆温梯度较大的情况下，其逆梯度输运区域要略微低一些，初步分析认为是由于逆温层对热对流的抑制造成的；对于逆温层高度不同的情况，高度越低的逆温层对逆梯度输运的抑制作用更明显。     

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

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

The Effect of Chemical Equilibrium on the Formation of Stable Stratification

The formation of a stable stratification, rather than a uniformly mixed state, is shown to depend critically on the equation of state (EOS) or mixing rule for the density as a function of the solution constitution. For instance, ideal solutions can form any stable stratification at hydrostatic equilibrium, although molecular diffusion will eventually make the concentration of solute uniform. Several other types of EOS are considered at hydrostatic equilibrium, resulting in two and three layer stratifications. In the presence of convective mixing, however, the EOS serves only as a guide, since the dynamics of, for example, turbulent dispersion, significantly alter the formation of stratifications. The case of recirculating flows, however, can be simply treated and results in two layer stratifications for an EOS which is superjacent to the ideal solution. The example of stratification formation in a Burt separator is given.     

Experimental study of stratified jet by simultaneous measurements of velocity and density fields

Stratified flows with small density difference commonly exist in geophysical and engineering applications, which often involve interaction of turbulence and buoyancy effect. A combined particle image velocimetry (PIV) and planar laser-induced fluorescence (PLIF) system is developed to measure the velocity and density fields in a dense jet discharged horizontally into a tank filled with light fluid. The illumination of PIV particles and excitation of PLIF dye are achieved by a dual-head pulsed Nd:YAG laser and two CCD cameras with a set of optical filters. The procedure for matching refractive indexes of two fluids and calibration of the combined system are presented, as well as a quantitative analysis of the measurement uncertainties. The flow structures and mixing dynamics within the central vertical plane are studied by examining the averaged parameters, turbulent kinetic energy budget, and modeling of momentum flux and buoyancy flux. At downstream, profiles of velocity and density display strong asymmetry with respect to its center. This is attributed to the fact that stable stratification reduces mixing and unstable stratification enhances mixing. In stable stratification region, most of turbulence production is consumed by mean-flow convection, whereas in unstable stratification region, turbulence production is nearly balanced by viscous dissipation. Experimental data also indicate that at downstream locations, mixing length model performs better in mixing zone of stable stratification regions, whereas in other regions, eddy viscosity/diffusivity models with static model coefficients represent effectively momentum and buoyancy flux terms. The measured turbulent Prandtl number displays strong spatial variation in the stratified jet.     

A Direct Numerical Study of Transport and Anisotropy in a Stably Stratified Turbulent Flow with Uniform Horizontal Shear

Direct numerical simulations of homogeneous turbulence in stably stratified shear flow have been performed to aid the understanding of turbulence and turbulent mixing in geophysical flow. Two cases are compared. In the first case, which has been studied in the past, the mean velocity has vertical shear and the mean density is vertically stably stratified. In the second case, which has not been studied systematically before, the mean velocity has horizontal shear and the mean density is again vertically stably stratified. The critical value of the gradient Richardson number, for which a constant turbulence level is obtained, is found to be an order of magnitude larger in the horizontal shear case. The turbulent transport coefficients of momentum and vertical mass transfer are also an order of magnitude larger in the horizontal shear case. The anisotropy of the turbulence intensities are found to be in the range expected of flows with mean shear with no major qualitative change in the range of Richardson numbers studied here. However, the anisotropy of the turbulent dissipation rate is strongly affected by stratification with the vertical component dominating the others. This revised version was published online in July 2006 with corrections to the Cover Date.     

Large eddy simulation of turbulent statistical and transport properties in stably stratified flows

Three dimensional large eddy simulation （LES） is performed in the investigation of stably stratified turbulence with a sharp thermal interface. Main results are focused on the turbulent characteristic scale, statistical properties, transport properties, and temporal and spatial evolution of the scalar field. Results show that the buoyancy scale increases first, and then goes to a certain constant value. The stronger the mean shear, the larger the buoyancy scale. The overturning scale increases with the flow, and the mean shear improves the overturning scale. The flatness factor of temperature departs from the Gaussian distribution in a fairly large region, and its statistical properties are clearly different from those of the velocity fluctuations in strong stratified cases. Turbulent mixing starts from small scale motions, and then extends to large scale motions.     

分层流体中栅格湍流的特性

本文考察了分层流体中栅格湍流衰减和演化过程的细节,通过对空间场的信息二维图象处理,获得了湍流动能、耗散率、功率谱及多种湍流尺度等主要湍流特征量。结果显示层结加速了湍流垂直动能的衰减,而水平动能的衰减几乎很少受影响。流场各特征量与前人结果也符合得较好。     

Analysis and Modeling of the Turbulent Diffusion of Turbulent Kinetic Energy in Natural Convection

Buoyant flows often contain regions with unstable and stable thermal stratification from which counter gradient turbulent fluxes are resulting, e.g. fluxes of heat or of any turbulence quantity. Basing on investigations in meteorology an improvement in the standard gradient-diffusion model for turbulent diffusion of turbulent kinetic energy is discussed. The two closure terms of the turbulent diffusion, the velocity-fluctuation triple correlation and the velocity-pressure fluctuation correlation, are investigated based on Direct Numerical Simulation (DNS) data for an internally heated fluid layer and for Rayleigh–Bénard convection. As a result it is decided to extend the standard gradient-diffusion model for the turbulent energy diffusion by modeling its closure terms separately. Coupling of two models leads to an extended RANS model for the turbulent energy diffusion. The involved closure term, the turbulent diffusion of heat flux, is studied based on its transport equation. This results in a buoyancy-extended version of the Daly and Harlow model. The models for all closure terms and for the turbulent energy diffusion are validated with the help of DNS data for internally heated fluid layers with Prandtl number Pr = 7 and for Rayleigh–Bénard convection with Pr = 0.71. It is found that the buoyancy-extended diffusion model which involves also a transport equation for the variance of the vertical velocity fluctuation gives improved turbulent energy diffusion data for the combined case with local stable and unstable stratification and that it allows for the required counter gradient energy flux.     

A von Neumann–Smagorinsky turbulent transport model for stratified shear flows

A simple subgrid turbulent diffusion model based on an analogy to the von Neumann–Richtmyer artificial viscosity is explored for use in modelling mixing in turbulent stratified shear flow. The model may be more generally applicable to multicomponent turbulent hydrodynamics and to subgrid turbulent transport of momentum, composition and energy. As in the case of the von Neumann artificial viscosity and many subgrid-scale models for large-eddy simulation, the turbulent diffusivity explicitly depends on the grid size and is not based on a quantitative model of the unresolved turbulence. In order to address the issue that it is often not known a priori when and where a flow will become turbulent, the turbulent diffusivity is set to zero when the flow is expected to be stable on the basis of a Richardson/Rayleigh–Taylor stability criterion, in analogy to setting the von Neumann artificial viscosity to zero in expanding flows. One-dimensional predictions of this model applied to a simple shear flow configuration are compared to those obtained using a K–ε model. The density and velocity profiles predicted by both models are shown to be very similar.     

稳定分层流密度界面处湍流混合与分形结构

采用室内混合箱研究稳定分层流(上层淡水、下层盐水) 无剪切密度界面处的湍流混合与分形结构. 湍流通过浸没在盐水层中的振动格栅产生, 密度界面结构通过在盐水层中添加荧光剂或染料可视化, 共进行12 组实验. 实验观测并记录了:(1) 淡、盐水密度界面距混合箱底部的平均高程(h);(2) 淡、盐水层的密度(ρ₀,ρ), (3) 淡、盐水密度界面. 其中, 淡、盐水密度界面通过照片、录像进行记录. 观测结果用于计算:(1) 盐水层密度;(2) 卷挟速度, (3) 整体理查孙数(Ri_o), (4) 二维、三维密度界面, (5) 二维、三维密度界面的分形维度. 结果分析发现:(1) 湍流卷挟率随Rio 增大而减小, 并且满足Ri_o的-3=2 或-7=4 幂律;表明随着湍流强度的减弱, 混合的速度也越来越缓慢;(2) 二维密度界面分形维度大于1, 三维密度界面分形维度大于2;表明二维、三维密度界面存在分形结构;(3) 分形维度随Ri_o的增大而减小;表明随着湍流强度的减弱, 密度界面也越来越趋于光滑.