On the correlation between enstrophy and energy dissipation rate in a turbulent wake

All three components of the vorticity fluctuation have been measured simultaneously in a turbulent wake using a new eight-sensor vorticity probe. The vorticity fluctuation spectra agree reasonably well with those from a direct numerical simulation of a turbulent channel flow at high wavenumbers. The similarity between the instantaneous energy dissipation rate ε and the instantaneous enstrophy ω² is examined using spectra and probability density functions. The correlation between ω² and ε is evaluated in some detail. The homogeneous value of ε is strongly correlated with ω². The full value of ε and, more especially its isotropic value, are less well correlated with the enstrophy. Conditional averaging indicates that high enstrophy regions are associated with high energy dissipation rate regions.     

The effect of fine structure on the stability of planar vortices

This study considers the linear, inviscid response to an external strain field of classes of planar vortices. The case of a Gaussian vortex has been considered elsewhere, and an enstrophy rebound phenomenon was noted: after the vortex is disturbed enstrophy feeds from the non-axisymmetric to mean flow. At the same time an irreversible spiral wind-up of vorticity fluctuations takes place. A top-hat or Rankine vortex, on the other hand, can support a non-decaying normal mode.In vortex dynamics processes such as stripping and collisions generate vortices with sharp edges and often with bands or rings of fine scale vorticity at their periphery, rather than smooth profiles. This paper considers the stability and response of a family of vortices that vary from a broad profile to a top-hat vortex. As the edge of the vortex becomes sharper, a quasi-mode emerges and vorticity winds up in a critical layer, at the radius where the angular velocity of the fluid matches that of a normal mode on a top-hat vortex. The decay rate of these quasi-modes is proportional to the vorticity gradient at the critical layer, in agreement with theory. As the vortex edge becomes sharper it is found that the rebound of enstrophy becomes stronger but slower.The stability and linear behaviour of coherent vortices is then studied for distributions which exhibit additional fine structure within the critical layer. In particular we consider vorticity profiles with ‘bumps’, ‘troughs’ or ‘steps’ as this fine structure. The modified evolution equation that governs the critical layer is studied using numerical simulations and asymptotic analysis. It is shown that depending on the form of the short-scale vorticity distribution, this can stabilise or destabilise quasi-modes, and it may also lead to oscillatory behaviour.     

Coherent vortical and straining structures in the finite wall-mounted square cylinder wake

Topological aspects of the turbulent wake of a finite, surface-mounted, square-cross-section cylinder of h/d = 4 are addressed by decomposing the velocity field into a quasi-periodic coherent part and the unresolved incoherent fluctuations. The three-dimensional large scale structure is educed through a reconstruction of planar phase-averaged PIV measurements using the simultaneously sampled surface pressure difference on opposing sides of the obstacle as a phase reference. A topological model for the vortex structure is educed and mean streamwise wake vorticity is explained in terms of the connections between initially vertical structures shed alternately from either side of the obstacle, rather than previously proposed ‘tip’ vortex structures generated at the obstacle free-end. The coherent structure educed accounts for a significant portion of the fluctuating energy in the wake. The turbulent field is further analyzed by finding Lagrangian straining structures that form by induction of the coherent vorticity field, and these structures are related to the energy transfer from the base phase-averaged flow since they act to stretch incoherent vorticity fluctuations in their neighbourhood.     

Education of coherent structures in a numerically simulated plane wake

Coherent structures (CS) are educed using a conditional sampling technique involving alignment of vorticity patches of largest size and strength; hence we educe dominant CS. A numerically simulated spatially evolving wake of a thick flat plate is used as the database, and the inflow condition for the simulated wake includes random velocity perturbations which emulate turbulent conditions at a plate exit in the laboratory. In addition to previously educed properties such as coherent vorticity and production, and incoherent Reynolds stress and turbulence intensity, other measures such as coherent pressure and passive scalar distributions are also studied. In spite of the geometry difference, the near-wake dynamics of the plate seem quite similar to that of a cylinder. For example, turbulence is mostly produced by vortex stretching of the ribs at the saddle and then advected to the structure center, where it accumulates, and is balanced by incoherent dissipation. The distribution of coherent passive scalar indicates that mixing occurs in the saddle regions and that the mixed fluid is advected into the structure center.     

