基于PIV技术分析颗粒在湍流边界层中的行为

采用粒子图像测速技术(particle image velocimetry,PIV)在平板湍流边界层内开展实验研究,对比颗粒相及单相液体的平均速度剖面、湍流强度、雷诺应力等湍流统计量,分析颗粒在湍流边界层中的行为.利用空间多尺度局部平均涡量的概念提取壁湍流发卡涡展向涡头(顺向涡)并统计其数量规律,得到不同法向位置处顺向涡周围流向脉动速度及流线的空间拓扑结构,比较分析顺向涡发展程度及周围的湍流相干结构.结果发现:与清水工况相比,颗粒相湍流边界层的缓冲层变薄、对数律区下移,湍流强度得到增强,雷诺应力在对数律区有所增大;颗粒的流向脉动速度在展向涡周围的分布与清水工况不同,颗粒能够被流体展向涡周围的猝发过程有效传递;颗粒相的顺向涡涡核较大,且随着法向位置的升高逐渐发展完整,涡和条带在流向上拉伸得更长;同时发现在两种工况下,顺向涡的左下方始终存在一个逆向涡,颗粒相逆向涡的形成弱于单相流体;两种工况下的顺向涡数量均随着法向位置的升高而减少,最后逐渐趋于稳定.      

基于PIV技术分析颗粒在湍流边界层中的行为

采用粒子图像测速技术(particle image velocimetry, PIV)在平板湍流边界层内开展实验研究,对比颗粒相及单相液体的平均速度剖面、湍流强度、雷诺应力等湍流统计量,分析颗粒在湍流边界层中的行为.利用空间多尺度局部平均涡量的概念提取壁湍流发卡涡展向涡头(顺向涡)并统计其数量规律,得到不同法向位置处顺向涡周围流向脉动速度及流线的空间拓扑结构,比较分析顺向涡发展程度及周围的湍流相干结构.结果发现:与清水工况相比,颗粒相湍流边界层的缓冲层变薄、对数律区下移,湍流强度得到增强,雷诺应力在对数律区有所增大;颗粒的流向脉动速度在展向涡周围的分布与清水工况不同,颗粒能够被流体展向涡周围的猝发过程有效传递;颗粒相的顺向涡涡核较大,且随着法向位置的升高逐渐发展完整,涡和条带在流向上拉伸得更长;同时发现在两种工况下,顺向涡的左下方始终存在一个逆向涡,颗粒相逆向涡的形成弱于单相流体;两种工况下的顺向涡数量均随着法向位置的升高而减少,最后逐渐趋于稳定.     

壁面加热湍流相干结构的子波分析实验研究

通过在平板湍流边界层沿流向固壁表面平行放置若干条通电加热的金属细丝,在平板表面形成沿展向周期性分布的温度场,利用该温度场引起的空气热对流,在湍流边界层近壁区域产生一组沿湍流边界层展向周期分布的大尺度流向涡结构,改变了平板湍流边界层中不同尺度结构及其能量分布。采用对壁湍流多尺度结构的子波分析表明,在湍流边界层近壁区域产生规则的流向涡结构将壁湍流各种尺度湍涡结构不规则的脉动有序地组织起来,抑制了壁湍流各种尺度湍涡结构脉动,特别抑制了能量最大尺度湍涡结构的脉动,减小由于湍流脉动引起的在湍流边界层法向和展向的动量和能量损耗,从而减小了湍流的阻力。     

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

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

EXPERIMENTAL STUDY ON QUANTITATIVE MEASUREMENT OF THREE-DIMENSIONAL STRUCTURE OF HAIRPIN VORTEX BY STEREO-PIV 1)

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

湍流边界层拟序结构的实验研究

