基于壁面主动变形的湍流减阻控制研究

采用谱方法, 对反向控制下壁面主动变形的槽道湍流进行了直接数值模 拟研究. 结果表明, 在壁面最大变形量小于5倍黏性尺度条件下, 压差阻力可略, 摩擦阻力 降低7.6%. 施加控制后, 湍流强度和雷诺应力受到明显抑制, 平均速度剖面对数区上移. 受壁面法向运动的影响, 条带结构强度减弱、尺度变大; 流向涡外移且强度减弱, 其倾斜和 抬起的角度均有不同程度的减小. 壁面变形呈现流向拉长的凹槽结构, 其平均间距 为90倍黏性尺度.     

Coherent structures and riblets

This work deals with the effect of the riblets on the coherent structures near the wall. The emphasis is put on the genesis of the quasi-streamwise vortices in the presence of the riblets. The quasi-streamwise vortices regenerate by the tilting of wall normal vorticity induced by prevailing structures. This requires a mechanism which leads to a temporal streamwise dependence near the elongated flow structures and to a subsequent formation of new wall normal vorticity. It is suggested here that the action of existing quasi-streamwise vortices on the sidewalls of wall normal vorticity may create a local, streamwise dependent spanwise velocity and therefore, a secondary wall normal vorticity field. A preliminary analysis of the set-up and the time and space development of this secondary three-dimensional flow associated with the regeneration mechanism, is given. An attempt is made, in order to explain the drag reduction performed by the riblets through an intermittent model, based on the protrusion height. Logical estimates of the amount of drag reduction are obtained. The differences between the mechanism suggested here and those based on forced control experiments are also discussed.     

Swirling-strength based large eddy simulation of turbulent flow around single square cylinder at low Reynolds numbers

In view of the fact that large scale vortices play the substantial role of momentum transport in turbulent flows, large eddy simulation(LES) is considered as a better simulation model. However, the sub-grid scale(SGS) models reported so far have not ascertained under what flow conditions the LES can lapse into the direct numerical simulation. To overcome this discrepancy, this paper develops a swirling strength based the SGS model to properly model the turbulence intermittency, with the primary characteristics that when the local swirling strength is zero, the local sub-grid viscosity will be vanished. In this paper, the model is used to investigate the flow characteristics of zero-incident incompressible turbulent flows around a single square cylinder(SC)at a low Reynolds number range Re ∈ [103, 104]. The flow characteristics investigated include the Reynolds number dependence of lift and drag coefficients, the distributions of time-spanwise averaged variables such as the sub-grid viscosity and the logarithm of Kolmogorov micro-scale to the base of 10 at Re = 2 500 and 104, the contours of spanwise and streamwise vorticity components at t = 170. It is revealed that the peak value of sub-grid viscosity ratio and its root mean square(RMS) values grow with the Reynolds number. The dissipation rate of turbulent kinetic energy is larger near the SC solid walls.The instantaneous factor of swirling strength intermittency(FSI) exhibits some laminated structure involved with vortex shedding.     

Drag reduction of turbulent channel flows over an anisotropic porous wall with reduced spanwise permeability

The direct numerical simulation(DNS) is carried out for the incompressible viscous turbulent flows over an anisotropic porous wall. Effects of the anisotropic porous wall on turbulence modifications as well as on the turbulent drag reduction are investigated. The simulation is carried out at a friction Reynolds number of 180, which is based on the averaged friction velocity at the interface between the porous medium and the clear fluid domain. The depth of the porous layer ranges from 0.9 to 54 viscous units. The permeability in the spanwise direction is set to be lower than the other directions in the present simulation. The maximum drag reduction obtained is about 15.3% which occurs for a depth of 9 viscous units. The increasing of drag is addressed when the depth of the porous layer is more than 25 wall units. The thinner porous layer restricts the spanwise extension of the streamwise vortices which suppresses the bursting events near the wall. However, for the thicker porous layer, the wall-normal fluctuations are enhanced due to the weakening of the wall-blocking effect which can trigger strong turbulent structures near the wall.     

Cross-Stream Vorticity Field Induced by Streamwise Vortices in Unbounded and Wall-Bounded Shear Flows

This computational study examines the unsteady cross-stream vorticity structures that form when one or more streamwise vortices are immersed in homogeneous and boundary-layer shear flows. A quasi-two-dimensional limit is considered in which the velocity and vorticity fields, while still possessing three nonzero components, have vanishing gradient in the streamwise direction. This idealization is suitable to applications such as streamwise vortices that occur along a ship hull or airplane fuselage and it can be used as an idealized representation of the quasi-streamwise vortices in the near-wall region of a turbulent boundary layer. In this quasi-two-dimensional idealization, the streamwise velocity has no effect on the cross-stream velocity associated with the vortex. However, the vortex acts to modify the cross-stream vorticity component, resulting in regions of the flow with strong deviations in streamwise velocity. This paper examines the complex structures that form as the cross-stream vorticity field is wrapped up by the vortex and the effect of these structures on the streamwise velocity field, first for vortices immersed in homogeneous shear flow and then for vortices immersed in a boundary layer along a flat wall. Received 2 January 2002 and accepted 13 August 2002 Published online 3 December 2002 RID="*" ID="*" This project was supported by the Office of Naval Research under Grant Number N00014-01-1-0015. Dr. Thomas Swain is the program manager. Communicated by T.B. Gatski     

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.     

