Vortex shedding from a square cylinder in presence of a moving wall

A numerical study on the flow past a square cylinder placed parallel to a wall, which is moving at the speed of the far field has been made. Flow has been investigated in the laminar Reynolds number (based on the cylinder length) range. We have studied the flow field for different values of the cylinder to wall separation length. The governing unsteady Navier–Stokes equations are discretized through the finite volume method on a staggered grid system. A SIMPLE type of algorithm has been used to compute the discretized equations iteratively. A shear layer of negative vortex generates along the surface of the wall, which influences the vortex shedding behind the cylinder. The flow‐field is distinct from the flow in presence of a stationary wall. An alternate vortex shedding occurs for all values of gap height in the unsteady regime of the flow. The strong positive vortex pushes the negative vortex upwards in the wake. The gap flow in the undersurface of the cylinder is strong and the velocity profile overshoots. The cylinder experiences a downward force for certain values of the Reynolds number and gap height. The drag and lift are higher at lower values of the Reynolds number. Copyright © 2005 John Wiley & Sons, Ltd.     

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.     

同轴圆柱间滑移旋转流动问题的同伦解析解

从理论上研究了具有延伸柱面的同轴圆柱间滑移流动问题.通过引入适当的相似变换将控制方程组转化为一类非线性边值问题, 利用同伦分析方法首次获得问题的近似解析解, 分析讨论滑移参数、雷诺数和内外筒半径比对流动的影响.结果表明: 滑移边界参数对剪切应力有较大的影响, 滑移边界小的流体对壁面和边界层内流场施加了更大的剪切力; 增大雷诺数Re, 能增大内筒的壁面剪切力和扭矩; 内外筒半径比对流动结构也有较大的影响, 内筒半径固定, 增大外筒半径能够减少内筒的纵向剪切力.      

Early stages of an impulsively started unsteady laminar flow past tapered trapezoidal cylinders

Characteristics of the developing recirculation region behind a tapered trapezoidal cylinder and its interaction with the separating shear layer from the leading edges were studied numerically for an impulsively started laminar flow. An unsteady stream function–vorticity formulation was used. The Reynolds numbers considered range from 25 to 1000. Pressure contours, surface pressure coefficient, wake length and drag coefficient were studied through the streamline flow field. Main flow and subflow regimes were identified by an analysis of the evolution of the flow characteristics. It was found that typically, for a given trapezoidal cylinder, flow starts with no separation. As time advances, the symmetrical standing zone of recirculation develops aft of the trapezoidal cylinder. The rate of growth in width, length and structure of the aft end eddies depends on the Reynolds number. In time, separated flow from the leading edges of the trapezoidal cylinder also develops and forms growing separation bubbles on the upper and lower inclined surfaces of the trapezoidal cylinder. As time advances, the separation bubbles on the upper and lower inclined surfaces of the cylinder grow towards the downstream regions and eventually merge with the swelling symmetrical eddies aft of the cylinder. This merging of the flows creates a complex flow regime with a disturbed tertiary flow zone near the merging junction. Eventually, depending on the Reynolds number and the tapered angle of the trapezoidal cylinder, the flow develops into a specific category of symmetrical standing recirculatory flow with its own distinct characteristics. Comparisons with the available results of other investigators showed very good agreement. © 1998 John Wiley & Sons, Ltd.     

Influence of the orientation of a square cylinder on the wake properties

An experimental investigation of flow around a square cylinder placed at various angles with respect to the approach fluid velocity is reported. The focus of the study is toward examining the sensitivity of the wake properties to the cylinder orientation and Reynolds number. Angles of incidence in the range of 0-60° and Reynolds numbers of 1340, 4990, and 9980 have been considered. Velocity measurements have been carried out using an X-wire hotwire anemometer. The Strouhal number and the drag coefficient of the cylinder have been computed from the wake measurements. Utilizing the velocity traces at distinct probe locations in the near and the far wake, statistical properties such as the RMS velocities and the spectra have been obtained. Results obtained in the present work revealed that for a cylinder with zero inclination, flow separates from the corners on the face exposed to the incoming flow. For inclinations greater than zero, the points of separation on the cylinder move downstream and the wake size increases, but the separated shear layer rolls up over a shorter distance. These factors lead to a reduced drag coefficient and a higher Strouhal number. The center-line recovery of the time-averaged velocity and the decay rates of velocity fluctuations depend on the Reynolds number. A marginal effect of the cylinder orientation is also seen.     

