POD analysis of a finite-length cylinder near wake

The near wake of a wall-mounted finite-length square cylinder with an aspect ratio of 7 is investigated based on the proper orthogonal decomposition (POD) of the PIV data measured in three spanwise planes, i.e., z/d = 6, 3.5 and 1.0, near the cylinder free end, mid-span and fixed end (wall), respectively. The Reynolds number based on free-stream velocity (U _∞) and cylinder width (d) is 9,300. A two-dimensional (2D) square cylinder wake is also measured and analyzed at the same Reynolds number for the purpose of comparison. The structures of various POD modes show marked differences between the two flows. While the coefficients, a ₁ and a ₂, of the POD modes 1 and 2 occur within an annular area centered at a ₁ = a ₂ = 0 in the 2D wake, their counterparts are scattered all over the entire circular plane at z/d = 1.0 and 3.5 of the finite-length cylinder wake. Flow at z/d = 6 is dominated by POD mode 1, which corresponds to symmetrical vortex shedding and accounts for 54.0 % of the total turbulent kinetic energy (TKE). On the other hand, the POD modes 1 and 2, corresponding to anti-symmetrical vortex shedding, are predominant, accounting for about 45.0 % of the total TKE, at z/d = 3.5 and 1. It has been found that the flow structure may be qualitatively and quantitatively characterized by the POD coefficients. For example, at z/d = 6, a larger a ₁ corresponds to a smaller length of flow reversal zone and a stronger downwash flow. At z/d = 3.5 and 1, two typical flow modes can be identified from a ₁ and a ₂. While large a ₁ and/or a ₂ correspond to anti-symmetrical vortex shedding, as in a 2D cylinder wake, small a ₁ and a ₂ lead to symmetrical vortex shedding. Any values between the large and small a ₁ and/or a ₂ correspond to the flow structure between these two typical flow modes. As such, the probability of occurrence of a flow structure may be determined from the distribution of the POD coefficients.     

Effects of large spanwise wavelength on the wake of a sinusoidal wavy cylinder

The wake of a sinusoidal wavy cylinder with a large spanwise wavelength λ/D_m (=3.79–7.57) and a constant wave amplitude a/D_m=0.152, where D_m is the mean diameter of the cylinder, is investigated using three dimensional (3D) large eddy simulation (LES) at a subcritical Reynolds number Re=3×10³, based on incoming free-stream velocity (U_∞) and D_m. Attention is paid to assimilating the effects of λ/D_m on the cylinder wake, including vortex shedding frequency, spanwise vortex formation length, streamwise velocity distribution, flow separation angle, 3D vortex structure, and turbulent kinetic energy (TKE) distribution. Based on the predominant role of λ/D_m in the near wake modification, three regimes are identified, i.e., regime I at λ/D_m<6.0, regime II at λ/D_m≈6.0 and regime III at λ/D_m>6.0. A dramatic decrease in fluid forces is observed at λ/D_m=6.06, about 16% and 93% reduction in time-averaged drag and fluctuating lift, respectively, compared to those of a smooth cylinder. We identified, for the first time, an optimum λ/D_m (=6.06) for the wavy cylinder with relatively large λ/D_m (>3.5) in the subcritical flow regime. The underlying mechanisms of force reduction are discussed, including the flow characteristics at the three λ/D_m regimes. A comparison is also made between the results of λ/D_m effects on the near wakes of a circular and a square cylinder.     

Role of the shear layer instability in the near wake behavior of two side-by-side circular cylinders