Three-component velocity field measurements of propeller wake using a stereoscopic PIV technique

A stereoscopic PIV (Particle Image Velocimetry) technique was used to measure the three-dimensional flow structure of the turbulent wake behind a marine propeller with five blades. The out-of-plane velocity component was determined using two CCD cameras with an angular displacement configuration. Four hundred instantaneous velocity fields were measured for each of four different blade phases, and ensemble averaged in order to find the spatial evolution of the propeller wake in the region from the trailing edge up to one propeller diameter (D) downstream. The influence of propeller loading conditions on the wake structure was also investigated by measuring the velocity fields at three advance ratios (J=0.59, 0.72 and 0.88). The phase-averaged velocity fields revealed that a viscous wake formed by the boundary layers developed along the blade surfaces. Tip vortices were generated periodically and the slipstream contracted in the near-wake region. The out-of-plane velocity component and strain rate had large values at the locations of the tip and trailing vortices. As the flow moved downstream, the turbulence intensity, the strength of the tip vortices, and the magnitude of the out-of-plane velocity component at trailing vortices all decreased due to effects such as viscous dissipation, turbulence diffusion, and blade-to-blade interaction.     

Planar velocity measurements of a two-dimensional compressible wake

The present study describes the application of particle image velocimetry (PIV) to investigate the compressible flow in the wake of a two-dimensional blunt base at a freestream Mach number M_X=2. The first part of the study addresses specific issues related to the application of PIV to supersonic wind tunnel flows, such as the seeding particle flow-tracing fidelity and the measurement spatial resolution. The seeding particle response is assessed through a planar oblique shock wave experiment. The measurement spatial resolution is enhanced by means of an advanced image-interrogation algorithm. In the second part, the experimental results are presented. The PIV measurements yield the spatial distribution of mean velocity and turbulence. The mean velocity distribution clearly reveals the main flow features such as expansion fans, separated shear layers, flow recirculation, reattachment, recompression and wake development. The turbulence distribution shows the growth of turbulent fluctuations in the separated shear layers up to the reattachment location. Increased velocity fluctuations are also present downstream of reattachment outside of the wake due to unsteady flow reattachment and recompression. The instantaneous velocity field is analyzed seeking coherent flow structures in the redeveloping wake. The instantaneous planar velocity and vorticity measurements return evidence of large-scale turbulent structures detected as spatially coherent vorticity fluctuations. The velocity pattern consistently shows large masses of fluid in vortical motion. The overall instantaneous wake flow is organized as a double row of counter-rotating structures. The single structures show vorticity contours of roughly elliptical shape in agreement with previous studies based on spatial correlation of planar light scattering. Peak vorticity is found to be five times higher than the mean vorticity value, suggesting that wake turbulence is dominated by the activity of large-scale structures. The unsteady behavior of the reattachment phenomenon is studied. Based on the instantaneous flow topology, the reattachment is observed to fluctuate mostly in the streamwise direction suggesting that the unsteady separation is dominated by a pumping-like motion.     

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

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

Dynamics of hairpin vortices generated by a mixing tab in a channel flow

To better understand mixing by hairpin vortices, time-series particle image velocimetry (PIV) was applied to the wake of a trapezoidal-shaped passive mixing tab mounted at the bottom of a square turbulent channel (Re _h =2,080 based on the tab height). Instantaneous velocity/vorticity fields were obtained in sequences of 10 Hz in the tab wake in the center plane (x–y) and in a plane (x–z) parallel to the wall. Periodically-shed hairpin vortices were clearly identified and seen to rise as they advected downstream. Experimental evidence shows that the vortex-induced ejection of the near-wall viscous fluid to the immediate upstream is important to the dynamics of hairpin vortices. It can increase the strength of the hairpin vortices in the near tab region and cause generation of secondary hairpin vortices further downstream when the hairpin heads are farther away from the wall. Measurements also reveal the existence of a type of new secondary vortice with the opposite-sign spanwise vorticity. The distribution of vortex loci in the x–y plane shows that the hairpin vortices and the reverse vortices are spatially segregated in distinct layers. Turbulence statistics, including mean velocity profiles, Reynolds stresses, and turbulent kinetic energy dissipation rate distributions, were obtained from the PIV data. These statistical quantities clearly reveal imprints of the identified vortex structures and provide insight into mixing effectiveness. Received: 24 February 2000/Accepted: 24 October 2000     

