Two-dimensional granular gas of inelastic spheres with multiplicative driving

Optimisation of the crystallisation process through staggered cooling in a nonvibrating granular system

We study experimentally the optimisation of the crystallisation process through a 2D-dimensional system of magnetic particles under an oscillating magnetic field. The time-dependent magnetic field fluidises the system, and by varying its magnitude, the effective temperature of the system is controlled. The system exhibits fluid-like behaviour when the effective temperature is high and evolves to a crystalline arrangement when the temperature is slowly lowered in linear cooling. We found that replacing the linear cooling path with a staggered one let us find the conditions to crystallise the system quicker than using linear cooling. We determine the minimum time required for the particles to find their minimum energy configuration at each temperature step. We found that the crystallisation time was considerably reduced using this method, which allowed us to optimise the crystallisation time.

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.

Anomalous flocking in nonpolar granular Brownian vibrators

Using Brownian vibrators, we investigated the structures and dynamics of quasi-2d granular materials, with packing fractions (ϕ) ranging from 0.111 to 0.832. Our observations revealed a remarkable large-scale flocking behavior in hard granular disk systems, encompassing four distinct phases: granular fluid, flocking fluid, poly-crystal, and crystal. Anomalous flocking emerges at ϕ = 0.317, coinciding with a peak in local density fluctuations, and ceased at ϕ = 0.713 as the system transitioned into a poly-crystal state. The poly-crystal and crystal phases resembled equilibrium hard disks, while the granular and flocking fluids differed significantly from equilibrium systems and previous experiments involving uniformly driven spheres. This disparity suggests that collective motion arises from a competition controlled by volume fraction, involving an active force and an effective attractive interaction resulting from inelastic particle collisions. Remarkably, these findings align with recent theoretical research on the flocking motion of spherical active particles without alignment mechanisms.

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.

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.

Glass- and crystal-forming model based on a granular two-dimensional system

We study a two-dimensional system of magnetic particles under an alternating magnetic field. The particles are settled on the surface of a negative lens where they tend to sediment toward the center due to gravity. The effective temperature is controlled by the intensity of the applied magnetic field. The system is cooled down from a gaslike state to a solidlike state at different rates. We observe that for some slow cooling rates the final configuration of system is a hexagonal compact arrange, while for the faster cooling rates the final configurations are glasslike states. We followed the time evolution of the system, which allows us to determine in detail changes in quantities such as the interparticle distance. We determine the glass transition temperature for different cooling rates, finding that such temperature increases as the cooling rate decreases, in contrast with some other glass-forming liquids.

Temperature distribution in driven granular mixtures does not depend on mechanism of energy dissipation

We study analytically and numerically the distribution of granular temperatures in granular mixtures for different dissipation mechanisms of inelastic inter-particle collisions. Both driven and force-free systems are analyzed. We demonstrate that the simplified model of a constant restitution coefficient fails to predict even qualitatively a granular temperature distribution in a homogeneous cooling state. At the same time we reveal for driven systems a stunning result - the distribution of temperatures in granular mixtures is universal. That is, it does not depend on a particular dissipation mechanism of inter-particles collisions, provided the size distributions of particles is steep enough. The results of the analytic theory are compared with simulation results obtained by the direct simulation Monte Carlo (DSMC). The agreement between the theory and simulations is perfect. The reported results may have important consequences for fundamental science as well as for numerous application, e.g. for the experimental modelling in a lab of natural processes.

Magnetic excitation of a granular gas as a bulk thermostat

A thermostat utilizing a varying magnetic field has been developed to agitate soft ferromagnetic particles in microgravity platforms for an investigation of an energy-dissipative granular gas. Although the method has experimentally realized a reasonably homogeneous spatial distribution of particles, the physics behind the magnetically excited particles has not been understood. Therefore, a numerical calculation based on the discrete element method is developed in this paper to explain the realization of homogeneously distributed particles. The calculation method allows considering inelastic and magnetic interactions between particles and tracking the motions due to those interactions during the excitation of the granular gas. The calculation results, compared with the experimental result, show that magnetic interactions between particles, a time-domain variation of magnetic-excitation directions, and random collisions of particles between each magnetic excitation contribute to distribute particles homogeneously.

Hydrodynamic granular segregation induced by boundary heating and shear

Segregation induced by a thermal gradient of an impurity in a driven low-density granular gas is studied. The system is enclosed between two parallel walls from which we input thermal energy to the gas. We study here steady states occurring when the inelastic cooling is exactly balanced by some external energy input (stochastic force or viscous heating), resulting in a uniform heat flux. A segregation criterion based on Navier-Stokes granular hydrodynamics is written in terms of the tracer diffusion transport coefficients, whose dependence on the parameters of the system (masses, sizes, and coefficients of restitution) is explicitly determined from a solution of the inelastic Boltzmann equation. The theoretical predictions are validated by means of Monte Carlo and molecular dynamics simulations, showing that Navier-Stokes hydrodynamics produces accurate segregation criteria even under strong shearing and/or inelasticity.

