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.     

Laminar flow and heat transfer in confined channel flow past square bars arranged side by side

A numerical investigation was conducted to analyze the unsteady laminar flow field and heat transfer characteristics in a plane channel with two square bars mounted side by side to the approaching flow. A finite volume technique is applied with a fine grid and time resolution. The transverse separation distance between the bars (G/d) is varied from 0 to 5, whereas the bar height to channel height is d/H=1/8, and the channel length is L=5H. Different flow regimes develop in the channel due the interaction between the two mounted square bars, steady flow, flow with vortex shedding synchronization either in phase or in anti-phase, or biased flow with low frequency modulation of vortex shedding are found. Results show that the pressure drop increase and heat transfer enhancement are strongly dependent of the transverse separation distance of the bars and the channel Reynolds number.     

Flow structures around an equilateral triangle arrangement of three spheres

This paper represents the results of an experimental study on the flow structure around a single sphere and three spheres in an equilateral-triangular arrangement. Flow field measurements were performed using a Particle Image Velocimetry (PIV) technique and dye visualization in an open water channel for a Reynolds number of Re = 5 × 10³ based on the sphere diameter. The distributions and flow features at the critical locations of the contours of the velocity fluctuations, the patterns of sectional streamlines, the vorticity contours, the turbulent kinetic energy, the Reynolds stress correlations and shedding frequency are discussed. The gap ratios (G/D) of the three spheres were varied in the range of 1.0 ⩽ G/D ⩽ 2.5 where G was the distance between the sphere centers, and D was the sphere diameter which was taken as 30 mm. Due to the interference of the shedding shear layers and the wakes, more complex features of the flow patterns can be found in the wake region of the two downstream spheres behind the leading sphere. For G/D = 1.25, a jet-like flow around the leading sphere through the gap between the two downstream spheres occurred, which significantly enhanced the wake region. It was observed that a continuous flow development involving shearing phenomena and the interactions of shedding vortices caused a high rate of fluctuations over the whole flow field although most of the time-averaged flow patterns were almost symmetric about the two downstream spheres.     

An Embedding Method for Bluff Body Flows: Interactions of Two Side-by-Side Cylinder Wakes

We introduce a simple method for the numerical simulation of bluff body flows where the solid object is represented by a distributed body force in the Navier–Stokes equations. The body force density is found at every time step to reduce the velocity within the computational cells occupied by the rigid body to a prescribed value. The method combines certain ideas from the immersed boundary method which was developed to treat biofluid mechanical flows and the volume-of-fluid method for simulating flows with fluid–fluid interfaces. The main advantage of this embedding method is that the computations can be effected on a regular Cartesian grid, without the need to fit the grid to the bluff body surfaces. Thus, flow past several complex bodies can be treated as easily as flow past a single body. The method is validated by reproducing well-established results for vortex shedding from a stationary cylinder. The flow past two side-by-side cylinders is then investigated. When the distance between the cylinders is small, they are seen to shed vortices in-phase, whereas for larger distances, the shedding occurs in anti-phase. For intermediate distances, various shedding patterns are observed, including quasi-periodic, asymmetric and chaotic regimes. Mean values and phase portraits associated with the cylinder lift and drag coefficients, as well as spectral analysis of the same data, are used to describe the flow. A transition diagram that can be compared with experiments or models outlines the various dynamical regimes as a function of the distance between the cylinders and the Reynolds number.     

串列双圆柱绕流问题的数值模拟

本文运用有限体积方法,对绕串列放置的双圆柱的二维不可压缩流动进行了数值计算。为研究两圆柱不同间距对圆柱相互作用和尾流特征的影响,选取间距比Ｌ／Ｄ（Ｌ为两圆柱中心间的距离,Ｄ为圆柱直径）在１．５～５．０之间每隔０．５共八个有代表性的间距进行了计算模拟。计算均在Ｒｅ＝２００条件下进行。计算结果表明：对该绕流问题,流动特征在很大程度上取决于间距的大小。且间距存在一临界值,间距比从小于临界值变化到大于临界     

Drag forces of interacting spheres in power-law fluids

