自由来流扰动引起的旁路转捩研究进展

本文回顾了自由来流扰动引起的边界层旁路转捩研究进展, 所讨论的自由来流扰动包含高湍流度自由来流和尾迹两大类. 首先介绍旁路转捩的定义、途经及其与规则转捩的区别; 继而讨论高湍流度自由来流引起旁路转捩的路径和机理, 涵盖感受性、瞬态增长、湍流斑、条带稳定性等旁路转捩的主要环节. 在此基础上介绍了尾迹引起旁路转捩的研究现状, 包括对流尾迹和大尺度钝体尾涡两种尾迹形式, 分析了转捩路径及机理, 最后给出研究展望.      

基于Stereo-PIV技术的三维发卡涡结构定量测量研究

发卡涡是湍流相干结构研究中最为关注的内容,实现发卡涡三维结构的定量测量并进行流体动力学分析,对深入研究湍流相干结构、实现湍流精准控制等具有重要意义.本研究通过对合成射流装置进行合理控制,使得层流边界层中产生了规则的人造发卡涡结构,进而用体视图像粒子测速仪(Stereo-PIV)锁相实验技术对发卡涡结构所在的三维空间流场进行了定量测量,并得到了一个完整周期内形成的发卡涡三维结构的空间流场. 结果发现,重构所得的三维发卡涡结构质量较高, 实验技术和方案具有可行性.发卡涡结构所在空间流场情况,符合目前人们对于发卡涡、高低速条带、喷射和扫掠事件的常规认识. 此外,对近壁二次流向涡、展向涡量集中区域的展向涡头和强剪切区域、与低速喷射流体相关的汇聚流动和发散流动等有了更细致的认识.同时, 也探讨了"基于二维脉动流场的相关特征去重构发卡涡三维流场"的可行性.为进一步定量探究发卡涡结构的形成演化、不同涡结构的融合及二次诱导等壁湍流相干结构问题提供思路.      

水下发射航行体尾涡不稳定性分析

在水下连续发射过程中前一发航行体尾流会对后一发航行体运动姿态稳定性产生流动干扰现象. 因此, 研究尾流中涡旋结构演变机理对解决多弹体水下连续发射流动干扰难题具有重要的意义. 本文采用改进型分离涡模型与能量方程, VOF多相流模型与重叠网格技术相结合方法, 对航行体水下发射尾流演变过程开展精细化模拟研究, 其中模拟结果和实验吻合度较好, 验证了本文数值方法的有效性. 以航行体尾流区域为重点研究对象, 分析了尾流区瞬态流场分布, 讨论了横流强度和雷诺数对尾涡结构演变以及脉动压力分布特性的影响. 结果表明: 由于尾流区高速流体核心区与低速自由流相互作用导致Kelvin-Helmholtz不稳定现象出现, 可以清晰地发现涡旋结构在剪切力的作用下发生脱落. 在横流条件下, 航行体尾端脱落的涡环与涡腿形成发卡涡, 而多个发卡涡沿轴向间隔排列组成发卡涡包存在于尾流中. 随着横流强度增大, 形成多级发卡涡包结构, 而导致脉动压力二次峰值均出现的主要原因是尾流涡旋流场演变引起的. 随着雷诺数的增大, 尾流中由圆柱形涡和U型涡组成的二次涡结构逐渐明显, 不稳定性加强.      

壁面湍流发卡涡包空间模态的TRPIV实验研究

运用高时间分辨率粒子图像测速(Time-resolved PIV简称TRPIV),测量得到平板湍流边界层流向/法向平面内瞬时速度矢量空间分布的时间序列;采用空间局部平均速度结构函数的概念,识别和提取湍流边界层中大尺度发卡涡包结构的空间特征。发现在湍流边界层中不同法向位置多个正负发卡涡包结构同时交替存在。这些分布在不同法向高度的发卡涡包结构之间通过倾斜的涡量剪切层相联系,构成了湍流边界层中内、外区紧密相连、相互作用的一种稳态的分布方式。     

圆柱尾迹涡的三维演化及结构特征

