Wakes of elliptical cylinders at low Reynolds number

The wakes of elliptical cylinders are numerically investigated at a Reynolds number Re_D = 150. ANSYS-Fluent, based on the finite volume method, is used to simulate two-dimensional Newtonian fluid flow. The cylinder cross-sectional aspect ratio (AR) is varied from 0.25 to 1.0 (circular cylinder), and the angle of attack (α) of the cylinder is changed as α = 0° – 90°. With the changes in AR and α, three distinct wake patterns (patterns I, II, III) are observed, associated with different characteristics of fluid forces. Steady wake (pattern I) is characterised by two steady bubbles forming behind the cylinder, occurring at AR < 0.37 and α < 2.5°. Time-mean drag and fluctuating lift coefficients are small. Pattern II refers to Karman wake followed by steady wake (AR ≥ 0.37 – 0.67, depending on α) with the Karman street transitioning to two steady shear layers downstream. An inflection angle α_i is identified where the time-mean drag of the elliptical cylinder is identical to that of a circular cylinder. Pattern III is the Karman wake followed by secondary wake (AR ≤ 0.67, α > 52°), where the Karman street forming behind the cylinder is modified to a secondary vortex street with a low frequency. The Time-mean drag coefficient is maximum for this pattern.     

Flow characteristics and flow-induced forces of a stationary and rotating triangular cylinder with different incidence angles at low Reynolds numbers

In this paper, the problem of two-dimensional fluid flow past a stationary and rotationally oscillating equilateral triangular cylinder with a variable incident angle, Reynolds number, oscillating amplitude, and oscillating frequency is numerically investigated. The computations are carried out by using a two-step Taylor-characteristic-based Galerkin (TCBG) algorithm. For the stationary cases, simulations are conducted at various incident angles of α=0.0–60.0° and Reynolds numbers of Re=50–160. For the oscillation cases, the investigations are done at various oscillating amplitudes of θ_max=7.5–30.0° and oscillating frequencies of F_s/F_o=0.5–3.0 considering two different incidence angles (α=0.0°, 60.0°) and three different Reynolds numbers (Re=50, 100, 150). The results show that the influences of key parameters (incidence angle, Reynolds number, oscillating amplitude, and oscillating frequency) are significant on the flow pattern and hydrodynamic forces. For the stationary cases, at smaller angle of incidence (α≤30.0°), Reynolds number has a large impact on the position of the separation points. When α is between 30.0° and 60.0°, it was found that the separation points are located at the rear corners. From a topological point of view, the diagram of flow pattern is summarized, including two distinct patterns, namely, main separation and vortex merging. A deep analysis of the influence of Reynolds number and incidence angles on the mean pressure coefficient along the triangular cylinder surface is presented. Additionally, for the oscillating cases, the lock-on phenomenon is captured. The dominant flow patterns are 2S mode and P+S mode in lock-on region at α=0.0°. It is found at α=60.0°, however, that the flow pattern is predominantly 2S mode. Furthermore, except for the case of F_s/F_o=2.0, the mean drag decreases as the oscillating amplitude increases for each Reynolds number at α=0.0°. At α=60.0°, the minimum mean drag for F_s/F_o=1.5 is lower than that for stationary case, and occurs at θ_max=15.0° (Re=100) and θ_max=22.5° (Re=150), respectively. Finally, the effect of Reynolds number on a rotational oscillation cylinder is elucidated.     

Direct numerical simulation of three-dimensional flow past a yawed circular cylinder of infinite length

Direct numerical simulation of flow past a stationary circular cylinder at yaw angles (α) in the range of 0–60° was conducted at Reynolds number of 1000. The three-dimensional (3-D) Navier–Stokes equations were solved using the Petrov–Galerkin finite element method. The transition of the flow from 2-D to 3-D was studied. The phenomena that were observed in flow visualization, such as the streamwise vortices, the vortex dislocation and the instability of the shear layer, were reproduced numerically. The effects of the yaw angle on wake structures, vortex shedding frequency and hydrodynamic forces of the cylinder were investigated. It was found that the Strouhal number at different yaw angles (α) follows the independence principle. The mean drag coefficient agrees well with the independence principle. It slightly increases with the increase of α and reaches a maximum value at α=60°, which is about 10% larger than that when α=0°. The root-mean-square (r.m.s.) values of the lift coefficient are noticeably dependent on α.     

