Vortex shedding and three-dimensional behaviour of flow past a cylinder confined in a channel

A numerical investigation of the flow past a circular cylinder centred in a two-dimensional channel of varying width is presented. For low Reynolds numbers, the flow is steady. For higher Reynolds numbers, vortices begin to shed periodically from the cylinder. In general, the Strouhal frequency of the shedding vortices increases with blockage ratio. In addition, a two-dimensional instability of the periodic vortex shedding is found, both empirically and by means of a Floquet stability analysis. The instability leads to a beating behaviour in the lift and drag coefficients of the cylinder, which occurs at a Reynolds number higher than the critical Reynolds number for the three-dimensional mode A-type instability, but lower than a Reynolds number for any mode B-type instability.     

Some new characteristics of the confined flow over circular cylinders at low reynolds numbers

This study summarises some new characteristics of the fluid flow over a confined circular cylinder at low Reynolds numbers. Results from both two- and three-dimensional direct numerical simulations are presented at blockage ratio between 0.1 and 0.9 and Reynolds number between 120 and 500. Floquet stability analysis of selected cases will also be presented. From the two-dimensional simulations, it is found that the fluctuating lift forces decreases with blockage ratio and becomes zero (where the flow is steady) at blockage ratio of approximately 0.7–0.8. Upon further increasing the blockage ratio to 0.9, the simulations show a dramatic increase in the fluctuating lift forces, nearly an order of magnitude greater than previously reported for an unconfined cylinder flow. It is also found that for blockage ratio of 0.5, there is a long term two-dimensional instability that becomes more prominent with increasing Reynolds number. This instability has a time scale of approximately 105 time units ($D/U_{\mathit{max}}$) at Reynolds number of 500. In addition, the transition between two- and three-dimensional flow at blockage ratios up to 0.5 is investigated. It is shown that the transition Reynolds number decreases with increasing blockage ratio. At high blockage ratio of 0.5, as we increase the Reynolds number, the transition to three-dimensional flow is shown to go from unsteady two-dimensional to steady three-dimensional before transitioning to unsteady three-dimensional flow.     

棱柱绕流流动的数值模拟

利用数值方法对长宽比为1／3, 1和3的棱柱绕流在雷诺数为100的非稳态流动特性进行了分析和研究。采用有限体积法对棱柱绕流的二维流动N-S方程进行离散求解,分析和研究了非稳态的棱柱绕流流场,升力系数,阻力系数和涡动特性,数值模拟的结果与相关文献的数据比较吻合。通过上述研究能够为了解棱柱绕流的非稳态流动特性提供有力的帮助。而对棱柱三维流动的模拟分析和对雷诺数的变化对棱柱流动特性的影响进行研究,将为掌握棱柱绕流的工程特性打下基础。     

Experimental investigation of three-dimensional flow instabilities in a rotating lid-driven cavity

The swirling flow between a rotating lid and a stationary cylinder is studied experimentally. The flow is governed by two parameters: the ratio of container height to disk radius, h, and the Reynolds number, Re, based on the disk angular velocity, cylinder radius and kinematic viscosity of the working liquid. For the first time, the onset of three-dimensional flow behavior is measured by combining the high spatial resolution of particle image velocimetry and the temporal accuracy of laser Doppler anemometry. A detailed mapping of the transition scenario from steady and axisymmetric flow to unsteady and three-dimensional flow is investigated for 1 ≥ h ≥ 3.5. The flow is characterized by the development of azimuthal modes of different wave numbers. A range of different modes is detected and critical Reynolds numbers and associated frequencies are identified. The results are compared to the numerical stability analysis of Gelfgat et al. (J Fluid Mech 438:363–377, 2001). In most cases, the measured onset of three-dimensionality is in good agreement with the numerical results and disagreements can be explained by bifurcations not accounted for by the numerical stability analysis.     

Topology of vortex creation in the cylinder wake

