Free Cooling of a Granular Gas of Rodlike Particles in Microgravity

Anomalous diffusion in polydisperse granular gases

We investigate both ensemble and time-averaged mean-squared displacements of particles in a polydisperse granular system in a homogeneous cooling state and derive rigorous analytical expressions valid both at short and long time scales. The discrepancies in ensemble and time-averaged mean-squared displacements indicate ergodicity breaking in granular systems consisting of an arbitrary number of species of different sizes and masses. We compare the results of our study with Monte Carlo simulations in terms of a powerful low-rank algorithm and find a nice agreement.

Onsager variational principle for granular fluids

Granular fluids, as defined by a collection of moving solid particles, is a paradigm of a dissipative system out of equilibrium. Inelastic collisions between particles is the source of dissipation, and is the origin of a transition from a gas to a liquidlike state. This transition can be triggered by an increase of the solid fraction. Moreover, in compartmentalized systems, this condensation is driving the entire granular fluid into a Maxwell demon phenomenon, localizing most of the grains into a specific compartment. Classical approaches fail to capture these phenomena, thus motivating many experimental and numerical works. Herein, we demonstrate that the Onsager variational principle is able to predict accurately the coexistence of gas-liquid states in granular systems, opening ways to model other phenomena observed in such dissipative systems like segregation or the jamming transition.

Rotational and translational motions in a homogeneously cooling granular gas

A granular gas composed of monodisperse spherical particles was studied in microgravity experiments in a drop tower. Translations and rotations of the particles were extracted from optical video data. Equipartition is violated, the rotational degrees of freedom were excited only to roughly 2/3 of the translational ones. After stopping the mechanical excitation, we observed granular cooling of the ensemble for a period of three times the Haff time, where the kinetic energy dropped to about 5% of its initial value. The cooling rates of all observable degrees of freedom were comparable, and the ratio of rotational and translational kinetic energies fluctuated around a constant value. The distributions of translational and rotational velocity components showed slight but systematic deviations from Gaussians at the start of cooling.

Cooling of a granular gas mixture in microgravity

Granular gases are fascinating non-equilibrium systems with interesting features such as spontaneous clustering and non-Gaussian velocity distributions. Mixtures of different components represent a much more natural composition than monodisperse ensembles but attracted comparably little attention so far. We present the observation and characterization of a mixture of rod-like particles with different sizes and masses in a drop tower experiment. Kinetic energy decay rates during granular cooling and collision rates were determined and Haff's law for homogeneous granular cooling was confirmed. Thereby, energy equipartition between the mixture components and between individual degrees of freedom is violated. Heavier particles keep a slightly higher average kinetic energy than lighter ones. Experimental results are supported by numerical simulations.

Diffusion of intruders in granular suspensions: Enskog theory and random walk interpretation

The Enskog kinetic theory is applied to compute the mean square displacement of impurities or intruders (modeled as smooth inelastic hard spheres) immersed in a granular gas of smooth inelastic hard spheres (grains). Both species (intruders and grains) are surrounded by an interstitial molecular gas (background) that plays the role of a thermal bath. The influence of the latter on the motion of intruders and grains is modeled via a standard viscous drag force supplemented by a stochastic Langevin-like force proportional to the background temperature. We solve the corresponding Enskog-Lorentz kinetic equation by means of the Chapman-Enskog expansion truncated to first order in the gradient of the intruder number density. The integral equation for the diffusion coefficient is solved by considering the first two Sonine approximations. To test these results, we also compute the diffusion coefficient from the numerical solution of the inelastic Enskog equation by means of the direct simulation Monte Carlo method. We find that the first Sonine approximation generally agrees well with the simulation results, although significant discrepancies arise when the intruders become lighter than the grains. Such discrepancies are largely mitigated by the use of the second Sonine approximation, in excellent agreement with computer simulations even for moderately strong inelasticities and/or dissimilar mass and diameter ratios. We invoke a random walk picture of the intruders' motion to shed light on the physics underlying the intricate dependence of the diffusion coefficient on the main system parameters. This approach, recently employed to study the case of an intruder immersed in a granular gas, also proves useful in the present case of a granular suspension. Finally, we discuss the applicability of our model to real systems in the self-diffusion case. We conclude that collisional effects may strongly impact the diffusion coefficient of the grains.

Enskog kinetic theory of binary granular suspensions: Heat flux and stability analysis of the homogeneous steady state

