A numerical code for the simulation of magma-rocks dynamics

Abstract

We present a numerical code for the simulation of the dynamics of compressible to incompressible, multicomponent ows, based on the _nite element algorithm by Hauke & Hughes (1998). Balance equations for mass, momentum, energy and composition are solved with space-time Galerkin least-squares and discontinuity-capturing stabilizing techniques. The code is used to study the dynamics of convection and mixing in magmatic systems such as replenishment of magma chambers and volcanic conduits, and it reveals the occurrence of previously not described processes. The uid-structure interaction of fully coupled magma-rock dynamics is being implemented by using the deforming-spatial domain method by Tezduyar (2006), that intrinsecally includes moving meshes.