正弦波沟槽对湍流边界层相干结构影响的TR-PIV实验研究

利用高时间分辨率粒子图像测速(time-resolved particle image velocimetry, TR-PIV)技术,在不同雷诺数下对光滑壁面和二维顺流向、三维正弦波(two/three dimensional, 2D/3D)沟槽壁面湍流边界层流场进行了实验测量,从不同沟槽对湍流边界层相干结构影响的角度分析了其减阻的机理.对比不同壁面的各阶统计量结果发现:沟槽降低了壁面摩擦阻力,存在减阻效果,正弦波沟槽的减阻率增大.运用相关函数、λ_(ci)检测准则等方法提取了不同壁面湍流边界层发卡涡和低速条带等典型相干结构的空间拓扑形态,结果表明:两种沟槽壁面的相干结构在流向和法向上的空间尺度均有不同程度的减小,且相干结构与主流之间的倾角趋于更小,流体在法向上的运动及结构的抬升受到明显抑制,发卡涡诱导喷射和扫掠的能力降低,从而影响了湍流中能量与动量的输运过程及湍流的自维持机制,且相比于2D沟槽, 3D正弦波沟槽作用效果更为明显.在同一雷诺数下,随着距离壁面法向位置的增加,不同壁面湍流边界层低速条带的展向间距都变宽;但同一法向位置处2D/3D沟槽壁面湍流边界层低速条带的间距与光滑壁面的相比更宽,沟槽的存在有效抑制了低速条带在展向上的运动,使得低速条带更稳定.     

湍流边界层内流向涡的数值模拟研究

采用准二维共振三波作为湍流边界层近壁区相干结构初值,用直接数值模拟方法计算了流动从二维结构发展到三维结构并且伴随流向涡生成的整个过程,分析结果显示流向涡对湍流动能和质量传输有着重要作用,是湍流边界层相干结构的重要特征和运动形式.     

湍流边界层中下扫流与“反发卡涡”

用氢气泡法观测湍流边界层的下扫流和有关的流动结构.实验中发现一种新型涡结构，它的特征与典型的发卡涡正好相反.发卡涡的头部指向下游，而它的头部指向上游； 发卡涡的两腿之间，由于涡的诱导产生上升流，而它则在其两腿之间，由于涡的诱导产生下扫流. 关键词： 湍流边界层 流动显示 流动结构 发卡涡     

边界层内涡旋-颗粒相互作用的直接数值模拟

本文采用直接数值模拟方法结合拉格朗日跟踪方法研究平板边界层内涡旋与重颗粒的相互作用。模拟结果显示了由半球状粗糙元诱导产生的发夹涡沿流向的演化过程。小惯性的重颗粒在低压区上边缘的聚集阻碍了涡量的聚集并破坏了发夹涡的产生。模拟结果还证实颗粒在边界层内的输运主要受扫过和喷出事件的影响,且在边界层近壁区域内,扫过事件输运颗粒的能力强于喷出事件的。     

湍流边界层近壁区多个相干结构的数值模拟

从流动稳定性理论中的一般共振三波概念出发,提出一种湍流边界层近壁区多个相干结构的理论模型,采用高精度差分格式和Fourier谱展开相结合的方法,求解三维不可压Navier-Stokes方程,直接数值模拟近壁湍流多个相干结构的演化问题.并将得到的湍流边界层近壁区多个相干结构的数值演化特性与实验观察到的特性进行了比较.     

新型大涡数值模拟亚格子模型及应用

基于湍流大小尺度间动量输运的结构函数方程，提出了一种新的湍流大涡模型(LES)亚格子涡粘模式．新亚格子涡粘系数正比于纵向速度增量的扭率，它表征大小尺度湍流间的能量输运和耗散之比．新模式通过各向同性湍流直接数值模拟数据库的检验，并用于槽道湍流的大涡模拟计算，将所得结果与DNS结果进行了比较．     

小尺度涡流发生器强化传热机理的研究

在矩形槽道内布置了不同高度的小尺度涡流发生器,采用湍流模型对流体在其中的流动与传热特性进行了数值模拟,并对流体在槽道中的传热和流动特性进行了对比研究.分析了小尺度涡流发生器诱导涡的特性及其对湍流相干结构的影响与作用,讨论了湍流相干结构对温度场的作用机理,解释了涡流发生器强化换热的机理.     

