Efficient 6D Vlasov simulation using the dynamical low-rank framework Ensign
arXiv:2110.13481 · doi:10.1016/j.cpc.2022.108489
Abstract
Running kinetic simulations using grid-based methods is extremely expensive due to the up to six-dimensional phase space. Recently, it has been shown that dynamical low-rank algorithms can drastically reduce the required computational effort, while still accurately resolving important physical features such as filamentation and Landau damping. In this paper, we propose a new second order projector-splitting dynamical low-rank algorithm for the full six-dimensional Vlasov--Poisson equations. An exponential integrator based Fourier spectral method is employed to obtain a numerical scheme that is CFL condition free but still fully explicit. The resulting method is implemented with the aid of Ensign, a software framework which facilitates the efficient implementation of dynamical low-rank algorithms on modern multi-core CPU as well as GPU based systems. Its usage and features are briefly described in the paper as well. The presented numerical results demonstrate that 6D simulations can be run on a single workstation and highlight the significant speedup that can be obtained using GPUs.
References in corpus (12)
- A semi-Lagrangian Vlasov solver in tensor train format
- A low-rank method for two-dimensional time-dependent radiation transport calculations
- An asymptotic-preserving dynamical low-rank method for the multi-scale multi-dimensional linear transport equation
- A mass, momentum, and energy conservative dynamical low-rank scheme for the Vlasov equation
- A low-rank projector-splitting integrator for the Vlasov--Maxwell equations with divergence correction
- A comparison of semi-Lagrangian discontinuous Galerkin and spline based Vlasov solvers in four dimensions
- A high-order / low-order (HOLO) algorithm for preserving conservation in time-dependent low-rank transport calculations
- A Low Rank Tensor Representation of Linear Transport and Nonlinear Vlasov Solutions and Their Associated Flow Maps
- Dynamical low-rank integrator for the linear Boltzmann equation: error analysis in the diffusion limit
- High performance computing aspects of a dimension independent semi-Lagrangian discontinuous Galerkin code
- An almost symmetric Strang splitting scheme for nonlinear evolution equations
- An almost symmetric Strang splitting scheme for the construction of high order composition methods