Numerical investigation of wake flow regimes behind a high-speed rotating circular cylinder in steady flow

2019 ◽  
Vol 878 ◽  
pp. 875-906
Author(s):  
Adnan Munir ◽  
Ming Zhao ◽  
Helen Wu ◽  
Lin Lu

Flow around a high-speed rotating circular cylinder for $Re\leqslant 500$ is investigated numerically. The Reynolds number is defined as $UD/\unicode[STIX]{x1D708}$ with $U$, $D$ and $\unicode[STIX]{x1D708}$ being the free-stream flow velocity, the diameter of the cylinder and the kinematic viscosity of the fluid, respectively. The aim of this study is to investigate the effect of a high rotation rate on the wake flow for a range of Reynolds numbers. Simulations are performed for Reynolds numbers of 100, 150, 200, 250 and 500 and a wide range of rotation rates from 1.6 to 6 with an increment of 0.2. Rotation rate is the ratio of the rotational speed of the cylinder surface to the incoming fluid velocity. A systematic study is performed to investigate the effect of rotation rate on the flow transition to different flow regimes. It is found that there is a transition from a two-dimensional vortex shedding mode to no vortex shedding mode when the rotation rate is increased beyond a critical value for Reynolds numbers between 100 and 200. Further increase in rotation rate results in a transition to three-dimensional flow which is characterized by the presence of finger-shaped (FV) vortices that elongate in the wake of the cylinder and very weak ring-shaped vortices (RV) that wrap the surface of the cylinder. The no vortex shedding mode is not observed at Reynolds numbers greater than or equal to 250 since the flow remains three-dimensional. As the rotation rate is increased further, the occurrence frequency and size of the ring-shaped vortices increases and the flow is dominated by RVs. The RVs become bigger in size and the flow becomes chaotic with increasing rotation rate. A detailed analysis of the flow structures shows that the vortices always exist in pairs and the strength of separated shear layers increases with the increase of rotation rate. A map of flow regimes on a plane of Reynolds number and rotation rate is presented.

2016 ◽  
Vol 792 ◽  
pp. 682-711 ◽  
Author(s):  
Michael O. John ◽  
Dominik Obrist ◽  
Leonhard Kleiser

The leading-edge boundary layer (LEBL) in the front part of swept airplane wings is prone to three-dimensional subcritical instability, which may lead to bypass transition. The resulting increase of airplane drag and fuel consumption implies a negative environmental impact. In the present paper, we present a temporal biglobal secondary stability analysis (SSA) and direct numerical simulations (DNS) of this flow to investigate a subcritical transition mechanism. The LEBL is modelled by the swept Hiemenz boundary layer (SHBL), with and without wall suction. We introduce a pair of steady, counter-rotating, streamwise vortices next to the attachment line as a generic primary disturbance. This generates a high-speed streak, which evolves slowly in the streamwise direction. The SSA predicts that this flow is unstable to secondary, time-dependent perturbations. We report the upper branch of the secondary neutral curve and describe numerous eigenmodes located inside the shear layers surrounding the primary high-speed streak and the vortices. We find secondary flow instability at Reynolds numbers as low as$Re\approx 175$, i.e. far below the linear critical Reynolds number$Re_{crit}\approx 583$of the SHBL. This secondary modal instability is confirmed by our three-dimensional DNS. Furthermore, these simulations show that the modes may grow until nonlinear processes lead to breakdown to turbulent flow for Reynolds numbers above$Re_{tr}\approx 250$. The three-dimensional mode shapes, growth rates, and the frequency dependence of the secondary eigenmodes found by SSA and the DNS results are in close agreement with each other. The transition Reynolds number$Re_{tr}\approx 250$at zero suction and its increase with wall suction closely coincide with experimental and numerical results from the literature. We conclude that the secondary instability and the transition scenario presented in this paper may serve as a possible explanation for the well-known subcritical transition observed in the leading-edge boundary layer.


Author(s):  
S. Nagaya ◽  
R. E. Baddour

CFD simulations of crossflows around a 2-D circular cylinder and the resulting vortex shedding from the cylinder are conducted in the present study. The capability of the CFD solver for vortex shedding simulation from a circular cylinder is validated in terms of the induced drag and lifting forces and associated Strouhal numbers computations. The validations are done for uniform horizontal fluid flows at various Reynolds numbers in the range 103 to 5×105. Crossflows around the circular cylinder beneath a free surface are also simulated in order to investigate the characteristics of the interaction between vortex shedding and a free surface at Reynolds number 5×105. The influence of the presence of the free surface on the vortex shedding due to the cylinder is discussed.


1990 ◽  
Vol 220 ◽  
pp. 459-484 ◽  
Author(s):  
H. M. Badr ◽  
M. Coutanceau ◽  
S. C. R. Dennis ◽  
C. Ménard

