Thermal convection in mono-disperse and bi-disperse granular gases: A simulation study

TL;DR

This simulation study investigates thermal convection in granular gases, revealing how parameters like inelasticity, gravity, and mixture composition influence the onset and nature of convection patterns.

Contribution

The paper introduces a new Event-Driven Molecular Dynamics simulation method to analyze convection in mono- and bi-disperse granular gases, aligning results with hydrodynamical predictions.

Findings

Increasing inelasticity causes a transition from no convection to one and then two convective rolls.

Higher gravity can suppress convection, leading to no flow.

Mixture properties and impurities significantly affect convection patterns.

Abstract

We present results of a simulation study of inelastic hard-disks vibrated in a vertical container. An Event-Driven Molecular Dynamics method is developed for studying the onset of convection. Varying the relevant parameters (inelasticity, number of layers at rest, intensity of the gravity) we are able to obtain a qualitative agreement of our results with recent hydrodynamical predictions. Increasing the inelasticity, a first continuous transition from the absence of convection to one convective roll is observed, followed by a discontinuous transition to two convective rolls, with hysteretic behavior. At fixed inelasticity and increasing gravity, a transition from no convection to one roll can be evidenced. If the gravity is further increased, the roll is eventually suppressed. Increasing the number of monolayers the system eventually localizes mostly at the bottom of the box: in this case multiple convective rolls as well as surface waves appear. We analyze the density and temperature fields and study the existence of symmetry breaking in these fields in the direction perpendicular to the injection of energy. We also study a binary mixture of grains with different properties (inelasticity or diameters). The effect of changing the properties of one of the components is analyzed, together with density, temperature and temperature ratio fields. Finally, the presence of a low-fraction of quasi-elastic impurities is shown to determine a sharp transition between convective and non-convective steady states.