应用无粘涡丝运动学理论和局部诱导近似(LIA)方法,以Lagrange观点跟踪涡丝在背景流场中运动,用数值方法研究了中等Re数(≈10~3)下圆柱分离尾迹中Kármán涡和涡辫区涡丝的三维演化的机制和动力学过程,及其结构特征。背景流场考虑为尾迹时间平均速度流场和Kármán涡街流场。初始展向小扰动为指数形式和谐波形式。结果指出:Kármán涡和涡辫区中的涡丝具有展向不稳定性,形成流向涡量。在尾迹的初期输运过程中,表现出有序的大、小尺度涡结构。并进一步分析了其产生的机理。     

轴对称钝物体的尾涡结构

本文对于圆球和70°角斜切尾钝头旋成体的尾涡结构作了流态显示和部分热线测量,流态显示结果表明,这类钝物体尾迹中的三维尾涡结构均呈现为一列相互锁定的“涡链”状结构,每一个涡元素自身并不闭合,呈“发卡”形,并与后继的涡相连接,对于70°斜切尾钝头旋成体模型,由于其斜切尾有固定三维尾涡面卷起和脱落方位的作用,故呈现出规则的三维尾涡结构。而圆球模型,由于涡面卷起和涡脱落的方位是随机旋转的,因此“涡链”呈现为被任意扭曲的千姿百态,流态照片各不相同,圆球尾迹的这种既有周期性又有随机旋转的特性,在热线测量的结果中也有所反映。     

圆柱涡激振动的结构-尾流振子耦合模型研究

建立了一个新的结构-尾流振子耦合模型. 流场近尾迹动力学特征被模化为非线性阻尼振子,采用van der Pol方程描述. 以控制体中结构与近尾迹流体间受力互为反作用关系来实现流固耦合. 采用该模型进行了二维结构涡激振动计算,得到了合理的振幅随来流流速的变化规律和共振幅值,并正确地预计了共振振幅值$A_{\max}^\ast$随着质量阻尼参数$\left( {m^\ast + C_A } \right)\zeta$的变化规律,给出了预测$A_{\max }^\ast$值的拟合公式. 采用该模型计算了三维柔性结构在均匀来流和简谐波形来流作用下的VIV响应. 结构在均匀来流作用下振动呈现由驻波向行波的变化过程, 并最后稳定为行波振动形态.在简谐波形来流作用下,结构呈现混合振动形态,幅值随时间呈周期变化.      

基于Navier-Stokes方程残差的隐式大涡模拟有限元模型

对应于湍流的大尺度与小尺度流场信息,本文在有限元的框架下,假设Navier-Stokes方程的解的形函数由大尺度和不可解尺度形函数叠加组成,引入对应的权函数,将Navier-Stokes方程的有限元变分形式分解为大尺度和不可解尺度系统.根据不可解尺度系统,构建基于Navier-Stokes大尺度方程残差的不可解尺度模型,将其代入Navier-Stokes方程的大尺度系统,进而数值求解大尺度系统得到Navier-Stokes方程的大尺度解.该方法无需像传统的大涡模拟方法那样对方程的解进行过滤,通过对形函数进行尺度分解实现解的尺度分解.本文使用该方法的自编程序代码开展了槽道湍流的数值模拟.通过与有限差分大涡模拟、DNS计算结果的比较,发现在使用较少网格情况下该方法预测的平均流向速度在近壁区与DNS数据吻合,在黏性外层略偏高;该方法对雷诺应力预测偏低导致从流向向垂向方向上湍动能输运略偏低.流向速度等值面图显示该方法有效捕捉到了大尺度旋涡结构;同时在近壁区可以观察到明显的低速条带结构.     

倾斜圆柱结构涡激振动研究进展

