高聚物减阻溶液对壁湍流输运过程的影响

基于相同雷诺数下清水和高分子聚合物溶液壁湍流的高时间分辨率粒子图像测速技术(time-resolved particle image velocimetry, TRPIV)的对比实验, 从高聚物溶液对湍流边界层动量能量输运影响的角度分析其减阻的机理. 对比两者的雷诺应力发现高聚物的存在抑制了湍流输运过程. 这一影响与高聚物对壁湍流中占主导地位的涡旋运动和低速条带等相干结构的作用密切相关. 运用条件相位平均、相关函数和线性随机估计(linear stochastic estimation, LSE)等方法, 分析提取了高聚物溶液流场中的发卡涡和发卡涡包等典型相干结构的空间拓扑形态. 相比于清水, 高聚物溶液中相干结构的流向尺度增大, 涡旋运动的发展及低速流体喷射的强度受到削弱, 表明了添加的高聚物阻碍了湍流原有的能量传递和自维持的机理. 正是通过影响相干结构, 高聚物抑制了湍流边界层中近壁区与外区之间的动量和能量输运, 使得湍流的无序性降低, 从而减小了湍流流动的阻力.     

横向粗糙元壁面槽道湍流的大涡模拟研究

采用大涡模拟和浸没边界法相结合对不同高度和不同间距横向粗糙元壁面槽道湍流进行了模拟,得到了光滑壁面和粗糙壁面湍流的流向平均速度分布,雷诺剪切应力,脉动速度均方根和近壁区拟序结构。结果发现横向粗糙元降低了流向平均速度,增大了流动阻力,粗糙壁面湍流的雷诺剪切应力大于光滑壁面。粗糙元降低了流向脉动速度,增强了展向和法向脉动速度。粗糙元高度越高,对湍流流动影响越大,而粗糙元间距对湍流统计特性的影响不大。粗糙壁面仍然存在着和光滑壁面类似的条带结构。     

沟槽壁面减阻机理实验研究

利用IFA300型热线风速仪,测量了光滑壁面和沟槽减阻壁面湍流边界层内的瞬时速度,利用自行设计的阻力天平仪测量了壁面摩擦力。得到了边界层无量纲速度分布和平均湍动能分布。对测得的脉动速度信号,利用离散正交小波变换按时间和尺度分解,得到各尺度分量的湍动能,并且发现其分布在湍流惯性区具有极大值。分析表明,当沟槽有减阻效果时,边界层内的平均湍动能减小,湍流惯性区各分量的湍动能极大值亦减小。     

不同粒径粗糙元对壁面湍流拟序结构的影响

采用氢气泡流动显示技术,研究了不同粒径粗糙元对壁湍流拟序结构的影响。实验中粗糙元布置在氢泡丝前,均匀排成一排,基于平均速度和水力直径的雷诺数从14300到48000变化,离散粗糙元直径分别为0.9、1.6、2.3、4、5和6mm,粗糙元间距为2 cm。实验发现:光滑壁面和粗糙元壁面近壁湍流条带数量和条带高度都随雷诺数增大而增加。相同雷诺数下,粗糙元壁面流动产生的条带数量和条带高度随粗糙元粒径增大而增加。粗糙元壁面湍流条带数量和条带高度均比光滑壁面大。该研究对壁湍流的流动控制工程应用具有参考价值。     

不同排列粗糙元对湍流拟序结构的影响

采用氢气泡流动显示技术,研究了壁面不同排列粗糙元对近壁湍流拟序结构的影响。实验中基于平均速度和水力直径的雷诺数分别为14300、31200、48000。氢泡丝上游布置的离散粗糙元直径为2．3、4和6mm,得到了光滑壁面和5种不同排列方式的离散粗糙元壁面湍流条带和湍流斑块特征。结果发现：相同雷诺数和相同排列方式下条带间距...     

八重准周期排列的短沟槽结构减阻机理分析

