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

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

Physics of multiple level hairpin vortex structures in turbulence

Previous experimental and numerical studies have revealed that the hairpin vortex is a basic flow element of transitional boundary layer. The hairpin vortex is believed to have legs, necks and a ring head. Based on our DNS study, the legs and the ring head are generated separately by different mechanisms. The legs function like an engine to generate low speed zones by rotation, create shear layers with surrounding high speed neighbor fluids, and further cause vortex ring formation through shear layer instability. In addition, the ring head is ?-shaped and separated from quasi-streamwise legs from the beginning. Contrary to the classical concept of "vortex breakdown", we believe transition from laminar flow to turbulence is a "buildup" process of multiple level vortical structures. The vortex rings of first level hairpins are mostly responsible for positive spikes, which cause new vorticity rollup, second level vortex leg formation and finally smaller second level vortex ring generation. The third and lower level vortices are generated following the same mechanism. In this paper, the physical process from ?-vortex to multi-level hairpin vortices is described in detail.     

Effects of dimples on laminar boundary layers

A series of direct numerical simulations of the flow past a flat plate with two and eight rows of dimples in a staggered arrangement is carried out. The Reynolds number based on the boundary layer thickness and freestream velocity near the inflow plane is 1000 and the dimples are spherical with a depth to diameter ratio of 0.1. The incoming flow is laminar and the boundary layer thickness before the dimples is half the dimple depth. At this low Reynolds number the flow is expected to remain laminar over a smooth flat plate. The presence of the dimples triggers instabilities that cause significant momentum transport. It is shown that the shear layer that forms as the flow separates over the first two rows of dimple becomes unstable and sheds coherent vortex sheets. The vortex sheets become unstable and are transformed into packets of horseshoe vortices. When these vortices evolve over a flat plate or over a series of dimples the flow dynamics are very different with important changes in momentum transport across the boundary layer.     

前缘曲率对三维边界层内被激发出非定常横流模态的影响研究

三维边界层感受性问题是三维边界层层流向湍流转捩的初始阶段,是实现三维边界层转捩预测与控制的关键环节.在高湍流度的环境下,非定常横流模态的失稳是导致三维边界层流动转捩的主要原因;但是,前缘曲率对三维边界层感受性机制作用的研究也是十分重要的课题之一.因此,本文采用直接数值模拟方法研究在自由来流湍流作用下具有不同椭圆形前缘三维（后掠翼平板）边界层内被激发出非定常横流模态的感受性机制;揭示不同椭圆形前缘曲率对三维边界层内被激发出非定常横流模态的扰动波波包传播速度、传播方向、分布规律、感受性系数以及分别提取获得一组扰动波的幅值、色散关系和增长率等关键因素的影响;建立在不同椭圆形前缘曲率情况下,三维边界层内被激发出非定常横流模态的感受性问题与自由来流湍流的强度和运动方向变化之间的内在联系;详细分析了不同强度各向异性的自由来流湍流在激发三维边界层感受性机制的物理过程中起着何种作用等.通过上述研究将有益于拓展和完善流动稳定性理论,为三维边界层内层流向湍流转捩的预测与控制提供依据.     

Convection and correlation of coherent structure in turbulent boundary layer using tomographic particle image velocimetry

The present experimental work focuses on a new model for space–time correlation and the scale-dependencies of convection velocity and sweep velocity in turbulent boundary layer over a flat wall. A turbulent boundary layer flow at Reθ= 2460 is measured by tomographic particle image velocimetry(tomographic PIV). It is demonstrated that arch, cane,and hairpin vortices are dominant in the logarithmic layer. Hairpins and hairpin packets are responsible for the elongated low-momentum zones observed in the instantaneous flow field. The conditionally-averaged coherent structures systemically illustrate the key roles of hairpin vortice in the turbulence dynamic events, such as ejection and sweep events and energy transport. The space–time correlations of instantaneous streamwise fluctuation velocity are calculated and confirm the new elliptic model for the space–time correlation instead of Taylor hypothesis. The convection velocities derived from the space–time correlation and conditionally-averaged method both suggest the scaling with the local mean velocity in the logarithmic layer. Convection velocity result based on Fourier decomposition(FD) shows stronger scale- dependency in the spanwise direction than in streamwise direction. Compared with FD, the proper orthogonal decomposition(POD) has a distinct distribution of convection velocity for the large- and small-scales which are separated in light of their contributions of turbulent kinetic energy.     

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

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