We analyze the topology of the two-dimensional flow around a circular cylinder at moderate Reynolds numbers in the regime where the vortex wake is created. A normal form for the stream function close to the cylinder is presented and used to predict the streamline pattern both in the steady and the periodic regime, where two different vortex shedding scenarios are identified. The theoretical predictions are verified numerically. For the vorticity, a very different topology occurs with infinite nested sequences of iso-curves moving downstream. General equations of motion for critical points are derived.     

Two-dimensional unsteady laminar flow of a power law fluid across a square cylinder

The two-dimensional and unsteady free stream flow of power law fluids past a long square cylinder has been investigated numerically in the range of conditions 60≤Re≤160 and 0.5≤n≤2.0. Over this range of Reynolds numbers, the flow is periodic in time. A semi-explicit finite volume method has been used on a non-uniform collocated grid arrangement to solve the governing equations. The global quantities such as drag coefficients, Strouhal number and the detailed kinematic variables like stream function, vorticity and so on, have been obtained for the above range of conditions. While, over this range of Reynolds number, the flow is known to be periodic in time for Newtonian fluids, a pseudo-periodic flow regime displaying more than one dominant frequency in the lift is observed for shear-thinning fluids. This seems to occur at Reynolds numbers of 120 and 140 for n=0.5 and 0.6, respectively. Broadly speaking, the smaller the value of the power law index, lower is the Reynolds number of the onset of the pseudo-periodic regime. This work is concerned only with the fully periodic regime and, therefore, the range of Reynolds numbers studied varies with the value of the power law index. Not withstanding this aspect, in particular here, the effects of Reynolds number and of the power law index have been elucidated in the unsteady laminar flow regime. The leading edge separation in shear-thinning fluids produces an increase in drag values with the increasing Reynolds number, while shear-thickening fluid behaviour delays this separation and shows the lowering of the drag coefficient with the Reynolds number. Also, the preliminary results suggest the transition from the steady to unsteady flow conditions to occur at lower Reynolds numbers in shear-thinning fluids than that in Newtonian fluids.     

Time resolved analysis of steady and oscillating flow in the upper human airways

In this experimental study a thorough analysis of the steady and unsteady flow field in a realistic transparent silicone lung model of the first bifurcation of the upper human airways will be presented. To determine the temporal evolution of the flow time-resolved particle-image velocimetry recordings were performed for a Womersley number range 3.3 ≤ α ≤ 5.8 and Reynolds numbers of Re _D = 1,050, 1,400, and 2,100. The results evidence a highly three-dimensional and asymmetric character of the velocity field in the upper human airways, in which the influence of the asymmetric geometry of the realistic lung model plays a significant role for the development of the flow field in the respiratory system. At steady inspiration, the flow shows independent of the Reynolds number a large zone with embedded counter-rotating vortices in the left bronchia ensuring a continuous streamwise transport into the lung. At unsteady flow the critical Reynolds number, which describes the onset of vortices in the first bifurcation, is increased at higher Womersley number and decreased at higher Reynolds number. At expiration the unsteady and steady flows are almost alike.     

Motion of a circular cylinder in a viscous liquid between parallel plates

The two-dimensional motion of a cylinder in a viscous fluid between two parallel walls of a vertical channel is studied. It is found that when the cylinder moves very closely along one of the channel walls, it always rotates in the direction opposite to that of contact rolling along the nearest wall. When the cylinder is away from the walls, its rotation depends on the Reynolds number of the flow. In this study two numerical methods were used. One is for the unsteady motion of a sedimenting cylinder initially released from a position close to one of the channel walls, where the Navier-Stokes equations are solved for the fluid and Newton's equations of motion are solved for the rigid cylinder. The other method is for the steady flow in which a cylinder is fixed in a uniform flow field where the channel walls are sliding past the cylinder at the speed of the approaching flow, or equivalently a cylinder is moving with a constant velocity in a quiescent fluid. The flow field, the drag, the side force (lift), and the torque experienced by the cylinder are studied in detail. The effects of the cylinder location in the channel, the size of the channel relative to the cylinder diameter, and the Reynolds number of the flow are examined. In the limit when the cylinder is translating very closely along one of the walls, the flow in the gap between the cylinder and the wall is solved analytically using lubrication theory, and the numerical solution in the other region is used to piece together the whole flow field.This research was supported by NSF DMR91-20668 through the Laboratory for Research on the Structure of Matter at the University of Pennsylvania and from the Research Foundation of the University of Pennsylvania.     

