Role of roughness on the hydrodynamic homogeneous base state of inelastic spheres

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.

Translational and rotational non-Gaussianities in homogeneous freely evolving granular gases

The importance of roughness in the modeling of granular gases has been increasingly considered in recent years. In this paper, a freely evolving homogeneous granular gas of inelastic and rough hard disks or spheres is studied under the assumptions of the Boltzmann kinetic equation. The homogeneous cooling state is studied from a theoretical point of view using a Sonine approximation, in contrast to a previous Maxwellian approach. A general theoretical description is done in terms of d_{t} translational and d_{r} rotational degrees of freedom, which accounts for the cases of spheres (d_{t}=d_{r}=3) and disks (d_{t}=2, d_{r}=1) within a unified framework. The non-Gaussianities of the velocity distribution function of this state are determined by means of the first nontrivial cumulants and by the derivation of non-Maxwellian high-velocity tails. The results are validated by computer simulations using direct simulation Monte Carlo and event-driven molecular dynamics algorithms.

Kinetic Theory and Memory Effects of Homogeneous Inelastic Granular Gases under Nonlinear Drag

We study a dilute granular gas immersed in a thermal bath made of smaller particles with masses not much smaller than the granular ones in this work. Granular particles are assumed to have inelastic and hard interactions, losing energy in collisions as accounted by a constant coefficient of normal restitution. The interaction with the thermal bath is modeled by a nonlinear drag force plus a white-noise stochastic force. The kinetic theory for this system is described by an Enskog-Fokker-Planck equation for the one-particle velocity distribution function. To get explicit results of the temperature aging and steady states, Maxwellian and first Sonine approximations are developed. The latter takes into account the coupling of the excess kurtosis with the temperature. Theoretical predictions are compared with direct simulation Monte Carlo and event-driven molecular dynamics simulations. While good results for the granular temperature are obtained from the Maxwellian approximation, a much better agreement, especially as inelasticity and drag nonlinearity increase, is found when using the first Sonine approximation. The latter approximation is, additionally, crucial to account for memory effects such as Mpemba and Kovacs-like ones.

Optimal Control of Uniformly Heated Granular Fluids in Linear Response

We present a detailed analytical investigation of the optimal control of uniformly heated granular gases in the linear regime. The intensity of the stochastic driving is therefore assumed to be bounded between two values that are close, which limits the possible values of the granular temperature to a correspondingly small interval. Specifically, we are interested in minimising the connection time between the non-equilibrium steady states (NESSs) for two different values of the granular temperature by controlling the time dependence of the driving intensity. The closeness of the initial and target NESSs make it possible to linearise the evolution equations and rigorously-from a mathematical point of view-prove that the optimal controls are of bang-bang type, with only one switching in the first Sonine approximation. We also look into the dependence of the optimal connection time on the bounds of the driving intensity. Moreover, the limits of validity of the linear regime are investigated.

Power statistics of Otto heat engines with the Mpemba effect

The Mpemba effect is a counterintuitive relaxation phenomenon whereby a system with a higher initial temperature may cool down to the thermal state faster than an identical system that was initially prepared at a lower temperature. Here, we investigate heat and work in a Markovian state transition system with cyclic switching hot-cold temperatures, which operates as an Otto heat engine working in long but finite time, either with or without the Mpemba effect. Under the condition of the periodic steady state having been reached, the time durations of the heating and cooling relaxation processes are determined by exploring a distance-from-equilibrium equivalent to the Kullback-Leibler divergence. We then numerically evaluate and compare the averages and variances of both the work and the power output of two scenarios with and without the Mpemba effect. The results show that the Markovian Mpemba effect can enhance the machine performance by significantly increasing the power output for a given efficiency without sacrificing the stability.

Hydrodynamics of granular gases of inelastic and rough hard disks or spheres. II. Stability analysis