Drag forces of interacting particles suspended in power-law fluid flows were investigated in this study. The drag forces of interacting spheres were directly measured by using a micro-force measuring system. The tested particles include a pair of interacting spheres in tandem and individual spheres in a cubic matrix of multi-sphere in flows with the particle Reynolds number from 0.7 to 23. Aqueous carboxymethycellulose (CMC) solutions and glycerin solutions were used as the fluid media in which the interacting spheres were suspended. The range of power-law index varied from 0.6 to 1.0. In conjunction to the drag force measurements, the flow patterns and velocity fields of power-law flows over a pair of interacting spheres were also obtained from the laser assisted flow visualization and numerical simulation.

Both experimental and computational results suggest that, while the drag force of an isolated sphere depends on the power-index, the drag coefficient ratio of an interacting sphere is independent from the power-law index but strongly depends on the separation distance and the particle Reynolds number. Our study also shows that the drag force of a particle in an assemblage is strongly positions dependent, with a maximum difference up to 38%.     

Numerical modeling of flow past a volumeless and thin rigid body using direct forcing immersed boundary method

The new capability has been added as the numerical method for modeling volumeless and thin rigid bodies to the direct forcing immersed boundary (DFIB) method. The DFIB approach is based on adding a virtual force to the Navier–Stokes equations of incompressible flow to account for the interaction between the fluid and structures. The volume of a solid function (VOS) identifies the stationary or moving solid structures in a given fluid domain. A new VOS-based algorithm was developed to identify thin, rigid structure boundary points in fluid flow and ensure that the fluid cannot cross through the boundary of a thin rigid structure while moving or stationary. The DFIB method was first validated in a three-dimensional (3D) turbulent flow over a circular cylinder. The large-eddy simulation simulated the turbulent flow scales. The proposed algorithm was tested using a 3D turbulent flow past a stationary and rotating Savonius wind turbine that functions as a thin, rigid body. The validation results showed that the selected DFIB approach, combined with the novel algorithm, could simulate a thin, volumeless, rigid structure that is stationary and rotating in incompressible turbulent flows. The current method is also applicable for two-way fluid-structure interaction problems.     

Dynamics of the flow around colliding spheres

An investigation of the flow resulting from the collision of two spheres at low Reynolds numbers is presented. Each sphere starts from rest and traverses a distance of 5 sphere diameters to the point of contact. Experimental and numerical results are compared for a symmetric collision; that is, a collision between two spheres of the same diameter and travelling with the same velocity. The flow consists of two axisymmetric recirculation zones which become a pair of colliding vortex rings, expanding radially from the collision point. Several examples of unbalanced collisions are also presented numerically, with one or both of the velocity and diameter of the spheres altered. These collisions break the symmetry, altering the post-collision expansion of the vortex rings.     

Vortex dynamics mechanisms of separated boundary layer in a highly loaded low pressure turbine cascade

This paper investigates the vortex dynamics in the suction-side boundary layer on an aero-engine low pressure turbine blade at two different Reynolds numbers at which short and long laminar separation bubbles occur. Different vortical patterns are observed and investigated through large eddy simulation (LES). The results show that at the higher Reynolds number, streamwise streaks exist upstream of separation line. These streaks initiate spanwise undulation in the form of vortex tubes, which roll-up and shed from the shear layer due to the Kelvin–Helmholtz instability. The vortex tubes alternately pair together and eventually distort and break down to small-scale turbulence structures near the mean reattachment location and convect into a fully turbulent boundary layer. At the lower Reynolds number, streamwise streaks are strong and the separated flow is unable to reattach to the blade surface immediately after transition to turbulence. Therefore, bursting of short bubbles into long bubbles can occur, and vortex tubes have larger diameters and cover a part of the blade span. In this case vortex pairing does not occur and vortex shedding process is promoted mainly by flapping phenomenon. Moreover, the results of dynamic mode decomposition (DMD) analysis show a breathing motion as a source of unsteadiness in the separation location, which is accompanied by the flapping phenomenon.     

