激波与转捩边界层干扰非定常特性数值分析

激波与边界层干扰的非定常问题是高速飞行器气动设计中基础研究内容之一.以往研究主要针对层流和湍流干扰,在分离激波低频振荡及其内在机理方面存在着上游机制和下游机制两类截然不同的理论解释.分析激波与转捩边界层干扰下非定常运动现象有助于进一步加深理解边界层状态以及分离泡结构对低频振荡特性的影响规律,为揭示其产生机理指出新的方向.采用直接数值模拟方法对来流马赫数2.9,24?压缩拐角内激波与转捩边界层干扰下激波的非定常运动特性进行了数值分析.通过在拐角上游平板特定的流向位置添加吹吸扰动激发流动转捩,使得进入拐角的边界层处于转捩初期阶段.在验证了计算程序可靠性的基础上,详细分析了转捩干扰下激波运动的间歇性和振荡特征,着重研究了分离泡展向三维结构对激波振荡特性的影响规律,最后还初步探索了转捩干扰下激波低频振荡产生的物理机制.研究结果表明:分离激波的非定常运动仍存在强间歇性和低频振荡特征,其时间尺度约为上游无干扰区内脉动信号特征尺度的10倍量级;分离泡展向三维结构不会对分离激波的低频振荡特征产生实质影响.依据瞬态脉动流场的低通滤波结果,转捩干扰下激波低频振荡的诱因来源于拐角干扰区下游,与流场中分离泡的收缩/膨胀运动存在一定的关联.     

管道凸起激波/湍流边界层干扰的数值模拟

采用基于部分平均的湍流方程组,模拟了管道凸起流动中两种不同强度的激波-湍流边界层干扰现象,对流场的时均参数与实验值进行了比较,计算得到的壁面压力分布、摩阻系数分布、边界层厚度和速度型与实验值比较吻合很好.结果表明基于部分平均的湍流方程组既能够准确模拟小分离的激波/湍流边界层干扰流动又能够准确模拟大分离的激波/湍流边界层干扰流动,较好地预测了典型的λ激波结构.     

高超声速激波湍流边界层干扰直接数值模拟研究

高超声速激波与湍流边界层干扰会导致飞行器表面出现局部热流峰值,严重影响飞行器气动性能和飞行安全. 针对高马赫数激波干扰问题,以往数值研究多采用雷诺平均方法,而在直接数值模拟方面的相关工作较为少见. 开展高超声速激波与湍流边界层干扰的直接数值模拟研究,有助于进一步提升对其复杂流动机理认识和理解,同时也将为现有湍流模型和亚格子应力模型的改进提供理论依据. 采用直接数值模拟方法对来流马赫数6.0,34°压缩拐角内激波与湍流边界层的干扰问题进行了研究. 基于雷诺应力各向异性张量,分析了高超声速湍流边界层在压缩拐角内的演化特性. 通过对湍动能输运方程的逐项分析,系统地研究了可压缩效应对湍动能及其输运的影响机制. 采用动态模态分解方法,探讨了干扰流场的非定常运动历程. 研究结果表明,随着湍流边界层往下游发展,近壁湍流的雷诺应力状态由两组元轴对称状态逐渐演化为两组元状态,外层区域则由轴对称膨胀趋近于各向同性. 干扰流场内存在强内在压缩性效应(声效应),其对湍动能输运的影响主要体现在压力--膨胀项,而对膨胀--耗散项影响较小. 高超声速下压缩拐角内的非定常运动仍存在以分离泡膨胀/收缩为特征的低频振荡特性,其物理机制与分离泡剪切层密切相关.      

高超声速激波湍流边界层干扰直接数值模拟研究

高超声速激波与湍流边界层干扰会导致飞行器表面出现局部热流峰值,严重影响飞行器气动性能和飞行安全.针对高马赫数激波干扰问题,以往数值研究多采用雷诺平均方法,而在直接数值模拟方面的相关工作较为少见.开展高超声速激波与湍流边界层干扰的直接数值模拟研究,有助于进一步提升对其复杂流动机理认识和理解,同时也将为现有湍流模型和亚格子应力模型的改进提供理论依据.采用直接数值模拟方法对来流马赫数6.0,34?压缩拐角内激波与湍流边界层的干扰问题进行了研究.基于雷诺应力各向异性张量,分析了高超声速湍流边界层在压缩拐角内的演化特性.通过对湍动能输运方程的逐项分析,系统地研究了可压缩效应对湍动能及其输运的影响机制.采用动态模态分解方法,探讨了干扰流场的非定常运动历程.研究结果表明,随着湍流边界层往下游发展,近壁湍流的雷诺应力状态由两组元轴对称状态逐渐演化为两组元状态,外层区域则由轴对称膨胀趋近于各向同性.干扰流场内存在强内在压缩性效应(声效应),其对湍动能输运的影响主要体现在压力-膨胀项,而对膨胀-耗散项影响较小.高超声速下压缩拐角内的非定常运动仍存在以分离泡膨胀/收缩为特征的低频振荡特性,其物理机制与分离泡剪切层密切相关.     

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

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