近壁区湍流脉动压力相关项的理论模型

提出在湍流边界层近壁区采用三维波的理论模型描述湍流相干结构,根据理论模型计算了Ｒｅｙｎｏｌｄｓ应力输运方程中的脉动速度与脉动压力梯度的相关项,理论计算结果与直接数值模拟（ＤＮＳ）符合很好。表明该理论方法不仅有益于对湍流机制的了解,而且可能为湍流的近壁模型化开辟一条新的途径。     

湍流射流与扩散火焰大涡拟序结构的波动特性研究

１前言湍流射流扩散燃烧方式提供了射流火焰与涡团相互作用的最基本形式,对研究在湍流射流剪切边界层内的反应物的卷吸混合、热量及动量的输运及湍流射流结构方面具有其特殊的意义。Ｋａｔｔａ［１］利用浮力与单步反应机理模型对Ｎ２－Ｈ２射流扩散火焰进行了直接数值模...     

部分相干涡旋光束在大气湍流中的远场传输特性

运用广义惠更斯-菲涅耳原理,详细研究了部分相干涡旋光束在湍流介质中的远场传输规律.研究表明,部分相干涡旋光束的光谱相干度及光强分布与光束的拓扑电荷数、空间相对相干长度及湍流介质的折射率结构常数等因素有关.在湍流介质中,光谱相干度存在相位奇点,并且随着空间相对相干长度的增大,相干涡旋逐渐演化为强度涡旋,而湍流介质的强弱对部分相干涡旋光束的影响则相反,随着湍流介质的折射率常数的增大,强度涡旋逐渐演化为相干涡旋.     

Origin of effectiveness degradation in active drag reduction control of turbulent channel flow at Reτ = 1000

Direct numerical simulations of turbulent channel flows are performed with opposition control at Re_τ = 180 and 1000. The drag reduction rate at the higher Reynolds number is reduced by 25% compared with that at the lower Reynolds number. In order to investigate the reason for the degradation of the control effectiveness, we examine the response of Reynolds stresses and coherent structures in both the outer and inner regions to the control and the role that large-scale motions play therein. In the outer region, the Reynolds stresses at different length scales are reduced at the same rate as the drag reduction rate, and conditionally averaged large-scale motions with spanwise scale larger than half channel width are still large-scale low-speed streaks flanked by a pair of large-scale counter-rotating streamwise vortices but with reduced velocity amplitudes. In the inner region, the effectiveness of the control in suppressing the turbulence deteriorates at the higher Reynolds number. In response to the superimposition effect of large-scale motions, the contribution to near-wall wall-parallel velocity fluctuations from large-scale motions becomes larger at the higher Reynolds number, while the suppression of large-scale motions by the control is weaker than that of near-wall coherent structures. In both controlled and uncontrolled cases, large-scale motions can modulate the amplitudes of near-wall coherent structures, and the attenuation of streamwise vortices by the control under large-scale high-speed streaks is significantly less effective than that under large-scale low-speed streaks. As a result, the effectiveness of control in suppressing near-wall coherent structures and Reynolds shear stresses becomes weaker at the higher Reynolds number. The quantitative analysis of the contributions to the drag reduction rate from outer and inner regions shows that the effectiveness of the control is mainly determined by the suppression degree of near-wall motions. Furthermore, budgets of streamwise enstrophy are analysed to reveal the interaction of large-scale motions with near-wall streamwise vorticity. The titling effect induced by large-scale motions is positive under large-scale high-speed streaks, but negative under large-scale low-speed streaks, which could be a possible way of large-scale motion to modulate streamwise vorticity. In the controlled cases, the positive titling effect induced by large-scale motions under large-scale high-speed streaks is even enhanced, while other terms in the budgets are reduced, which could explain the degradation of control effectiveness in suppressing near-wall streamwise vortices under large-scale high-speed streaks. Therefore, the loss in the drag reduction rate at the higher Reynolds number is due to the weakened control effectiveness on near-wall coherent structures, which are exposed to the modulation effect of large-scale motions.     

Dynamics of hairpin vortices and polymer-induced turbulent drag reduction

