The effect of slip distribution on flow past a circular cylinder

A slip boundary has been shown to have a significant impact on flow past bluff bodies. In this work and using a circular cylinder as a model system, the effects of various slip configurations on the passing flow are investigated. A theoretical analysis using matched-asymptotic expansion is first performed in the small-Reynolds number regime following Stokes and Oseen. A slip boundary condition is shown to lead to only higher-order effects (~1/ln(Re)) on the resulting drag coefficient. For higher Reynolds numbers (100–500), the effects of five types of symmetric slip boundary conditions, namely, no slip, fore-side slip, aft-side slip, flank slip, and all slip on the flow field and pertinent parameters are investigated with numerical simulations. Detailed results on the flow structure and force distribution are presented. Flank slip is found to have the best effect for drag reduction with comparable coverage of slip area. For asymmetric slip distributions, torque and lift are found to generally occur.     

Creeping motion of a slip spherical particle in a circular cylindrical pore

A combined analytical–numerical study for the creeping flow caused by a spherical fluid or solid particle with a slip-flow surface translating in a viscous fluid along the centerline of a circular cylindrical pore is presented. To solve the axisymmetric Stokes equations for the fluid velocity field, a general solution is constructed from the superposition of the fundamental solutions in both cylindrical and spherical coordinate systems. The boundary conditions are enforced first at the pore wall by the Fourier transforms and then on the particle surface by a collocation technique. Numerical results for the hydrodynamic drag force acting on the particle are obtained with good convergence for various values of the relative viscosity or slip coefficient of the particle, the slip parameter of the pore wall, and the ratio of radii of the particle and pore. For the motion of a fluid sphere along the axis of a cylindrical pore, our drag results are in good agreement with the available solutions in the literature. As expected, the boundary-corrected drag force for all cases is a monotonic increasing function of the ratio of particle-to-pore radii, and approaches infinity in the limit. Except for the case that the cylindrical pore is hardly slip and the value of the ratio of particle-to-pore radii is close to unity, the drag force exerted on the particle increases monotonically with an increase in its relative viscosity or with a decrease in its slip coefficient for a constant ratio of radii. In a comparison for the pore shape effect on the axial translation of a slip sphere, it is found that the particle in a circular cylindrical pore in general acquires a lower hydrodynamic drag than in a spherical cavity, but this trend can be reversed for the case of highly slippery particles and pore walls.     

Very slow flow of Bingham viscoplastic fluid around a circular cylinder

Numerical simulations have been used to study the flow of a Bingham viscoplastic fluid around a circular cylinder in an infinite medium with negligible inertia effects. Papanastasiou's regularisation technique has been adopted to approximate the model. The case corresponding to preponderant plasticity effects has been particularly studied and convergence of the solutions examined in detail. The flow kinematics and stresses have been determined. The rigid zones have been identified and characterised. At large Oldroyd numbers, when plasticity effects become preponderant, a viscoplastic boundary layer appears around the cylinder. The characteristics of this viscoplastic boundary layer are quantified. The results are compared with existing theoretical results, concerning particularly the predictions of the viscoplastic boundary layer theory and the plasticity theory.     

Slow motion of a slip spheroid along its axis of revolution

A combined analytical and numerical study of the Stokes flow caused by a rigid spheroidal particle translating along its axis of revolution in a viscous fluid is presented. The fluid is allowed to slip at the surface of the particle. The general solution for the stream function in prolate and oblate spheroidal coordinates can be expressed in an infinite-series form of semi-separation of variables. The slip boundary condition incorporating the shear stress at the particle surface is applied to this general solution to determine its unknown coefficients of the leading orders. The solution of these coefficients can be either numerical results obtained from a boundary-collocation method or explicit formulas derived analytically. The drag force exerted on the spheroidal particle by the fluid is evaluated with good convergence behavior for various values of the slip parameter and aspect ratio of the particle. The agreement between our hydrodynamic drag results and the relevant numerical solutions obtained previously using a singularity method is excellent. Although the drag force acting on the translating spheroid normalized by that on a corresponding sphere with equal equatorial radius increases monotonically with an increase in the axial-to-radial aspect ratio for a no-slip spheroid, it decreases monotonically as this aspect ratio increases for a perfect-slip spheroid. The normalized drag force exerted on a spheroid with a given surface slip coefficient in between the no-slip and perfect-slip limits is not a monotonic function of its aspect ratio. For a spheroid with a fixed aspect ratio, its drag force is a monotonically decreasing function of the slip coefficient of the particle.     