圆柱结构涡激振动现象在生活中十分常见,如海洋工程中的管道、土木工程中的高耸建筑、桥梁斜拉索,核工程中的热交换器等频繁受到涡激振动影响,诱发结构的疲劳损伤,甚至破坏失效.现阶段,人们对垂直来流作用下圆柱结构涡激振动机理已有较为全面的认识.然而,当圆柱倾斜置于流场中,结构后缘的尾流形态与垂直放置差异显著,结构与流体的耦合作用机理更为复杂.为简化倾斜圆柱涡激振动问题,提出了不相关原则,来流速度被分解为垂直圆柱结构轴向和平行圆柱结构轴向的两个速度分量,仅考虑垂直结构轴向速度分量的影响,忽略平行结构轴向速度分量的影响.近年来,针对倾斜圆柱涡激振动及不相关原则的适用性,出现了大量实验和数值模拟研究成果.为了深化对倾斜圆柱结构涡激振动相关机理的认知,本文全面阐述了倾斜圆柱结构涡激振动响应规律、尾迹流场模式和流体力特性等方面的研究进展,分析了不相关原则的适用范围,探讨了倾斜圆柱结构涡激振动抑制措施,并对今后该领域的研究进行了力所能及的展望.     

流向强迫振荡圆柱绕流的涡脱落模态分析

通过求解采用ALE方法描述的运动坐标系Navier-Stokes方程组,分析均匀来流下雷诺 数为150的静止和流向振荡的圆柱绕流. 主要研究了强迫振荡频率和较大振幅比 (A/D=0.3-1.2)对圆柱升力、阻力变化特性以及涡脱落模态的影响. 研究表 明,流向振荡圆柱绕流存在多种涡脱落模态,如对称S以及反对称A-I, A-III, A-IV等多种形式;比较研究结果,拓展了各模态下对应的锁定区域,并将其分为5个 子区;A-I模态中圆柱受力较以前所知更复杂;通过分析计算结果,发现最大加速度 比Af_{c}^{2}/Df_{s0}^{2}可能是涡脱落模态(尤其是对称S模态)最有效的控制参数.     

Statistical scale of hairpin packets in the later stage of bypass transition induced by cylinder wake

The hairpin packet's structure and its statistical scale in the later stage of bypass transition induced by a cylinder wake are investigated by time-resolved particle image velocimetry from the side and top view, respectively. Linear stochastic estimation is used to achieve the conditionally averaged velocity fields. For the side view case, the conditionally averaged structure consists of a series of swirling motions located along a line inclining at a large angle (18°) from the wall and a low-speed region occupied by the cylinder wake appearing right above them. In the (x, z)-plane at the wall-normal height y/???=?0.2, the dominant structures are shown to be the large-scale regions of low momentum elongated almost over 3?? along the streamwise. The low-speed regions are consistently bordered by small-scale counter-rotating vortice pairs organized in the streamwise with a statistical spanwise width of 0.55??. The results suggest that in the later part of the transitional zone, the upward induction of the cylinder wake enhances both the wall-normal and spanwise extent of the hairpin packets.     

A bypass wake induced laminar/turbulent transition

The process of laminar to turbulent transition induced by a von Karman vortex street wake, was studied for the case of a flat plate boundary layer. The boundary layer developed under zero pressure gradient conditions. The vortex street was generated by a cylinder positioned in the free stream. An X-type hot-wire probe located in the boundary layer, measured the streamwise and normal to the wall velocity components. The measurements covered two areas; the region of transition onset and development and the region where the wake and the boundary layer merged producing a turbulent flow. The evolution of Reynolds stresses and rms-values of velocity fluctuations along the transition region are presented and discussed. From the profiles of the Reynolds stress and the mean velocity profile, a ‘negative' energy production region along the transition region, was identified. A quadrant splitting analysis was applied to the instantaneous Reynolds stress signals. The contributions of the elementary coherent structures to the total Reynolds stress were evaluated, for several x-positions of the near wall region. Distinct regions in the streamwise and normal to the wall directions were identified during the transition.     

TR PIV INVESTIGATION ON THE LAGRANGIAN HAIRPIN VORTICES INDUCED BY A CUBIC DISTURBANCE

