ترغب بنشر مسار تعليمي؟ اضغط هنا

The geometric theory of charge conservation in particle-in-cell simulations

68   0   0.0 ( 0 )
 نشر من قبل Alexander S. Glasser
 تاريخ النشر 2019
  مجال البحث فيزياء
والبحث باللغة English




اسأل ChatGPT حول البحث

In recent years, several gauge-symmetric particle-in-cell (PIC) methods have been developed whose simulations of particles and electromagnetic fields exactly conserve charge. While it is rightly observed that these methods gauge symmetry gives rise to their charge conservation, this causal relationship has generally been asserted via ad hoc derivations of the associated conservation laws. In this work, we develop a comprehensive theoretical grounding for charge conservation in gauge-symmetric Lagrangian and Hamiltonian PIC algorithms. For Lagrangian variational PIC methods, we apply Noethers second theorem to demonstrate that gauge symmetry gives rise to a local charge conservation law as an off-shell identity. For Hamiltonian splitting methods, we show that the momentum map establishes their charge conservation laws. We define a new class of algorithms -- gauge-compatible splitting methods -- that exactly preserve the momentum map associated with a Hamiltonian systems gauge symmetry -- even after time discretization. This class of algorithms affords splitting schemes a decided advantage over alternative Hamiltonian integrators. We apply this general technique to design a novel, explicit, symplectic, gauge-compatible splitting PIC method, whose momentum map yields an exact local charge conservation law. Our study clarifies the appropriate initial conditions for such schemes and examines their symplectic reduction.



قيم البحث

اقرأ أيضاً

The paper provides a tutorial to the conceptual layout of a self-consistently coupled Particle-In-Cell/Test-Particle model for the kinetic simulation of sputtering transport in capacitively coupled plasmas at low gas pressures. It explains when a kin etic approach is actually needed and which numerical concepts allow for the inherent nonequilibrium behavior of the charged and neutral particles. At the example of a generic sputtering discharge both the fundamentals of the applied Monte Carlo methods as well as the conceptual details in the context of the sputtering scenario are elaborated on. Finally, two in the context of sputtering transport simulations often exploited assumptions, namely on the energy distribution of impinging ions as well as on the test particle approach, are validated for the proposed example discharge.
We construct a particle integrator for nonrelativistic particles by means of the splitting method based on the exact flow of the equation of motion of particles in the presence of constant electric and magnetic field. This integrator is volume-preser ving similar to the standard Boris integrator and is suitable for long-term integrations in particle-in-cell simulations. Numerical tests reveal that it is significantly more accurate than previous volume-preserving integrators with second-order accuracy. For example, in the $E times B$ drift test, this integrator is more accurate than the Boris integrator and the integrator based on the exact solution of gyro motion by three and two orders of magnitude, respectively. In addition, we derive approximate integrators that incur low computational cost and high-precision integrators displaying fourth- to tenth-order accuracy with the aid of the composition method. These integrators are also volume-preserving. It is also demonstrated that the Boris integrator is equivalent to the simplest case of the approximate integrators derived in this study.
This paper discusses temporally continuous and discrete forms of the speed-limited particle-in-cell (SLPIC) method first treated by Werner et al. [Phys. Plasmas 25, 123512 (2018)]. The dispersion relation for a 1D1V electrostatic plasma whose fast pa rticles are speed-limited is derived and analyzed. By examining the normal modes of this dispersion relation, we show that the imposed speed-limiting substantially reduces the frequency of fast electron plasma oscillations while preserving the correct physics of lower-frequency plasma dynamics (e.g. ion acoustic wave dispersion and damping). We then demonstrate how the timestep constraints of conventional electrostatic particle-in-cell methods are relaxed by the speed-limiting approach, thus enabling larger timesteps and faster simulations. These results indicate that the SLPIC method is a fast, accurate, and powerful technique for modeling plasmas wherein electron kinetic behavior is nontrivial (such that a fluid/Boltzmann representation for electrons is inadequate) but evolution is on ion timescales.
Upon inclusion of collisions, the speed-limited particle-in-cell (SLPIC) simulation method successfully computed the Paschen curve for argon. The simulations modelled an electron cascade across an argon-filled capacitor, including electron-neutral io nization, electron-neutral elastic collisions, electron-neutral excitation, and ion-induced secondary-electron emission. In electrical breakdown, the timescale difference between ion and electron motion makes traditional particle-in-cell (PIC) methods computationally slow. To decrease this timescale difference and speed up computation, we used SLPIC, a time-domain algorithm that limits the speed of the fastest electrons in the simulation. The SLPIC algorithm facilitates a straightforward, fully-kinetic treatment of dynamics, secondary emission, and collisions. SLPIC was as accurate as PIC, but ran up to 200 times faster. SLPIC accurately computed the Paschen curve for argon over three orders of magnitude in pressure.
181 - A. Bret , M.E. Dieckmann 2010
Particle-in-cell (PIC) simulations are widely used as a tool to investigate instabilities that develop between a collisionless plasma and beams of charged particles. However, even on contemporary supercomputers, it is not always possible to resolve t he ion dynamics in more than one spatial dimension with such simulations. The ion mass is thus reduced below 1836 electron masses, which can affect the plasma dynamics during the initial exponential growth phase of the instability and during the subsequent nonlinear saturation. The goal of this article is to assess how far the electron to ion mass ratio can be increased, without changing qualitatively the physics. It is first demonstrated that there can be no exact similarity law, which balances a change of the mass ratio with that of another plasma parameter, leaving the physics unchanged. Restricting then the analysis to the linear phase, a criterion allowing to define a maximum ratio is explicated in terms of the hierarchy of the linear unstable modes. The criterion is applied to the case of a relativistic electron beam crossing an unmagnetized electron-ion plasma.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا