Eventos Anais de eventos
ENCIT 2022
19th Brazilian Congress of Thermal Sciences and Engineering
Evaluation of a Parallel High-Order CFD Solver Developed Using Chapel
Submission Author:
Fábio Malacco Moreira , SP
Co-Authors:
Fábio Malacco Moreira, João Luiz F. Azevedo
Presenter: Fábio Malacco Moreira
doi://10.26678/ABCM.ENCIT2022.CIT22-0540
Abstract
The Flux Reconstruction (FR) method is a recent approach to compact high-order spatial discretizations that unifies schemes such as the Discontinuous Galerkin (DG), Spectral Volume (SFV), and Spectral Difference (SD) methods into a single mathematical framework. It also allows for a wide range of new scheme variants to be developed inside this framework and it is generally of simpler computational implementation than its predecessors. The Chapel programming language is also a recent development that integrates parallel programming directives and controls natively into the language. It uses a Partitioned Global Address Space (PGAS) approach in which the data communication in a distributed memory cluster is hidden from the programmer, while still allowing fine control over data locality. The objective of this effort is to build a competitively performant, shared-memory parallel, 2-D Euler flow solver in Chapel using the FR method. While 2nd-order Finite Volume (FV) solvers are common and sufficient for many applications, especially steady-state flows, high-order methods are particularly well suited for problems for which spatial resolution is critical, such as Large Eddy Simulations (LES) and vortex dominated flows. The Chapel programming language offers an opportunity to simplify software development, reducing entry barriers for new researchers, increasing productivity, and, therefore, reducing time spent on computational optimization. Such a programming environment should also allow for easier scaling of the resulting code to a distributed memory solver in the future. The solver developed will be validated through comparisons of flow solutions to some standard CFD test cases, such as the Ringleb flow and airfoil calculations. These comparisons will take into consideration both solution accuracy and computational cost. The work demonstrates both the capability of the FR method in accurately resolving critical 1-D and 2-D benchmark problems and of the Chapel programming language as a viable alternative to the traditional MPI approach. The present effort can be seen as a proof of concept to allow for further development of more complex, 3-D solvers in the future.
Keywords
CFD, Chapel, High-order methods, Flux Reconstruction, Parallel Programming
