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

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

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

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

FLOW STRUCTURE IN THE TURBULENT BOUNDARY LAYER OVER DIRECTIONAL RIBLETS SURFACES

本文采用时间解析的二维粒子图像测速技术,对零压力梯度光滑以及汇聚和发散沟槽表面平板湍流边界层统计特性和流动结构进行了研究.结果表明在垂直于汇聚和发散沟槽表面的对称平面内,相对于光滑壁面,发散沟槽壁面使当地边界层厚度、壁面摩擦阻力、湍流脉动、雷诺应力等明显减小;而汇聚沟槽壁面对湍流边界层特性和流动结构的影响正好相反,汇聚沟槽使壁面流体有远离壁面向上运动的趋势,因而导致边界层厚度增加了约43%;同时,在汇聚沟槽表面情况下流向大尺度相干结构更容易形成,这对减阻是不利的.此外,顺向涡数量在湍流边界层的对数区均存在一个极大值,发散沟槽表面所对应的极大值位置更靠近沟槽壁面,而在汇聚沟槽表面则有远离壁面的趋势,由顺向涡诱导产生的较强的喷射和扫掠运动会在湍流边界层中产生较强的剪切作用,顺向涡数量的减少是发散沟槽壁面当地摩擦阻力降低的主要原因.     

沟槽方向对湍流边界层流动结构影响的实验研究

本文采用时间解析的二维粒子图像测速技术,对零压力梯度光滑以及汇聚和发散沟槽表面平板湍流边界层统计特性和流动结构进行了研究.结果表明在垂直于汇聚和发散沟槽表面的对称平面内,相对于光滑壁面,发散沟槽壁面使当地边界层厚度、壁面摩擦阻力、湍流脉动、雷诺应力等明显减小;而汇聚沟槽壁面对湍流边界层特性和流动结构的影响正好相反,汇聚沟槽使壁面流体有远离壁面向上运动的趋势,因而导致边界层厚度增加了约43%;同时,在汇聚沟槽表面情况下流向大尺度相干结构更容易形成,这对减阻是不利的.此外,顺向涡数量在湍流边界层的对数区均存在一个极大值,发散沟槽表面所对应的极大值位置更靠近沟槽壁面,而在汇聚沟槽表面则有远离壁面的趋势,由顺向涡诱导产生的较强的喷射和扫掠运动会在湍流边界层中产生较强的剪切作用,顺向涡数量的减少是发散沟槽壁面当地摩擦阻力降低的主要原因.     

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

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

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

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

纵向沟槽壁面湍流边界层内类开尔文-亥姆霍兹涡结构的流动显示

纵向沟槽壁面的湍流边界层,当沟槽的脊-脊内尺度无量纲展向间距s+ 在一定范围内,与光滑壁面湍流边界层相比,具有减阻效应;并在s+ 约为17 个黏性长度单位时减阻效果达到最优,之后其减阻趋弱,直至增阻;其原因可能是沟槽壁面湍流边界层由于“开尔文-亥姆霍兹” 不稳定性而产生的一种“类开尔文-亥姆霍兹”展向涡结构. 实验采用烟雾流动显示技术,首次在风洞中显示了这种“类开尔文-亥姆霍兹” 展向涡结构,确认了其存在性,并在概念上简单勾勒了其结构模型.      

水下湍流减阻途径分析

综述了壁湍流边界层分层结构及近壁区湍流猝发过程,对现有典型水下湍流减阻技术进行了较深入地总结.根据各类减阻方法的特点对水下湍流减阻方法进行了分类,并分析了各类减阻方法对壁湍流流场的影响特征,总结了维持或延长层流、干扰湍流"猝发"、壁面隔离、增大湍流阻尼和改变壁面物性等5种水下湍流减阻可行途径,为开展水下湍流减阻研究及发现新的减阻方法提供参考.     

基于单个压电振子的湍流边界层主动控制

利用安装在壁面上的单个压电振子周期振荡,采用开环主动控制方案,实现了对平板湍流边界层相干结构猝发的主动控制和壁湍流减阻.根据不同的输入电压幅值和频率,完成了10种工况的实验.在压电振子下游2mm处,用热线风速仪和迷你热线单丝探针,精细测量湍流边界层不同法向位置瞬时流向速度信号的时间序列,分析了在Re_?=2183压电振子振动对湍流边界层平均速度剖面、减阻率和相干结构猝发过程的影响.实验结果表明,施加控制的工况使平均速度剖面对数律层上移,产生减阻效果;压电振子振幅越大,减阻率越高,减阻效果越明显;通过对施加控制前后流向瞬时速度的多尺度湍涡结构脉动动能的尺度分析,当压电振子振动频率与壁湍流能量最大尺度的猝发频率相近时,减阻率达到最大,为25%,说明控制壁湍流能量最大尺度相干结构的猝发是实现壁湍流减阻的关键;通过对比相干结构猝发的流向速度分量条件相位平均波形,发现施加控制的工况中相干结构猝发流向速度分量的波形幅值明显降低,且流向速度在扫掠后期高速阶段迅速衰减,缩短了高速流体的下扫过程,说明压电振子的振动能抑制相干结构的高速流体下扫过程,减弱高速流体与壁面的强烈剪切过程,并使近壁区域相干结构的振幅显著减弱,迁移速度加快,从而减小壁面摩擦阻力.      

