A low-rank projector-splitting integrator for the Vlasov--Maxwell equations with divergence correction
arXiv:1902.00424 · doi:10.1016/j.jcp.2019.109063
Abstract
The Vlasov--Maxwell equations are used for the kinetic description of magnetized plasmas. As they are posed in an up to 3+3 dimensional phase space, solving this problem is extremely expensive from a computational point of view. In this paper, we exploit the low-rank structure in the solution of the Vlasov equation. More specifically, we consider the Vlasov--Maxwell system and propose a dynamic low-rank integrator. The key idea is to approximate the dynamics of the system by constraining it to a low-rank manifold. This is accomplished by a projection onto the tangent space. There, the dynamics is represented by the low-rank factors, which are determined by solving lower-dimensional partial differential equations. The proposed scheme performs well in numerical experiments and succeeds in capturing the main features of the plasma dynamics. We demonstrate this good behavior for a range of test problems. The coupling of the Vlasov equation with the Maxwell system, however, introduces additional challenges. In particular, the divergence of the electric field resulting from Maxwell's equations is not consistent with the charge density computed from the Vlasov equation. We propose a correction based on Lagrange multipliers which enforces Gauss' law up to machine precision.
References in corpus (8)
- Time integration of tensor trains
- A Hamiltonian splitting for the Vlasov-Maxwell system
- A semi-Lagrangian Vlasov solver in tensor train format
- Block-Structured Adaptive Mesh Refinement Algorithms for Vlasov Simulation
- A dynamical adaptive tensor method for the Vlasov-Poisson system
- High performance computing aspects of a dimension independent semi-Lagrangian discontinuous Galerkin code
- Relativistic Vlasov-Maxwell modelling using finite volumes and adaptive mesh refinement
- An asymptotic preserving scheme for the relativistic Vlasov--Maxwell equations in the classical limit
Cited by in corpus (16)
- 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 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
- A quantum-inspired method for solving the Vlasov-Poisson equations
- A robust and conservative dynamical low-rank algorithm
- A low-rank power iteration scheme for neutron transport criticality problems
- Semi-Lagrangian Vlasov simulation on GPUs
- A Parallel Low-Rank Solver for the Six-Dimensional Vlasov-Maxwell Equations
- Quantized tensor networks for solving the Vlasov-Maxwell equations
- A sweep-based low-rank method for the discrete ordinate transport equation
- Efficient 6D Vlasov simulation using the dynamical low-rank framework Ensign
- Dynamical low-rank approximation of the Vlasov-Poisson equation with piecewise linear spatial boundary
- An efficient dynamical low-rank algorithm for the Boltzmann-BGK equation close to the compressible viscous flow regime
- Dynamical low-rank tensor approximations to high-dimensional parabolic problems: existence and convergence of spatial discretizations