The Spherical Shell NSinker benchmark

The Spherical Shell NSinker benchmark#

This folder contains the “Spherical Shell NSinker” benchmark, which is loosely based on the NSinker benchmark defined in

Weighted BFBT Preconditioner for Stokes Flow Problems with Highly Heterogeneous Viscosity, Johann Rudi, Georg Stadler, Omar Ghattas, SIAM Journal of Scientific Computing, Vol 39, Issue 5, 2017. https://doi.org/10.1137/16M108450X

but adapted to a spherical shell geometry. See the “nsinker” benchmark in this repository for the version in a unit box. The motivation for creating this benchmark was the desire to test linear solver performance with a spherical geometry, which is often used in practice and that turns out to be more challenging.

This benchmark is defined and used in

Efficient distributed matrix-free multigrid methods on locally refined meshes for FEM computations, Peter Munch, Timo Heister, Laura Prieto Saavedra, Martin Kronbichler, ACM Transactions on Parallel Computing, 2023, 10, 1-38. https://doi.org/10.1145/3580314

We consider a spherical shell with inner radius 0.54 and outer radius 1.0 with N spherical, high-viscosity, high-density sinking spheres at random (but predetermined in the code) locations in the domain. Additionally, a transition zone with distance 0.89583 from the origin (or about 660km depth if this is assumed to be a non-dimensionalized Earth-like computation) with a sharp viscosity jump can be enabled.