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

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

固体颗粒对沟槽湍流边界层影响的实验研究

本文采用粒子图像测速技术(particles image velocimetry, PIV)研究固体颗粒对放置在平板湍流边界层中的平壁和沟槽壁面减阻效果的影响. 实验对清水和加入粒径为155 μm聚苯乙烯颗粒的流法向二维速度场信息进行采集, 对不同工况下的平均速度剖面、雷诺应力和湍流度等统计量进行对比, 分析流体在边界层中的行为. 运用空间局部平均结构函数提取了不同工况湍流边界层喷射?扫掠行为的空间拓扑结构并进行比较. 结果发现, 在不同的壁面条件下, 粒子加入后的对数律区中无量纲速度均略大于清水组, 雷诺切应力有所降低, 湍流度有所减弱. 对于不同流场速度下的沟槽而言, 颗粒的加入均降低了壁面附近的阻力, 而颗粒单独作用于光滑壁面的减阻效果并不明显. 加入粒子后的相干结构数目有所增加, 法向脉动速度下降. 沟槽壁面附近的相干结构数目有所增加, 法向脉动速度在自由来流速度较大时有所上升, 在速度较小时有所下降. 这表明不同减阻状况下的沟槽均能将大涡破碎成更多的涡, 并且粒子的加入强化了这种破碎作用.      

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

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

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

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

FLOW STRUCTURE IN THE TURBULENT BOUNDARY LAYER OVER DIRECTIONAL RIBLETS SURFACES

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

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

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

壁湍流多尺度相干结构复涡黏模型的实验研究

在湍流相干结构动力学方程中,非相干结构成分对相干结构贡献的雷诺应力的模型为涡黏性模型,即涡黏性系数乘以相干结构平均速度变形率的形式.基于非相干结构成分对相干结构贡献的雷诺应力与相干结构速度变形率之间存在相位差的事实,在理论上提出了非相干结构成分对相干结构贡献的雷诺应力复涡黏性模型的假设.应用热线测速技术,在低速风洞中对湍流边界层非相干结构成分对相干结构贡献的雷诺应力与相干结构法向速度变形率之间的相位关系进行了实验测量.通过分析湍流相干结构猝发过程中非相干结构成分对相干结构贡献的雷诺应力与相干结构速度变形率之间的相位关系,研究了相干结构雷诺应力分量与流向速度法向梯度之间的相位差沿湍流边界层法向的变化规律,肯定了湍流相干结构复涡黏性系数模型的合理性.     

壁湍流猝发过程中速度分量的相位差对雷诺应力影响的实验研究

用IFA300恒温热线风速仪和×形二分量热线探针,以采样间隔小于最小湍流时间尺度的分辨率,精细测量了风洞中平板湍流边界层不同法向位置的瞬时流向、展向速度分量的时间序列信号.用子波分析辨识壁湍流相干结构猝发事件的能量最大准则,确定壁湍流相干结构猝发事件的时间尺度;用条件相位平均技术提取了相干结构猝发过程中流向、展向脉动速度分量条件相位平均波形,用互相关方法研究了相干结构猝发过程中流向、展向脉动速度分量条件相位平均波形的相位差关系及其对雷诺应力的影响,发现在缓冲层和对数律区,展向脉动速度与流向脉动速度的条件相位平均波形具有不同的相位;当两者相位基本一致时,雷诺应力达到正的最大值,此时湍流相干结构的产生非常活跃;当两者相位差分别集中在90°和270°附近时,雷诺应力的幅值减小并接近于零,此时湍流相干结构的产生和猝发都得到了抑制.     

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

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

输气管道壁面涂料减阻机理的实验研究

用IFA-300热线风速仪以高于对应最小湍流时间尺度的分辨率精细测量了风洞中不同壁面涂料的管道湍流边界层不同法向位置流向速度分量的时间序列信号,利用湍流边界层近壁区域对数律平均速度剖面与壁面摩擦速度、流体黏性系数等内尺度物理量的关系和壁面摩擦速度与壁面摩擦切应力的关系,在准确测量湍流边界层近壁区域对数律平均速度剖面的基础上,间接测量湍流边界层的壁面摩擦阻力．对不同壁面涂料的壁湍流脉动速度信号用子波分析进行多尺度分解,用子波系数的瞬时强度因子和平坦因子检测管道湍流边界层中的多尺度相干结构,提取不同尺度相干结构的条件相位平均波形,对比研究输气管道壁面涂料的减阻机理．     

The effect of streamwise vortices on the turbulence structure of a separating boundary layer

