An Energy- and Charge-conserving, Implicit, Electrostatic Particle-in-Cell Algorithm

Details An Energy- and Charge-conserving, Implicit, Electrostatic Particle-in-Cell Algorithm An Energy- and Charge-conserving, Implicit, Electrostatic Particle-in-Cell Algorithm

2011-01-19

arxiv.org/abs/1101.3701v1

J. Comput. Phys. 230, 18 (2011), p. 7018

Physics

Computational Physics

29 pages, 8 figures, submitted to Journal of Computational Physics

Scientific paper

10.1016/j.jcp.2011.05.031

This paper discusses a novel fully implicit formulation for a 1D electrostatic particle-in-cell (PIC) plasma simulation approach. Unlike earlier implicit electrostatic PIC approaches (which are based on a linearized Vlasov-Poisson formulation), ours is based on a nonlinearly converged Vlasov-Amp\`ere (VA) model. By iterating particles and fields to a tight nonlinear convergence tolerance, the approach features superior stability and accuracy properties, avoiding most of the accuracy pitfalls in earlier implicit PIC implementations. In particular, the formulation is stable against temporal (CFL) and spatial (aliasing) instabilities. It is charge- and energy-conserving to numerical roundoff for arbitrary implicit time steps. While momentum is not exactly conserved, errors are kept small by an adaptive particle sub-stepping orbit integrator, which is instrumental to prevent particle tunneling. The VA model is orbit-averaged along particle orbits to enforce an energy conservation theorem with particle sub-stepping. As a result, very large time steps, constrained only by the dynamical time scale of interest, are possible without accuracy loss. Algorithmically, the approach features a Jacobian-free Newton-Krylov solver. A main development in this study is the nonlinear elimination of the new-time particle variables (positions and velocities). Such nonlinear elimination, which we term particle enslavement, results in a nonlinear formulation with memory requirements comparable to those of a fluid computation, and affords us substantial freedom in regards to the particle orbit integrator. Numerical examples are presented that demonstrate the advertised properties of the scheme. In particular, long-time ion acoustic wave simulations show that numerical accuracy does not degrade even with very large implicit time steps, and that significant CPU gains are possible.

Affiliated with

Also associated with

No associations

Rating

An Energy- and Charge-conserving, Implicit, Electrostatic Particle-in-Cell Algorithm does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.

If you have personal experience with An Energy- and Charge-conserving, Implicit, Electrostatic Particle-in-Cell Algorithm, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and An Energy- and Charge-conserving, Implicit, Electrostatic Particle-in-Cell Algorithm will most certainly appreciate the feedback.

Rate now

Comments { 0 }

Profile ID: LFWR-SCP-O-157023