Three-dimensional simulations of flow past two circular cylinders in side-by-side arrangements at right and oblique attacks

The flow past two identical circular cylinders in side-by-side arrangements at right and oblique attack angles is numerically investigated by solving the three-dimensional Navier–Stokes equations using the Petrov–Galerkin finite element method. The study is focused on the effect of flow attack angle and gap ratio between the two cylinders on the vortex shedding flow and the hydrodynamic forces of the cylinders. For an oblique flow attack angle, the Reynolds number based on the velocity component perpendicular to the cylinder span is defined as the normal Reynolds number Re_N and that based on the total velocity is defined as the total Reynolds number Re_T. Simulations are conducted for two Reynolds numbers of Re_N=500 and Re_T=500, two flow attack angles of α=0° and 45° and four gap ratios of G/D=0.5, 1, 3 and 5. The biased gap flow for G/D=0.5 and 1 and the flip-flopping bistable gap flow for G/D=1 are observed for both α=0° and 45°. For a constant normal Reynolds number of Re_N=500, the mean drag and lift coefficients at α=0° are very close to those at α=45°. The difference between the root mean square (RMS) lift coefficient at α=0° and that at α=45° is about 20% for large gap ratios of 3 and 5. From small gap ratios of 0.5 and 1, the RMS lift coefficients at α=0° and 45° are similar to each other. The present simulations show that the agreement in the force coefficients between the 0° and 45° flow attack angles for a constant normal Reynolds number is better than that for a constant total Reynolds number. This indicates that the normal Reynolds number should be used in the implementation of the independence principle (i.e., the independence of the force coefficients on the flow attack angle). The effect of Reynolds number on the bistable gap flow is investigated by simulating the flow for Re_N=100–600, α=0° and 45° and G/D=1. Flow for G/D=1 is found to be two-dimensional at Re_N=100 and weak three-dimensional at Re_N=200. While well defined biased flow can be identified for Re_N=300–600, the gap flow for Re_N=100 and 200 changes its biased direction too frequently to allow stable biased flow to develop.     

Flow around four cylinders arranged in a square configuration

This paper presents an experimental study of the flow around four circular cylinders arranged in a square configuration. The Reynolds number was fixed at Re=8000, the pitch-to-diameter ratio between adjacent cylinders was varied from P/D=2 to 5 and the incidence angle was changed from α=0° (in-line square configuration) to 45° (diamond configuration) at an interval of 7.5°. The flow field was measured using digital Particle Image Velocimetry (PIV) to examine the vortex shedding characteristics of the cylinders, together with direct measurement of fluid dynamic forces (lift and drag) on each cylinder using a piezoelectric load cell. Depending on the pitch ratio, the flow could be broadly classified as shielding regime (P/D≤2), shear layer reattachment regime (2.5≤P/D≤3.5) and vortex impinging regime (P/D≥4). However, this classification is valid only in the case that the cylinder array is arranged nearly in-line with the free stream (α≈0°), because the flow is also sensitive to α. As α increases from 0° to 45°, each cylinder experiences a transition of vortex shedding pattern from a one-frequency mode to a two-frequency mode. The flow interference among the cylinders is complicated, which could be non-synchronous, quasi-periodic or synchronized with a definite phase relationship with other cylinders depending on the combined value of α and P/D. The change in vortex pattern is also reflected by some integral parameters of the flow such as force coefficients, power spectra and Strouhal numbers.     

Control of two-degree-of-freedom vortex-induced vibrations of a circular cylinder using a pair of synthetic jets at low Reynolds number: Influence of position angle and momentum coefficient

