Single-particle motion in three-dimensional vibrofluidized granular beds

Application of Positron Emission Particle Tracking (PEPT) for the evaluation of powder behaviour in an incline linear blender for Continuous Direct Compression (CDC)

Positron Emission Particle Tracking (PEPT) is a non-invasive measurement technique which offers the ability to track the motion of individual particles with high temporal and spatial resolution, and thus build up an understanding of the bulk behaviour of a system from its microscopic (particle level) dynamics. Using this measurement technique, we have developed a series of novel metrics to better understand the behaviours of powders during the steady-state operation of a continuous blender system. Results are presented concerning the response of particle motion to processing parameters (mixing blade configuration and RPM), quantifying the motion in terms of predicted mixing performance. It was found that both increasing rpm and increasing hold-up mass (by selecting fewer transport blades and more mixing blades) provided improved mixing conditions. Interestingly, under specific conditions, there is evidence of convection-like mixing occurring at the interface of the transport and mixing region. This suggests the existence of a potential 'folding region' whereby powder is transported up the barrel (and away from the powder bulk bed) before being reconstituted back into the bulk mass. The results also provide valuable experimental data for the development, calibration and validation of future Discrete Element Method (DEM) simulations.

Coffee bean particle motion in a rotating drum measured using Positron Emission Particle Tracking (PEPT)

Physicochemical transformation of coffee during roasting depends on the applied time-temperature profile (i.e., rate of heat transfer), with heat transfer phenomena governed by particle dynamics. Positron Emission Particle Tracking (PEPT), a non-invasive imaging technique, was used here to characterise the granular flow of coffee in a real, pilot-scale rotating drum roaster. The experimental study established the impact of drum speed, batch size and bean density (i.e., roast degree) on the system's particle dynamics. Particle motion data revealed two distinct regions: (i) a disperse (low occupancy, high velocity) region of in-flight particles and (ii) a dense (high occupancy, low velocity) bean bed. Implications of these results for heat transfer suggest that controlling drum speed for different density coffees will provide roaster operators with a tool to modulate conductive heat transfer from the heated drum to the bean bed. These comprehensive data thus inform roasting best practices and support the development of physics-driven models coupling heat and mass transfer to particle dynamics.

Recent advances in positron emission particle tracking: a comparative review

Positron emission particle tracking (PEPT) is a technique which allows the high-resolution, three-dimensional imaging of particulate and multiphase systems, including systems which are large, dense, and/or optically opaque, and thus difficult to study using other methodologies. In this work, we bring together researchers from the world's foremost PEPT facilities not only to give a balanced and detailed overview and review of the technique but, for the first time, provide a rigorous, direct, quantitative assessment of the relative strengths and weaknesses of all contemporary PEPT methodologies. We provide detailed explanations of the methodologies explored, including also interactive code examples allowing the reader to actively explore, edit and apply the algorithms discussed. The suite of benchmarking tests performed and described within the document is made available in an open-source repository for future researchers.

Positron Emission Particle Tracking of Granular Flows

Positron emission particle tracking (PEPT) is a noninvasive technique capable of imaging the three-dimensional dynamics of a wide variety of powders, particles, grains, and/or fluids. The PEPT technique can track the motion of particles with high temporal and spatial resolution and can be used to study various phenomena in systems spanning a broad range of scales, geometries, and physical states. We provide an introduction to the PEPT technique, an overview of its fundamental principles and operation, and a brief review of its application to a diverse range of scientific and industrial systems.

New Insight into Pseudo-Thermal Convection in Vibrofluidised Granular Systems

