Moment preserving constrained resampling with applications to particle-in-cell methods
From MaRDI portal
Publication:778280
DOI10.1016/J.JCP.2020.109317zbMATH Open1435.76100arXiv1702.05198OpenAlexW3006460058MaRDI QIDQ778280FDOQ778280
Authors: D. Kharzeev 
Publication date: 2 July 2020
Published in: Journal of Computational Physics (Search for Journal in Brave)
Abstract: In simulations of partial differential equations using particle-in-cell (PIC) methods, it is often advantageous to resample the particle distribution function to increase simulation accuracy, reduce compute cost, and/or avoid numerical instabilities. We introduce an algorithm for particle resampling called Moment Preserving Contrained Resampling (MPCR). The general algorithm partitions the system space into smaller subsets and is designed to conserve any number of particle and grid quantities with a high degree of accuracy (i.e. machine accuracy). The resampling scheme can be integrated into any PIC code. The advantages of MPCR, including performance, accuracy, and stability, are presented by examining several numerical tests, including a use-case study in gyrokinetic fusion plasma simulations. The tests demonstrate that while the computational cost of MPCR is negligible compared to the nascent particle evolution in PIC methods, periodic particle resampling yields a significant improvement in the accuracy and stability of the results.
Full work available at URL: https://arxiv.org/abs/1702.05198
Recommendations
- Method to increase the simulation speed of particle-in-cell (PIC) code
- Particle merging algorithm for PIC codes
- A Particle-in-cell Method with Adaptive Phase-space Remapping for Kinetic Plasmas
- On algorithms implementing the particles method on unstructured grids
- A particle-in-cell method for the simulation of plasmas based on an unconditionally stable field solver
Stochastic analysis applied to problems in fluid mechanics (76M35) Ionized gas flow in electromagnetic fields; plasmic flow (76X05)
Cites Work
- A numerically stable dual method for solving strictly convex quadratic programs
- Title not available (Why is that?)
- Title not available (Why is that?)
- On a finite-element method for solving the three-dimensional Maxwell equations
- An improved control-variate scheme for particle-in-cell simulations with collisions
- Dynamic and selective control of the number of particles in kinetic plasma simulations
- A new method for coalescing particles in PIC codes.
- Particle rezoning for multidimensional kinetic particle-in-cell simulations
- A new hybrid-Lagrangian numerical scheme for gyrokinetic simulation of tokamak edge plasma
- Adaptive particle management in a particle-in-cell code
- High-order nodal discontinuous Galerkin particle-in-cell method on unstructured grids
- Two statistical particle split and merge methods for particle-in-cell codes
- Particle merging algorithm for PIC codes
- Controlling the weights of simulation particles: adaptive particle management using \(k\)-d trees
- A finite volume method for the approximation of Maxwell's equations in two space dimensions on arbitrary meshes
- Control of the number of particles in fluid and MHD particle in cell methods
- A split control variate scheme for PIC simulations with collisions
- Nonlinear gyrokinetic equations for tokamak microturbulence
- Higher-dimensional Voronoi diagrams in linear expected time
- An Efficient Method for Weighted Sampling without Replacement
- Electromagnetic gyrokinetic PIC simulation with an adjustable control variates method
- A maximum entropy method for particle filtering
- Optimal SIR algorithm vs. fully adapted auxiliary particle filter: a non asymptotic analysis
- Tokamaks. With contributions from D. J. Campbell, J. W. Connor, R. D. Gill, N. C. Hawkes, J. Hugill, C. N. Lashmore-Davies, G. M. McCracken, H. R. Wilson, A. E. Costley, R. J. Hastie, A. Hermann, A. S. Kaye, B. Lloyd, G. F. Matthews, J. J. O'Rourke, C. M. Roach, S. E. Sharapov, D. F. Start, B. J. D. Tubbing and D. J. Ward.
- A 4th-order particle-in-cell method with phase-space remapping for the Vlasov-Poisson equation
Cited In (4)
- An unsupervised machine-learning checkpoint-restart algorithm using Gaussian mixtures for particle-in-cell simulations
- Hedging direct simulation Monte Carlo bets via event splitting
- FLEKS: a flexible particle-in-cell code for multi-scale plasma simulations
- A performant energy-conserving particle reweighting method for particle-in-cell simulations
Uses Software
This page was built for publication: Moment preserving constrained resampling with applications to particle-in-cell methods
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q778280)
