Efficient parallel 3D computation of the compressible Euler equations with an invariant-domain preserving second-order finite-element scheme

arXiv2007.00094MaRDI QIDQ6344148FDOQ6344148

Authors: Matthias Maier, Martin Kronbichler 

Publication date: 30 June 2020

Abstract: We discuss the efficient implementation of a high-performance second-order collocation-type finite-element scheme for solving the compressible Euler equations of gas dynamics on unstructured meshes. The solver is based on the convex limiting technique introduced by Guermond et al. (SIAM J. Sci. Comput. 40, A3211-A3239, 2018). As such it is invariant-domain preserving, i.e., the solver maintains important physical invariants and is guaranteed to be stable without the use of ad-hoc tuning parameters. This stability comes at the expense of a significantly more involved algorithmic structure that renders conventional high-performance discretizations challenging. We develop an algorithmic design that allows SIMD vectorization of the compute kernel, identify the main ingredients for a good node-level performance, and report excellent weak and strong scaling of a hybrid thread/MPI parallelization.

Has companion code repository: https://github.com/conservation-laws/ryujin

This page was built for publication: Efficient parallel 3D computation of the compressible Euler equations with an invariant-domain preserving second-order finite-element scheme

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6344148)