Utilising a combination of experimental results obtained via positron emission particle tracking (PEPT) and numerical simulations, we study the influence of a system’s geometric and elastic properties on the convective behaviours of a dilute, vibrofluidised granular assembly. Through the use of a novel, ‘modular’ system geometry, we demonstrate the existence of several previously undocumented convection-inducing mechanisms and compare their relative strengths across a broad, multi-dimensional parameter space, providing criteria through which the dominant mechanism within a given system – and hence its expected dynamics – may be predicted. We demonstrate a range of manners through which the manipulation of a system’s geometry, material properties and imposed motion may be exploited in order to induce, suppress, strengthen, weaken or even invert granular convection. The sum of our results demonstrates that boundary-layer effects due to wall (in)elasticity or directional impulses due to ‘rough’ boundaries exert only a secondary influence on the system’s behaviour. Rather, the direction and strength of convective motion is predominantly determined by the energy flux in the vicinity of the system’s lateral boundaries, demonstrating unequivocally that pseudo-thermal granular convection is decidedly a collective phenomenon.

Real-time probing of granular dynamics with magnetic resonance

Granular dynamics govern earthquakes, avalanches, and landslides and are of fundamental importance in a variety of industries ranging from energy to pharmaceuticals to agriculture. Nonetheless, our understanding of the underlying physics is poor because we lack spatially and temporally resolved experimental measurements of internal grain motion. We introduce a magnetic resonance imaging methodology that provides internal granular velocity measurements that are four orders of magnitude faster compared to previous work. The technique is based on a concerted interplay of scan acceleration and materials engineering. Real-time probing of granular dynamics is explored in single- and two-phase systems, providing fresh insight into bubble dynamics and the propagation of shock waves upon impact of an intruder. We anticipate that the methodology outlined here will enable advances in understanding the propagation of seismic activity, the jamming transition, or the rheology and dynamics of dense suspensions.

Dynamic X-ray radiography reveals particle size and shape orientation fields during granular flow

When granular materials flow, the constituent particles segregate by size and align by shape. The impacts of these changes in fabric on the flow itself are not well understood, and thus novel non-invasive means are needed to observe the interior of the material. Here, we propose a new experimental technique using dynamic X-ray radiography to make such measurements possible. The technique is based on Fourier transformation to extract spatiotemporal fields of internal particle size and shape orientation distributions during flow, in addition to complementary measurements of velocity fields through image correlation. We show X-ray radiography captures the bulk flow properties, in contrast to optical methods which typically measure flow within boundary layers, as these are adjacent to any walls. Our results reveal the rich dynamic alignment of particles with respect to streamlines in the bulk during silo discharge, the understanding of which is critical to preventing destructive instabilities and undesirable clogging. The ideas developed in this paper are directly applicable to many other open questions in granular and soft matter systems, such as the evolution of size and shape distributions in foams and biological materials.

Modifying self-assembly and species separation in three-dimensional systems of shape-anisotropic particles

The behaviors of large, dynamic assemblies of macroscopic particles are of direct relevance to geophysical and industrial processes and may also be used as easily studied analogs to micro- or nano-scale systems, or model systems for microbiological, zoological, and even anthropological phenomena. We study vibrated mixtures of elongated particles, demonstrating that the inclusion of differing particle "species" may profoundly alter a system's dynamics and physical structure in various diverse manners. The phase behavior observed suggests that our system, despite its athermal nature, obeys a minimum free energy principle analogous to that observed for thermodynamic systems. We demonstrate that systems of exclusively spherical objects, which form the basis of numerous theoretical frameworks in many scientific disciplines, represent only a narrow region of a wide, multidimensional phase space. Thus, our results raise significant questions as to whether such models can accurately describe the behaviors of systems outside this highly specialized case.

Dynamic self-assembly of non-Brownian spheres studied by molecular dynamics simulations

Granular self-assembly of confined non-Brownian spheres under gravity is studied by molecular dynamics simulations. Starting from a disordered phase, dry or cohesive spheres organize, by vibrational annealing, into body-centered-tetragonal or face-centered-cubic structures, respectively. During the self-assembling process, isothermal and isodense points are observed. The existence of such points indicates that both granular temperature and packing fraction undergo an inversion process that may be in the core of crystal nucleation. Around the isothermal point, a sudden growth of granular clusters having the maximum coordination number takes place, indicating the outcome of a first-order phase transition. We propose a heuristic equation that successfully describes the dynamic evolution of the local packing fraction in terms of the local granular temperature, along the entire crystallization process.

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.

Maximizing energy transfer in vibrofluidized granular systems