INTERFERENCE OF SHOCK AND BOUNDARY LAYER IN SUPERSONIC FLOW

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

高超声速流场激光测速技术研究进展

高超声速气流条件下飞行器内/外部流动中存在强湍流及脉动、边界层转捩、激波-边界层干扰和高温真实气体效应等耦合效应,表征该非定常流动现象对飞行器气动力、气动热以及目标光电特性等产生的影响是高超声速流动研究中的前沿课题.速度作为表征流动过程最重要的参数之一,准确的速度测量对于深入理解上述复杂流动-传输机理以及高超声速飞行器设计具有重要指导意义.文章针对高超声速流场速度测量中几种常用的非接触式激光测试技术进行了综述,主要包括基于空间法的粒子图像测速,基于激光吸收光谱、激光诱导荧光和瑞利散射的多普勒测速,基于飞行时间法的分子标记测速,以及基于流场折射率的聚焦激光差分干涉测速技术.首先简要介绍每种激光测速技术的基本原理,然后进一步介绍该技术在高超声速自由流、层/湍流边界层、激波/边界层干扰、尾流或其他复杂流动区域的速度及其脉动度测量等方面的典型应用,分析各种技术环境适用性及面临的局限性和挑战.最后对基于激光技术的高超声速流场速度测量进行了总结及发展趋势展望.     

TRPIV MEASUREMENT OF DRAG-REDUCTION IN THE TURBULENT BOUNDARY LAYER OVER RIBLETS PLATE

采用高时间分辨率粒子图像测速技术对沟槽壁面平板湍流边界层速度矢量场的时间序列及其统计量进行了实验测量,讨论了在同一来流速度下沟槽壁面对平均速度剖面﹑雷诺切应力及湍流强度的影响. 用流向速度分量的多尺度空间局部平均结构函数辨识壁湍流多尺度相干结构,用条件采样和相位平均技术提取壁湍流多尺度相干结构喷射和扫掠事件的脉动速度、展向涡量的二维空间拓扑形态. 结果表明,与同材料光滑壁面对比,沟槽壁面实现了10.73%的摩阻减小量;沟槽壁面湍流边界层湍流强度及雷诺切应力皆比光滑平板湍流边界层对应统计量小,说明沟槽壁面有效降低了湍流边界层内流体的脉动. 通过比较壁湍流相干结构猝发事件各脉动速度分量与展向涡量的空间分布特征,肯定了沟槽壁面的减阻效果,发现沟槽壁面通过抑制相干结构猝发事件实现减阻.     

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

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

湍流边界层等动量区演化机理的实验研究

等动量区是瞬时流场中流体动量接近的局部区域,其生成和分布与相干结构密切相关. 对等动量区的研究有助于更深入认识湍流边界层相干结构,但目前对其演化过程还缺乏实验支持和机理分析. 设计并使用移动式高时间分辨率粒子图像测速技术(TRPIV)系统对光滑平板湍流边界层进行了跟踪测量,用滤波方式对数据进行降噪,结合对直接数值模拟数据的插值结果,获得脉动速度信号. 使用改进方法去掉非湍流的影响,检测边界层内的等动量区,得到其数量的时间序列,结合流向速度概率密度函数分布的变化,分析得出了等动量区的数量在大的时间尺度下从一个稳态到另一个稳态的阶梯状变化特点. 分解不同尺度的脉动速度,对大尺度和小尺度脉动信号进行条件平均,发现大尺度脉动对等动量区数量变化起主要作用,表现为不同速度流体通过发生不同猝发事件改变流向速度概率密度函数分布. 分析流向大尺度脉动空间分布的变化,发现等动量区内常含有多个大尺度脉动区域,不同区域的扩张、收缩、分裂、合并影响流向速度的集中程度,进而导致等动量区数量的变化.      

TRPIV investigation of space-time correlation in turbulent flows over flat and wavy walls