On streamwise vortices in Large Eddy Simulations of initially laminar plane mixing layers

This paper details the influence of the nature of imposed inflow fluctuations on Large Eddy Simulations of a spatially developing turbulent mixing layer originating from laminar boundary layers. A simulation with imposed white-noise random fluctuations, commonly used in numerical simulations, produces mean-flow statistics that agree well with reference experimental data. Whilst flow visualisation images show evidence for streamwise vorticity in this simulation, quantitative statistics do not reveal the presence of statistically stationary streamwise vortices. A further simulation that uses physically-correlated inflow fluctuations also produces good mean-flow statistical agreement with reference data. From secondary shear stress contours it can be inferred that this simulation does, however, predict the presence of statistically stationary streamwise vortices. The properties of the streamwise vortices are in good agreement with experimental data. The data presented here indicate that, even for initially laminar conditions, plane mixing layer simulations require accurate physically correlated inflow conditions in order to reproduce the flow features found experimentally.     

Experiments in Turbulent Pipe Flow with Polymer Additives at Maximum Drag Reduction

In this paper we report on (two-component) LDV experiments in a fully developed turbulent pipe flow with a drag-reducing polymer (partially hydrolyzed polyacrylamide) dissolved in water. The Reynolds number based on the mean velocity, the pipe diameter and the local viscosity at the wall is approximately 10000. We have used polymer solutions with three different concentrations which have been chosen such that maximum drag reduction occurs. The amount of drag reduction found is 60–70%. Our experimental results are compared with results obtained with water and with a very dilute solution which exhibits only a small amount of drag reduction. We have focused on the observation of turbulence statistics (mean velocities and turbulence intensities) and on the various contributions to the total shear stress. The latter consists of a turbulent, a solvent (viscous) and a polymeric part. The polymers are found to contribute significantly to the total stress. With respect to the mean velocity profile we find a thickening of the buffer layer and an increase in the slope of the logarithmic profile. With respect to the turbulence statistics we find for the streamwise velocity fluctuations an increase of the root mean square at low polymer concentration but a return to values comparable to those for water at higher concentrations. The root mean square of the normal velocity fluctuations shows a strong decrease. Also the Reynolds (turbulent) shear stress and the correlation coefficient between the stream wise and the normal components are drastically reduced over the entire pipe diameter. In all cases the Reynolds stress stays definitely non-zero at maximum drag reduction. The consequence of the drop of the Reynolds stress is a large polymer stress, which can be 60% of the total stress. The kinetic-energy balance of the mean flow shows a large transfer of energy directly to the polymers instead of the route by turbulence. The kinetic energy of the turbulence suggests a possibly negative polymeric dissipation of turbulent energy. This revised version was published online in July 2006 with corrections to the Cover Date.     

Riblet drag reduction in mild adverse pressure gradients: A numerical investigation

Riblet films are a passive method of turbulent boundary layer control that can reduce viscous drag. They have been studied with great detail for over 30 years. Although common riblet applications include flows with Adverse Pressure Gradients (APG), nearly all research thus far has been performed in channel flows. Recent research has provided motivation to study riblets in more complicated turbulent flows with claims that riblet drag reduction can double in mild APG common to airfoils at moderate angles of attack. Therefore, in this study, we compare drag reduction by scalloped riblet films between riblets in a zero pressure gradient and those in a mild APG using high-resolution large eddy simulations. In order to gain a fundamental understanding of the relationship between drag reduction and pressure gradient, we simulated several different riblet sizes that encompassed a broad range of s⁺ (riblet width in wall units), similarly to many previously published experimental studies. We found that there was only a slight improvement in drag reduction for riblets in the mild APG. We also observed that peak values of streamwise turbulence intensity, turbulent kinetic energy, and streamwise vorticity scale with riblet width. Primary Reynolds shear stresses and turbulence kinetic energy production however scale with the ability of the riblet to reduce skin-friction.     

Direct simulation of breakdown to turbulence following oblique instability waves in a supersonic boundary layer