Numerical analysis of unsteady cavitating turbulent flow and shedding horse-shoe vortex structure around a twisted hydrofoil

Cavitating turbulent flow around hydrofoils was simulated using the Partially-Averaged Navier–Stokes (PANS) method and a mass transfer cavitation model with the maximum density ratio (ρ_l/ρ_v,clip) effect between the liquid and the vapor. The predicted cavity length and thickness of stable cavities as well as the pressure distribution along the suction surface of a NACA66(MOD) hydrofoil compare well with experimental data when using the actual maximum density ratio (ρ_l/ρ_v,clip = 43391) at room temperature. The unsteady cavitation patterns and their evolution around a Delft twisted hydrofoil were then simulated. The numerical results indicate that the cavity volume fluctuates dramatically as the cavitating flow develops with cavity growth, destabilization, and collapse. The predicted three dimensional cavity structures due to the variation of attack angle in the span-wise direction and the shedding cycle as well as its frequency agree fairly well with experimental observations. The distinct side-lobes of the attached cavity and the shedding U-shaped horse-shoe vortex are well captured. Furthermore, it is shown that the shedding horse-shoe vortex includes a primary U-shaped vapor cloud and two secondary U-shaped vapor clouds originating from the primary shedding at the cavity center and the secondary shedding at both cavity sides. The primary shedding is related to the collision of a radially-diverging re-entrant jet and the attached cavity surface, while the secondary shedding is due to the collision of side-entrant jets and the radially-diverging re-entrant jet. The local flow fields show that the interaction between the circulating flow and the shedding vapor cloud may be the main mechanism producing the cavitating horse-shoe vortex. Two side views described by iso-surfaces of the vapor volume fraction for a 10% vapor volume, and a non-dimensional Q-criterion equal to 200 are used to illustrate the formation, roll-up and transport of the shedding horse-shoe vortex. The predicted height of the shedding horse-shoe vortex increases as the vortex moves downstream. It is shown that the shape of the horse-shoe vortex for the non-dimensional Q-criterion is more complicated than that of the 10% vapor fraction iso-surface and is more consistent with the experiments. Further, though the time-averaged lift coefficient predicted by the PANS calculation is about 12% lower than the experimental value, it is better than other predictions based on RANS solvers.     

Study on flow past two spheres in tandem arrangement using a local mesh refinement virtual boundary method

A local mesh refinement virtual boundary method based on a uniform grid is designed to study the transition between the flow patterns of two spheres in tandem arrangement for Re=250. For a small gap (L/D=1.5), the flow field is axisymmetric. As the spacing ratio increases to 2.0, the pressure gradient induces the circumferential fluid motion and a plane‐symmetric flow is constructed through a regular bifurcation. For L/D?2.5, the vortices are periodically shed from the right sphere, but the planar symmetry remains. The case for L/D=3.0 is picked up to give a detail investigation for the unsteady flow. The shedding frequency of vortical structure from the upper side of the right sphere is found to be double of the frequency of the lower side. With the flow spectra of various gaps given, the underlying competitive mechanism between the two shedding frequencies is studied and a critical spacing gap is revealed. Copyright © 2005 John Wiley & Sons, Ltd.     

Separation efficiency of liquid–solid undergoing vibration based on breakage of liquid bridge

Vibrating separation is a significant method for liquid–solid separation. A typical example is the vibrating screen to dewater wet granular matter. The properties of granular matter and the vibrating parameters significantly affect the separation efficiency. This study investigates the effect of vibration parameters in separation based on the breakage of large-scale liquid bridge numerically by using a calibrated simulation model. Through analysing the simulation results, the liquid bridge shape and the volume between two sphere particles for various particle sizes and particle distances were studied in the static condition under the effect of gravity. The results show a general reducing trend of liquid bridge volume when the radius ratio of two particles increases, particularly when the ratio increases to 5. Additionally, a set of vibrating motion was applied to the liquid bridge in the simulation model. A group of experiments were also performed to validate the simulation model with vibration. Then, the effect of vibrating peak acceleration, distance between spheres and radius on the separation efficiency which was reflected by the residual water were investigated. It is found that separation efficiency increased obviously with the peak acceleration and the increase slowed down after the peak acceleration over 1 m/s².     