Subcritical spatial transition of swept Hiemenz flow

It is known from experimental investigations that the leading-edge boundary layer of a swept wing exhibits transition to turbulence at subcritical Reynolds numbers, i.e. at Reynolds numbers which lie below the critical Reynolds number predicted by linear stability theory. In the present work, we investigate this subcritical transition process by direct numerical simulations of a swept Hiemenz flow in a spatial setting. The laminar base flow is perturbed upstream by a pair of stationary counter-rotating vortex-like disturbances. This perturbation generates high- and low-speed streaks by a non-modal growth mechanism. Further downstream, these streaky structures exhibit a strong instability to secondary perturbations which leads to a breakdown to turbulence.The observed transition mechanism has strong similarities to by-pass transition mechanisms found for two-dimensional boundary layers. It can be shown that transition strongly depends on the amplitude of the primary perturbation as well as on the frequency of the secondary perturbation.     

A numerical study of steady and unsteady viscoelastic flow past bounded cylinders

We consider two-dimensional, inertia-free, flow of a constant-viscosity viscoelastic fluid obeying the FENE-CR equation past a cylinder placed symmetrically in a channel, with a blockage ratio of 0.5. Through numerical simulations we show that the flow becomes unsteady when the Deborah number (using the usual definition) is greater than De ≈ 1.3, for an extensibility parameter of the model of L² = 144. The transition from steady to unsteady flow is characterised by a small pulsating recirculation zone of size approximately equal to 0.15 cylinder radius attached to the downstream face of the cylinder. There is also a rise in drag coefficient, which shows a sinusoidal variation with time. The results suggest a possible triggering mechanism leading to the steady three-dimensional Gortler-type vortical structures, which have been observed in experiments of the flow of a viscoelastic fluid around cylinders. The results reveal that the reason for failure of the search for steady numerical solutions at relatively high Deborah numbers is that the two-dimensional flow separates and eventually becomes unsteady. For a lower extensibility parameter, L² = 100, a similar recirculation is formed given rise to a small standing eddy behind the cylinder which becomes unsteady and pulsates in time for Deborah numbers larger than De ≈ 4.0–4.5.     

Numerical simulation of Reynolds number effects on velocity shear flow around a circular cylinder

Three-dimensional Direct Numerical Simulation (DNS) and Large Eddy Simulation (LES) are performed to investigate the shear effects on flow around a circular cylinder at Reynolds numbers of Re=60–1000. The shear parameter, β, which is based on the velocity gradient, cylinder diameter and upstream mean velocity at the center plane of the cylinder, varies from 0 to 0.30. Variations of Strouhal number, drag and lift coefficients, and unsteady wake structures with shear parameter are studied, along with their dependence on Reynolds number. The presented simulation provides detailed information for the flow field around a circular cylinder in shear flow. This study shows that the Strouhal number exhibits no significant variation with shear parameter. The stagnation point moves to the high-velocity side almost linearly with shear parameter, and this result mainly influences the aerodynamic forces acting on a circular cylinder in shear flow. Both the Reynolds number and shear parameter influence the movement of the stagnation point and separation point. Mode A wake instability is suppressed into parallel vortex shedding mode at a certain shear parameter. The lift force increases with increasing shear parameter and acts from the high-velocity side to the low-velocity side. In addition, a simple method to estimate the lift force coefficient in shear flow is provided.     

A review of rotating cylinder wake transitions

