The Virgo Consortium

The AREPO Code

Kelvin Helmholtz instability

The Kelvin Helmholtz instability calculated with the moving-mesh method. The figure shows the density field with the Voronoi mesh overlaid in black.


The AREPO code is based on a moving unstructured mesh defined by the Voronoi tessellation of a set of discrete points . The mesh is used to solve the hyperbolic conservation laws of ideal hydrodynamics with a finite volume approach, based on a second-order unsplit Godunov scheme with an exact Riemann solver.

This method is fully Galilean-invariant, unlike ordinary Eulerian codes, a property that is of significant importance for cosmological simulations where highly supersonic bulk flows are common. In addition, it can adjust its spatial resolution automatically and continuously, and hence inherits the principal advantage of SPH for simulations of cosmological structure growth.

Deep wedge of LCDM

Time evolution of the mixing of two fluids induced by the motion of a solid object. This test illustrates the ability of AREPO to cope with arbitrarily curved, moving boundary conditions.


The initial paper for “The AREPO code” has been published as Springel 2010. For more information and relevant publications please visit HITS’s project description page or contact Volker Springel.