In order to understand the effects of synthetic jets on the active control of two-degree-of-freedom (2DOF) vortex-induced vibrations (VIVs) of a circular cylinder, a series of numerical simulation were carried out at Reynolds number of 150. The synthetic jet excitation frequency was fixed at five times of the natural frequency of the cylinder in still water. The influence of two key parameters of synthetic jets, the position angle (α) and the momentum coefficient (Cu), on 2DOF VIVs was analyzed. Results indicated that both in-flow oscillation and cross-flow oscillation can be suppressed when the synthetic jets with sufficient momentum coefficient were positioned at the circular cylinder's leeward side (0° ≤ α ≤ 75°). When Cu = 4, 15° ≤ α ≤ 60°, the reductions of cross-flow and in-flow oscillation amplitudes were all larger than 99% and 70%, respectively. Besides, the in-flow oscillation frequency was locked-in to the excitation frequency of synthetic jets when the in-flow oscillation was effectively suppressed. A symmetric wake can be observed when the cross-flow oscillation was completely suppressed, and the 2P+2S vortex pattern can be observed when Cu = 4, α =165° and 180°.     

3-D measurements of separated flow over rigid plates with spanwise periodic cambering at low Reynolds number

Measurements of the flow field around a flat plate and rigid plates with spanwise periodic cambering were performed using volumetric three-component velocimetry (V3V) at a Reynolds numbers of 28,000 at α=12° where the flow is fully separated. The Reynolds normal and shear stresses, and the streamwise, spanwise and normal components of the vorticity vector are investigated for three-dimensionality. Flow features are discussed in context of the periodic cambering and corresponding aerodynamic force measurements. The periodic cambering results in spanwise variation in the reversed-flow region, Reynolds stresses and spanwise vorticity. These spanwise variations are induced by streamwise and normal vortices of opposite directions of rotation. Moreover, measurements were carried out for the cambered plates at α=8°, where a long separation bubble exists, to further understand the behavior of the streamwise and normal vortices. These vortices become more organized and increase in strength and size at the lower angle of attack. It is also speculated that these vortices contribute to the increase in lift at and beyond the onset of stall angle of attack.     

Dynamic interaction of fixed dual spheres for several configurations and inflow conditions

The changes in force characteristics as well as the shedding patterns for various dual sphere configurations are studied. The Reynolds numbers considered are 300, 600 and two different inflow conditions are used: steady and pulsating. The sphere formations are defined by the separation distance D₀ between the spheres and the angle between the line connecting the centres of the spheres and the main flow direction, γ. The position of one of the spheres is varied in the range 0°–90° using a 15° increment. Two separation distances are studied; 1.5D and 3D. The method used for the simulations is the Volume of Solid (VOS) approach, a method based on Volume of Fluid (VOF). A major conclusion from this work is that the sphere interaction alters the wake dynamics by obstructing the vortex shedding (generating a steady wake or a wake with lower Strouhal number) and by changing the direction of the lift force so that it in most cases is directed in the plane containing the sphere centres. The results also show that changing the inflow condition gives the same relative change in drag and lift as for a single sphere. The drag is substantially reduced by placing the sphere downstream in a tandem arrangement and slightly increased in a side-by-side arrangement. However, the effect is decreased by increasing separation distance and increasing Reynolds number.     

Scaling of flow separation on a pitching low aspect ratio plate

We use two different dye injection approaches, in two different water tunnels, to visualize the formation and subsequent evolution of leading-edge vortices and related separated structures, for a pitching low aspect ratio plate. The motion is a smoothed linear pitch ramp from 0° to 40° incidence, brief hold, and return to 0°, executed at reduced pitch rates ranging from 0.1 to 0.35 and about various pivot locations. All cases evince a leading edge vortex with pronounced axial flow, which leads to formation of large-scale, three-dimensional flow structures, culminating in a large vortical structure centered at the wing symmetry plane. Pitch is also compared to plunge, where the plunge-induced angle of attack is taken as the geometric pitch incidence angle, ignoring pitch-rate effects. At successively increasing values of convective time C/U, the three-dimensional patterns of the flow structure are remarkably similar for the pitching and plunging motions. The similarity of these patterns persists, though they are shifted in time, for variation of either the location of the pitching axis or the dimensionless pitch rate.     