Vortex dynamics of a sphere wake in proximity to a wall

Toward getting the vortex dynamics characteristics and wake structure of a sphere in proximity to a wall, the effect of a proximal flat plate on the wake of a stationary sphere is investigated via direct numerical simulation. The vortex shedding process and the significant variation of the wake structure are described in detail. The drag coefficient reduces and the wake structure of the sphere becomes complex due to the combined effect of the wake flow and the wall. A jet flow forms between the sphere and the flat plate, which suppresses the vortex separation on the bottom of the sphere. The asymmetric distributions of the coherent structures and the recirculation region behind the sphere are discussed. Besides vortex shedding patterns, the time-averaged velocity distribution, vortex dynamics, distribution regularities of turbulent kinetic energy and enstrophy are investigated.     

PIV analysis of flow around a container ship model with a rotating propeller

The flow characteristics of the propeller wake behind a container ship model with a rotating propeller were investigated using a two-frame PIV (Particle Image Velocimetry) technique. Ensemble-averaged mean velocity fields were measured at four different blade phases and ensemble-averaged to investigate the flow structure in the near-wake region. The mean velocity fields in longitudinal planes show that a velocity deficit is formed in the regions near the blade tips and hub. As the flow develops in the downstream direction, the trailing vortices formed behind the propeller hub move upward slightly due to the presence of the hull wake and free surface. Interaction between the bilge vortices and the incoming flow around the hull causes the flow structure to be asymmetric. Contour plots of the vorticity give information on the radial distribution of the loading on the blades. The radial velocity profiles fluctuate to a greater extent under the heavy (J=0.59) and light loading (J=0.88) conditions than under the design loading condition (J=0.72). The turbulence intensity has large values around the tip and trailing vortices. As the wake develops in the downstream direction, the strength of the vorticity diminishes and the turbulence intensity increases due to turbulent diffusion and active mixing between the tip vortices and the adjacent wake flow.     

Coherent Vortex Simulation (CVS), A Semi-Deterministic Turbulence Model Using Wavelets

In the spirit of Ha Minh's semi-deterministic model, we propose a new method for computing fully-developed turbulent flows, called Coherent Vortex Simulation (CVS). It is based on the observation that turbulent flows contain both an organized part, the coherent vortices, and a random part, the incoherent background flow. The separation into coherent and incoherent contributions is done using the wavelet coefficients of the vorticity field and the Biot–Savart kernel to reconstruct the coherent and incoherent velocity fields. The evolution of the coherent part is computed using a wavelet basis, adapted at each time step to resolve the regions of strong gradients, while the incoherent part is discarded during the flow evolution, which models turbulent dissipation. The CVS method is similar to LES, but it uses nonlinear multiscale band-pass filters, which depend on the instantaneous flow realization, while LES uses linear low-pass filters, which do not adapt to the flow evolution. As example, we apply the CVS method to compute a time developing two-dimensional mixing layer and a wavelet forced two-dimensional homogeneous isotropic flow. We also demonstrate how walls or obstacles can be taken into account using penalization and compute a two-dimensional flow past an array of cylinders. Finally, we perform the same segmentation into coherent and incoherent components in a three-dimensional homogeneous isotropic turbulent flow. We show that the coherent components correspond to vortex tubes, which exhibit non-Gaussian statistics and long-range correlation, with the same k ^−5/3power-law energy spectrum as the total flow. In contrast, the incoherent components correspond to an homogeneous random background flow which does not contain organized structures and presents an energy equipartition together with a Gaussian PDF of velocity. This justifies their elimination during the CVS computation to model turbulent dissipation. This revised version was published online in August 2006 with corrections to the Cover Date.     

Vortices,Complex Flows and Inertial Particles