The Enskog kinetic theory of multicomponent granular suspensions employed previously [Gómez González, Khalil, and Garzó, Phys. Rev. E 101, 012904 (2020)2470-004510.1103/PhysRevE.101.012904] is considered further to determine the four transport coefficients associated with the heat flux. These transport coefficients are obtained by solving the Enskog equation by means of the application of the Chapman-Enskog method around the local version of the homogeneous state. Explicit forms of the heat flux transport coefficients are provided in steady-state conditions by considering the so-called second Sonine approximation to the distribution function of each species. Their quantitative variation on the control parameters of the mixture (masses and diameters, coefficients of restitution, concentration, volume fraction, and the background temperature) is demonstrated and the results show that in general the dependence of the heat flux transport coefficients on inelasticity is clearly different from that found in the absence of the gas phase (dry granular mixtures). As an application of the general results, the stability of the homogeneous steady state is analyzed by solving the linearized Navier-Stokes hydrodynamic equations. The linear stability analysis (which holds for wavelengths long compared with the mean free path) shows that the transversal and longitudinal modes are always stable with respect to long-enough wavelength excitations. This conclusion agrees with previous results derived for monocomponent and (dilute) bidisperse granular suspensions but contrasts with the instabilities found in previous works in dry (no gas phase) granular mixtures.

Granular cooling of ellipsoidal particles in microgravity

A three-dimensional granular gas of ellipsoids is established by exposing the system to the microgravity environment of the International Space Station. We use two methods to measure the dynamics of the constituent particles and report the long-time development of the granular temperature until no further particle movement is detectable. The resulting cooling behavior can be well described by Haff’s cooling law with time scale τ. Different analysis methods show evidence of particle clustering towards the end of the experiment. By using the kinetic theory for ellipsoids we compare the translational energy dissipation of individual collision events with the overall cooling time scale τ. The difference from this comparison indicates how energy is distributed in different degrees of freedom including both translation and rotation during the cooling.

Gravity Governs Shear Localization in Confined Dense Granular Flows

The prediction of flow profiles of slowly sheared granular materials is a major geophysical and industrial challenge. Understanding the role of gravity is particularly important for future planetary exploration in varying gravitational environments. Using the principle of minimization of energy dissipation, and combining experiments and variational analysis, we disentangle the contributions of the gravitational acceleration, confining pressure, and layer thickness on shear strain localization induced by moving fault boundaries at the bottom of a granular layer. The flow profile is independent of the gravity for geometries with a free top surface. However, under a confining pressure or if the sheared layer withstands the weight of the upper layers, increasing gravity promotes the transition from closed shear zones buried in the bulk to open ones that intersect the top surface. We show that the center position and width of the shear zone and the axial angular velocity at the top surface follow universal scaling laws when properly scaled by the gravity, applied pressure, and layer thickness. Our finding that the flow profiles lie on a universal master curve opens the possibility to predict the quasistatic shear flow of granular materials in varying gravitational environments.

Direct observation of crystal nucleation and growth in a quasi-two-dimensional nonvibrating granular system

We study a quasi-two-dimensional macroscopic system of magnetic spherical particles settled on a shallow concave dish under a temporally oscillating magnetic field. The system reaches a stationary state where the energy losses from collisions and friction with the concave dish surface are compensated by the continuous energy input coming from the oscillating magnetic field. Random particle motions show some similarities with the motions of atoms and molecules in a glass or a crystal-forming fluid. Because of the curvature of the surface, particles experience an additional force toward the center of the concave dish. When decreasing the magnetic field, the effective temperature is decreased and diffusive particle motion slows. For slow cooling rates we observe crystallization, where the particles organize into a hexagonal lattice. We study the birth of the crystalline nucleus and the subsequent growth of the crystal. Our observations support nonclassical theories of crystal formation. Initially a dense amorphous aggregate of particles forms, and then in a second stage this aggregate rearranges internally to form the crystalline nucleus. As the aggregate grows, the crystal grows in its interior. After a certain size, all the aggregated particles are part of the crystal and after that crystal growth follows the classical theory for crystal growth.

Visual analysis of density and velocity profiles in dense 3D granular gases

Granular multiparticle ensembles are of interest from fundamental statistical viewpoints as well as for the understanding of collective processes in industry and in nature. Extraction of physical data from optical observations of three-dimensional (3D) granular ensembles poses considerable problems. Particle-based tracking is possible only at low volume fractions, not in clusters. We apply shadow-based and feature-tracking methods to analyze the dynamics of granular gases in a container with vibrating side walls under microgravity. In order to validate the reliability of these optical analysis methods, we perform numerical simulations of ensembles similar to the experiment. The simulation output is graphically rendered to mimic the experimentally obtained images. We validate the output of the optical analysis methods on the basis of this ground truth information. This approach provides insight in two interconnected problems: the confirmation of the accuracy of the simulations and the test of the applicability of the visual analysis. The proposed approach can be used for further investigations of dynamical properties of such media, including the granular Leidenfrost effect, granular cooling, and gas-clustering transitions.

Particle Dynamics at the Onset of the Granular Gas-Liquid Transition

