A Multigrid reduction framework for domains with symmetries
Divergence constraints are present in the governing equations of numerous physical phenomena, and they usually lead to a Poisson equation whose solution represents a bottleneck in many simulation codes. Algebraic multigrid (AMG) is arguably the most powerful preconditioner for Poisson’s equation, an...
| Autores: | , , , |
|---|---|
| Formato: | artículo |
| Fecha de publicación: | 2024 |
| País: | España |
| Recursos: | Universitat Politècnica de Catalunya (UPC) |
| Repositorio: | UPCommons. Portal del coneixement obert de la UPC |
| Idioma: | inglés |
| OAI Identifier: | oai:upcommons.upc.edu:2117/426497 |
| Acesso em linha: | https://hdl.handle.net/2117/426497 https://dx.doi.org/10.1137/24M1638513 |
| Access Level: | acceso abierto |
| Palavra-chave: | AMG Multigrid reduction SpMM Spatial symmetries Poisson’s equation Àrees temàtiques de la UPC::Enginyeria mecànica::Mecànica de fluids Àrees temàtiques de la UPC::Física::Termodinàmica |
| Resumo: | Divergence constraints are present in the governing equations of numerous physical phenomena, and they usually lead to a Poisson equation whose solution represents a bottleneck in many simulation codes. Algebraic multigrid (AMG) is arguably the most powerful preconditioner for Poisson’s equation, and its effectiveness results from the complementary roles played by the smoother, responsible for damping high-frequency error components, and the coarse-grid correction, which in turn reduces low-frequency modes. This work presents several strategies to make AMG more compute-intensive by leveraging reflection, translational, and rotational symmetries. AMGR, our final proposal, does not require boundary conditions to be symmetric, therefore applying to a broad range of academic and industrial configurations. It is based on a multigrid reduction framework that introduces an aggressive coarsening to the multigrid hierarchy, reducing the memory footprint, setup, and application costs of the top-level smoother. While preserving AMG’s excellent convergence, AMGR allows one to replace the standard sparse matrix-vector product with the more compute-intensive sparse matrix-matrix product, yielding significant accelerations. Numerical experiments on industrial CFD applications demonstrated up to 70% speed-ups when solving Poisson’s equation with AMGR instead of AMG. Additionally, strong and weak scalability analyses revealed no significant degradation. |
|---|