20世纪60年代后, 先后从流动显示发现了快慢斑、猝发、上升流、下扫流和多种涡结构等湍流边界层的拟序结构. 它们对湍流边界层的摩阻、传热传质和湍动能的产生等特性有重要影响. 涡结构是上述拟序结构的核心, 它影响其它拟序结构的发展和演变. 发卡涡通常被认为是基本涡结构. 发卡涡等涡结构的再生, 是湍流边界层拟序结构能够自持续的必要的因素.壁面低速流上升产生猝发, 是湍流边界层湍能的主要来源; 条件采样是测量猝发频率和其它拟序结构出现频率的重要手段. 流动显示对湍流边界层拟序结构作了大量定性观察, 有许多减阻和增加传热率等应用性研究在此基础上发展起来. 80年代后, 出现了测量湍流边界层的瞬时流速矢量场的多热线法和PIV技术, 三维PIV技术可望将来为湍流边界层的实验研究带来重大进展. 本文评述了流动显示法、多热线法和PIV技术的优点和不足之处, 以及它们在对湍流边界层拟序结构的研究中的贡献.      

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

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

逆向涡对超疏水壁面减阻影响的TRPIV实验研究

超疏水壁面由于具有减阻和自清洁功能而成为国内外减阻和海洋防污等研究领域的热点之一,而20世纪湍流中相干结构的发现为湍流的控制指出新的方向,尤其近壁区涡结构对摩擦阻力贡献很大.利用高时间分辨率粒子图像测速技术,研究了超疏水壁面(SH)以及亲水壁面(PH)湍流边界层中正负展向涡的空间分布特征,研究逆向涡对超疏水壁面近壁区流动结构的影响和超疏水壁面的减阻机理.首先利用空间多尺度局部平均涡量的概念提取壁湍流发卡涡展向涡头(顺向涡)和逆向涡,实现了准确识别涡心并排除小尺度涡的干扰;然后根据检测到的顺向涡和逆向涡流线分布图,发现逆向涡始终处于正向涡的上游和下方,并且对正向涡的进一步发展起抑制作用;最后对两种壁面边界层中逆向涡数量以及出现概率进行对比,发现具有减阻效果的超疏水壁面边界层中出现更多逆向涡.说明逆向涡可抑制上方顺向涡与壁面的强烈剪切,并使靠近壁面的流体加速,从而产生减阻效果;超疏水壁面中涡结构具有更大的β角,使其更好地阻碍了发卡涡头附近强烈的喷射和扫略;超疏水壁面逆向涡出现概率明显大于亲水壁面.这些结果表明:超疏水壁面表现出的减阻特性(Reδ≈13 500,减阻5.8%)与两板产生逆向涡的差异有关.     

金属丝壁面加热湍流标度律的实验研究

通过在固壁表面的平板湍流边界层沿流向平行放置若干通电加热的金属细丝,在平板表面形成沿展向周期性分布的温度场,利用该温度场引起的空气热对流,在湍流边界层近壁区域产生一组沿湍流边界层展向周期分布的流向涡结构。对壁湍流小尺度结构标度律统计特性的研究表明,金属丝加热后形成的规则流向涡结构将壁湍流各种尺度湍涡结构不规则的脉动有序地组织起来,增强了湍流小尺度结构的层次结构相似性,减小了壁湍流中小尺度结构的间歇性和奇异性,抑制了壁湍流中奇异的湍涡结构。     

圆柱尾迹影响旁路转捩末期发卡涡涡包的研究

基于TR-PIV技术, 通过侧视和俯视两种情况对圆柱尾迹影响下旁路转 捩末期发卡涡涡包的结构及特征尺寸进行了实验研究. 结合二维空间子波变换和\lambda _{ci}准则, 运用线性随机估计方法对速度信号进行条件平均. 在侧视情况下, 条件平均结 果显示, 在边界层中一系列发卡涡涡头与壁面构成17^{\circ}的倾角, 并且被尾迹涡所占据的低 速区域出现在涡包上方的主流区中. 在俯视的结果中, 沿流向方向拉伸(流向尺度 3\delta, 展向尺度0.55δ)的低速条带结构出现在法向高 度为y/δ =0.2的流向-展向平面中, 并且在该低速条带的两侧对称地出 现了沿流向分布的反向旋转的涡结构. 可以得出: 在圆柱尾迹影响下旁路转捩的末期, 由于 尾迹涡诱导作用的影响, 发卡涡涡包在形态上显示出了更大尺度的特征.     

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     

