Volume 10, Issue 2
Modeling Magma Dynamics with a Mixed Fourier Collocation - Discontinuous Galerkin Method

Alan R. Schiemenz ,  Marc A. Hesse and Jan S. Hesthaven


Commun. Comput. Phys., 10 (2011), pp. 433-452.

Preview Full PDF BiBTex 82 375
  • Abstract

A high-order discretization consisting of a tensor product of the Fourier collocation and discontinuous Galerkin methods is presented for numerical modeling of magma dynamics. The physical model is an advection-reaction type system consisting of two hyperbolic equations and one elliptic equation. The high-order solution basis allows for accurate and efficient representation of compaction-dissolution waves that are predicted from linear theory. The discontinuous Galerkin method provides a robust and efficient solution to the eigenvalue problem formed by linear stability analysis of the physical system. New insights into the processes of melt generation and segregation, such as melt channel bifurcation, are revealed from two-dimensional time-dependent simulations.

  • History

Published online: 2011-10

  • Keywords

  • AMS Subject Headings

  • Cited by