It has been known for over six decades that the dissolution of minute amounts of high molecular weight polymers in wall-bounded turbulent flows results in a dramatic reduction in turbulent skin friction by up to 70%. First principles simulations of turbulent flow of model polymer solutions can predict the drag reduction (DR) phenomenon. However, the essential dynamical interactions between the coherent structures present in turbulent flows and polymer conformation field that lead to DR are poorly understood. We examine this connection via dynamical simulations that track the evolution of hairpin vortices, i.e., counter-rotating pairs of quasistreamwise vortices whose nonlinear autogeneration and growth, decay and breakup are centrally important to turbulence stress production. The results show that the autogeneration of new vortices is suppressed by the polymer stresses, thereby decreasing the turbulent drag.     

A new mixed subgrid-scale model for large eddy simulation of turbulent drag-reducing flows of viscoelastic fluids

A mixed subgrid-scale(SGS) model based on coherent structures and temporal approximate deconvolution(MCT) is proposed for turbulent drag-reducing flows of viscoelastic fluids. The main idea of the MCT SGS model is to perform spatial filtering for the momentum equation and temporal filtering for the conformation tensor transport equation of turbulent flow of viscoelastic fluid, respectively. The MCT model is suitable for large eddy simulation(LES) of turbulent dragreducing flows of viscoelastic fluids in engineering applications since the model parameters can be easily obtained. The LES of forced homogeneous isotropic turbulence(FHIT) with polymer additives and turbulent channel flow with surfactant additives based on MCT SGS model shows excellent agreements with direct numerical simulation(DNS) results. Compared with the LES results using the temporal approximate deconvolution model(TADM) for FHIT with polymer additives, this mixed SGS model MCT behaves better, regarding the enhancement of calculating parameters such as the Reynolds number.For scientific and engineering research, turbulent flows at high Reynolds numbers are expected, so the MCT model can be a more suitable model for the LES of turbulent drag-reducing flows of viscoelastic fluid with polymer or surfactant additives.     

Direct numerical simulation of a fully developed compressible wall turbulence over a wavy wall

Compressible turbulent channel flow over a wavy surface is investigated by direct numerical simulations using high-resolution finite difference schemes. The Reynolds number considered in the present paper is 3380 based on the bulk velocity, the channel half-width and the kinetic viscosity at the wall. Four test cases are simulated and analysed at Ma_m = 0.33, 0.8, 1.2, 1.5 based on the bulk velocity and the speed of sound at the wall. We mainly focus on the curvature and the Mach number effects on the compressible turbulent flows. Numerical results show that although the wavy wall has effects on the mean and fluctuation quantities, log law still exists in the distribution of the wave-averaged streamwise velocity if the roughness effects are taken into consideration in the scaling of it. Near-wall streaks are broken by the wavy surface and near-wall quasi-streamwise vortices mostly begin at the upslope of the wave and pass over the crest of it. The wavy wall makes the turbulence more active and the flow easier to be blended. From the viewpoint of turbulent kinetic budgets, curvature effects strengthen both the diffusion terms and the dissipation terms. At the same time, they change the properties of the compressibility-related terms and promote more inner energy transferring into turbulent kinetic energy. As the Mach number increases, the reattachment of the mean flow is delayed, which indicates the mean separation bubble becomes larger. Concerning the near-wall coherent structures, the vortices are more sparsely distributed with the increasing of the Mach number. For the supersonic cases, shock waves appear. Though they have little effects on the mean turbulent quantities, they change the structures of the flow fields and induce local separations at the upper wall of the channel.     

Scalar dispersion by coherent structures in uniformly sheared flow generated in a water tunnel

In order to investigate the role of coherent structures as mechanisms of scalar dispersion, we studied measurements of a passive scalar plume released in a uniformly sheared turbulent flow generated in a water tunnel. The flow had homogeneous turbulence properties in the measurement domain and contained hairpin vortices similar to those in boundary layers, and so was an ideal test bed to study the effects of coherent structures on turbulent dispersion, free from the effects of inhomogeneities or boundaries. Measurements of the velocity and concentration fields were acquired simultaneously using stereo particle image velocimetry and planar laser-induced fluorescence. We found that dye was preferentially located far away from vortices and was less likely to appear in close proximity to vortices, which is attributed to the high dissipation at the periphery of the vortices. However, we also found that dye was not directly correlated with the uniform momentum zones in the flow, suggesting a more complex relationship exists between these zones, the locations of vortices, and dye transport. Considering scalar flux events rather than simply the presence of dye as our condition of interest, a conditional eddy analysis demonstrated that hairpin vortices are responsible for the large scalar flux events as well as the large Reynolds stress events in the flow. The fact that the Reynolds stress was correlated with the scalar flux further confirmed that coherent structures are dominant mechanisms for scalar transport. Furthermore, we found that the scalar flux vector was preferentially inclined by 155° and ?25° with respect to the streamwise direction, and was thus approximately orthogonal to the planes of the legs of the most common upright and inverted hairpin structures in the flow. These findings demonstrate that coherent structures play an important and intricate role in turbulent diffusion.     

