skip to main content


Title: Transport of anisotropic particles under waves
Using a numerical model, we analyse the effects of shape on both the orientation and transport of anisotropic particles in wavy flows. The particles are idealized as prolate and oblate spheroids, and we consider the regime of small Stokes and particle Reynolds numbers. We find that the particles preferentially align into the shear plane with a mean orientation that is solely a function of their aspect ratio. This alignment, however, differs from the Jeffery orbits that occur in the residual shear flow (that is, the Stokes drift velocity field) in the absence of waves. Since the drag on an anisotropic particle depends on its alignment with the flow, this preferred orientation determines the effective drag on the particles, which in turn impacts their net downstream transport. We also find that the rate of alignment of the particles is not constant and depends strongly on their initial orientation; thus, variations in initial particle orientation result in dispersion of anisotropic-particle plumes. We show that this dispersion is a function of the particle’s eccentricity and the ratio of the settling and wave time scales. Due to this preferential alignment, we find that a plume of anisotropic particles in waves is on average transported farther but dispersed less than it would be if the particles were randomly oriented. Our results demonstrate that accurate prediction of the transport of anisotropic particles in wavy environments, such as microplastic particles in the ocean, requires the consideration of these preferential alignment effects.  more » « less
Award ID(s):
1706586
NSF-PAR ID:
10062299
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Journal of Fluid Mechanics
Volume:
837
ISSN:
0022-1120
Page Range / eLocation ID:
320 to 340
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. Particle-laden flows of sedimenting solid particles or droplets in a carrier gas have strong inter-phase coupling. Even at low particle volume fractions, the two-way coupling can be significant due to the large particle to gas density ratio. In this semi-dilute regime, the slip velocity between phases leads to sustained clustering that strongly modulates the overall flow. The analysis of perturbations in homogeneous shear reveals the process by which clusters form: (i) the preferential concentration of inertial particles in the stretching regions of the flow leads to the formation of highly concentrated particle sheets, (ii) the thickness of the latter is controlled by particle-trajectory crossing, which causes a local dispersion of particles, (iii) a transverse Rayleigh–Taylor instability, aided by the shear-induced rotation of the particle sheets towards the gravity normal direction, breaks the planar structure into smaller clusters. Simulations in the Euler–Lagrange formalism are compared to Euler–Euler simulations with the two-fluid and anisotropic-Gaussian methods. It is found that the two-fluid method is unable to capture the particle dispersion due to particle-trajectory crossing and leads instead to the formation of discontinuities. These are removed with the anisotropic-Gaussian method which derives from a kinetic approach with particle-trajectory crossing in mind. 
    more » « less
  3. Particulate matter in the environment, such as sediment, marine debris and plankton, is transported by surface waves. The transport of these inertial particles is different from that of fluid parcels described by Stokes drift. In this study, we consider the transport of negatively buoyant particles that settle in flow induced by surface waves as described by linear wave theory in arbitrary depth. We consider particles that fall under both a linear drag regime in the low Reynolds number limit and in a nonlinear drag regime in the transitional Reynolds number range. Based on an analysis of typical applications, we find that the nonlinear regime is the most widely applicable. From an expansion in the particle Stokes number, we find kinematic expressions for inertial particle motion in waves, and from a multiscale expansion in the dimensionless wave amplitude, we find expressions for the wave-averaged drift velocities. These drift velocities are analogous to Stokes drift and can be used in large-scale models that do not resolve surface waves. We find that the horizontal drift velocity is reduced relative to the Stokes drift of fluid parcels and that the vertical drift velocity is enhanced relative to the particle terminal settling velocity. We also demonstrate that a cloud of settling particles released simultaneously will disperse in the horizontal direction. Finally, we discuss the accuracy of our expressions by comparing against numerical simulations, which show excellent agreement, and against experimental data, which show the same trends. 
    more » « less
  4. Hypothesis: The dip coating of suspensions made of monodisperse non-Brownian spherical particles dispersed in a Newtonian fluid leads to different coating regimes depending on the ratio of the particle diameter to the thickness of the film entrained on the substrate. In particular, dilute particles dispersed in the liquid are entrained only above a threshold value of film thickness. In the case of anisotropic particles, in particular fibers, the smallest characteristic dimension will control the entrainment of the particle. Furthermore, it is possible to control the orientation of the anisotropic particles depending on the substrate geometry. In the thick film regime, the Landau-Levich-Derjaguin model remains valid if one account for the change in viscosity. Experiment: To test the hypotheses, we performed dip-coating experiments with dilute suspensions of non-Brownian fibers with different length-to-diameter aspect ratios. We characterize the number of fibers entrained on the surface of the substrate as a function of the withdrawal velocity, allowing us to estimate a threshold capillary number below which all the particles remain in the liquid bath. Besides, we measure the angular distribution of the entrained fibers for two different substrate geometries: flat plates and cylindrical rods. We then measure the film thickness for more concentrated fiber suspensions. Findings: The entrainment of the fibers on a flat plate and a cylindrical rod is primarily controlled by the smaller characteristic length of the fibers: their diameter. At first order, the entrainment threshold scales similarly to that of spherical particles. The length of the fibers only appears to have a minor influence on the entrainment threshold. No preferential alignment is observed for non-Brownian fibers on a flat plate, except for very thin films, whereas the fibers tend to align themselves along the axis of a cylindrical rod for a large enough ratio of the fiber length to the radius of the cylindrical rod. The Landau-Levich-Derjaguin law is recovered for more concentrated suspension by introducing an effective capillary number accounting for the change in viscosity. 
    more » « less
  5. At low Reynolds numbers, axisymmetric ellipsoidal particles immersed in a shear flow undergo periodic tumbling motions known as Jeffery orbits, with the orbit determined by the initial orientation. Understanding this motion is important for predicting the overall dynamics of a suspension. While slender fibres may follow Jeffery orbits, many such particles in nature are neither straight nor rigid. Recent work exploring the dynamics of curved or elastic fibres have found Jeffery-like behaviour along with chaotic orbits, decaying orbital constants and cross-streamline drift. Most work focuses on particles with reflectional symmetry; we instead consider the behaviour of a composite asymmetric slender body made of two straight rods, suspended in a two-dimensional shear flow, to understand the effects of the shape on the dynamics. We find that for certain geometries the particle does not rotate and undergoes persistent drift across streamlines, the magnitude of which is consistent with other previously identified forms of cross-streamline drift. For this class of particles, such geometry-driven cross-streamline motion may be important in giving rise to dispersion in channel flows, thereby potentially enhancing mixing. 
    more » « less