We study experimentally the dynamical behavior of few large tracer particles placed in a quasi-2D granular "gas" made of many small beads in a low-gravity environment. Multiple inelastic collisions transfer momentum from the uniaxially driven gas to the tracers whose velocity distributions are studied through particle tracking. Analyzing these distributions for an increasing system density reveals that translational energy equipartition is reached at the onset of the gas-liquid granular transition corresponding to the emergence of local clusters. The dynamics of a few tracer particles thus appears as a simple and accurate tool to detect this transition. A model is proposed for describing accurately the formation of local heterogeneities.

Cluster dynamics in dense granular gases of rod-like particles

Granular gases are interesting multiparticle systems which, irrespective of the apparent simplicity of particle interactions, exhibit a rich scenario of so far only little understood features. We have numerically investigated a dense granular gas composed of frictional spherocylinders which are excited mechanically by lateral vibrating container walls. This study was stimulated by experiments in microgravity on parabolic flights. The formation of spatial inhomogeneities (clusters) was observed in a region near the corners of the container, about halfway from the excitation plates. The particles in the clusters show a tendency to align parallel to the container walls, seemingly increasing the stabilizing effect of friction. The simulation results provide hints that the phase difference of the vibrations of the two excitation walls might affect the cluster dynamics.

Magnetically heated granular gas in a low-gravity environment

Magnetic forces are used to heat up thousands of spherical particles under low-gravity. This long range external excitation, combined with the induced particle-particle interactions, results in a homogeneous spatial distribution of the particles. Comparisons with predictions of kinetic theories can hence be carried out. Haff’s cooling law is verified qualitatively, while the measured cooling time scale is quantitatively different from the prediction. The high velocity tail of the velocity distribution during homogeneous cooling state (HCS) is measured, while the expected cluster formation after HCS can not be verified by our experiment.

Continuously heated granular gas of elongated particles

Some years ago, Harth et al. experimentally explored the steady state dynamics of a heated granular gas of rod-like particles in microgravity [K. Harth et al. Phys. Rev. Lett. 110, 144102 (2013)]. Here, we report numerical results that quantitatively reproduce their experimental findings and provide additional insight into the process. A system of sphero-cylinders is heated by the vibration of three flat side walls, resulting in one symmetrically heated direction, one non-symmetrically heated direction, and one non-heated direction. In the non-heated direction, the speed distribution follows a stretched exponential distribution $$p(\upsilon )\, \propto \,{\rm{exp}}\left( { - {{\left( {{{\left| \upsilon \right|} \mathord{\left/ {\vphantom {{\left| \upsilon \right|} C}} \right. \kern-\nulldelimiterspace} C}} \right)}^{1.5}}} \right)$$. In the symmetrically heated direction, the velocity statistics at low speeds is similar but it develops pronounced exponential tails at high speeds. In the non-symmetrically heated direction (not accessed experimentally), the distribution also follows $$p(\upsilon )\, \propto \,{\rm{exp}}\left( { - {{\left( {{{\left| \upsilon \right|} \mathord{\left/ {\vphantom {{\left| \upsilon \right|} C}} \right. \kern-\nulldelimiterspace} C}} \right)}^{1.5}}} \right)$$ , but the velocity statistics of rods moving toward the vibrating wall resembles the indirectly excited direction, whereas the velocity statistics of those moving away from the wall resembles the direct excited direction.

Crowding-Enhanced Diffusion: An Exact Theory for Highly Entangled Self-Propelled Stiff Filaments

We study a strongly interacting crowded system of self-propelled stiff filaments by event-driven Brownian dynamics simulations and an analytical theory to elucidate the intricate interplay of crowding and self-propulsion. We find a remarkable increase of the effective diffusivity upon increasing the filament number density by more than one order of magnitude. This counterintuitive "crowded is faster" behavior can be rationalized by extending the concept of a confining tube pioneered by Doi and Edwards for highly entangled, crowded, passive to active systems. We predict a scaling theory for the effective diffusivity as a function of the Péclet number and the filament number density. Subsequently, we show that an exact expression derived for a single self-propelled filament with motility parameters as input can predict the nontrivial spatiotemporal dynamics over the entire range of length and timescales. In particular, our theory captures short-time diffusion, directed swimming motion at intermediate times, and the transition to complete orientational relaxation at long times.

Granular flow of cylinder-like particles in a cylindrical hopper under external pressure based on DEM simulations