Recent work on the flow past a rotating cylinder is reviewed and further investigated at low Reynolds numbers. The various two- and three-dimensional transitions that occur as the rotation rate is increased are detailed. Two steady states, steady state I and steady state II, are identified based on the physical characteristics of the wake and the drag force on the body. Steady state I occurs at lower rotation rates, while state steady state II occurs at higher rotation rates. Linear stability analysis shows that two three-dimensional modes become unstable on steady state I and steady state II. Floquet stability analysis of the unsteady base flows that occur at very low rotation rates shows the presence of five three-dimensional modes. The curves of marginal stability are presented, followed by a comparison of numerical simulations to their experimentally obtained counterparts. Furthermore, the spatio-temporal characteristics of each mode and the likely underlying physical mechanisms are briefly discussed.     

Unsteady characteristics of the cavity flow around a rough circular cylinder

Effects of surface roughness on the unsteady cavitating flow around a two-dimensional circular cylinder were experimentally investigated at Reynolds numbers from 1.36 × 10⁵ to 1.78 × 10⁵. Two patterns of surface roughness were investigated, a double-cut pattern and a single-cut pattern. The cavity elongates with an increase of the surface roughness, especially in supercavitating flow. However, for some roughness parameters tested, the cavity length exhibits an extreme decrease. In a particular case of the double-cut pattern, there exists the minimum cavity behind the cylinder.     

Anisotropic Organised Eddy Simulation for the prediction of non-equilibrium turbulent flows around bodies

The unsteady turbulent flow around bodies at high Reynolds number is predicted by an anisotropic eddy-viscosity model in the context of the Organised Eddy Simulation (OES). A tensorial eddy-viscosity concept is developed to reinforce turbulent stress anisotropy, that is a crucial characteristic of non-equilibrium turbulence in the near-region. The theoretical aspects of the modelling are investigated by means of a phase-averaged PIV in the flow around a circular cylinder at Reynolds number 1.4×10⁵. A pronounced stress–strain misalignment is quantified in the near-wake region of the detached flow, that is well captured by a tensorial eddy-viscosity concept. This is achieved by modelling the turbulence stress anisotropy tensor by its projection onto the principal matrices of the strain-rate tensor. Additional transport equations for the projection coefficients are derived from a second-order moment closure scheme. The modification of the turbulence length scale yielded by OES is used in the Detached Eddy Simulation hybrid approach. The detached turbulent flows around a NACA0012 airfoil (2-D) and a circular cylinder (3-D) are studied at Reynolds numbers 10⁵ and 1.4×10⁵, respectively. The results compared to experimental ones emphasise the predictive capabilities of the OES approach concerning the flow physics capture for turbulent unsteady flows around bodies at high Reynolds numbers.     

Effect of the pattern of initial perturbations on the steady pipe flow regime

The influence of the inlet flow formation mode on the steady flow regime in a circular pipe has been investigated experimentally. For a given inlet flow formation mode the Reynolds number Re* at which the transition from laminar to turbulent steady flow occurred was determined. With decrease in the Reynolds number the difference between the resistance coefficients for laminar and turbulent flows decreases. At a Reynolds number approximately equal to 1000 the resistance coefficients calculated from the Hagen-Poiseuille formula for laminar steady flow and from the Prandtl formula for turbulent steady flow are equal. Therefore, we may assume that at Re > 1000 steady pipe flow can only be laminar and in this case it is meaningless to speak of a transition from one steady pipe flow regime to the other. The previously published results [1–9] show that the Reynolds number at which laminar goes over into turbulent steady flow decreases with increase in the intensity of the inlet pulsations. However, at the highest inlet pulsation intensities realized experimentally, turbulent flow was observed only at Reynolds numbers higher than a certain value, which in different experiments varied over the range 1900–2320 [10]. In spite of this scatter, it has been assumed that in the experiments a so-called lower critical Reynolds number was determined, such that at higher Reynolds numbers turbulent flow can be observed and at lower Reynolds numbers for any inlet perturbations only steady laminar flow can be realized. In contrast to the lower critical Reynolds number, the Re* values obtained in the present study, were determined for given (not arbitrary) inlet flow formation modes. In this study, it is experimentally shown that the Re* values depend not only on the pipe inlet pulsation intensity but also on the pulsation flow pattern. This result suggests that in the previous experiments the Re* values were determined and that their scatter is related with the different pulsation flow patterns at the pipe inlet. The experimental data so far obtained are insufficient either to determine the lower critical Reynolds number or even to assert that this number exists for a pipe at all.     

