Three-dimensional transition in the wake of a transversely oscillating cylinder

2007 ◽  
Vol 577 ◽  
pp. 79-104 ◽  
Author(s):  
JUSTIN S. LEONTINI ◽  
M. C. THOMPSON ◽  
K. HOURIGAN

A Floquet stability analysis of the transition to three-dimensionality in the wake of a cylinder forced to oscillate transversely to the free stream has been undertaken. The effect of varying the oscillation amplitude is determined for a frequency of oscillation close to the natural shedding frequency. The three-dimensional modes that arise are identified, and the effect of the oscillation amplitude on their structure and growth rate quantified.It is shown that when the two-dimensional wake is in the 2S configuration (which is similar to the Kármán vortex street), the three-dimensional modes that arise are similar in nature and symmetry structure to the modes in the wake of a fixed cylinder. These modes are known as modes A, B and QP and occur in this order with increasing Re. However, increasing the amplitude of oscillation causes the critical Reynolds number for mode A to increase significantly, to the point where mode B becomes critical before mode A. The critical wavelength for mode A is also affected by the oscillation, becoming smaller with increasing amplitude. Elliptic instability theory is shown also to predict this trend, providing further support that mode A primarily arises as a result of an elliptic instability.At higher oscillation amplitudes, the spatio-temporal symmetry of the two-dimensional wake changes and it takes on the P + S configuration, with a pair of vortices on one side of the wake and a single vortex on the other side, for each oscillation cycle. With the onset of this configuration, modes A, B and QP cease to exist. It is shown that two new three-dimensional modes arise from this base flow, which we call modes SL and SS. Both of these modes are subharmonic, repeating over two base-flow periods. Also, either mode can be the first to become critical, depending on the amplitude of oscillation of the cylinder.The emergence of these two new modes, as well as the reversal of the order of inception of the three-dimensional modes A and B, leads to the observation that for an oscillating cylinder wake there are four different modes that can lead the transition to three-dimensionality, depending on the amplitude of oscillation. Therefore this type of flow provides a good example for studying the effect of mode-order inception on the path taken to turbulence in bluff-body wakes.For the range of amplitudes studied, the maximum Re value for which the flow remains two-dimensional is 280.

2013 ◽  
Vol 736 ◽  
pp. 316-350 ◽  
Author(s):  
J. A. Bourgeois ◽  
B. R. Noack ◽  
R. J. Martinuzzi

AbstractWe experimentally investigate the three-dimensional wake behind a finite wall-mounted square cylinder at $\mathit{Re}= 12\hspace{0.167em} 000$ and aspect ratio of 4. Focus is placed on the base flow and oscillatory fluctuation. Time-resolved three-dimensional velocity fields are constructed from high-frame-rate particle image velocimetry (PIV) and simultaneously recorded surface pressure measurements. All three velocity components are resolved in a rectangular near-wake region by two orthogonal dense arrays of parallel PIV planes. A key enabler is a generalized phase average incorporating a slowly varying base flow, a variable oscillation amplitude and higher harmonics. These generalizations reduce the instantaneous residual 30 % below those of a traditional phase average. Moreover, the resolved variations reveal analytical constraints of the mean flow and oscillation levels, such as the mean-field paraboloid. The proposed methodology for generalized phase averaging and for construction of three-dimensional velocity fields from two-dimensional PIV data is applicable to a large class of turbulent flows with oscillatory dynamics.


2013 ◽  
Vol 727 ◽  
pp. 346-378 ◽  
Author(s):  
Yongyun Hwang ◽  
Jinsung Kim ◽  
Haecheon Choi

AbstractControlling vortex shedding using spanwise-varying passive or active actuation (namely three-dimensional control) has recently been reported as a very efficient method for regulating two-dimensional bluff-body wakes. However, the mechanism of how the designed controller regulates vortex shedding is not clearly understood. To understand this mechanism, we perform a linear stability analysis of two-dimensional wakes, the base flow of which is modified with a given spanwise waviness. Absolute and convective instabilities of the spanwise wavy base flows are investigated using Floquet theory. Two types of base-flow modification are considered: varicose and sinuous. Both of these modifications attenuate absolute instability of two-dimensional wakes. In particular, the varicose modification is found to be much more effective in the attenuation than the sinuous one, and its spanwise lengths resulting in maximum attenuation show good agreement with those in three-dimensional controls. The physical mechanism of the stabilization is found to be associated with the formation of streamwise vortices from tilting of two-dimensional Kármán vortices and the subsequent tilting of these streamwise vortices by the spanwise shear in the base flow. Finally, the sensitivity of absolute instability to spanwise wavy base-flow modification is investigated. It is shown that absolute instability of two-dimensional wakes is much less sensitive to spanwise wavy base-flow modification than to two-dimensional modification. This suggests that the high efficiency observed in several three-dimensional controls is not due to the sensitive response of the wake instability to the spanwise waviness in the base flow.


