Related papers: A Local Macroscopic Conservative (LoMaC) low rank …
We extend our previous work on a semi-Lagrangian adaptive rank (SLAR) integrator, in the finite difference framework for nonlinear Vlasov-Poisson systems, to the general high-order tensor setting. The proposed scheme retains the high-order…
In this paper, we consider a finite difference grid-based semi-Lagrangian approach in solving the Vlasov-Poisson (VP) system. Many of existing methods are based on dimensional splitting, which decouples the problem into solving linear…
Tensor methods are among the most prominent tools for the numerical solution of high-dimensional problems where functions of multiple variables have to be approximated. These methods exploit the tensor structure of function spaces and apply…
High-resolution simulations of particle-based kinetic plasma models typically require a high number of particles and thus often become computationally intractable. This is exacerbated in multi-query simulations, where the problem depends on…
Density varies spatiotemporally in low Mach number flows. Hence, incompressibility cannot be assumed, and the density must be accurately solved. Various methods have been proposed to analyze low Mach number flows, but their energy…
A popular approach to analyze the dynamics of high-dimensional many-body systems, such as macromolecules, is to project the trajectories onto a space of slowly-varying collective variables, where subsequent analyses are made, such as…
Recent development of structure-preserving geometric particle-in-cell (PIC) algorithms for Vlasov-Maxwell systems is summarized. With the arriving of 100 petaflop and exaflop computing power, it is now possible to carry out direct…
Lorentz invariant structure-preserving algorithms possess reference-independent secular stability, which is vital for simulating relativistic multi-scale dynamical processes. The splitting method has been widely used to construct…
We propose a class of Particle-In-Cell (PIC) methods for the Vlasov-Poisson system with a strong and inhomogeneous external magnetic field with fixed direction, where we focus on the motion of particles in the plane orthogonal to the…
Many applications in data science and scientific computing involve large-scale datasets that are expensive to store and compute with, but can be efficiently compressed and stored in an appropriate tensor format. In recent years, randomized…
In this paper, a novel high order semi-Lagrangian (SL) spectral volume (SV) method is proposed and studied for nonlinear Vlasov-Poisson (VP) simulations via operator splitting. The proposed algorithm combines both advantages of…
Despite advancements in high-performance computing and modern numerical algorithms, computational cost remains prohibitive for multi-query kinetic plasma simulations. In this work, we develop data-driven reduced-order models (ROMs) for…
We construct high order symmetric volume-preserving methods for the relativistic dynamics of a charged particle by the splitting technique with processing. Via expanding the phase space to include time $t$, we give a more general…
The Nordstr\"om-Vlasov system is a Lorentz invariant model for a self-gravitating collisionless gas. We establish suitable a-priori-bounds on the solutions of this system, which together with energy estimates and the smoothing effect of…
In this work, we aim at constructing numerical schemes, that are as efficient as possible in terms of cost and conservation of invariants, for the Vlasov--Fokker--Planck system coupled with Poisson or Amp\`ere equation. Splitting methods…
Fragmentation methods applied to multireference wave functions constitute a road towards the application of highly accurate ab initio wave function calculations to large molecules and solids. However, it is important for reproducibility and…
We show that nonrelativsitic scaling of the collisionless Vlasov-Maxwell system implies the existence of a formal invariant slow manifold in the infinite-dimensional Vlasov-Maxwell phase space. Vlasov-Maxwell dynamics restricted to the slow…
We introduce an algorithm to find possible constants of motion for a given replicator equation. The algorithm is inspired by Dirac geometry and a Hamiltonian description for the replicator equations with such constants of motion, up to a…
We present an efficient and accurate energy-conserving implicit particle-in-cell~(PIC) algorithm for the electrostatic Vlasov system, with particular emphasis on its high robustness for simulating complex plasma systems with multiple…
This paper investigates nonlinear Landau damping in the 3D Vlasov-Poisson (VP) system. We study the asymptotic stability of the Poisson equilibrium $\mu(v)=\frac{1}{\pi^2(1+|v|^2)^2}$ under small perturbations. Building on the foundational…