Granular flow is widely found in nature or industrial production. Although the external driving force significantly affects the dynamic behavior of a granular system, a large number of numerical simulations have been conducted to study granular flows driven by gravity. In this study, a superquadric equation was used to construct spherical and cylindrical elements, and the flow processes of granular materials under external pressure were simulated by the discrete element method. To examine the validity of the DEM model, the Janssen effect of spherical particles, the static packing of cylindrical particles and the flow process of spherical particles under external pressure are simulated and compared with the previous experimental and theoretical results. Subsequently, the effects of blockiness, orifice diameter, and particle friction on the flow characteristics are investigated. Results show that the flow rate of spherical particles increases as the external pressure and opening diameter increase or the particle friction decreases. However, the flow rate of cylindrical particles decreases as the blockiness parameter increases, and the external pressure has little effect on the flow rate of the cylindrical particles when the blockiness parameter is greater than 4. Furthermore, the external pressure causes a change in the flow pattern of granular systems. In a gravity-driven granular flow, cylindrical particles appear in funnel flow, and spherical particles in both mass and funnel flows. In a pressure-driven granular flow, spherical particles appear in mass flow, and cylindrical particles in both mass and funnel flows. The critical height of the transition between mass and funnel flows decreases with increasing external pressure and eventually reaches a steady state. Meanwhile, the critical height increases with the blockiness parameter, which indicates that more cylindrical than spherical particles appear in funnel flow. Finally, the basic flow characteristics of granular materials under external pressure are further analyzed by the velocity uniformity index, the normal contact force between particles, and the bottom pressure. Overall, the numerical results are useful for understanding the changes in the flow characteristics of spherical and cylindrical granular materials under external pressure, and further provide guidance for the appropriate design and optimization of cylindrical hoppers.

Velocity Distribution of a Homogeneously Cooling Granular Gas

In contrast to molecular gases, granular gases are characterized by inelastic collisions and require therefore permanent driving to maintain a constant kinetic energy. The kinetic theory of granular gases describes how the average velocity of the particles decreases after the driving is shut off. Moreover, it predicts that the rescaled particle velocity distribution will approach a stationary state with overpopulated high-velocity tails as compared to the Maxwell-Boltzmann distribution. While this fundamental theoretical result was reproduced by numerical simulations, an experimental confirmation is still missing. Using a microgravity experiment that allows the spatially homogeneous excitation of spheres via magnetic fields, we confirm the theoretically predicted exponential decay of the tails of the velocity distribution.

Crystallization via shaking in a granular gas with van der Waals interactions

We investigate the effect of van der Waals forces on a collection of granular particles by means of molecular dynamics simulations of a vibrated system in three dimensions. The van der Waals interactions introduce two phase coexistences: one between a random close packing and a gas and a second between a polycrystalline dense state and a gas, where the dense, disordered component crystallizes when the driving amplitude exceeds a threshold value. The region of stability of the ordered state in the nonequilibrium phase diagram grows in size as the Hamaker constant increases or the degree of dissipation increases.

Inertial delay of self-propelled particles

The motion of self-propelled massive particles through a gaseous medium is dominated by inertial effects. Examples include vibrated granulates, activated complex plasmas and flying insects. However, inertia is usually neglected in standard models. Here, we experimentally demonstrate the significance of inertia on macroscopic self-propelled particles. We observe a distinct inertial delay between orientation and velocity of particles, originating from the finite relaxation times in the system. This effect is fully explained by an underdamped generalisation of the Langevin model of active Brownian motion. In stark contrast to passive systems, the inertial delay profoundly influences the long-time dynamics and enables new fundamental strategies for controlling self-propulsion in active matter.

Interplay between polydispersity, inelasticity, and roughness in the freely cooling regime of hard-disk granular gases

A polydisperse granular gas made of inelastic and rough hard disks is considered. Focus is laid on the kinetic-theory derivation of the partial energy production rates and the total cooling rate as functions of the partial densities and temperatures (both translational and rotational) and of the parameters of the mixture (masses, diameters, moments of inertia, and mutual coefficients of normal and tangential restitution). The results are applied to the homogeneous cooling state of the system and the associated nonequipartition of energy among the different components and degrees of freedom. It is found that disks typically present a stronger rotational-translational nonequipartition but a weaker component-component nonequipartition than spheres. A noteworthy "mimicry" effect is unveiled, according to which a polydisperse gas of disks having common values of the coefficient of restitution and of the reduced moment of inertia can be made indistinguishable from a monodisperse gas in what concerns the degree of rotational-translational energy nonequipartition. This effect requires the mass of a disk of component i to be approximately proportional to 2σ_{i}+〈σ〉, where σ_{i} is the diameter of the disk and 〈σ〉 is the mean diameter.

An instrument for studying granular media in low-gravity environment

A new experimental facility has been designed and constructed to study driven granular media in a low-gravity environment. This versatile instrument, fully automatized, with a modular design based on several interchangeable experimental cells, allows us to investigate research topics ranging from dilute to dense regimes of granular media such as granular gas, segregation, convection, sound propagation, jamming, and rheology-all without the disturbance by gravitational stresses active on Earth. Here, we present the main parameters, protocols, and performance characteristics of the instrument. The current scientific objectives are then briefly described and, as a proof of concept, some first selected results obtained in low gravity during parabolic flight campaigns are presented.