Extensive characterisation of a high Reynolds number decelerating boundary layer using advanced optical metrology

Over the last years, the observation of large-scale structures in turbulent boundary layer flows has stimulated intense experimental and numerical investigations. Nevertheless, partly due to the lack of comprehensive experimental data at sufficiently high Reynolds number, our understanding of turbulence near walls, especially in decelerating situations, is still quite limited. The aim of the present contribution is to combine the equipment and skills of several teams to perform a detailed characterisation of a large-scale turbulent boundary layer under adverse pressure gradient. Extensive particle image velocimetry (PIV) measurements are performed, including a set-up with 16 sCMOS cameras allowing the characterisation of the boundary layer on 3.5 m, stereo PIV and high resolution near wall measurements. In this paper, detailed statistics are presented and discussed, boundary conditions are carefully characterised, making this experiment a challenging test case for numerical simulation.     

确定分布的展向Lorentz力调制下的槽道湍流涡结构

采用直接数值模拟方法,对槽道湍流中确定分布的Lorentz力的流动控制与减阻问题进行研究.讨论了Lorentz力作用于槽道湍流后,流场的特性和涡结构的特性,并对此类Lorentz力对槽道湍流的控制与减阻机理进行了讨论.研究发现:1)Lorentz力诱导的层流流场壁面附近存在梯度极大的展向速度剪切层,该剪切层容易形成流向涡结构;2)在给定合适参数的确定分布的Lorentz力作用下,湍流流场仅剩周期分布的准流向涡;3)与未控制流场相比,控制后的流场中,准流向涡的抬升高度大大降低,从而减小猝发强度,使壁面阻力下降.     

A numerical analysis on flow field and heat transfer by interaction between a pair of vortices in rectangular channel flow

This paper is a numerical study of the effect of flow field and heat transfer created by interactions between a pair of vortices generated by a vortex generator in a rectangular channel flow. In order to analyze the vortices produced by the vortex generator, the pseudo-compressibility method is introduced into the Navier–Strokes (NS) equation of a three-dimensional unsteady, incompressible viscous flow. A two-layer k–ε turbulence model is used on the flat plate three-dimensional turbulence boundary to predict the turbulence characteristics of the vortices. The computational results accurately predict the vortex characteristics, which are related to Reynolds stress, turbulent kinetic energy, and flow field. Also, in the prediction of thermal boundary layers, skin friction characteristics, and heat transfers, the present results are reasonably close to the experimental results obtained by other researchers.     

Constrained large-eddy simulation of supersonic turbulent boundary layer over a compression ramp

A supersonic turbulent boundary layer over a compression ramp is numerically investigated using the constrained large-eddy simulation (CLES) method. The compression corner is characterised by a deflection angle of 24°. The free-stream Mach number is Ma_∞ = 2.9, and the Reynolds number based on the momentum thickness of inlet boundary layer is Re_θ = 2300. The mean and statistical quantities, such as mean velocity, wall pressure and Reynolds stresses, are thoroughly analysed and compared with those from traditional large-eddy simulation (LES), experimental measurement and direct numerical simulation (DNS). It turns out that CLES can predict the friction coefficient, wall-pressure distribution, size of separation bubble, Reynolds stresses, etc. more accurately than traditional LES, and the results are in reasonable agreement with the experimental and/or DNS data. Also discussed are the effects of specific parameterisations of the Reynolds constraint and interfacial positions separating the constrained and unconstrained regions on the performance of the CLES method.     

Spatial resolution analysis of planar PIV measurements to characterise vortices in turbulent flows