A high Reynolds number flat plate turbulent boundary layer is investigated in a wind-tunnel experiment. The flow is subjected to an adverse pressure gradient which is strong enough to generate a weak separation bubble. This experimental study attempts to shed some new light on separation control by means of streamwise vortices with emphasize on the change in the boundary layer turbulence structure. In the present case, counter-rotating and initially non-equidistant streamwise vortices become and remain equidistant and confined within the boundary layer, contradictory to the prediction by inviscid theory. The viscous diffusion cause the vortices to grow, the swirling velocity component to decrease and the boundary layer to develop towards a two-dimensional state. At the position of the eliminated separation bubble the following changes in the turbulence structure were observed. The anisotropy state in the near-wall region is unchanged, which indicates that it is determined by the presence of the wall rather than the large scale vortices. However, the turbulence in the outer part of the boundary layer becomes overall more isotropic due to an increased wall-normal mixing and a significantly decreased production of streamwise fluctuations. The turbulent kinetic energy is decreased as a consequence of the latter. Despite the complete change in mean flow, the spatial turbulence structure and the anisotropy state, the process of transfer of turbulent kinetic energy to the spanwise fluctuating component seems to be unchanged. Local regions of anisotropy are strongly connected to maxima in the turbulent production. For example, at spanwise positions in between those of symmetry, the spanwise gradient of the streamwise velocity cause significant production of turbulent fluctuations. Transport of turbulence in the spanwise direction occurs in the same direction as the rotation of the vortices.     

Direct numerical simulation of a particle-laden flow in a flat plate boundary layer

In this paper, a direct numerical simulation of particle-laden flow in a flat plate boundary layer is performed, using the Eulerian–Lagrangian point-particle approach. This is, as far as we know, the first simulation of a particle-laden spatially-developing turbulent boundary layer with two-way coupling. A local minimum of the particle number density is observed in the close vicinity of the wall. The present simulation results indicate that the inertial particles displace the quasi-streamwise vortices towards the wall, which, in turn, enhance the mean streamwise fluid velocity. As a result, the skin-friction coefficient is increased whereas the boundary layer integral thicknesses are reduced. The presence of particles augments the streamwise fluctuating velocity in the near-wall region but attenuates it in the outer layer. Nevertheless, the wall-normal and spanwise velocity fluctuations are significantly damped, and so is the Reynolds stress. In addition, the combined effect of a reduced energy production and an increased viscous dissipation leads to the attenuation of the turbulent kinetic energy.     

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

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

Tomographic PIV investigation of coherent structures in a turbulent boundary layer flow

Tomographic particle image velocimetry was used to quantitatively visualize the three-dimensional coherent structures in the logarithmic region of the turbulent boundary layer in a water tunnel.The Reynolds number based on momentum thickness is Reθ = 2 460.The instantaneous velocity fields give evidence of hairpin vortices aligned in the streamwise direction forming very long zones of low speed fluid,which is flanked on either side by highspeed ones.Statistical support for the existence of hairpins is given by conditional averaged eddy within an increasing spanwise width as the distance from the wall increases,and the main vortex characteristic in different wall-normal regions can be reflected by comparing the proportion of ejection and its contribution to Reynolds stress with that of sweep event.The pre-multiplied power spectra and two-point correlations indicate the presence of large-scale motions in the boundary layer,which are consistent with what have been termed very large scale motions(VLSMs).The three dimen-sional spatial correlations of three components of velocity further indicate that the elongated low-speed and highspeed regions will be accompanied by a counter-rotating roll modes,as the statistical imprint of hairpin packet structures,all of which together make up the characteristic of coherent structures in the logarithmic region of the turbulent boundary layer(TBL).     

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.     

Numerical simulation on evolution of subharmonic low-speed streaks in minimal channel turbulent flow

The evolution of low-speed streaks in the turbulent boundary layer of the minimum channel flow unit at a low Reynolds number is simulated by the direct numer- ical simulation （DNS） based on the standard Fourier-Chebyshev spectral method. The subharmonic sinuous （SS） mode for two spanwise-aligned low-speed streaks is excited by imposing the initial perturbations. The possibilities and the physical realities of the turbulent sustaining in the minimal channel unit are examined. Based on such a flow field environment, the evolution of the low-speed streaks during a cycle of turbulent sus- taining, including lift-up, oscillation, and breakdown, is investigated. The development of streamwise vortices and the dynamics of vortex structures are examined. The results show that the vortices generated from the same streak are staggered along the streamwise direction, while the vortices induced by different streaks tilt toward the normal direction due to the mutual induction effect. It is the spatial variations of the streamwise vortices that cause the lift-up of the streaks. By resolving the transport dynamics of enstrophy, the strength of the vortices is found to continuously grow in the logarithmic layer through the vortex stretching mechanism during the evolution of streaks. The enhancement of the vortices contributes to the spanwise oscillation and the following breakdown of the low-speed streaks.     

雷诺应力各向异性涡黏模型的层析TRPIV测量

利用层析TRPIV测量水洞中平板湍流边界层3D-3C速度场的高分辨率时间序列数据库. 提出了空间局部平均多尺度速度结构函数的新概念, 描述湍流多尺度涡结构的空间拉伸、压缩、剪切变形和旋转. 用空间局部平均多尺度速度结构函数对湍流脉动速度进行了空间多尺度分解. 用空间流向局部平均多尺度速度结构函数, 根据湍流多尺度涡结构在流向的拉伸和压缩物理特征, 提出了新的湍流相干结构条件采样方法, 检测并提取了层析TRPIV数据中相干结构“喷射”和“扫掠”事件中的脉动速度、平均速度变形率、雷诺应力等物理量的空间拓扑形态. 通过研究平均速度变形率各分量与雷诺应力各分量之间的空间相位差异,肯定了壁湍流相干结构雷诺应力各向异性复涡黏模型的合理性.