Using discrete particle simulations validated by experimental data acquired using the positron emission particle tracking technique, we study the efficiency of energy transfer from a vibrating wall to a system of discrete, macroscopic particles. We demonstrate that even for a fixed input energy from the wall, energy conveyed to the granular system under excitation may vary significantly dependent on the frequency and amplitude of the driving oscillations. We investigate the manner in which the efficiency with which energy is transferred to the system depends on the system variables and determine the key control parameters governing the optimization of this energy transfer. A mechanism capable of explaining our results is proposed, and the implications of our findings in the research field of granular dynamics as well as their possible utilization in industrial applications are discussed.

Competition between geometrically induced and density-driven segregation mechanisms in vibrofluidized granular systems

The behaviors of granular systems are sensitive to a wide variety of particle properties, including size, density, elasticity, and shape. Differences in any of these properties between particles in a granular mixture may lead to segregation, or "demixing," a process of great industrial relevance. Despite the known influence of particle geometry in granular systems, a considerable fraction of research into these systems concerns only uniformly spherical particles. We address, for the case of vertically vibrated granular systems, the important question of whether the introduction of differing particle geometries entirely invalidates our existing knowledge based on purely spherical granulates, or whether current models may simply be adapted to account for the effects of particle shape. We demonstrate that while shape effects can indeed influence the dynamical and segregative behaviors of a granular system, the segregative mechanisms associated with particle geometry are decidedly secondary to those related to particle density. The relevant control parameters determining the extent of geometrically induced segregation are established. Finally, a manner in which shape effects may be accounted for in simulations utilizing purely spherical particles is proposed.

Low-frequency oscillations and convective phenomena in a density-inverted vibrofluidized granular system

Low-frequency oscillations (LFOs) are thought to play an important role in the transition between the Leidenfrost and convective states of a vibrated granular bed. This work details the experimental observation of LFOs, which are found to be consistently present for a range of driving frequencies and amplitudes, with particles of varying material and using containers of differing material properties. The experimentally acquired results show a close qualitative and quantitative agreement with both theory and simulations across the range of parameters tested. Strong agreement between experimental and simulation results was also observed when investigating the influence of sidewall dissipation on LFOs and vertical density profiles. This paper additionally provides evidence of two phenomena present in the Leidenfrost state: a circulatory motion over extended time periods in near-crystalline configurations, and a Leidenfrost-like state in which the dense upper region displays an unusual inverse thermal convection.

Center of mass scaling in three-dimensional binary granular systems

Using a combination of experimental results acquired through positron emission particle tracking and simulational results obtained via the discrete particle method, we determine a scaling relationship for the center of mass height of a vibrofluidized three-dimensional, bidisperse granular system. We find the scaling to be dependent on the characteristic velocity with which the system is driven, the depth of the granular bed, and the elasticities of the particles involved, as well as the degree of segregation exhibited by the system and the ratio of masses between particle species. The scaling is observed to be robust over a significant range of system parameters.

Dynamics of a first-order transition to an absorbing state

A granular system confined in a quasi-two-dimensional box that is vertically vibrated can transit to an absorbing state in which all particles bounce vertically in phase with the box, with no horizontal motion. In principle, this state can be reached for any density lower than the one corresponding to one complete monolayer, which is then the critical density. Below this critical value, the transition to the absorbing state is of first order, with long metastable periods, followed by rapid transitions driven by homogeneous nucleation. Molecular dynamics simulations and experiments show that there is a dramatic increase on the metastable times far below the critical density; in practice, it is impossible to observe spontaneous transitions close to the critical density. This peculiar feature is a consequence of the nonequilibrium nature of this first-order transition to the absorbing state. A Ginzburg-Landau model, with multiplicative noise, describes qualitatively the observed phenomena and explains the macroscopic size of the critical nuclei. The nuclei become of small size only close to a second critical point where the active phase becomes unstable via a saddle node bifurcation. It is only close to this second critical point that experiments and simulations can evidence spontaneous transitions to the absorbing state while the metastable times grow dramatically moving away from it.

