The effects of acceleration in jets: kinematics of the near field vortices

Direct and large-eddy simulations (DNS/LES) of accelerating round jets are used to analyze the effects of acceleration on the kinematics of vortex rings in the near field of the jet (x/D < 12). The acceleration is obtained by increasing the nozzle jet velocity with time, in a previously established (steady) jet, and ends once the inlet jet velocity is equal to twice its initial value. Several acceleration rates (α = 0.02–0.6) and Reynolds numbers (Re _D = 500–20000) were simulated. Acceleration maps were used to make a detailed study of the kinematics of vortex rings in accelerating jets. One of the effects of the acceleration is to cause a number of new primary and secondary vortex merging events that are absent from steady jets. As the acceleration rate α increases, both the number of primary merging events between rings and the axial position where these take place decreases. The statistics for the speed of the starting ring that forms at the start of the acceleration phase for each simulation, agree well with the statistics for the “front” speed observed by Zhang and Johari (Phys Fluids 8:2185–2195, 1996). Acceleration maps and flow visualizations show that during the acceleration phase the near field coherent vortices become smaller and are formed at an higher frequency than in the steady jet, and their (mean) shedding frequency increases linearly with the acceleration rate. Finally, it was observed that the acceleration decreases the spreading rate of the jet, in agreement with previous experimental works.      

Flow around a circular cylinder—structure of the near wake shear layer

The separated shear layer in the near wake of a circular cylinder was investigated using a single hot wire probe, with special attention given to the shear layer instability characteristics. Without end plates to force parallel vortex shedding, the critical Reynolds number for the onset of the instability was 740. The present data, together with all previously published data, show that the ratio of the instability frequency f_sl to the vortex shedding frequency f_v varies as Re^0.65, which is in agreement with the Re^0.67 dependence obtained by Prasad and Williamson [1997, J Fluid Mech 333:375–402]. However, the distribution of f_sl/f_v and the spectra of the longitudinal velocity fluctuation (u) suggest that, on either side of Re=5,000, the shear layer exhibits lower and upper subcritical regimes, in support of the observations by Norberg [1987, publication no. 87/2, Chalmers University of Technology, Sweden] and Prasad and Williamson [1997, J Fluid Mech 343:235–265]. The spectra of u provide strong evidence for the occurrence of vortex pairing in wake shear layers, suggesting that the near wake develops in a similar manner to a mixing layer.     

Vortex shedding in cylinder flow of shear-thinning fluids: I. Identification and demarcation of flow regimes

An experimental study on the flow of non-Newtonian fluids around a cylinder was undertaken to identify and delimit the various shedding flow regimes as a function of adequate non-dimensional numbers. The measurements of vortex shedding frequency and formation length (l_f) were carried out by laser-Doppler anemometry in Newtonian fluids and in aqueous polymer solutions of CMC and tylose. These were shear thinning and elastic at weight concentrations ranging from 0.1 to 0.6%. The 10 and 20 mm diameter cylinders (D) used in the experiments had aspect ratios of 12 and 6 and blockage ratios of 5 and 10%, respectively. The Reynolds number (Re^*) was based on a characteristic shear rate of U_∞/(2D) and ranged from 50 to 9×10³ thus encompassing the laminar shedding, the transition and shear-layer transition regimes. Increasing fluid elasticity reduced the various critical Reynolds numbers (Re_etr^*, Re_lf^*, Re_bbp^*) and narrowed the extent of the transition regime. For the 0.6% tylose solution the transition regime was even suppressed. On the other end, pseudoplasticity was found to be indirectly responsible for the observed reduction in Re_otr^*: it increases the Strouhal number which in turn increases the vortex filaments, precursors of the transition regime. Elasticity was better quantified by the elasticity number Re′/We than by the Weissenberg number. This elasticity number involves the calculation of the viscosity at a high characteristic shear rate, typical of the boundary layer, rather than at the average value (U_∞/(2D)) used for the Reynolds number, Re^*.     

