Motion of spheroid particles in shear flow with inertia

2014 ◽  
Vol 749 ◽  
pp. 145-166 ◽  
Author(s):  
Wenbin Mao ◽  
Alexander Alexeev
Keyword(s):  
Reynolds Number ◽  
Angular Velocity ◽  
Shear Flow ◽  
Rotation Period ◽  
Reynolds Numbers ◽  
Stokes Number ◽  
Simple Shear Flow ◽  
Fluid Inertia ◽  
Particle Rotation ◽  
Particle Inertia

AbstractIn this article, we investigate the motion of a solid spheroid particle in a simple shear flow. Using a lattice Boltzmann method, we examine individual effects of fluid inertia and particle rotary inertia as well as their combination on the dynamics and trajectory of spheroid particles at low and moderate Reynolds numbers. The motion of a single spheroid is shown to be dependent on the particle Reynolds number, particle aspect ratio, particle initial orientation and the Stokes number. Spheroids with random initial orientations are found to drift to stable orbits influenced by fluid inertia and/or particle inertia. Specifically, prolate spheroids drift towards the tumbling mode of motion, whereas oblate spheroids drift to the rolling mode. The rotation period and the variation of angular velocity of tumbling spheroids decrease as Stokes number increases. With increasing Reynolds number, both the maximum and minimum values of angular velocity decrease, whereas the particle rotation period increases. We show that particle inertia does not affect the hydrodynamic torque on the particle. We also demonstrate that superposition can be used to estimate the combined effect of fluid inertia and particle inertia on the dynamics of spheroid particles at sufficiently low Reynolds numbers.

Effect of fluid inertia on the dynamics and scaling of neutrally buoyant particles in shear flow

2013 ◽  
Vol 738 ◽  
pp. 563-590 ◽  
Author(s):  
T. Rosén ◽  
F. Lundell ◽  
C. K. Aidun
Keyword(s):  
Reynolds Number ◽  
Steady State ◽  
Aspect Ratio ◽  
Shear Flow ◽  
Stokes Number ◽  
Spheroidal Particle ◽  
Fluid Inertia ◽  
Prolate Spheroidal ◽  
Particle Inertia ◽  
Neutrally Buoyant

AbstractThe basic dynamics of a prolate spheroidal particle suspended in shear flow is studied using lattice Boltzmann simulations. The spheroid motion is determined by the particle Reynolds number (${\mathit{Re}}_{p} $) and Stokes number ($\mathit{St}$), estimating the effects of fluid and particle inertia, respectively, compared with viscous forces on the particle. The particle Reynolds number is defined by ${\mathit{Re}}_{p} = 4G{a}^{2} / \nu $, where $G$ is the shear rate, $a$ is the length of the spheroid major semi-axis and $\nu $ is the kinematic viscosity. The Stokes number is defined as $\mathit{St}= \alpha \boldsymbol{\cdot} {\mathit{Re}}_{p} $, where $\alpha $ is the solid-to-fluid density ratio. Here, a neutrally buoyant prolate spheroidal particle ($\mathit{St}= {\mathit{Re}}_{p} $) of aspect ratio (major axis/minor axis) ${r}_{p} = 4$ is considered. The long-term rotational motion for different initial orientations and ${\mathit{Re}}_{p} $ is explained by the dominant inertial effect on the particle. The transitions between rotational states are subsequently studied in detail in terms of nonlinear dynamics. Fluid inertia is seen to cause several bifurcations typical for a nonlinear system with odd symmetry around a double zero eigenvalue. Particle inertia gives rise to centrifugal forces which drives the particle to rotate with the symmetry axis in the flow-gradient plane (tumbling). At high ${\mathit{Re}}_{p} $, the motion is constrained to this planar motion regardless of initial orientation. At a certain critical Reynolds number, ${\mathit{Re}}_{p} = {\mathit{Re}}_{c} $, a motionless (steady) state is created through an infinite-period saddle-node bifurcation and consequently the tumbling period near the transition is scaled as $\vert {\mathit{Re}}_{p} - {\mathit{Re}}_{c} {\vert }^{- 1/ 2} $. Analyses in this paper show that if a transition from tumbling to steady state occurs at ${\mathit{Re}}_{p} = {\mathit{Re}}_{c} $, then any parameter $\beta $ (e.g. confinement or particle spacing) that influences the value of ${\mathit{Re}}_{c} $, such that ${\mathit{Re}}_{p} = {\mathit{Re}}_{c} $ as $\beta = {\beta }_{c} $, will lead to a period that scales as $\vert \beta - {\beta }_{c} {\vert }^{- 1/ 2} $ and is independent of particle shape or any geometric aspect ratio in the flow.