2004 ◽  
Vol 126 (5) ◽  
pp. 861-870 ◽  
Author(s):  
A. Thakur ◽  
X. Liu ◽  
J. S. Marshall

An experimental and computational study is performed of the wake flow behind a single yawed cylinder and a pair of parallel yawed cylinders placed in tandem. The experiments are performed for a yawed cylinder and a pair of yawed cylinders towed in a tank. Laser-induced fluorescence is used for flow visualization and particle-image velocimetry is used for quantitative velocity and vorticity measurement. Computations are performed using a second-order accurate block-structured finite-volume method with periodic boundary conditions along the cylinder axis. Results are applied to assess the applicability of a quasi-two-dimensional approximation, which assumes that the flow field is the same for any slice of the flow over the cylinder cross section. For a single cylinder, it is found that the cylinder wake vortices approach a quasi-two-dimensional state away from the cylinder upstream end for all cases examined (in which the cylinder yaw angle covers the range 0⩽ϕ⩽60°). Within the upstream region, the vortex orientation is found to be influenced by the tank side-wall boundary condition relative to the cylinder. For the case of two parallel yawed cylinders, vortices shed from the upstream cylinder are found to remain nearly quasi-two-dimensional as they are advected back and reach within about a cylinder diameter from the face of the downstream cylinder. As the vortices advect closer to the cylinder, the vortex cores become highly deformed and wrap around the downstream cylinder face. Three-dimensional perturbations of the upstream vortices are amplified as the vortices impact upon the downstream cylinder, such that during the final stages of vortex impact the quasi-two-dimensional nature of the flow breaks down and the vorticity field for the impacting vortices acquire significant three-dimensional perturbations. Quasi-two-dimensional and fully three-dimensional computational results are compared to assess the accuracy of the quasi-two-dimensional approximation in prediction of drag and lift coefficients of the cylinders.


2009 ◽  
Vol 643 ◽  
pp. 349-362 ◽  
Author(s):  
DAVID LO JACONO ◽  
JUSTIN S. LEONTINI ◽  
MARK C. THOMPSON ◽  
JOHN SHERIDAN

A study of the flow past an oscillatory rotating cylinder has been conducted, where the frequency of oscillation has been matched to the natural frequency of the vortex street generated in the wake of a stationary cylinder, at Reynolds number 300. The focus is on the wake transition to three-dimensional flow and, in particular, the changes induced in this transition by the addition of the oscillatory rotation. Using Floquet stability analysis, it is found that the fine-scale three-dimensional mode that typically dominates the wake at a Reynolds number beyond that at the second transition to three-dimensional flow (referred to as mode B) is suppressed for amplitudes of rotation beyond a critical amplitude, in agreement with past studies. However, the rotation does not suppress the development of three-dimensionality completely, as other modes are discovered that would lead to three-dimensional flow. In particular, the longer-wavelength mode that leads the three-dimensional transition in the wake of a stationary cylinder (referred to as mode A) is left essentially unaffected at low amplitudes of rotation. At higher amplitudes of oscillation, mode A is also suppressed as the two-dimensional near wake changes in character from a single- to a double-row wake; however, another mode is predicted to render the flow three-dimensional, dubbed mode D (for double row). This mode has the same spatio-temporal symmetries as mode A.


2018 ◽  
Vol 841 ◽  
pp. 636-653
Author(s):  
Ting-Yueh Chang ◽  
Falin Chen ◽  
Min-Hsing Chang

A three-dimensional linear stability analysis is carried out for a convecting layer in which both the temperature and solute distributions are linear in the horizontal direction. The three-dimensional results show that, for $Le=3$ and 100, the most unstable mode occurs invariably as the longitudinal mode, a vortex roll with its axis perpendicular to the longitudinal plane, suggesting that the two-dimensional results are sufficient to illustrate the stability characteristics of the convecting layer. Two-dimensional results show that the stability boundaries of the transverse mode (a vortex roll with its axis perpendicular to the transverse plane) and the longitudinal modes are virtually overlapped in the regime dominated by thermal diffusion and the regime dominated by solute diffusion, while these two modes hold a significant difference in the regime the salt-finger instability prevails. More precisely, the instability area in terms of thermal Grashof number $Gr$ and solute Grashof number $Gs$ is larger for the longitudinal mode than the transverse mode, implying that, under any circumstance, the longitudinal mode is always more unstable than the transverse mode.


