A scalable high-order multigrid-FFT Poisson solver for unbounded domains on adaptive multiresolution grids

Fuente: arXiv
Gespeichert in:
Bibliographische Detailangaben
Hauptverfasser: Poncelet, Gilles, Lambrechts, Jonathan, Gillis, Thomas, Chatelain, Philippe
Format: Preprint
Veröffentlicht: 2025
Schlagworte:
Online-Zugang:
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
_version_ 1866908701689905152
author Poncelet, Gilles
Lambrechts, Jonathan
Gillis, Thomas
Chatelain, Philippe
author_facet Poncelet, Gilles
Lambrechts, Jonathan
Gillis, Thomas
Chatelain, Philippe
contents Multigrid solvers are among the most efficient methods for solving the Poisson equation, which is ubiquitous in computational physics. For example, in the context of incompressible flows, it is typically the costliest operation. The present document expounds upon the implementation of a flexible multigrid solver that is capable of handling any type of boundary conditions within murphy, a multiresolution framework for solving partial differential equations (PDEs) on collocated adaptive grids. The utilization of a Fourier-based direct solver facilitates the attainment of flexibility and enhanced performance by accommodating any combination of unbounded and semi-unbounded boundary conditions. The employment of high-order compact stencils contributes to the reduction of communication demands while concurrently enhancing the accuracy of the system. The resulting solver is validated against analytical solutions for periodic and unbounded domains. In conclusion, the solver has been demonstrated to demonstrate scalability to 16,384 cores within the context of leading European high-performance computing infrastructures.
format Preprint
id arxiv_https___arxiv_org_abs_2512_08555
institution arXiv
publishDate 2025
record_format arxiv
spellingShingle A scalable high-order multigrid-FFT Poisson solver for unbounded domains on adaptive multiresolution grids
Poncelet, Gilles
Lambrechts, Jonathan
Gillis, Thomas
Chatelain, Philippe
Numerical Analysis
Distributed, Parallel, and Cluster Computing
35J99, 35-04, 65N06, 65N50, 65N55, 65N80, 68W15
Multigrid solvers are among the most efficient methods for solving the Poisson equation, which is ubiquitous in computational physics. For example, in the context of incompressible flows, it is typically the costliest operation. The present document expounds upon the implementation of a flexible multigrid solver that is capable of handling any type of boundary conditions within murphy, a multiresolution framework for solving partial differential equations (PDEs) on collocated adaptive grids. The utilization of a Fourier-based direct solver facilitates the attainment of flexibility and enhanced performance by accommodating any combination of unbounded and semi-unbounded boundary conditions. The employment of high-order compact stencils contributes to the reduction of communication demands while concurrently enhancing the accuracy of the system. The resulting solver is validated against analytical solutions for periodic and unbounded domains. In conclusion, the solver has been demonstrated to demonstrate scalability to 16,384 cores within the context of leading European high-performance computing infrastructures.
title A scalable high-order multigrid-FFT Poisson solver for unbounded domains on adaptive multiresolution grids
topic Numerical Analysis
Distributed, Parallel, and Cluster Computing
35J99, 35-04, 65N06, 65N50, 65N55, 65N80, 68W15
url https://arxiv.org/abs/2512.08555