Effect of tine rake angle and aspect ratio on soil failure patterns in dry loam soil

Studies were conducted in a laboratory glass-sided soil bin with dry compact loam soil (c = 0.02 kPa, Φ = 20° and cone index 210 kPa) with the specific objective of observing the effect of flat tine rake angle and aspect ratio on soil failure patterns. The tine was moved in the soil in a quasi-static condition and soil failure patterns were observed through a glass window. Tine rake angles of 50°, 90° and 130° were used while aspect ratio effects were studied by varying both width and depth of the tine. Individual effects of width and depth were investigated by maintaining a constant aspect ratio of 2.0 but varying width and depth. Results obtained indicated that soil failure patterns are affected by tine design parameters. Soil failure patterns were observed to be of progressive shear type in all cases. For 50° rake angle tines, the patterns consisted of inclined shear lines starting from the tine tip and gradually moving upwards towards the horizontal soil surface, intersecting it at an average failure angle of 32°. In the case of 90° rake angle tines, the inclined shear surface was at a distance from the tine tip whereas, for 130° rake angle tines, prismatic-shaped stationary soil wedges were formed adjacent to the tine. Failure angles for the 90° and 130° rake angle tines were almost the same as those for 50° rake angle tines. The results of this study also indicated that aspect ratio alone cannot account for changes in soil failure patterns, their corresponding soil reactions, forward rupture or surcharge profiles. The effects are mainly due to the individual changes in width and depth. There were no distinct zones as described in the passive soil pressure theory. Soil failures were in regular cycles resulting in corresponding variations in the soil reactions on the tines.     

Influence of inlet shear on the 3-D flow past a square cylinder at moderate Reynolds number

Unsteady three-dimensional (3-D) numerical simulations of linear shear flow past a square cylinder at moderate Reynolds number (Re=200) are performed. The shear parameter (K) considered in this study is varied as 0.0, 0.1, and 0.2. For the uniform flow (K=0.0) case, the chosen Re falls in the transition Reynolds number range. The low frequency force pulsations of square cylinder transition phenomena are observed to decrease with increasing shear parameter. The evolution of streamwise vortical structures indicates a mode A spanwise instability in the uniform flow. Unlike in uniform flow, mixed mode A and mode B spanwise instability is observed in the case of a shear flow. The autocorrelation function of the lift and the drag coefficients is improved for any particular separation distance with increasing K.     

Mixing efficiency of Y-type micromixers with different angles

The angle effect (α = +30°, ±60°, ±90°, and ±120°) of Y-type micromixers has been extensively studied using micro Laser Induced Fluorescence (μLIF) and micro Particle Image Velocimetry (μPIV) optical techniques to visualize and quantify concentration and velocity fields in order to examine mixing performance. An optimal α was obtained based upon mixing efficiency and mixing length. In addition, the effect of Joule heating on mixing was studied and discussed here.     

Analytical solutions for tapering quadratic and cubic rat-holes in highly frictional granular solids