2019 ◽  
Vol 874 ◽  
pp. 299-338 ◽  
Author(s):  
Rui Wang ◽  
Yan Bao ◽  
Dai Zhou ◽  
Hongbo Zhu ◽  
Huan Ping ◽  
...  

In this paper, instabilities in the flow over a circular cylinder of diameter $D$ with dual splitter plates attached to its rear surface are numerically investigated using the spectral element method. The key parameters are the splitter plate length $L$, the attachment angle $\unicode[STIX]{x1D6FC}$ and the Reynolds number $Re$. The presence of the plates was found to significantly modify the flow topology, leading to substantial changes in both the primary and secondary instabilities. The results showed that the three instability modes present in the bare circular cylinder wake still exist in the wake of the present configurations and that, in general, the occurrences of modes A and B are delayed, while the onset of mode QP is earlier in the presence of the splitter plates. Furthermore, two new synchronous modes, referred to as mode A$^{\prime }$ and mode B$^{\prime }$, are found to develop in the wake. Mode A$^{\prime }$ is similar to mode A but with a quite long critical wavelength. Mode B$^{\prime }$ shares the same spatio-temporal symmetries as mode B but has a distinct spatial structure. With the exception of the case of $L/D=0.25$, mode A$^{\prime }$ persists for all configurations investigated here and always precedes the transition through mode A. The onset of mode B$^{\prime }$ occurs for $\unicode[STIX]{x1D6FC}>20^{\circ }$ with $L/D=1.0$ and for $L/D>0.5$ with $\unicode[STIX]{x1D6FC}=60^{\circ }$. The characteristics of all the transition modes are analysed, and their similarities and differences are discussed in detail in comparison with the existing modes. In addition, the physical mechanism responsible for the instability mode B$^{\prime }$ is proposed. The weakly nonlinear feature of mode B$^{\prime }$, as well as that of mode A$^{\prime }$, is assessed by employing the Landau model. Finally, selected three-dimensional simulations are performed to confirm the existence of these two new modes and to investigate the nonlinear evolution of the three-dimensional modes.


2017 ◽  
Vol 822 ◽  
pp. 813-847 ◽  
Author(s):  
Azan M. Sapardi ◽  
Wisam K. Hussam ◽  
Alban Pothérat ◽  
Gregory J. Sheard

This study seeks to characterise the breakdown of the steady two-dimensional solution in the flow around a 180-degree sharp bend to infinitesimal three-dimensional disturbances using a linear stability analysis. The stability analysis predicts that three-dimensional transition is via a synchronous instability of the steady flows. A highly accurate global linear stability analysis of the flow was conducted with Reynolds number $\mathit{Re}<1150$ and bend opening ratio (ratio of bend width to inlet height) $0.2\leqslant \unicode[STIX]{x1D6FD}\leqslant 5$. This range of $\mathit{Re}$ and $\unicode[STIX]{x1D6FD}$ captures both steady-state two-dimensional flow solutions and the inception of unsteady two-dimensional flow. For $0.2\leqslant \unicode[STIX]{x1D6FD}\leqslant 1$, the two-dimensional base flow transitions from steady to unsteady at higher Reynolds number as $\unicode[STIX]{x1D6FD}$ increases. The stability analysis shows that at the onset of instability, the base flow becomes three-dimensionally unstable in two different modes, namely a spanwise oscillating mode for $\unicode[STIX]{x1D6FD}=0.2$ and a spanwise synchronous mode for $\unicode[STIX]{x1D6FD}\geqslant 0.3$. The critical Reynolds number and the spanwise wavelength of perturbations increase as $\unicode[STIX]{x1D6FD}$ increases. For $1<\unicode[STIX]{x1D6FD}\leqslant 2$ both the critical Reynolds number for onset of unsteadiness and the spanwise wavelength decrease as $\unicode[STIX]{x1D6FD}$ increases. Finally, for $2<\unicode[STIX]{x1D6FD}\leqslant 5$, the critical Reynolds number and spanwise wavelength remain almost constant. The linear stability analysis also shows that the base flow becomes unstable to different three-dimensional modes depending on the opening ratio. The modes are found to be localised near the reattachment point of the first recirculation bubble.


2015 ◽  
Vol 25 (09) ◽  
pp. 1530022
Author(s):  
N. Nabatian ◽  
N. W. Mureithi