运用时分粒图仪 (time-resolved particle imagevelocimetry)对平板边界层中立方体扰动诱导产生的发卡涡涡结构进行了实验测量, 得到了立方体扰动尾迹区中展向平面内的速度矢量分布.基于拉格朗日坐标系, 运用李亚普诺夫指数对流场中的相干结构进行了辨识, 比较欧拉坐标系下的结果发现拉格朗日体系下李亚普诺夫指数场不仅有效地包含欧拉坐标系下涡强所表征的发卡涡涡头信息外, 并且能够呈现出涡颈在空间中的分布, 能够完全刻画出发卡涡涡结构在展向平面内的信息, 还发现李亚普诺夫指数对流场拟序结构的表征有很好的鲁棒性.另外, 李亚普诺夫指数场的互相关分析表明诱导产生的发卡涡涡结构随着向下游的传播, 其涡头不断抬升, 发卡涡涡结构与壁面的倾角也逐渐变大, 有效地刻画了拟序结构的空间演化过程.     

Wake structure of a finite circular cylinder of small aspect ratio

The wake of a finite circular cylinder of small aspect ratio was studied with a seven-hole probe and thermal anemometry. The cylinder was mounted normal to a ground plane and was partially immersed in a turbulent boundary layer. The time-averaged velocity and streamwise vorticity fields showed the development of the tip vortex structures, the extent of the near-wake recirculation zone, the downwash phenomenon and base vortex structures within the boundary layer. The wake structure and power spectra were similar for cylinder aspect ratios of 5 to 9, but a distinctly different behaviour was observed for an aspect ratio of 3.     

Wake development in staggered short cylinder arrays within a channel

Staggered arrays of short cylinders, known as pin?Cfins, are commonly used as a heat exchange method in many applications such as cooling electronic equipment and cooling the trailing edge of gas turbine airfoils. This study investigates the near wake flow as it develops through arrays of staggered pin fins. The height-to-diameter ratio was unity while the transverse spacing was kept constant at two cylinder diameters. The streamwise spacing was varied between 3.46 and 1.73 cylinder diameters. For each geometric arrangement, experiments were conducted at Reynolds numbers of 3.0e3 and 2.0e4 based on cylinder diameter and velocity through the minimum flow area of the array. Time-resolved flowfield measurements provided insight into the dependence of row position, Reynolds number, and streamwise spacing. Decreasing streamwise spacing resulted in increased Strouhal number as the near wake length scales were confined. In the first row of the bundle, low Reynolds number flows were mainly shear-layer-driven while high Reynolds number flows were dominated by periodic vortex shedding. The level of velocity fluctuations increased for cases having stronger vortex shedding. The effect of streamwise spacing was most apparent in the reduction of velocity fluctuations in the wake when the spacing between rows was reduced from 2.60 diameters to 2.16 diameters.     

DIRECT NUMERICAL SIMULATIONS ON THE TURBULENT FLOW PAST A CONFINED CIRCULAR CYLINDER WITH THE INFLUENCE OF THE STREAMWISE MAGNETIC FIELDS1)

绕流是托卡马克装置中液态包层内常见的流动形态,对流场与热量分布有着重要的影响.本文通过直接数值模拟(DNS),研究了不同磁场强度下$Re=3900$的圆柱绕流,分析了磁场强度对于湍流尾迹的影响.无磁场情况下,直接数值模拟的结果与前人的实验及模拟结果吻合很好.圆柱下游的尾迹中,随着流向距离的增大, 流向速度剖面逐渐从U型进化呈V型, 并慢慢趋于平缓,这表明尾迹中的流动结构受圆柱影响逐渐减小.圆柱后方两侧的剪切层中,由于Kelvin-Helmholtz不稳定性的影响,可以清晰地看到小尺度剪切层涡的脱落.通过对无磁场的计算结果施加流向磁场,本文计算了哈特曼数($Ha$)分别为20, 40和80的工况,以研究磁场效应对于湍流的影响.结果表明磁场较弱时,流动依然呈三维湍流状态.随着磁场增强, 近圆柱尾流区受磁场抑制明显,回流区被拉长,剪切层失稳位置向下游转移.圆柱后方的涡结构由于受到竖直方向洛伦兹力的挤压作用,随着哈特曼数的增加尾迹区域逐渐变窄.相比于无磁场情况的涡结构,由于磁场的耗散作用,相应的涡结构尺度变小.该研究不仅扩展了现有磁场下湍流运动的参数范围,对于液态包层的设计及安全运行同样具有重要的理论指导意义和工程应用价值.     

流向磁场作用下圆柱绕流的直接数值模拟

