Rigorous charge conservation for local electromagnetic field solvers

doi:10.1016/0010-4655(92)90169-Y

Computer Physics Communications

Volume 69, Issues 2–3, March–April 1992, Pages 306-316

https://doi.org/10.1016/0010-4655(92)90169-Y Get rights and content

Abstract

In this paper we present a rigorous method for finding the electric field in two-and-one-half and three dimensional electromagnetic plasma simulations without resorting to computationally expensive transforms. A finite grid interpretation of the divergence equation $∇· J = -ϖρ/ϖt$ is offered which allows the current density and thus new local electric and magnetic field strengths to be determined directly from knowledge of charge motion.

References (7)

B. Marder
J. Comput. Phys.
(1987)
E. Lindman
J. Comput. Phys.
(1975)
O. Buneman et al.
Relativistic Plasmas

There are more references available in the full text version of this article.

Cited by (398)

GPIC: A set of high-efficiency CUDA Fortran code using gpu for particle-in-cell simulation in space physics
2024, Computer Physics Communications
Previous implementation of the full kinetic particle-in-cell (PIC) simulation on the Graphical Processing Unit (GPU) device has shown advantages compared to traditional Central Processing Unit (CPU) computing. However, conventional PIC simulations using GPU computing have faced limitations in terms of low performance when the simulation box size or the number of particles per cell increases. In this study, we boost computing efficiency by designing novel schemes of kernels with the combination of numerical and technical aspects. The time expenses of the global data transferring and duplicate data fetching processes have been massively reduced by the utilization of the on-chip memory. The reduction treatment on the particles, the graded current computing scheme, and the 2.5D thread launch strategy are designed to accelerate the simulation iterations. The new scheme can reach up to about 5.5 times the acceleration ratio than the old one and attain the highest to 734 times faster than the program using CPUs. Our new scheme can realize large-scale PIC simulations with both decent performances on the legacy GPUs and the increasing trend of the acceleration ratio on the lasted GPUs.
An energy-conserving Fourier particle-in-cell method with asymptotic-preserving preconditioner for Vlasov-Ampère system with exact curl-free constraint
2023, Journal of Computational Physics
We present an efficient and accurate energy-conserving implicit particle-in-cell (PIC) algorithm for the electrostatic Vlasov system, with particular emphasis on its high robustness for simulating complex plasma systems with multiple physical scales. This method consists of several indispensable elements: (i) the reformulation of the original Vlasov-Poisson system into an equivalent Vlasov-Ampère system with divergence-/curl-free constraints; (ii) a novel structure-preserving Fourier spatial discretization, which exactly preserves these constraints at the discrete level; (iii) a preconditioned Anderson-acceleration algorithm for the solution of the highly nonlinear system; and (iv) a linearized and uniform approximation of the implicit Crank-Nicolson scheme for various Debye lengths, based on the generalized Ohm's law, which serves as an asymptotic-preserving preconditioner for the proposed method. Numerical experiments are conducted, and comparisons are made among the proposed energy-conserving scheme, the classical leapfrog scheme, and a Strang operator-splitting scheme to demonstrate the superiority of the proposed method, especially for plasma systems crossing physical scales.
An implicit particle code with exact energy and charge conservation for electromagnetic studies of dense plasmas
2023, Journal of Computational Physics
A collisional particle code based on implicit energy- and charge-conserving methods is presented. A modified version of the particle-suppressed Jacobian-Free Newton-Krylov method that can enhance the solver efficiency is introduced. Mathematically, it is shown that this new approach can be viewed as a fixed-point iteration method for the particle positions. The model can exactly conserve global energy and local charge and can efficiently use time steps larger than the plasma period. The algorithm's ability to simulate dense plasmas accurately and efficiently is quantified by simulating the dynamic compression of a plasma slab via a magnetic piston in 1D planar geometry.
EZ: An efficient, charge conserving current deposition algorithm for electromagnetic particle-in-cell simulations
2023, Computer Physics Communications
We present EZ, a novel current deposition algorithm for particle-in-cell (PIC) simulations. EZ calculates the current density on the electromagnetic grid due to macro-particle motion within a time step by solving the continuity equation of electrodynamics. Being a charge conserving hybridization of Esirkepov's method and ZigZag, we refer to it as “EZ” as shorthand for “Esirkepov meets ZigZag”. Simulations of a warm, relativistic plasma with PIConGPU show that EZ achieves the same level of charge conservation as the commonly used method by Esirkepov, yet reaches higher performance for macro-particle assignment-functions up to third-order. In addition to a detailed description of the functioning of EZ, reasons for the expected and observed performance increase are given, and guidelines for its implementation aiming at highest performance on GPUs are provided.
A pseudospectral implicit particle-in-cell method with exact energy and charge conservation
2023, Computer Physics Communications
The standard particle-in-cell (PIC) method employs explicit finite-difference (FD) schemes, such as the leap-frog scheme, for spatial and temporal integrations. In this study, we introduce a novel approach by employing a pseudospectral method for solving the Poisson equation and adopting a fully implicit time integration technique to achieve exact energy conservation. The pseudospectral field solver offers the advantage of spectral accuracy in obtaining field solutions. Previous investigations on implicit time integration of PIC FD equations have demonstrated the capability of enforcing precise energy exchange between the field and particles, resulting in energy-conserving schemes. We prove that the exact energy conservation property can be extended to the pseudospectral scheme. Additionally, we provide a solution to ensure a pseudospectral charge continuity equation. To validate our methodology, we implement the new scheme in a 2D electrostatic PIC code and present numerical examples that confirm the theoretical findings.
Verification and benchmarking relativistic electron beam transport through a background gas
2023, Computer Physics Communications
It is necessary to establish confidence in high-consequence codes containing an extensive suite of physics algorithms in the regimes of interest. Verification problems allow code developers to assess numerical accuracy and increase confidence that specific sets of model physics were implemented correctly in the code. The two main verification techniques are code verification and solution verification. In this work, we present verification problems that can be used in other codes to increase confidence in simulations of relativistic beam transport. Specifically, we use the general plasma code EMPIRE to model and compare with the analytical solution to the evolution of the outer radial envelope of a relativistic charged particle beam. We also outline a benchmark test of a relativistic beam propagating through a vacuum and pressurized gas cell, and present the results between EMPIRE and the hybrid code GAZEL. Further, we discuss the subtle errors that were caught with these problems and detail lessons learned.

View all citing articles on Scopus

View full text

Rigorous charge conservation for local electromagnetic field solvers

Abstract

J. Comput. Phys.

J. Comput. Phys.

Relativistic Plasmas