Observations of meandering superstructures in the roughness sublayer of a turbulent boundary layer

The existence of meandering superstructures in a turbulent boundary layer overlying irregular roughness is explored with high-frame-rate particle-image velocimetry measurements within the roughness sublayer at moderate Reynolds number. Elongated streamwise fields of view reconstructed from this data using a Taylor’s hypothesis approach revealed spanwise-meandering regions of low-streamwise-momentum fluid that extended several boundary-layer thicknesses in the streamwise direction. Such signatures are consistent with previous observations of superstructures in smooth-wall turbulence. Counter-rotating wall-normal vortex cores were found to reside along the spanwise boundaries of these superstructure patterns, consistent with the spatial characteristics of hairpin vortex packets.     

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.     

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.     

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.     

Vortex ring/viscous wall layer interaction model of the turbulence production process near walls

An experimental simulation of the interaction of vortex ring-like eddies with the sublayer of a turbulent boundary layer is investigated. An artificially generated vortex ring interacting with a Stokes' layer enables investigation of the interaction with reproducible initial conditions and in the absence of background turbulence. All of the observed features in the turbulent boundary layer production process such as the streaky structure, the pockets, the hairpin vortices, streak lift-up, oscillation, and breakup, have been observed to form. The model shows us that hairpin vortices can pinchoff and reconnect forming new vortex ring-like eddies. Interestingly, the model includes interactions that occur with low probability in the turbulent boundary layer, but which contribute significantly to transport, and may be the events most readily controllable.     

DPIV, HPTV and visualization study of a vortex ring–moving wall interaction

This paper presents an experimental study of the interaction between a vortex ring and a moving wall. This type of flow can be considered as modeling, in a simplified way, the interaction between a "typical eddy" and the viscous sublayer of a turbulent boundary layer. In the present study, the vortex ring is considered as a three-dimensional (3D) perturbation of a viscous Stokes layer. The interaction was first characterized by visualization. To obtain quantitative information, digital particle image velocimetry (DPIV) and holographic particle-tracking velocimetry (HPTV) were used. These different techniques led to a precise and detailed characterization of the vortex ring alone and of an interaction in which a hairpin vortex is generated in the Stokes layer. The results obtained show a good similarity between the observed vortex ring and the Oseen model. They also validate the Stokes layer model and show that in the present conditions, the hairpin vorticity is comparable to that of the Stokes layer. The holographic study, which was undertaken to obtain full 3D three-component (3D3C) velocity maps, showed the present limitations of HPTV.     

Time-resolved PIV measurements of the interaction of polystyrene beads with near-wall-coherent structures in a turbulent channel flow

Time-resolved PIV measurements were performed in a dilute particle-laden flow tracking near-neutrally buoyant polystyrene beads and the velocity field of a near wall turbulent boundary layer. Data were taken in a vertical light sheet aligned in the streamwise direction at the center of a horizontal, closed loop, transparent square water channel facility. In addition, low speed measurements were performed characterizing the effects of the dispersed phase on mean and turbulence flow quantities. Reynolds shear stress slightly differed from clear water conditions whereas fluid mean and rms values were not affected. A case study for several beads revealed a clear relation between their movement and near-wall coherent structures. Several structures having 2D vorticity signatures of near-wall hairpin vortices and hairpin packets, directly affected bead movement. A statistical analysis showed that the mean streamwise velocity of ascending beads lagged behind the mean fluid velocity and bead rms values were higher than fluid ones. Particle Reynolds numbers based on the magnitude of the instantaneous relative velocity vector peaked near the wall; values not exceeding 100, too low for vortex shedding to occur. Quadrant analysis showed a clear preference for ascending beads to reside in ejections while for descending beads the preference for sweeps was less.