Non-linear dynamics of semi-dilute polydisperse polymer solutions in microfluidics: effects of flow geometry

The non-linear dynamics of a semi-dilute (c/c* = 15) polydisperse polyethylene oxide (PEO) solution in microfluidics are studied experimentally using benchmark contraction–expansion flow geometries with three contraction–expansion ratios (4:1:4, 8:1:8 and 16:1:16) and two narrow channel lengths (L _c/D _h = 53 and 5.3, where L _c is the length of the narrow channel and D _h is its hydraulic diameter). Complex flows over a range of elasticity numbers (El), Weissenberg numbers (Wi) and Reynolds numbers (Re) are characterized using micro-particle image velocimetry ( \upmu\upmu-PIV) and pressure drop measurements. The evolution of vortex formation and dynamics has been visualized through a step-flow-rate experiment. Various flow dynamics regimes have been quantified and are presented in a Wi–Re diagram. The experimental results reveal that the contraction ratio can result in qualitatively different vortex dynamics of semi-dilute polymer solutions in microfluidics, whereas the length of the narrow channel merely affects the dynamics at a quantitative level. A single elasticity number, if defined by the size of the narrow channel, is not sufficient to account for the effects of contraction ratio on the non-linear vortex dynamics.     

The effect of chamber length and Reynolds number on jet precession

The jet axial velocity field exiting from a nozzle/chamber configuration with an expansion ratio of 5 is investigated using Stereo-PIV for a range of chamber lengths and Reynolds (Re) numbers. The jet can exit the chamber in axial jet (AJ) mode with the maximum velocity near the chamber axis or precessing jet (PJ) mode with the maximum velocity near the chamber wall and rotating or precessing about the chamber axis. Algorithms were developed to determine the jet mode from exit conditions and allow conditional averaging of the velocity field in PJ mode. The probability of the jet in PJ mode was found to be a strong function of chamber length, L/D and only a mild function of Re for Re > 10,000. High precession probability was found for chambers of length in the range 2 < L/D < 2.75 for all cases for Re > 10,000. An abrupt reduction in precession probability occurred for chamber lengths L/D~3. For increasing chamber lengths, an increase in precession probability was observed. The ratio of entrainment-into-the-chamber of surrounding fluid to jet exit fluid was found not to be a function of Re or jet mode (AJ or PJ) but only a function of L/D. A maximum ratio entrainment-into-the-chamber was observed to occur in the range 2 < L/D < 2.5. Conditionally averaged velocity profiles also showed the exiting jet to be a strong function of L/D and with only a mild effect of Re for all cases of Re > 10,000.     

Vortex characteristics in laminar cavity flow at very low Mach number

In the present paper, a laminar cavity is analysed at very low Mach numbers. The characteristics of core-vortices are proposed and commented. The experiments were performed in an open subsonic wind tunnel using particle image velocimetry (PIV). A rectangular cavity with a length-to-depth ratio of 4 was used (shallow and open type). Three different Reynolds numbers, based on cavity depth and free stream velocity, were examined (Re_h=4,000, 9,000 and 13,000). The upstream boundary layer was investigated using classical hot-wire anemometry and was found to be laminar. For each Reynolds number, a total of 1,000 vectors fields were acquired. The results are given in terms of conventional quantities (mean flow velocity, turbulence characteristics, Reynolds shear stress) and also in terms of vortex characteristics (such as probability density function of vortex location, vortex size and vortex circulation). Some of these vortex characteristics are then proposed in a local averaged presentation. The extraction of vortices from instantaneous flow fields has been done through the use of a home-made algorithm based on continuous wavelet analysis.     

Effect of local forcing on a turbulent boundary layer