Conditions for the stability under linear perturbations around the homogeneous cooling state are studied for dilute granular gases of inelastic and rough hard disks or spheres with constant coefficients of normal (α) and tangential (β) restitution. After a formally exact linear stability analysis of the Navier-Stokes-Fourier hydrodynamic equations in terms of the translational (d_{t}) and rotational (d_{r}) degrees of freedom, the transport coefficients derived in the companion paper [A. Megías and A. Santos, "Hydrodynamics of granular gases of inelastic and rough hard disks or spheres. I. Transport coefficients" Phys. Rev. E 104, 034901 (2021)10.1103/PhysRevE.104.034901] are employed. Known results for hard spheres [Garzó, Santos, and Kremer, Phys. Rev. E 97, 052901 (2018)10.1103/PhysRevE.97.052901] are recovered by setting d_{t}=d_{r}=3, while novel results for hard disks (d_{t}=2, d_{r}=1) are obtained. In the latter case, a high-inelasticity peculiar region in the (α,β) parameter space is found, inside which the critical wave number associated with the longitudinal modes diverges. Comparison with event-driven molecular dynamics simulations for dilute systems of hard disks at α=0.2 shows that this theoretical region of absolute instability may be an artifact of the extrapolation to high inelasticity of the approximations made in the derivation of the transport coefficients, although it signals a shrinking of the conditions for stability. In the case of moderate inelasticity (α=0.7), however, a good agreement between the theoretical predictions and the simulation results is found.

Hydrodynamics of granular gases of inelastic and rough hard disks or spheres. I. Transport coefficients

The transport coefficients for dilute granular gases of inelastic and rough hard disks or spheres with constant coefficients of normal (α) and tangential (β) restitution are obtained in a unified framework as functions of the number of translational (d_{t}) and rotational (d_{r}) degrees of freedom. The derivation is carried out by means of the Chapman-Enskog method with a Sonine-like approximation in which, in contrast to previous approaches, the reference distribution function for angular velocities does not need to be specified. The well-known case of purely smooth d-dimensional particles is recovered by setting d_{t}=d and formally taking the limit d_{r}→0. In addition, previous results [G. M. Kremer, A. Santos, and V. Garzó, Phys. Rev. E 90, 022205 (2014)10.1103/PhysRevE.90.022205] for hard spheres are reobtained by taking d_{t}=d_{r}=3, while novel results for hard-disk gases are derived with the choice d_{t}=2, d_{r}=1. The singular quasismooth limit (β→-1) and the conservative Pidduck's gas (α=β=1) are also obtained and discussed.

Linear response in the uniformly heated granular gas

We analyze the linear response properties of the uniformly heated granular gas. The intensity of the stochastic driving fixes the value of the granular temperature in the nonequilibrium steady state reached by the system. Here, we investigate two specific situations. First, we look into the "direct" relaxation of the system after a single (small) jump of the driving intensity. This study is carried out by two different methods. Not only do we linearize the evolution equations around the steady state, but we also derive generalized out-of-equilibrium fluctuation-dissipation relations for the relevant response functions. Second, we investigate the behavior of the system in a more complex experiment, specifically a Kovacs-like protocol with two jumps in the driving. The emergence of an anomalous Kovacs response is explained in terms of the properties of the direct relaxation function: it is the second mode changing sign at the critical value of the inelasticity that demarcates anomalous from normal behavior. The analytical results are compared with numerical simulations of the kinetic equation, and a good agreement is found.

Pseudo-two-dimensional dynamics in a system of macroscopic rolling spheres

We study in this work the dynamics of a collection of identical hollow spheres (ping-pong balls) that rest on a horizontal metallic grid. Fluidization is achieved by means of a turbulent air current coming from below. The upflow is adjusted so that the particles do not levitate over the grid, resulting in quasi-two-dimensional dynamics. We show that the behavior of diffusion and correlations in this system is particularly rich. Noticeably as well (and related to the complex dynamical behavior), a variety of phases appear, with important peculiarities with respect to analogous setups. We observe gas, liquid, glass, and hexagonal crystal phases. Most notably, we show that the melting of the hexagonal crystal occurs in coexistence with a liquid phase. This strikingly differs from the corresponding transition in a purely two-dimensional systems of air-fluidized disks, for which no phase coexistence has been reported in the literature.

Large Mpemba-like effect in a gas of inelastic rough hard spheres

We report the emergence of a giant Mpemba effect in the uniformly heated gas of inelastic rough hard spheres: The initially hotter sample may cool sooner than the colder one, even when the initial temperatures differ by more than one order of magnitude. In order to understand this behavior, it suffices to consider the simplest Maxwellian approximation for the velocity distribution in a kinetic approach. The largeness of the effect stems from the fact that the rotational and translational temperatures, which obey two coupled evolution equations, are comparable. Our theoretical predictions agree very well with molecular dynamics and direct simulation Monte Carlo data.

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.

Incorporation of velocity-dependent restitution coefficient and particle surface friction into kinetic theory for modeling granular flow cooling

