Modern inertial microfluidics routinely employs oscillatory flows around localized solid features or microbubbles for controlled, specific manipulation of particles, droplets, and cells. It is shown that theories of inertial effects that have been state of the art for decades miss major contributions and strongly underestimate forces on small suspended objects in a range of practically relevant conditions. An analytical approach is presented that derives a complete set of inertial forces and quantifies them in closed form as easy-to-use equations of motion, spanning the entire range from viscous to inviscid flows. The theory predicts additional attractive contributions toward oscillating boundaries, even for density-matched particles, a previously unexplained experimental observation. The accuracy of the theory is demonstrated against full-scale, three-dimensional direct numerical simulations throughout its range.
more »
« less
Density-contrast induced inertial forces on particles in oscillatory flows
Oscillatory flows have become an indispensable tool in microfluidics, inducing inertial effects for displacing and manipulating fluid-borne objects in a reliable, controllable and label-free fashion. However, the quantitative description of such effects has been confined to limit cases and specialized scenarios. Here we develop an analytical formalism yielding the equation of motion of density-mismatched spherical particles in oscillatory background flows, generalizing previous work. Inertial force terms are systematically derived from the geometry of the flow field together with analytically known Stokes number dependences. Supported by independent, first-principles direct numerical simulations, we find that these forces are important even for nearly density-matched objects such as cells or bacteria, enabling their fast displacement and separation. Our formalism thus consistently incorporates particle inertia into the Maxey–Riley equation, and in doing so provides a generalization of Auton's modification to added mass, as well as recovering the description of acoustic radiation forces on particles as a limiting case.
more »
« less
- Award ID(s):
- 1846752
- PAR ID:
- 10586102
- Publisher / Repository:
- Journal of Fluid Mechanics
- Date Published:
- Journal Name:
- Journal of Fluid Mechanics
- Volume:
- 985
- ISSN:
- 0022-1120
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
The drag force on a spherical intruder in dense granular shear flows is studied using discrete element method simulations. Three regimes of the intruder dynamics are observed depending on the magnitude of the drag force (or the corresponding intruder velocity) and the flow inertial number: a fluctuation-dominated regime for small drag forces; a viscous regime for intermediate drag forces; and an inertial (cavity formation) regime for large drag forces. The transition from the viscous regime (linear force-velocity relation) to the inertial regime (quadratic force-velocity relation) depends further on the inertial number. Despite these distinct intruder dynamics, we find a quantitative similarity between the intruder drag in granular shear flows and the Stokesian drag on a sphere in a viscous fluid for intruder Reynolds numbers spanning five orders of magnitude. Beyond this first-order description, a modified Stokes drag model is developed that accounts for the secondary dependence of the drag coefficient on the inertial number and the intruder size and density ratios. When the drag model is coupled with a segregation force model for intruders in dense granular flows, it is possible to predict the velocity of gravity-driven segregation of an intruder particle in shear flow simulations.more » « less
-
The Stokes equation describes the motion of fluids when inertial forces are negligible compared with viscous forces. In this article, we explore the consequence of parity-violating and non-dissipative (i.e. odd) viscosities on Stokes flows in three dimensions. Parity-violating viscosities are coefficients of the viscosity tensor that are not invariant under mirror reflections of space, while odd viscosities are those which do not contribute to dissipation of mechanical energy. These viscosities can occur in systems ranging from synthetic and biological active fluids to magnetized and rotating fluids. We first systematically enumerate all possible parity-violating viscosities compatible with cylindrical symmetry, highlighting their connection to potential microscopic realizations. Then, using a combination of analytical and numerical methods, we analyse the effects of parity-violating viscosities on the Stokeslet solution, on the flow past a sphere or a bubble and on many-particle sedimentation. In all the cases that we analyse, parity-violating viscosities give rise to an azimuthal flow even when the driving force is parallel to the axis of cylindrical symmetry. For a few sedimenting particles, the azimuthal flow bends the trajectories compared with a traditional Stokes flow. For a cloud of particles, the azimuthal flow impedes the transformation of the spherical cloud into a torus and the subsequent breakup into smaller parts that would otherwise occur. The presence of azimuthal flows in cylindrically symmetric systems (sphere, bubble, cloud of particles) can serve as a probe for parity-violating viscosities in experimental systems.more » « less
-
The oscillatory dynamics of nanoelectromechanical systems (NEMS) is at the heart of many emerging applications in nanotechnology. For common NEMS, such as beams and strings, the oscillatory dynamics is formulated using a dissipationless wave equation derived from elasticity. Under a harmonic ansatz, the wave equation gives an undamped free vibration equation; solving this equation with the proper boundary conditions provides the undamped eigenfunctions with the familiar standing wave patterns. Any harmonically driven solution is expressible in terms of these undamped eigenfunctions. Here, we show that this formalism becomes inconvenient as dissipation increases. To this end, we experimentally map out the position- and frequency-dependent oscillatory motion of a NEMS string resonator driven linearly by a non-symmetric force at one end at different dissipation limits. At low dissipation (high Q factor), we observe sharp resonances with standing wave patterns that closely match the eigenfunctions of an undamped string. With a slight increase in dissipation, the standing wave patterns become lost, and waves begin to propagate along the nanostructure. At large dissipation (low Q factor), these propagating waves become strongly attenuated and display little, if any, resemblance to the undamped string eigenfunctions. A more efficient and intuitive description of the oscillatory dynamics of a NEMS resonator can be obtained by superposition of waves propagating along the nanostructure.more » « less
-
Understanding the orientation dynamics of anisotropic colloidal particles is important for suspension rheology and particle self-assembly. However, even for the simplest case of dilute suspensions in shear flow, the orientation dynamics of non-spherical Brownian particles are poorly understood. Here we analytically calculate the time-dependent orientation distributions for non-spherical axisymmetric particles confined to rotate in the flow–gradient plane, in the limit of small but non-zero Brownian diffusivity. For continuous shear, despite the complicated dynamics arising from the particle rotations, we find a coordinate change that maps the orientation dynamics to a diffusion equation with a remarkably simple ratio of the enhanced rotary diffusivity to the zero shear diffusion: $$D_{eff}^{r}/D_{0}^{r}=(3/8)(p-1/p)^{2}+1$$ , where $$p$$ is the particle aspect ratio. For oscillatory shear, the enhanced diffusion becomes orientation dependent and drastically alters the long-time orientation distributions. We describe a general method for solving the time-dependent oscillatory shear distributions and finding the effective diffusion constant. As an illustration, we use this method to solve for the diffusion and distributions in the case of triangle-wave oscillatory shear and find that they depend strongly on the strain amplitude and particle aspect ratio. These results provide new insight into the time-dependent rheology of suspensions of anisotropic particles. For continuous shear, we find two distinct diffusive time scales in the rheology that scale separately with aspect ratio $$p$$ , as $$1/D_{0}^{r}p^{4}$$ and as $$1/D_{0}^{r}p^{2}$$ for $$p\gg 1$$ . For oscillatory shear flows, the intrinsic viscosity oscillates with the strain amplitude. Finally, we show the relevance of our results to real suspensions in which particles can rotate freely. Collectively, the interplay between shear-induced rotations and diffusion has rich structure and strong effects: for a particle with aspect ratio 10, the oscillatory shear intrinsic viscosity varies by a factor of $${\approx}2$$ and the rotational diffusion by a factor of $${\approx}40$$ .more » « less
An official website of the United States government