Numerical solutions for impulsively started and decelerated viscous flow past a circular cylinder

A finite difference study of the unsteady two-dimensional flow past a circular cylinder has been conducted using vorticity and streamfunction as the dependent variables. The two cases considered were impulsively started and decelerated flows. The impulsively started problem was considered to validate the method and has yielded results which agree quite closely with existing results from both calculations and experiments. The decelerated flow analysis produced results which can be explained in terms of induced velocity effects from existing wake vortices for both suddenly stopped and uniformly decelerated flows.     

Numerical investigation of steady suction control of flow around a circular cylinder

This paper numerically investigates the effectiveness of the control of steady suction on a stationary circular cylinder with several isolated suction holes on the surface at a subcritical Reynolds number. The control effectiveness as a function of the azimuthal position, spanwise spacing and suction flow rate of the suction holes on the control of the aerodynamic forces on the cylinder and the suppression of alternate vortex shedding are taken into account. The study of the azimuthal location of the suction holes indicates that azimuthal angles of θ=90° and 270°, which are close to the separation point, provide the most substantial decreases in the aerodynamic forces. When restricted to the most effective azimuthal angle, a remarkable control effectiveness can be achieved when the axial spacing between two neighboring suction holes is less than a minimum value even under a small suction momentum coefficient. However, if the axial spacing exceeds the minimum spacing, the control effectiveness will not be saturated even under a very large suction momentum coefficient. Thus, the cause of the effective aerodynamic force control is suggested to be a result of obvious three-dimensional phenomenon in the near wake, which is characterized by the generation of a convergent flow between two neighboring suction hole sections and a stronger, larger three-dimensional vortex pair adjacent to the convergent flow. It has been suggested that this strongly three-dimensional flow pattern is induced by the strong interaction between two neighboring but counter-rotating three-dimensional vortices separately produced by two neighboring suction holes. Moreover, the effects of such three-dimensional flow patterns are investigated in detail based on variations in the flow field and sectional aerodynamic forces in different cross sections. Finally, the upper limit of the axial spacing between two neighboring suction holes to form such a three-dimensional flow pattern is suggested to be between 0.75 D and 1.5 D when the suction flow rate exceeds a certain value.     

Slow motions of a solid spherical particle close to a viscous interface

In order to investigate the hydrodynamic interaction between an interface and a spherical particle and its dependence on the type of interface, it is essential to compute the drag and torque exerted on the sphere in the vicinity of the interface. In this paper, the problem of all slow elementary motions (relative translation and rotation) and stationary movement of a spherical particle next to a solid, viscous or free interface is considered. For low capillary numbers and different values of surface dilatational and shear viscosities in a curvilinear co-ordinate system of revolution with bicylindrical co-ordinates in meridian planes, the problem reduces from three to two dimensions. The model equations and boundary conditions, which contain second-order derivatives of the velocities, transform to an equivalent well-defined system of second-order partial differential equations which is solved numerically for medium and small values of the dimensionless distance to the interface. Very good agreement with the asymptotic equation for a translating sphere close to a solid interface could be achieved. The numerical results reveal in all cases the strong influence of the surface viscosity on the motion of the solid sphere. For small distances from the interface, the drag and torque coefficients change significantly depending on the surface viscosity.     

Vortex shedding around a near-wall circular cylinder induced by a solitary wave

This study developed a two-dimensional generalized vortex method to analyze the shedding of vortices and the hydrodynamic forces resulting from a solitary wave passing over a submerged circular cylinder placed near a flat seabed. Numerical results for validation are compared with other numerical and experimental results, and satisfactory agreement is found. A series of simulations were performed to study the effects of gap-to-diameter ratio and incident wave height on vorticity pattern as well as the forces exerted on the cylinder. The range of the heights of incident waves is from 0.3h to 0.7h, where h is the still water depth. The range of the gap-to-diameter ratios is from 0.1 to 0.8. The results indicate that the flow pattern and the pressure distribution change significantly because of the close proximity of the seabed where the vorticity flux on the seabed-side surface of the cylinder is suppressed. Placing the cylinder nearer the seabed increases the drag and the positive lift on the cylinder. When the gap-to-diameter ratio increases, the pattern of vortices changes because of the interaction between the main recirculation zone and the shear layers separated from the gap. The maxima of drag, lift and total force increase linearly with the height of the incident wave.     