Numerical analysis of shock wave and supersonic turbulent boundary interaction between adiabatic and cold walls

Direct numerical simulations of shock wave and supersonic turbulent boundary layer interaction in a 24° compression ramp with adiabatic and cold-wall temperatures are conducted. The wall temperature effects on turbulence structures and shock motions are investigated. The results are validated against previous experimental and numerical data. The effects of wall cooling on boundary layer characteristics are analysed. Statistical data show that wall cooling has a significant effect on the logarithmic region of mean velocity profile downstream the interaction region. Moreover, the influence of wall temperature on Reynolds stress anisotropy is mainly limited in the near-wall region and has little change on the outer layer. As the wall temperature decreases, the streamwise coherency of streaks increases. Based on the analysis of instantaneous Lamb vector divergence, the momentum transport between small-scale vortices on cold-wall condition is significantly enhanced. In addition, spectral analysis of wall pressure signals indicates that the location of peak of low-frequency energy shifts toward higher frequencies in cold case. Furthermore, the dynamic mode decomposition results reveal two characteristic modes, namely a low-frequency mode exhibiting the breathing motion of separation bubble and a high-frequency mode associated with the propagation of instability waves above separation bubble. The shape of dynamic modes is not sensitive to wall temperature.     

Influence of polymer additives on turbulent energy cascading in forced homogeneous isotropic turbulence studied by direct numerical simulations

Direct numerical simulations(DNS) were performed for the forced homogeneous isotropic turbulence(FHIT) with/without polymer additives in order to elaborate the characteristics of the turbulent energy cascading influenced by drag-reducing effects.The finite elastic non-linear extensibility-Peterlin model(FENE-P) was used as the conformation tensor equation for the viscoelastic polymer solution.Detailed analyses of DNS data were carried out in this paper for the turbulence scaling law and the topological dynamics of FHIT as well as the important turbulent parameters,including turbulent kinetic energy spectra,enstrophy and strain,velocity structure function,small-scale intermittency,etc.A natural and straightforward definition for the drag reduction rate was also proposed for the drag-reducing FHIT based on the decrease degree of the turbulent kinetic energy.It was found that the turbulent energy cascading in the FHIT was greatly modified by the drag-reducing polymer additives.The enstrophy and the strain fields in the FHIT of the polymer solution were remarkably weakened as compared with their Newtonian counterparts.The small-scale vortices and the small-scale intermittency were all inhibited by the viscoelastic effects in the FHIT of the polymer solution.However,the scaling law in a fashion of extended self-similarity for the FHIT of the polymer solution,within the presently simulated range of Weissenberg numbers,had no distinct differences compared with that of the Newtonian fluid case.     

Mechanism of polymer drag reduction using a low-dimensional model

Using a retarded-motion expansion to describe the polymer stress, we derive a low-dimensional model to understand the effects of polymer elasticity on the self-sustaining process that maintains the coherent wavy streamwise vortical structures underlying wall-bounded turbulence. Our analysis shows that at small Weissenberg numbers, Wi, elasticity enhances the coherent structures. At higher Wi, however, polymer stresses suppress the streamwise vortices (rolls) by calming down the instability of the streaks that regenerates the rolls. We show that this behavior can be attributed to the nonmonotonic dependence of the biaxial extensional viscosity on Wi, and identify it as the key rheological property controlling drag reduction.     

近壁区理论耗散率k-ε模型的研究

提出在湍流边界层近壁区采用共振三波的理论模型描述湍流相干结构,根据理论模型计算了ε的分布。并且在传统k-ε模式基础上依照理论ε值计算了平均速度分布。在粘性作用层理论值与直接数值模拟符合很好。表明该理论方法不仅有益于对湍流机制的了解,而且可能为湍流的近壁模型化开辟一条新的途径。