A Multi Level Multi Domain Method for Particle In Cell Plasma Simulations
arXiv:1201.6208 · doi:10.1016/j.jcp.2012.12.028
Abstract
A novel adaptive technique for electromagnetic Particle In Cell (PIC) plasma simulations is presented here. Two main issues are identified in designing adaptive techniques for PIC simulation: first, the choice of the size of the particle shape function in progressively refined grids, with the need to avoid the exertion of self-forces on particles, and, second, the necessity to comply with the strict stability constraints of the explicit PIC algorithm. The adaptive implementation presented responds to these demands with the introduction of a Multi Level Multi Domain (MLMD) system (where a cloud of self-similar domains is fully simulated with both fields and particles) and the use of an Implicit Moment PIC method as baseline algorithm for the adaptive evolution. Information is exchanged between the levels with the projection of the field information from the refined to the coarser levels and the interpolation of the boundary conditions for the refined levels from the coarser level fields. Particles are bound to their level of origin and are prevented from transitioning to coarser levels, but are repopulated at the refined grid boundaries with a splitting technique. The presented algorithm is tested against a series of simulation challenges.
References in corpus (1)
Cited by in corpus (10)
- Relativistic magnetic reconnection in collisionless ion-electron plasmas explored with particle-in-cell simulations
- Arbitrarily high order Convected Scheme solution of the Vlasov-Poisson system
- Note on quantitatively correct simulations of the kinetic beam-plasma instability
- Simflowny 3: An upgraded platform for scientific modelling and simulation
- The muphyII Code: Multiphysics Plasma Simulation on Large HPC Systems
- A new method to dispatch split particles in Particle-In-Cell codes
- Parallelization of an implicit algorithm for multi-dimensional particle-in-cell simulations
- Comparison of Integrators for Electromagnetic Particle in Cell Methods: Algorithms and Applications
- Optimal Landau-type closure parameters for two-fluid simulations of plasma turbulence at kinetic scales
- An exact particle splitting algorithm for PIC codes