Simulations of Droplet Collisions in Shear Flow

2012 ◽  
Author(s):  
Orest Shardt ◽  
J. J. Derksen ◽  
Sushanta K. Mitra
Keyword(s):  
Reynolds Number ◽  
Shear Flow ◽  
Capillary Number ◽  
Three Dimensional ◽  
Reynolds Numbers ◽  
Simple Shear Flow ◽  
Number Range ◽  
Droplet Collisions ◽  
Critical Capillary Number ◽  
Boltzmann Method

When droplets collide in a shear flow, they may coalesce or remain separate after the collision. At low Reynolds numbers, droplets coalesce when the capillary number does not exceed a critical value. We present three-dimensional simulations of droplet coalescence in a simple shear flow. We use a free-energy lattice Boltzmann method (LBM) and study the collision outcome as a function of the Reynolds and capillary numbers. We study the Reynolds number range from 0.2 to 1.4 and capillary numbers between 0.1 and 0.5. We determine the critical capillary number for the simulations (0.19) and find that it is does not depend on the Reynolds number. The simulations are compared with experiments on collisions between confined droplets in shear flow. The critical capillary number in the simulations is about a factor of 25 higher than the experimental value.

On the rotation of a circular porous particle in 2D simple shear flow with fluid inertia

2016 ◽  
Vol 808 ◽  
Author(s):  
Chenggong Li ◽  
Mao Ye ◽  
Zhongmin Liu
Keyword(s):  
Angular Velocity ◽  
Shear Flow ◽  
Simple Shear ◽  
Momentum Equation ◽  
Confinement Effect ◽  
Porous Particle ◽  
Simple Shear Flow ◽  
Rotational Behaviour ◽  
Two Dimensional ◽  
Fluid Inertia

We investigate numerically the rotational behaviour of a circular porous particle suspended in a two-dimensional (2D) simple shear flow with fluid inertia at particle shear Reynolds number up to 108. We use the volume-averaged macroscopic momentum equation to formulate the flow field inside and outside the moving porous particle, which is solved by a modified single relaxation time lattice Boltzmann method. The effects of fluid inertia, confinement of the bounding walls, and permeability of the particle are studied. Our two-dimensional simulation results confirm that the permeability has little effect on the rotation of a porous particle in unbounded shear flow without fluid inertia (Masoud, Stone & Shelley, J. Fluid Mech., vol. 733, 2013, R6), but also suggest that the role of permeability cannot be neglected when the confinement effect is significant, or the fluid inertia is not negligible. The fluid inertia and the confined walls have similar effects on the rotation of a porous particle as that on a solid impermeable particle. The angular velocity decays with an increase in fluid inertia, and the confinement effect suppresses the angular velocity to a shear rate ratio below 0.5. A simple scaling argument based on the balance of torque exerted by fluid flows adjacent to the two bounding walls and that due to the flow recirculation can explain our results.

Three-dimensional calculations of the simple shear flow around a single particle between two moving walls

1995 ◽  
Vol 283 ◽  
pp. 273-285 ◽  
Author(s):  
H. Nirschl ◽  
H. A. Dwyer ◽  
V. Denk
Keyword(s):  
Reynolds Number ◽  
Shear Flow ◽  
Simple Shear ◽  
Particle Surface ◽  
Three Dimensional ◽  
Reynolds Numbers ◽  
Rigid Sphere ◽  
Simple Shear Flow ◽  
Number Range ◽  
Grid Scheme

Three-dimensional solutions have been obtained for the steady simple shear flow over a spherical particle in the intermediate Reynolds number range 0.1 [les ] Re [les ] 100. The shear flow was generated by two walls which move at the same speed but in opposite directions, and the particle was located in the middle of the gap between the walls. The particle-wall interaction is treated by introducing a fully three-dimensional Chimera or overset grid scheme. The Chimera grid scheme allows each component of a flow to be accurately and efficiently treated. For low Reynolds numbers and without any wall influence we have verified the solution of Taylor (1932) for the shear around a rigid sphere. With increasing Reynolds numbers the angular velocity for zero moment for the sphere decreases with increasing Reynolds number. The influence of the wall has been quantified with the global particle surface characteristics such as net torque and Nusselt number. A detailed analysis of the influence of the wall distance and Reynolds number on the surface distributions of pressure, shear stress and heat transfer has also been carried out.