Numerical simulation and experimental visualization of the influence of the deformation frequency of a radially deforming circular cylinder impulsively started on cylinder wake

A periodic superimposed motion may notably influence the flow structure and the development of the convective heat transfer relative to non‐deformable case. In particular, a radial deformation of a circular cylinder, may cause a possible synchronization with the cylinder wake, which is itself periodic when Vortex Street takes place. This synchronization phenomenon, often called ‘lock‐in’, may cause undesirable effects, but may also constitute a way of controlling the wake development. Body deformability may be used as wake control device that would favourably affect the interplay of primary and secondary vorticities, thus reducing the drag coefficient. These numerical and experimental studies are done herein for a Reynolds number equal to 23500. The problem is resolved by using the Navier–Stokes equations in the vorticity‐stream function form. The vorticity transport equation is solved by a second‐order finite difference method in both directions of the domains. The Poisson equation for the stream‐function is solved by a SOR method. The advance in time is achieved by a second‐order Adams–Bashforth scheme. The effect of turbulence is represented by eddy viscosity ν_t, which is determined by a sub‐grid‐scale model. In the present study, we use a Smagorinsky model. Copyright © 2003 John Wiley & Sons, Ltd.     

Wall proximity effects on the effectiveness of upstream control rod

Two dimensional flow over a circular cylinder with an upstream control rod of same diameter is simulated in unbound condition and in wall bounded conditions. The cylinders are placed at various heights from the wall and the inter-distance between cylinders is also varied. The control rod is subjected to different rotation rates. It is found that, in unbound condition, rotating the control rod decreases the critical pitch length (S/D_cr) and increases the drag and Strouhal number of the main cylinder. In presence of plane wall, the shielding provided by the separated shear layers from the control rod in cavity regime is deteriorated due to deflection of shear layers which results in higher drag and large fluctuation of lift coefficient. However, in wake impingement regime, the binary vortices from the control rod are weakened due to diffusion of vorticity and hence, the main cylinder experiences a lower drag and small lift fluctuations than that of unbound condition. The critical height of vortex suppression (H/D_cr) is higher in cavity regime than that of wake impingement regime due to the single extended-bluff body like configuration. The rotation of control rod energizes the wall boundary layer and increases the critical height of vortex suppression. Increasing the rotational rate of control rod decreases the drag force and reduces the amplitude of lift fluctuation. Analysis of the wall shear stress distribution reveals that it suffers a sudden drop at moderate height where the normal Karman vortex shedding changes to irregular shedding consisting of single row of negative vortices. Modal structures obtained from dynamic mode decomposition (DMD) reveal that the flow structures behind the main cylinder are suppressed due to wall and the flow is dominated by the wake of control rod.     

Orthogonal wavelet analysis of turbulent wakes behind various bluff bodies

The multi-scale structures of turbulent wakes generated by three kinds of bluff body, i.e. circular cylinder, square cylinder and compound of cylinder and square (CS) cylinders, have been experimentally investigated in this paper. Firstly, the instantaneous velocity fields and vorticity were measured by the high-speed PIV technique in a circulating water channel. The instantaneous streamlines and corresponding normalized vorticity contours are obtained at a Reynolds number of 5600. Then one- and two-dimensional wavelet multi-resolution technique was used to analyze the instantaneous velocities and vorticity measured by the high-speed PIV. The turbulence structures were separated into a number of subsets based on their central frequencies, which are linked with the turbulence scales. The instantaneous vorticity and Reynolds shear stresses of various scales were examined and compared between the three generators. It is found that the large-scale turbulent structure makes the largest contribution to the vorticity and Reynolds shear stresses for the three wake generators and exhibits a strong dependence upon the initial conditions or the wake generators. The large-scale vorticity and the sizes of vortex in the circular and square cylinders are larger than those in the CS cylinder wake. The contributions to the Reynolds shear stresses from the large-scale turbulent structures account for 90-96% to the measured maximum Reynolds shear stresses for the three wakes. However, the small-scale structures make less contribution to the vorticity and Reynolds shear stresses.     