The properties of vortical structures at high Reynolds number in uniform flows and near rigid boundaries are reviewed. New properties are derived by analysing the dynamics of the main flow features and the related integral constraints, including the relations between mean swirl and bulk speed, the relative level of internal fluctuations to bulk properties, and connections between the steadiness and topology of the structures. A crucial property that determines energy dissipation and the transport of inertial particles (with finite fall speed) is the variation across the structure of the ratio of the mean strain rate (Σ) to the mean vorticity (Ω). It is shown how, once such particles are entrained into the vortical regions of a coherent structure, they can be transported over significant distances even as the vortices grow and their internal structure is distorted by internal turbulence, swirling motions and the presence of rigid boundaries. However if the vortex is strongly distorted by a straining motion so that Σ is greater than Ω, the entrained particles are ejected quite rapidly. These mechanisms are consistent with previous studies of entrained and sedimenting particles in disperse two phase flows over flat surfaces, and over bluff obstacles and dunes. They are also tested in more detail here through laboratory observations and measurements of 50–200-μm particles entrained into circular and non-circular vortices moving first into still air and then onto rigid surfaces placed parallel and perpendicular to the direction of motion of the vortices.     

Vortex motion in a rotating barotropic fluid layer

To study vortex motion and the mechanisms of geostrophic adjustment (i.e. the equilibrium between pressure gradient and Coriolis force, which leads to the weakening of inertio-gravity waves) in large scale geophysical flows, we simulate the dynamics of a shallow-water layer in uniform rotation, without any forcing other than the initial injection of energy and potential enstrophy. Such a flow generates inertio-gravity waves which interact with the rotational eddies. We found that both inertio-gravity waves and rotation reduce the non-linear interactions between vortices, namely the condensation of the vorticity field into isolated coherent vortices, corresponding to the inverse rotational energy cascade, and the associated production of vorticity filaments, due to the direct potential enstrophy cascade. Rotation also inhibits the direct inertio-gravitational energy cascade for scales larger than the Rossby deformation radius. Therefore, if inertio-gravity waves are initially excited at large enough scales, they will remain trapped there due to rotation and there will be no geostrophic adjustment. On the contrary, if inertio-gravity waves are only present at scales smaller than the Rossby deformation radius, which are insensitive to the effect of rotation, they will non-linearly interact and cascade towards the dissipative scales, leaving the flow in geostrophic equilibrium.     

Partial roll-up of a viscoelastic Kármán street

The instantaneous velocity field of a circular cylinder wake is built using a PIV technique when a small amount of viscoelastic liquid is introduced through the cylinder. It is shown that the viscoelastic fluid slows down the vorticity concentration and produces a street of partially rolled-up vortices. The underlying mechanism appears very analogous to that of a surface tension in the Kelvin–Helmholtz instability. The partial roll-up is studied in terms of the Weiss determinant. This quantity is a local measurement of the spatial separation between the strain and vorticity. In the viscoelastic wake, the Weiss determinant reaches much lower values than in the Newtonian wake. This result shows that the elasticity prevents the clear separation between vorticity and strain during the roll-up process. Since the Weiss determinant is directly related to the pressure field, it suggests that elasticity can drastically modify the pressure levels even when vorticity and strain levels are unaffected.     

Comparison of flow structures in the downstream region of a cylinder and sphere

An experimental investigation of flow structures downstream of a circular cylinder and sphere immersed in a free-stream flow is performed for Re = 5000 and 10,000 using qualitative and quantitative flow visualization techniques. The obtained results are presented in terms of time-averaged velocity vectors, patterns of streamlines, vorticity, Reynolds stress correlations and turbulent kinetic energy distributions. Flow data reveal that the size of wake flow region, the location of singular and double points, the peak values of turbulence quantities, such as Reynolds stress correlations, vorticity fluctuations and turbulent kinetic energy vary as a function of models’ geometry and Reynolds Numbers. The concentration of small scale vortices is more dominant in the wake of the sphere than that of the cylinder. The maximum value of turbulent kinetic energy (TKE) occurs close to the saddle point for the cylinder case while two maximum values of TKE occur along shear layers for the sphere one because of the 3-D flow behavior.     