HYDRODYNAMIC RESISTANCE EFFECT OF FLUID LAYER BETWEEN TWO IMMERSED APPROACHING PARTICLES

1. Introduction The mechanisms of impact and rebound of solid parti- cles in particulate flow systems are of interest over a wide range of application areas such as fluidized beds, pneu- matic transport, filtration processes, erosion and pollution control of suspended particles. In many cases, the colli- sions of particles against themselves and against walls may affect the properties of the mixture. Efforts have been made to describe the fundamental mechanics of particle collisions. The conta…     

Wake development in staggered short cylinder arrays within a channel

Staggered arrays of short cylinders, known as pin?Cfins, are commonly used as a heat exchange method in many applications such as cooling electronic equipment and cooling the trailing edge of gas turbine airfoils. This study investigates the near wake flow as it develops through arrays of staggered pin fins. The height-to-diameter ratio was unity while the transverse spacing was kept constant at two cylinder diameters. The streamwise spacing was varied between 3.46 and 1.73 cylinder diameters. For each geometric arrangement, experiments were conducted at Reynolds numbers of 3.0e3 and 2.0e4 based on cylinder diameter and velocity through the minimum flow area of the array. Time-resolved flowfield measurements provided insight into the dependence of row position, Reynolds number, and streamwise spacing. Decreasing streamwise spacing resulted in increased Strouhal number as the near wake length scales were confined. In the first row of the bundle, low Reynolds number flows were mainly shear-layer-driven while high Reynolds number flows were dominated by periodic vortex shedding. The level of velocity fluctuations increased for cases having stronger vortex shedding. The effect of streamwise spacing was most apparent in the reduction of velocity fluctuations in the wake when the spacing between rows was reduced from 2.60 diameters to 2.16 diameters.     

Experimental wall correction factors of single solid spheres in triangular and square cylinders,and parallel plates

To examine wall effects, the steady state settling velocities of single solid spheres of various diameters were measured in triangular and square cylinders, and between parallel plates, which were filled up with highly viscous aqueous milet jelley solutions. The wall effect for the steady state motion of single spheres in a narrow space has to date been correlated by the wall correction factor. The wall correction factors of single solid spheres were correlated in such apparatus in the creeping flow region.     

Pipe flow with solid particles fixed in space

A group of solid particles were hung by slender rods in a pipe to make a model of two-phase flow of coarse particles. Pressure gradient and velocities were measured for different types of the models. The drag on the particles (spheres) were obtained from measurements of pressure gradient with some assumptions. The results are summarized as follows. (1) Mean velocities of fluid are lower in the central part of the pipe than in the circumferential part. Turbulence is remarkably increased by particles. The spectrum distribution of turbulent velocity becomes flatter. These results are similar to the gas-solid flow of coarse particles in a vertical pipe. (2) At a large Reynolds number, the drag coefficient per one sphere in the group is larger than that of a single isolated sphere in a uniform flow. When the spheres are arranged along the same line in the longitudinal direction, the drag coefficient becomes smaller as the longitudinal distance between the spheres is shortened.     

Wall Region Porosity Distributions for Packed Beds of Uniform Spheres with Modified and Unmodified Walls

A mathematical derivation of the porosity (local void fraction) distribution near the walls of packed beds of uniform spheres is presented. This investigation supports the study of methods of reducing or eliminating the so-called wall effect, or bypass flow, which accompanies the increase in porosity when spheres come in contact with a solid boundary. Limiting the amount of bypass flow is important in certain applications such as in packed bed nuclear reactors where bypass flow allows some coolant to avoid the high power density core region. Four basic porosity distributions are determined. The first investigates the case where spheres against a flat wall are packed in the tightest possible packing density. This density is then reduced by changing the sphere spacing until the minimum porosity matches that obtained experimentally. In the other cases, the effect of various ways of embedding spheres in the wall on the wall region porosity is examined. By partially embedding spheres in the wall, the porosity at the wall is reduced and the most direct cause of the bypass flow is thereby eliminated. The porosity is found by evaluating the ratio of the solid area to total area in a plane which is parallel to the wall. The local porosity is derived as a function of distance from the wall in the region within one-half a sphere diameter from the wall. The average porosity of the wall region is also calculated. This research has application to flow situations such as packed bed chemical reactors, pebble bed nuclear reactors and flow in packed beds.     

