Log In Sign Up

Adaptive symplectic model order reduction of parametric particle-based Vlasov-Poisson equation

by   Jan S. Hesthaven, et al.

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 a set of parameters. In this work, we derive reduced order models for the semi-discrete Hamiltonian system resulting from a geometric particle-in-cell approximation of the parametric Vlasov-Poisson equations. Since the problem's non-dissipative and highly nonlinear nature makes it reducible only locally in time, we adopt a nonlinear reduced basis approach where the reduced phase space evolves in time. This strategy allows a significant reduction in the number of simulated particles, but the evaluation of the nonlinear operators associated with the Vlasov-Poisson coupling remains computationally expensive. We propose a novel reduction of the nonlinear terms that combines adaptive parameter sampling and hyper-reduction techniques to address this. The proposed approach allows decoupling the operations having a cost dependent on the number of particles from those that depend on the instances of the required parameters. In particular, in each time step, the electric potential is approximated via dynamic mode decomposition (DMD) and the particle-to-grid map via a discrete empirical interpolation method (DEIM). These approximations are constructed from data obtained from a past temporal window at a few selected values of the parameters to guarantee a computationally efficient adaptation. The resulting DMD-DEIM reduced dynamical system retains the Hamiltonian structure of the full model, provides good approximations of the solution, and can be solved at a reduced computational cost.


page 22

page 26


Gradient-preserving hyper-reduction of nonlinear dynamical systems via discrete empirical interpolation

This work proposes a hyper-reduction method for nonlinear parametric dyn...

Symplectic model reduction methods for the Vlasov equation

Particle-based simulations of the Vlasov equation typically require a la...

Sparse Grids based Adaptive Noise Reduction strategy for Particle-In-Cell schemes

We propose a sparse grids based adaptive noise reduction strategy for el...

Reduced-order modeling for Ablowitz-Ladik equation

In this paper, reduced-order models (ROMs) are constructed for the Ablow...

Efficient and Accurate Adaptive Resolution for Weakly-Compressible SPH

In this paper we propose an accurate, and computationally efficient meth...

Bayesian phase estimation with adaptive grid refinement

We introduce a novel Bayesian phase estimation technique based on adapti...

Sensitivity-driven adaptive sparse stochastic approximations in plasma microinstability analysis

Quantifying uncertainty in predictive simulations for real-world problem...