The unsteady flow past a circular cylinder which starts translating and rotating impulsively from rest in a viscous fluid is investigated both theoretically and experimentally in the Reynolds number range 103 [les ] R [les ] 104 and for rotational to translational surface speed ratios between 0.5 and 3. The theoretical study is based on numerical solutions of the two-dimensional unsteady Navier–Stokes equations while the experimental investigation is based on visualization of the flow using very fine suspended particles. The object of the study is to examine the effect of increase of rotation on the flow structure. There is excellent agreement between the numerical and experimental results for all speed ratios considered, except in the case of the highest rotation rate. Here three-dimensional effects become more pronounced in the experiments and the laminar flow breaks down, while the calculated flow starts to approach a steady state. For lower rotation rates a periodic structure of vortex evolution and shedding develops in the calculations which is repeated exactly as time advances. Another feature of the calculations is the discrepancy in the lift and drag forces at high Reynolds numbers resulting from solving the boundary-layer limit of the equations of motion rather than the full Navier–Stokes equations. Typical results are given for selected values of the Reynolds number and rotation rate.


2017 ◽  
Vol 829 ◽  
pp. 486-511 ◽  
Author(s):  
K. W. L. Wong ◽  
J. Zhao ◽  
D. Lo Jacono ◽  
M. C. Thompson ◽  
J. Sheridan

While flow-induced vibration of bluff bodies has been extensively studied over the last half-century, only limited attention has been given to flow-induced vibration of elastically mounted rotating cylinders. Since recent low-Reynolds-number numerical work suggests that rotation can enhance or suppress the natural oscillatory response, the former could find applications in energy harvesting and the latter in vibration control. The present experimental investigation characterises the dynamic response and wake structure of a rotating circular cylinder undergoing vortex-induced vibration at a low mass ratio ($m^{\ast }=5.78$) over the reduced velocity range leading to strong oscillations. The experiments were conducted in a free-surface water channel with the cylinder vertically mounted and attached to a motor that provided constant rotation. Springs and an air-bearing system allow the cylinder to undertake low-damped transverse oscillations. Under cylinder rotation, the normalised frequency response was found to be comparable to that of a freely vibrating non-rotating cylinder. At reduced velocities consistent with the upper branch of a non-rotating transversely oscillating cylinder, the maximum oscillation amplitude increased with non-dimensional rotation rate up to $\unicode[STIX]{x1D6FC}\approx 2$. Beyond this, there was a sharp decrease in amplitude. Notably, this critical value corresponds approximately to the rotation rate at which vortex shedding ceases for a non-oscillating rotating cylinder. Remarkably, at $\unicode[STIX]{x1D6FC}=2$ there was approximately an 80 % increase in the peak amplitude response compared to that of a non-rotating cylinder. The observed amplitude response measured over the Reynolds-number range of ($1100\lesssim Re\lesssim 6300$) is significantly different from numerical predictions and other experimental results recorded at significantly lower Reynolds numbers.


2018 ◽  
Vol 847 ◽  
pp. 664-699 ◽  
Author(s):  
J. Zhao ◽  
D. Lo Jacono ◽  
J. Sheridan ◽  
K. Hourigan ◽  
M. C. Thompson

This study experimentally investigates the in-line flow-induced vibration (FIV) of an elastically mounted circular cylinder under forced axial rotation in a free stream. The present experiments characterise the structural vibration, fluid forces and wake structure of the fluid–structure system at a low mass ratio (the ratio of the total mass to the displaced fluid mass) over a wide parameter space spanning the reduced velocity range $5\leqslant U^{\ast }\leqslant 32$ and the rotation rate range $0\leqslant \unicode[STIX]{x1D6FC}\leqslant 3.5$, where $U^{\ast }=U/(\,f_{nw}D)$ and $\unicode[STIX]{x1D6FC}=|\unicode[STIX]{x1D6FA}|D/(2U)$, with $U$ the free-stream velocity, $D$ the cylinder outer diameter, $f_{nw}$ the natural frequency of the system in quiescent water and $|\unicode[STIX]{x1D6FA}|$ the angular velocity of the cylinder rotation. The corresponding Reynolds number (defined by $Re=UD/\unicode[STIX]{x1D708}$, with $\unicode[STIX]{x1D708}$ the kinematic viscosity of the fluid) was varied over the interval $1349\leqslant Re\leqslant 8624$, where it is expected that the FIV response is likely to be relatively insensitive to the Reynolds number. The fluid–structure system was modelled using a low-friction air-bearing system in conjunction with a free-surface water-channel facility. Three vibration regions that exhibited vortex-induced vibration (VIV) synchronisation, rotation-induced galloping and desynchronised responses were observed. In both the VIV synchronisation and rotation-induced galloping regions, significant cylinder vibration was found to be correlated with wake–body synchronisation within the rotation rate range $2.20\lesssim \unicode[STIX]{x1D6FC}\lesssim 3.15$. Of significant interest, the frequency of the streamwise fluid force could be modulated by the imposed rotation to match that of the transverse lift force, resulting in harmonic synchronisation. Measurements using the particle image velocimetry (PIV) technique were performed to identify the wake structure. Interestingly, the imposed rotation can cause regular vortex shedding in in-line FIV at rotation rates that see suppression of the Bénard–von-Kármán vortex shedding in the case of a rigidly mounted cylinder ($\unicode[STIX]{x1D6FC}\gtrsim 1.75$). There is a monotonic increase in the drag coefficient with rotation rate beyond $\unicode[STIX]{x1D6FC}=2$ for a non-oscillating rotating cylinder. This suggests that the mechanism for sustaining the large rotation-induced galloping oscillations at higher $\unicode[STIX]{x1D6FC}$ is due to a combination of aerodynamic forcing from the locked induced vortex shedding associated with the oscillations, assisted by aerodynamic forcing, evaluated using quasi-steady theory.


