3D Numerical Modelling of Compressible Coupled Magma/Mantle Dynamics With Adaptive Mesh Refinement
Duration: 2 hours 22 mins
Share this media item:
Embed this media item:
Embed this media item:
About this item
| Description: |
Heister, T - Dannberg, J (Texas A&M University)
Friday 10th June 2016 - 09:00 to 11:00 |
|---|
| Created: | 2016-06-16 10:24 |
|---|---|
| Collection: | Melt in the Mantle |
| Publisher: | Isaac Newton Institute |
| Copyright: | Heister, T - Dannberg, J |
| Language: | eng (English) |
| Distribution: |
World
(downloadable)
|
| Explicit content: | No |
| Aspect Ratio: | 16:9 |
| Screencast: | No |
| Bumper: | UCS Default |
| Trailer: | UCS Default |
| Abstract: | Juliane Dannberg (dannberg@math.tamu.edu), Ryan Grove, Timo Heister (heister@clemson.edu)
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. |
|---|---|
Available Formats
| Format | Quality | Bitrate | Size | |||
|---|---|---|---|---|---|---|
| MPEG-4 Video | 640x360 | 1.94 Mbits/sec | 2.03 GB | View | Download | |
| WebM | 640x360 | 576.07 kbits/sec | 599.14 MB | View | Download | |
| iPod Video | 480x270 | 524.1 kbits/sec | 545.09 MB | View | Download | |
| MP3 | 44100 Hz | 250.93 kbits/sec | 260.98 MB | Listen | Download | |
| Auto * | (Allows browser to choose a format it supports) | |||||