设计构建了一排和三排阵列的八重准周期短条纹沟槽减阻结构,以及作为对比研究的无序和周期结构,并采用雷诺Navier-Stokes方程和RANGκ-ε湍流模型,系统计算了这些结构表面的湍流边界层状态和总应力,模拟结果显示:八重准周期沟槽结构相对于周期和无序结构具有更优的减阻效应,且为三排阵列时的减阻效果明显优于单排阵列,这一结果得到了减阻实验的验证,通过分析比较不同结构的流体边界层特性发现,八重准周期结构可有效抑制附面层的涡强度,减小湍流耗散速率,保持流体条纹相的稳定性,结合二维光栅的夫琅禾费衍射波模型分析表明,八重准周期结构可减弱衍射谱在大角度方向上的谱强度,揭示出该结构抑制流体相干扰动波扩展的物理机制,并与流场分析结果相符合。     

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

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

粗糙壁面湍流边界层流动和拟序结构的大涡模拟研究

采用LES方法和虚拟粗糙元模型,,模拟了Rer=180有粗糙壁面的充分发展槽道湍流流动.从流向平均速度、脉动速度均方根、雷诺剪切应力、脉动速度相关,以及瞬时近壁准流向涡结构等多个方面,对粗糙壁面湍流边界层流动和拟序结构进行了分析.结果表明,粗糙壁面处雷诺剪切应力、脉动速度均方根以及近壁条带平均距离和准流向涡尺度都大于光滑壁面,且随着粗糙元高度的增加这种趋势更加明显.本工作为进一步研究颗粒在近壁区域的弥散规律奠定了基础.     

壁面换热和压力梯度条件下边界层演化的LES研究

本文采用大涡模拟(LES)方法对壁面换热边界条件下有压力梯度的分离边界层流动问题进行了数值模拟研究,对不同壁面热边界条件下分离边界层进行了统计特性及瞬态流场分析.结果表明:与绝热壁面条件相比,等温冷壁面条件(0.8倍来流总温)下附着边界层内的速度剖面比绝热壁面更加饱满,边界层形状因子降低、分离泡也被显著抑制,其时均分离泡尺寸无论是在流向还是法向均更小;等温冷壁面条件下分离剪切层展向涡脱落的位置提前,展向涡卷起的频率降低,分离剪切层内相干结构尺度整体减小,促使分离边界层转捩加速和再附位置前移。     

基于蚯蚓背孔射流的仿生射流表面减阻性能研究

为了减小流体对固体壁面的阻力, 基于蚯蚓生物学特征, 对蚯蚓背孔射流特性进行分析, 建立仿蚯蚓背孔射流的仿生射流表面计算模型, 采用SST k-ω 湍流模型对仿生射流表面的减阻特性进行数值模拟, 同时对数值模拟结果进行实验验证, 并以此研究了仿蚯蚓背孔射流表面的减阻机理.结果表明, 在一定条件下, 仿蚯蚓背孔射流的仿生射流表面具有较好的减阻效果; 在同一射流方向角下, 随着射流速度的增加, 减阻率逐渐增大; 在同一射流速度下, 随着射流方向角的增加, 减阻率呈先减小后增大的变化趋势; 数值模拟与实验均在射流速度为1 m·s^-1、射流方向角为-30°时达到最大, 分别为8.69%, 7.86%; 射流表面改变了原有光滑壁面的边界层结构, 对壁面边界层进行了有效的控制, 减小了壁面的剪应力, 降低了壁面边界层的速度.     

Coherent structures over riblets in turbulent boundary layer studied by combining time-resolved particle image velocimetry(TRPIV),proper orthogonal decomposition(POD),and finite-time Lyapunov exponent(FTLE)