柔性壁面湍流减阻机理与多尺度相干结构控制的实验研究

用热膜测速仪以高于对应最小湍流时间尺度的分辨率,精细测量了水槽中刚性壁面和柔性壁面平板湍流边界层不同法向位置流向速度分量的时间序列信号,利用湍流边界层近壁区域的对数律平均速度剖面与壁面摩擦速度、流体粘性系数、壁面摩擦切应力等内尺度物理量的关系,在准确测量湍流边界层近壁区域对数律平均速度剖面的基础上,通过非线性迭代求解壁面摩擦速度以及湍流边界层壁面摩擦切应力.结果表明柔性壁的湍流边界层速度分布在对数律层有所上移,缓冲层增厚,说明柔性壁面具有一定的减阻作用. 利用湍流多尺度局部平均结构函数的概念和多尺度局部平均结构函数的瞬时强度因子、平坦因子检测多尺度相干结构及其间歇性的方法,提取了湍流边界层多尺度相干结构的条件相位平均波形.对比研究了刚性壁面和柔性壁面平板湍流边界层近壁区域多尺度相干结构的条件相位平均波形及其间歇性的统计特征,分析了柔性壁面具有减阻作用的物理机理.     

The turbulence vorticity as a window to the physics of friction-drag reduction by oscillatory wall motion

DNS data for channel flow, subjected to spanwise (in-plane) wall oscillations at a friction Reynolds number of 1025, are used to examine the turbulence interactions that cause the observed substantial reduction in drag provoked by the forcing. Following a review of pertinent interactions between the forcing-induced unsteady Stokes strain and the Reynolds stresses, identified in previous work by the present authors, attention is focused on the equations governing the components of the enstrophy, with particular emphasis placed on the wall-normal and the spanwise components. The specific objective is to study the mechanisms by which the Stokes strain modifies the enstrophy field, and thus the turbulent stresses. As such, the present analysis sheds fresh light on the drag-reduction processes, illuminating the interactions from a different perspective than that analysed in previous work. The investigation focuses on the periodic rise and fall in the drag and phase-averaged properties during the actuation cycle at sub-optimal actuation conditions, in which case the drag oscillates by around ±2% around the time-averaged 20% drag-reduction margin. The results bring out the important role played by specific strain-related production terms in the enstrophy-component equations, and also identifies vortex tilting/stretching in regions of high skewness as being responsible for the observed strong increase in the spanwise enstrophy components during the drag-reduction phase. Simultaneously, the wall-normal enstrophy component, closely associated with near-wall streak intensity, diminishes, mainly as a result of a reduction in a production term that involves the correlation between wall-normal vorticity fluctuations and the spanwise derivative of wall-normal-velocity fluctuations, which pre-multiplies the streamwise shear strain.     

Heat transfer enhancement in a turbulent natural convection boundary layer along a vertical flat plate

An experimental study on heat transfer enhancement for a turbulent natural convection boundary layer in air along a vertical flat plate has been performed by inserting a long flat plate in the spanwise direction (simple heat transfer promoter) and short flat plates aligned in the spanwise direction (split heat transfer promoter) with clearances into the near-wall region of the boundary layer. For a simple heat transfer promoter, the heat transfer coefficients increase by a peak value of approximately 37% in the downstream region of the promoter compared with those in the usual turbulent natural convection boundary layer. It is found from flow visualization and simultaneous measurements of the flow and thermal fields with hot- and cold-wires that such increase of heat transfer coefficients is mainly caused by the deflection of flows toward the outer region of the boundary layer and the invasion of low-temperature fluids from the outer region to the near-wall region with large-scale vortex motions riding out the promoter. However, heat transfer coefficients for a split heat transfer promoter exhibit an increase in peak value of approximately 60% in the downstream region of the promoter. Flow visualization and PIV measurements show that such remarkable heat transfer enhancement is attributed to longitudinal vortices generated by flows passing through the clearances of the promoter in addition to large-scale vortex motions riding out the promoter. Consequently, it is concluded that heat transfer enhancement of the turbulent natural convection boundary layer can be substantially achieved in a wide area of the turbulent natural convection boundary layer by employing multiple column split heat transfer promoters. It may be expected that the heat transfer enhancement in excess of approximately 40% can be accomplished by inserting such promoters.     