New exact analytical solutions are presented for both stress and velocity fields for a Coulomb–Mohr granular solid assuming non-dilatant double-shearing theory. The solutions determined apply to highly frictional materials for which the angle of internal friction φ is assumed equal to 90°. This major assumption is made primarily to facilitate exact analytical solutions, and it is discussed at length in the Introduction, both in the context of real materials which exhibit large angles of internal friction, and in the context of using the solutions derived here as the leading term in a regular perturbation solution involving powers of 1−sinφ. The analytical velocity fields so obtained are illustrated graphically by showing the direction of the principal stress as compared to the streamlines. The stress solutions are also exploited to determine the static stress distribution for a granular material contained within vertical boundaries and a horizontal base, which is assumed to have an infinitesimal central outlet through which material flows until a rat-hole of parabolic or cubic profile is obtained, and no further flow takes place. A rat-hole is a stable structure that may form in storage hoppers and stock-piles, preventing any further flow of material. Here we consider the important problems of two-dimensional parabolic rat-holes of profile y=ax², and three-dimensional cubic rat-holes of profile z=ar³, which are both physically realistic in practice. Analytical solutions are presented for both two and three-dimensional rat-holes for the case of a highly frictional granular solid, which is stored at rest between vertical walls and a horizontal rigid plane, and which has an infinitesimal central outlet. These solutions are bona fide exact solutions of the governing equations for a Coulomb–Mohr granular solid, and satisfy exactly the free surface condition along the rat-hole surface, but approximate frictional conditions along the containing boundaries. The analytical solutions presented here constitute the only known solutions for any realistic rat-hole geometry, other than the classical solution which applies to a perfectly vertical cylindrical cavity.     

Simulation of non-Newtonian ink transfer between two separating plates for gravure-offset printing

The inks used in gravure-offset printing are non-Newtonian fluids with higher viscosities and lower surface tensions compared to Newtonian fluids. This paper examines the transfer of a non-Newtonian ink between two parallel plates when the top plate is moved upward with a constant velocity while the bottom plate is held fixed. Numerical simulations were carried out using the Carreau model to explore the behavior of a non-Newtonian ink in gravure-offset printing. The volume of fluid (VOF) model was adopted to demonstrate the stretching and break-up behaviors of the ink. The results indicate that the ink transfer ratio is greatly influenced by the contact angle, especially the contact angle at the upper plate (α). For lower values of α, oscillatory or unstable behavior of the position of minimum thickness of the ink between the two parallel plates during the stretching period is observed. This oscillation gradually diminishes as the contact angle at the upper plate is increased. Moreover, the number of satellite droplets increases as the velocity of the upper plate is increased. The surface tension of the conductive ink shows a positive impact on the ink transfer ratio to the upper plate. Indeed, the velocity of the upper plate has a significant influence on the ink transfer in gravure-offset printing when the Capillary number (Ca) is greater than 1 and the surface tension dominates over the ink transfer process when Ca is less than 1.     

Experimental study of pitching and plunging airfoils at low Reynolds numbers

Measurements of the unsteady flow structure and force time history of pitching and plunging SD7003 and flat plate airfoils at low Reynolds numbers are presented. The airfoils were pitched and plunged in the effective angle of attack range of 2.4°–13.6° (shallow-stall kinematics) and ?6° to 22° (deep-stall kinematics). The shallow-stall kinematics results for the SD7003 airfoil show attached flow and laminar-to-turbulent transition at low effective angle of attack during the down stroke motion, while the flat plate model exhibits leading edge separation. Strong Re-number effects were found for the SD7003 airfoil which produced approximately 25 % increase in the peak lift coefficient at Re = 10,000 compared to higher Re flows. The flat plate airfoil showed reduced Re effects due to leading edge separation at the sharper leading edge, and the measured peak lift coefficient was higher than that predicted by unsteady potential flow theory. The deep-stall kinematics resulted in leading edge separation that led to formation of a large leading edge vortex (LEV) and a small trailing edge vortex (TEV) for both airfoils. The measured peak lift coefficient was significantly higher (~50 %) than that for the shallow-stall kinematics. The effect of airfoil shape on lift force was greater than the Re effect. Turbulence statistics were measured as a function of phase using ensemble averages. The results show anisotropic turbulence for the LEV and isotropic turbulence for the TEV. Comparison of unsteady potential flow theory with the experimental data showed better agreement by using the quasi-steady approximation, or setting C(k) = 1 in Theodorsen theory, for leading edge–separated flows.     

Numerical simulation of shear effect on vortex shedding behind a square cylinder

