@Article{CiCP-18-1,
author = {},
title = {Modeling 3D Magma Dynamics Using a Discontinuous Galerkin Method},
journal = {Communications in Computational Physics},
year = {2015},
volume = {18},
number = {1},
pages = {230--246},
abstract = {<p style="text-align: justify;">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.</p>},
issn = {1991-7120},
doi = {https://doi.org/10.4208/cicp.090314.151214a},
url = {https://global-sci.com/article/80318/modeling-3d-magma-dynamics-using-a-discontinuous-galerkin-method}
}