Author(s):  
A. Inasawa ◽  
K. Toda ◽  
M. Asai

Disturbance growth in the wake of a circular cylinder moving at a constant acceleration is examined experimentally. The cylinder is installed on a carriage moving in the still air. The results show that the critical Reynolds number for the onset of the global instability leading to a self-sustained wake oscillation increases with the magnitude of acceleration, while the Strouhal number of the growing disturbance at the critical Reynolds number is not strongly dependent on the magnitude of acceleration. It is also found that with increasing the acceleration, the Ka´rma´n vortex street remains two-dimensional even at the Reynolds numbers around 200 where the three-dimensional instability occurs to lead to the vortex dislocation in the case of cylinder moving at constant velocity or in the case of cylinder wake in the steady oncoming flow.


Crystals ◽  
2019 ◽  
Vol 9 (4) ◽  
pp. 217 ◽  
Author(s):  
Jia-Jia Yu ◽  
Lu Zhang ◽  
Ting Shen ◽  
Li Zhang ◽  
You-Rong Li

A series of three-dimensional numerical simulations were performed to understand the thermal-solutal capillary-buoyancy flow of Ge1-xSix melts during Czochralski crystal growth with a rotating crystal or crucible. The crystal and crucible rotation Reynolds numbers in this work are 0∼3.5 × 103 (0∼4.4 rpm) and 0∼−2.4 × 103 (0∼−1.5 rpm), respectively. Simulation results show that if the thermal capillary Reynolds number is relatively low, the flow will be steady and axisymmetric, even though the crystal or crucible rotates at a constant rate. The critical thermal capillary Reynolds number for the initiation of the three-dimensional oscillatory flow is larger than that of pure fluids. As the crystal or crucible rotation rate increases, the critical thermal capillary Reynolds number first increases and then decreases. The dominant flow pattern after the flow destabilization is azimuthal traveling waves. Furthermore, a reversed evolution from the oscillatory spoke pattern to traveling waves appears in the melt. Once the crystal or crucible rotation rate is relatively large, the traveling waves respectively evolve to rotating waves at the crystal rotation and a spindle-like pattern at the crucible rotation. In addition, the maximum amplitude of solute concentration oscillation on the free surface initially decreases, but finally rises with the crystal or crucible rotation rate increasing.


2000 ◽  
Vol 123 (2) ◽  
pp. 347-358 ◽  
Author(s):  
P. Bagchi ◽  
M. Y. Ha ◽  
S. Balachandar

Direct numerical solution for flow and heat transfer past a sphere in a uniform flow is obtained using an accurate and efficient Fourier-Chebyshev spectral collocation method for Reynolds numbers up to 500. We investigate the flow and temperature fields over a range of Reynolds numbers, showing steady and axisymmetric flow when the Reynolds number is less than 210, steady and nonaxisymmetric flow without vortex shedding when the Reynolds number is between 210 and 270, and unsteady three-dimensional flow with vortex shedding when the Reynolds number is above 270. Results from three-dimensional simulation are compared with the corresponding axisymmetric simulations for Re>210 in order to see the effect of unsteadiness and three-dimensionality on heat transfer past a sphere. The local Nusselt number distribution obtained from the 3D simulation shows big differences in the wake region compared with axisymmetric one, when there exists strong vortex shedding in the wake. But the differences in surface-average Nusselt number between axisymmetric and three-dimensional simulations are small owing to the smaller surface area associated with the base region. The shedding process is observed to be dominantly one-sided and as a result axisymmetry of the surface heat transfer is broken even after a time-average. The one-sided shedding also results in a time-averaged mean lift force on the sphere.


Sign in / Sign up

Export Citation Format

Share Document