Modeling 3D Magma Dynamics Using a Discontinuous Galerkin Method

Modeling 3D Magma Dynamics Using a Discontinuous Galerkin Method

Year:    2015

Communications in Computational Physics, Vol. 18 (2015), Iss. 1 : pp. 230–246

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.

You do not have full access to this article.

Already a Subscriber? Sign in as an individual or via your institution

Journal Article Details

Publisher Name:    Global Science Press

Language:    English

DOI:    https://doi.org/10.4208/cicp.090314.151214a

Communications in Computational Physics, Vol. 18 (2015), Iss. 1 : pp. 230–246

Published online:    2015-01

AMS Subject Headings:    Global Science Press

Copyright:    COPYRIGHT: © Global Science Press

Pages:    17

Keywords:   

  1. A hybridized discontinuous Galerkin method for a linear degenerate elliptic equation arising from two-phase mixtures

    Kang, Shinhoo | Bui-Thanh, Tan | Arbogast, Todd

    Computer Methods in Applied Mechanics and Engineering, Vol. 350 (2019), Iss. P.315

    https://doi.org/10.1016/j.cma.2019.03.018 [Citations: 5]
  2. Melt Segregation and Depletion During Ascent of Buoyant Diapirs in Subduction Zones

    Zhang, Nan | Behn, Mark D. | Parmentier, E. Marc | Kincaid, Christopher

    Journal of Geophysical Research: Solid Earth, Vol. 125 (2020), Iss. 2

    https://doi.org/10.1029/2019JB018203 [Citations: 12]
  3. Multilevel and local time-stepping discontinuous Galerkin methods for magma dynamics

    Tirupathi, S. | Hesthaven, J. S. | Liang, Y. | Parmentier, M.

    Computational Geosciences, Vol. 19 (2015), Iss. 4 P.965

    https://doi.org/10.1007/s10596-015-9514-7 [Citations: 7]
  4. A scalable exponential-DG approach for nonlinear conservation laws: With application to Burger and Euler equations

    Kang, Shinhoo | Bui-Thanh, Tan

    Computer Methods in Applied Mechanics and Engineering, Vol. 385 (2021), Iss. P.114031

    https://doi.org/10.1016/j.cma.2021.114031 [Citations: 1]
  5. A high‐order numerical study of reactive dissolution in an upwelling heterogeneous mantle: 2. Effect of shear deformation

    Baltzell, Conroy | Parmentier, E. M. | Liang, Yan | Tirupathi, Seshu

    Geochemistry, Geophysics, Geosystems, Vol. 16 (2015), Iss. 11 P.3855

    https://doi.org/10.1002/2015GC006038 [Citations: 12]