This paper is concerned with the numerical simulation of the flow structure around a square cylinder in a uniform shear flow. The calculations were conducted by solving the unsteady 2D Navier–Stokes equations with a finite difference method. The effect of the shear parameter K of the approaching flow on the vortex-shedding Strouhal number and the force coefficients acting on the square cylinder is investigated in the range K=0·0–0·25 at various Reynolds numbers from 500 to 1500. The computational results are compared with some existing experimental data and previous studies. The effect of shear rate on the Strouhal number and the force acting on the cylinder has a tendency to reduce the oscillation. The Strouhal number, mean drag and amplitude of the fluctuating force tend to decrease as the shear rate increases, but show no significant change at low shear rate. Increasing the Reynolds number decreases the Strouhal number and increases the force acting on the cylinder. At high shear rate the shedding frequencies of the fluctuating drag and lift coefficients are identical. © 1997 John Wiley & Sons, Ltd.     

The flexibility effect of a plate according to various angles of attack in a free-stream

This paper presents a computational fluid–structure interaction analysis for a flexible plate in a free-stream to investigate the effects of flexibility and angle of attack on force generation. A Lattice Boltzmann Method with an immersed boundary technique using a direct forcing scheme model of the fluid is coupled to a finite element model with rectangular bending elements. We investigated the effects of various angles of attack of a flexible plate fixed at one of the end edges in a free-stream at a Reynolds number of 5000, which represents the wing flapping condition of insects and small birds in nature. The lift of the flexible plate is maintained at the large angle of attack, whereas the rigid plate shows the largest lift at angles of attack around 30–40° and then drastic reductions in the lift at the large angle of attack. If we consider the efficiency as the lift divided by the drag, the flexible plate shows better efficiency at angles of attack greater than 30° compared to the rigid plate. The better performance of the flexible plate at large angles of attack comes from the deformation of the plate, which produces an interaction between the trailing edge vortex and the short edge vortex. The horseshoe-shaped vortex produced by a large vortex interaction at the trailing edge side has an important role in increasing the lift, and the small projection area due to the deformation reduces the drag. Furthermore, we investigate the role of flexibility on the lift and the drag force of the rectangular plate in a free-stream as the Reynolds number increases. Whenever a large vortex interaction at the trailing edge side is shown, the efficiency of the rectangular plate is improved. Especially, the flexible plate shows better efficiency as the Reynolds number increases regardless of the angle of attack.     

Effects of boundary slippage on thin-film lubrication between two nonparallel plane plates

Hydrodynamic lubrications between two planeplates with an intersection angle θ have been investigated using the boundary slippage theory,and relations are obtainedbetween dimensionless pressures and coordinate x,betweenbearing capacity,friction force,friction coefficient and dimensionless slipping size factor.The results show that bearing capacity of two plane plates without boundary slippagesignificantly increases with increasing intersection angle θwhen 0 θ 1,whereas decreases with increasing intersection angle θ when θ 1.The results also show thatnegative pressure occurs in fluid entrance region and bearingcapacity decreases,and friction force and friction coefficientincrease with the increase of dimensionless slipping size factor.     

History force on a sphere in a weak linear shear flow

A numerical study of history forces acting on a spherical particle in a linear shear flow, over a range of finite Re, is presented. In each of the cases considered, the particle undergoes rapid acceleration from Re₁ to Re₂ over a short-time period. After acceleration, the particle is maintained at Re₂ in order to allow for clean extraction of drag and lift kernels. Good agreement is observed between current drag kernel results and previous investigations. Furthermore, ambient shear is found to have little influence on the drag kernel. The lift kernel is observed to be oscillatory, which translates to a non-monotonic change in lift force to the final steady state. In addition, strong dependence on the start and end conditions of acceleration is observed. Unlike drag, the lift history kernel scales linearly with Reynolds number and shear rate. This behavior is consistent with a short-time inviscid evolution. A simple expression for the lift history kernel is presented.