A semi-Lagrangian Vlasov solver in tensor train format

Abstract

In this article, we derive a semi-Lagrangian scheme for the solution of the Vlasov equation represented as a low-parametric tensor. Grid-based methods for the Vlasov equation have been shown to give accurate results but their use has mostly been limited to simulations in two dimensional phase space due to extensive memory requirements in higher dimensions. Compression of the solution via high-order singular value decomposition can help in reducing the storage requirements and the tensor train (TT) format provides efficient basic linear algebra routines for low-rank representations of tensors. In this paper, we develop interpolation formulas for a semi-Lagrangian solver in TT format. In order to efficiently implement the method, we propose a compression of the matrix representing the interpolation step and an efficient implementation of the Hadamard product. We show numerical simulations for standard test cases in two, four and six dimensional phase space. Depending on the test case, the memory requirements reduce by a factor $10^2-10^3$ in four and a factor $10^5-10^6$ in six dimensions compared to the full-grid method.

References in corpus (2)

• Alternating minimal energy methods for linear systems in higher dimensions. Part II: Faster algorithm and application to nonsymmetric systems
• Low-rank approximation in the numerical modeling of the Farley-Buneman instability in ionospheric plasma

Cited by in corpus (27)

• Vlasov methods in space physics and astrophysics
• A low-rank method for two-dimensional time-dependent radiation transport calculations
• 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 Low Rank Tensor Representation of Linear Transport and Nonlinear Vlasov Solutions and Their Associated Flow Maps
• Stable interpolation with isotropic and anisotropic Gaussians using Hermite generating function
• A quantum-inspired method for solving the Vlasov-Poisson equations
• A massively parallel semi-Lagrangian solver for the six-dimensional Vlasov-Poisson equation
• Dynamical low-rank integrator for the linear Boltzmann equation: error analysis in the diffusion limit
• A dynamical adaptive tensor method for the Vlasov-Poisson system
• Tensor methods for the Boltzmann-BGK equation
• Parallel numerical tensor methods for high-dimensional PDEs
• A Parallel Low-Rank Solver for the Six-Dimensional Vlasov-Maxwell Equations
• Semi-Lagrangian Vlasov simulation on GPUs
• Quantized tensor networks for solving the Vlasov-Maxwell equations
• Numerical solution of the Boltzmann equation with S-model collision integral using tensor decompositions
• Efficient 6D Vlasov simulation using the dynamical low-rank framework Ensign
• An Energy Conserving Vlasov Solver That Tolerates Coarse Velocity Space Resolutions: Simulation of MMS Reconnection Events
• Solving high-dimensional parabolic PDEs using the tensor train format
• Discontinuous Galerkin schemes for a class of Hamiltonian evolution equations with applications to plasma fluid and kinetic problems
• An efficient dynamical low-rank algorithm for the Boltzmann-BGK equation close to the compressible viscous flow regime
• Sparse Grid Discontinuous Galerkin Methods for the Vlasov-Maxwell System
• Sylvester-Preconditioned Adaptive-Rank Implicit Time Integrators for Advection-Diffusion Equations with Variable Coefficients
• Adaptive hyper-reduction of non-sparse operators: application to parametric particle-based kinetic plasma models
• Neural semi-Lagrangian method for high-dimensional advection-diffusion problems
• Optimal Landau-type closure parameters for two-fluid simulations of plasma turbulence at kinetic scales
• Dynamical Tensor Train Approximation for Kinetic Equations