The present experimental work is devoted to in- vestigate a new space-time correlation model for the turbulent boundary layer over a flat and a wavy walls. A turbulent boundary layer flow at Reo = 2460 is measured by tomographic time-resolved particle image velocimetry （Tomo-TRPIV）. The space-time correlations of instantaneous streamwise fluctuation velocity are calculated at 3 different wall-normal locations in logarithmic layer. It is found that the scales of coherent structure increase with moving far away from the wall. The growth of scales is a manifestation of the growth of prevalent coherent structures in the turbulent boundary layer like hairpin vortex or hairpin packets when they lift up. The resulting contours of the space-time correlation exhibit elliptic-like shapes rather than straight lines. It is suggested that, instead of Taylor hypothesis, the elliptic model of the space-time correlation is valid for the wallbounded turbulent flow over either a flat wall or a wavy wall. The elliptic iso-correlation curves have a uniform preferred orientation whose slope is determined by the convection velocity. The convection velocity derived from the space-time correlation represents the velocity at which the large-scale eddies carry small-scale eddies. The sweep velocity rep- resents the distortions of the small-scale eddies and is intimately associated with the fluctuation velocity in the logarithmic layer of turbulent boundary layers. The nondimensionalized correlation curves confirm that the elliptic model is more proper for approximating the space-time correlation than Taylor hypothesis, because the latter can not embody the small-scale motions which have non-negligible distortions. A second flow over a wavy wall is also recorded using TRPIV. Due to the combined effect of shear layers and the adverse pressure gradient, the space-time correlation does not show an elliptic-like shape at some specific heights over the wavy wall, but in the outer region of the wavy wallbounded flow, the elliptic model remai     

Correlation study in shock wave�Cturbulent boundary layer interaction

Shock wave–turbulent boundary layer interaction is a critical problem in aircraft design. Therefore, a thorough understanding of the processes occurring in such flows is necessary. The most important task is to study the unsteady phenomena, in particular, the low-frequency ones, for this interaction. An experimental study of separated flow has been performed in the zone of interaction of the incident oblique shock wave with a turbulent boundary layer at Mach 2. Two-point correlation data in the separation zone and the upstream flow were obtained and showed that low-frequency oscillations of the reflected shock waves are related to pulsations in the inflow turbulent boundary layer.     

On the use of a two-layer model for large-eddy simulations of supersonic boundary layers with separation

The two-layer modeling approach has become one of the most promising and successful methodology for simulating turbulent boundary layers in the past ten years. In the present study, a mixed wall model for large-eddy simulations (LES) of high-speed flows is proposed which combine two approaches; the thin-Boundary Layer Equations (TBLE) model of Kawai and Larsson (1994) and the analytical wall-layer model of Duprat et al. (2011) for streamwise pressure gradients. The new hybrid model has been efficiently implemented into a three-dimensional compressible LES solver and validated against DNS of a spatially-evolving supersonic boundary layer (BL) under moderate and strong pressure gradients, before being employed for the prediction of nozzle flow separations at different flow conditions, ranging from weakly to highly over-expanded regimes. A good agreement is obtained in terms of mean and fluctuating quantities compared to the DNS results. Particularly, the current wall-modeled LES results are found to perfectly match the DNS data of supersonic BL with/out pressure gradient. It is also shown that the model can account for the effect of the large-scale turbulent motions of the outer layer, indicating a good interaction between the inner and the outer part of the wall layer. In terms of simulations costs and improvements of computing power, the obtained results highlight the capability of the current wall-modeling LES strategy in saving a considerable amount of computational time compared to the wall-resolved LES counterpart, allowing to push further the simulations limits. Furthermore, the application of these computationally low-costly LES simulations to nozzle flow separation allows to clearly identify the origin of the shock unsteadiness, and the existence of broadband and energetically-significant low-frequency oscillations (LFO) in the vicinity of the separation region.     

Numerical investigation of a spatially developing turbulent natural convection boundary layer along a vertical heated plate

Large-eddy simulation (LES) on a spatially developing natural convection boundary layer along a vertical heated plate was conducted. The heat transfer rate, friction velocity, mean velocity and temperature, and second-order turbulent properties both in the wall-normal and the stream-wise direction showed reasonable agreement with the findings of past experiments. The spectrum of velocity and temperature fluctuation showed a -2/3-power decay slope and -2-power decay slope respectively. Quadrant analysis revealed the inclination on Q1 and Q3 in the Reynolds stress and turbulent heat flux, changing their contribution along the distance from the plate surface. Following the convention, we defined the threshold region where the stream-wise mean velocity takes local maximum, the inner layer which is closer to the plate than the threshold region, the outer layer which is farther to the plate than the threshold region. The space correlation of stream-wise velocity tilted the head toward the wall in the propagating direction in the outer layer; on the other hand, the correlated motion had little inclination in the threshold region. The time history of the second invariant of gradient tensor Q revealed that the vortex strength oscillates both in the inner and the outer layers in between the laminar and the transition region. In the turbulent region, the vortex was often dominant in the outer layer. Instantaneous three-dimensional visualization of Q revealed the existence of high-speed fluid parcels associated with arch-shape vortices. These results were considered as an intrinsic structure in the outer layer, which is symmetrical to the structure of canonical smooth/rough wall bounded layer flow in forced convection.     

