scholarly journals Particle-in-cell simulation of two stream instability in the non-extensive statistics

2014 ◽  
Vol 32 (3) ◽  
pp. 399-407 ◽  
Author(s):  
Mohammad Ghorbanalilu ◽  
Elahe Abdollahzadeh ◽  
S.H. Ebrahimnazhad Rahbari

AbstractWe have performed extensive one dimensional particle-in-cell (PIC) simulations to explore generation of electrostatic waves driven by two-stream instability (TSI) that arises due to the interaction between two symmetric counterstreaming electron beams. The electron beams are considered to be cold, collisionless and magnetic-field-free in the presence of neutralizing background of static ions. Here, electrons are described by the non-extensive q-distributions of the Tsallis statistics. Results shows that the electron holes structures are different for various q values such that: (i) for q > 1 cavitation of electron holes are more visible and the excited waves were more strong (ii) for q < 1 the degree of cavitation decreases and for q = 0.5 the holes are not distinguishable. Furthermore, time development of the velocity root-mean-square (VRMS) of electrons for different q-values demonstrate that the maximum energy conversion is increased upon increasing the non-extensivity parameter q up to the values q > 1. The normalized total energy history for a arbitrary entropic index q = 1.5, approves the energy conserving in our PIC simulation.

2019 ◽  
Vol 236 ◽  
pp. 153-163 ◽  
Author(s):  
Diego Gonzalez-Herrero ◽  
Alfredo Micera ◽  
Elisabetta Boella ◽  
Jaeyoung Park ◽  
Giovanni Lapenta

2021 ◽  
Author(s):  
Seiji Zenitani ◽  
Tsunehiko Kato

&lt;div&gt; &lt;div&gt; &lt;div&gt; &lt;p&gt;&amp;#160;Particle-in-cell (PIC) simulation has long been used in theoretical plasma physics. In PIC simulation, the Boris solver is the de-facto standard for solving particle motion, and it has been used over a half century. Meanwhile, there is a continuous demand for better particle solvers. In this contribution, we introduce a family of Boris-type schemes for integrating the motion of charged particles. We call the new solvers the multiple Boris solvers. The new solvers essentially repeat the standard two-step procedure multiple times in the Lorentz-force part, and we derive a single-step form for arbitrary subcycle number &lt;em&gt;n&lt;/em&gt;. The new solvers give &lt;em&gt;n&lt;sup&gt;2&lt;/sup&gt;&lt;/em&gt; times smaller errors, allow larger timesteps, but they are computationally affordable for moderate &lt;em&gt;n&lt;/em&gt;. The multiple Boris solvers also reduce a numerical error in long-term plasma motion in a relativistic magnetized flow.&lt;/p&gt; &lt;/div&gt; &lt;/div&gt; &lt;/div&gt;&lt;p&gt;Reference:&lt;/p&gt;&lt;ul&gt;&lt;li&gt;S. Zenitani &amp; T. N. Kato, &lt;em&gt;Multiple Boris integrators for particle-in-cell simulation&lt;/em&gt;, Comput. Phys. Commun. &lt;strong&gt;247&lt;/strong&gt;, 106954, doi:10.1016/j.cpc.2019.106954 (2020)&lt;/li&gt; &lt;/ul&gt;


1993 ◽  
Vol 182 (4-6) ◽  
pp. 426-432 ◽  
Author(s):  
M. Fivaz ◽  
A. Fasoli ◽  
K. Appert ◽  
F. Skiff ◽  
T.M. Tran ◽  
...  

1999 ◽  
Vol 6 (3/4) ◽  
pp. 211-219 ◽  
Author(s):  
L. Muschietti ◽  
I. Roth ◽  
R. E. Ergun ◽  
C. W. Carlson

