Saved in:
Bibliographic Details
Main Authors: Böhm, Fabian, Kohl, Nils, Köstler, Harald, Rüde, Ulrich
Format: Preprint
Published: 2025
Subjects:
Online Access:https://arxiv.org/abs/2511.13109
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866911269783601152
author Böhm, Fabian
Kohl, Nils
Köstler, Harald
Rüde, Ulrich
author_facet Böhm, Fabian
Kohl, Nils
Köstler, Harald
Rüde, Ulrich
contents We propose a robust, adaptive coarse-grid correction scheme for matrix-free geometric multigrid targeting PDEs with strongly varying coefficients. The method combines uniform geometric coarsening of the underlying grid with heterogeneous coarse-grid operators: Galerkin coarse grid approximation is applied locally in regions with large coefficient gradients, while lightweight, direct coarse grid approximation is used elsewhere. This selective application ensures that local Galerkin operators are computed and stored only where necessary, minimizing memory requirements while maintaining robust convergence. We demonstrate the method on a suite of sinker benchmark problems for the generalized Stokes equation, including grid-aligned and unaligned viscosity jumps, smoothly varying viscosity functions with large gradients, and different viscosity evaluation techniques. We analytically quantify the solver's memory consumption and demonstrate its efficiency by solving Stokes problems with $10^{10}$ degrees of freedom, viscosity jumps of $10^{6}$ magnitude, and more than 100{,}000 parallel processes.
format Preprint
id arxiv_https___arxiv_org_abs_2511_13109
institution arXiv
publishDate 2025
record_format arxiv
spellingShingle Large-scale Multigrid with Adaptive Galerkin Coarsening
Böhm, Fabian
Kohl, Nils
Köstler, Harald
Rüde, Ulrich
Performance
76M10, 65N55, 65M60, 35R05, 86A17, 35Q30, 76D07
We propose a robust, adaptive coarse-grid correction scheme for matrix-free geometric multigrid targeting PDEs with strongly varying coefficients. The method combines uniform geometric coarsening of the underlying grid with heterogeneous coarse-grid operators: Galerkin coarse grid approximation is applied locally in regions with large coefficient gradients, while lightweight, direct coarse grid approximation is used elsewhere. This selective application ensures that local Galerkin operators are computed and stored only where necessary, minimizing memory requirements while maintaining robust convergence. We demonstrate the method on a suite of sinker benchmark problems for the generalized Stokes equation, including grid-aligned and unaligned viscosity jumps, smoothly varying viscosity functions with large gradients, and different viscosity evaluation techniques. We analytically quantify the solver's memory consumption and demonstrate its efficiency by solving Stokes problems with $10^{10}$ degrees of freedom, viscosity jumps of $10^{6}$ magnitude, and more than 100{,}000 parallel processes.
title Large-scale Multigrid with Adaptive Galerkin Coarsening
topic Performance
76M10, 65N55, 65M60, 35R05, 86A17, 35Q30, 76D07
url https://arxiv.org/abs/2511.13109