Kinetic theory (KT) has been successfully used to model rapid granular flows in which particle interactions are frictionless and near elastic. However, it fails when particle interactions become frictional and inelastic. For example, the KT is not able to accurately predict the free cooling process of a vibrated granular medium that consists of inelastic frictional particles under microgravity. The main reason that the classical KT fails to model these flows is due to its inability to account for the particle surface friction and its inelastic behavior, which are the two most important factors that need be considered in modeling collisional granular flows. In this study, we have modified the KT model that is able to incorporate these two factors. The inelasticity of a particle is considered by establishing a velocity-dependent expression for the restitution coefficient based on many experimental studies found in the literature, and the particle friction effect is included by using a tangential restitution coefficient that is related to the particle friction coefficient. Theoretical predictions of the free cooling process by the classical KT and the improved KT are compared with the experimental results from a study conducted on an airplane undergoing parabolic flights without the influence of gravity [Y. Grasselli, G. Bossis, and G. Goutallier, Europhys. Lett. 86, 60007 (2009)10.1209/0295-5075/86/60007]. Our results show that both the velocity-dependent restitution coefficient and the particle surface friction are important in predicting the free cooling process of granular flows; the modified KT model that integrates these two factors is able to improve the simulation results and leads to better agreement with the experimental results.

Energy nonequipartition in gas mixtures of inelastic rough hard spheres: The tracer limit

The dynamical properties of a tracer or impurity particle immersed in a host gas of inelastic and rough hard spheres in the homogeneous cooling state is studied. Specifically, the breakdown of energy equipartition as characterized by the tracer/host ratios of translational and rotational temperatures is analyzed by exploring a wide spectrum of values of the control parameters of the system (masses, moments of inertia, sizes, and coefficients of restitution). Three complementary approaches are considered. On the theoretical side, the Boltzmann and Boltzmann-Lorentz equations (both assuming the molecular chaos ansatz) are solved by means of a multitemperature Maxwellian approximation for the velocity distribution functions. This allows us to obtain explicit analytical expressions for the temperature ratios. On the computational side, two different techniques are used. First, the kinetic equations are numerically solved by the direct simulation Monte Carlo (DSMC) method. Second, molecular dynamics simulations for dilute gases are performed. Comparison between theory and simulations shows a general good agreement. This means that (i) the impact of the molecular chaos ansatz on the temperature ratios is not significant (except at high inelasticities and/or big impurities) and (ii) the simple Maxwellian approximation yields quite reliable predictions.

Transport coefficients of a granular gas of inelastic rough hard spheres

The Boltzmann equation for inelastic and rough hard spheres is considered as a model of a dilute granular gas. In this model, the collisions are characterized by constant coefficients of normal and tangential restitution, and hence the translational and rotational degrees of freedom are coupled. A normal solution to the Boltzmann equation is obtained by means of the Chapman-Enskog method for states near the homogeneous cooling state. The analysis is carried out to first order in the spatial gradients of the number density, the flow velocity, and the granular temperature. The constitutive equations for the momentum and heat fluxes and for the cooling rate are derived, and the associated transport coefficients are expressed in terms of the solutions of linear integral equations. For practical purposes, a first Sonine approximation is used to obtain explicit expressions of the transport coefficients as nonlinear functions of both coefficients of restitution and the moment of inertia. Known results for purely smooth inelastic spheres and perfectly elastic and rough spheres are recovered in the appropriate limits.

Higher-order effects on orientational correlation and relaxation dynamics in homogeneous cooling of a rough granular gas

The orientational or angular correlation between the directions of the translational and rotational motions is analyzed theoretically for the homogeneous cooling state of a rough granular gas. The dynamical equations are derived using an approximate form of the single-particle distribution function that incorporates angular correlations. The goal is to assess the effects of higher-order angular corrections for which both quadratic- and quartic-order terms (in translational and rotational velocities of particles) are retained in the perturbation expansion of the distribution function. We show that higher-order corrections can markedly affect the steady-state orientational correlation when the normal restitution coefficient is moderate or small, and this effect is more prominent for nearly smooth particles. The transient evolution of orientational correlation is found to be significantly affected by higher-order terms. In particular the higher-order orientational correlations can dominate over the leading-order contribution during short times even in the quasi-elastic limit, although the steady correlation remains unaffected by such corrections in the same limit. The buildup of correlations during the transient stage seems to be closely tied to the evolution of the ratio between the rotational and translational temperatures. It is demonstrated that the transient dynamics of the temperature ratio and its steady state remain insensitive to higher-order angular correlation.

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.