Comparison of flow structures in the downstream region of a cylinder and sphere

An experimental investigation of flow structures downstream of a circular cylinder and sphere immersed in a free-stream flow is performed for Re = 5000 and 10,000 using qualitative and quantitative flow visualization techniques. The obtained results are presented in terms of time-averaged velocity vectors, patterns of streamlines, vorticity, Reynolds stress correlations and turbulent kinetic energy distributions. Flow data reveal that the size of wake flow region, the location of singular and double points, the peak values of turbulence quantities, such as Reynolds stress correlations, vorticity fluctuations and turbulent kinetic energy vary as a function of models’ geometry and Reynolds Numbers. The concentration of small scale vortices is more dominant in the wake of the sphere than that of the cylinder. The maximum value of turbulent kinetic energy (TKE) occurs close to the saddle point for the cylinder case while two maximum values of TKE occur along shear layers for the sphere one because of the 3-D flow behavior.     

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.     

Asymmetric vortex shedding flow past an inclined flat plate at high incidence

This paper reports an experimental investigation of the vortex shedding wake behind a long flat plate inclined at a small angle of attack to a main flow stream. Detailed velocity fields are obtained with particle-image velocimetry (PIV) at successive phases in a vortex shedding cycle at three angles of attack, α=20°, 25° and 30°, at a Reynolds number Re≈5,300. Coherent patterns and dynamics of the vortices in the wake are revealed by the phase-averaged PIV vectors and derived turbulent properties. A vortex street pattern comprising a train of leading edge vortices alternating with a train of trailing edge vortices is found in the wake. The trailing edge vortex is shed directly from the sharp trailing edge while there are evidences that the formation and shedding of the leading edge vortex involve a more complicated mechanism. The leading edge vortex seems to be shed into the wake from an axial location near the trailing edge. After shedding, the vortices are convected downstream in the wake with a convection speed roughly equal to 0.8 the free-stream velocity. On reaching the same axial location, the trailing edge vortex, as compared to the leading edge vortex, is found to possess a higher peak vorticity level at its centre and induce more intense fluid circulation and Reynolds stresses production around it. It is found that the results at the three angles of attack can be collapsed into similar trends by using the projected plate width as the characteristic length of the flow.     

Large Eddy Simulation of boundary layer transition over an isolated ramp-type micro roughness element

Boundary layer transition over an isolated surface roughness element is investigated by means of numerical simulation. Large Eddy Simulation (LES) flow-modeling approach is employed to study flow characteristics and transition phenomenon past a roughness element immersed within an incoming developing boundary layer, at a height-based Reynolds number of 1170. LES numerical results are compared to experimental data from literature showing the time-averaged velocity distribution, the velocity fluctuation statistics and the instantaneous flow topology.Despite slight difference in the intensity of streamwise velocity fluctuations, the present LES results and experimental data show very good agreement. The mean flow visualization shows streamwise counter-rotating vortices pairs formation downstream of the obstacle. The primary pair induces an upwash motion and a momentum deficit that creates a Kelvin-Helmholtz type flow instability. The instantaneous flow topology reveals the formation of coherent K-H vortices downstream that produce turbulent fluctuations in the wake of the roughness element. These vortices are streched and lifted up when moving downstream. The velocity fluctuations results show that the onset of the turbulence is dominated by the energy transfer of large-scale vortices.     

An experimental study on unsteady turbulent near wake of a rectangular cylinder in channel flow

 Unsteady turbulent near wake of a rectangular cylinder in channel flow has been studied experimentally with a laser Doppler velocimetry (LDV). The time-averaged and phase-averaged statistics were measured for the cylinders having various width-to-height ratios, b/h. It is shown that the turbulent intensities on the centerline of the channel have their maxima near the rear stagnation point of a recirculation region. The contours of coherent vorticity and streamline reproduce clearly the shed vortices from the cylinder observed by the flow visualization. The characteristics of the flow field, which depends on b/h, are discussed and the significant contribution of the coherent structure to the flow field is clarified. Moreover, the turbulent kinetic energy budget has been examined. Received: 19 January 1998/Accepted: 21 July 1998