Capturing patterns and symmetries in chaotic granular flow

Segregation patterns in three-dimensional granular flows

Flow of size-bidisperse particle mixtures in a spherical tumbler rotating alternately about two perpendicular axes produces segregation patterns that track the location of nonmixing islands predicted by a dynamical systems approach. To better understand the paradoxical accumulation of large particles in regions defined by barriers to transport, we perform discrete element method (DEM) simulations to visualize the three-dimensional structure of the segregation patterns and track individual particles. Our DEM simulations and modeling results indicate that segregation pattern formation in the biaxial spherical tumbler is due to the interaction of size-driven radial segregation with the weak spanwise component of the advective surface flow. Specifically, we find that after large particles segregate to the surface, slow axial drift in the flowing layer, which is inherent to spherical tumblers, is sufficient to drive large particles across nominal transport barriers and into nonmixing islands predicted by an advective flow model in the absence of axial drift. Axial drift alters the periodic dynamics of nonmixing islands, turning them into "sinks" where large particles accumulate even in the presence of collisional diffusion. Overall, our results indicate that weak perturbation of chaotic flow has the potential to alter key dynamical system features (e.g., transport barriers), which ultimately can result in unexpected physical phenomena.

Particle capture in a model chaotic flow

To better understand and optimize the capture of passive scalars (particles, pollutants, greenhouse gases, etc.) in complex geophysical flows, we study capture in the simpler, but still chaotic, time-dependent double-gyre flow model. For a range of model parameters, the domain of the double-gyre flow consists of a chaotic region, characterized by rapid mixing, interspersed with nonmixing islands in which particle trajectories are regular. Capture units placed within the domain remove all particles that cross their perimeters without altering the velocity field. To predict the capture capability of a unit at an arbitrary location, we characterize the trajectories of a uniformly seeded ensemble of particles as chaotic or nonchaotic, and then use them to determine the spatially resolved fraction of time that the flow is chaotic. With this information, we can predict where to best place units for maximum capture. We also examine the time dependence of the capture process, and demonstrate that there can be a trade-off between the amount of material captured and the capture rate.

Pattern formation in a fully three-dimensional segregating granular flow

Segregation patterns of size-bidisperse particle mixtures in a fully three-dimensional flow produced by alternately rotating a spherical tumbler about two perpendicular axes are studied over a range of particle sizes and volume ratios using both experiments and a continuum model. Pattern formation results from the interaction of size segregation with chaotic regions and nonmixing islands of the flow. Specifically, large particles in the flowing surface layer are preferentially deposited in nonmixing islands despite the effects of collisional diffusion and chaotic transport. The protocol-dependent structure of the unstable manifolds of the flow surrounding the nonmixing islands provides further insight into why certain segregation patterns are more robust than others.

Modeling Segregation in Granular Flows

Accurate continuum models of flow and segregation of dense granular flows are now possible. This is the result of extensive comparisons, over the last several years, of computer simulations of increasing accuracy and scale, experiments, and continuum models, in a variety of flows and for a variety of mixtures. Computer simulations-discrete element methods (DEM)-yield remarkably detailed views of granular flow and segregation. Conti-nuum models, however, offer the best possibility for parametric studies of outcomes in what could be a prohibitively large space resulting from the competition between three distinct driving mechanisms: advection, diffusion, and segregation. We present a continuum transport equation-based framework, informed by phenomenological constitutive equations, that accurately predicts segregation in many settings, both industrial and natural. Three-way comparisons among experiments, DEM, and theory are offered wherever possible to validate the approach. In addition to the flows and mixtures described here, many straightforward extensions of the framework appear possible.

Forced axial segregation in axially inhomogeneous rotating systems

Controlling segregation is both a practical and a theoretical challenge. Using a novel drum design comprising concave and convex geometry, we explore, through the application of both discrete particle simulations and positron emission particle tracking, a means by which radial size segregation may be used to drive axial segregation, resulting in an order of magnitude increase in the rate of separation. The inhomogeneous drum geometry explored also allows the direction of axial segregation within a binary granular bed to be controlled, with a stable, two-band segregation pattern being reliably and reproducibly imposed on the bed for a variety of differing system parameters. This strong banding is observed to persist even in systems that are highly constrained in the axial direction, where such segregation would not normally occur. These findings, and the explanations provided of their underlying mechanisms, could lead to radical new designs for a broad range of particle processing applications but also may potentially prove useful for medical and microflow applications.

Chaotic mixing via streamline jumping in quasi-two-dimensional tumbled granular flows