On the interaction between two fixed spherical particles

The variation of the drag (C_D) and lift coefficients (C_L) of two fixed solid spherical particles placed at different positions relative each other is studied. Simulations are carried out for particle Reynolds numbers of 50, 100 and 200 and the particle position is defined by the angle between the line connecting the centers of the particles and the free-stream direction (α) and the separation distance (d₀) between the particles. The flow around the particles is simulated using two different methods; the Lattice Boltzmann Method (LBM), using two different computational codes, and a conventional finite difference approach, where the Volume of Solid Method (VOS) is used to represent the particles. Comparisons with available numerical and experimental data show that both methods can be used to accurately resolve the flow field around particles and calculate the forces the particles are subjected to. Independent of the Reynolds number, the largest change in drag, as compared to the single particle case, occurs for particles placed in tandem formation. Compared to a single particle, the drag reduction for the secondary particle in tandem arrangement is as high as 60%, 70% and 80% for Re = 50, 100 and 200, respectively. The development of the recirculation zone is found to have a significant influence on the drag force. Depending on the flow situation in-between the particles for various particle arrangements, attraction and repulsion forces are detected due to low and high pressure regions, respectively. The results show that the inter-particle forces are not negligible even under very dilute conditions.     

Steady shear rheology of dilute polystyrene particle gels

This paper describes an experimental study on dispersions of monodisperse polystyrene (PS) spheres with a typical radius of 1 μm, dispersed in an electrolyte at high ionic strength, screening the electrostatic repulsion. These suspensions gelate at rest even at low volume fractions of PS particles. The density of the particles is matched with the solvent by using deuterium oxide for volume fractions φ≤0.117. Steady-state flow curves, viscosity as a function of shear rate, are measured and reported for 0.014<φ<0.322. The measured flow curves are analyzed on the basis of two models: 1. In the giant floc model (van Diemen and Stein 1983, 1984; Schreuder et al. 1986, 1987; Laven et al. 1988), at low shear rates, the shear is not distributed homogeneously but is limited to certain shear planes; the energy dissipation during steady flow is due primarily to overcoming the viscous drag on the suspended particles during motion caused by encounters of particles in the shear planes. Though this model was developed for higher solid volume fractions (0.35–0.425), we found that it also describes the rheology of dilute particle gels for 0.15≤φ≤0.3, using the same values for the parameters in the model as in the high solid volume fraction region. For φ<0.15, the model also describes the data if the fraction of distance by which a moving particle entrains its neighbors, is assumed to increase in this φ region. 2. The model of de Rooij (de Rooij et al. 1993, 1994) considers aggregates in shear flow to be monodisperse impermeable spheres with a fractal structure. The permeability is taken into account by considering a hydrodynamic radius smaller than the gyration radius in the Krieger-Dougherty expression for the hydrodynamic contribution to the viscosity. Through the use of a yield criterion the aggregate radius is modeled as a function of shear rate. We found that the model describes our experimental results, with a combination of parameter values used already by de Rooij, but only for φ<0.15. Received: 7 May 1998 Accepted: 22 December 1998     

Pore-Scale Modeling of Viscous Flow and Induced Forces in Dense Sphere Packings

We propose a method for effectively upscaling incompressible viscous flow in large random polydispersed sphere packings: the emphasis of this method is on the determination of the forces applied on the solid particles by the fluid. Pore bodies and their connections are defined locally through a regular Delaunay triangulation of the packings. Viscous flow equations are upscaled at the pore level, and approximated with a finite volume numerical scheme. We compare numerical simulations of the proposed method to detailed finite element simulations of the Stokes equations for assemblies of 8–200 spheres. A good agreement is found both in terms of forces exerted on the solid particles and effective permeability coefficients.