Mixed convection boundary layer flow of a viscoelastic fluid over a horizontal circular cylinder

The steady mixed convection boundary layer flow of a viscoelastic fluid over a horizontal circular cylinder in a stream flowing vertically upwards is numerically studied for both cases of heated and cooled cylinders. The governing partial differential equations are transformed into dimensionless forms using an appropriate transformation and then solved numerically using the Keller-box method. The comparison between the solutions obtained and those for a Newtonian fluid is found to be very good. Effects of the mixed convection and elasticity parameters on the skin friction and heat transfer coefficients for a fluid having the Prandtl number equal to one are also discussed. It is found that for some values of the viscoelastic parameter and some negative values of the mixed convection parameter (opposing flow) the boundary layer separates from the cylinder. Heating the cylinder delays separation and can, if the cylinder is warm enough, suppress the separation completely. Similar to the case of a Newtonian fluid, cooling the cylinder brings the separation point nearer to the lower stagnation point. However, for a sufficiently cold cylinder there will not be a boundary layer.     

Numerical simulation of two-degree-of-freedom vortex-induced vibration of a circular cylinder close to a plane boundary

Two-degree-of-freedom vortex-induced vibrations (VIV) of a circular cylinder close to a plane boundary are investigated numerically. The Reynolds-Averaged Navier-Stokes (RANS) equations are solved using the Arbitrary Lagrangian Eulerian (ALE) scheme with a k-ω turbulence model closure. The numerical model is validated against experimental data of VIV of a cylinder in uniform flow and VIV of a cylinder close to a plane boundary at low mass ratios. The numerical results of the vibration mode, vibration amplitude and frequency agree well with the experimental data. VIV of a circular cylinder close to a plane boundary is simulated with a mass ratio of 2.6 and gap ratios of e/D=0.002 and 0.3 (gap ratio is defined as the ratio of gap between the cylinder and the bed (e) to cylinder diameter (D)). Simulations are carried out for reduced velocities ranging from 1 to 15 and Reynolds numbers ranging from 1000 to 15 000. It is found that vortex-induced vibrations occur even if the initial gap ratio is as small as e/D=0.002, although reported research indicated that vortex shedding behind a fixed circular cylinder is suppressed at small gap ratios (e/D<0.3 or 0.2). It was also found that vibration amplitudes are dependant on the bouncing back coefficient when the cylinder hits the plane boundary. Three vortex shedding modes are identified according to the numerical results: (i) single-vortex mode where the vortices are only shed from the top of the cylinder; (ii) vortex-shedding-after-bounce-back mode; (iii) vortex-shedding-before-bounce-back mode. It was found that the vortex shedding mode depends on the reduced velocity.     

Amplification/suppression of flow-induced motions of an elastically mounted circular cylinder by attaching tripping wires

This study has experimentally studied the effect of tripping wires on the vibration of a circular cylinder subject to flows. We placed a pair of tripping wires onto the surface of a circular cylinder symmetrically about the stagnation point, and submerged the cylinder clamped or elastically mounted. The test for the clamped cylinder (hereinafter stationary cylinder test) covered a wide domain of angular position from 15° to 165° to find what angular positions enhance or weaken hydrodynamic forces acting on the cylinder. The test for the moving cylinder used four angular positions: 60°, 75°, 105° and 120° to examine how the tripping wires amplify or suppress the flow-induced vibration of the cylinder. The stationary cylinder test revealed that with the angular position ranging from 20° to 52.5° or greater than 97.5° the hydrodynamic forces diminish and with the angular position ranging from 52.5° to 97.5° the forces increase. Particularly, positioning the tripping wires at 75° maximally increases the lift and drag coefficients by 63% and 44%, respectively, in comparison with the coefficients produced without the tripping wires, and these coefficients are maximally reduced by 67% and 20%, respectively, by positioning them at 112.5°. The moving cylinder test has elucidated that the angular positions at 60° and 75° intensify the vibration of the cylinder, involving a monotonic increase in the amplitudes with an increase in reduced velocity and wide lock-in range. With the tripping wires attached at 105° and 120° the vibrations of the cylinder are considerably suppressed, particularly when the angular position is 120° the vibration completely disappears.     