Effects of packing density on the segregative behaviors of granular systems

We present results concerning the important role of system packing in the processes of density- and inelasticity-induced segregation in vibrofluidized binary granular beds. Data are acquired through a combination of experimental results acquired from positron emission particle tracking and simulations performed using the discrete particle method. It is found that segregation due to inelasticity differences between particle species is most pronounced in moderately dense systems, yet still exerts a significant effect in all but the highest density systems. Results concerning segregation due to disparities in particles' material densities show that the maximal degree to which a system can achieve segregation is directly related to the density of the system, while the rate at which segregation occurs shows an inverse relation. Based on this observation, a method of minimizing the time and energy requirements associated with producing a fully segregated system is proposed.

Energy non-equipartition in strongly convective granular systems

Using positron emission particle tracking, the effects of convective motion and the resulting segregative behaviour on the partition of kinetic energy between the components of a bidisperse granular system are, for the first time, systematically investigated. It is found that the distribution of energy between the two system components, which are equal in size but differ in their material properties, is strongly dependent on the degree of segregation observed in the granular bed. The results obtained demonstrate that the difference in energy obtained by dissimilar particle species is not an innate property of the materials in question, but can in fact be altered through variation of the relevant system parameters. The existence of a relationship between the convective and segregative properties of a granular system and the degree of energy equipartition within the system implies the possibility of extending existing theory into the convective regime. Thus, our findings represent an incremental step towards the definition of a granular analogy to temperature that can be applied to more generalised systems and, through this, an improved understanding of inhomogeneous granular systems in general.

Influence of thermal convection on density segregation in a vibrated binary granular system

Using a combination of experimental results and discrete particle method simulations, the role of buoyancy-driven convection in the segregative behavior of a three-dimensional, binary granular system is investigated. A relationship between convective motion and segregation intensity is presented, and a qualitative explanation for this behavior is proposed. This study also provides an insight into the role of diffusive behavior in the segregation of a granular bed in the convective regime. The results of this work strongly imply the possibility that, for an adequately fluidized granular bed, the degree of segregation may be indirectly controlled through the adjustment of the system's driving parameters, or the dissipative properties of the system's side-boundaries.

Thermal convection and temperature inhomogeneity in a vibrofluidized granular bed: the influence of sidewall dissipation

Using a vertically vibrated, fully three-dimensional granular system, we investigate the impact of dissipative interactions between the particles in the system and the vertical sidewalls bounding it. We find that sidewall dissipation influences various properties of the bed including, but not limited to, the spatial distribution of granular temperatures, the functional form of velocity distributions, and the strength of convection. Simple, monotonic relationships are observed for all the aforementioned properties, including a striking linear relationship between convection strength and wall dissipation. We conclude that sidewall effects are not limited to the vicinity of the walls themselves, but extend into the bulk of the system and hence must be considered even in relatively wide, three-dimensional systems. We also propose the possibility of using the alteration of sidewall material as a method of "tuning" certain system parameters in situations where changing the bulk properties or driving parameters of a granular system may be undesirable.

Boltzmann statistics in a three-dimensional vibrofluidized granular bed: idealizing the experimental system

The importance of the method of energy injection into a vertically vibrated granular bed was investigated using positron emission particle tracking. A comparison was made between two experimental systems. The first was driven by a flat base, the second by a base comprising a monolayer of spheres, each free to move in three dimensions but constrained to the bottom of the system, effectively randomizing energy transfer into the bed. The latter system exhibited more Gaussian velocity distributions, a more isotropic temperature, and a better approximation of molecular chaos than the former, behavior more akin to the idealized situation of white noise forcing than is generally observed in an experimental system.

Three-dimensional simulations of a vertically vibrated granular bed including interstitial air

We present a numerical study of the effect of interstitial air on a vertically vibrated granular bed within one period of oscillation. We use a three-dimensional molecular-dynamics simulation including air phenomenologically. The simulations are validated with experiments made with spherical glass beads in a rectangular container. After validation, results are reported for a granular column of 9000 grains and approximately 50 layers deep (at rest), agitated with a sinusoidal excitation with maximal acceleration 4.7g at 11.7 Hz. We report the evolution of density, granular temperature, and coordination number within a vibration cycle, and the effect of interstitial air on those parameters. In three-dimensional computer simulations we found that the presence of interstitial air can promote the collective motion of the granular material as a whole.