Numerical study of vortex shedding from a rotating cylinder immersed in a uniform flow field

A numerical study is made of the unsteady two‐dimensional, incompressible flow past an impulsively started translating and rotating circular cylinder. The Reynolds number (Re) and the rotating‐to‐translating speed ratio (α) are two controlled parameters, and the influence of their different combinations on vortex shedding from the cylinder is investigated by the numerical scheme sketched below. Associated with the streamfunction (ψ)–vorticity (ω) formulation of the Navier–Stokes equations, the Poisson equation for ψ is solved by a Fourier/finite‐analytic, separation of variable approach. This approach allows one to attenuate the artificial far‐field boundary, and also yields a global conditioning on the wall vorticity in response to the no‐slip condition. As for the vorticity transport equation, spatial discretization is done by means of finite difference in which the convection terms are handled with the aid of an ENO (essentially non‐oscillatory)‐like data reconstruction process. Finally, the interior vorticity is updated by an explicit, second‐order Runge–Kutta method. Present computations fall into two categories. One with Re=10³ and α≤3; the other with Re=10⁴ and α≤2. Comparisons with other numerical or physical experiments are included. Copyright © 2000 John Wiley & Sons, Ltd.     

Suppression of vortex shedding for flow around a circular cylinder using optimal control

Adjoint formulation is employed for the optimal control of flow around a rotating cylinder, governed by the unsteady Navier–Stokes equations. The main objective consists of suppressing Karman vortex shedding in the wake of the cylinder by controlling the angular velocity of the rotating body, which can be constant in time or time‐dependent. Since the numerical control problem is ill‐posed, regularization is employed. An empirical logarithmic law relating the regularization coefficient to the Reynolds number was derived for 60?Re?140. Optimal values of the angular velocity of the cylinder are obtained for Reynolds numbers ranging from Re=60 to Re=1000. The results obtained by the computational optimal control method agree with previously obtained experimental and numerical observations. A significant reduction of the amplitude of the variation of the drag coefficient is obtained for the optimized values of the rotation rate. Copyright © 2002 John Wiley & Sons, Ltd.     

Flow past a square cylinder at low Reynolds numbers

Results are presented for the flow past a stationary square cylinder at zero incidence for Reynolds number, Re ? 150. A stabilized finite‐element formulation is employed to discretize the equations of incompressible fluid flow in two‐dimensions. For the first time, values of the laminar separation Reynolds number, Re_s, and separation angle, θ_s, at Re_s are predicted. Also, the variation of θ_s with Re is presented. It is found that the steady separation initiates at Re = 1.15. Contrary to the popular belief that separation originates at the rear sharp corners, it is found to originate from the base point, i.e. θ_s=180^° at Re = Re_s. For Re > 5, θ_s approaches the limit of 135 ^°. The length of the separation bubble increases approximately linearly with increasing Re. The drag coefficient varies as Re^?0.66. Flow characteristics at Re ? 40 are also presented for elliptical cylinders of aspect ratios 0.2, 0.5, 0.8 and 1 (circle) having the same characteristic dimension as the square and major axis oriented normal to the free‐stream. Compared with a circular cylinder, the flow separates at a much lower Re from a square cylinder leading to the formation of a bigger wake (larger bubble length and width). Consequently, at a given Re, the drag on a square cylinder is more than the drag of a circular cylinder. This suggests that a cylinder with square section is more bluff than the one with circular section. Among all the cylinder shapes studied, the square cylinder with sharp corners generates the largest amount of drag. Copyright © 2010 John Wiley & Sons, Ltd.     