An experimental study is performed to analyze flow structures behind local suction and blowing in a flat-plate turbulent boundary layer. The local forcing is given to the boundary layer flow by means of a sinusoidally oscillating jet issuing from a thin spanwise slot at the wall. The Reynolds number based on the momentum thickness is about Re _θ =1700. The effects of local forcing are scrutinized by altering the forcing frequency (0.011 ≤ f ⁺≤ 0.044). The forcing amplitude is fixed at A ₀=0.4. It is found that a small local forcing reduces the skin friction and the skin friction reduction increases with the forcing frequency. A phase-averaging technique is employed to capture the large-scale vortex evolution. An organized spanwise vortical structure is generated by the local forcing. The cross-sectional area of vortex and the time fraction of vortex are examined by changing the forcing frequency. An investigation of the random fluctuation components reveals that turbulent energy is concentrated near the center of vortical structures. Received: 17 March 2000/Accepted: 3 April 2001     

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     

An experimental investigation of flow past a dual step cylinder

Flow development in the wake of a dual step cylinder has been investigated experimentally using Laser Doppler Velocimetry and flow visualization. The dual step cylinder model is comprised of a large diameter cylinder (D) mounted at the mid-span of a small diameter cylinder (d). The experiments have been performed for a Reynolds number (Re _D ) of 1,050, a diameter ratio (D/d) of 2, and a range of large cylinder aspect ratios (L/D). The results show that the flow development is highly dependent on L/D. The following four distinct flow regimes can be identified based on vortex dynamics in the wake of the large cylinder: (1) for L/D ≥ 15, three vortex shedding cells form in the wake of the large cylinder, one central cell bounded by two cells of lower frequency, (2) for 8 < L/D ≤ 14, a single vortex shedding cell forms in the wake of the large cylinder, (3) for 2 < L/D ≤ 6, vortex shedding from the large cylinder is highly three-dimensional. When spanwise vortices are shed, they deform substantially and attain a hairpin shape in the near wake, (4) for 0.2 ≤ L/D ≤ 1, the large cylinder induces vortex dislocations between small cylinder vortices. The results show that for Regimes I to III, on the average, the frequency of vortex shedding in the large cylinder wake decreases with L/D, which is accompanied by a decrease in coherence of the shed vortices. In Regime IV, small cylinder vortices connect across the large cylinder wake, but these connections are interrupted by vortex dislocations. With decreasing L/D, the frequency of dislocations decreases and the dominant frequency in the large cylinder wake increases toward the small cylinder shedding frequency.     

Performance of a 12-sensor vorticity probe in the near field of a rectangular turbulent jet

The design and operational characteristics of a 12-sensor hot wire probe for three-dimensional velocity–vorticity measurements in turbulent flow fields is described and discussed. The performance of the probe is investigated in comparison with X-sensor probe measurements in the near field of a rectangular turbulent jet with aspect ratio 6. Measurements have been conducted at Reynolds number Re _D = 21,000 at nozzle distances of x/D = 1, 3, 6 and 11, where D is the width of the nozzle. The results obtained with the 12-sensor probe compare well to the results of the X-sensor probe. Distributions of mean and fluctuating velocity–vorticity fields are presented and discussed. Among the results the most prominent is the experimental confirmation of the high levels of fluctuating vorticity in the shear layers.     

Mechanical forcing of the wake of a flat plate

We report experimental results of the forced wake of a thin symmetric flat plate, placed parallel to an uniform air stream, in the range of thickness-based Reynolds number 50< Re _e <200. External wake forcing was introduced by small harmonic oscillations of a moving flap, placed at the trailing-edge of the flat plate. When the flap remains in a fixed horizontal position, the mean velocity profiles obtained by hot wire measurements, for different Reynolds numbers, are self similar. In the presence of harmonic forcing, within a certain range of the forcing frequency, the mean velocity profiles change and coherent structures are formed in the wake. Two independent flow-type resonances were observed: (i) when the inverse of the forcing frequency matches the flight time of the fluid particles along the flap. (ii) when the forcing frequency of the flap equals one half of the vortex shedding frequency of the flat plate and flap system. Implications of the two observed resonances on the wake structure are important. The first resonance (i) is associated to a wide but less intense (energy fluctuations) wake flow and the second resonance (ii) generates a thin but intense resultant wake flow.     

