Skip to content
Preprint

High-Order Structure-Preserving SBP Finite Difference Methods for the Vlasov-Maxwell System on Matrix-Free GPUs

Sep 2026 · 0 citations · 24 references
Mathematics Computer Science Physics

TL;DR

A high-order, stable summation-by-parts (SBP) finite difference method for solving the Vlasov-Maxwell system in a 2D2V phase space and proves that the fully discrete scheme exactly conserves mass and preserves momentum up to truncation error.

Abstract

In this paper, we present a high-order, stable summation-by-parts (SBP) finite difference method for solving the Vlasov-Maxwell system in a 2D2V phase space. Central SBP operators for the advection terms are not stable when the solution becomes non-smooth and fine-scale filamentary structures develop, as is typical in high-dimensional Vlasov-Maxwell simulations. To address this issue, the method is stabilized using high-order upwind SBP operators. High-order explicit Runge-Kutta methods are employed for time integration. We prove that the fully discrete scheme exactly conserves mass and preserves momentum up to truncation error. Furthermore, we present a matrix-free implementation of the method on modern GPU architectures. A range of challenging benchmark problems is solved to demonstrate the accuracy, robustness, and performance of the proposed scheme.

View source

Similar papers

Preprint Sep 2026

A second-order structure- and positivity-preserving convex limiting method for the Vlasov equations

In this paper, we introduce a novel second-order, positivity-preserving finite element method for the Vlasov equations using a convex limiting algorithm. The method employs strong-stability-preserving (SSP) Runge-Kutta time integration and a tensor-product construction of the phase-space mesh for efficient high-dimensi...

K. Kormann, Murtazo Nazarov, Wen Junjie · 1 citation
Preprint Aug 2026

Sweep-based, implicit solutions of the multidimensional BGK equation on unstructured grids

We present a nodal discontinuous Galerkin method for solving the Bhatnagar-Gross-Krook (BGK) kinetic equation on multi-dimensional, unstructured grids. The method uses implicit, sweep-based solvers and a moment-preserving projection of the Maxwellian source to enable high-order accuracy in time while avoiding restricti...

T. Evans, Ryan S. Glasby, Cory D. Hauck et al. · 1 citation
Preprint Sep 2026

Mixed-precision GPU algorithms for efficient turbulent flow simulations with Raviart-Thomas finite elements

A solution framework with fast matrix-free operator evaluation for all ingredients, combined with multigrid solvers for the Poisson problem, and a robust mixed-precision algorithmic framework to generate accurate initial guesses for the iterative linear solvers are developed.

Ivan Prusak, Enes Mustafa Soydan, Ivan Pribec et al. · 0 citations
Preprint Sep 2026

An energy stable and accuracy-preserving finite volume scheme based on the SAV method with application to wall-distance computation

A novel semi-implicit second-order finite volume scheme integrating the scalar auxiliary variable (SAV) approach is proposed for solving the pseudo-time Eikonal equation in wall-distance computation. Unconditional energy stability under zero boundary conditions is rigorously proved, eliminating the dependence of the ti...

Xiao-Rui Xu, Qian Wang · 0 citations
Preprint Sep 2026

Positivity-preserving semi-Lagrangian discontinuous Galerkin methods for multi-species Vlasov-Amp\`ere models of plasma

We develop a high-order and positivity-preserving kinetic solver for the multi-species Vlasov-Amp\`ere system that is built within a semi-Lagrangian discontinuous Galerkin (SLDG) framework. The method is designed for efficiency and accuracy across a broad range of conditions, including stiff, multi-scale regimes with s...

D. Gambo, James A. Rossmanith · 0 citations
Sep 2026

An efficient spectral element phase-field method for incompressible multiphase flow simulation

This study develops a graphics processing unit (GPU)-accelerated spectral element phase-field method for efficient and accurate simulation of incompressible two-phase flows. The objective is to reduce the cost of the repeated elliptic subproblems that dominate Navier–Stokes–Cahn–Hilliard simulations while retaining hig...

Yao Xiao, Li-Ai Han, Ji Yao et al. · 0 citations

We use cookies to run the site and, with your consent, for analytics and to show ads. See our Cookie Policy.