绕流是托卡马克装置中液态包层内常见的流动形态,对流场与热量分布有着重要的影响.本文通过直接数值模拟(DNS),研究了不同磁场强度下$Re=3900$的圆柱绕流,分析了磁场强度对于湍流尾迹的影响.无磁场情况下,直接数值模拟的结果与前人的实验及模拟结果吻合很好.圆柱下游的尾迹中,随着流向距离的增大, 流向速度剖面逐渐从U型进化呈V型, 并慢慢趋于平缓,这表明尾迹中的流动结构受圆柱影响逐渐减小.圆柱后方两侧的剪切层中,由于Kelvin-Helmholtz不稳定性的影响,可以清晰地看到小尺度剪切层涡的脱落.通过对无磁场的计算结果施加流向磁场,本文计算了哈特曼数($Ha$)分别为20, 40和80的工况,以研究磁场效应对于湍流的影响.结果表明磁场较弱时,流动依然呈三维湍流状态.随着磁场增强, 近圆柱尾流区受磁场抑制明显,回流区被拉长,剪切层失稳位置向下游转移.圆柱后方的涡结构由于受到竖直方向洛伦兹力的挤压作用,随着哈特曼数的增加尾迹区域逐渐变窄.相比于无磁场情况的涡结构,由于磁场的耗散作用,相应的涡结构尺度变小.该研究不仅扩展了现有磁场下湍流运动的参数范围,对于液态包层的设计及安全运行同样具有重要的理论指导意义和工程应用价值.      

Absolute and convective instabilities of two-dimensional bluff body wakes in ground effect

Local and global instabilities are investigated of wakes of general two-dimensional bluff bodies placed near and parallel to a plane boundary or ground. A spatio-temporal linear stability analysis is first applied to a four-parameter family of local wake profiles to investigate the fundamental local stability characteristics of the wake in ground effect. The analysis shows significant dependencies of the stability characteristics of the wake on the distance from the wake centreline to the ground (normalised by the wake width), and also on the velocity ratio of the near- and far-ground sides of the wake. The analysis is then compared with earlier experiments on a circular cylinder to examine, according to the transition scenario of the steep global modes, the streamwise variation of the local stability characteristics of the wake in ground effect. The comparison indicates that the near wake region of the cylinder changes from being absolutely unstable to being convectively unstable when the cylinder comes down into the near-ground range in which the von Kármán-type vortex shedding from the cylinder is suppressed, being qualitatively consistent with the transition scenario for general wake-type flows. A possible explanation is also given for the counter-intuitive relation between the thickness of the boundary layer on the ground and the critical gap distance for the cessation of the von Kármán-type vortex shedding in ground effect.     

Flow above the free end of a surface-mounted finite-height circular cylinder: A review

The wake of a surface-mounted finite-height circular cylinder and the associated vortex patterns are strongly dependent on the cylinder aspect ratio and the thickness of the boundary layer on the ground plane relative to the dimensions of the cylinder. Above a critical aspect ratio, the mean wake is characterized by streamwise tip vortex structures and Kármán vortex shedding from the sides of the cylinder. Below a critical aspect ratio, a unique mean wake structure is observed. Recent experimental studies in the literature that used phase-averaged techniques, as well as recent numerical simulations, have led to an improved physical understanding of the near-wake vortex flow patterns. However, the flow above the free end of the finite circular cylinder, and its relationship to the near wake, has not been systematically studied. The effects of aspect ratio and boundary layer thickness on the free-end flow field are also not completely understood, nor has the influence of Reynolds number on the free-end flow field been fully explored. Common features associated with the free end include separation from the leading edge, a mean recirculation zone containing a prominent cross-stream arch (or mushroom) vortex, and reattachment onto the free-surface. Other flow features that remain to be clarified include a separation bubble near the leading edge, one or two cross-stream vortices within this separation bubble, the origins of the streamwise tip or trailing vortices, and various critical points in the near-surface flow topology. This paper reviews the current understanding of the flow above the free end of a surface-mounted finite-height circular cylinder, with a focus on models of the flow field, surface oil flow visualization studies, pressure and heat flux distributions on the free-end surface, measurements of the local velocity field, and numerical simulations, found in the literature.