The effects of spatial resolution of planar particle image velocimetry (PIV) on vortex size, swirling strength, circulation and population density characterisation are analysed using a series of experimental and numerical databases. The databases comprise a PIV database of an adverse-pressure-gradient turbulent boundary layer (APG TBL), a PIV database of a zero-pressure-gradient (ZPG) TBL in streamwise-wall-normal planes and streamwise-wall-normal slices of a direct numerical simulation (DNS) of a ZPG TBL. The effects of interrogation window and mesh sizes on the vortex parameters are analysed in the outer region of these flows using different qualitative and quantitative approaches. The quantitative analysis mainly capitalises on the possibility of mimicking the PIV data-sets with the DNS one. These approaches allow us to not only isolate the effects of mesh size and the interrogation window size but also to deduce the combined effects of other measurement errors in PIV. Typical values of mesh size and interrogation window size (0.01–0.03 of the boundary layer thickness) and typical levels of measurement uncertainties have significant effects on the vortex parameters. Moreover, each PIV error source affects the vortex parameters in different and frequently opposite manners. Hence, an optimal selection of measurement parameters such as the interrogation window size is indispensable in order to minimise the effects of spatial resolution and other measurement errors on the vortex parameters. Guidelines are presented in the Conclusions section of this paper. Finally, it is found that all the vortex parameters, when averaged across the outer region, are reasonably comparable in the ZPG and APG TBLs despite the fact that these are very different flows.     

Decaying two-dimensional turbulence in a circular container

We present direct numerical simulations of two-dimensional decaying turbulence at initial Reynolds number 5 x 10(4) in a circular container with no-slip boundary conditions. Starting with random initial conditions the flow rapidly exhibits self-organization into coherent vortices. We study their formation and the role of the viscous boundary layer on the production and decay of integral quantities. The no-slip wall produces vortices which are injected into the bulk flow and tend to compensate the enstrophy dissipation. The self-organization of the flow is reflected by the transition of the initially Gaussian vorticity probability density function (PDF) towards a distribution with exponential tails. Because of the presence of coherent vortices the pressure PDF become strongly skewed with exponential tails for negative values.     

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

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

超声速湍流边界层湍流统计与湍流结构分析

可压缩边界层转捩问题与湍流问题一直是制约高超声速飞行器发展的关键基础问题,也是近年来流体力学领域研究的热点问题.采用直接数值模拟方法,获得了空间发展的Ma=2.25超声速湍流边界层流场,通过对湍流边界层的发展状态进行评估,得出有效的Reynolds数Reθ 范围约为2600～4600.对壁面摩阻系数开展了分解,获得了各...     

Direct numerical simulation of the flow around a wall-mounted square cylinder under various inflow conditions

The flow around a wall-mounted square cylinder of side d is investigated by means of direct numerical simulation (DNS). The effect of inflow conditions is assessed by considering two different cases with matching momentum-thickness Reynolds numbers Re_θ ? 1000 at the obstacle: the first case is a fullyturbulent zero pressure gradient boundary layer, and the second one is a laminar boundary layer with prescribed Blasius inflow profile further upstream. An auxiliary simulation carried out with the pseudo-spectral Fourier–Chebyshev code SIMSON is used to obtain the turbulent time-dependent inflow conditions which are then fed into the main simulation where the actual flow around the cylinder is computed. This main simulation is performed, for both laminar and turbulent-inflows, with the spectral-element method code Nek5000. In both cases the wake is completely turbulent, and we find the same Strouhal number St ? 0.1, although the two wakes exhibit structural differences for x > 3d downstream of the cylinder. Transition to turbulence is observed in the laminar-inflow case, induced by the recirculation bubble produced upstream of the obstacle, and in the turbulent-inflow simulation the streamwise fluctuations modulate the horseshoe vortex. The wake obtained in our laminar-inflow case is in closer agreement with reference particle image velocimetry measurements of the same geometry, revealing that the experimental boundary layer was not fully turbulent in that dataset, and highlighting the usefulness of DNS to assess the quality of experimental inflow conditions.     

Anisotropic Reynolds stress tensor representation in shear flows using DNS and experimental data

