Hostname: page-component-78c5997874-ndw9j Total loading time: 0 Render date: 2024-11-10T14:49:50.247Z Has data issue: false hasContentIssue false
  • English
  • Français

Modeling 3D Magma Dynamics Using a Discontinuous Galerkin Method

Part of: Partial differential equations, initial value and time-dependent initial-boundary value problems

Published online by Cambridge University Press:  03 July 2015

Seshu Tirupathi ,
Jan S. Hesthaven  and
Yan Liang
Seshu Tirupathi*
Affiliation:
Division of Applied Mathematics, Brown University, Providence, RI 02912, USA IBM Research, Dublin, Damastown Industrial Park, Dublin 15, Ireland
Jan S. Hesthaven
Affiliation:
Division of Applied Mathematics, Brown University, Providence, RI 02912, USA EPFL-SB-MATHICSE-MCSS, Ecole Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
Yan Liang
Affiliation:
Department of Geological Sciences, Brown University, Providence, RI 02912, USA
*
*Corresponding author. Email addresses: seshutir@ie.ibm.com (S. Tirupathi), jan.hesthaven@epfl.ch (J. S. Hesthaven), Yan_Liang@brown.edu (Y. Liang)
Get access

Abstract

Discontinuous Galerkin (DG) and matrix-free finite element methods with a novel projective pressure estimation are combined to enable the numerical modeling of magma dynamics in 2D and 3D using the library deal.II. The physical model is an advection-reaction type system consisting of two hyperbolic equations to evolve porosity and soluble mineral abundance at local chemical equilibrium and one elliptic equation to recover global pressure. A combination of a discontinuous Galerkin method for the advection equations and a finite element method for the elliptic equation provide a robust and efficient solution to the channel regime problems of the physical system in 3D. A projective and adaptively applied pressure estimation is employed to significantly reduce the computational wall time without impacting the overall physical reliability in the modeling of important features of melt segregation, such as melt channel bifurcation in 2D and 3D time dependent simulations.

Keywords

magma dynamicsdiscontinuous Galerkin methodfinite element methodmatrix-free method

MSC classification

Secondary: 65Z05: Applications to physics 65M22: Solution of discretized equations 65M60: Finite elements, Rayleigh-Ritz and Galerkin methods, finite methods
Type
Research Article
Information
Communications in Computational Physics , Volume 18 , Issue 1 , July 2015 , pp. 230 - 246
Copyright
Copyright © Global-Science Press 2015 

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]Abreu, E., Douglas, J., Furtado, F. and Pereira, F., Operator splitting for three-phase flow in heterogeneous porous media, Commun. Comput. Phys., 6 (2009), 72–84.CrossRefGoogle Scholar
[2]Aharonov, E., Whitehead, J., Kelemen, P. B. and Spiegelman, M. S, Channeling instability of upwelling melt in the mantle, J. Geophys. Res., 100 (1995), 20433–20450.Google Scholar
[3]Bangerth, W., Hartmann, R. and Kanschat, G., deal.II: A general-purpose object-oriented finite element library, ACM Transactions on Mathematical Software, 33 (2007), 1–27.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]Chueh, C., Djilali, N. and Bangerth, W., An h-adaptive operator splitting method for two-phase flow in 3D heterogeneous porous media, SIAM J. Sci. Comput., 35 (2013), B149–B175.Google Scholar
[7]Hesse, M.A., Schiemenz, A. R., Liang, Y. and Parmentier, E.M., Compaction-dissolution waves in an upwelling mantle column, Geophys. J. Int., 187 (2011), 1057–1075.CrossRefGoogle Scholar
[8]Hesthaven, J. S. and Warburton, T., Nodal discontinuous Galerkin methods: Algorithms, analysis, and applications, Springer, 2008.Google Scholar
[9]Kelemen, P. B., The origin of the land under the sea, Scientific American, 300 (2009), 52–57.CrossRefGoogle ScholarPubMed
[10]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
[11]Kronbichler, M. and Kormann, K., A generic interface for parallel cell-based finite element operator application, Fluids Comput., 63 (2012), 135–147.Google Scholar
[12]Liang, Y., Schiemenz, A. R., Hesse, M., Parmentier, E. M. and Hesthaven, J.S., High-porosity channels for melt migration in the mantle: Top is the dunite and bottom is the harzburgite and lherzolite, Geophys. Res. Lett., 37 (2010), L15306.Google Scholar
[13]Liang, Y., Schiemenz, A. R. and Hesse, M. A.Waves, channels, and the preservation of chemical heterogeneities during melt migration in the mantle, Geophys. Res. Lett., 38 (2011), L20308.Google Scholar
[14]Saad, Y., Iterative methods for sparse linear systems (2nd edition), SIAM Publishing, 2003.Google Scholar
[15]Schiemenz, A. R., Hesse, M. A. and Hesthaven, J. S.Modeling magma dynamics with a mixed Fourier collocation – discontinuous Galerkin method, Commun. Comput. Phys., 10 (2011), 433–452.CrossRefGoogle Scholar
[16]Schiemenz, A. R., Liang, Y. and Parmentier, E.M., A high-order numerical study of reactive dissolution in an upwelling heterogeneous mantle: I. Channelization, channel lithology, and channel geometry, Geophys. J. Int., 186 (2011), 641–664.Google Scholar
[17]Spiegelman, M. and Kelemen, P. B.Extreme chemical variability as a consequence of channelized melt transport, Geochem. Geophys. Geosyst., 4 (2003), 1055.Google Scholar
[18]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.Google Scholar