Fully implicit Particle-in-cell algorithms for multiscale plasma simulation
- Los Alamos National Laboratory
The outline of the paper is as follows: Particle-in-cell (PIC) methods for fully ionized collisionless plasmas, explicit vs. implicit PIC, 1D ES implicit PIC (charge and energy conservation, moment-based acceleration), and generalization to Multi-D EM PIC: Vlasov-Darwin model (review and motivation for Darwin model, conservation properties (energy, charge, and canonical momenta), and numerical benchmarks). The author demonstrates a fully implicit, fully nonlinear, multidimensional PIC formulation that features exact local charge conservation (via a novel particle mover strategy), exact global energy conservation (no particle self-heating or self-cooling), adaptive particle orbit integrator to control errors in momentum conservation, and canonical momenta (EM-PIC only, reduced dimensionality). The approach is free of numerical instabilities: ωpeΔt >> 1, and Δx >> λD. It requires many fewer dofs (vs. explicit PIC) for comparable accuracy in challenging problems. Significant CPU gains (vs explicit PIC) have been demonstrated. The method has much potential for efficiency gains vs. explicit in long-time-scale applications. Moment-based acceleration is effective in minimizing NFE, leading to an optimal algorithm.
- Research Organization:
- Los Alamos National Laboratory (LANL)
- Sponsoring Organization:
- LDRD
- DOE Contract Number:
- AC52-06NA25396
- OSTI ID:
- 1196172
- Report Number(s):
- LA-UR-15-25379
- Country of Publication:
- United States
- Language:
- English
Similar Records
Multirate Particle-in-Cell Time Integration Techniques of Vlasov-Maxwell Equations for Collisionless Kinetic Plasma Simulations
A curvilinear, fully implicit, conservative electromagnetic PIC algorithm in multiple dimensions