Reduction in drag and vortex shedding frequency through porous sheath around a circular cylinder

A numerical study on the laminar vortex shedding and wake flow due to a porous‐wrapped solid circular cylinder has been made in this paper. The cylinder is horizontally placed, and is subjected to a uniform cross flow. The aim is to control the vortex shedding and drag force through a thin porous wrapper around a solid cylinder. The flow field is investigated for a wide range of Reynolds number in the laminar regime. The flow in the porous zone is governed by the Darcy–Brinkman–Forchheimer extended model and the Navier–Stokes equations in the fluid region. A control volume approach is adopted for computation of the governing equations along with a second‐order upwind scheme, which is used to discretize the convective terms inside the fluid region. The inclusion of a thin porous wrapper produces a significant reduction in drag and damps the oscillation compared with a solid cylinder. Dependence of Strouhal number and drag coefficient on porous layer thickness at different Reynolds number is analyzed. The dependence of Strouhal number and drag on the permeability of the medium is also examined. Copyright © 2009 John Wiley & Sons, Ltd.     

Calculation of oseen flows past a circular cylinder at low reynolds numbers

The velocity, pressure, vorticity and streamfunction are computed in the Oseen hydrodynamic field of an unbounded fluid past a circular cylinder in the Reynolds Number range going from 0.4 to 12. The boundary condition is satisfied by means of the method of least squares that determines suitable coefficients for Faxén series. Particular investigation is made of the wake region in which calculations are made of flow patterns, velocity and vorticity distributions. It is shown that, attached vortices arise at the rear of the cylinder at Reynolds Number Re=3.025. Calculated drag coefficients are in good agreement with known results of the works of several authors up to a Reynolds Number of 20.     

圆角化对方柱气动性能影响的流场机理

方形截面柱体的圆角化处理是常用的流动控制方法,但其流场作用机理尚未被澄清.采用大涡模拟方法,在雷诺数为2.2$\times$10$^{4}$时,考虑风攻角的影响,对均匀流作用下的标准方柱和圆角方柱的气动性能和流场特性进行了研究,定量分析了圆角化气动措施和风攻角变化对分离泡特性的影响规律,从流场角度澄清了圆角化气动措施对方柱气动性能的影响机理.研究表明:与标准方柱相比,圆角方柱的表面风压、气动力和涡脱强度呈整体下降的趋势,但圆角方柱的斯特劳哈尔数更高;圆角方柱的"分离泡流态'发生在更小的风攻角范围内,分离泡的出现会进一步造成方柱的尾流变窄,涡脱强度减弱;随着风攻角的增大,分离泡的长度会逐渐减小直至消失,分离泡的中心会逐渐向方柱前角(迎风向)和方柱壁面移动;与标准方柱相比,圆角方柱的气流发生初次分离的位置向下游移动,分离后的剪切层更贴近方柱,因而更易发生再附现象;方柱尾流宽度的减小和涡脱强度的减弱是导致圆角方柱气动力减小和斯特劳哈尔数增大的主要原因.      

Direct simulation of multiphase flows with modeling of dynamic interface contact angle

We describe a modeling technique for dynamic contact angle between a phase interface and a solid wall using a generalized Navier boundary condition in the context of a front-tracking-based multiphase method. The contact line motion is determined by the generalized Navier slip boundary condition in order to eliminate the infinite shear stress at the contact line. Applying this slip boundary condition only to the interface movement with various slip ratios shows good agreement with experimental results compared to allowing full fluid slip along the solid surface. The interface slip model performs well on grid convergence tests using both the slip ratio and slip length models. A detailed energy analysis was performed to identify changes in kinetic, surface, and potential energies as well as viscous and contact line dissipation with time. A friction coefficient for contact line dissipation was obtained based on the other computed energy terms. Each energy term and the friction coefficient were compared for different grid resolutions. The effect of varying the slip ratio as well as the contact angle distribution versus contact line speed was analyzed. The behavior of drop impact on a solid wall with different advancing and receding angles was investigated. Finally, the proposed dynamic contact model was extended to three dimensions for large-scale parallel calculations. The impact of a droplet on a solid cylinder was simulated to demonstrate the capabilities of the proposing formulation on general solid structures. Widely different contact angles were tested and showed distinctive characteristic behavior clearly.     

