Comparison Between Algebraic and Matrix-free Geometric Multigrid for a Stokes Problem on Adaptive Meshes with Variable Viscosity

07/15/2019
by   Thomas C. Clevenger, et al.
0

Problems arising in Earth's mantle convection involve finding the solution to Stokes systems with large viscosity contrasts. These systems contain localized features which, even with adaptive mesh refinement, result in linear systems that can be on the order of 100+ million unknowns. One common approach for preconditioning to the velocity block of these systems is to apply an Algebraic Multigrid (AMG) v-cycle (as is done in the ASPECT software, for example), however, with AMG, robustness can be difficult with respect to problem size and number of parallel processes. Additionally, we see an increase in iteration counts with adaptive refinement when using AMG. In contrast, the Geometric Multigrid (GMG) method, by using information about the geometry of the problem, should offer a more robust option. Here we present a matrix-free GMG v-cycle which works on adaptively refined, distributed meshes, and we will compare it against the current AMG preconditioner (Trilinos ML) used in the ASPECT software. We will demonstrate the robustness of GMG with respect to problem size and show scaling up to 24576 cores and 2.2B unknowns. All computations are run using the open source, finite element library deal.II.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
09/30/2021

The Scott-Vogelius Method for Stokes Problem on Anisotropic Meshes

This paper analyzes the Scott-Vogelius divergence-free element pair on a...
research
04/05/2019

A Flexible, Parallel, Adaptive Geometric Multigrid method for FEM

We present data structures and implementation details of a geometric mul...
research
03/23/2022

Efficient distributed matrix-free multigrid methods on locally refined meshes for FEM computations

This work studies three multigrid variants for matrix-free finite-elemen...
research
03/19/2017

Coupling parallel adaptive mesh refinement with a nonoverlapping domain decomposition solver

We study the effect of adaptive mesh refinement on a parallel domain dec...
research
10/05/2020

hr-adaptivity for nonconforming high-order meshes with the target matrix optimization paradigm

We present an hr-adaptivity framework for optimization of high-order mes...
research
11/04/2020

Local Coarsening Algorithms on Adaptively Refined Meshes in 2D and Their Efficient Implementation in MATLAB

Adaptive meshing includes local refinement as well as coarsening of mesh...
research
10/24/2019

The deal.II finite element library: design, features, and insights

deal.II is a state-of-the-art finite element library focused on generali...

Please sign up or login with your details

Forgot password? Click here to reset