Measurement and analysis of unsteady flow structures in rotor blade wakes

Time-variant data are obtained to investigate the exit flow field from a rotor in a research compressor. In the free-stream region, the instantaneous data are analogous to one another and to the ensemble averaged free-stream results. However, in the wake region, some of the instantaneous signals are similar to one another and to the ensemble averaged wake, but others differ significantly. These variations in the instantaneous data are interpreted and shown to be due to a vortex street structure in the wake. This is accomplished by: (1) developing a mathematical model of the rotor blade exit flow field based on a wake vortex street structure analogous to the unsteady flow field behind bluff bodies due to classical von Karman vortex shedding; and (2) correlating predictions of both the ensemble averaged and instantaneous rotor blade exit flow fields as well as the velocity probability density distributions from this vortex wake flow field model with the corresponding data. The correlation of the ensemble averaged rotor blade exit flow fields is very good and the flow angle distribution correlation excellent. The predicted instantaneous rotor blade exit flow field exhibits many of the flow features found in the data. Also, the probability density distributions for the data and the vortex wake flow field model are analogous to one another.List of symbols N number of rotor revolutions - S _w rotor blade wake width - S _x vortex core horizontal spacing - S _y vortex core vertical spacing - u velocity component parallel to vortex street motion - v velocity component normal to vortex street motion - W instantaneous relative velocity - W _i velocity induced by vortex street - W free-stream relative velocity - W _s velocity of vortex street - x coordinate parallel to vortex street motion - y coordinate normal to vortex street motion - free-stream relative flow angle - _inst instantaneous relative flow angle - vortex strength     

Effects of riblet bends on pipe flows

Four riblet bends were tested to investigate the effects of riblets on pipe flows including the secondary flow on the Reynolds numbers; Re _D =6×10³–4×10⁴. The pressure gradients on the smooth pipe downstream from the riblet bends were measured, and also the pressure losses of the bends only were measured. All riblet bends reduced the pressure gradient on the smooth pipe downstream from them, which means a drag reduction. Two of the riblet bends showed the maximum drag reduction of about 4 percent at Re _D = 6500; this reduction rate was significant considering the uncertainty of the present experiments. Since the pressure losses of these two riblet bends were almost identical to that of the smooth bend at Re _D = 6500, they could cause a net drag reduction of about 4 percent on the piping system including these bends at that Reynolds number. Furthermore, the velocity profiles measured by LDV indicated that the secondary flow becomes weaker downstream from the riblet bends when a drag reduction is recognized there.Nomenclature D pipe diameter - D ₀ the distance from the valley to the valley passing through the pipe center - H height of groove - P nondimensional static pressure (p/it/(U ₀ ² ):p is gauge pressure) - dP/dX nondimensional pressure gradient - Rc curvature of bend - Re _D Reynolds number based on bulk velocity and pipe diameter - s spacing of groove - U mean streamwise velocity along the horizontal diameter - U ₀ bulk velocity - V mean vertical velocity along the horizontal diameter - x streamwise direction along the pipe axis (see Fig. 1) - X nondimensionalx (=x/D) - y radial direction in the horizontal plane which is perpendicular to the plane including the bend (see Fig. 1) - yUV swirl intensity (nondimensional swirl intensity:yUV/(DU ₀ ² ))     

Turbulent wake and vortex shedding for a stack partially immersed in a turbulent boundary layer

