A robust matrix-free approach for large-scale non-isothermal high-contrast viscosity Stokes flow on blended domains with applications to geophysics

We consider a compressible Stokes problem in the quasi-stationary case coupled with a time dependent advection-diffusion equation with special emphasis on high viscosity contrast geophysical mantle convection applications. In space, we use a P2-P1 Taylor--Hood element which is generated by a blending approach to account for the non-planar domain boundary without compromising the stencil data structure of uniformly refined elements. In time, we apply an operator splitting approach for the temperature equation combining the BDF2 method for diffusion and a particle method for advection, resulting in an overall second order scheme. Within each time step, a stationary Stokes problem with a high viscosity contrast has to be solved for which we propose a matrix-free, robust and scalable iterative solver based on Uzawa type block preconditioners, polynomial Chebyshev smoothers and a BFBT type Schur complement approximation. Our implementation is using a hybrid hierarchical grid approach allowing for massively parallel, high resolution Earth convection simulations.