Numerical Simulation of Two-Dimensional Drops Suspended in Simple Shear Flow at Nonzero Reynolds Numbers

2011 ◽  
Vol 133 (3) ◽  
Author(s):  
S. Mortazavi ◽  
Y. Afshar ◽  
H. Abbaspour
Keyword(s):  
Reynolds Number ◽  
Shear Flow ◽  
Normal Stress ◽  
Reynolds Numbers ◽  
Normal Stress Difference ◽  
Two Dimensions ◽  
Simple Shear Flow ◽  
Stress Difference ◽  
Areal Fraction ◽  
Small Reynolds Numbers

The motion of deformable drops suspended in a linear shear flow at nonzero Reynolds numbers is studied by numerical simulations in two dimensions. It is found that a deformable drop migrates toward the center of the channel in agreement with experimental findings at small Reynolds numbers. However, at relatively high Reynolds numbers (Re=80) and small deformation, the drop migrates to an equilibrium position off the centerline. Suspension of drops at a moderate areal fraction (φ=0.44) is studied by simulations of 36 drops. The flow is studied as a function of the Reynolds number and a shear thinning behavior is observed. The results for the normal stress difference show oscillations around a mean value at small Reynolds numbers, and it increases as the Reynolds number is raised. Simulations of drops at high areal fraction (φ=0.66) show that if the Capillary number is kept constant, the effective viscosity does not change in the range of considered Reynolds numbers (0.8–80). The normal stress difference is also a weak function of the Reynolds number. It is also found that similar to flows of granular materials, suspension of drops at finite Reynolds numbers shows the same trend for the density and fluctuation energy distribution across the channel.

Pair-sphere trajectories in finite-Reynolds-number shear flow

2008 ◽  
Vol 596 ◽  
pp. 413-435 ◽  
Author(s):  
PANDURANG M. KULKARNI ◽  
JEFFREY F. MORRIS
Keyword(s):  
Reynolds Number ◽  
Shear Flow ◽  
Fluid Motion ◽  
Reynolds Numbers ◽  
Simple Shear Flow ◽  
Rigid Spheres ◽  
Newtonian Dynamics ◽  
Flow Reynolds Number ◽  
Boltzmann Method ◽  
Neutrally Buoyant

The pair trajectories of neutrally buoyant rigid spheres immersed in finite-inertia simple-shear flow are described. The trajectories are obtained using the lattice-Boltzmann method to solve the fluid motion, with Newtonian dynamics describing the sphere motions. The inertia is characterized by the shear-flow Reynolds number ${\it Re} \,{=}\,\rho\dot{\gamma}a^2/\mu$, where μ and ρ are the viscosity and density of the fluid respectively, $\dot{\gamma}$ is the shear rate and a is the radius of the larger of the pair of spheres in the case of unequal sizes; the majority of results presented are for pairs of equal radii. Reynolds numbers of 0 ≤ Re ≤ 1 are considered with a focus on inertia at Re = O(0.1). At finite inertia, the topology of the pair trajectories is altered from that predicted at Re = 0, as closed trajectories found in Stokes flow vanish and two new forms of trajectories are observed. These include spiralling and reversing trajectories in addition to largely undisturbed open trajectories. For Re = O(0.1), the limits of the various regions in pair space yielding open, reversing and spiralling trajectories are roughly defined.

scholarly journals Reynolds-number dependence of turbulence enhancement on collision growth

2016 ◽  
Vol 16 (19) ◽  
pp. 12441-12455 ◽  
Author(s):  
Ryo Onishi ◽  
Axel Seifert
Keyword(s):  
Reynolds Number ◽  
Isotropic Turbulence ◽  
Reynolds Numbers ◽  
Stokes Number ◽  
Homogeneous Isotropic Turbulence ◽  
Collision Efficiency ◽  
Cloud Droplets ◽  
Clustering Effect ◽  
Reynolds Number Dependence ◽  
Coagulation Kernel