The late stages of transition to turbulence in a Mach two boundary layer are investigated by direct numerical simulation of the compressible Navier-Stokes equations. The primary instability at this Mach number consists of oblique waves, which are known to form a pattern of quasi-streamwise vortices. It is found that breakdown does not follow immediately from these vortices, which decay in intensity. The generation of new vortices is observed by following the evolution of the pressure and vorticity in the simulation, and analysed by consideration of vorticity stretching. It is found that the slight inclined and skewed nature of the quasi-streamwise vortices leads to a production of oppositely signed streamwise vorticity, which serves as a strong localised forcing of the shear layer alongside the original vortices, formed by convection and stretching of spanwise vorticity. The shear layer rolls up into many new vortices, and is followed by a sharp increase in the energy of higher frequencies and in the skin friction.     

Relationship between wall shear stresses and streamwise vortices in turbulent flows over wavy boundaries

The relationship between wall shear stresses and near-wall streamwise vortices is investigated via a direct numerical simulation(DNS) of turbulent flows over a wavy boundary with traveling-wave motion. The results indicate that the wall shear stresses are still closely related to the near-wall streamwise vortices in the presence of a wave. The wave age and wave phase significantly affect the distribution of a two-point correlation coefficient between the wall shear stresses and streamwise vorticity. For the slow wave case of c/Um = 0.14, the correlation is attenuated above the leeward side while the distribution of correlation function is more elongated and also exhibits a larger vertical extent above the crest. With respect to the fast wave case of c/U_m=1.4, the distribution of the correlation function is recovered in a manner similar to that in the flat-wall case. In this case, the maximum correlation coefficient exhibits only slight differences at different wave phases while the vertical distribution of the correlation function depends on the wave phase.     

Time resolved flow-field measurements of a turbulent mixing layer over a rectangular cavity

High Reynolds number, low Mach number, turbulent shear flow past a rectangular, shallow cavity has been experimentally investigated with the use of dual-camera cinematographic particle image velocimetry (CPIV). The CPIV had a 3 kHz sampling rate, which was sufficient to monitor the time evolution of large-scale vortices as they formed, evolved downstream and impinged on the downstream cavity wall. The time-averaged flow properties (velocity and vorticity fields, streamwise velocity profiles and momentum and vorticity thickness) were in agreement with previous cavity flow studies under similar operating conditions. The time-resolved results show that the separated shear layer quickly rolled-up and formed eddies immediately downstream of the separation point. The vortices convect downstream at approximately half the free-stream speed. Vorticity strength intermittency as the structures approach the downstream edge suggests an increase in the three-dimensionality of the flow. Time-resolved correlations reveal that the in-plane coherence of the vortices decays within 2–3 structure diameters, and quasi-periodic flow features are present with a vortex passage frequency of ~1 kHz. The power spectra of the vertical velocity fluctuations within the shear layer revealed a peak at a non-dimensional frequency corresponding to that predicted using linear, inviscid instability theory.     

纵向涡对近壁湍流影响的激光测试

用激光测试方法详细地研究了单涡、上洗双涡、下洗双涡对近壁湍流的影响．在涡的下洗区，流体流速高，湍流度低，流体流向壁面；在涡的上洗区，流体流速低，湍流度高，流体远离壁面     

Spatial characterization of the numerically simulated vorticity fields of a flow in a flume

The topology of large scale structures in a turbulent boundary layer is investigated numerically. Spatial characteristics of the large scale structure are presented through an original method, proper orthogonal decomposition (POD) of the three-dimensional vorticity fields. The DNS results, obtained by Tiselj et al. [23] for a fully developed turbulent flow in a flume, are used in the present work to analyze coherent structures with the proposed methodology. In contrast to the reconstruction methods that use instantaneous flow quantities, this approach utilizes the whole dataset of the numerical simulation. The analysis uses one thousand 3D vorticity fields from 50000 time steps of the simulation for the Reynolds number of 2600 (the turbulent Reynolds number Re^*=171). The computational domain is 2146×171×537 wall units and the grid resolution is 128×65×72 points (in streamwise, wall-normal and spanwise directions, respectively). Experimental results obtained by using particle image velocimetry (PIV) in a fully developed turbulent boundary layer in a flume, which were analyzed with the same statistical characterization method, are in agreement with the DNS analysis: the dominant vortical structure appears to have a longitudinal streamwise orientation, an inclination angle of about 8°, streamwise length of several hundred wall units, and a distance between the neighboring structures of about 100 wall units in the spanwise direction. PACS 47.27.Nz, 47.54+r     

Temporal development of turbulent boundary layers with embedded streamwise vortices

