Very fine sand is prepared in a well-defined and fully decompactified state by letting gas bubble through it. After turning off the gas stream, a steel ball is dropped on the sand. On impact of the ball, sand is blown away in all directions ("splash") and an impact crater forms. When this cavity collapses, a granular jet emerges and is driven straight into the air. A second jet goes downwards into the air bubble entrained during the process, thus pushing surface material deep into the ground. The air bubble rises slowly towards the surface, causing a granular eruption. In addition to the experiments and the discrete particle simulations we present a simple continuum theory to account for the void collapse leading to the formation of the upward and downward jets.

The stability of two neighboring surface nanobubbles on a chemically heterogeneous surface is studied by molecular dynamics (MD) simulations of binary mixtures consisting of Lennard-Jones (LJ) particles. A diffusion equation-based stability analysis suggests that two nanobubbles sitting next to each other remain stable, provided the contact line is pinned, and that their radii of curvature are equal. However, many experimental observations seem to suggest some long-term kind of ripening or shrinking of the surface nanobubbles. In our MD simulations we find that the growth/dissolution of the nanobubbles can occur due to the transfer of gas particles from one nanobubble to another along the solid substrate. That is, if the interaction between the gas and the solid is strong enough, the solid–liquid interface can allow for the existence of a “tunnel” which connects the liquid–gas interfaces of the two nanobubbles to destabilize the system. The crucial role of the gas–solid interaction energy is a nanoscopic element that hitherto has not been considered in any macroscopic theory of surface nanobubbles and may help to explain experimental observations of the long-term ripening.

Buoyancy driven granular convection is studied for a shallow, vertically shaken granular bed in a quasi 2D container. Starting from the granular Leidenfrost state, in which a dense particle cluster floats on top of a dilute gaseous layer of fast particles (Meerson et al. in Phys Rev Lett 91:024301, 2003; Eshuis et al. in Phys Rev Lett 95:258001, 2005), we witness the emergence of counter-rotating convection rolls when the shaking strength is increased above a critical level. This resembles the classical onset of convection-at a critical value of the Rayleigh number-in a fluid heated from below. The same transition, even quantitatively, is seen in molecular dynamics simulations, and explained by a hydrodynamic-like model in which the granular material is treated as a continuum. The critical shaking strength for the onset of granular convection is accurately reproduced

The dissolution of a multicomponent nanodrop in a sparingly miscible liquid is studied by molecular dynamics (MD) simulations. We studied both binary and ternary systems, in which nanodroplets are formed from one and two components, respectively. Whereas for a single-component droplet the dissolution can easily be calculated, the situation is more complicated for a multicomponent drop, as the interface concentrations of the drop constituents depend on the drop composition, which changes with time. In this study, the variation of the interface concentration with the drop composition is determined from independent 'numerical experiments', which are then used in the theoretical model for the dissolution dynamics of a multicomponent drop. The MD simulations reveal that when the interaction strengths between the drop constituents and the surrounding bulk liquid are significantly different, the concentration of the more soluble component near the drop interface may become larger than in the drop bulk. This effect is the larger the smaller the drop radius. While the present study is limited to binary and ternary systems, the same method can be easily extended to a larger number of components.

Discrete element model is used to simulate the flow of liquid-granule mixtures in an inclined channel containing a linear contraction. All the relevant particle/particle and particle/fluid interactions are included in the numerical model. The presence of the contraction induces different steady morphologies of the solid phase or the mixture depending on whether closed or open channels are used. These flows behave quite differently depending on the upstream Froude number and the contraction size ratio. The model is first validated by comparing with the existing results for dry granular (glass particles) chute flows (Vreman et al., 2007). Then simulations of a chute of glass particles in water flowing in a closed channel are compared to the dry granular case. With the same solid flux at the inlet, the hydrodynamic forces in the liquid-solid mixture induce higher particle solid volume fractions in the part of the flow containing the solid phase. The streamwise particle velocity (resp. depth of the solid phase) has the same evolution along the channel with smaller (larger) values than in the dry granular flow case.