Opin vísindi

A FFT-accelerated multi-block finite-difference solver for massively parallel simulations of incompressible flows

A FFT-accelerated multi-block finite-difference solver for massively parallel simulations of incompressible flows


Title: A FFT-accelerated multi-block finite-difference solver for massively parallel simulations of incompressible flows
Author: Costa, Pedro Simoes
Date: 2022-02
Language: English
Scope:
University/Institute: University of Iceland
Series: Computer Physics Communications; 271()
ISSN: 0010-4655
DOI: https://doi.org/10.1016/j.cpc.2021.108194
Subject: Tölvunarfræði; Vélbúnaður (tölvur); Computational fluid dynamics; Direct numerical simulation; Fast Poisson solver; High-performance computing; Multi-block solver; Hardware and Architecture; Physics and Astronomy (all)
URI: https://hdl.handle.net/20.500.11815/3016

Show full item record

Citation:

Costa , P S 2022 , ' A FFT-accelerated multi-block finite-difference solver for massively parallel simulations of incompressible flows ' , Computer Physics Communications , vol. 271 , 108194 . https://doi.org/10.1016/j.cpc.2021.108194

Abstract:

We present a multi-block finite-difference solver for massively parallel Direct Numerical Simulations (DNS) of incompressible flows. The algorithm combines the versatility of a multi-block solver with the method of eigenfunctions expansions, to speedup the solution of the pressure Poisson equation. This is achieved by employing FFT-based transforms along one homogeneous direction, which effectively reduce the problem complexity at a low cost. These FFT-based expansions are implemented in a framework that unifies all valid combinations of boundary conditions for this type of method. Subsequently, a geometric multigrid solver is employed to solve the reduced Poisson equation in a multi-block geometry. Particular care was taken here, to guarantee the parallel performance of the multigrid solver when solving the reduced linear systems equations. We have validated the overall numerical algorithm and assessed its performance in a massively parallel setting. The results show that 2- to 8-fold reductions in computational cost may be easily achieved when exploiting FFT-accelerated for the solution of the Poisson equation. The solver, SNaC, has been made freely available and open-source under the terms of an MIT license.

Description:

© 2021 Elsevier B.V. This work was supported by the University of Iceland Recruitment Fund grant No. 1515-151341, TURBBLY.

Files in this item

This item appears in the following Collection(s)