The Lagrangian hydrodynamics code MAGMA2

Understand

We present the methodology and performance of the new Lagrangian hydrodynamics code MAGMA2, a Smoothed Particle Hydrodynamics code that benefits from a number of non-standard enhancements.

• By default it uses high-order smoothing kernels and wherever gradients are needed, they are calculated via accurate matrix inversion techniques, but a more conventional formulation with kernel gradients has also been implemented for comparison purposes.
• We also explore a matrix inversion formulation of SPH with a symmetrisation in the particle indices that is not frequently used.
• We find interesting advantages of this formulation in some of the tests, for example, a substantial reduction of surface tension effects for non-ideal particle setups and more accurate peak densities in Sedov blast waves.