High performance computing aspects of a dimension independent semi-Lagrangian discontinuous Galerkin code
From MaRDI portal
(Redirected from Publication:339403)
parallelizationhigh-performance computingVlasov-Poisson equationssemi-Lagrangian discontinuous Galerkin methods
Parallel numerical computation (65Y05) Vlasov equations (35Q83) Finite element, Rayleigh-Ritz and Galerkin methods for initial value and initial-boundary value problems involving PDEs (65M60) Software, source code, etc. for problems pertaining to numerical analysis (65-04) Software, source code, etc. for problems pertaining to partial differential equations (35-04)
Abstract: The recently developed semi-Lagrangian discontinuous Galerkin approach is used to discretize hyperbolic partial differential equations (usually first order equations). Since these methods are conservative, local in space, and able to limit numerical diffusion, they are considered a promising alternative to more traditional semi-Lagrangian schemes (which are usually based on polynomial or spline interpolation). In this paper, we consider a parallel implementation of a semi-Lagrangian discontinuous Galerkin method for distributed memory systems (so-called clusters). Both strong and weak scaling studies are performed on the Vienna Scientific Cluster 2 (VSC-2). In the case of weak scaling, up to 8192 cores, we observe a parallel efficiency above 0.89 for both two and four dimensional problems. Strong scaling results show good scalability to at least 1024 cores (we consider problems that can be run on a single processor in reasonable time). In addition, we study the scaling of a two dimensional Vlasov--Poisson solver that is implemented using the framework provided. All of the simulation are conducted in the context of worst case communication overhead; i.e., in a setting where the CFL number increases linearly with the problem size. The framework introduced in this paper facilitates a dimension independent implementation (based on C++ templates) of scientific codes using both an MPI and a hybrid approach to parallelization. We describe the essential ingredients of our implementation.
Recommendations
- A performance comparison of semi-Lagrangian discontinuous Galerkin and spline based Vlasov solvers in four dimensions
- Discontinuous Galerkin semi-Lagrangian method for Vlasov-Poisson
- A high order semi-Lagrangian discontinuous Galerkin method for Vlasov-Poisson simulations without operator splitting
- Efficient implementation of ADER discontinuous Galerkin schemes for a scalable hyperbolic PDE engine
- Performance and scalability improvements for discontinuous Galerkin solutions to conservation laws on unstructured grids
Cites work
- scientific article; zbMATH DE number 852535 (Why is no real title available?)
- scientific article; zbMATH DE number 3345258 (Why is no real title available?)
- A parallel Vlasov solver based on local cubic spline interpolation on patches
- A splitting approach for the Kadomtsev-Petviashvili equation
- Conservative semi-Lagrangian schemes for Vlasov equations
- Convergence analysis of a discontinuous Galerkin/Strang splitting approximation for the Vlasov-Poisson equations
- Discontinuous Galerkin semi-Lagrangian method for Vlasov-Poisson
- Numerical study of blow up and stability of solutions of generalized Kadomtsev-Petviashvili equations
- Numerical study of oscillatory regimes in the Kadomtsev-Petviashvili equation
- On the error propagation of semi-Lagrange and Fourier methods for advection problems
- Positivity preserving semi-Lagrangian discontinuous Galerkin formulation: theoretical analysis and application to the Vlasov-Poisson system
Cited in
(14)- The numerical flow iteration for the Vlasov-Poisson equation
- A Low-Rank Projector-Splitting Integrator for the Vlasov--Poisson Equation
- A hybrid semi-Lagrangian DG and ADER-DG solver on a moving mesh for Liouville's equation of geometrical optics
- A robust second-order low-rank BUG integrator based on the midpoint rule
- Dynamical Low-Rank Integrator for the Linear Boltzmann Equation: Error Analysis in the Diffusion Limit
- Efficient 6D Vlasov simulation using the dynamical low-rank framework \texttt{Ensign}
- Suppressing instability in a Vlasov-Poisson system by an external electric field through constrained optimization
- A low-rank projector-splitting integrator for the Vlasov-Maxwell equations with divergence correction
- A performance comparison of semi-Lagrangian discontinuous Galerkin and spline based Vlasov solvers in four dimensions
- A split step Fourier/discontinuous Galerkin scheme for the Kadomtsev-Petviashvili equation
- Semi-Lagrangian Vlasov simulation on GPUs
- Exponential methods for solving hyperbolic problems with application to collisionless kinetic equations
- A quasi-conservative Dynamical Low-rank algorithm for the Vlasov equation
- Performance and scalability improvements for discontinuous Galerkin solutions to conservation laws on unstructured grids
This page was built for publication: High performance computing aspects of a dimension independent semi-Lagrangian discontinuous Galerkin code
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q339403)