Fluid-dynamic forces and wake frequencies on a tilted rectangular cylinder near a solid wall

We investigate the impact of different boundary conditions on the flow field developing around a tilted rectangular cylinder with two different values of the aspect ratio (l/s=3 and 4). We are mainly interested in analyzing the changes in force coefficients and in the vortex shedding Strouhal number when the cylinder is placed at various distances from a bottom wall and different values of attack angle. The angle of attack ranges between −30° and +30° and the cylinder elevation above the bottom wall is varied between almost zero and 5 times the thickness of the cylinder. A large body of experimental results is related to the small elevation conditions at different attack angles, where the presence of the wall has a non-negligible effect on the behavior of the force coefficients and Strouhal number of the vortex shedding.     

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.     

An immersed interface method for flow past circular cylinder in the vicinity of a plane moving wall

Two‐dimensional flows past a stationary circular cylinder near a plane boundary are numerically simulated using an immersed interface method with second‐order accuracy. Instead of a fixed wall, a moving wall with no‐slip boundary is considered to avoid the complex involvement of the boundary layer and to focus only on the shear‐free wall proximity effects for investigating the force dynamics and flow fields. To analyze the convergence and accuracy of our implementation, numerical studies have been first performed on a simple test problem of rotational flow, where the second order of convergence is confirmed through numerical experiments and an optimal range of relative grid‐match ratio of Lagrangian to Eulerian grid sizes has been recommended. By comparing the force quantities and the Strouhal number, the accuracy of this method has been demonstrated on the flow past a stationary isolated cylinder. The cylinder is then put in proximity to the wall to investigate the shear‐free wall proximity effects in the low Reynolds number regime (20≤Re≤200). The gap ratio, e/D, where e denotes the gap between the cylinder and the moving wall and D denotes the diameter of the cylinder, is taken from 0.10 to 2.00 to determine the critical gap ratio, (e/D)_critical, for the alternate vortex shedding, where the fluid forces, flow fields and the streamwise velocity profiles are studied. One of the key findings is that the (e/D)_critical for the alternate vortex shedding decreases as the Reynolds number increases. We also find that, in this low Reynolds number regime, the mean drag coefficient increases and peaks at e/D = 0.5 with the increase of e/D and keeps decreasing gently from e/D = 0.5 to e/D = 2.0, while the mean lift coefficient decreases monotonically with the increase of e/D. New correlations are then proposed for computing force coefficients as a function of Re and e/D for a cylinder in the vicinity of a moving plane wall. Copyright © 2015 John Wiley & Sons, Ltd.     

Numerical simulation of laminar flow past a circular cylinder with slip conditions

The no‐slip condition is an assumption that cannot be derived from first principles and a growing number of literatures replace the no‐slip condition with partial‐slip condition, or Navier‐slip condition. In this study, the influence of partial‐slip boundary conditions on the laminar flow properties past a circular cylinder was examined. Shallow‐water equations are solved by using the finite element method accommodating SU/PG scheme. Four Reynolds numbers (20, 40, 80, and 100) and six slip lengths were considered in the numerical simulation to investigate the effects of slip length and Reynolds number on characteristic parameters such as wall vorticity, drag coefficient, separation angle, wake length, velocity distributions on and behind the cylinder, lift coefficient, and Strouhal number. The simulation results revealed that as the slip length increases, the drag coefficient decreases since the frictional component of drag is reduced, and the shear layer developed along the cylinder surface tends to push the separation point away toward the rear stagnation point so that it has larger separation angle than that of the no‐slip condition. The length of the wake bubble zone was shortened by the combined effects of the reduced wall vorticity and wall shear stress which caused a shift of the reattachment point closer to the cylinder. The frequency of the asymmetrical vortex formation with partial slip velocity was increased due to the intrinsic inertial effect of the Navier‐slip condition. Copyright © 2011 John Wiley & Sons, Ltd.     

