Multilevel and Local Timestepping Discontinuous Galerkin Methods for Magma Dynamics

Discontinuous Galerkin (DG) method is presented for numerical modeling of melt migration in a chemically reactive and viscously deforming upwelling mantle column. DG methods for both advection and elliptic equations provide a robust and efficient solution to the problems of melt migration in the asthenospheric upper mantle. Assembling and solving the elliptic equation is the major bottleneck in these computations. To address this issue, adaptive mesh refinement and local timestepping methods have been proposed to significantly improve the computational wall time. The robustness of DG methods is demonstrated through two benchmark problems by modeling detailed structure of high-porosity dissolution channels and compaction-dissolution waves.


Published in:
Computational Geosciences, 19, 4, 965-978
Year:
2015
Publisher:
Dordrecht, Springer
ISSN:
1420-0597
Keywords:
Laboratories:




 Record created 2014-03-07, last modified 2018-03-17

Postprint:
Download fulltext
PDF

Rate this document:

Rate this document:
1
2
3
 
(Not yet reviewed)