This paper presents three-dimensional (3D) direct numerical simulations (DNS) of flow past a circular cylinder over a range of Reynolds number ($Re$) up to 300. The gradual wake transition process from mode A* (i.e. mode A with large-scale vortex dislocations) to mode B is well captured over a range of $Re$ from 230 to 260. The mode swapping process is investigated in detail with the aid of numerical flow visualization. It is found that the mode B structures in the transition process are developed based on the streamwise vortices of mode A or A* which destabilize the braid shear layer region. For each case within the transition range, the transient mode swapping process consists of dislocation and non-dislocation cycles. With the increase of $Re$, it becomes more difficult to trigger dislocations from the pure mode A structure and form a dislocation cycle, and each dislocation stage becomes shorter in duration, resulting in a continuous decrease in the probability of occurrence of mode A* and a continuous increase in the probability of occurrence of mode B. The occurrence of mode A* results in a relatively strong flow three-dimensionality. A critical condition is confirmed at approximately $Re=265{-}270$, where the weakest flow three-dimensionality is observed, marking a transition from the disappearance of mode A* to the emergence of increasingly disordered mode B structures.

With the aid of direct numerical simulation, this paper presents a detailed investigation on the flow around a finite square cylinder at a fixed aspect ratio (AR) of 4 and six Reynolds numbers (Re = 50, 100, 150, 250, 500, and 1000). It is found that the mean streamwise vortex structure is also affected by Re, apart from the AR value. Three types of mean streamwise vortices have been identified and analyzed in detail, namely, “Quadrupole Type” at Re = 50 and Re = 100, “Six-Vortices Type” at Re = 150 and Re = 250, and “Dipole Type” at Re = 500 and Re = 1000. It is the first time that the “Six-Vortices Type” mean streamwise vortices are reported, which is considered as a transitional structure between the other two types. Besides, three kinds of spanwise vortex-shedding models have been observed in this study, namely, “Hairpin Vortex Model” at Re = 150, “C and Reverse-C and Hairpin Vortex Model (Symmetric Shedding)” at Re = 250, and “C and Reverse-C and Hairpin Vortex Model (Symmetric/Antisymmetric Shedding)” at Re = 500 and Re = 1000. The newly proposed “C and Reverse-C and Hairpin Vortex Model” shares some similarities with “Wang’s Model” [H. F. Wang and Y. Zhou, “The finite-length square cylinder near wake,” J. Fluid Mech. 638, 453–490 (2009)] but differs in aspects such as the absence of the connection line near the free-end and the “C-Shape” vortex structure in the early stage of the formation of the spanwise vortex.

This paper presents the response and the wake modes of a freely vibrating D-section prism with varying angles of attack ( $\alpha = 0^\circ \text {--}180^\circ$ ) and reduced velocity ( $U^* = 2\text {--}20$ ) by a numerical investigation. The Reynolds number, based on the effective diameter, is fixed at 100. The vibration of the prism is allowed only in the transverse direction. We found six types of response with increasing angle of attack: typical vortex-induced vibration (VIV) at $\alpha = 0^\circ \text {--}35^\circ$ ; extended VIV at $\alpha = 40^\circ \text {--}65^\circ$ ; combined VIV and galloping at $\alpha = 70^\circ \text {--}80^\circ$ ; narrowed VIV at $\alpha = 85^\circ \text {--}150^\circ$ ; transition response, from narrowed VIV to pure galloping, at $\alpha = 155^\circ \text {--}160^\circ$ ; and pure galloping at $\alpha = 165^\circ \text {--}180^\circ$ . The typical and narrowed VIVs are characterized by linearly increasing normalized vibration frequency with increasing $U^*$ , which is attributed to the stationary separation points of the boundary layer. On the other hand, in the extended VIV, the vortex shedding frequency matches the natural frequency in a large $U^*$ range with increasing $\alpha$ generally. The galloping is characterized by monotonically increasing amplitude with enlarging $U^*$ , with the largest amplitude being $A^* = 3.2$ . For the combined VIV and galloping, the vibration amplitude is marginal in the VIV branch while it significantly increases with $U^*$ in the galloping branch. In the transition from narrowed VIV to pure galloping, the vibration frequency shows a galloping-like feature, but the amplitude does not monotonically increase with increasing $U^*$ . Moreover, a partition of the wake modes in the $U^*$ – $\alpha$ parametric plane is presented, and the flow physics is elucidated through time variations of the displacement, drag and lift coefficients and vortex dynamics. The angle-of-attack range of galloping is largely predicted by performing a quasi-steady analysis of the galloping instability. Finally, the effects of $m^*$ and ${\textit {Re}}$ , the roles of afterbody and the roles of separation point in determining vibration responses and vortex shedding frequency are further discussed.

