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

Författare: Rosén, T., Lundell, F., Aidun, C. K.
Dokumenttyp: Artikel
Tillstånd: Publicerad
Tidskrift: Journal of Fluid Mechanics
Volym: 738   563-590
År: 2013


The 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 (Re p) and Stokes number (St), estimating the effects of fluid and particle inertia, respectively, compared with viscous forces on the particle. The particle Reynolds number is defined by Re_p=4Ga^2/ν, where G is the shear rate, a is the length of the spheroid major semi-axis and ν is the kinematic viscosity. The Stokes number is defined as St=α⋅Re p, where α is the solid-to-fluid density ratio. Here, a neutrally buoyant prolate spheroidal particle (St=Re p) of aspect ratio (major axis/minor axis) r p=4 is considered. The long-term rotational motion for different initial orientations and 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 Rep, the motion is constrained to this planar motion regardless of initial orientation. At a certain critical Reynolds number, Re p=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 |Re p-Re c| -1/2. Analyses in this paper show that if a transition from tumbling to steady state occurs at Re p=Re c, then any parameter β (e.g. confinement or particle spacing) that influences the value of Re c, such that Re p=Re c as β=β c, will lead to a period that scales as |β-β c| -1/2 and is independent of particle shape or any geometric aspect ratio in the flow.