Investigation of turbulence structures in a drag-reduced turbulent channel flow with surfactant additive by stereoscopic particle image velocimetry

In the present study, we employed stereoscopic particle image velocimetry (PIV) to investigate the characteristics of turbulence structures in a drag-reduced turbulent channel flow with addition of surfactant. The tested drag-reducing fluid was a CTAC/NaSal/Water (CTAC: cetyltrimethyl ammonium chloride; NaSal: sodium salicylate) system at 25°C. The weight concentration of CTAC was 30 ppm. Stereoscopic PIV measurement was performed for a water flow (Re=1.1×10⁴) and a CTAC solution flow (Re=1.5×10⁴ with 54% drag reduction) in both the streamwise–spanwise and wall-normal-spanwise planes, respectively. The three-dimensionality of hairpin vortex structures in the near-wall region for wall-bounded turbulent flow was reproduced by conditionally averaging the stereoscopic two-dimensional-three-component velocity fields. A series of wall-normal vortex cores were found to align with the near-wall low-speed streaks with opposite vorticity signals at both sides of the streaks and with the vorticity decreased on average by about one order of magnitude in CTAC solution flow compared with water flow; the spanwise spacing between the near-wall low-speed streaks in the solution flow is increased by about 46%. The streamwise vorticity of the vortex cores appearing in the wall-normal-spanwise plane was also decreased by the use of drag-reducing surfactant additives.     

Tomographic PIV investigation on coherent vortex structures over shark-skin-inspired drag-reducing riblets

Nature has shown us that the microstructure of the skin of fast-swimming sharks in the ocean can reduce the skin friction drag due to the well-known shark-skin effect.In the present study,the effect of shark-skin-inspired riblets on coherent vortex structures in a turbulent boundary layer(TBL) is investigated.This is done by means of tomographic particle image velocimetry(TPIV) measurements in channel fl ws over an acrylic plate of drag-reducing riblets at a friction Reynolds number of 190.The turbulent fl ws over drag-reducing riblets are verifie by a planar time-resolved particle image velocimetry(TRPIV) system initially,and then the TPIV measurements are performed.Two-dimensional(2D) experimental results with a dragreduction rate of around 4.81% are clearly visible over triangle riblets with a peak-to-peak spacing s+of 14,indicating from the drag-reducing performance that the buffer layer within the TBL has thickened;the logarithmic law region has shifted upward and the Reynolds shear stress decreased.A comparison of the spatial topological distributions of the spanwise vorticity of coherent vortex structures extracted at different wall-normal heights through the improved quadrant splitting method shows that riblets weaken the amplitudesof the spanwise vorticity when ejection(Q2) and sweep(Q4) events occur at the near wall,having the greatest effect on Q4 events in particular.The so-called quadrupole statistical model for coherent structures in the whole TBL is verified Meanwhile,their spatial conditional-averaged topological shapes and the spatial scales of quadrupole coherent vortex structures as a whole in the overlying turbulent fl w over riblets are changed,suggesting that the riblets dampen the momentum and energy exchange between the regions of near-wall and outer portion of the TBL by depressing the bursting events(Q2 and Q4),thereby reducing the skin friction drag.     

Skin-friction drag reduction in turbulent channel flow based on streamwise shear control

It is known that stretching and intensification of a hairpin vortex by mean shear play an important role to create a hairpin vortex packet, which generates the large Reynolds shear stress associated with skin-friction drag in wall-bounded turbulent flows. In order to suppress the mean shear at the wall for high efficient drag reduction (DR), in the present study, we explore an active flow control concept using streamwise shear control (SSC) at the wall. The longitudinal control surface is periodically spanwise-arranged with no-control surface while varying the structural spacing, and an amplitude parameter for imposing the strength of the actuating streamwise velocity at the wall is introduced to further enhance the skin-friction DR. Significant DR is observed with an increase in the two parameters with an accompanying reduction of the Reynolds stresses and vorticity fluctuations, although a further increase in the parameters amplifies the turbulence activity in the near-wall region. In order to study the direct relationship between turbulent vortical structures and DR under the SSC, temporal evolution with initial eddies extracted by conditional averages for Reynolds-stress-maximizing Q2 events are examined. It is shown that the generation of new vortices is dramatically inhibited with an increase in the parameters throughout the flow, causing fewer vortices to be generated under the control. However, when the structural spacing is sufficiently large, the generation of new vortex is not suppressed over the no-control surface in the near-wall region, resulting in an increase of the second- and fourth-quadrant Reynolds shear stresses. Although strong actuating velocity intensifies the near-wall turbulence, the increase in the turbulence activity is attributed to the generation of counter-clockwise near-wall vortices by the increased vortex transport.     