Linear hydrodynamics for driven granular gases

We study the dynamics of a granular gas heated by a stochastic thermostat. From a Boltzmann description, we derive the hydrodynamic equations for small perturbations around the stationary state that is reached in the long time limit. Transport coefficients are identified as Green-Kubo formulas obtaining explicit expressions as a function of the inelasticity and the spatial dimension.

Universal reference state in a driven homogeneous granular gas

We study the dynamics of a homogeneous granular gas heated by a stochastic thermostat, in the low density limit. It is found that, before reaching the stationary regime, the system quickly "forgets" the initial condition and then evolves through a universal state that does not only depend on the dimensionless velocity, but also on the instantaneous temperature, suitably renormalized by its steady state value. We find excellent agreement between the theoretical predictions at the Boltzmann equation level for the one-particle distribution function and the direct Monte Carlo simulations. We conclude that at variance with the homogeneous cooling phenomenology, the velocity statistics should not be envisioned as a single-parameter, but as a two-parameter scaling form, keeping track of the distance to stationarity.

Effect of Coulomb friction on orientational correlation and velocity distribution functions in a sheared dilute granular gas

From particle simulations of a sheared frictional granular gas, we show that the Coulomb friction can have dramatic effects on orientational correlation as well as on both the translational and angular velocity distribution functions even in the Boltzmann (dilute) limit. The dependence of orientational correlation on friction coefficient (μ) is found to be nonmonotonic, and the Coulomb friction plays a dual role of enhancing or diminishing the orientational correlation, depending on the value of the tangential restitution coefficient (which characterizes the roughness of particles). From the sticking limit (i.e., with no sliding contact) of rough particles, decreasing the Coulomb friction is found to reduce the density and spatial velocity correlations which, together with diminished orientational correlation for small enough μ, are responsible for the transition from non-gaussian to gaussian distribution functions in the double limit of small friction (μ→0) and nearly elastic particles (e→1). This double limit in fact corresponds to perfectly smooth particles, and hence the maxwellian (gaussian) is indeed a solution of the Boltzmann equation for a frictional granular gas in the limit of elastic collisions and zero Coulomb friction at any roughness. The high-velocity tails of both distribution functions seem to follow stretched exponentials even in the presence of Coulomb friction, and the related velocity exponents deviate strongly from a gaussian with increasing friction.

Model for density waves in gravity-driven granular flow in narrow pipes

A gravity-driven flow of grains through a narrow pipe in vacuum is studied by means of a one-dimensional model with two coefficients of restitution. Numerical simulations show clearly how density waves form when a strikingly simple criterion is fulfilled: that dissipation due to collisions between the grains and the walls of the pipe is greater per collision than that which stems from collisions between particles. Counterintuitively, the highest flow rate is observed when the number of grains per density wave grows large. We find strong indication that the number of grains per density wave always approaches a constant as the particle number tends to infinity, and that collapse to a single wave, which was often observed also in previous simulations, occurs because the number of grains is insufficient for multiple wave formation.

Constitutive relations in dense granular flows

We use simulations to investigate constitutive relations in dry granular flow. Our system is comprised of polydisperse sets of spherical grains falling down a vertical chute under the influence of gravity. Three phases or states of granular matter are observed: a free-fall dilute granular gas region at the top of the chute, a granular fluid in the middle and then a glassy region at the bottom. We examine a complete closed set of constitutive relations capable of describing the local stresses, heat flow, and dissipation in the different regions. While the pressure can be reasonably described by hard sphere gas models, the transport coefficients cannot. Transport coefficients such as viscosity and heat conductivity increase with decreasing temperature in the fluid and glassy phases. The glass exhibits signs of a finite yield stress and we show that the static sand pile is a limit of our glassy state.

Velocity fluctuations in dense granular flows

We use simulations to investigate velocity fluctuations in dry granular flow. Our system is comprised of mono- and polydisperse sets of spherical grains falling down a vertical chute under the influence of gravity. We find three different classes of velocity distributions depending on factors such as the local density. The class of the velocity distribution depends on whether the grains are in a free-fall, fluid, or glassy state. The analytic form of the distributions match those that have been found by other authors in fairly diverse systems. Here, we have all three present in a single system in steady state. Power-law tails that match recent experiments are also found but in a transition area suggesting they may be an artifact of crossover from one class of velocity distribution to another. We find evidence that the transition from one class to another may correspond to a second order dynamical phase transition in the limit that the vertical flow speed goes to zero.