Sinusoidally oscillatory flow around four circular cylinders in an in-line square arrangement is numerically investigated at Keulegan–Carpenter numbers ($\mathit{KC}$) ranging from 1 to 12 and at Reynolds numbers ($\mathit{Re}$) from 20 to 200. A set of flow patterns is observed and classified based on known oscillatory flow regimes around a single cylinder. These include six types of reflection symmetry regimes to the axis of flow oscillation, two types of spatio-temporal symmetry regimes and a series of symmetry-breaking flow patterns. In general, at small gap distances, the four structures behave more like a single body, and the flow fields therefore resemble those around a single cylinder with a large effective cylinder diameter. With increasing gap distance, flow structures around each individual cylinder in the array start to influence the overall flow patterns, and the flow field shows a variety of symmetry and asymmetry patterns as a result of vortex and shear layer interactions. The characteristics of hydrodynamic forces on individual cylinders as well as on the cylinder group are also examined. It is found that the hydrodynamic forces respond in a similar manner to the flow field to the cylinder proximity and wake interference.

No abstract

Two-dimensional (2D) and three-dimensional (3D) instabilities in the wake of a circular cylinder placed near to a moving wall are investigated using direct numerical simulation (DNS). The study covers a parameter space spanning a non-dimensional gap ratio ($G^{\ast }$) between 0.1 to 19.5 and Reynolds number ($Re$) up to 300. Variations in the flow characteristics with $Re$ and $G^{\ast }$ are studied, and their correlations with the hydrodynamic forces on the cylinder are investigated. It is also found that the monotonic increase of the critical $Re$ for 2D instability ($Re_{cr2D}$) with decreasing $G^{\ast }$ is influenced by variations in the mean flow rate around the cylinder, the confinement of the near-wake flow by the plane wall and the characteristics of the shear layer formed above the moving wall directly below the cylinder. The first factor destabilizes the wake flow at a moderate $G^{\ast }$ while the latter two factors stabilize the wake flow with decreasing $G^{\ast }$. In terms of 3D instability, the flow transition sequence of ‘2D steady $\rightarrow$ 3D steady $\rightarrow$ 3D unsteady’ for small gap ratios is analysed at $G^{\ast }=0.2$. It is found that the 3D steady and 3D unsteady flows are triggered by Mode C instability due to wall proximity. However, the Mode C structure is not sustained indefinitely, since interference with the shear layer leads to other 3D steady and unsteady flow structures.

Three-dimensional (3-D) wake transition for flow past a square cylinder aligned with sides perpendicular and parallel to the approaching flow is investigated using direct numerical simulation. The secondary wake instability, namely a Mode A instability, occurs at a Reynolds number ($Re$) of 165.7. A gradual wake transition from Mode A* (i.e. Mode A with vortex dislocations) to Mode B is observed over a range of $Re$ from 185 to 210, within which the probability of occurrence of vortex dislocations decreases monotonically with increasing $Re$. The characteristics of the Strouhal–Reynolds number relationship are analysed. At the onset of Mode A*, a sudden drop of the 3-D Strouhal number from its two-dimensional counterpart is observed, which is due to the subcritical nature of the Mode A* instability. A continuous 3-D Strouhal–Reynolds number curve is observed over the mode swapping regime, since Mode A* and Mode B have extremely close vortex shedding frequencies and therefore only a single merged peak is observed in the frequency spectrum. The existence of hysteresis for the Mode A and Mode B wake instabilities is examined. The unconfined Mode A and Mode B wake instabilities are hysteretic and non-hysteretic, respectively. However, a spanwise confined Mode A could be non-hysteretic. It is proposed that the existence of hysteresis at a wake instability can be identified by examining the sudden/gradual variation of the 3-D flow properties at the onset of the wake instability, with sudden and gradual variations corresponding to hysteretic (subcritical) and non-hysteretic (supercritical) flows, respectively.

The Honji instability is studied using direct numerical simulations of sinusoidal oscillatory flow around a circular cylinder. The three-dimensional Navier–Stokes equations are solved by a finite element method at a relatively small value of the Keulegan–Carpenter number KC. The generation and subsequent development of Honji vortices are discussed over a range of frequency parameters by means of flow visualization. It is found that the spacing between Honji vortices is only weakly dependent on the frequency of oscillation, but is strongly correlated to KC because it is the terms within the governing equation containing KC that dominate the three-dimensional features of the flow. An empirical relationship between KC and the spacing between neighbouring vortices is proposed. The three-dimensional steady streaming structure within the vortices is identified and it is found that at high frequencies the steady streaming is two-dimensional although the instantaneous flow structure is itself fully three-dimensional.