Tensorial decompositions and projections are used to study the performance of algebraic non-linear models and predict the anisotropy of the Reynolds stresses. Direct numerical simulation (DNS) data for plane channel flows at friction Reynolds number (Re_τ = 180, 395, 590, 1000), and for the boundary layer using both DNS (Re_τ = 359, 830, 1271) and experimental data (Re_τ = 2680, 3891, 4941, 7164) are used to build and evaluate the models. These data are projected into tensorial basis formed from the symmetric part of mean velocity gradient and non-persistence-of-straining tensor. Six models are proposed and their performances are investigated. The scalar coefficients for these six different levels of approximations of the Reynolds stress tensor are derived, and made dimensionless using the classical turbulent scales, the kinetic turbulent energy (κ) and its dissipation rate (ε). The dimensionless coefficients are then coupled with classical wall functions. One model is selected by comparing the predicted Reynolds stress components with experimental and DNS data, presenting a good prediction for the shear and normal Reynolds stresses.     

Three-Dimensional Stability of Burgers Vortices

Burgers vortices are explicit stationary solutions of the Navier-Stokes equations which are often used to describe the vortex tubes observed in numerical simulations of three-dimensional turbulence. In this model, the velocity field is a two-dimensional perturbation of a linear straining flow with axial symmetry. The only free parameter is the Reynolds number Re = Γ/ν, where Γ is the total circulation of the vortex and ν is the kinematic viscosity. The purpose of this paper is to show that Burgers vortices are asymptotically stable with respect to small three-dimensional perturbations, for all values of the Reynolds number. This general result subsumes earlier studies by various authors, which were either restricted to small Reynolds numbers or to two-dimensional perturbations. Our proof relies on the fact that the linearized operator at Burgers vortex has a simple and very specific dependence upon the axial variable. This allows to reduce the full linearized equations to a vectorial two-dimensional problem, which can be treated using an extension of the techniques developed in earlier works. Although Burgers vortices are found to be stable for all Reynolds numbers, the proof indicates that perturbations may undergo an important transient amplification if Re is large, a phenomenon that was indeed observed in numerical simulations.     

Numerical investigation of separation-induced transition in a lowpressure turbine cascade in a low-disturbance environment

In this paper,the separation-induced transition in an LPT(low-pressure turbine)cascade is investigated at low Reynolds number with DNS(direct numerical simulation).The transition process is accurately predicted giving good agreements between the DNS and experimental results.To illustrate the secondary instability of separation-induced transition in a low-disturbance environment,the results are comprehensively analyzed in both Fourier space and physical space.It is illustrated that the effect of hyperbolic instability dominates around the saddle point of hyperbolic streamlines.This instability mechanism is responsible for the emergence of the streamwise vortices in the braid region.Elongated and intensified because of the“stretching”effect of the background flow,these vortices become the most noticeable characteristic of the flow field.Fundamental modes of small spanwise wavelength are excited in the braid region,so as some low-frequency modes.The elliptical instability plays a minor role than hyperbolic instability.It is also observed that the fundamental mode with a larger spanwise wavelength is unstable in the vortex core which is associated with the deformation of the vortex core via elliptical instability.There is no convincing evidence for the existence of subharmonic instability.     

A study on drag reduction of a rotationally oscillating circular cylinder at low Reynolds number

The flow field around a rotationally oscillating circular cylinder in a uniform flow is studied by using a particle image velocimetry to understand the mechanism of drag reduction and the corresponding suppression of vortex shedding in the cylinder wake at low Reynolds number. Experiments are conducted on the flow around the circular cylinder under rotational oscillation at forcing Strouhal number 1, rotational amplitude 2 and Reynolds number 2,000. It is found from the flow measurement by PIV that the width of the wake is narrowed and the velocity fluctuations are reduced by the rotational oscillation of the cylinder, which results in the drag reduction rate of 30%. The mechanism of drag reduction is studied by phase-averaged PIV measurement, which indicates the formation of periodic small-scale vortices from both sides of the cylinder. It is found from the cross-correlation measurement between the velocity fluctuations that the large-scale structure of vortex shedding is almost removed in the cylinder wake, when the small-scale vortices are generated at the unstable frequency of shear layer by the influence of rotational oscillation.