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

Modeling Magma Dynamics with a Mixed Fourier Collocation - Discontinuous Galerkin Method

Alan R. Schiemenz 1*, Marc A. Hesse 2, Jan S. Hesthaven 3

1 Department of Geological Sciences and Division of Applied Mathematics, Brown University, Providence, RI 02910, USA; and Department of Earth and Environmental Sciences, Geophysics, Munich University, Theresienstr. 41, 80333 Munich, Germany.
2 Department of Geological Sciences, University of Texas at Austin, Austin, TX 78712, USA.
3 Division of Applied Mathematics, Brown University, Providence, RI 02910, USA.

Received 3 February 2010; Accepted (in revised version) 24 September 2010
Available online 28 April 2011


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.

AMS subject classifications: 65Z05, 65M22, 65M70, 65T40, 65T50

Notice: Undefined variable: pac in /var/www/html/issue/abstract/readabs.php on line 164
Key words: Magma dynamics, discontinuous Galerkin method, Fourier collocation method, numerical simulations, linear stability analysis.

*Corresponding author.
Email: (A. R. Schiemenz), (M. A. Hesse), (J. S. Hesthaven)

The Global Science Journal