A blender in a granular medium

Опубликовано: 15 Март 2025
на канале: Nils Berglund
1,017
45

This simulation of a granular medium features a "blender", composed of three blades angled at 120 degrees. The blender starts at rest, to give the medium a chance to reach an equilibrium configuration. Then its angular speed increases, before decreasing again. The angle is given in degrees, and the angular speed in radians per simulation second.
The particles in this simulation interact via a Lennard-Jones potential, and are subject to a viscous damping term, gravity directed downwards, and to a one-sided harmonic repulsion from the blender. There are periodic boundary conditions.
This simulation has two parts, showing the evolution with two different color schemes:
Initial y coordinate: 0:00
Direction of velocity: 2:08
In the first part, each particle is given a color at the beginning of the simulation, which is a function of its y-coordinate. The particle then keeps the same color throughout the simulation.In the second part, the hue of the particles depends on the direction of their velocity, and the luminosity depends on their speed.
To save on computation time, particles are placed into a "hash grid", each cell of which contains between 3 and 10 particles. Then only the influence of other particles in the same or neighboring cells is taken into account for each particle.
The Lennard-Jones potential is strongly repulsive at short distance, and mildly attracting at long distance. It is widely used as a simple yet realistic model for the motion of electrically neutral molecules. The force results from the repulsion between electrons due to Pauli's exclusion principle, while the attractive part is a more subtle effect appearing in a multipole expansion. For more details, see https://en.wikipedia.org/wiki/Lennard...

Render time: 1 hour 13 minutes
Color scheme: Part 1 - HSV/Jet
Part 2 - Twilight by Bastian Bechtold
https://github.com/bastibe/twilight

Music: "Simmering" by Asher Fulero‪@AsherFulero‬

Current version of the C code used to make these animations:
https://github.com/nilsberglund-orlea...
https://www.idpoisson.fr/berglund/sof...
Some outreach articles on mathematics:
https://images.math.cnrs.fr/_Berglund...
(in French, some with a Spanish translation)

#granular_material #lennardjones #LennardJones #molecular_dynamics #phonon