NMR measurements and hydrodynamic simulations of phase-resolved velocity distributions within a three-dimensional vibrofluidized granular bed

We report the results of nuclear magnetic resonance imaging experiments on vertically vibrated granular beds of mustard grains. A novel spin-echo velocity profiling technique was developed that allows granular temperature, mean velocity and packing fraction distributions within the three-dimensional cell to be measured as a function of both vertical position and vibration phase. Bimodal velocity distributions were observed at certain portions of the vibration cycle, and in general the ability to acquire time-resolved data demonstrated the significant distortions to the velocity distributions and the systematic errors in calculated temperature distributions that may arise with time-averaged measurements. The experimental behaviour was compared with predictions from a time-varying one-dimensional hydrodynamic model using the experimental parameters as input to the code. In both cases, damping of longitudinal sound waves was linked to significant volume heating effects, which contrasts with the usual heat transport mechanism (i.e. diffusion from the boundaries) currently assumed in most steady-state models. This leads to a new explanation for the counterintuitive upturn in granular temperature in vibrofluidized granular beds, based on amplification and damping of sound waves in the high-altitude region.

Fluidization of a vertically vibrated two-dimensional hard sphere packing: a granular meltdown

We report measurements of the fluidization process in vertically vibrated two-dimensional granular packings. An initially close packed granular bed is exposed to sinusoidal container oscillations with gradually increasing amplitude. At first the particles close to the free surface become mobile. When a critical value of the forcing strength is reached the remaining crystal suddenly breaks up and the bed fluidizes completely. This transition leads to discontinuous changes in the density distribution and in the root mean square displacement of the individual particles. Likewise the vertical center of mass coordinate increases by leaps and bounds at the transition. It turns out that the maximum container velocity v0 is the crucial driving parameter determining the state of a fully fluidized system. For particles of various sizes the transition to full fluidization occurs at the same value of v 2 0/gd, where d is the particle diameter and g is the gravitational acceleration. A discontinuous fluidization transition is only observed when the particles are highly elastic.

Shock wave propagation in vibrofluidized granular materials

Shock wave formation and propagation in vertically vibrated quasi-two-dimensional granular materials are studied by digital high speed photography. Steep density and temperature wave fronts form at the bottom of the granular layer when the layer collides with vibrating plate. Then the fronts propagate upwards through the layer. The temperature front is always in the transition region between the upward and downward granular flows. The effects of driving parameters and particle number on the shock are also explored.

Convection in vibrated annular granular beds

The response to vibration of a granular bed, consisting of a standard cylindrical geometry but with the addition of a dissipative cylindrical inner wall, has been investigated both experimentally (using positron emission particle tracking) and numerically (using hard sphere molecular dynamics simulation). The packing fraction profiles and granular temperature distributions (in both vertical and horizontal directions) were determined as a function of height and distance from the axis. The two sets of results were in reasonable agreement. The molecular dynamics simulations were used to explore the behavior of the granular bed in the inner wall-outer wall coefficient of restitution phase space. It was observed that one could control the direction of the toroidal convection rolls by manipulating the relative dissipation at the inner and outer walls via the coefficients of restitution, and with several layers of grains it was seen that double convection rolls could also be formed, a result that was subsequently confirmed experimentally.

Kink-induced transport and segregation in oscillated granular layers

We use experiments and molecular dynamics simulations of vertically oscillated granular layers to study horizontal particle segregation induced by a kink (a boundary between domains oscillating out of phase). Counterrotating convection rolls carry the larger particles in a bidisperse layer along the granular surface to a kink, where they become trapped. The convection originates from avalanches that occur inside the layer, along the interface between solidified and fluidized grains. The position of a kink can be controlled by modulation of the container frequency, making possible systematic harvesting of the larger particles.