Interaction between two circular cylinders in slow flow of Bingham viscoplastic fluid

The interaction between two circular cylinders was studied in the slow flow of a Bingham viscoplastic fluid in an infinite medium without any inertia effects. The configuration studied is that in which the flow direction is parallel to the centre line of the cylinders. Finite-element numerical simulations were used with an approximation by Papanastasiou's regularisation method. The case of high yield stress effect was particularly examined. The convergence of the solutions was examined in detail. Changes in the rigid zones, kinematics and stresses were determined in relation to the degree of interaction, which is a function of the distance between the cylinders and the effect of yield stress. The results compared with the case of a single cylinder show that yield stress reduces interaction effects. The transition between configurations with interacting cylinders and configuration with isolated cylinders was examined as a function of the effect of yield stress. Correlations were proposed for the drag coefficient and the stability criterion when the cylinders are interacting.     

Flow-induced vibration of a circular cylinder subjected to wake interference at low Reynolds number

Two- and three-dimensional numerical simulations of the flow around two circular cylinders in tandem arrangements are performed. The upstream cylinder is fixed and the downstream cylinder is free to oscillate in the transverse direction, in response to the fluid loads. The Reynolds number is kept constant at 150 for the two-dimensional simulations and at 300 for the three-dimensional simulations, and the reduced velocity is varied by changing the structural stiffness. The in-line centre-to-centre distance is varied from 1.5 to 8.0 diameters, and the results are compared to that of a single isolated flexible cylinder with the same structural characteristics, m^?=2.0 and ζ=0.007. The calculations show that significant changes occur in the dynamic behaviour of the cylinders, when comparing the flow around the tandem arrangements to that around an isolated cylinder: for the tandem arrangements, the lock-in boundaries are wider, the maximum displacement amplitudes are greater and the amplitudes of vibration for high reduced velocities, outside the lock-in, are very significant. The main responsible for these changes appears to be the oscillatory flow in the gap between the cylinders.     

Spanwise correlations in the wake of a circular cylinder and a trapezoid placed inside a circular pipe

The spanwise correlation of a circular cylinder and a trapezoidal bluff body placed inside a circular pipe in fully developed turbulent regime is studied using hotwire anemometer. The present configuration possesses complex fluid structure interaction owing to the following features: high blockage effect; low aspect ratio of the body; upstream turbulence and interaction of axisymmetric flow with a two dimensional bluff body. The spatial correlation of such configuration is seldom reported in the literature. Results are presented for Reynolds number of Re_D=1×10⁵. Three different blockage ratios (0.14, 0.19 and 0.28) are considered in the present study. Correlation coefficient is observed to improve with increase in blockage ratio. Compared to a circular cylinder, a trapezoidal bluff body possesses high correlation length. The near wall effects tend to increase the phase drift, which is reflected in low correlation coefficients close to the pipe wall. The results show that the simultaneous effect of curvature, low aspect ratio and upstream turbulence reduces the correlation coefficients significantly as compared to unconfined and confined (parallel channel) flows. The low frequency modulations with a circular cylinder are higher for lower blockage ratios. The three-dimensionality of vortex shedding for trapezoid with a blockage ratio of 0.28 was observed to be lower compared to circular cylinder and all other blockage ratios. Low frequency modulations were found to be responsible for weak vortex shedding from a circular cylinder compared to a trapezoidal bluff body. The vortex shedding is observed to be nearly two dimensional in case of a trapezoidal bluff body of blockage ratio 0.28.     

Low Reynolds number hydrodynamic interaction of a solid particle with a planar wall

The slow viscous flow problem of an arbitrary solid particle in motion near a planar wall is recast into a boundary integral formulation. The present formulation employs the Green function appropriate to the planar wall problem and is developed in sufficient generality to allow calculations for arbitrary particles in any base flow which satisfies Stokes equations and no-slip on the wall. The resulting integral equations are easily discretized and solved for the particle surface tractions. Calculations are performed for axisymmetric motions of a variety of ellips?ids near the planar wall. Agreement with existing theory is excellent.