Reconstruction of velocity fields from wall pressure measurements in a shock wave/turbulent boundary layer interaction

We present here experimental results in a shock wave/turbulent boundary layer interaction at Mach number of 2.3 impinged by an oblique shock wave, with a deflection angle of 9.5°, as installed in the supersonic wind tunnel of the IUSTI laboratory, France. For such a shock intensity, strong unsteadiness are developing inside the separated zone involving very low frequencies associated with reflected shock motions.The present work consists in simultaneous PIV velocity fields and unsteady wall pressure measurements. The wall pressure and PIV measurements were used to characterize the pressure distribution at the wall in an axial direction, and the flow field associated. These results give access for the first time to the spatial-time correlation between wall pressure and velocity in a shock wave turbulent boundary layer interaction and show the feasibility of such coupling techniques in compressible flows. Linear Stochastic Estimation (LSE) coupled with Proper Orthogonal Decomposition (POD) has been applied to these measurements, and first results are presented here, showing the ability of these techniques to reproduce both the unsteady breathing of the recirculating bubble at low frequency and the Kelvin–Helmholtz instabilities developing at moderate frequency.     

Error of large-eddy simulation in the wall pressure fluctuation of a turbulent channel flow

We analyze the error of large-eddy simulation(LES) in wall pressure fluctuation of a turbulent channel flow. To separate different sources of the error, we conduct both direct numerical simulations(DNS)and LES, and apply an explicit filter on DNS data to obtain filtered DNS(FDNS) data. The error of LES is consequently decomposed into two parts: The first part is the error of FDNS with respect to DNS,which quantifies the influence of the filter operation. The second part is the difference between LES and FDNS induced by the error of LES in velocity field. By comparing the root-mean-square value and the wavenumber-frequency spectrum of the wall pressure fluctuation, it is found that the inaccuracy of the velocity fluctuations is the dominant source that induces the error of LES in the wall pressure fluctuation.The present study provides a basis on future LES studies of the wall pressure fluctuation.     

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.     

PIV study on a shock-induced separation in a transonic flow

A transonic interaction between a steady shock wave and a turbulent boundary layer in a Mach 1.4 channel flow is experimentally investigated by means of particle image velocimetry (PIV). In the test section, the lower wall is equipped with a contour profile shaped as a bump allowing flow separation. The transonic interaction, characterized by the existence in the outer flow of a lambda shock pattern, causes the separation of the boundary layer, and a low-speed recirculating bubble is observed downstream of the shock foot. Two-component PIV velocity measurements have been performed using an iterative gradient-based cross-correlation algorithm, providing high-speed and flexible calculations, instead of the classic multi-pass processing with FFT-based cross-correlation. The experiments are performed discussing all the hypotheses linked to the experimental set-up and the technique of investigation such as the two-dimensionality assumption of the flow, the particle response assessment, the seeding system, and the PIV correlation uncertainty. Mean velocity fields are presented for the whole interaction with particular attention for the recirculating bubble downstream of the detachment, especially in the mixing layer zone where the effects of the shear stress are most relevant. Turbulence is discussed in details, the results are compared to previous study, and new results are given for the turbulent production term and the return to isotropy mechanism. Finally, using different camera lens, a zoom in the vicinity of the wall presents mean and turbulent velocity fields for the incoming boundary layer.     

DNS and scaling laws from new symmetry groups of ZPG turbulent boundary layer flow

The Lie group, or symmetry approach, developed by Oberlack (see e.g. Oberlack [26] and references therein) is used to derive new scaling laws for various quantities of a zero pressure gradient turbulent boundary layer flow. The approach unifies and extends the work done by Oberlack for the mean velocity of stationary parallel turbulent shear flows. From the two-point correlation (TPC) equations the knowledge of the symmetries allows us to derive a variety of invariant solutions (scaling laws) for turbulent flows, one of which is the new exponential mean velocity profile that is found in the mid-wake region of flat-plate boundary layers. Further, a third scaling group was found in the TPC equations for the one-dimensional turbulent boundary layer. This is in contrast to the Navier–Stokes and Euler equations, which have one and two scaling groups, respectively. The present focus is on the exponential law in the outer region of turbulent boundary layer corresponding new scaling laws for one- and two-point correlation functions. A direct numerical simulation (DNS) of a flat plate turbulent boundary layer with zero pressure gradient was performed at two different Reynolds numbers Re=750,2240. The Navier–Stokes equations were numerically solved using a spectral method with up to 140 million grid points. The results of the numerical simulations are compared with the new scaling laws. TPC functions are presented. The numerical simulation shows good agreement with the theoretical results, however only for a limited range of applicability. PACS 02.20.-a, 47.11.+j, 47.27.Nz, 47.27.Eq