2025/12/09 by Poncelet, Gilles, Lambrechts, Jonathan, Gillis, Thomas +1
Engineering · #35-04 #35J99 #65N06 #65N50 #65N55 #65N80 #68W15 #Advanced Numerical Methods in Computational Mathematics #Computational Fluid Dynamics and Aerodynamics #Distributed #FOS: Computer and information sciences #FOS: Mathematics #Lattice Boltzmann Simulation Studies #Numerical Analysis (math.NA) #Parallel #and Cluster Computing (cs.DC)
paper · doi:10.48550/arxiv.2512.08555
openalex publication_date 2025/12/09 · openalex created_date 2025/12/11 · openalex updated_date 2026/07/28
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.