Hostname: page-component-8448b6f56d-cfpbc Total loading time: 0 Render date: 2024-04-20T00:29:41.689Z Has data issue: false hasContentIssue false
  • English
  • Français

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

Published online by Cambridge University Press:  20 August 2015

Alan R. Schiemenz ,
Marc A. Hesse  and
Jan S. Hesthaven
Alan R. Schiemenz*
Affiliation:
Department of Geological Sciences, Brown University, Providence, RI 02910, USA Division of Applied Mathematics, Brown University, Providence, RI 02910, USA
Marc A. Hesse*
Affiliation:
Department of Geological Sciences, University of Texas at Austin, Austin, TX 78712, USA
Jan S. Hesthaven*
Affiliation:
Division of Applied Mathematics, Brown University, Providence, RI 02910, USA
*
Corresponding author.Email:alan.schiemenz@geophysik.uni-muenchen.de
Email:mhesse@jsg.utexas.edu
Email:Jan.Hesthaven@brown.edu
Get access

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.

Keywords

65Z0565M2265M7065T4065T50Magma dynamicsdiscontinuous Galerkin methodFourier collocation methodnumerical simulationslinear stability analysis
Type
Research Article
Information
Communications in Computational Physics , Volume 10 , Issue 2 , August 2011 , pp. 433 - 452
Copyright
Copyright © Global Science Press Limited 2011

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

[1] Aharonov, E., Whitehead, J., Kelemen, P. B. and Spiegelman, M., Channeling instability of upwelling melt in the mantle, J. Geophys. Res., 100 (1995), 433–450.Google Scholar
[2] Arnold, D. N., Brezzi, F., Cockburn, B. and Marini, L. D., Unified analysis of discontinuous galerkin methods for elliptic problems, SIAM J. Numer. Anal., 39 (2002), 1749–1779.Google Scholar
[3] Bercovici, D., Ricard, Y. and Schubert, G., A two-phase model for compaction and damage 1. General theory, J. Geophys. Res., 106 (2001), 8887–8906.Google Scholar
[4] Buneman, O., Diagnosing oscillatory growth or decay, J. Comput. Phys., 29 (1978), 295–296.Google Scholar
[5] Carpenter, M. H. and Kennedy, C. A., Fourth-order 2N-storage Runge-Kutta schemes, NASA TM 109112, NASA Langley Research Center, 1994.Google Scholar
[6] Hesse, M.A., Schiemenz, A. R., Liang, Y. and Parmentier, E.M., Compaction-dissolution waves in an upwelling mantle column, Geophys. J. Int., 2011 (submitted).CrossRefGoogle Scholar
[7] Hesthaven, J. S., Gottlieb, S. and Gottlieb, D., Spectral Methods for Time-Dependent Problems, Cambridge University Press, 2006.Google Scholar
[8] Hesthaven, J. S. and Warburton, T., Nodal Discontinuous Galerkin Methods: Algorithms, Analysis, and Applications, Springer, 2008.Google Scholar
[9] Kelemen, P. B., Hirth, G., Shimizu, N., Spiegelman, M. and Dick, H. J. B., A review of melt migration processes in the adiabatically upwelling mantle beneath oceanic spreading ridges, Phil. Trans. R. Soc. Lond., 355 (1997), 282–318.Google Scholar
[10] Saad, Y., Iterative Methods for Sparse Linear Systems (2nd edition), SIAM Publishing, 2003.CrossRefGoogle Scholar
[11] Schiemenz, A. R., Liang, Y. and Parmentier, E. M., A high-order numerical study of reactive dissolution in an upwelling heterogenous mantle: I. Channelization, channel lithology, and channel geometry, Geophys. J. Int., 2011. DOI: 10.1111/j.1365-246X.2011.05065.x (in press).Google Scholar
[12] Spiegelman, M. and Kelemen, P. B., Extreme chemical variability as a consequence of channelized melt transport, Geochem. Geophys. Geosyst., 4(7) (2003), 1055.Google Scholar
[13] Spiegelman, M., Kelemen, P. B. and Aharonov, E., Causes and consequences of flow organization during melt transport: The reaction infiltration instability in compactible media, J. Geo. Res., 106 (2001), 2061–2077.CrossRefGoogle Scholar