Abstract. Recent observations from satellites crossing regions of magnetic-field-aligned electron streams reveal solitary potential structures that move at speeds much greater than the ion acoustic/thermal velocity. The structures appear as positive potential pulses rapidly drifting along the magnetic field, and are electrostatic in their rest frame. We interpret them as BGK electron holes supported by a drifting population of trapped electrons. Using Laplace transforms, we analyse the behavior of one phase-space electron hole. The resulting potential shapes and electron distribution functions are self-consistent and compatible with the field and particle data associated with the observed pulses. In particular, the spatial width increases with increasing amplitude. The stability of the analytic solution is tested by means of a two-dimensional particle-in-cell simulation code with open boundaries. We consider a strongly magnetized parameter regime in which the bounce frequency of the trapped electrons is much less than their gyrofrequency. Our investigation includes the influence of the ions, which in the frame of the hole appear as an incident beam, and impinge on the BGK potential with considerable energy. The nonlinear structure is remarkably resilient


Author(s):  
Xiaojuan Wang ◽  
Zhanghu Hu ◽  
Younian Wang

Abstract A two-dimensional(2D) electromagnetic particle-in-cell(PIC) simulation model is proposed to study the density evolution and collective stopping of electron beams in background plasmas. We show here the formation of the multi-layer structure of the relativistic electron beam in the plasma due to the different betatron frequency from the beam front to the beam tail. Meanwhile, the nonuniformity of the longitudinal wakefield is the essential reason for the multilayer structure formation in beam phase space. The influences of beam parameters (beam radius and transverse density profile) on the formation of the multi-layer structure and collective stopping in background plasmas are also considered.


1994 ◽  
Vol 12 (2) ◽  
pp. 273-282 ◽  
Author(s):  
Glenn Joyce ◽  
Jonathan Krall ◽  
Steven Slinker

ELBA is a three-dimensional, particle-in-cell, simulation code that has been developed to study the propagation and transport of relativistic charged particle beams. The code is particularly suited to the simulation of relativistic electron beams propagating through collisionless or slightly collisional plasmas or through external electric or magnetic fields. Particle motion is followed via a coordinate “window” in the laboratory frame that moves at the speed of light. This scheme allows us to model only the immediate vicinity of the beam. Because no information can move in the forward direction in these coordinates, particle and field data can be handled in a simple way that allows for very large scale simulations. A mapping scheme has been implemented that, with corrections to Maxwell's equations, allows the inclusion of bends in the simulation system.


2016 ◽  
Vol 34 (2) ◽  
pp. 284-293 ◽  
Author(s):  
E. Khalilzadeh ◽  
J. Yazdanpanah ◽  
J. Jahanpanah ◽  
A. Chakhmachi

AbstractIn this paper, a heuristic algorithm based on particle-in-cell (PIC) simulation is introduced to investigate the harmonic generation during the ionization and formation of plasma by a non-relativistic laser field when it propagates through hydrogen atoms. The harmonic generation is considered for the radiative recombination of an ionized electron with its nearest ion. The ionization algorithm is improved by considering the Stark effect and nonzero velocity for ionized electrons. Energy conservation is evaluated during the recombination process. In our code, for the first time, Maxwell's equations are integrated for harmonic fields in a separate mesh using the artificial recombination current as a source term. The simulation results are then used to illustrate the intensity spectrum of generated fields. It is shown that the initial momentum of ionized electrons affects the harmonic spectrum because the energy of radiated photons varies with the electron energy.


Photonics ◽  
2021 ◽  
Vol 8 (6) ◽  
pp. 172
Author(s):  
Vladimir Annenkov ◽  
Evgeny Berendeev ◽  
Evgeniia Volchok ◽  
Igor Timofeev

Based on particle-in-cell simulations, we propose to generate sub-nanosecond pulses of narrowband terahertz radiation with tens of MW power using unique properties of kiloampere relativistic (2 MeV) electron beams produced by linear induction accelerators. Due to small emittance of such beams, they can be focused into millimeter and sub-millimeter spots comparable in sizes with the wavelength of THz radiation. If such a beam is injected into a plasma, it becomes unstable against the two-stream instability and excites plasma oscillations that can be converted to electromagnetic waves at the plasma frequency and its harmonics. It is shown that several radiation mechanisms with high efficiency of power conversion (∼1%) come into play when the radial size of the beam–plasma system becomes comparable with the wavelength of the emitted waves.


Sign in / Sign up

Export Citation Format

Share Document