University of Cambridge > > Isaac Newton Institute Seminar Series > 3D Numerical Modelling of Compressible Coupled Magma/Mantle Dynamics With Adaptive Mesh Refinement

3D Numerical Modelling of Compressible Coupled Magma/Mantle Dynamics With Adaptive Mesh Refinement

Add to your list(s) Download to your calendar using vCal

If you have a question about this talk, please contact INI IT.

MIMW03 - From the continuum to the tectonic: the magma/mantle dynamics of planet earth

Juliane Dannberg (, Ryan Grove, Timo Heister (

Melt generation and migration are important processes for the evolution of the Earth's interior and impact the global convection of the mantle.
While they have been the subject of numerous investigations, the typical time and length-scales of melt transport are vastly different from global mantle convection, which determines where melt is generated. This makes it difficult to study mantle convection and melt migration in a unified framework. In addition, modelling magma dynamics poses the challenge of highly non-linear and spatially variable material properties, in particular the viscosity.

Here, we present our extension of the community mantle convection code ASPECT , which adds the equations of two-phase flow of melt and solid, as an example for how these challenges can be addressed. First, We will analyse well-posedness, existence, and uniqueness of the problem. Then we will discuss the correct way to do a stable higher order finite element discretization. Finally, the resulting linear system is solved with an iterated solver preconditioned by a Schur complement-based block preconditioner. We demonstrate that applying adaptive mesh refinement to this type of problem is particularly advantageous, as the resolution can be increased in mesh cells where melt is present and viscosity gradients are high, whereas a lower resolution is sufficient in regions without melt. Together with a high-performance, massively parallel implementation, this allows for high resolution, 3d, compressible, global mantle convection simulations coupled with melt migration.

We present benchmarks of our solver to confirm the theoretical results, and apply our software to large-scale 3d simulations of melting and melt transport in mantle plumes interacting with the lithosphere to show robustness and parallel scalability of the linear solver. Our model incorporates the individual compressibilities of the solid and the fluid phase in addition to compaction, and we demonstrate that including these effects can change melt volumes by more than 20%. Moreover, we show how including melting, melt migration and freezing of melt in global convection models can influence convection patterns and the distribution of chemical heterogeneities in the mantle.

Our model of magma dynamics provides a framework for modelling processes on different scales and investigating links between processes occurring in the deep mantle and melt generation and migration.

This talk is part of the Isaac Newton Institute Seminar Series series.

Tell a friend about this talk:

This talk is included in these lists:

Note that ex-directory lists are not shown.


© 2006-2022, University of Cambridge. Contact Us | Help and Documentation | Privacy and Publicity