Simulations of granular media

This page shows the animated part of my thesis. The content is the following:

Other stuff:

I present below some animations obtained in the same manner as a cartoon. We used the technique of simulation called distinct (or discrete) element method (DEM). This means that every grain is modeled separately. For the first five types of experience, we used the rigid-particle method to manage collisions.


Convection in granular media

The following animations (Mpeg format) are simulations of an experience lasting one minute in reality. Computational time is on an average of 22 hours to simulate 60 seconds. You will see various circular motions called convections produced by vibrating the box containing the grains.

Click here (1.5 Mb) or here (1.2 Mb) to see an animation of the convection induced by high amplitude and low frequency vertical vibrations. If we shake the box with low amplitude and high frequency vertical vibrations, the result is quite different (1.2 Mb).

Click here (1.2 Mb) to see an animation of the convection induced by horizontal vibrations.

Click here or here (1.2 Mb) to see an animation of the convection induced by vertical AND horizontal vibrations.

Click here (1.2 Mb) to see an animation of the convection induced by vertical vibrations. The vertical walls go up and down and induce a shear.

Click here (1.2 Mb) to see an animation of the convection induced by vertical vibrations. The vertical walls go up and down in the opposite direction of the previous simulation.

This simulation is the same as the first one, but with friction equal to 0. You will be convinced that friction has a great influence when looking at here (1.7 Mb). However, circular motions appear when vibrations are weak, as you can see here (1.1 Mb).

Some typical streams


Segregation induced by vibrations

When we shake a granular medium, large grains go up. You have certainly noticed this phenomenon by looking at the previous animations concerning convection. Now we increase the difficulty by modeling grains with polygons. Click here (0.5 Mb) to see a simulation of this phenomenon.


Flows

We simulated flows of media composed of particles having varying properties in order to see their influence on the flow time.

Click here (0.3 Mb) to see an animation.

In the experiment below, we compared four media composed. Curves depict the percentage of grains passed in terms of time. P means sharp-pointed aspect, R round aspect, e elastic (coefficient of restitution=0.9, coefficient of friction=0.1), r rugose (0.5, 0.5). Incline of walls is 30 degrees with respect to vertical.

Clearly "elastic" grains flow more quickly than rugose ones. Shape has no effect. On the other hand, the speed of flow does not depend on the number of grains still in the hopper : the flow is linear. "Flat" regions indicate temporary bottlenecks.

Click here (0.3 Mb) to see what happens when we turn the gravity vector on 45 degrees.


Regular structures

What happens in a regular structure when we remove some grains ? Here is our most beautiful animation (0.7 Mb).


Landslide

We simulated a landslide (1.1 Mb).


Penetration

We simulated the penetration of a body inside a granular medium. Contrary to the previous simulations, we used here the soft-particle model, more precisely Cundall's model. Click here (3.1 Mb) to see an experience lasting 0.7 seconds in reality. The large disc is falling down from a height of 23 m.

This method has several advantages. For example we can visualize forces between the grains, as you can see on this image (44 Ko).

The manner grains are arranged plays a pominent part in dissipation of energy, as we can see it on three examples : when the medium is amorphous (183 Ko), when all grains have the same size and are arranged in alternate rows (93 Ko) or arranged on a square mesh 188 (Ko).



Last updated on May 23, 1996. All pages copyright © D. Muller 1995. All Rights Reserved.