Turbulence properties in the cylinder wake at high Reynolds numbers

The present contribution analyses the turbulence properties in unsteady flows around bluff body wakes and provides a database for improvement and validation of turbulence models, concerning the present class of nonequilibrium flows. The flow around a circular cylinder with a low aspect ratio and a high blockage coefficient is investigated. This confined environment is used in order to allow direct comparisons with realisable 3-D Navier–Stokes computations avoiding ‘infinite’ conditions. The flow is investigated in the beginning of the critical regime at Reynolds number 140 000. The analysis is carried out by means of 2-D PIV, of 3-C PIV and of high-frequency 2-D PIV. The experimental analysis contributes to confirm the validity of advanced statistical turbulence modelling for unsteady flows around bodies.     

A thermal LBM-LES model in body-fitted coordinates: Flow and heat transfer around a circular cylinder in a wide Reynolds number range

In order to apply the lattice Boltzmann method (LBM) for modeling passive heat transfer at high Reynolds numbers, a number of models were proposed by introducing the large eddy simulation (LES) into the LBM framework to improve numerical stability. Our study finds that the generalized form of interpolation-supplemented LBM (GILBM), likewise, can locally modify the dimensionless relaxation time, thus enhancing the numerical stability at high Reynolds numbers. Given additional advantages of the GILBM in dealing with complicated geometries and improving computing accuracy, a thermal LBM-LES model in body-fitted coordinates is established in this paper. Numerical validation is performed by investigating the flow and heat transfer around a circular cylinder over a wide range of Reynolds numbers. The obtained results agree well with both experimental and numerical data in the previous work. Meanwhile, the effects of Reynolds number and Prandtl number on thermodynamic features of flow past a circular cylinder are revealed. It is found out that when the Reynolds number exceeds the critical value, the local Nusselt number fluctuates rapidly in a specific region of the rear cylinder surface affected by the Prandtl number. In the near-wake region, the temperature field exhibits significant dependence on the Prandtl number at moderate Reynolds numbers, while such effects turn to be slight at high Reynolds numbers.     

Numerical study of a separated-reattaching flow

The flow over a normal plate with a long, downstream splitter plate is numerically studied with fine spatial and temporal resolutions. The time-dependent two-dimensional Navier-Stokes equations are integrated in time using a high-order, upwind-biased finite-difference scheme. Calculations are performed at several Reynolds numbers to study the evolution of the unsteady nature of the flow. Steady separated flow has been observed for Reynolds numbers up to 150, after which unsteady vortical structures are seen to develop from the shear layer. The time-mean flow characteristics in the steady and unsteady regimes are described. The calculations are seen to agree fairly well with experimental data at high Reynolds numbers.This work was supported by the Office of Naval Research under Grant N-00014-92-J-1640.     

Linear stability analysis of flow in a periodically grooved channel

We have conducted the linear stability analysis of flow in a channel with periodically grooved parts by using the spectral element method. The channel is composed of parallel plates with rectangular grooves on one side in a streamwise direction. The flow field is assumed to be two‐dimensional and fully developed. At a relatively small Reynolds number, the flow is in a steady‐state, whereas a self‐sustained oscillatory flow occurs at a critical Reynolds number as a result of Hopf bifurcation due to an oscillatory instability mode. In order to evaluate the critical Reynolds number, the linear stability theory is applied to the complex laminar flow in the periodically grooved channel by constituting the generalized eigenvalue problem of matrix form using a penalty‐function method. The critical Reynolds number can be determined by the sign of a linear growth rate of the eigenvalues. It is found that the bifurcation occurs due to the oscillatory instability mode which has a period two times as long as the channel period. Copyright © 2003 John Wiley & Sons, Ltd.