The two-dimensional numerical simulation of the flow over a cylinder forced to oscillate in the streamwise direction for Re = 200 is performed in CFX ANSYS. The controlled-vibration comprises of prescribed inline vibration from displacement amplitude-to-cylinder diameter A/D = 0.05 up to 0.5 with the excitation frequency ratios of 1, 1.5 and 2 including the harmonic and superharmonic excitation regions. The immersed boundary method is used to model the cylinder oscillation. Modal decomposition of the transverse velocity field via the proper orthogonal decomposition (POD) method is applied to uncover the interaction of symmetric and antisymmetric modes of the near wake. A model using the first two POD modes is developed based on symmetry group equivariance. The model predicts the mode interactions and bifurcated solution branches for all cases, and is shown to be in good agreement with numerical as well as previous experimental results. Lock-on is determined for a range of values of the oscillation amplitudes and frequency ratios. It is shown that the lock-on range widens with increasing nondimensional oscillation amplitude. The asymmetric 2S, P + S and symmetric pattern S with symbol S for a single vortex and P for a vortex pair shed per cycle, as well as a regime in which vortex formation is not synchronized with cylinder motion are observed in the cylinder wake depending on the combination of oscillation amplitude and frequency ratio. The frequency ratio variation from 1 to 2 leads to the switching from asymmetric to symmetric modes. The symmetric flow pattern corresponds to a near zero lift coefficient on the cylinder.


2008 ◽  
Vol 612 ◽  
pp. 237-260 ◽  
Author(s):  
JÉRÔME FONTANE ◽  
LAURENT JOLY

We perform a three-dimensional stability analysis of the Kelvin–Helmholtz (KH) billow, developing in a shear layer between two fluids with different density. We begin with two-dimensional simulations of the temporally evolving mixing layer, yielding the unsteady base flow fields. The Reynolds number is 1500 while the Schmidt and Froude numbers are infinite. Then exponentially unstable modes are extracted from a linear stability analysis performed at the saturation of the primary mode kinetic energy. The spectrum of the least stable modes exhibits two main classes. The first class comprises three-dimensional core-centred and braid-centred modes already present in the homogeneous case. The baroclinic vorticity concentration in the braid lying on the light side of the KH billow turns the flow into a sharp vorticity ridge holding high shear levels. The hyperbolic modes benefit from the enhanced level of shear in the braid whereas elliptic modes remain quite insensitive to the modifications of the base flow. In the second class, we found typical two-dimensional modes resulting from a shear instability of the curved vorticity-enhanced braid. For a density contrast of 0.5, the wavelength of the two-dimensional instability is about ten times shorter than that of the primary wave. Its amplification rate competes well against those of the hyperbolic three-dimensional modes. The vorticity-enhanced braid thus becomes the preferred location for the development of secondary instabilities. This stands as the key feature of the transition of the variable-density mixing layer. We carry out a fully resolved numerical continuation of the nonlinear development of the two-dimensional braid-mode. Secondary roll-ups due to a small-scale Kelvin–Helmholtz mechanism are promoted by the underlying strain field and develop rapidly in the compression part of the braid. Originally analysed by Reinoud et al. (Phys. Fluids, vol. 12, 2000, p. 2489) from two-dimensional non-viscous numerical simulations, this instability is shown to substantially increase the mixing.


1989 ◽  
Vol 200 ◽  
pp. 189-216 ◽  
Author(s):  
Arnon Chait ◽  
Seppo A. Korpela

The multicellular flow between two vertical parallel plates is numerically simulated using a time-splitting pseudospectral method. The steady flow of air, and the time-periodic flow of oil (Prandtl numbers of 0.71 and 1000, respectively) are investigated and descriptions of these flows using both physical and spectral approaches are presented. The details of the time dependency of the flow and temperature fields of oil are shown, and the dynamics of the process is discussed. The spectral transfer of energy among the axial modes comprising the flow is explored. The spectra of kinetic energy and thermal variance for air are found to be smooth and viscously dominated. Similar spectra for oil are bumpier, and the dynamics of the time-dependent flow are determined to be confined to the lower end of the spectrum alone.The three-dimensional linear stability of the multicellular flow of air is parametrically studied. The domain of stable two-dimensional cellular motion was found to be constrained by the Eckhaus instability and by two types of monotone instability. The two-dimensional multicellular flow is unstable above a Grashof number of about 8550 (with the critical Grashof number for the base flow being 8037). Therefore the flow of air in a sufficiently tall vertical enclosure should be considered to be three-dimensional for most practical applications.


Sign in / Sign up

Export Citation Format

Share Document