New numerical approaches for modeling thermochemical convection in a compositionally stratified fluid

TL;DR

This paper introduces two novel numerical algorithms, a Discontinuous Galerkin method with a Bound Preserving limiter and a Volume-of-Fluid approach, to improve the sharpness of compositional interfaces in mantle convection simulations.

Contribution

It presents two new algorithms for more accurately modeling sharp compositional boundaries in thermochemical convection, addressing limitations of existing methods.

Findings

New algorithms maintain sharper compositional boundaries.

Comparison shows improved interface sharpness over traditional methods.

Algorithms integrated into open-source ASPECT code.

Abstract

Seismic imaging of the mantle has revealed large and small scale heterogeneities in the lower mantle; specifically structures known as large low shear velocity provinces (LLSVP) below Africa and the South Pacific. Most interpretations propose that the heterogeneities are compositional in nature, differing in composition from the overlying mantle, an interpretation that would be consistent with chemical geodynamic models. Numerical modeling of persistent compositional interfaces presents challenges, even to state-of-the-art numerical methodology. For example, some numerical algorithms for advecting the compositional interface cannot maintain a sharp compositional boundary as the fluid migrates and distorts with time dependent fingering due to the numerical diffusion that has been added in order to maintain the upper and lower bounds on the composition variable and the stability of the advection method. In this work we present two new algorithms for maintaining a sharper computational boundary than the advection methods that are currently openly available to the computational mantle convection community; namely, a Discontinuous Galerkin method with a Bound Preserving limiter and a Volume-of-Fluid interface tracking algorithm. We compare these two new methods with two approaches commonly used for modeling the advection of two distinct, thermally driven, compositional fields in mantle convection problems; namely, an approach based on a high-order accurate finite element method advection algorithm that employs an artificial viscosity technique to maintain the upper and lower bounds on the composition variable as well as the stability of the advection algorithm and the advection of particles that carry a scalar quantity representing the location of each compositional field. All four of these algorithms are implemented in the open source FEM code ASPECT.