Effect of fluid inertia on the dynamics and scaling of neutrally buoyant particles in shear flow
2014 (English)In: Journal of Fluid Mechanics, ISSN 0022-1120, E-ISSN 1469-7645, Vol. 738, 563-590 p.Article in journal (Refereed) Published
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)/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 St = alpha . Re-p, where alpha 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 Re-p, 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 vertical bar Re-p - Re-c vertical bar(-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 beta (e. g. confinement or particle spacing) that influences the value of Re-c, such that Re-p = Re-c as beta = beta(c), will lead to a period that scales as vertical bar beta - beta c vertical bar(-1/2) and is independent of particle shape or any geometric aspect ratio in the flow.
Place, publisher, year, edition, pages
2014. Vol. 738, 563-590 p.
complex fluids, multiphase and particle-laden flows, nonlinear dynamical systems
Engineering and Technology
IdentifiersURN: urn:nbn:se:kth:diva-140130DOI: 10.1017/jfm.2013.599ISI: 000328486400024ScopusID: 2-s2.0-84910138134OAI: oai:DiVA.org:kth-140130DiVA: diva2:689665
QC 201401212014-01-212014-01-172016-09-29Bibliographically approved