TY - JOUR T1 - Modeling Magma Dynamics with a Mixed Fourier Collocation - Discontinuous Galerkin Method AU - Alan R. Schiemenz, Marc A. Hesse & Jan S. Hesthaven JO - Communications in Computational Physics VL - 2 SP - 433 EP - 452 PY - 2011 DA - 2011/10 SN - 10 DO - http://doi.org/10.4208/cicp.030210.240910a UR - https://global-sci.org/intro/article_detail/cicp/7449.html KW - AB -
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.