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

ECsim-CYL: Energy Conserving Semi-Implicit particle in cell simulation in axially symmetric cylindrical coordinates

178   0   0.0 ( 0 )
 نشر من قبل Giovanni Lapenta
 تاريخ النشر 2018
  مجال البحث فيزياء
والبحث باللغة English




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

Based on the previously developed Energy Conserving Semi Implicit Method (ECsim) code, we present its cylindrical implementation, called ECsim-CYL, to be used for axially symmetric problems. The main motivation for the development of the cylindrical version is to greatly improve the computational speed by utilizing cylindrical symmetry. The ECsim-CYL discretizes the field equations in two-dimensional cylindrical coordinates using the finite volume method . For the particle mover, it uses a modification of ECsims mover for cylindrical coordinates by keeping track of all three components of velocity vectors, while only keeping radial and axial coordinates of particle positions. In this paper, we describe the details of the algorithm used in the ECsim-CYL and present a series of tests to validate the accuracy of the code including a wave spectrum in a homogeneous plasmas inside a cylindrical waveguide and free expansion of a spherical plasma ball in vacuum. The ECsim-CYL retains the stability properties of ECsim and conserves the energy within machine precision, while accurately describing the plasma behavior in the test cases.



قيم البحث

اقرأ أيضاً

The recently developed energy conserving semi-implicit method (ECsim) for PIC simulation is applied to multiple scale problems where the electron-scale physics needs to be only partially retained and the interest is on the macroscopic or ion-scale pr ocesses. Unlike hybrid methods, the ECsim is capable of providing kinetic electron information, such as wave-electron interaction (Landau damping or cyclotron resonance) and non-Maxwellian electron velocity distributions. However, like hybrid, the ECsim does not need to resolve all electron scales, allowing time steps and grid spacing orders of magnitude larger than in explicit PIC schemes. The additional advantage of the ECsim is that the stability at large scale is obtained while conserving energy exactly. Three examples are presented: ion acoustic waves, electron acoustic instability and reconnection processes.
We present in this work the implementation of the Energy Conserving Semi-Implicit Method in a parallel code called ECsim. This new code is a three-dimensional, fully electromagnetic particle in cell (PIC) code. It is written in C/C++ and uses MPI to allow massive parallelization. ECsim is unconditionally stable in time, eliminates the finite grid instability, has the same cycle scheme as the explicit code with a computational cost comparable to other semi-implicit PIC codes. All this features make it a very valuable tool to address situations which have not been possible to analyze until now with other PIC codes. In this work, we show the details of the algorithm implementation and we study its performance in different systems. ECsim is compared with another semi-implicit PIC code with different time and spectral resolution, showing its sability to address situations where other codes fail.
158 - Tsunehiko N. Kato 2013
When a charged particle moves through a plasma at a speed much higher than the thermal velocity of the plasma, it is subjected to the force of the electrostatic field induced in the plasma by itself and loses its energy. This process is well-known as the stopping power of a plasma. In this paper we show that the same process works in particle-in-cell (PIC) simulations as well and the energy loss rate of fast particles due to this process is mainly determined by the number of plasma electrons contained in the electron skin depth volume. However, since there are generally very few particles in that volume in PIC simulations compared with real plasmas, the energy loss effect can be exaggerated significantly and can affect the results. Therefore, especially for the simulations that investigate the particle acceleration processes, the number of particles used in the simulations should be chosen large enough to avoid this artificial energy loss.
We present a new technique for transferring momentum and velocity between particles and grid with Particle-In-Cell (PIC) calculations which we call Affine-Particle-In-Cell (APIC). APIC represents particle velocities as locally affine, rather than loc ally constant as in traditional PIC. We show that this representation allows APIC to conserve linear and angular momentum across transfers while also dramatically reducing numerical diffusion usually associated with PIC. Notably, conservation is achieved with lumped mass, as opposed to the more commonly used Fluid Implicit Particle (FLIP) transfers which require a full mass matrix for exact conservation. Furthermore, unlike FLIP, APIC retains a filtering property of the original PIC and thus does not accumulate velocity modes on particles as FLIP does. In particular, we demonstrate that APIC does not experience velocity instabilities that are characteristic of FLIP in a number of Material Point Method (MPM) hyperelasticity calculations. Lastly, we demonstrate that when combined with the midpoint rule for implicit update of grid momentum that linear and angular momentum are exactly conserved.
We construct Boris-type schemes for integrating the motion of charged particles in particle-in-cell (PIC) simulation. The new solvers virtually combine the 2-step Boris procedure arbitrary n times in the Lorentz-force part, and therefore we call them the multiple Boris solvers. Using Chebyshev polynomials, a one-step form of the new solvers is provided. The new solvers give n^2 times smaller errors, allow larger timesteps, and have a long-term stability. We present numerical tests of the new solvers, in comparison with other particle integrators.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

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