Vortex shedding and aerodynamic forces on a circular cylinder in linear shear flow at subcritical Reynolds number

Vortex shedding and aerodynamic forces on a circular cylinder in a linear shear flow with its axis normal to the plane of the velocity shear profile at subcritical Reynolds number are investigated experimentally. 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.27. The Strouhal number has no significant variation with the shear parameter. The time-mean base pressure increases and the fluctuating component of the base pressure decreases significantly with increasing shear parameter. Vortex shedding is suppressed by the velocity shear. Dislocation of the stagnation point takes place and this influences the pressure distribution around the cylinder together with the velocity shear. A mean lift force arises in the shear flow due to asymmetry of the pressure distribution, and it acts from the high velocity side to the low velocity side. In addition, the lift coefficient increases and the drag coefficient decreases with increasing shear parameter.     

大雷诺数线性剪切流绕圆柱的分离流动

本文采用离散涡方法计算了大雷诺数下线性剪切流绕圆柱的分离流动,计算结果给出阻力系数、横向力系数、边界层分离点、斯特罗哈数及尾流中旋涡的分布图案,计算数据与已有的实验数据相吻合。     

大雷诺数线性剪切流绕圆柱的分离流动

本文采用离散涡方法计算了大雷诺数下线性剪切流绕圆柱的分离流动,计算结果给出阻力系数、横向力系数、边界层分离点、斯特罗哈数及尾流中旋涡的分布图案,计算数据与已有的实验数据相吻合。     

Quadratic programming algorithm for wall slip and free boundary pressure condition

Wall slip is often observed in a highly sheared fluid film in a solid gap. This makes a difficulty in mathematical analysis for the hydrodynamic effect because fluid velocity at the liquid–solid interfaces is not known a priori. If the gap has a convergent–divergent wedge, a free boundary pressure condition, i.e. Reynolds pressure boundary condition, is usually used in the outlet zone in numerical solution. This paper, based on finite element method and parametric quadratic programming technique, gives a numerical solution technique for a coupled boundary non‐linearity of wall slip and free boundary pressure condition. It is found that the numerical error decreases with the number of elements in a negative power law having an index larger than 2. Our method does not need an iterative process and can simultaneously gives rise to fluid film pressure distribution, wall slip velocity and surface shear stress. Wall slip always decreases the hydrodynamic pressure. Large wall slip even causes a null hydrodynamic pressure in a pure sliding solid gap. Copyright © 2005 John Wiley & Sons, Ltd.     

Flow around a porous cylinder subject to continuous suction or blowing

In the present experimental investigation the surface pressure distribution, vortex shedding frequency, and the wake flow behind a porous circular cylinder are studied when continuous suction or blowing is applied through the cylinder walls. It is found that even moderate levels of suction/blowing (5% of the oncoming streamwise velocity) have a large impact on the flow around the cylinder. Suction delays separation contributing to a narrower wake width, and a corresponding reduction of drag, whereas blowing shows the opposite behaviour. Both uniform suction and blowing display unexpected flow features which are analysed in detail. Suction shows a decrease of the turbulence intensity throughout the whole wake when compared with the natural case, whilst blowing only shows an effect up to five diameters downstream of the cylinder. The drag on the cylinder is shown to increase linearly with the blowing rate, whereas for suction there is a drastic decrease at a specific suction rate. This is shown to be an effect of the separation point moving towards the rear part of the cylinder, similar to what happens when transition to turbulence occurs in the boundary layer on a solid cylinder. The suction/blowing rate can empirically be represented by an effective Reynolds number for the solid cylinder, and an analytical expression for this Reynolds number representation is proposed and verified. Flow visualizations expose the complexity of the flow field in the near wake of the cylinder, and image averaging enables the retrieval of quantitative information, such as the vortex formation length.