Orientational correlation and velocity distributions in uniform shear flow of a dilute granular gas

Using particle simulations of the uniform shear flow of a rough dilute granular gas, we show that the translational and rotational velocities are strongly correlated in direction, but there is no orientational correlation-induced singularity at perfectly smooth (beta=-1) and rough (beta=1) limits for elastic collisions (e=1); both the translational and rotational velocity distribution functions remain close to a Gaussian for these two limiting cases. Away from these two limits, the orientational as well as spatial velocity correlations are responsible for the emergence of non-Gaussian high-velocity tails. The tails of both distribution functions follow stretched exponentials, with the exponents depending on normal (e) and tangential (beta) restitution coefficients.

Boltzmann equation for dissipative gases in homogeneous states with nonlinear friction

Combining analytical and numerical methods, we study within the framework of the homogeneous nonlinear Boltzmann equation a broad class of models relevant for the dynamics of dissipative fluids, including granular gases. We use the method presented in a previous paper [J. Stat. Phys. 124, 549 (2006)] and extend our results to a different heating mechanism: namely, a deterministic nonlinear friction force. We derive analytically the high-energy tail of the velocity distribution and compare the theoretical predictions with high-precision numerical simulations. Stretched exponential forms are obtained when the nonequilibrium steady state is stable. We derive subleading corrections and emphasize their relevance. In marginal stability cases, power-law behaviors arise, with exponents obtained as the roots of transcendental equations. We also consider some simple Bhatnagar-Gross-Krook models, driven by similar heating devices, to test the robustness of our predictions.

Velocity distribution and the effect of wall roughness in granular Poiseuille flow

From event-driven simulations of a gravity-driven channel flow of inelastic hard disks, we show that the velocity distribution function remains close to a Gaussian for a wide range densities (even when the Knudsen number is of order 1) if the walls are smooth and the particle collisions are nearly elastic. For dense flows, a transition from a Gaussian to a power-law distribution for the high-velocity tails occurs with increasing dissipation in the center of the channel, irrespective of wall roughness. For a rough wall, the near-wall distribution functions are distinctly different from those in the bulk, even in the quasielastic limit.

Local heat flux and energy loss in a two-dimensional vibrated granular gas

We performed event-driven simulations of a two-dimensional granular gas between two vibrating walls and directly measured the local heat flux and local energy dissipation in the stationary state. Describing the local heat flux as a function of the coordinate in the direction perpendicular to the driving walls, we test a generalization of Fourier's law, q(x)=-kappa inverted delta T(x)+mu inverted delta rho(x), by relating the local heat flux to the local gradients of the temperature and density. This ansatz accounts for the fact that heat flux can also be generated by density gradients, not only by temperature gradients. Assuming the transport coefficients kappa and mu to be independent of x, we check the validity of this assumption and test the generalized Fourier law in the simulations. Both kappa and mu are determined for different system parameters, in particular, for a wide range of coefficients of restitution. We also compare our numerical results to existing hydrodynamic theories. Agreement is found for kappa for very small inelasticities only, i.e., when the gradients are small. Beyond this region, kappa and mu exhibit a striking nonmonotonic behavior. This may hint that hydrodynamics to Navier-Stokes order cannot be applied to moderately inelastic vibrated systems.

Segregation mechanisms in a numerical model of a binary granular mixture

A simple phenomenological numerical model of a binary granular mixture is developed and investigated numerically. We attempt to model a recently reported experimental system where a horizontally vibrated binary monolayer was found to exhibit a transition from a mixed to a segregated state as the filling fraction of the mixture was increased. This numerical model is found to reproduce much of the experimentally observed behavior, most importantly the transition from the mixed to the segregated state. We use the numerical model to investigate granular segregation mechanisms and explain the experimentally observed behavior.

Local equation of state and velocity distributions of a driven granular gas

We present event-driven simulations of a granular gas of inelastic hard disks with incomplete normal restitution in two dimensions between vibrating walls (without gravity). We measure hydrodynamic quantities such as the stress tensor, density and temperature profiles, as well as velocity distributions. Relating the local pressure to the local temperature and local density, we construct a local constitutive equation. For strong inelasticities the local constitutive relation depends on global system parameters, like the volume fraction and the aspect ratio. For moderate inelasticities the constitutive relation is approximately independent of the system parameters and can hence be regarded as a local equation of state, even though the system is highly inhomogeneous with heterogeneous temperature and density profiles arising as a consequence of energy injection. With respect to local velocity distributions we find that they do not scale with the square root of the local granular temperature. Moreover the high-velocity tails are different for the distribution of the x and the y components of the velocity, and even depend on the position in the sample, the global volume fraction, and the coefficient of restitution.