Time-resolved particle image velocimetry(TRPIV) experiments are performed to investigate the coherent structure's performance of riblets in a turbulent boundary layer(TBL) at a friction Reynolds number of 185. To visualize the energetic large-scale coherent structures(CSs) over a smooth surface and riblets, the proper orthogonal decomposition(POD) and finite-time Lyapunov exponent(FTLE) are used to identify the CSs in the TBL. Spatial-temporal correlation is implemented to obtain the characters and transport properties of typical CSs in the FTLE fields. The results demonstrate that the generic flow structures, such as hairpin-like vortices, are also observed in the boundary layer flow over the riblets, consistent with its smooth counterpart. Low-order POD modes are more sensitive to the riblets in comparison with the high-order ones,and the wall-normal movement of the most energy-containing structures are suppressed over riblets. The spatial correlation analysis of the FTLE fields indicates that the evolution process of the hairpin vortex over riblets are inhibited. An apparent decrease of the convection velocity over riblets is noted, which is believed to reduce the ejection/sweep motions associated with high shear stress from the viscous sublayer. These reductions exhibit inhibition of momentum transfer among the structures near the wall in the TBL flows.     

Effects of single synthetic jet on turbulent boundary layer

The turbulent boundary layer (TBL) is actively controlled by the synthetic jet generated from a circular hole. According to the datasets of velocity fields acquired by a time-resolved particle image velocimetry (TR-PIV) system, the average drag reduction rate of 6.2% in the downstream direction of the hole is obtained with control. The results of phase averaging show that the synthetic jet generates one vortex pair each period and the consequent vortex evolves into hairpin vortex in the environment with free-stream, while the reverse vortex decays rapidly. From the statistical average, it can be found that a low-speed streak is generated downstream. Induced by the two vortex legs, the fluid under them converges to the middle. The drag reduction effect produced by the synthetic jet is local, and it reaches a maximum value at x⁺=400, where the drag reduction rate reaches about 12.2%. After the extraction of coherent structure from the spatial two-point correlation analysis, it can be seen that the synthetic jet suppresses the streamwise scale and wall-normal scale of the large scale coherent structure, and slightly weakens the spanwise motion to achieve the effect of drag reduction.     

Observation of vortex packets in direct numerical simulation of fully turbulent channel flow

A number of experimental studies have inferred the existence of packets of inclined, hairpinlike vortices in wall turbulence on the basis of observations made in two-dimensional x−y planes using visualization and particle image velocimetry (PIV). However, there are very few observations of hairpins in existing three-dimensional studies made using direct numerical simulation (DNS), and no such study claims to have revealed packets. We demonstrate, for the first time, the existence of hairpin vortex packets in DNS of turbulent flow. The vortex packet structure found in the present study at low Reynolds number,Re _t=300, is consistent with and substantiates the observations and the results from twodimensional PIV measurements at higher Reynolds numbers in channel, pipe and boundary layer flows. Thus, the evidence supports the view that vortex packets are a universal feature of wall turbulence, independent of effects due to boundary layer trips or critical conditions in the aforementioned numerical studies. Visualization of the DNS velocity field and vortices also shows the close association of hairpin packets with long low-momentum streaks and the regions of high Reynolds shear stress.     

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.     

A Reynolds stress model for turbulent flows of viscoelastic fluids

A second-order closure is developed for predicting turbulent flows of viscoelastic fluids described by a modified generalised Newtonian fluid model incorporating a nonlinear viscosity that depends on a strain-hardening Trouton ratio as a means to handle some of the effects of viscoelasticity upon turbulent flows. Its performance is assessed by comparing its predictions for fully developed turbulent pipe flow with experimental data for four different dilute polymeric solutions and also with two sets of direct numerical simulation data for fluids theoretically described by the finitely extensible nonlinear elastic – Peterlin model. The model is based on a Newtonian Reynolds stress closure to predict Newtonian fluid flows, which incorporates low Reynolds number damping functions to properly deal with wall effects and to provide the capability to handle fluid viscoelasticity more effectively. This new turbulence model was able to capture well the drag reduction of various viscoelastic fluids over a wide range of Reynolds numbers and performed better than previously developed models for the same type of constitutive equation, even if the streamwise and wall-normal turbulence intensities were underpredicted.     

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.