Mechanism of drag reduction by spanwise oscillating Lorentz force in turbulent channel flow

Numerical simulations and experimental research are both carried out to investigate the controlled effect of spanwise oscillating Lorentz force on a turbulent channel flow. The variations of the streaks and the skin friction drag are obtained through the PIV system and the drag measurement system, respectively. The flow field in the near-wall region is shown through direct numerical simulations utilizing spectral method. The experimental results are consistent with the numerical simulation results qualitatively, and both the results indicate that the streaks are tilted into the spanwise direction and the drag reduction utilizing spanwise oscillating Lorentz forces can be realized. The numerical simulation results reveal more detail of the drag reduction mechanism which can be explained, since the spanwise vorticity generated from the interaction between the induced Stokes layer and intrinsic turbulent flow in the near-wall region can make the longitudinal vortices tilt and oscillate, and leads to turbulence suppression and drag reduction.     

Vortex reynolds number in turbulent boundary layers

The effects of vortex Reynolds number on the statistics of turbulence in a turbulent boundary layer have been investigated. Vortex Reynolds number is defined as the ratio of circulation around the vortex structure to the fluid viscosity. The vortex structure of the outer region was modeled and a full numerical simulation was then conducted using a high-order spectral method. A unit domain of the outer region of a turbulent boundary layer was assumed to be composed of essentially three elements: a wall, a Blasius mean shear, and an elliptic vortex inclined at 45° to the flow direction. The laminar base-flow Reynolds number is roughly in the same range as that of a turbulent boundary layer based on eddy viscosity, and the vortex-core diameter based on the boundary-layer thickness is nearly the same as the maximum mixing length in a turbulent boundary layer. The computational box size, namely, 500, 150, and 250 wall units in the streamwise, surface-normal, and spanwise directions, respectively, is approximately the same as the measured quasi-periodic spacings of the near-wall turbulence-producing events in a turbulent boundary layer. The effects of vortex Reynolds number and the signs of the circulation on the moments of turbulence were examined. The signs mimic the ejection and sweep types of organized motions of a turbulent boundary layer. A vortex Reynolds number of 200 describes the turbulence moments in the outer layer reasonably well.     

Streamwise Vortices and Velocity Streaks in a Locally Drag-Reduced Turbulent Boundary Layer

This work aims to understand the changes associated with the near-wall streaky structures in a turbulent boundary layer (TBL) where the local skin-friction drag is substantially reduced. The Reynolds number is R e _?? = 1000 based on the momentum thickness or R e _τ = 440 based on the friction velocity of the uncontrolled flow. The TBL is perturbed via a local surface oscillation produced by an array of spanwise-aligned piezo-ceramic (PZT) actuators and measurements are made in two orthogonal planes using particle image velocimetry (PIV). Data analyses are conducted using the vortex detection, streaky structure identification, spatial correlation and proper orthogonal decomposition (POD) techniques. It is found that the streaky structures are greatly modified in the near-wall region. Firstly, the near-wall streamwise vortices are increased in number and swirling strength but decreased in size, and are associated with greatly altered velocity correlations. Secondly, the velocity streaks grow in number and strength but contract in width and spacing, exhibiting a regular spatial arrangement. Other aspects of the streaky structures are also characterized; they include the spanwise gradient of the longitudinal fluctuating velocity and both streamwise and spanwise integral length scales. The POD analysis indicates that the turbulent kinetic energy of the streaky structures is reduced. When possible, our results are compared with those obtained by other control techniques such as a spanwise-wall oscillation, a spanwise oscillatory Lorentz force and a transverse traveling wave.     

壁湍流相干结构和减阻控制机理

剪切湍流中相干结构的发现是上世纪湍流研究的重大进展之一,这些大尺度的相干运动在湍流的动力学过程中起重要作用,也为湍流的控制指出了新的方向.壁湍流高摩擦阻力的产生与近壁区流动结构密切相关,基于近壁区湍流动力学过程的减阻控制方案可以有效降低湍流的摩擦阻力,但是随着雷诺数的升高, 这些控制方案的有效性逐渐降低.近年来研究发现, 在高雷诺数情况下外区存在大尺度的相干运动,这种大尺度运动对近壁区湍流和壁面摩擦阻力的产生有重要影响,为高雷诺数湍流减阻控制策略的设计提出了新的挑战.该文将对壁湍流相干结构的研究历史加以简单的回顾,重点介绍近壁区相干结构及其控制机理、近年来高雷诺数外区大尺度运动的研究进展,在此基础上提出高雷诺数减阻控制研究的关键科学问题.