Cluster growth in two- and three-dimensional granular gases

Dissipation in granular media leads to interesting phenomena such as cluster formation and crystallization in nonequilibrium dynamical states. The freely cooling system is examined concerning the energy decay and the cluster evolution in time, both in two and three dimensions. We also suggest an interpretation of the three-dimensional cluster growth in terms of percolation theory, but this point deserves further study.

Granular fluid thermostated by a bath of elastic hard spheres

The homogeneous steady state of a fluid of inelastic hard spheres immersed in a bath of elastic hard spheres kept at equilibrium is analyzed by means of the first Sonine approximation to the (spatially homogeneous) Enskog-Boltzmann equation. The temperature of the granular fluid relative to the bath temperature and the kurtosis of the granular distribution function are obtained as functions of the coefficient of restitution, the mass ratio, and a dimensionless parameter beta measuring the cooling rate relative to the friction constant. Comparison with recent results obtained from an iterative numerical solution of the Enskog-Boltzmann equation [Biben et al., Physica A 310, 308 (2002)] shows an excellent agreement. Several limiting cases are also considered. In particular, when the granular particles are much heavier than the bath particles (but have a comparable size and number density), it is shown that the bath acts as a white noise external driving. In the general case, the Sonine approximation predicts the lack of a steady state if the control parameter beta is larger than a certain critical value beta(c) that depends on the coefficient of restitution and the mass ratio. However, this phenomenon appears outside the expected domain of applicability of the approximation.

Power-law distribution of pressure fluctuations in multiphase flow

Bubbling fluidized beds are granular systems, in which a deep layer of particles is set in motion by a vertical gas stream, with the excess gas rising as bubbles through the bed. We show that pressure fluctuations in such a system have non-Gaussian statistics. The probability density function has a power-law drop-off and is very well represented by a Tsallis distribution. Its shape is explained through the folding of the Gaussian distribution of pressure fluctuations produced by a monodisperse set of bubbles, onto the actual distribution of bubble sizes in the bed, assuming that bubbles coalesce via a Smoluchowski-type aggregation process. Therefore, the Tsallis statistics arise as a result of bubble polydispersity, rather than system nonextensivity.

Computing a non-Maxwellian velocity distribution from first principles

We investigate a family of single-particle anomalous velocity distribution by solving a particular class of stochastic Liouville equations. The stationary state is obtained analytically and the Maxwell-Boltzmann distribution is reobtained in a particular limit. We discuss the comparison with other different methods to obtain the stationary state. Extensions when the models cannot be solved in an exact way are also pointed out in connection with the one-ficton approximation.

Fluctuation-dissipation relations in driven granular gases

We study the dynamics of a two-dimensional driven inelastic gas, by means of direct simulation Monte Carlo techniques, i.e., under the assumption of molecular chaos. Under the effect of a uniform stochastic driving in the form of a white noise plus a friction term, the gas is kept in a nonequilibrium steady state characterized by fractal density correlations and non-Gaussian distributions of velocities; the mean-squared velocity, that is the so-called granular temperature, is lower than the bath temperature. We observe that a modified form of the Kubo relation, which relates the autocorrelation and the linear response for the dynamics of a system at equilibrium, still holds for the off equilibrium, though stationary, dynamics of the systems under investigation. Interestingly, the only needed modification to the equilibrium Kubo relation is the replacement of the equilibrium temperature with an effective temperature, which results equal to the global granular temperature. We present two independent numerical experiments, where two different observables are studied: (a) the staggered density current, whose response to an impulsive shear is proportional to its autocorrelation in the unperturbed system and (b) the response of a tracer to a small constant force, switched on at time t(w), which is proportional to the mean-square displacement in the unperturbed system. Both measures confirm the validity of Kubo's formula, provided that the granular temperature is used as the proportionality factor between response and autocorrelation, at least for not too large inelasticities.

Randomly driven granular fluids: collisional statistics and short scale structure

We present a molecular-dynamics and kinetic theory study of granular material, modeled by inelastic hard disks, fluidized by a random driving force. The focus is on collisional averages and short-distance correlations in the nonequilibrium steady state, in order to analyze in a quantitative manner the breakdown of molecular chaos, i.e., factorization of the two-particle distribution function, f((2))(x(1),x(2)) approximately chif((1))(x(1))f((1))(x(2)) in a product of single-particle ones, where x(i)=[r(i),v(i)] with i=1,2 and chi represents the position correlation. We have found that molecular chaos is only violated in a small region of the two-particle phase space [x(1),x(2)], where there is a predominance of grazing collisions. The size of this singular region grows with increasing inelasticity. The existence of particle- and noise-induced recollisions magnifies the departure from mean-field behavior. The implications of this breakdown in several physical quantities are explored.