We study, numerically and analytically, the singular limit of a vanishing flowing layer in tumbled granular flows in quasi-two-dimensional rotating containers. The limiting behavior is found to be identical under the two versions of the kinematic continuum model of such flows, and the transition to the limiting dynamics is analyzed in detail. In particular, we formulate the no-shear-layer dynamical system as a piecewise isometry. It is shown how such a discontinuous map, through the concordant mechanism of streamline jumping, leads to the physical mixing of granular matter. The dependence of the dynamics of Lagrangian particle trajectories on the tumbler fill fraction is also established through Poincaré sections, and, in the special case of a half-full tumbler, chaotic behavior is shown to disappear completely in the singular limit. At other fill levels, stretching in the sense of shear strain is replaced by spreading due to streamline jumping. Finally, we use finite-time Lyapunov exponents to establish the manifold structure and understand "how chaotic" the limiting piecewise isometry is.

Unsteady granular flows in a rotating tumbler

The characteristics of steady granular flow in quasi-two-dimensional rotating tumblers have been thoroughly investigated and are fairly well understood. However, unsteady time-varying flow has not been studied in detail. The linear response of granular flow in quasi-two-dimensional rotating tumblers is presented for periodic forcing protocols via sinusoidal variation in the rotational speed of the tumbler and for step changes in rotational speed. Variations in the tumbler radius, particle size, and forcing frequency are explored. Similarities to steady flow include the fastest flow occurring at the free surface of the flowing layer and an instantaneous approximately linear velocity profile through the depth. The flowing layer depth varies by 2-5 particle diameters between minimum and maximum rotation rates. However, unsteady forcing also causes the flow to exhibit dynamic properties. For periodic rotational speeds, the phase lag of the flowing layer depth increases linearly with increasing input forcing frequency up to nearly 2.0 rad over 0-20 cycles per tumbler revolution. The amplitude responses of the velocity and shear rate show a resonance behavior unique to the system level parameters. The phase lag of all flow properties appears to be related to the number of particle contacts from the edge of the rotating tumbler. Characterization via step changes in rotational speed shows dynamic properties of overshoot (up to 35%) and rise times on the order of 0.2-0.7 s. The results suggest that the unsteady granular flow analysis may be beneficial for characterizing the "flowability" and "rheology" of granular materials based on particle size, moisture content, or other properties.

Granular flow in rotating cylinders with noncircular cross sections

An experimental and theoretical study is carried out of the flow of granular material in cylinders with different cross-sectional shapes rotated about their axes. The flow of particles in such geometries is confined to a shallow layer at the free surface. The length and thickness of the layer shrink and expand periodically with rotation of the cylinder, resulting in chaotic advection and improved mixing of passive tracers. Experimental results obtained by flow visualization are reported for quasi-two-dimensional mixers half filled with glass beads. A depth-averaged flow model to predict the time-varying layer thickness profile is presented, along with a perturbation solution in terms of a small parameter k , which is the ratio of the maximum layer thickness to the half length of the layer (L) , at the cross-section orientation when the length is minimum. To the lowest order [O(k 0)] , the model predicts that the layer profiles scaled with L(theta) at different mixer orientation angles (theta) are identical and the same as that for a circle. The measured layer thickness profiles averaged over different orientations of noncircular mixers match reasonably well with the theory, but the standard deviations are larger for the noncircular cylinders compared to the circle. The O(k) perturbation solution and the full theory both predict that the scaled layer thickness varies periodically; the deviations are proportional to the rate of change of the length with orientation. The perturbation solution gives results close to those from the numerical solution except at cylinder orientations when the length of the flowing layer changes sharply. The measured variation of the scaled midlayer thickness with orientation for all geometries is well predicted by the theory.

Radial granular segregation under chaotic flow in two-dimensional tumblers

An initially well mixed granular material composed of two distinct subclasses of particles, small and large or light and heavy, segregates radially into stable lobed patterns when rotated in various quasi-two-dimensional, regular polygonal tumblers. The patterns are highly sensitive to the time-periodic flow, which in turn depends critically on the fill fraction and container shape. Simulations of a simple model reproduce the segregation patterns observed in experiment. Kolmogorov-Arnol'd-Moser (KAM) regions in Poincaré plots of the velocity field used to model the flow attract smaller (denser) particles and their spatial symmetries mirror those of the segregation patterns, suggesting that competition between the driving forces for radial segregation (percolation and buoyancy) and those for chaotic mixing plays a key role.