Thermal convection in mono-disperse and bi-disperse granular gases: A   simulation study
  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. 
Daniela Paolotti, Alain Barrat, Umberto Marini Bettolo Marconi and   Andrea Puglisi