Wakes, and their interaction behind two parallel cylinders lying in a plane perpendicular to the flow, have been investigated experimentally in the sub-critical Reynolds number regime. The experiments were performed in a water channel using laser Doppler velocimetry. The gap between the two cylinders was less than the cylinder diameter, a geometry referred to as strong interaction configuration. In this case the blockage is strong and a gap-jet appears between the cylinders. Two flow regimes of the near wake region have been identified: one below a critical Reynolds number Re _c ]1000;1700[, where the gap jet is stably deflected to one side and the double near-wake becomes asymmetric; the other, above Re _c, where the gap-jet deflection is unstable and a random flopping phenomenon takes place. When Re<Re _c, two different Strouhal numbers are identified, related to the Kármán vortex shedding behind each cylinder. When Re>Re _c, a third frequency appears in the near wake, related to the development of Kelvin-Helmholtz vortices in the separated shear layer of the cylinders [Prasad A, Williamson CHK (1997) J Fluid Mech 333:375]. The observed flopping behavior is attributed to the birth of these Kelvin-Helmholtz instabilities and their intermittent nature. Further downstream, beyond about five cylinder diameters, the random flopping flow phenomena disappear while a slightly asymmetric single wake persists. It is characterized by a Strouhal number St=0.13, a value that one would normally measure behind a single cylinder of twice its diameter.     

On the flow behaviour of confined finite-length wavy cylinders

Confined aspect-ratio of 6 wavy cylinders with a mean blockage-ratio of 0.5 were studied using time-resolved particle-image velocimetry at a sub-critical Reynolds number of 2700. Wavelengths and wave amplitudes of 2–4 and 0.1–0.3 mean diameters respectively were investigated. Results show that vortices are generally shed from the wavy cylinder and channel walls regularly, reminiscent of the unsteady symmetric flow configuration in confined non-wavy cylinders. Furthermore, vortex formation lengths for confined wavy cylinders are generally shorter than their unconfined counterparts, though their variations with respect to geometrical changes remain consistent with unconfined flow conditions. Gross cross-stream flow behaviour does not differ significantly between confined and unconfined wavy cylinders, indicating that finite-length effects are independent of the present confinement. Confined wavy cylinder wake regions are more sensitive towards geometrical changes and a combination of small wavelength and large wave amplitude leads to significant suppression of coherent cylinder and wall vortex-shedding. This is supported by phase-averaged flow reconstructions derived from Proper Orthogonal Decomposition analysis. Lastly, larger wave amplitudes lead to redistributions of dominant flow energy further downstream and to higher mode numbers, which suggests a causal link to the formation of stronger and more coherent streamwise vortices.     

Vortex and wake-induced vibrations of a tandem arrangement of two flexible circular cylinders with near wake interference

Results showing the dynamic response of a tandem arrangement of two vertical high aspect ratio (length over diameter) and low mass ratio (mass over mass of displaced fluid) flexible cylinders vibrating at low mode number are presented in this paper. Two circular cylinder models were aligned with the flow, so the downstream or trailing cylinder was immersed in the wake of the leading one. Centre-to-centre distances from 2 to 4 diameters were studied. The models were very similar in design, with external diameters of 16 mm and a total length of 1.5 m. Reynolds numbers up to 12 000 were achieved with reduced velocities, based on the fundamental natural frequency of the downstream cylinder in still water, up to 16. The trailing model had a mass ratio of 1.8 with a combined mass-damping parameter of 0.049, whilst the corresponding figures for the leading cylinder were 1.45 and 0.043, respectively. The dynamic response of the trailing model has been analysed by studying cross-flow and in-line amplitudes, dominant frequencies and modal amplitudes. The dynamic response of the leading one is analysed by means of its cross-flow amplitudes and dominant frequencies and it is also related to the motion of the trailing cylinder by studying the synchronisation between their instantaneous cross-flow motions. Planar digital particle image velocimetry (DPIV) was used to visualise the wake. Different response regimes have been identified based on the type of oscillations exhibited by the cylinders: vortex-induced (VIV), wake-induced (WIV) or combinations of both.     

The near wake structure of a square cylinder

The near wake structure of a square cross section cylinder in flow perpendicular to its length was investigated experimentally over a Reynolds number (based on cylinder width) range of 6700–43,000. The wake structure and the characteristics of the instability wave, scaling on θ at separation, were strongly dependent on the incidence angle () of the freestream velocity. The nondimensional frequency (St_θ) of the instability wave varied within the range predicted for laminar instability frequencies for flat plate wakes, jets and shear layers. For = 22.5°, the freestream velocity was accelerated over the side walls and the deflection of the streamlines (from both sides of the cylinder) towards the center line was higher compared to the streamlines for = 0°. This caused the vortices from both sides of the cylinder to merge by x/d 2, giving the mean velocity distribution typical of a wake profile. For = 0°, the vortices shed from both sides of the cylinder did not merge until x/d 4.5. The separation boundary layer for all cases was either transitional or turbulent, yet the results showed good qualitative, and for some cases even quantitative, agreement with linearized stability results for small amplitude disturbances waves in laminar separation layers.     

An experimental study of interacting coherent structures in the turbulent near wake behind a pair of square cylinders

An experimental investigation (LDA technique) of coherent structures in the turbulent plane near wake behind a pair of square cylinders in side-by-side arrangement is presented with particular emphasis put on the use of 2D invariant structural measures.     

Obtaining phase averaged turbulence properties in the near wake of a circular cylinder at high Reynolds number using POD

The flow past a circular cylinder at high Reynolds number is studied by means of PIV, 3C-PIV and Time-Resolved PIV techniques. One of the goals of this study was to allow comparisons with numerical simulations on a domain identical to that of the experiment. For this reason, the cylinder was placed in a confined environment, with a high blockage and a low aspect ratio, thereby allowing computations on a mesh of reasonable size, and avoiding “infinite conditions”. This paper deals with the decomposition of the flow in a coherent and random parts. To this aim, phase averaged quantities were first obtained using the wall pressure signal on the cylinder as a trigger signal. This was achieved using both conditional sampling and LSE with similar results. This decomposition is then analysed using the Time Resolved PIV measurements, as well as by comparison of the contributions of the organised and turbulent fluctuations to the time-independent Reynolds stress tensor with those estimated from velocity spectra by interpolation and integration of the continuous part. In agreement with other studies, it is found that the contribution of the turbulent motion is overestimated as a result of the occurence of phase jitter between the trigger and velocity signal. A POD analysis was then performed to extract the coherent motion and to compare this decomposition with that obtained by phase averaging. Similarly to the phase averaging, the POD allows the decomposition of the time-independent stress tensor as the sum of two contributions corresponding to the first N modes, and the rest of the modes. This decomposition is then analysed by comparing these contributions to those obtained from the velocity spectra, according to the value N chosen. It is found that these contributions are strongly dependent on N, and the contribution of the first modes greatly overestimate the coherent motion if N is too large. In order to obtain a good decomposition of the flow in coherent and random parts, the difficulty in this case lies in the choice of the modes. Finally, the POD coefficients of the first two modes are used instead of the pressure signal to determine the phase of the vortex shedding, and the phase averaging is reconsidered. It is found that the phase averaged vortices are less smeared by the averaging process, the turbulent stresses better follow the evolution of the vortices, and the contributions of both coherent and turbulent fluctuations are found to agree well with those evaluated from the velocity spectra. This enhancement is obtained because the phase angle is determined directly from the velocity fields to be averaged, thereby reducing the phase-jitter effect. A comparison with a detached eddy simulation is also briefly shown and demonstrates the high level of agreement obtainable between simulation and experiment, as well as confirming the enhancement of the phase averaging using this procedure.     

Strouhal numbers in the wake of two inline cylinders

The dominant vortex frequencies f _s in the wake of two tandem circular cylinders of identical diameter d have been measured simultaneously using two hot wires placed behind each cylinder. Measurements were conducted over the Reynolds number Re (U d/, where U and are the free-stream velocity and the kinematic viscosity of fluid, respectively) =800–4.2×10⁴ and the cylinder centre-to-centre spacing L/d=1–15. The Strouhal number St (f _s d/U ) exhibits a strong dependence on L/d and Re. For L/d<(L/d)_c, which is a critical value and ranges between 3.5 and 5, there is no vortex street formed in the gap between the cylinders, and St measured behind the downstream cylinder drops rapidly for increasing L/d. For L/d>(L/d)_c, co-shedding occurs, that is, vortices are shed from the upstream as well as the downstream cylinder, and their frequencies are found to be identical. St climbs with increasing L/d, approaching a constant between 0.18 and 0.22 for L/d>10. The St–Re relationship is classified into four categories, based on their behaviours, which are associated with distinct flow physics—category 1: for 1L/d<2, the shear layers separated from the upstream cylinder roll up behind the downstream cylinder; category 2: for 2L/d3, there is a transition from the shear layer rollup behind to reattachment on the downstream cylinder; category 3: for 3<L/d5, transition from the reattachment to co-shedding regime occurs at a critical Reynolds number; and category 4: the flow for L/d>5 is characterized by co-shedding only. The present measurements reconfirm the previous observation of a bi-stable flow at the transition from the reattachment to co-shedding regime. It is found for the first time that another bi-stable flow occurs at the transition from category 1 to 2, that is, the stable reattachment co-exists with the stable rollup (behind the downstream cylinder) of shear layers separating from the upstream cylinder, resulting in two distinct vortex-shedding frequencies even at the same Re and L/d. The St behaviour is further discussed along with flow visualization using the laser-induced fluorescence technique.     

Elastic analyses of heterogeneous hollow cylinders

Two different kinds of heterogeneous elastic hollow cylinders are studied in the present paper. One is a multi-layered cylinder with different values in different layers for both elastic modulus and Poisson’s ratio. Another is an elastic hollow cylinder with continuously graded material properties. By introducing two recursive algorithms, the extrusion stresses between two neighbor layers in the multi-layered cylinder submitted to uniform pressures on the inner and outer surfaces can be simply determined. Then the exact solutions of the multi-layered structure can be found based on Lamé’s solution. For the hollow cylinder with continuously graded properties, the displacement method is used. Both Whittaker equation and hyper-geometric equation are derived and successfully solved, and then the exact solutions are found. The results obtained in the present paper are compared with the numerical solutions and good agreements are found. At the end of the present paper, some inherent properties of these two different kinds of heterogeneous elastic hollow cylinders are presented and discussed. The results obtained in the present paper are useful in the design and analysis for composites reinforced by unidirectional fiber layers.     

Laminar heat transfer characteristics of internally finned tube with sinusoidal wavy fin

Comparative numerical study of laminar heat transfer characteristics of annular tubes with sinusoidal wavy fins has been conducted both experimentally and numerically with Re = 299–1,475. The uniform heat flux is imposed on the tube outside wall surface. Two tube materials (copper and stainless steel) are considered. It is found that the fluid temperature profile is not linear but convex along the flow direction due to the axial heat conduction in tube wall, and the effects of axial heat conduction on the heat transfer decreases with an increase in Reynolds number or decrease in tube wall thermal conductivity. The axial distributions of local Nusselt number could reach periodically fully developed after 3–5 cycles. The convectional data reduction method based on the traditional method should be improved for tube with high thermal conductivity or low Reynolds numbers, Otherwise, the heat transfer performance of internally finned tube may be underestimated.     

Unsteady motion of a profile near a wavy shield

The investigation of the motion of a profile near a plane shield and near a wavy shield is one of the ways of approximately taking into account the effect of the surface of a fluid on the characteristics of surface wings. The steady-state characteristics of the profiles near a plane shield have been investigated in a number of theoretical and experimental studies [1, 2]. In the present study we consider the unsteady motion of a thin profile near plane and wavy shields in an ideal incompressible fluid at rest at infinity, and when there is a fluid flow.Translated from Izvestiya Akadeinii Nauk SSSR. Mekhanika Zhidkosti i Gaza, No. 3, pp. 10–16, May–June, 1973.     

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.     

The flow regime in the turbulent near wake of a hemisphere

The work presented is a wind tunnel study of the near wake region behind a hemisphere immersed in three different turbulent boundary layers. In particular, the effect of different boundary layer profiles on the generation and distribution of near wake vorticity and on the mean recirculation region is examined. Visualization of the flow around a hemisphere has been undertaken, using models in a water channel, in order to obtain qualitative information concerning the wake structure.List of symbols C _p pressure coefficient, - D diameter of hemisphere - n vortex shedding frequency - p pressure on model surface - p ₀ static pressure - Re Reynolds number, - St Strouhal number, - U, V, W local mean velocity components - mean freestream velocity inX direction - U ^* shear velocity, - u, v, w velocity fluctuations inX, Y andZ directions - X Cartesian coordinate in longitudinal direction - Y Cartesian coordinate in lateral direction - Z Cartesian coordinate in direction perpendicular to the wall - it* boundary layer displacement thickness, - diameter of model surface roughness - elevation angleI - O boundary layer momentum thickness, - _w wall shearing stress - dynamic viscosity of fluid - density of fluid - streamfunction - _x longitudinal component of vorticity, - _y lateral component of vorticity, - _z vertical component of vorticity, This paper was presented at the Ninth symposium on turbulence, University of Missouri-Rolla, October 1–3, 1984     

Experimental study on wake flow structures of screen cylinders using PIV

Experiments were conducted in a water flume using Particle Image Velocimetry (PIV) to study the evolution of the vortical structures in the wakes of four types of screen cylinders at a Reynolds number of about 3200. The results were compared with that of a bare cylinder. The screen cylinders were made of stainless steel screen meshes of various porosities (37%, 48%, 61% and 67%) rolled into cylindrical shapes. Smoke wire flow visualisations in a wind tunnel were also conducted in support of the PIV tests. Depending on the porosity of the screen mesh, two vortex formation mechanisms for the screen cylinder wakes were identified. One was associated with wake instability and the other was associated with shear-layer (Kelvin-Helmholtz) convective instability which involved merging through pairing and tripling of small-scale vortices within the shear layers. The former was responsible for the formation of large-scale vortices in the bare cylinder and the screen cylinder wakes with 37% and 48% porosities, while the latter was responsible for the screen cylinder wakes with 61% and 67% porosities. The results also showed that with increasing porosity, the vortex formation region was extended farther downstream and the Reynolds shear stress, the Turbulent Kinetic Energy (TKE) and vortex intensity were decreased constantly.     

The velocity field of the turbulent very near wake of a circular cylinder

Hot-wire measurements were conducted in the very near wake (x/d10) of a circular cylinder at a Reynolds number based on cylinder diameter, Re _d of 3900. Measurements of the streamwise velocity component with the use of single sensor hot-wire probes were found to be inaccurate for such flowfields where high flow angles are present. An X-array probe provided detailed streamwise and lateral velocity component statistics. Frequency spectra of these two velocity components are also presented. Measurements with a 4-sensor hot-wire probe confirmed that the very near wake region is dominantly two-dimensional, thus validating the accuracy of the present X-array data.This study has been funded by the NASA-Ames University Consortium Cooperative Agreement, NCC2-5003. We wish to thank Patrick Beaudan for providing us with the LES results for comparison and Parviz Moin for his interest in and encouragement of this experiment to provide validation data for the LES. We also wish to thank loseph Murray for his help with the look-up-table data reduction program.     

Strain gradient near interface of coaxial cylinders in torsion

Classical elastoplastic theory predicts that the rotation angle near an interface between two mismatched materials is discontinuous under shear. The strain gradient effects, however, can be significant within a narrow region near the interface. This can be shown by application of the strain gradient plasticity. The matching expansion method was used to obtain asymptotic results. Comparison is then made with those found numerically for the interface torsion problem of a two-layered cylindrical tube. The strain gradient plasticity theory solution differs from that of the classical elastoplastic theory solution, depending on the properties aside from the interface behavior and the loading mode. A failure criterion is also proposed that accounts for the strain gradients.     

Influence of wall proximity on characteristics of wake behind a square cylinder: PIV measurements and POD analysis

Influence of wall proximity on characteristics of the wake behind a two-dimensional square cylinder was experimentally studied in the present work. A low-speed recirculation water channel was established for the experiment; the Reynolds number based on the free-stream velocity and cylinder width (D) was kept at Re_D = 2250. Four cases with different gap width, e.g., G/D = 0.1, 0.2, 0.4 and 0.8, were chosen for comparison. Two experimental techniques, e.g., the standard PIV with high image-density CCD camera and TR-PIV with a high-speed camera were employed in measuring the wake field, enabling a comprehensive view of the time-averaged wake pattern at high spatial resolution and the instantaneous flow field at high temporal resolution, respectively. For the four cases, the difference in spatial characteristics of the wake in the vicinity of the plane wall was analyzed in terms of the time-averaged quantities measured by the standard PIV, e.g., the streamline pattern, the vector field, the streamwise velocity fluctuation intensity and the reverse-flow intermittency. The proper orthogonal decomposition (POD) method was extensively used to decompose the TR-PIV measurements, giving a close-up view of the energetic POD modes buried in the wake. The low-order flow model of the wake at G/D = 0.8 and 0.4 was constructed by using the linear combination of the first two POD modes and the time-mean flow field, which reflected well the vortex shedding process in the sense of the phase-dependent patterns. The intermittent appearance of the weakly separated region near the wall was found at G/D = 0.4. On going from G/D = 0.8 to 0.4, the remarkable variation of the instantaneous wake in the longitudinal direction confirmed that the wall constraint stretches the vortices in the plane of the wall and transfers the energy to the longitudinal component at the expense of the lateral one.     

Shear flow induced vibrations of long slender cylinders with a wake oscillator model

A time domain model is presented to study the vibrations of long slender cylinders placed in shear flow. Long slender cylinders such as risers and tension legs are widely used in the field of ocean engineering. They are subjected to vortex-induced vibrations(VIV) when placed within a transverse incident flow. A three dimensional model coupled with wake oscillators is formulated to describe the response of the slender cylinder in cross-flow and in-line directions. The wake oscillators are distributed along the cylinder and the vortex-shedding frequency is derived from the local current velocity. A non-linear fiuid force model is accounted for the coupled effect between cross-flow and in-line vibrations. The comparisons with the published experimental data show that the dynamic features of VIV of long slender cylinder placed in shear flow can be obtained by the proposed model,such as the spanwise average displacement,vibration frequency,dominant mode and the combination of standing and traveling waves. The simulation in a uniform flow is also conducted and the result is compared with the case of nonuniform flow. It is concluded that the flow shear characteristic has significantly changed the cylinder vibration behavior.