Abstract. This study investigates the Reynolds-number dependence of turbulence enhancement on the collision growth of cloud droplets. The Onishi turbulent coagulation kernel proposed in Onishi et al. (2015) is updated by using the direct numerical simulation (DNS) results for the Taylor-microscale-based Reynolds number (Reλ) up to 1140. The DNS results for particles with a small Stokes number (St) show a consistent Reynolds-number dependence of the so-called clustering effect with the locality theory proposed by Onishi et al. (2015). It is confirmed that the present Onishi kernel is more robust for a wider St range and has better agreement with the Reynolds-number dependence shown by the DNS results. The present Onishi kernel is then compared with the Ayala–Wang kernel (Ayala et al., 2008a; Wang et al., 2008). At low and moderate Reynolds numbers, both kernels show similar values except for r2 ∼ r1, for which the Ayala–Wang kernel shows much larger values due to its large turbulence enhancement on collision efficiency. A large difference is observed for the Reynolds-number dependences between the two kernels. The Ayala–Wang kernel increases for the autoconversion region (r1, r2 < 40 µm) and for the accretion region (r1 < 40 and r2 > 40 µm; r1 > 40 and r2 < 40 µm) as Reλ increases. In contrast, the Onishi kernel decreases for the autoconversion region and increases for the rain–rain self-collection region (r1, r2 > 40 µm). Stochastic collision–coalescence equation (SCE) simulations are also conducted to investigate the turbulence enhancement on particle size evolutions. The SCE with the Ayala–Wang kernel (SCE-Ayala) and that with the present Onishi kernel (SCE-Onishi) are compared with results from the Lagrangian Cloud Simulator (LCS; Onishi et al., 2015), which tracks individual particle motions and size evolutions in homogeneous isotropic turbulence. The SCE-Ayala and SCE-Onishi kernels show consistent results with the LCS results for small Reλ. The two SCE simulations, however, show different Reynolds-number dependences, indicating possible large differences in atmospheric turbulent clouds with large Reλ.

On a Circular Cylinder in a Uniform Planar Shear Flow at Subcritical Reynolds Number

2002 ◽  
Author(s):  
D. Sumner ◽  
O. O. Akosile
Keyword(s):  
Reynolds Number ◽  
Circular Cylinder ◽  
Shear Flow ◽  
Lift Force ◽  
Experimental Studies ◽  
Reynolds Numbers ◽  
Low Velocity ◽  
Static Pressure Distribution ◽  
Planar Shear ◽  
The Mean

An experimental investigation was conducted of a circular cylinder immersed in a uniform planar shear flow, where the approach velocity varies across the diameter of the cylinder. The study was motivated by some apparent discrepancies between numerical and experimental studies of the flow, and the general lack of experimental data, particularly in the subcritical Reynolds number regime. Of interest was the direction and origin of the steady mean lift force experienced by the cylinder, which has been the subject of contradictory results in the literature, and for which measurements have rarely been reported. The circular cylinder was tested at Reynolds numbers from Re = 4.0×104 − 9.0×104, and the dimensionless shear parameter ranged from K = 0.02 − 0.07, which corresponded to a flow with low to moderate shear. The results showed that low to moderate shear has no appreciable influence on the Strouhal number, but has the effect of lowering the mean drag coefficient. The circular cylinder develops a small steady mean lift force directed towards the low-velocity side, which is attributed to an asymmetric mean static pressure distribution on its surface. The reduction in the mean drag force, however, cannot be attributed solely to this asymmetry.

Gravitational and zero-drag motion of a spheroid adjacent to an inclined plane at low Reynolds number

1994 ◽  
Vol 268 ◽  
pp. 267-292 ◽  
Author(s):  
Richard Hsu ◽  
Peter Ganatos
Keyword(s):  
Reynolds Number ◽  
Angular Velocity ◽  
Shear Flow ◽  
Equilibrium Position ◽  
Prolate Spheroid ◽  
Inclined Plane ◽  
Gravitational Settling ◽  
Angular Velocities ◽  
Planar Wall ◽  
Neutrally Buoyant

The first highly accurate solutions for the resistance tensor of an oblate or prolate spheroid moving near a planar wall obtained by Hsu & Ganatos (1989) are used to compute the translational and angular velocities and trajectories of a neutrally buoyant spheroid in shear flow and the gravitational settling motion of a non-neutrally buoyant spheroid adjacent to an inclined plane. The neutrally buoyant spheroid in shear flow undergoes a periodical motion toward and away from the wall as it continually tumbles forward. For some orientation angles it is found that the wall actually enhances the angular velocity of the particle. For certain inclinations a spheroid settling under gravity near an inclined plane reaches an equilibrium position, after which it translates parallel to the wall without rotation.

Sign in / Sign up

Export Citation Format

Share Document