The effect of the jet-to-cross-flow velocity ratio, R, on the turbulent wake and Kármán vortex shedding for a cylindrical stack of aspect ratio AR=9 was investigated in a low-speed wind tunnel using thermal anemometry. The cross-flow Reynolds number was Re_D=2.3×10⁴, the jet Reynolds number ranged from Re_d=7.6×10³ to 4.7×10⁴, and R was varied from 0 to 3. The stack was partially immersed in a flat-plate turbulent boundary layer, with a boundary layer thickness-to-stack-height ratio of δ/H=0.5 at the location of the stack. From the behaviour of the turbulent wake and the vortex shedding, the flow around the stack could be classified into three regimes depending on the value of R, which were the downwash (R<0.7), cross-wind-dominated (0.7R<1.5), and jet-dominated (R1.5) flow regimes. Each flow regime had a distinct structure to the mean velocity (streamwise and wall-normal directions), turbulence intensity (streamwise and wall-normal directions), and Reynolds shear stress fields, as well as the variation of the Strouhal number and the power spectrum along the stack height.     

Measuring velocity fields in wind tunnels with holographic interferometry

This paper investigates the feasibility of using holographic interferometry in wind tunnel flows for measuring velocity fields rather than density or temperature fields. First results were obtained in a vortex street behind a cylinder at Re=190(U _∞=0.7 m/s). The light scattered from an illuminated fluid plane was holographically recorded twice with the same reference beam. Using a time interval of 10 μs, local fluid displacements smaller than a few microns were recorded. The holographic plate was placed in front and as close as possible to the fluid plane. The interferograms obtained from the hologram reconstruction give information about one velocity component, at 45° with the illuminated plane. The alignment of the cylinder axis with this 45° direction provided definite confirmation about the vortex street having a non-negligible axial velocity. The constant velocity fluid region has proven to be very useful for quantifying the velocity information contained in the interferogram. Received: 8 November 1999/Accepted: 14 March 2000     

Passive control of wake flow by two small control cylinders at Reynolds number 80

Passive control of the wake behind a circular cylinder in uniform flow is studied by numerical simulation at Re_D=80. Two small control cylinders are placed symmetrically along the separating shear layers at various stream locations. In the present study, the detailed flow mechanisms that lead to a significant reduction in the fluctuating lift but maintain the shedding vortex street are clearly revealed. When the stream locations lie within 0.8≤X_C/D≤3.0, the alternate shedding vortex street remains behind the control cylinders. In this case, the symmetric standing eddies immediately behind the main cylinder and the downstream delay of the shedding vortex street are the two primary mechanisms that lead to a 70–80% reduction of the fluctuating lift on the main cylinder. Furthermore, the total drag of all the cylinders still has a maximum 5% reduction. This benefit is primarily attributed to the significant reduction of the pressure drag on the main cylinder. Within X_C/D>3.0, the symmetry of the standing eddy breaks down and the staggered vortex street is similar to that behind a single cylinder at the same Reynolds number. In the latter case, the mean pressure drag and the fluctuating lift coefficients on the main cylinder will recover to the values of a single cylinder.     

Vortex mechanism of heat transfer enhancement in a channel with spherical and oval dimples

Vortex mechanism of heat transfer enhancement in a narrow channel with dimples has been investigated numerically using LES and URANS methods. The flow separation results in a formation of vortex structures which significantly enhance heat transfer on dimpled surfaces leading to a small increase in pressure loss. The heat transfer can be significantly increased by rounding the dimple edge and use of oval dimples. To get a deep insight into flow physics LES is performed for single phase flow in a channel with a spherical dimple. The instantaneous vortex formation and separation are investigated in and around the dimple area. Considered are Reynolds numbers (based on dimple print diameter) Re_D = 20,000 and Re_D = 40,000 the depth to print diameter ratio of Δ = 0.26. Frequency analysis of LES data revealed the presence of dominating frequencies in unsteady flow oscillations. Direct analysis of the flow field revealed the presence of coherent vortex structure inclined to the mean flow. The structure changes its orientation in time causing the long period oscillations with opposite-of-phase motion. Three dimensional proper orthogonal decomposition (POD) analysis is carried out on LES pressure and velocity fields to identify spatio-temporal structures hidden in the random fluctuations. Tornado-like spatial POD structures have been determined inside dimples.     