The interaction of streamwise vortices with turbulent boundary layer has been investigated using large-eddy simulation. The initial conditions are a pair of counterrotating Oseen vortices with flow between them directed toward the wall (common-flow-down), superimposed on various instantaneous realizations of a turbulent boundary layer. The time development of the vortices and their interaction with the boundary layer are studied by integrating the filtered Navier-Stokes equations in time. The most important effects of the vortices on the boundary layer are the thinning of the boundary layer between vortices (downwash region) and the thickening of the boundary layer in the upwash region. The vortices first move toward the wall as a result of the self-induced velocity, and then apart from each other because of the image vortices due to the solid wall. The Reynolds stress profiles highlight the highly three-dimensional structure of the turbulent boundary layer modified by the vortices. The presence of significant turbulent activity near the vortex center and in the upwash region suggests that localized instability mechanisms in addition to the convection of turbulent energy by the secondary flow are responsible for this effect. High levels of turbulent kinetic energy and secondary stresses in the vicinity of the vortex center are also observed. The numerical results show good agreement with experimental results.This work was supported by the Office of Naval Research under Grant N00014-89-J-1638. Computer time was supplied by the San Diego Supercomputing Center.     

基于涡旋识别方法的高速列车尾涡结构的讨论

利用改进型延迟分离涡模拟方法对缩尺比例1:30的高速列车简化模型的绕流流场进行数值计算,主要针对近尾流区的涡旋结构展开具体讨论. 通过不同的涡旋识别方法,发现在尾涡结构中,高涡量的强涡旋主要聚集于尾车附近,而涡量较低但处于相对稳定状态的涡旋分布在大部分尾流空间中. 对此,主要基于最新提出的涡旋定义及其物理意义认为,由于边界层在尾部发生的流动分离,剪切变形以及高涡量的扩散对强涡旋的形成发挥着重要的作用,而涡旋会被较强的剪切旋转拉伸,使得局部复杂的流动表现出突出的湍流特性;另一方面,尽管涡强度明显下降,但是在强剪切应变迅速衰减的情况下,流向涡核中的涡旋涡量是主要的,此时,在较接近地面的情况下,流体微团以涡核为中心的旋转运动使得涡旋与地面之间的相互作用成为主导的流动机制. 虽然涡强度会相对缓慢地衰减,但是从湍流能量产生的角度,该机制对涡旋的自维持发挥重要的作用,从而使尾涡结构能够相对稳定地存在于尾流流动中.      

Using digital holographic microscopy for simultaneous measurements of 3D near wall velocity and wall shear stress in a turbulent boundary layer

A digital holographic microscope is used to simultaneously measure the instantaneous 3D flow structure in the inner part of a turbulent boundary layer over a smooth wall, and the spatial distribution of wall shear stresses. The measurements are performed in a fully developed turbulent channel flow within square duct, at a moderately high Reynolds number. The sample volume size is 90 × 145 × 90 wall units, and the spatial resolution of the measurements is 3–8 wall units in streamwise and spanwise directions and one wall unit in the wall-normal direction. The paper describes the data acquisition and analysis procedures, including the particle tracking method and associated method for matching of particle pairs. The uncertainty in velocity is estimated to be better than 1 mm/s, less than 0.05% of the free stream velocity, by comparing the statistics of the normalized velocity divergence to divergence obtained by randomly adding an error of 1 mm/s to the data. Spatial distributions of wall shear stresses are approximated with the least square fit of velocity measurements in the viscous sublayer. Mean flow profiles and statistics of velocity fluctuations agree very well with expectations. Joint probability density distributions of instantaneous spanwise and streamwise wall shear stresses demonstrate the significance of near-wall coherent structures. The near wall 3D flow structures are classified into three groups, the first containing a pair of counter-rotating, quasi streamwise vortices and high streak-like shear stresses; the second group is characterized by multiple streamwise vortices and little variations in wall stress; and the third group has no buffer layer structures.     

On the etiology of shear layer vortices

In an effort to isolate the mechanism by which streamwise structures form in turbulent wall layers, evolution equations were derived for the streamwise velocity and vorticity perturbations about a mean turbulent fully developed channel flow. The stability of these equations, which take their most concise form when derived from the Generalized Lagrangian mean equations of Andrews and McIntyre, are studied assuming normal modes and infinitesimal disturbances. The resulting stability diagram yields, inter alia, the spanwise periodicity of the resulting structures, which we term shear layer vortices. If streaks are thought of as the footprints of these vortices, we then have a formal way of determining the spacing of streaks. The first three modes of instability are determined; at the first not just two vortices form per period, but four. It is also evident that an intense local shear layer forms about the plane in which the convection velocity equals the mean Eulerian velocity.Dedicated to Professor J.L. Lumley on the occasion of his 60th birthday.This work was supported in part by the U.S. Office of Naval Research under SRO IV Grant No. N00014-85-K-0172 and in part by the National Science Foundation Grant CTS-9008477.     

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.     

附面层中马蹄形旋涡的实验研究

本文用烟流显示和热线方法研究附面层中马蹄形旋涡的形态、生成频率、配对和涡系之间的相互作用。