Smoke-wire flow visualization of the near-wake region behind a circular disc at low Reynolds numbers

The flow structures in the near field of the unducted wake region behind a circular disc for annular flow at low Reynolds numbers were studied by smoke-wire flow visualization technique. A twisted-dual-wire was employed to perform the time evolving visualization. Three typical characteristic flow modes: Q-tip, open-top toroid, and closed toroid, were identified in the near disc region. For Reynolds number between 130 and 390, the Q-tip flow mode which subject to a periodic up-down oscillatory motion was observed. The open-top toroid mode which experiences the expelling vortex shedding was found for Reynolds number between 390 and 455. The free separation surface turns around and merges to the central axisymmetric axis to form the conventionally observed toroidal recirculation bubble for Reynolds number higher than 455. The closed toroid mode exhibits both expelling and shear-layer vortex sheddings. With the identified flow modes at low Reynolds numbers, the recirculation contours, recirculation length, and the shedding frequency in each mode were measured and discussed.List of symbols B.R. blockage ratio (=D ² /D _a ² ) - D _a outer diameter of annular jet, 30 mm - D diameter of circular disc, 20 mm - f frequency of vortex shedding, Hz - L _r axial length of recirculation zone - R radius of circular disc, 10 mm - u _a average exit velocity of annular jet - ₀ stream function with value of zero - mass density of annular flow - u average axial velocity - r radial coordinate, originated from center of circular disk - r ₀ radial coordinate of the boundary of the recirculation zone - Re _a Reynolds number of annular jet based on the disc diameter - Z axial coordinate, originated from center of circular disk - w _max maximum half-width of the recirculation zone - St Strouhal number (=fD/D _a )     

Scale-adaptive simulation of flow past wavy cylinders at a subcritical Reynolds number

The Xu & Yan scale-adaptive simulation (XYSAS) model is employed to simulate the flows past wavy cylinders at Reynolds number 8 × 10 3.This approach yields results in good agreement with experimental measurements.The mean flow field and near wake vortex structure are replicated and compared with that of a corresponding circular cylinder.The effects of wavelength ratios λ/D m from 3 to 7,together with the amplitude ratios a /D m of 0.091 and 0.25,are fully investigated.Owing to the wavy configuration,a maximum reduction of Strouhal number and root-meansquare (r.m.s) fluctuating lift coefficients are up to 50% and 92%,respectively,which means the vortex induced vibration (VIV) could be effectively alleviated at certain larger values of λ/D m and a /D m.Also,the drag coefficients can be reduced by 30%.It is found that the flow field presents contrary patterns with the increase of λ/D m.The free shear layer becomes much more stable and rolls up into mature vortex only further downstream when λ/D m falls in the range of 5-7.The amplitude ratio a /D m greatly changes the separation line,and subsequently influences the wake structures.     

Control of vortex shedding in an axisymmetric bluff body wake

In the present experimental study the effect of a control disc mounted at the rear of an axisymmetric blunt-based body of revolution, first studied by Mair, is investigated in the Reynolds number range 3×10³≤Re_D≤5×10⁴ . As the distance of the control disc from the blunt base is increased, four vortex shedding regimes are identified: at small distances there is no effect, then a sharp increase of vortex shedding activity and total drag is observed, followed by an interval with reduced activity and drag and finally at large distances a regime where the flow around the main body and disc become essentially independent, i.e. where the drag forces of the two elements become additive. The near and far wake velocity fields corresponding to the different regimes are documented with time- and phase-averaged hot-wire and LDA measurements, with spectral analysis of the data and with flow visualizations of the near wake. The results are used to develop an improved understanding of the instability mechanism leading to high vortex shedding activity.