Concurrent particle-in-cell plasma simulation on a multi-transputer parallel computer

International Nuclear Information System (INIS)

Khare, A.N.; Jethra, A.; Patel, Kartik

1992-01-01

This report describes the parallelization of a Particle-in-Cell (PIC) plasma simulation code on a multi-transputer parallel computer. The algorithm used in the parallelization of the PIC method is described. The decomposition schemes related to the distribution of the particles among the processors are discussed. The implementation of the algorithm on a transputer network connected as a torus is presented. The solutions of the problems related to global communication of data are presented in the form of a set of generalized communication functions. The performance of the program as a function of data size and the number of transputers show that the implementation is scalable and represents an effective way of achieving high performance at acceptable cost. (author). 11 refs., 4 figs., 2 tabs., appendices

Study of effect of grain size on dust charging in an RF plasma using three-dimensional PIC-MCC simulations

International Nuclear Information System (INIS)

Ikkurthi, V. R.; Melzer, A.; Matyash, K.; Schneider, R.

2008-01-01

A 3-dimensional Particle-Particle Particle-Mesh (P 3 M) code is applied to study the charging process of micrometer size dust grains confined in a capacitive RF discharge. In our model, particles (electrons and ions) are treated kinetically (Particle-in-Cell with Monte Carlo Collisions (PIC-MCC)). In order to accurately resolve the plasma particles' motion close to the dust grain, the PIC technique is supplemented with Molecular Dynamics (MD), employing an an analytic electrostatic potential for the interaction with the dust grain. This allows to self-consistently resolve the dust grain charging due to absorption of plasma electrons and ions. The charging of dust grains confined above lower electrode in a capacitive RF discharge and its dependence on the size and position of the dust is investigated. The results have been compared with laboratory measurements

Implementation of a 3D plasma particle-in-cell code on a MIMD parallel computer

International Nuclear Information System (INIS)

Liewer, P.C.; Lyster, P.; Wang, J.

1993-01-01

A three-dimensional plasma particle-in-cell (PIC) code has been implemented on the Intel Delta MIMD parallel supercomputer using the General Concurrent PIC algorithm. The GCPIC algorithm uses a domain decomposition to divide the computation among the processors: A processor is assigned a subdomain and all the particles in it. Particles must be exchanged between processors as they move. Results are presented comparing the efficiency for 1-, 2- and 3-dimensional partitions of the three dimensional domain. This algorithm has been found to be very efficient even when a large fraction (e.g. 30%) of the particles must be exchanged at every time step. On the 512-node Intel Delta, up to 125 million particles have been pushed with an electrostatic push time of under 500 nsec/particle/time step

Particle-in-Cell Modeling of Magnetized Argon Plasma Flow Through Small Mechanical Apertures

Energy Technology Data Exchange (ETDEWEB)

Adam B. Sefkow and Samuel A. Cohen

2009-04-09

Motivated by observations of supersonic argon-ion flow generated by linear helicon-heated plasma devices, a three-dimensional particle-in-cell (PIC) code is used to study whether stationary electrostatic layers form near mechanical apertures intersecting the flow of magnetized plasma. By self-consistently evaluating the temporal evolution of the plasma in the vicinity of the aperture, the PIC simulations characterize the roles of the imposed aperture and applied magnetic field on ion acceleration. The PIC model includes ionization of a background neutral-argon population by thermal and superthermal electrons, the latter found upstream of the aperture. Near the aperture, a transition from a collisional to a collisionless regime occurs. Perturbations of density and potential, with mm wavelengths and consistent with ion acoustic waves, propagate axially. An ion acceleration region of length ~ 200-300 λD,e forms at the location of the aperture and is found to be an electrostatic double layer, with axially-separated regions of net positive and negative charge. Reducing the aperture diameter or increasing its length increases the double layer strength.

Particle-in-Cell Modeling of Magnetized Argon Plasma Flow Through Small Mechanical Apertures

International Nuclear Information System (INIS)

Sefkow, Adam B.; Cohen, Samuel A.

2009-01-01

Motivated by observations of supersonic argon-ion flow generated by linear helicon-heated plasma devices, a three-dimensional particle-in-cell (PIC) code is used to study whether stationary electrostatic layers form near mechanical apertures intersecting the flow of magnetized plasma. By self-consistently evaluating the temporal evolution of the plasma in the vicinity of the aperture, the PIC simulations characterize the roles of the imposed aperture and applied magnetic field on ion acceleration. The PIC model includes ionization of a background neutral-argon population by thermal and superthermal electrons, the latter found upstream of the aperture. Near the aperture, a transition from a collisional to a collisionless regime occurs. Perturbations of density and potential, with mm wavelengths and consistent with ion acoustic waves, propagate axially. An ion acceleration region of length ∼ 200-300 λ D,e forms at the location of the aperture and is found to be an electrostatic double layer, with axially-separated regions of net positive and negative charge. Reducing the aperture diameter or increasing its length increases the double layer strength

Towards the optimization of a gyrokinetic Particle-In-Cell (PIC) code on large-scale hybrid architectures

International Nuclear Information System (INIS)

Ohana, N; Lanti, E; Tran, T M; Brunner, S; Hariri, F; Villard, L; Jocksch, A; Gheller, C

2016-01-01

With the aim of enabling state-of-the-art gyrokinetic PIC codes to benefit from the performance of recent multithreaded devices, we developed an application from a platform called the “PIC-engine” [1, 2, 3] embedding simplified basic features of the PIC method. The application solves the gyrokinetic equations in a sheared plasma slab using B-spline finite elements up to fourth order to represent the self-consistent electrostatic field. Preliminary studies of the so-called Particle-In-Fourier (PIF) approach, which uses Fourier modes as basis functions in the periodic dimensions of the system instead of the real-space grid, show that this method can be faster than PIC for simulations with a small number of Fourier modes. Similarly to the PIC-engine, multiple levels of parallelism have been implemented using MPI+OpenMP [2] and MPI+OpenACC [1], the latter exploiting the computational power of GPUs without requiring complete code rewriting. It is shown that sorting particles [3] can lead to performance improvement by increasing data locality and vectorizing grid memory access. Weak scalability tests have been successfully run on the GPU-equipped Cray XC30 Piz Daint (at CSCS) up to 4,096 nodes. The reduced time-to-solution will enable more realistic and thus more computationally intensive simulations of turbulent transport in magnetic fusion devices. (paper)

Deploying electromagnetic particle-in-cell (EM-PIC) codes on Xeon Phi accelerators boards

Science.gov (United States)

Fonseca, Ricardo

2014-10-01

The complexity of the phenomena involved in several relevant plasma physics scenarios, where highly nonlinear and kinetic processes dominate, makes purely theoretical descriptions impossible. Further understanding of these scenarios requires detailed numerical modeling, but fully relativistic particle-in-cell codes such as OSIRIS are computationally intensive. The quest towards Exaflop computer systems has lead to the development of HPC systems based on add-on accelerator cards, such as GPGPUs and more recently the Xeon Phi accelerators that power the current number 1 system in the world. These cards, also referred to as Intel Many Integrated Core Architecture (MIC) offer peak theoretical performances of >1 TFlop/s for general purpose calculations in a single board, and are receiving significant attention as an attractive alternative to CPUs for plasma modeling. In this work we report on our efforts towards the deployment of an EM-PIC code on a Xeon Phi architecture system. We will focus on the parallelization and vectorization strategies followed, and present a detailed performance evaluation of code performance in comparison with the CPU code.

The effects of particle recycling on the divertor plasma: A particle-in-cell with Monte Carlo collision simulation

Science.gov (United States)

Chang, Mingyu; Sang, Chaofeng; Sun, Zhenyue; Hu, Wanpeng; Wang, Dezhen

2018-05-01

A Particle-In-Cell (PIC) with Monte Carlo Collision (MCC) model is applied to study the effects of particle recycling on divertor plasma in the present work. The simulation domain is the scrape-off layer of the tokamak in one-dimension along the magnetic field line. At the divertor plate, the reflected deuterium atoms (D) and thermally released deuterium molecules (D2) are considered. The collisions between the plasma particles (e and D+) and recycled neutral particles (D and D2) are described by the MCC method. It is found that the recycled neutral particles have a great impact on divertor plasma. The effects of different collisions on the plasma are simulated and discussed. Moreover, the impacts of target materials on the plasma are simulated by comparing the divertor with Carbon (C) and Tungsten (W) targets. The simulation results show that the energy and momentum losses of the C target are larger than those of the W target in the divertor region even without considering the impurity particles, whereas the W target has a more remarkable influence on the core plasma.

Computer-assisted Particle-in-Cell code development

International Nuclear Information System (INIS)

Kawata, S.; Boonmee, C.; Teramoto, T.; Drska, L.; Limpouch, J.; Liska, R.; Sinor, M.

1997-12-01

This report presents a new approach for an electromagnetic Particle-in-Cell (PIC) code development by a computer: in general PIC codes have a common structure, and consist of a particle pusher, a field solver, charge and current density collections, and a field interpolation. Because of the common feature, the main part of the PIC code can be mechanically developed on a computer. In this report we use the packages FIDE and GENTRAN of the REDUCE computer algebra system for discretizations of field equations and a particle equation, and for an automatic generation of Fortran codes. The approach proposed is successfully applied to the development of 1.5-dimensional PIC code. By using the generated PIC code the Weibel instability in a plasma is simulated. The obtained growth rate agrees well with the theoretical value. (author)

A particle-in-cell method for modeling small angle Coulomb collisions in plasmas

International Nuclear Information System (INIS)

Parker, S.E.

1989-01-01

We propose a computational method to self-consistently model small angle collisional effects. This method may be added to standard Particle-In-Cell (PIC) plasma simulations to include collisions, or as an alternative to solving the Fokker-Planck (FP) equation using finite difference methods. The distribution function is represented by a large number of particles. The particle velocities change due to the drag force, and the diffusion in velocity is represented by a random process. This is similar to previous Monte-Carlo methods except we calculate the drag force and diffusion tensor self- consistently. The particles are weighted to a grid in velocity space and associated ''Poisson equations'' are solved for the Rosenbluth potentials. The motivation is to avoid the very time consuming method of Coulomb scattering pair by pair. First the approximation for small angle Coulomb collisions is discussed. Next, the FP-PIC collision method is outlined. Then we show a test of the particle advance modeling an electron beam scattering off a fixed ion background. 4 refs

Comparison of several algorithms of the electric force calculation in particle plasma models

International Nuclear Information System (INIS)

Lachnitt, J; Hrach, R

2014-01-01

This work is devoted to plasma modelling using the technique of molecular dynamics. The crucial problem of most such models is the efficient calculation of electric force. This is usually solved by using the particle-in-cell (PIC) algorithm. However, PIC is an approximative algorithm as it underestimates the short-range interactions of charged particles. We propose a hybrid algorithm which adds these interactions to PIC. Then we include this algorithm in a set of algorithms which we test against each other in a two-dimensional collisionless magnetized plasma model. Besides our hybrid algorithm, this set includes two variants of pure PIC and the direct application of Coulomb's law. We compare particle forces, particle trajectories, total energy conservation and the speed of the algorithms. We find out that the hybrid algorithm can be a good replacement of direct Coulomb's law application (quite accurate and much faster). It is however probably unnecessary to use it in practical 2D models.

Status and future plans for open source QuickPIC

Science.gov (United States)

An, Weiming; Decyk, Viktor; Mori, Warren

2017-10-01

QuickPIC is a three dimensional (3D) quasi-static particle-in-cell (PIC) code developed based on the UPIC framework. It can be used for efficiently modeling plasma based accelerator (PBA) problems. With quasi-static approximation, QuickPIC can use different time scales for calculating the beam (or laser) evolution and the plasma response, and a 3D plasma wake field can be simulated using a two-dimensional (2D) PIC code where the time variable is ξ = ct - z and z is the beam propagation direction. QuickPIC can be thousand times faster than the normal PIC code when simulating the PBA. It uses an MPI/OpenMP hybrid parallel algorithm, which can be run on either a laptop or the largest supercomputer. The open source QuickPIC is an object-oriented program with high level classes written in Fortran 2003. It can be found at https://github.com/UCLA-Plasma-Simulation-Group/QuickPIC-OpenSource.git

Recent advances in the modeling of plasmas with the Particle-In-Cell methods

Science.gov (United States)

Vay, Jean-Luc; Lehe, Remi; Vincenti, Henri; Godfrey, Brendan; Lee, Patrick; Haber, Irv

2015-11-01

The Particle-In-Cell (PIC) approach is the method of choice for self-consistent simulations of plasmas from first principles. The fundamentals of the PIC method were established decades ago but improvements or variations are continuously being proposed. We report on several recent advances in PIC related algorithms, including: (a) detailed analysis of the numerical Cherenkov instability and its remediation, (b) analytic pseudo-spectral electromagnetic solvers in Cartesian and cylindrical (with azimuthal modes decomposition) geometries, (c) arbitrary-order finite-difference and generalized pseudo-spectral Maxwell solvers, (d) novel analysis of Maxwell's solvers' stencil variation and truncation, in application to domain decomposition strategies and implementation of Perfectly Matched Layers in high-order and pseudo-spectral solvers. Work supported by US-DOE Contracts DE-AC02-05CH11231 and the US-DOE SciDAC program ComPASS. Used resources of NERSC, supported by US-DOE Contract DE-AC02-05CH11231.

Particle acceleration in near critical density plasma

International Nuclear Information System (INIS)

Gu, Y.J.; Kong, Q.; Kawata, S.; Izumiyama, T.; Nagashima, T.

2013-01-01

Charged particle acceleration schemes driven by ultra intense laser and near critical density plasma interactions are presented. They include electron acceleration in a plasma channel, ion acceleration by the Coulomb explosion and high energy electron beam driven ion acceleration. It is found that under the near critical density plasma both ions and electrons are accelerated with a high acceleration gradient. The electron beam containing a large charge quantity is accelerated well with 23 GeV/cm. The collimated ion bunch reaches 1 GeV. The investigations and discussions are based on 2.5D PIC (particle-in-cell) simulations. (author)

Apar-T: code, validation, and physical interpretation of particle-in-cell results

Science.gov (United States)

Melzani, Mickaël; Winisdoerffer, Christophe; Walder, Rolf; Folini, Doris; Favre, Jean M.; Krastanov, Stefan; Messmer, Peter

2013-10-01

We present the parallel particle-in-cell (PIC) code Apar-T and, more importantly, address the fundamental question of the relations between the PIC model, the Vlasov-Maxwell theory, and real plasmas. First, we present four validation tests: spectra from simulations of thermal plasmas, linear growth rates of the relativistic tearing instability and of the filamentation instability, and nonlinear filamentation merging phase. For the filamentation instability we show that the effective growth rates measured on the total energy can differ by more than 50% from the linear cold predictions and from the fastest modes of the simulation. We link these discrepancies to the superparticle number per cell and to the level of field fluctuations. Second, we detail a new method for initial loading of Maxwell-Jüttner particle distributions with relativistic bulk velocity and relativistic temperature, and explain why the traditional method with individual particle boosting fails. The formulation of the relativistic Harris equilibrium is generalized to arbitrary temperature and mass ratios. Both are required for the tearing instability setup. Third, we turn to the key point of this paper and scrutinize the question of what description of (weakly coupled) physical plasmas is obtained by PIC models. These models rely on two building blocks: coarse-graining, i.e., grouping of the order of p ~ 1010 real particles into a single computer superparticle, and field storage on a grid with its subsequent finite superparticle size. We introduce the notion of coarse-graining dependent quantities, i.e., quantities depending on p. They derive from the PIC plasma parameter ΛPIC, which we show to behave as ΛPIC ∝ 1/p. We explore two important implications. One is that PIC collision- and fluctuation-induced thermalization times are expected to scale with the number of superparticles per grid cell, and thus to be a factor p ~ 1010 smaller than in real plasmas, a fact that we confirm with

High-Fidelity RF Gun Simulations with the Parallel 3D Finite Element Particle-In-Cell Code Pic3P

Energy Technology Data Exchange (ETDEWEB)

Candel, A; Kabel, A.; Lee, L.; Li, Z.; Limborg, C.; Ng, C.; Schussman, G.; Ko, K.; /SLAC

2009-06-19

SLAC's Advanced Computations Department (ACD) has developed the first parallel Finite Element 3D Particle-In-Cell (PIC) code, Pic3P, for simulations of RF guns and other space-charge dominated beam-cavity interactions. Pic3P solves the complete set of Maxwell-Lorentz equations and thus includes space charge, retardation and wakefield effects from first principles. Pic3P uses higher-order Finite Elementmethods on unstructured conformal meshes. A novel scheme for causal adaptive refinement and dynamic load balancing enable unprecedented simulation accuracy, aiding the design and operation of the next generation of accelerator facilities. Application to the Linac Coherent Light Source (LCLS) RF gun is presented.

An EDDY/particle-in-cell simulation of erosion of plasma facing walls bombarded by a collisional plasma

International Nuclear Information System (INIS)

Inai, Kensuke; Ohya, Kaoru

2011-01-01

To investigate the erosion of a plasma-facing wall intersecting an oblique magnetic field, we performed a kinetic particle-in-cell (PIC) simulation of magnetized plasma, in which collision processes between charged and neutral particles were taken into account. Sheath formation and local physical quantities, such as the incident angle and energy distributions of plasma ions at the wall, were examined at a plasma density of 10 18 m -3 , a temperature of 10 eV, and a magnetic field strength of 5 T. The erosion rate of a carbon wall was calculated using the ion-solid interaction code EDDY. At a high neutral density (>10 20 m -3 ), the impact energy of the ions dropped below the threshold for physical sputtering, so that the sputtering yield was drastically decreased and wall erosion was strongly suppressed. Sputter erosion was also suppressed when the angle of the magnetic field with respect to the surface normal was sufficiently large. (author)

Studies of Particle Wake Potentials in Plasmas

Science.gov (United States)

Ellis, Ian; Graziani, Frank; Glosli, James; Strozzi, David; Surh, Michael; Richards, David; Decyk, Viktor; Mori, Warren

2011-10-01

Fast Ignition studies require a detailed understanding of electron scattering, stopping, and energy deposition in plasmas with variable values for the number of particles within a Debye sphere. Presently there is disagreement in the literature concerning the proper description of these processes. Developing and validating proper descriptions requires studying the processes using first-principle electrostatic simulations and possibly including magnetic fields. We are using the particle-particle particle-mesh (PPPM) code ddcMD and the particle-in-cell (PIC) code BEPS to perform these simulations. As a starting point in our study, we examine the wake of a particle passing through a plasma in 3D electrostatic simulations performed with ddcMD and with BEPS using various cell sizes. In this poster, we compare the wakes we observe in these simulations with each other and predictions from Vlasov theory. Prepared by LLNL under Contract DE-AC52-07NA27344 and by UCLA under Grant DE-FG52-09NA29552.

Particle-in-cell simulations of fast magnetic field penetration into plasmas due to the Hall electric field

International Nuclear Information System (INIS)

Swanekamp, S.B.; Grossmann, J.M.; Fruchtman, A.; Oliver, B.V.; Ottinger, P.F.

1996-01-01

Particle-in-cell (PIC) simulations are used to study the penetration of magnetic field into plasmas in the electron-magnetohydrodynamic (EMHD) regime. These simulations represent the first definitive verification of EMHD with a PIC code. When ions are immobile, the PIC results reproduce many aspects of fluid treatments of the problem. However, the PIC results show a speed of penetration that is between 10% and 50% slower than predicted by one-dimensional fluid treatments. In addition, the PIC simulations show the formation of vortices in the electron flow behind the EMHD shock front. The size of these vortices is on the order of the collisionless electron skin depth and is closely coupled to the effects of electron inertia. An energy analysis shows that one-half the energy entering the plasma is stored as magnetic field energy while the other half is shared between internal plasma energy (thermal motion and electron vortices) and electron kinetic energy loss from the volume to the boundaries. The amount of internal plasma energy saturates after an initial transient phase so that late in time the rate that magnetic energy increases in the plasma is the same as the rate at which kinetic energy flows out through the boundaries. When ions are mobile it is observed that axial magnetic field penetration is followed by localized thinning in the ion density. The density thinning is produced by the large electrostatic fields that exist inside the electron vortices which act to reduce the space-charge imbalance necessary to support the vortices. This mechanism may play a role during the opening process of a plasma opening switch. copyright 1996 American Institute of Physics

Two-way coupling of magnetohydrodynamic simulations with embedded particle-in-cell simulations

Science.gov (United States)

Makwana, K. D.; Keppens, R.; Lapenta, G.

2017-12-01

We describe a method for coupling an embedded domain in a magnetohydrodynamic (MHD) simulation with a particle-in-cell (PIC) method. In this two-way coupling we follow the work of Daldorff et al. (2014) [19] in which the PIC domain receives its initial and boundary conditions from MHD variables (MHD to PIC coupling) while the MHD simulation is updated based on the PIC variables (PIC to MHD coupling). This method can be useful for simulating large plasma systems, where kinetic effects captured by particle-in-cell simulations are localized but affect global dynamics. We describe the numerical implementation of this coupling, its time-stepping algorithm, and its parallelization strategy, emphasizing the novel aspects of it. We test the stability and energy/momentum conservation of this method by simulating a steady-state plasma. We test the dynamics of this coupling by propagating plasma waves through the embedded PIC domain. Coupling with MHD shows satisfactory results for the fast magnetosonic wave, but significant distortion for the circularly polarized Alfvén wave. Coupling with Hall-MHD shows excellent coupling for the whistler wave. We also apply this methodology to simulate a Geospace Environmental Modeling (GEM) challenge type of reconnection with the diffusion region simulated by PIC coupled to larger scales with MHD and Hall-MHD. In both these cases we see the expected signatures of kinetic reconnection in the PIC domain, implying that this method can be used for reconnection studies.

SHARP: A Spatially Higher-order, Relativistic Particle-in-cell Code

Energy Technology Data Exchange (ETDEWEB)

Shalaby, Mohamad; Broderick, Avery E. [Department of Physics and Astronomy, University of Waterloo, 200 University Avenue West, Waterloo, ON, N2L 3G1 (Canada); Chang, Philip [Department of Physics, University of Wisconsin-Milwaukee, 1900 E. Kenwood Boulevard, Milwaukee, WI 53211 (United States); Pfrommer, Christoph [Leibniz-Institut für Astrophysik Potsdam (AIP), An der Sternwarte 16, D-14482 Potsdam (Germany); Lamberts, Astrid [Theoretical Astrophysics, California Institute of Technology, Pasadena, CA 91125 (United States); Puchwein, Ewald, E-mail: mshalaby@live.ca [Institute of Astronomy and Kavli Institute for Cosmology, University of Cambridge, Madingley Road, Cambridge, CB3 0HA (United Kingdom)

2017-05-20

Numerical heating in particle-in-cell (PIC) codes currently precludes the accurate simulation of cold, relativistic plasma over long periods, severely limiting their applications in astrophysical environments. We present a spatially higher-order accurate relativistic PIC algorithm in one spatial dimension, which conserves charge and momentum exactly. We utilize the smoothness implied by the usage of higher-order interpolation functions to achieve a spatially higher-order accurate algorithm (up to the fifth order). We validate our algorithm against several test problems—thermal stability of stationary plasma, stability of linear plasma waves, and two-stream instability in the relativistic and non-relativistic regimes. Comparing our simulations to exact solutions of the dispersion relations, we demonstrate that SHARP can quantitatively reproduce important kinetic features of the linear regime. Our simulations have a superior ability to control energy non-conservation and avoid numerical heating in comparison to common second-order schemes. We provide a natural definition for convergence of a general PIC algorithm: the complement of physical modes captured by the simulation, i.e., those that lie above the Poisson noise, must grow commensurately with the resolution. This implies that it is necessary to simultaneously increase the number of particles per cell and decrease the cell size. We demonstrate that traditional ways for testing for convergence fail, leading to plateauing of the energy error. This new PIC code enables us to faithfully study the long-term evolution of plasma problems that require absolute control of the energy and momentum conservation.

An arbitrary curvilinear-coordinate method for particle-in-cell modeling

International Nuclear Information System (INIS)

Fichtl, C A; Finn, J M; Cartwright, K L

2012-01-01

A new approach to kinetic simulation of plasmas in complex geometries, based on the particle-in-cell (PIC) simulation method, is explored. In the two-dimensional (2D) electrostatic version of our method, called the arbitrary curvilinear-coordinate PIC method, all essential PIC operations are carried out in 2D on a uniform grid on the unit square logical domain, and mapped to a nonuniform boundary-fitted grid on the physical domain. As the resulting logical grid equations of motion are not separable, we have developed an extension of the semi-implicit modified leapfrog integration technique to preserve the symplectic nature of the logical grid particle mover. A generalized, curvilinear-coordinate formulation of Poisson's equations to solve for the electrostatic fields on the uniform logical grid is also developed. By our formulation, we compute the plasma charge density on the logical grid based on the particles' positions on the logical domain. That is, the plasma particles are weighted to the uniform logical grid and the self-consistent mean electrostatic fields obtained from the solution of the logical grid Poisson equation are interpolated to the particle positions on the logical grid. This process eliminates the complexity associated with the weighting and interpolation processes on the nonuniform physical grid and allows us to run the PIC method on arbitrary boundary-fitted meshes. (paper)

Study on Characteristics of Constricted DC Plasma Using Particle-In-Cell Simulator

International Nuclear Information System (INIS)

Jo, Jong Gap; Park, Yeong Shin; Hwang, Yong Seok

2010-01-01

In dc glow discharge, when anode size is smaller than cathode, very small and bright plasma ball occurs in front of anode. This plasma is called constricted dc plasma and characterized by a high plasma density in positive glow, so called plasma ball, compared to the conventional dc plasma. For the reason, this plasma is utilized to ion or electron beam sources since the beam currents are enhanced by the dense anode glow. However, correlations between characteristics of the plasma (plasma density, electron temperature and space potential) and discharge conditions (anode size, discharge voltage, discharge current, pressure) have been a little investigated definitely clear in previous study because of the trouble of a diagnosis. The plasma ball which is the most essential part of the constricted plasma is too small to diagnose precisely without disturbing plasma. Therefore, we tried to analyze the constricted plasma through computer simulation with Particle-In-Cell (PIC) code. In this study, simulation result of constricted dc plasma as well as conventional dc glow discharge will be addressed and compared with each others

Low frequency, electrodynamic simulation of kinetic plasmas with the DArwin Direct Implicit Particle-In-Cell (DADIPIC) method

International Nuclear Information System (INIS)

Gibbons, M.R.

1995-06-01

This dissertation describes a new algorithm for simulating low frequency, kinetic phenomena in plasmas. DArwin Direct Implicit Particle-in-Cell (DADIPIC), as its name implies, is a combination of the Darwin and direct implicit methods. One of the difficulties in simulating plasmas lies in the enormous disparity between the fundamental scale lengths of a plasma and the scale lengths of the phenomena of interest. The objective is to create models which can ignore the fundamental constraints without eliminating relevant plasma properties. Over the past twenty years several PIC methods have been investigated for overcoming the constraints on explicit electrodynamic PIC. These models eliminate selected high frequency plasma phenomena while retaining kinetic phenomena at low frequency. This dissertation shows that the combination of Darwin and Direct Implicit allows them to operate better than they have been shown to operate in the past. Through the Darwin method the hyperbolic Maxwell's equations are reformulated into a set of elliptic equations. Propagating light waves do not exist in the formulation so the Courant constraint on the time step is eliminated. The Direct Implicit method is applied only to the electrostatic field with the result that electrostatic plasma oscillations do not have to be resolved for stability. With the elimination of these constraints spatial and temporal discretization can be much larger than that possible with explicit, electrodynamic PIC. The code functions in a two dimensional Cartesian region and has been implemented with all components of the particle velocities, the E-field, and the B-field. Internal structures, conductors or dielectrics, may be placed in the simulation region, can be set at desired potentials, and driven with specified currents

Global Particle-in-Cell Simulations of Mercury's Magnetosphere

Science.gov (United States)

Schriver, D.; Travnicek, P. M.; Lapenta, G.; Amaya, J.; Gonzalez, D.; Richard, R. L.; Berchem, J.; Hellinger, P.

2017-12-01

Spacecraft observations of Mercury's magnetosphere have shown that kinetic ion and electron particle effects play a major role in the transport, acceleration, and loss of plasma within the magnetospheric system. Kinetic processes include reconnection, the breakdown of particle adiabaticity and wave-particle interactions. Because of the vast range in spatial scales involved in magnetospheric dynamics, from local electron Debye length scales ( meters) to solar wind/planetary magnetic scale lengths (tens to hundreds of planetary radii), fully self-consistent kinetic simulations of a global planetary magnetosphere remain challenging. Most global simulations of Earth's and other planet's magnetosphere are carried out using MHD, enhanced MHD (e.g., Hall MHD), hybrid, or a combination of MHD and particle in cell (PIC) simulations. Here, 3D kinetic self-consistent hybrid (ion particle, electron fluid) and full PIC (ion and electron particle) simulations of the solar wind interaction with Mercury's magnetosphere are carried out. Using the implicit PIC and hybrid simulations, Mercury's relatively small, but highly kinetic magnetosphere will be examined to determine how the self-consistent inclusion of electrons affects magnetic reconnection, particle transport and acceleration of plasma at Mercury. Also the spatial and energy profiles of precipitating magnetospheric ions and electrons onto Mercury's surface, which can strongly affect the regolith in terms of space weathering and particle outflow, will be examined with the PIC and hybrid codes. MESSENGER spacecraft observations are used both to initiate and validate the global kinetic simulations to achieve a deeper understanding of the role kinetic physics play in magnetospheric dynamics.

Acceleration of PIC simulation with GPU

International Nuclear Information System (INIS)

Suzuki, Junya; Shimazu, Hironori; Fukazawa, Keiichiro; Den, Mitsue

2011-01-01

Particle-in-cell (PIC) is a simulation technique for plasma physics. The large number of particles in high-resolution plasma simulation increases the volume computation required, making it vital to increase computation speed. In this study, we attempt to accelerate computation speed on graphics processing units (GPUs) using KEMPO, a PIC simulation code package. We perform two tests for benchmarking, with small and large grid sizes. In these tests, we run KEMPO1 code using a CPU only, both a CPU and a GPU, and a GPU only. The results showed that performance using only a GPU was twice that of using a CPU alone. While, execution time for using both a CPU and GPU is comparable to the tests with a CPU alone, because of the significant bottleneck in communication between the CPU and GPU. (author)

A 3D gyrokinetic particle-in-cell simulation of fusion plasma microturbulence on parallel computers

Science.gov (United States)

Williams, T. J.

1992-12-01

One of the grand challenge problems now supported by HPCC is the Numerical Tokamak Project. A goal of this project is the study of low-frequency micro-instabilities in tokamak plasmas, which are believed to cause energy loss via turbulent thermal transport across the magnetic field lines. An important tool in this study is gyrokinetic particle-in-cell (PIC) simulation. Gyrokinetic, as opposed to fully-kinetic, methods are particularly well suited to the task because they are optimized to study the frequency and wavelength domain of the microinstabilities. Furthermore, many researchers now employ low-noise delta(f) methods to greatly reduce statistical noise by modelling only the perturbation of the gyrokinetic distribution function from a fixed background, not the entire distribution function. In spite of the increased efficiency of these improved algorithms over conventional PIC algorithms, gyrokinetic PIC simulations of tokamak micro-turbulence are still highly demanding of computer power--even fully-vectorized codes on vector supercomputers. For this reason, we have worked for several years to redevelop these codes on massively parallel computers. We have developed 3D gyrokinetic PIC simulation codes for SIMD and MIMD parallel processors, using control-parallel, data-parallel, and domain-decomposition message-passing (DDMP) programming paradigms. This poster summarizes our earlier work on codes for the Connection Machine and BBN TC2000 and our development of a generic DDMP code for distributed-memory parallel machines. We discuss the memory-access issues which are of key importance in writing parallel PIC codes, with special emphasis on issues peculiar to gyrokinetic PIC. We outline the domain decompositions in our new DDMP code and discuss the interplay of different domain decompositions suited for the particle-pushing and field-solution components of the PIC algorithm.

Collisional particle-in-cell modeling for energy transport accompanied by atomic processes in dense plasmas

Energy Technology Data Exchange (ETDEWEB)

Mishra, R.; Beg, F. N. [Center for Energy Research, University of California, San Diego, California 92093 (United States); Leblanc, P.; Sentoku, Y. [Department of Physics, University of Nevada, Reno, Nevada 89557 (United States); Wei, M. S. [General Atomics, San Diego, California 92121 (United States)

2013-07-15

Fully relativistic collisional Particle-in-Cell (PIC) code, PICLS, has been developed to study extreme energy density conditions produced in intense laser-solid interaction. Recent extensions to PICLS, such as the implementation of dynamic ionization, binary collisions in a partially ionized plasma, and radiative losses, enhance the efficacy of simulating intense laser plasma interaction and subsequent energy transport in resistive media. Different ionization models are introduced and benchmarked against each other to check the suitability of the model. The atomic physics models are critical to determine the energy deposition and transport in dense plasmas, especially when they consist of high Z (atomic number) materials. Finally we demonstrate the electron transport simulations to show the importance of target material on fast electron dynamics.

Plume expansion of a laser-induced plasma studied with the particle-in-cell method

DEFF Research Database (Denmark)

Ellegaard, O.; Nedelea, T.; Schou, Jørgen

2002-01-01

energy as well as electron energy. We have estimated the time constant for energy transfer between the electrons and the ions. The scaling of these processes is given by a single parameter determined by the Debye length obtained from the electron density in the plasma outside the surface. (C) 2002......The initial stage of laser-induced plasma plume expansion from a solid in vacuum and the effect of the Coulomb field have been studied. We have performed a one-dimensional numerical calculation by mapping the charge on a computational grid according to the particle-in-cell (PIC) method of Birdsall...... et al. It is assumed that the particle ablation from a surface with a fixed temperature takes place as a pulse, i.e. within a finite period of time. A number of characteristic quantities for the plasma plume are compared with similar data for expansion of neutrals as well as fluid models: Density...

Gyrokinetic particle-in-cell simulations of plasma microturbulence on advanced computing platforms

International Nuclear Information System (INIS)

Ethier, S; Tang, W M; Lin, Z

2005-01-01

Since its introduction in the early 1980s, the gyrokinetic particle-in-cell (PIC) method has been very successfully applied to the exploration of many important kinetic stability issues in magnetically confined plasmas. Its self-consistent treatment of charged particles and the associated electromagnetic fluctuations makes this method appropriate for studying enhanced transport driven by plasma turbulence. Advances in algorithms and computer hardware have led to the development of a parallel, global, gyrokinetic code in full toroidal geometry, the gyrokinetic toroidal code (GTC), developed at the Princeton Plasma Physics Laboratory. It has proven to be an invaluable tool to study key effects of low-frequency microturbulence in fusion plasmas. As a high-performance computing applications code, its flexible mixed-model parallel algorithm has allowed GTC to scale to over a thousand processors, which is routinely used for simulations. Improvements are continuously being made. As the US ramps up its support for the International Tokamak Experimental Reactor (ITER), the need for understanding the impact of turbulent transport in burning plasma fusion devices is of utmost importance. Accordingly, the GTC code is at the forefront of the set of numerical tools being used to assess and predict the performance of ITER on critical issues such as the efficiency of energy confinement in reactors

Particle Acceleration in Pulsar Wind Nebulae: PIC Modelling

Science.gov (United States)

Sironi, Lorenzo; Cerutti, Benoît

We discuss the role of PIC simulations in unveiling the origin of the emitting particles in PWNe. After describing the basics of the PIC technique, we summarize its implications for the quiescent and the flaring emission of the Crab Nebula, as a prototype of PWNe. A consensus seems to be emerging that, in addition to the standard scenario of particle acceleration via the Fermi process at the termination shock of the pulsar wind, magnetic reconnection in the wind, at the termination shock and in the Nebula plays a major role in powering the multi-wavelength signatures of PWNe.

Merging for Particle-Mesh Complex Particle Kinetic Modeling of the Multiple Plasma Beams

Science.gov (United States)

Lipatov, Alexander S.

2011-01-01

We suggest a merging procedure for the Particle-Mesh Complex Particle Kinetic (PMCPK) method in case of inter-penetrating flow (multiple plasma beams). We examine the standard particle-in-cell (PIC) and the PMCPK methods in the case of particle acceleration by shock surfing for a wide range of the control numerical parameters. The plasma dynamics is described by a hybrid (particle-ion-fluid-electron) model. Note that one may need a mesh if modeling with the computation of an electromagnetic field. Our calculations use specified, time-independent electromagnetic fields for the shock, rather than self-consistently generated fields. While a particle-mesh method is a well-verified approach, the CPK method seems to be a good approach for multiscale modeling that includes multiple regions with various particle/fluid plasma behavior. However, the CPK method is still in need of a verification for studying the basic plasma phenomena: particle heating and acceleration by collisionless shocks, magnetic field reconnection, beam dynamics, etc.

Nonequilibrium Gyrokinetic Fluctuation Theory and Sampling Noise in Gyrokinetic Particle-in-cell Simulations

International Nuclear Information System (INIS)

Krommes, John A.

2007-01-01

The present state of the theory of fluctuations in gyrokinetic (GK) plasmas and especially its application to sampling noise in GK particle-in-cell (PIC) simulations is reviewed. Topics addressed include the Δf method, the fluctuation-dissipation theorem for both classical and GK many-body plasmas, the Klimontovich formalism, sampling noise in PIC simulations, statistical closure for partial differential equations, the theoretical foundations of spectral balance in the presence of arbitrary noise sources, and the derivation of Kadomtsev-type equations from the general formalism

Nonequilibrium Gyrokinetic Fluctuation Theory and Sampling Noise in Gyrokinetic Particle-in-cell Simulations

Energy Technology Data Exchange (ETDEWEB)

John A. Krommes

2007-10-09

The present state of the theory of fluctuations in gyrokinetic GK plasmas and especially its application to sampling noise in GK particle-in-cell PIC simulations is reviewed. Topics addressed include the Δf method, the fluctuation-dissipation theorem for both classical and GK many-body plasmas, the Klimontovich formalism, sampling noise in PIC simulations, statistical closure for partial differential equations, the theoretical foundations of spectral balance in the presence of arbitrary noise sources, and the derivation of Kadomtsev-type equations from the general formalism.

Molecular Dynamic Studies of Particle Wake Potentials in Plasmas

Science.gov (United States)

Ellis, Ian; Graziani, Frank; Glosli, James; Strozzi, David; Surh, Michael; Richards, David; Decyk, Viktor; Mori, Warren

2010-11-01

Fast Ignition studies require a detailed understanding of electron scattering, stopping, and energy deposition in plasmas with variable values for the number of particles within a Debye sphere. Presently there is disagreement in the literature concerning the proper description of these processes. Developing and validating proper descriptions requires studying the processes using first-principle electrostatic simulations and possibly including magnetic fields. We are using the particle-particle particle-mesh (P^3M) code ddcMD to perform these simulations. As a starting point in our study, we examined the wake of a particle passing through a plasma. In this poster, we compare the wake observed in 3D ddcMD simulations with that predicted by Vlasov theory and those observed in the electrostatic PIC code BEPS where the cell size was reduced to .03λD.

Charge-conserving FEM-PIC schemes on general grids

International Nuclear Information System (INIS)

Campos Pinto, M.; Jund, S.; Salmon, S.; Sonnendruecker, E.

2014-01-01

Particle-In-Cell (PIC) solvers are a major tool for the understanding of the complex behavior of a plasma or a particle beam in many situations. An important issue for electromagnetic PIC solvers, where the fields are computed using Maxwell's equations, is the problem of discrete charge conservation. In this article, we aim at proposing a general mathematical formulation for charge-conserving finite-element Maxwell solvers coupled with particle schemes. In particular, we identify the finite-element continuity equations that must be satisfied by the discrete current sources for several classes of time-domain Vlasov-Maxwell simulations to preserve the Gauss law at each time step, and propose a generic algorithm for computing such consistent sources. Since our results cover a wide range of schemes (namely curl-conforming finite element methods of arbitrary degree, general meshes in two or three dimensions, several classes of time discretization schemes, particles with arbitrary shape factors and piecewise polynomial trajectories of arbitrary degree), we believe that they provide a useful roadmap in the design of high-order charge-conserving FEM-PIC numerical schemes. (authors)

Particle-in-cell Simulations with Kinetic Electrons

International Nuclear Information System (INIS)

Lewandowski, J.L.V.

2004-01-01

A new scheme, based on an exact separation between adiabatic and nonadiabatic electron responses, for particle-in-cell (PIC) simulations of drift-type modes is presented. The (linear and nonlinear) elliptic equations for the scalar fields are solved using a multi-grid solver. The new scheme yields linear growth rates in excellent agreement with theory and it is shown to conserve energy well into the nonlinear regime. It is also demonstrated that simulations with few electrons are reliable and accurate, suggesting that large-scale, PIC simulations with electron dynamics in toroidal geometry (e.g., tokamaks and stellarators plasmas) are within reach of present-day massively parallel supercomputers

Realistic PIC modelling of laser-plasma interaction: a direct implicit method with adjustable damping and high order weight functions

International Nuclear Information System (INIS)

Drouin, M.

2009-11-01

This research thesis proposes a new formulation of the relativistic implicit direct method, based on the weak formulation of the wave equation which is solved by means of a Newton algorithm. The first part of this thesis deals with the properties of the explicit particle-in-cell (PIC) methods: properties and limitations of an explicit PIC code, linear analysis of a numerical plasma, numerical heating phenomenon, interest of a higher order interpolation function, and presentation of two applications in high density relativistic laser-plasma interaction. The second and main part of this report deals with adapting the direct implicit method to laser-plasma interaction: presentation of the state of the art, formulating of the direct implicit method, resolution of the wave equation. The third part concerns various numerical and physical validations of the ELIXIRS code: case of laser wave propagation in vacuum, demonstration of the adjustable damping which is a characteristic of the proposed algorithm, influence of space-time discretization on energy conservation, expansion of a thermal plasma in vacuum, two cases of plasma-beam unsteadiness in relativistic regime, and then a case of the overcritical laser-plasma interaction

The use of electromagnetic particle-in-cell codes in accelerator applications

International Nuclear Information System (INIS)

Eppley, K.

1988-12-01

The techniques developed for the numerical simulation of plasmas have numerous applications relevant to accelerators. The operation of many accelerator components involves transients, interactions between beams and rf fields, and internal plasma oscillations. These effects produce non-linear behavior which can be represented accurately by particle in cell (PIC) simulations. We will give a very brief overview of the algorithms used in PIC Codes. We will examine the range of parameters over which they are useful. We will discuss the factors which determine whether a two or three dimensional simulation is most appropriate. PIC codes have been applied to a wide variety of diverse problems, spanning many of the systems in a linear accelerator. We will present a number of practical examples of the application of these codes to areas such as guns, bunchers, rf sources, beam transport, emittance growth and final focus. 8 refs., 8 figs., 2 tabs

Plasma density enhancement in atmospheric-pressure dielectric-barrier discharges by high-voltage nanosecond pulse in the pulse-on period: a PIC simulation

International Nuclear Information System (INIS)

Sang Chaofeng; Sun Jizhong; Wang Dezhen

2010-01-01

A particle-in-cell (PIC) plus Monte Carlo collision simulation is employed to investigate how a sustainable atmospheric pressure single dielectric-barrier discharge responds to a high-voltage nanosecond pulse (HVNP) further applied to the metal electrode. The results show that the HVNP can significantly increase the plasma density in the pulse-on period. The ion-induced secondary electrons can give rise to avalanche ionization in the positive sheath, which widens the discharge region and enhances the plasma density drastically. However, the plasma density stops increasing as the applied pulse lasts over certain time; therefore, lengthening the pulse duration alone cannot improve the discharge efficiency further. Physical reasons for these phenomena are then discussed.

Plasma density enhancement in atmospheric-pressure dielectric-barrier discharges by high-voltage nanosecond pulse in the pulse-on period: a PIC simulation

Science.gov (United States)

Sang, Chaofeng; Sun, Jizhong; Wang, Dezhen

2010-02-01

A particle-in-cell (PIC) plus Monte Carlo collision simulation is employed to investigate how a sustainable atmospheric pressure single dielectric-barrier discharge responds to a high-voltage nanosecond pulse (HVNP) further applied to the metal electrode. The results show that the HVNP can significantly increase the plasma density in the pulse-on period. The ion-induced secondary electrons can give rise to avalanche ionization in the positive sheath, which widens the discharge region and enhances the plasma density drastically. However, the plasma density stops increasing as the applied pulse lasts over certain time; therefore, lengthening the pulse duration alone cannot improve the discharge efficiency further. Physical reasons for these phenomena are then discussed.

PIConGPU - How to build one of the fastest GPU particle-in-cell codes in the world

Energy Technology Data Exchange (ETDEWEB)

Burau, Heiko; Debus, Alexander; Helm, Anton; Huebl, Axel; Kluge, Thomas; Widera, Rene; Bussmann, Michael; Schramm, Ulrich; Cowan, Thomas [HZDR, Dresden (Germany); Juckeland, Guido; Nagel, Wolfgang [TU Dresden (Germany); ZIH, Dresden (Germany); Schmitt, Felix [NVIDIA (United States)

2013-07-01

We present the algorithmic building blocks of PIConGPU, one of the fastest implementations of the particle-in-cell algortihm on GPU clusters. PIConGPU is a highly-scalable, 3D3V electromagnetic PIC code that is used in laser plasma and astrophysical plasma simulations.

Particle-in-cell plasma simulation codes on the connection machine

International Nuclear Information System (INIS)

Walker, D.W.

1991-01-01

Methods for implementing three-dimensional, electromagnetic, relativistic PIC plasma simulation codes on the Connection Machine (CM-2) are discussed. The gather and scatter phases of the PIC algorithm involve indirect indexing of data, which results in a large amount of communication on the CM-2. Different data decompositions are described that seek to reduce the amount of communication while maintaining good load balance. These methods require the particles to be spatially sorted at the start of each time step, which introduced another form of overhead. The different methods are implemented in CM Fortran on the CM-2 and compared. It was found that the general router is slow in performing the communication in the gather and scatter steps, which precludes an efficient CM Fortran implementation. An alternative method that uses PARIS calls and the NEWS communication network to pipeline data along the axes of the VP set is suggested as a more efficient algorithm

Load management strategy for Particle-In-Cell simulations in high energy particle acceleration

Energy Technology Data Exchange (ETDEWEB)

Beck, A., E-mail: beck@llr.in2p3.fr [Laboratoire Leprince-Ringuet, École polytechnique, CNRS-IN2P3, Palaiseau 91128 (France); Frederiksen, J.T. [Niels Bohr Institute, University of Copenhagen, Blegdamsvej 17, 2100 København Ø (Denmark); Dérouillat, J. [CEA, Maison de La Simulation, 91400 Saclay (France)

2016-09-01

In the wake of the intense effort made for the experimental CILEX project, numerical simulation campaigns have been carried out in order to finalize the design of the facility and to identify optimal laser and plasma parameters. These simulations bring, of course, important insight into the fundamental physics at play. As a by-product, they also characterize the quality of our theoretical and numerical models. In this paper, we compare the results given by different codes and point out algorithmic limitations both in terms of physical accuracy and computational performances. These limitations are illustrated in the context of electron laser wakefield acceleration (LWFA). The main limitation we identify in state-of-the-art Particle-In-Cell (PIC) codes is computational load imbalance. We propose an innovative algorithm to deal with this specific issue as well as milestones towards a modern, accurate high-performance PIC code for high energy particle acceleration.

The Darwin direct implicit particle-in-cell (DADIPIC) method for simulation of low frequency plasma phenomena

International Nuclear Information System (INIS)

Gibbons, M.R.; Hewett, D.W.

1995-01-01

We describe a new algorithm for simulating low frequency, kinetic phenomena in plasma. Darwin direct implicit particle-in-cell (DADIPIC), as its name implies, is a combination of the Darwin and direct implicit methods. Through the Darwin method the hyperbolic Maxwell's equations are reformulated into a set of elliptic equations. Propagating light waves do not exist in the formulation so the Courant constraint on the time step is eliminated. The direct implicit method is applied only to the electrostatic field with the result that electrostatic plasma oscillations do not have to be resolved for stability. With the elimination of these constraints spatial and temporal discretization can be much larger than that possible with explicit, electrodynamic PIC. We discuss the algorithms for pushing the particles and solving the fields in 2D cartesian geometry. We also detail boundary conditions for conductors and dielectrics. Finally, we present two test cases, electron cyclotron waves and collisionless heating in inductively coupled plasmas. For these test cases DADIPIC shows agreement with analytic kinetic theory and good energy conservation characteristics. 33 refs., 7 figs., 2 tabs

Nonlinear PIC simulation in a Penning trap

International Nuclear Information System (INIS)

Lapenta, G.; Delzanno, G.L.; Finn, J. M.

2002-01-01

We study the nonlinear dynamics of a Penning trap plasma, including the effect of the finite length and end curvature of the plasma column. A new cylindrical PIC code, called KANDINSKY, has been implemented by using a new interpolation scheme. The principal idea is to calculate the volume of each cell from a particle volume, in the same manner as it is done for the cell charge. With this new method, the density is conserved along streamlines and artificial sources of compressibility are avoided. The code has been validated with a reference Eulerian fluid code. We compare the dynamics of three different models: a model with compression effects, the standard Euler model and a geophysical fluid dynamics model. The results of our investigation prove that Penning traps can really be used to simulate geophysical fluids

Second-order particle-in-cell (PIC) computational method in the one-dimensional variable Eulerian mesh system

International Nuclear Information System (INIS)

Pyun, J.J.

1981-01-01

As part of an effort to incorporate the variable Eulerian mesh into the second-order PIC computational method, a truncation error analysis was performed to calculate the second-order error terms for the variable Eulerian mesh system. The results that the maximum mesh size increment/decrement is limited to be α(Δr/sub i/) 2 where Δr/sub i/ is a non-dimensional mesh size of the ith cell, and α is a constant of order one. The numerical solutions of Burgers' equation by the second-order PIC method in the variable Eulerian mesh system wer compared with its exact solution. It was found that the second-order accuracy in the PIC method was maintained under the above condition. Additional problems were analyzed using the second-order PIC methods in both variable and uniform Eulerian mesh systems. The results indicate that the second-order PIC method in the variable Eulerian mesh system can provide substantial computational time saving with no loss in accuracy

On the numerical dispersion of electromagnetic particle-in-cell code: Finite grid instability

International Nuclear Information System (INIS)

Meyers, M.D.; Huang, C.-K.; Zeng, Y.; Yi, S.A.; Albright, B.J.

2015-01-01

The Particle-In-Cell (PIC) method is widely used in relativistic particle beam and laser plasma modeling. However, the PIC method exhibits numerical instabilities that can render unphysical simulation results or even destroy the simulation. For electromagnetic relativistic beam and plasma modeling, the most relevant numerical instabilities are the finite grid instability and the numerical Cherenkov instability. We review the numerical dispersion relation of the Electromagnetic PIC model. We rigorously derive the faithful 3-D numerical dispersion relation of the PIC model, for a simple, direct current deposition scheme, which does not conserve electric charge exactly. We then specialize to the Yee FDTD scheme. In particular, we clarify the presence of alias modes in an eigenmode analysis of the PIC model, which combines both discrete and continuous variables. The manner in which the PIC model updates and samples the fields and distribution function, together with the temporal and spatial phase factors from solving Maxwell's equations on the Yee grid with the leapfrog scheme, is explicitly accounted for. Numerical solutions to the electrostatic-like modes in the 1-D dispersion relation for a cold drifting plasma are obtained for parameters of interest. In the succeeding analysis, we investigate how the finite grid instability arises from the interaction of the numerical modes admitted in the system and their aliases. The most significant interaction is due critically to the correct representation of the operators in the dispersion relation. We obtain a simple analytic expression for the peak growth rate due to this interaction, which is then verified by simulation. We demonstrate that our analysis is readily extendable to charge conserving models

On the numerical dispersion of electromagnetic particle-in-cell code: Finite grid instability

Science.gov (United States)

Meyers, M. D.; Huang, C.-K.; Zeng, Y.; Yi, S. A.; Albright, B. J.

2015-09-01

The Particle-In-Cell (PIC) method is widely used in relativistic particle beam and laser plasma modeling. However, the PIC method exhibits numerical instabilities that can render unphysical simulation results or even destroy the simulation. For electromagnetic relativistic beam and plasma modeling, the most relevant numerical instabilities are the finite grid instability and the numerical Cherenkov instability. We review the numerical dispersion relation of the Electromagnetic PIC model. We rigorously derive the faithful 3-D numerical dispersion relation of the PIC model, for a simple, direct current deposition scheme, which does not conserve electric charge exactly. We then specialize to the Yee FDTD scheme. In particular, we clarify the presence of alias modes in an eigenmode analysis of the PIC model, which combines both discrete and continuous variables. The manner in which the PIC model updates and samples the fields and distribution function, together with the temporal and spatial phase factors from solving Maxwell's equations on the Yee grid with the leapfrog scheme, is explicitly accounted for. Numerical solutions to the electrostatic-like modes in the 1-D dispersion relation for a cold drifting plasma are obtained for parameters of interest. In the succeeding analysis, we investigate how the finite grid instability arises from the interaction of the numerical modes admitted in the system and their aliases. The most significant interaction is due critically to the correct representation of the operators in the dispersion relation. We obtain a simple analytic expression for the peak growth rate due to this interaction, which is then verified by simulation. We demonstrate that our analysis is readily extendable to charge conserving models.

LPIC++. A parallel one-dimensional relativistic electromagnetic particle-in-cell code for simulating laser-plasma-interaction

International Nuclear Information System (INIS)

Lichters, R.; Pfund, R.E.W.; Meyer-ter-Vehn, J.

1997-08-01

The code LPIC++ presented here, is based on a one-dimensional, electromagnetic, relativistic PIC code that has originally been developed by one of the authors during a PhD thesis at the Max-Planck-Institut fuer Quantenoptik for kinetic simulations of high harmonic generation from overdense plasma surfaces. The code uses essentially the algorithm of Birdsall and Langdon and Villasenor and Bunemann. It is written in C++ in order to be easily extendable and has been parallelized to be able to grow in power linearly with the size of accessable hardware, e.g. massively parallel machines like Cray T3E. The parallel LPIC++ version uses PVM for communication between processors. PVM is public domain software, can be downloaded from the world wide web. A particular strength of LPIC++ lies in its clear program and data structure, which uses chained lists for the organization of grid cells and enables dynamic adjustment of spatial domain sizes in a very convenient way, and therefore easy balancing of processor loads. Also particles belonging to one cell are linked in a chained list and are immediately accessable from this cell. In addition to this convenient type of data organization in a PIC code, the code shows excellent performance in both its single processor and parallel version. (orig.)

On the elimination of numerical Cerenkov radiation in PIC simulations

International Nuclear Information System (INIS)

Greenwood, Andrew D.; Cartwright, Keith L.; Luginsland, John W.; Baca, Ernest A.

2004-01-01

Particle-in-cell (PIC) simulations are a useful tool in modeling plasma in physical devices. The Yee finite difference time domain (FDTD) method is commonly used in PIC simulations to model the electromagnetic fields. However, in the Yee FDTD method, poorly resolved waves at frequencies near the cut off frequency of the grid travel slower than the physical speed of light. These slowly traveling, poorly resolved waves are not a problem in many simulations because the physics of interest are at much lower frequencies. However, when high energy particles are present, the particles may travel faster than the numerical speed of their own radiation, leading to non-physical, numerical Cerenkov radiation. Due to non-linear interaction between the particles and the fields, the numerical Cerenkov radiation couples into the frequency band of physical interest and corrupts the PIC simulation. There are two methods of mitigating the effects of the numerical Cerenkov radiation. The computational stencil used to approximate the curl operator can be altered to improve the high frequency physics, or a filtering scheme can be introduced to attenuate the waves that cause the numerical Cerenkov radiation. Altering the computational stencil is more physically accurate but is difficult to implement while maintaining charge conservation in the code. Thus, filtering is more commonly used. Two previously published filters by Godfrey and Friedman are analyzed and compared to ideally desired filter properties

Exactly energy conserving semi-implicit particle in cell formulation

International Nuclear Information System (INIS)

Lapenta, Giovanni

2017-01-01

We report a new particle in cell (PIC) method based on the semi-implicit approach. The novelty of the new method is that unlike any of its semi-implicit predecessors at the same time it retains the explicit computational cycle and conserves energy exactly. Recent research has presented fully implicit methods where energy conservation is obtained as part of a non-linear iteration procedure. The new method (referred to as Energy Conserving Semi-Implicit Method, ECSIM), instead, does not require any non-linear iteration and its computational cycle is similar to that of explicit PIC. The properties of the new method are: i) it conserves energy exactly to round-off for any time step or grid spacing; ii) it is unconditionally stable in time, freeing the user from the need to resolve the electron plasma frequency and allowing the user to select any desired time step; iii) it eliminates the constraint of the finite grid instability, allowing the user to select any desired resolution without being forced to resolve the Debye length; iv) the particle mover has a computational complexity identical to that of the explicit PIC, only the field solver has an increased computational cost. The new ECSIM is tested in a number of benchmarks where accuracy and computational performance are tested. - Highlights: • We present a new fully energy conserving semi-implicit particle in cell (PIC) method based on the implicit moment method (IMM). The new method is called Energy Conserving Implicit Moment Method (ECIMM). • The novelty of the new method is that unlike any of its predecessors at the same time it retains the explicit computational cycle and conserves energy exactly. • The new method is unconditionally stable in time, freeing the user from the need to resolve the electron plasma frequency. • The new method eliminates the constraint of the finite grid instability, allowing the user to select any desired resolution without being forced to resolve the Debye length. • These

Exactly energy conserving semi-implicit particle in cell formulation

Energy Technology Data Exchange (ETDEWEB)

Lapenta, Giovanni, E-mail: giovanni.lapenta@kuleuven.be

2017-04-01

We report a new particle in cell (PIC) method based on the semi-implicit approach. The novelty of the new method is that unlike any of its semi-implicit predecessors at the same time it retains the explicit computational cycle and conserves energy exactly. Recent research has presented fully implicit methods where energy conservation is obtained as part of a non-linear iteration procedure. The new method (referred to as Energy Conserving Semi-Implicit Method, ECSIM), instead, does not require any non-linear iteration and its computational cycle is similar to that of explicit PIC. The properties of the new method are: i) it conserves energy exactly to round-off for any time step or grid spacing; ii) it is unconditionally stable in time, freeing the user from the need to resolve the electron plasma frequency and allowing the user to select any desired time step; iii) it eliminates the constraint of the finite grid instability, allowing the user to select any desired resolution without being forced to resolve the Debye length; iv) the particle mover has a computational complexity identical to that of the explicit PIC, only the field solver has an increased computational cost. The new ECSIM is tested in a number of benchmarks where accuracy and computational performance are tested. - Highlights: • We present a new fully energy conserving semi-implicit particle in cell (PIC) method based on the implicit moment method (IMM). The new method is called Energy Conserving Implicit Moment Method (ECIMM). • The novelty of the new method is that unlike any of its predecessors at the same time it retains the explicit computational cycle and conserves energy exactly. • The new method is unconditionally stable in time, freeing the user from the need to resolve the electron plasma frequency. • The new method eliminates the constraint of the finite grid instability, allowing the user to select any desired resolution without being forced to resolve the Debye length. • These

Development of in-situ visualization tool for PIC simulation

International Nuclear Information System (INIS)

Ohno, Nobuaki; Ohtani, Hiroaki

2014-01-01

As the capability of a supercomputer is improved, the sizes of simulation and its output data also become larger and larger. Visualization is usually carried out on a researcher's PC with interactive visualization software after performing the computer simulation. However, the data size is becoming too large to do it currently. A promising answer is in-situ visualization. For this case a simulation code is coupled with the visualization code and visualization is performed with the simulation on the same supercomputer. We developed an in-situ visualization tool for particle-in-cell (PIC) simulation and it is provided as a Fortran's module. We coupled it with a PIC simulation code and tested the coupled code on Plasma Simulator supercomputer, and ensured that it works. (author)

Efficient particle-in-cell simulation of auroral plasma phenomena using a CUDA enabled graphics processing unit

Science.gov (United States)

Sewell, Stephen

This thesis introduces a software framework that effectively utilizes low-cost commercially available Graphic Processing Units (GPUs) to simulate complex scientific plasma phenomena that are modeled using the Particle-In-Cell (PIC) paradigm. The software framework that was developed conforms to the Compute Unified Device Architecture (CUDA), a standard for general purpose graphic processing that was introduced by NVIDIA Corporation. This framework has been verified for correctness and applied to advance the state of understanding of the electromagnetic aspects of the development of the Aurora Borealis and Aurora Australis. For each phase of the PIC methodology, this research has identified one or more methods to exploit the problem's natural parallelism and effectively map it for execution on the graphic processing unit and its host processor. The sources of overhead that can reduce the effectiveness of parallelization for each of these methods have also been identified. One of the novel aspects of this research was the utilization of particle sorting during the grid interpolation phase. The final representation resulted in simulations that executed about 38 times faster than simulations that were run on a single-core general-purpose processing system. The scalability of this framework to larger problem sizes and future generation systems has also been investigated.

Hybrid-PIC Computer Simulation of the Plasma and Erosion Processes in Hall Thrusters

Science.gov (United States)

Hofer, Richard R.; Katz, Ira; Mikellides, Ioannis G.; Gamero-Castano, Manuel

2010-01-01

HPHall software simulates and tracks the time-dependent evolution of the plasma and erosion processes in the discharge chamber and near-field plume of Hall thrusters. HPHall is an axisymmetric solver that employs a hybrid fluid/particle-in-cell (Hybrid-PIC) numerical approach. HPHall, originally developed by MIT in 1998, was upgraded to HPHall-2 by the Polytechnic University of Madrid in 2006. The Jet Propulsion Laboratory has continued the development of HPHall-2 through upgrades to the physical models employed in the code, and the addition of entirely new ones. Primary among these are the inclusion of a three-region electron mobility model that more accurately depicts the cross-field electron transport, and the development of an erosion sub-model that allows for the tracking of the erosion of the discharge chamber wall. The code is being developed to provide NASA science missions with a predictive tool of Hall thruster performance and lifetime that can be used to validate Hall thrusters for missions.

The Multipole Plasma Trap-PIC Modeling Results

Science.gov (United States)

Hicks, Nathaniel; Bowman, Amanda; Godden, Katarina

2017-10-01

A radio-frequency (RF) multipole structure is studied via particle-in-cell computer modeling, to assess the response of quasi-neutral plasma to the imposed RF fields. Several regimes, such as pair plasma, antimatter plasma, and conventional (ion-electron) plasma are considered. In the case of equal charge-to-mass ratio of plasma species, the effects of the multipole field are symmetric between positive and negative particles. In the case of a charge-to-mass disparity, the multipole RF parameters (frequency, voltage, structure size) may be chosen such that the light species (e.g. electrons) is strongly confined, while the heavy species (e.g. positive ions) does not respond to the RF field. In this case, the trapped negative space charge creates a potential well that then traps the positive species. 2D and 3D particle-in-cell simulations of this concept are presented, to assess plasma response and trapping dependences on multipole order, consequences of the formation of an RF plasma sheath, and the effects of an axial magnetic field. The scalings of trapped plasma parameters are explored in each of the mentioned regimes, to guide the design of prospective experiments investigating each. Supported by U.S. NSF/DOE Partnership in Basic Plasma Science and Engineering Grant PHY-1619615.

Observation and particle simulation of vaporized W, Mo, and Be in PISCES-B plasma for vapor-shielding studies

Directory of Open Access Journals (Sweden)

K. Ibano

2017-08-01

Full Text Available Interactions of Tungsten (W, Molybdenum (Mo, and Beryllium (Be vapors with a steady-state plasma were studied by the PISCES-B liner plasma experiments as well as Particle-In-Cell (PIC simulations for the understanding of vapor-shielding phenomena. Effective cooling of the plasma by laser-generated Be vapor was observed in PISCES-B. On the other hand, no apparent cooling was observed for W and Mo vapors. The PIC simulation explains these experimental observations of the difference between low-Z and high-Z vapors. Decrease of electron temperature due to the vapor ejection was observed in case of a simulation of the Be vapor. As for the W vapor, it was found that the plasma cooling is localized only near the wall at a higher electron density plasma (∼1019m−3. On the other hand, the appreciable plasma cooling can be observed in a lower density plasma (∼1018m−3 for the W vapor.

New methods in plasma simulation

International Nuclear Information System (INIS)

Mason, R.J.

1990-01-01

The development of implicit methods of particle-in-cell (PIC) computer simulation in recent years, and their merger with older hybrid methods have created a new arsenal of simulation techniques for the treatment of complex practical problems in plasma physics. The new implicit hybrid codes are aimed at transitional problems that lie somewhere between the long time scale, high density regime associated with MHD modeling, and the short time scale, low density regime appropriate to PIC particle-in-cell techniques. This transitional regime arises in ICF coronal plasmas, in pulsed power plasma switches, in Z-pinches, and in foil implosions. Here, we outline how such a merger of implicit and hybrid methods has been carried out, specifically in the ANTHEM computer code, and demonstrate the utility of implicit hybrid simulation in applications. 25 refs., 5 figs

Numerical heating in Particle-In-Cell simulations with Monte Carlo binary collisions

Science.gov (United States)

Alves, E. Paulo; Mori, Warren; Fiuza, Frederico

2017-10-01

The binary Monte Carlo collision (BMCC) algorithm is a robust and popular method to include Coulomb collision effects in Particle-in-Cell (PIC) simulations of plasmas. While a number of works have focused on extending the validity of the model to different physical regimes of temperature and density, little attention has been given to the fundamental coupling between PIC and BMCC algorithms. Here, we show that the coupling between PIC and BMCC algorithms can give rise to (nonphysical) numerical heating of the system, that can be far greater than that observed when these algorithms operate independently. This deleterious numerical heating effect can significantly impact the evolution of the simulated system particularly for long simulation times. In this work, we describe the source of this numerical heating, and derive scaling laws for the numerical heating rates based on the numerical parameters of PIC-BMCC simulations. We compare our theoretical scalings with PIC-BMCC numerical experiments, and discuss strategies to minimize this parasitic effect. This work is supported by DOE FES under FWP 100237 and 100182.

The Particle-in-Cell and Kinetic Simulation Software Center

Science.gov (United States)

Mori, W. B.; Decyk, V. K.; Tableman, A.; Fonseca, R. A.; Tsung, F. S.; Hu, Q.; Winjum, B. J.; An, W.; Dalichaouch, T. N.; Davidson, A.; Hildebrand, L.; Joglekar, A.; May, J.; Miller, K.; Touati, M.; Xu, X. L.

2017-10-01

The UCLA Particle-in-Cell and Kinetic Simulation Software Center (PICKSC) aims to support an international community of PIC and plasma kinetic software developers, users, and educators; to increase the use of this software for accelerating the rate of scientific discovery; and to be a repository of knowledge and history for PIC. We discuss progress towards making available and documenting illustrative open-source software programs and distinct production programs; developing and comparing different PIC algorithms; coordinating the development of resources for the educational use of kinetic software; and the outcomes of our first sponsored OSIRIS users workshop. We also welcome input and discussion from anyone interested in using or developing kinetic software, in obtaining access to our codes, in collaborating, in sharing their own software, or in commenting on how PICKSC can better serve the DPP community. Supported by NSF under Grant ACI-1339893 and by the UCLA Institute for Digital Research and Education.

A Generalized Weight-Based Particle-In-Cell Simulation Scheme

International Nuclear Information System (INIS)

Lee, W.W.; Jenkins, T.G.; Ethier, S.

2010-01-01

A generalized weight-based particle simulation scheme suitable for simulating magnetized plasmas, where the zeroth-order inhomogeneity is important, is presented. The scheme is an extension of the perturbative simulation schemes developed earlier for particle-in-cell (PIC) simulations. The new scheme is designed to simulate both the perturbed distribution ((delta)f) and the full distribution (full-F) within the same code. The development is based on the concept of multiscale expansion, which separates the scale lengths of the background inhomogeneity from those associated with the perturbed distributions. The potential advantage for such an arrangement is to minimize the particle noise by using (delta)f in the linear stage stage of the simulation, while retaining the flexibility of a full-F capability in the fully nonlinear stage of the development when signals associated with plasma turbulence are at a much higher level than those from the intrinsic particle noise.

Role of positive ions on the surface production of negative ions in a fusion plasma reactor type negative ion source--Insights from a three dimensional particle-in-cell Monte Carlo collisions model

Science.gov (United States)

Fubiani, G.; Boeuf, J. P.

2013-11-01

Results from a 3D self-consistent Particle-In-Cell Monte Carlo Collisions (PIC MCC) model of a high power fusion-type negative ion source are presented for the first time. The model is used to calculate the plasma characteristics of the ITER prototype BATMAN ion source developed in Garching. Special emphasis is put on the production of negative ions on the plasma grid surface. The question of the relative roles of the impact of neutral hydrogen atoms and positive ions on the cesiated grid surface has attracted much attention recently and the 3D PIC MCC model is used to address this question. The results show that the production of negative ions by positive ion impact on the plasma grid is small with respect to the production by atomic hydrogen or deuterium bombardment (less than 10%).

Role of positive ions on the surface production of negative ions in a fusion plasma reactor type negative ion source—Insights from a three dimensional particle-in-cell Monte Carlo collisions model

International Nuclear Information System (INIS)

Fubiani, G.; Boeuf, J. P.

2013-01-01

Results from a 3D self-consistent Particle-In-Cell Monte Carlo Collisions (PIC MCC) model of a high power fusion-type negative ion source are presented for the first time. The model is used to calculate the plasma characteristics of the ITER prototype BATMAN ion source developed in Garching. Special emphasis is put on the production of negative ions on the plasma grid surface. The question of the relative roles of the impact of neutral hydrogen atoms and positive ions on the cesiated grid surface has attracted much attention recently and the 3D PIC MCC model is used to address this question. The results show that the production of negative ions by positive ion impact on the plasma grid is small with respect to the production by atomic hydrogen or deuterium bombardment (less than 10%)

2D arc-PIC code description: methods and documentation

CERN Document Server

Timko, Helga

2011-01-01

Vacuum discharges are one of the main limiting factors for future linear collider designs such as that of the Compact LInear Collider. To optimize machine efficiency, maintaining the highest feasible accelerating gradient below a certain breakdown rate is desirable; understanding breakdowns can therefore help us to achieve this goal. As a part of ongoing theoretical research on vacuum discharges at the Helsinki Institute of Physics, the build-up of plasma can be investigated through the particle-in-cell method. For this purpose, we have developed the 2D Arc-PIC code introduced here. We present an exhaustive description of the 2D Arc-PIC code in two parts. In the first part, we introduce the particle-in-cell method in general and detail the techniques used in the code. In the second part, we provide a documentation and derivation of the key equations occurring in the code. The code is original work of the author, written in 2010, and is therefore under the copyright of the author. The development of the code h...

The relaxation of plasmas with dust particles

International Nuclear Information System (INIS)

Chutov, Yu.I.; Kravchenko, A.Yu.; Schram, P.P.J.M.

1997-01-01

Various parameters of relaxing plasmas with dust particles including the electron and ion energy distributions function are numerically simulated at various parameters of the dust particles using the PIC method and taking into account the dynamics of the dust particle charge without the assumption about the equilibrium of electrons and ions. Coulomb collisions are taken into account in the framework of the method of stochastic differential equations. The relaxation of bounded plasma clouds expanding into a vacuum as well as the relaxation of a uniform plasma, in which dust particles appear at some initial time, are investigated. The obtained results show that the relaxation of plasmas can be accompanied by a deviation of the ion distribution function from equilibrium as well as a change of the mean energy of electrons and ions because of the dependence of the collection of electrons and ions by dust particles on their energy. (author)

A fully kinetic, self-consistent particle simulation model of the collisionless plasma--sheath region

International Nuclear Information System (INIS)

Procassini, R.J.; Birdsall, C.K.; Morse, E.C.

1990-01-01

A fully kinetic particle-in-cell (PIC) model is used to self-consistently determine the steady-state potential profile in a collisionless plasma that contacts a floating, absorbing boundary. To balance the flow of particles to the wall, a distributed source region is used to inject particles into the one-dimensional system. The effect of the particle source distribution function on the source region and collector sheath potential drops, and particle velocity distributions is investigated. The ion source functions proposed by Emmert et al. [Phys. Fluids 23, 803 (1980)] and Bissell and Johnson [Phys. Fluids 30, 779 (1987)] (and various combinations of these) are used for the injection of both ions and electrons. The values of the potential drops obtained from the PIC simulations are compared to those from the theories of Emmert et al., Bissell and Johnson, and Scheuer and Emmert [Phys. Fluids 31, 3645 (1988)], all of which assume that the electron density is related to the plasma potential via the Boltzmann relation. The values of the source region and total potential drop are found to depend on the choice of the electron source function, as well as the ion source function. The question of an infinite electric field at the plasma--sheath interface, which arises in the analyses of Bissell and Johnson and Scheuer and Emmert, is also addressed

PIC simulation of electron acceleration in an underdense plasma

Directory of Open Access Journals (Sweden)

S Darvish Molla

2011-06-01

Full Text Available One of the interesting Laser-Plasma phenomena, when the laser power is high and ultra intense, is the generation of large amplitude plasma waves (Wakefield and electron acceleration. An intense electromagnetic laser pulse can create plasma oscillations through the action of the nonlinear pondermotive force. electrons trapped in the wake can be accelerated to high energies, more than 1 TW. Of the wide variety of methods for generating a regular electric field in plasmas with strong laser radiation, the most attractive one at the present time is the scheme of the Laser Wake Field Accelerator (LWFA. In this method, a strong Langmuir wave is excited in the plasma. In such a wave, electrons are trapped and can acquire relativistic energies, accelerated to high energies. In this paper the PIC simulation of wakefield generation and electron acceleration in an underdense plasma with a short ultra intense laser pulse is discussed. 2D electromagnetic PIC code is written by FORTRAN 90, are developed, and the propagation of different electromagnetic waves in vacuum and plasma is shown. Next, the accuracy of implementation of 2D electromagnetic code is verified, making it relativistic and simulating the generating of wakefield and electron acceleration in an underdense plasma. It is shown that when a symmetric electromagnetic pulse passes through the plasma, the longitudinal field generated in plasma, at the back of the pulse, is weaker than the one due to an asymmetric electromagnetic pulse, and thus the electrons acquire less energy. About the asymmetric pulse, when front part of the pulse has smaller time rise than the back part of the pulse, a stronger wakefield generates, in plasma, at the back of the pulse, and consequently the electrons acquire more energy. In an inverse case, when the rise time of the back part of the pulse is bigger in comparison with that of the back part, a weaker wakefield generates and this leads to the fact that the electrons

A gridding method for object-oriented PIC codes

International Nuclear Information System (INIS)

Gisler, G.; Peter, W.; Nash, H.; Acquah, J.; Lin, C.; Rine, D.

1993-01-01

A simple, rule-based gridding method for object-oriented PIC codes is described which is not only capable of dealing with complicated structures such as multiply-connected regions, but is also computationally faster than classical gridding techniques. Using, these smart grids, vacant cells (e.g., cells enclosed by conductors) will never have to be stored or calculated, thus avoiding the usual situation of having to zero electromagnetic fields within conductors after valuable cpu time has been spent in calculating the fields within these cells in the first place. This object-oriented gridding technique makes use of encapsulating characteristics of actual physical objects (particles, fields, grids, etc.) in C ++ classes and supporting software reuse of these entities through C ++ class inheritance relations. It has been implemented in the form of a simple two-dimensional plasma particle-in-cell code, and forms the initial effort of an AFOSR research project to develop a flexible software simulation environment for particle-in-cell algorithms based on object-oriented technology

Electromagnetic particle-in-cell (PIC) method for modeling the formation of metal surface structures induced by femtosecond laser radiation

Energy Technology Data Exchange (ETDEWEB)

Djouder, M. [Laboratoire de Physique et Chimie Quantique, Université Mouloud Mammeri de Tizi-ouzou, BP 17 RP, 15000 Tizi-Ouzou (Algeria); Lamrous, O., E-mail: omarlamrous@mail.ummto.dz [Laboratoire de Physique et Chimie Quantique, Université Mouloud Mammeri de Tizi-ouzou, BP 17 RP, 15000 Tizi-Ouzou (Algeria); Mitiche, M.D. [Laboratoire de Physique et Chimie Quantique, Université Mouloud Mammeri de Tizi-ouzou, BP 17 RP, 15000 Tizi-Ouzou (Algeria); Itina, T.E. [Laboratoire Hubert Curien, UMR CNRS 5516/Université Jean Monnet, 18 rue de Professeur Benoît Lauras, 42000 Saint-Etienne (France); Zemirli, M. [Laboratoire de Physique et Chimie Quantique, Université Mouloud Mammeri de Tizi-ouzou, BP 17 RP, 15000 Tizi-Ouzou (Algeria)

2013-09-01

The particle in cell (PIC) method coupled to the finite-difference time-domain (FDTD) method is used to model the formation of laser induced periodic surface structures (LIPSS) at the early stage of femtosecond laser irradiation of smooth metal surface. The theoretical results were analyzed and compared with experimental data taken from the literature. It was shown that the optical properties of the target are not homogeneous and the ejection of electrons is such that ripples in the electron density were obtained. The Coulomb explosion mechanism was proposed to explain the ripples formation under the considered conditions.

Electromagnetic particle-in-cell (PIC) method for modeling the formation of metal surface structures induced by femtosecond laser radiation

International Nuclear Information System (INIS)

Djouder, M.; Lamrous, O.; Mitiche, M.D.; Itina, T.E.; Zemirli, M.

2013-01-01

The particle in cell (PIC) method coupled to the finite-difference time-domain (FDTD) method is used to model the formation of laser induced periodic surface structures (LIPSS) at the early stage of femtosecond laser irradiation of smooth metal surface. The theoretical results were analyzed and compared with experimental data taken from the literature. It was shown that the optical properties of the target are not homogeneous and the ejection of electrons is such that ripples in the electron density were obtained. The Coulomb explosion mechanism was proposed to explain the ripples formation under the considered conditions.

Canonical symplectic particle-in-cell method for long-term large-scale simulations of the Vlasov–Maxwell equations

Energy Technology Data Exchange (ETDEWEB)

Qin, Hong; Liu, Jian; Xiao, Jianyuan; Zhang, Ruili; He, Yang; Wang, Yulei; Sun, Yajuan; Burby, Joshua W.; Ellison, Leland; Zhou, Yao

2015-12-14

Particle-in-cell (PIC) simulation is the most important numerical tool in plasma physics. However, its long-term accuracy has not been established. To overcome this difficulty, we developed a canonical symplectic PIC method for the Vlasov-Maxwell system by discretising its canonical Poisson bracket. A fast local algorithm to solve the symplectic implicit time advance is discovered without root searching or global matrix inversion, enabling applications of the proposed method to very large-scale plasma simulations with many, e.g. 10(9), degrees of freedom. The long-term accuracy and fidelity of the algorithm enables us to numerically confirm Mouhot and Villani's theory and conjecture on nonlinear Landau damping over several orders of magnitude using the PIC method, and to calculate the nonlinear evolution of the reflectivity during the mode conversion process from extraordinary waves to Bernstein waves.

Review of particle-in-cell modeling for the extraction region of large negative hydrogen ion sources for fusion

Science.gov (United States)

Wünderlich, D.; Mochalskyy, S.; Montellano, I. M.; Revel, A.

2018-05-01

Particle-in-cell (PIC) codes are used since the early 1960s for calculating self-consistently the motion of charged particles in plasmas, taking into account external electric and magnetic fields as well as the fields created by the particles itself. Due to the used very small time steps (in the order of the inverse plasma frequency) and mesh size, the computational requirements can be very high and they drastically increase with increasing plasma density and size of the calculation domain. Thus, usually small computational domains and/or reduced dimensionality are used. In the last years, the available central processing unit (CPU) power strongly increased. Together with a massive parallelization of the codes, it is now possible to describe in 3D the extraction of charged particles from a plasma, using calculation domains with an edge length of several centimeters, consisting of one extraction aperture, the plasma in direct vicinity of the aperture, and a part of the extraction system. Large negative hydrogen or deuterium ion sources are essential parts of the neutral beam injection (NBI) system in future fusion devices like the international fusion experiment ITER and the demonstration reactor (DEMO). For ITER NBI RF driven sources with a source area of 0.9 × 1.9 m2 and 1280 extraction apertures will be used. The extraction of negative ions is accompanied by the co-extraction of electrons which are deflected onto an electron dump. Typically, the maximum negative extracted ion current is limited by the amount and the temporal instability of the co-extracted electrons, especially for operation in deuterium. Different PIC codes are available for the extraction region of large driven negative ion sources for fusion. Additionally, some effort is ongoing in developing codes that describe in a simplified manner (coarser mesh or reduced dimensionality) the plasma of the whole ion source. The presentation first gives a brief overview of the current status of the ion

Rotating structures in low temperature magnetized plasmas - Insight from particle simulations

Directory of Open Access Journals (Sweden)

Jean-Pierre eBoeuf

2014-12-01

Full Text Available The EXB configuration of various low temperature plasma devices is often responsible for the formation of rotating structures and instabilities leading to anomalous electron transport across the magnetic field. In these devices, electrons are strongly magnetized while ions are weakly or not magnetized and this leads to specific physical phenomena that are not present in fusion plasmas where both electrons and ions are strongly magnetized. In this paper we describe basic phenomena involving rotating plasma structures in simple configurations of low temperature EXB plasma devices on the basis of PIC-MCC (Particle-In-Cell Monte Carlo Collisions simulations. We focus on three examples: rotating electron vortices and rotating spokes in cylindrical magnetrons, and azimuthal electron-cyclotron drift instability in Hall thrusters. The simulations are not intended to give definite answers to the many physics issues related to low temperature EXB plasma devices but are used to illustrate and discuss some of the basic questions that need further studies.

Sparse grid techniques for particle-in-cell schemes

Science.gov (United States)

Ricketson, L. F.; Cerfon, A. J.

2017-02-01

We propose the use of sparse grids to accelerate particle-in-cell (PIC) schemes. By using the so-called ‘combination technique’ from the sparse grids literature, we are able to dramatically increase the size of the spatial cells in multi-dimensional PIC schemes while paying only a slight penalty in grid-based error. The resulting increase in cell size allows us to reduce the statistical noise in the simulation without increasing total particle number. We present initial proof-of-principle results from test cases in two and three dimensions that demonstrate the new scheme’s efficiency, both in terms of computation time and memory usage.

Final Report for grant ER54958, 'Gyrokinetic Particle Simulation of Turbulent Transport in Burning Plasmas'

International Nuclear Information System (INIS)

Decyk, Viktor K.

2011-01-01

widely used in plasma modeling in DOE, not only in areas in fusion energy as exemplified by GTC, but in high energy physics, plasma accelerators, ICF, and other areas. In 2010, about 12% of the INCITE grants in DOE were devoted to PIC codes. We began by developing a simple 2D electrostatic PIC code for the NVIDIA Tesla C1060 GPU based on one of the codes from the UPIC Framework. The major new feature of this code was the implementation of a streaming algorithm, where the two major data elements (particles and fields) are read only once each time step with an optimal memory access pattern (unit stride). To achieve this, particles need to be ordered by cell and we developed a particle reordering scheme that worked effectively on this hardware. The first results used global memory only and achieved a speedup of 13 compared to a 2.67 GHz Intel Nehalem processor, and were published in the ICAP conference proceedings in 2009. In the next version, we parameterized the code to make it adaptable to different architectures. The reordering algorithm was generalized to allow more than one grid per sorting cell and more than one sorting cell per thread. We also added support for shared memory. The four tunable parameters were defined as follows: lth, the number of tightly coupled threads, ngpx and ngpy, the number of grids in a sorting cell, and ngpt, the number of sorting cells assigned to a thread. Increasing the number of grids per sorting cell reduced the cost of particle reordering, but it could increase the particle processing time because more shared memory was required. For the NVIDIA C1060, the optimal parameters were lth=32, ngpx=2, ngpy=3, ngpt = 1. Speedups of 15-25 were obtained compared to the Intel Nehalem, depending on plasma temperature. Details about the algorithm and performance results were published in 2011. The electrostatic PIC code was very simple, with low computational intensity (few floating point operations per memory access). Codes with higher

Modeling of the plasma generated in a rarefied hypersonic shock layer

International Nuclear Information System (INIS)

Farbar, Erin D.; Boyd, Iain D.

2010-01-01

In this study, a rigorous numerical model is developed to simulate the plasma generated in a rarefied, hypersonic shock layer. The model uses the direct simulation Monte Carlo (DSMC) method to treat the particle collisions and the particle-in-cell (PIC) method to simulate the plasma dynamics in a self-consistent manner. The model is applied to compute the flow along the stagnation streamline in front of a blunt body reentering the Earth's atmosphere at very high velocity. Results from the rigorous DSMC-PIC model are compared directly to the standard DSMC modeling approach that uses the ambipolar diffusion approximation to simulate the plasma dynamics. It is demonstrated that the self-consistent computation of the plasma dynamics using the rigorous DSMC-PIC model captures many physical phenomena not accurately predicted by the standard modeling approach. These computations represent the first assessment of the validity of the ambipolar diffusion approximation when predicting the rarefied plasma generated in a hypersonic shock layer.

Multi-dimensional PIC-simulations of parametric instabilities for shock-ignition conditions

Directory of Open Access Journals (Sweden)

Riconda C.

2013-11-01

Full Text Available Laser-plasma interaction is investigated for conditions relevant for the shock-ignition (SI scheme of inertial confinement fusion using two-dimensional particle-in-cell (PIC simulations of an intense laser beam propagating in a hot, large-scale, non-uniform plasma. The temporal evolution and interdependence of Raman- (SRS, and Brillouin- (SBS, side/backscattering as well as Two-Plasmon-Decay (TPD are studied. TPD is developing in concomitance with SRS creating a broad spectrum of plasma waves near the quarter-critical density. They are rapidly saturated due to plasma cavitation within a few picoseconds. The hot electron spectrum created by SRS and TPD is relatively soft, limited to energies below one hundred keV.

Electron acceleration in the Solar corona - 3D PiC code simulations of guide field reconnection

Science.gov (United States)

Alejandro Munoz Sepulveda, Patricio

2017-04-01

The efficient electron acceleration in the solar corona detected by means of hard X-ray emission is still not well understood. Magnetic reconnection through current sheets is one of the proposed production mechanisms of non-thermal electrons in solar flares. Previous works in this direction were based mostly on test particle calculations or 2D fully-kinetic PiC simulations. We have now studied the consequences of self-generated current-aligned instabilities on the electron acceleration mechanisms by 3D magnetic reconnection. For this sake, we carried out 3D Particle-in-Cell (PiC) code numerical simulations of force free reconnecting current sheets, appropriate for the description of the solar coronal plasmas. We find an efficient electron energization, evidenced by the formation of a non-thermal power-law tail with a hard spectral index smaller than -2 in the electron energy distribution function. We discuss and compare the influence of the parallel electric field versus the curvature and gradient drifts in the guiding-center approximation on the overall acceleration, and their dependence on different plasma parameters.

PIC simulations of conical magnetically insulated transmission line with LTD generator: Transition from self-limited to load-limited flow

Science.gov (United States)

Liu, Laqun; Wang, Huihui; Guo, Fan; Zou, Wenkang; Liu, Dagang

2017-04-01

Based on the 3-dimensional Particle-In-Cell (PIC) code CHIPIC3D, with a new circuit boundary algorithm we developed, a conical magnetically insulated transmission line (MITL) with a 1.0-MV linear transformer driver (LTD) is explored numerically. The values of switch jitter time of LTD are critical parameters for the system, which are difficult to be measured experimentally. In this paper, these values are obtained by comparing the PIC results with experimental data of large diode-gap MITL. By decreasing the diode gap, we find that all PIC results agree well with experimental data only if MITL works on self-limited flow no matter how large the diode gap is. However, when the diode gap decreases to a threshold, the self-limited flow would transfer to a load-limited flow. In this situation, PIC results no longer agree with experimental data anymore due to the anode plasma expansion in the diode load. This disagreement is used to estimate the plasma expansion speed.

Analysis of instability growth and collisionless relaxation in thermionic converters using 1-D PIC simulations

International Nuclear Information System (INIS)

Kreh, B.B.

1994-12-01

This work investigates the role that the beam-plasma instability may play in a thermionic converter. The traditional assumption of collisionally dominated relaxation is questioned, and the beam-plasma instability is proposed as a possible dominant relaxation mechanism. Theory is developed to describe the beam-plasma instability in the cold-plasma approximation, and the theory is tested with two common Particle-in-Cell (PIC) simulation codes. The theory is first confirmed using an unbounded plasma PIC simulation employing periodic boundary conditions, ES1. The theoretically predicted growth rates are on the order of the plasma frequencies, and ES1 simulations verify these predictions within the order of 1%. For typical conditions encountered in thermionic converters, the resulting growth period is on the order of 7 x 10 -11 seconds. The bounded plasma simulation PDP1 was used to evaluate the influence of finite geometry and the electrode boundaries. For this bounded plasma, a two-stream interaction was supported and resulting in nearly complete thermalization in approximately 5 x 10 -10 seconds. Since the electron-electron collision rate of 10 9 Hz and the electron atom collision rate of 10 7 Hz are significantly slower than the rate of development of these instabilities, the instabilities appear to be an important relaxation mechanism

photon-plasma: A modern high-order particle-in-cell code

International Nuclear Information System (INIS)

Haugbølle, Troels; Frederiksen, Jacob Trier; Nordlund, Åke

2013-01-01

We present the photon-plasma code, a modern high order charge conserving particle-in-cell code for simulating relativistic plasmas. The code is using a high order implicit field solver and a novel high order charge conserving interpolation scheme for particle-to-cell interpolation and charge deposition. It includes powerful diagnostics tools with on-the-fly particle tracking, synthetic spectra integration, 2D volume slicing, and a new method to correctly account for radiative cooling in the simulations. A robust technique for imposing (time-dependent) particle and field fluxes on the boundaries is also presented. Using a hybrid OpenMP and MPI approach, the code scales efficiently from 8 to more than 250.000 cores with almost linear weak scaling on a range of architectures. The code is tested with the classical benchmarks particle heating, cold beam instability, and two-stream instability. We also present particle-in-cell simulations of the Kelvin-Helmholtz instability, and new results on radiative collisionless shocks

Ion drag force on dust grains in the magnetized edge plasma

International Nuclear Information System (INIS)

Matyash, K.; Schneider, R.; Ikkurthi, V.R.; Melzer, A.

2009-01-01

A 3-dimensional Particle-Particle Particle-Mesh (P3M) code [K. Matyash, R. Schneider, F. Taccogna, D. Tskhakaya, J. Nucl. Mater. 363-365 (2007) 458] is applied to simulate a small-size (smaller than a Debye length) spherical dust grain confined in the magnetized plasma near the material wall of a tokamak. Plasma particles (electrons and ions) are treated kinetically (Particle-in-Cell with Monte Carlo Collisions (PIC MCC)), which allows to resolve self-consistently the electrostatic sheath in front of the wall. In order to describe accurately the plasma particles' motion close to the dust grain, the PIC technique is supplemented with Molecular Dynamics (MD), employing an analytic electrostatic potential for the interaction with the dust grain. The charging of a spherical, conducting dust grain confined in the sheath potential close to the wall of a tokamak is simulated. A magnetic field normal to the wall was investigated. The ion drag force resulting from dust grain collisions with the streaming ions is calculated. This force is critical for a realistic description of the dust particle dynamics and transport in fusion plasmas.

SMILEI: A collaborative, open-source, multi-purpose PIC code for the next generation of super-computers

Science.gov (United States)

Grech, Mickael; Derouillat, J.; Beck, A.; Chiaramello, M.; Grassi, A.; Niel, F.; Perez, F.; Vinci, T.; Fle, M.; Aunai, N.; Dargent, J.; Plotnikov, I.; Bouchard, G.; Savoini, P.; Riconda, C.

2016-10-01

Over the last decades, Particle-In-Cell (PIC) codes have been central tools for plasma simulations. Today, new trends in High-Performance Computing (HPC) are emerging, dramatically changing HPC-relevant software design and putting some - if not most - legacy codes far beyond the level of performance expected on the new and future massively-parallel super computers. SMILEI is a new open-source PIC code co-developed by both plasma physicists and HPC specialists, and applied to a wide range of physics-related studies: from laser-plasma interaction to astrophysical plasmas. It benefits from an innovative parallelization strategy that relies on a super-domain-decomposition allowing for enhanced cache-use and efficient dynamic load balancing. Beyond these HPC-related developments, SMILEI also benefits from additional physics modules allowing to deal with binary collisions, field and collisional ionization and radiation back-reaction. This poster presents the SMILEI project, its HPC capabilities and illustrates some of the physics problems tackled with SMILEI.

Wavenumber spectrum of whistler turbulence: Particle-in-cell simulation

International Nuclear Information System (INIS)

Saito, S.; Gary, S. Peter; Narita, Y.

2010-01-01

The forward cascade of decaying whistler turbulence is studied in low beta plasma to understand essential properties of the energy spectrum at electron scales, by using a two-dimensional electromagnetic particle-in-cell (PIC) simulation. This simulation demonstrates turbulence in which the energy cascade rate is greater than the dissipation rate at the electron inertial length. The PIC simulation shows that the magnetic energy spectrum of forward-cascaded whistler turbulence at electron inertial scales is anisotropic and develops a very steep power-law spectrum which is consistent with recent solar wind observations. A comparison of the simulated spectrum with that predicted by a phenomenological turbulence scaling model suggests that the energy cascade at the electron inertial scale depends on both magnetic fluctuations and electron velocity fluctuations, as well as on the whistler dispersion relation. Thus, not only kinetic Alfven turbulence but also whistler turbulence may explain recent solar wind observations of very steep magnetic spectra at short scales.

Plasma non-uniformity in a symmetric radiofrequency capacitively-coupled reactor with dielectric side-wall: a two dimensional particle-in-cell/Monte Carlo collision simulation

Science.gov (United States)

Liu, Yue; Booth, Jean-Paul; Chabert, Pascal

2018-02-01

A Cartesian-coordinate two-dimensional electrostatic particle-in-cell/Monte Carlo collision (PIC/MCC) plasma simulation code is presented, including a new treatment of charge balance at dielectric boundaries. It is used to simulate an Ar plasma in a symmetric radiofrequency capacitively-coupled parallel-plate reactor with a thick (3.5 cm) dielectric side-wall. The reactor size (12 cm electrode width, 2.5 cm electrode spacing) and frequency (15 MHz) are such that electromagnetic effects can be ignored. The dielectric side-wall effectively shields the plasma from the enhanced electric field at the powered-grounded electrode junction, which has previously been shown to produce locally enhanced plasma density (Dalvie et al 1993 Appl. Phys. Lett. 62 3207-9 Overzet and Hopkins 1993 Appl. Phys. Lett. 63 2484-6 Boeuf and Pitchford 1995 Phys. Rev. E 51 1376-90). Nevertheless, enhanced electron heating is observed in a region adjacent to the dielectric boundary, leading to maxima in ionization rate, plasma density and ion flux to the electrodes in this region, and not at the reactor centre as would otherwise be expected. The axially-integrated electron power deposition peaks closer to the dielectric edge than the electron density. The electron heating components are derived from the PIC/MCC simulations and show that this enhanced electron heating results from increased Ohmic heating in the axial direction as the electron density decreases towards the side-wall. We investigated the validity of different analytical formulas to estimate the Ohmic heating by comparing them to the PIC results. The widespread assumption that a time-averaged momentum transfer frequency, v m , can be used to estimate the momentum change can cause large errors, since it neglects both phase and amplitude information. Furthermore, the classical relationship between the total electron current and the electric field must be used with caution, particularly close to the dielectric edge where the (neglected

A new approach to theoretical investigations of high harmonics generation by means of fs laser interaction with overdense plasma layers. Combining particle-in-cell simulations with machine learning

International Nuclear Information System (INIS)

Mihailescu, A.

2016-01-01

Within the past decade, various experimental and theoretical investigations have been performed in the field of high-order harmonics generation (HHG) by means of femtosecond ( fs ) laser pulses interacting with laser produced plasmas. Numerous potential future applications thus arise. Beyond achieving higher conversion efficiency for higher harmonic orders and hence harmonic power and brilliance, there are more ambitious scientific goals such as attaining shorter harmonic wavelengths or reducing harmonic pulse durations towards the attosecond and even the zeptosecond range. High order harmonics are also an attractive diagnostic tool for the laser-plasma interaction process itself. Particle-in-Cell (PIC) simulations are known to be one of the most important numerical instruments employed in plasma physics and in laser-plasma interaction investigations. The novelty brought by this paper consists in combining the PIC method with several machine learning approaches. For predictive modelling purposes, a universal functional approximator is used, namely a multi-layer perceptron (MLP), in conjunction with a self-organizing map (SOM). The training sets have been retrieved from the PIC simulations and also from the available literature in the field. The results demonstrate the potential utility of machine learning in predicting optimal interaction scenarios for gaining higher order harmonics or harmonics with particular features such as a particular wavelength range, a particular harmonic pulse duration or a certain intensity. Furthermore, the author will show how machine learning can be used for estimations of electronic temperatures, proving that it can be a reliable tool for obtaining better insights into the fs laser interaction physics.

Particle-in-Cell Laser-Plasma Simulation on Xeon Phi Coprocessors

OpenAIRE

Surmin, I. A.; Bastrakov, S. I.; Efimenko, E. S.; Gonoskov, A. A.; Korzhimanov, A. V.; Meyerov, I. B.

2015-01-01

This paper concerns development of a high-performance implementation of the Particle-in-Cell method for plasma simulation on Intel Xeon Phi coprocessors. We discuss suitability of the method for Xeon Phi architecture and present our experience of porting and optimization of the existing parallel Particle-in-Cell code PICADOR. Direct porting with no code modification gives performance on Xeon Phi close to 8-core CPU on a benchmark problem with 50 particles per cell. We demonstrate step-by-step...

Performance of particle in cell methods on highly concurrent computational architectures

International Nuclear Information System (INIS)

Adams, M.F.; Ethier, S.; Wichmann, N.

2009-01-01

Particle in cell (PIC) methods are effective in computing Vlasov-Poisson system of equations used in simulations of magnetic fusion plasmas. PIC methods use grid based computations, for solving Poisson's equation or more generally Maxwell's equations, as well as Monte-Carlo type methods to sample the Vlasov equation. The presence of two types of discretizations, deterministic field solves and Monte-Carlo methods for the Vlasov equation, pose challenges in understanding and optimizing performance on today large scale computers which require high levels of concurrency. These challenges arises from the need to optimize two very different types of processes and the interactions between them. Modern cache based high-end computers have very deep memory hierarchies and high degrees of concurrency which must be utilized effectively to achieve good performance. The effective use of these machines requires maximizing concurrency by eliminating serial or redundant work and minimizing global communication. A related issue is minimizing the memory traffic between levels of the memory hierarchy because performance is often limited by the bandwidths and latencies of the memory system. This paper discusses some of the performance issues, particularly in regard to parallelism, of PIC methods. The gyrokinetic toroidal code (GTC) is used for these studies and a new radial grid decomposition is presented and evaluated. Scaling of the code is demonstrated on ITER sized plasmas with up to 16K Cray XT3/4 cores.

Performance of particle in cell methods on highly concurrent computational architectures

International Nuclear Information System (INIS)

Adams, M F; Ethier, S; Wichmann, N

2007-01-01

Particle in cell (PIC) methods are effective in computing Vlasov-Poisson system of equations used in simulations of magnetic fusion plasmas. PIC methods use grid based computations, for solving Poisson's equation or more generally Maxwell's equations, as well as Monte-Carlo type methods to sample the Vlasov equation. The presence of two types of discretizations, deterministic field solves and Monte-Carlo methods for the Vlasov equation, pose challenges in understanding and optimizing performance on today large scale computers which require high levels of concurrency. These challenges arises from the need to optimize two very different types of processes and the interactions between them. Modern cache based high-end computers have very deep memory hierarchies and high degrees of concurrency which must be utilized effectively to achieve good performance. The effective use of these machines requires maximizing concurrency by eliminating serial or redundant work and minimizing global communication. A related issue is minimizing the memory traffic between levels of the memory hierarchy because performance is often limited by the bandwidths and latencies of the memory system. This paper discusses some of the performance issues, particularly in regard to parallelism, of PIC methods. The gyrokinetic toroidal code (GTC) is used for these studies and a new radial grid decomposition is presented and evaluated. Scaling of the code is demonstrated on ITER sized plasmas with up to 16K Cray XT3/4 cores

FLIP-MHD: A particle-in-cell mehtod for magnetohydrodynamics

International Nuclear Information System (INIS)

Brackbill, J.U.

1990-01-01

A particle-in-cell (PIC) method, FLIP is extended to magnetohydrodynamic (MHD) flow in two dimensions. Particles are used to reduce computational diffusion of the magnetic field. FLIP is an extension of ''classical'' PIC, where particles have mass, but every other property of the fluid is stored on a grid. In FLIP, particles have every property of the fluid, so that they provide a complete Lagrangian description not only to resolve contact discontinuities but also to reduce computational diffusion of linear and angular momentum. The interactions among the particles are calculated on a grid, for convenience and economy. The present study extends FLIP to MHD, by including information about the magnetic field among the attributes of the particles. 6 refs

PIC Modeling of Argon Plasma Flow in MNX

Science.gov (United States)

Cohen, Samuel; Sefkow, Adam

2007-11-01

A linear helicon-heated plasma device - the Magnetic Nozzle Experiment (MNX) at the Princeton Plasma Physics Laboratory - is used for studies of the formation of strong electrostatic double layers near mechanical and magnetic apertures and the acceleration of plasma ions into supersonic directed beams. In order to characterize the role of the aperture and its involvement with ion acceleration, detailed particle-in-cell simulations are employed to study the effects of the surrounding boundary geometry on the plasma dynamics near the aperture region, within which the transition from a collisional to collisionless regime occurs. The presence of a small superthermal electron population is examined, and the model includes a background neutral population which can be ionized by energetic electrons. By self-consistently evaluating the temporal evolution of the plasma in the vicinity of the aperture, the formation mechanism of the double layer is investigated.

Output power fluctuations due to different weights of macro particles used in particle-in-cell simulations of Cerenkov devices

International Nuclear Information System (INIS)

Bao, Rong; Li, Yongdong; Liu, Chunliang; Wang, Hongguang

2016-01-01

The output power fluctuations caused by weights of macro particles used in particle-in-cell (PIC) simulations of a backward wave oscillator and a travelling wave tube are statistically analyzed. It is found that the velocities of electrons passed a specific slow-wave structure form a specific electron velocity distribution. The electron velocity distribution obtained in PIC simulation with a relative small weight of macro particles is considered as an initial distribution. By analyzing this initial distribution with a statistical method, the estimations of the output power fluctuations caused by different weights of macro particles are obtained. The statistical method is verified by comparing the estimations with the simulation results. The fluctuations become stronger with increasing weight of macro particles, which can also be determined reversely from estimations of the output power fluctuations. With the weights of macro particles optimized by the statistical method, the output power fluctuations in PIC simulations are relatively small and acceptable.

Energy loss of heavy ion beams in plasma

Energy Technology Data Exchange (ETDEWEB)

Okada, T; Hotta, T [Tokyo Univ. of Agriculture and Technology, Koganei (Japan). Faculty of Technology

1997-12-31

The energy loss of heavy-ion beams (HIB) is studied by means of Vlasov theory and Particle-in-Cell (PIC) simulations in a plasma. The interaction of HIB with a plasma is of central importance for inertial confinement fusion (ICF). A number of studies on the HIB interaction with target plasma have been published. It is important for heavy-ion stopping that the effects of the non-linear interaction of HIB within the Vlasov theory are included. Reported are results of a numerical study of nonlinear effects to the stopping power for HIB in plasma. It is shown that the PIC simulations of collective effects of the stopping power are in a good agreement with the Vlasov theory. (author). 2 tabs., 1 fig., 5 refs.

An incompressible two-dimensional multiphase particle-in-cell model for dense particle flows

Energy Technology Data Exchange (ETDEWEB)

Snider, D.M. [SAIC, Albuquerque, NM (United States); O`Rourke, P.J. [Los Alamos National Lab., NM (United States); Andrews, M.J. [Texas A and M Univ., College Station, TX (United States). Dept. of Mechanical Engineering

1997-06-01

A two-dimensional, incompressible, multiphase particle-in-cell (MP-PIC) method is presented for dense particle flows. The numerical technique solves the governing equations of the fluid phase using a continuum model and those of the particle phase using a Lagrangian model. Difficulties associated with calculating interparticle interactions for dense particle flows with volume fractions above 5% have been eliminated by mapping particle properties to a Eulerian grid and then mapping back computed stress tensors to particle positions. This approach utilizes the best of Eulerian/Eulerian continuum models and Eulerian/Lagrangian discrete models. The solution scheme allows for distributions of types, sizes, and density of particles, with no numerical diffusion from the Lagrangian particle calculations. The computational method is implicit with respect to pressure, velocity, and volume fraction in the continuum solution thus avoiding courant limits on computational time advancement. MP-PIC simulations are compared with one-dimensional problems that have analytical solutions and with two-dimensional problems for which there are experimental data.

2D PIC simulations for an EN discharge with magnetized electrons and unmagnetized ions

Science.gov (United States)

Lieberman, Michael A.; Kawamura, Emi; Lichtenberg, Allan J.

2009-10-01

We conducted 2D particle-in-cell (PIC) simulations for an electronegative (EN) discharge with magnetized electrons and unmagnetized ions, and compared the results to a previously developed 1D (radial) analytical model of an EN plasma with strongly magnetized electrons and weakly magnetized ions [1]. In both cases, there is a static uniform applied magnetic field in the axial direction. The 1D radial model mimics the wall losses of the particles in the axial direction by introducing a bulk loss frequency term νL. A special (desired) solution was found in which only positive and negative ions but no electrons escaped radially. The 2D PIC results show good agreement with the 1D model over a range of parameters and indicate that the analytical form of νL employed in [1] is reasonably accurate. However, for the PIC simulations, there is always a finite flux of electrons to the radial wall which is about 10 to 30% of the negative ion flux.[4pt] [1] G. Leray, P. Chabert, A.J. Lichtenberg and M.A. Lieberman, J. Phys. D, accepted for publication 2009.

Loading relativistic Maxwell distributions in particle simulations

Science.gov (United States)

Zenitani, S.

2015-12-01

In order to study energetic plasma phenomena by using particle-in-cell (PIC) and Monte-Carlo simulations, we need to deal with relativistic velocity distributions in these simulations. However, numerical algorithms to deal with relativistic distributions are not well known. In this contribution, we overview basic algorithms to load relativistic Maxwell distributions in PIC and Monte-Carlo simulations. For stationary relativistic Maxwellian, the inverse transform method and the Sobol algorithm are reviewed. To boost particles to obtain relativistic shifted-Maxwellian, two rejection methods are newly proposed in a physically transparent manner. Their acceptance efficiencies are 􏰅50% for generic cases and 100% for symmetric distributions. They can be combined with arbitrary base algorithms.

High-performance modeling of plasma-based acceleration and laser-plasma interactions

Science.gov (United States)

Vay, Jean-Luc; Blaclard, Guillaume; Godfrey, Brendan; Kirchen, Manuel; Lee, Patrick; Lehe, Remi; Lobet, Mathieu; Vincenti, Henri

2016-10-01

Large-scale numerical simulations are essential to the design of plasma-based accelerators and laser-plasma interations for ultra-high intensity (UHI) physics. The electromagnetic Particle-In-Cell (PIC) approach is the method of choice for self-consistent simulations, as it is based on first principles, and captures all kinetic effects, and also scale favorably to many cores on supercomputers. The standard PIC algorithm relies on second-order finite-difference discretization of the Maxwell and Newton-Lorentz equations. We present here novel formulations, based on very high-order pseudo-spectral Maxwell solvers, which enable near-total elimination of the numerical Cherenkov instability and increased accuracy over the standard PIC method for standard laboratory frame and Lorentz boosted frame simulations. We also present the latest implementations in the PIC modules Warp-PICSAR and FBPIC on the Intel Xeon Phi and GPU architectures. Examples of applications will be given on the simulation of laser-plasma accelerators and high-harmonic generation with plasma mirrors. Work supported by US-DOE Contracts DE-AC02-05CH11231 and by the European Commission through the Marie Slowdoska-Curie fellowship PICSSAR Grant Number 624543. Used resources of NERSC.

Fluctuations and transport in fusion plasmas. Final report

International Nuclear Information System (INIS)

Gould, R.W.; Liewer, P.C.

1995-01-01

The energy confinement in tokamaks in thought to be limited by transport caused by plasma turbulence. Three dimensional plasma particle-in-cell (PIC) codes are used to model the turbulent transport in tokamaks to attempt to understand this phenomena so that tokamaks can be made more efficient. Presently, hundreds of hours of Cray time are used to model these experiments and much bigger and longer runs are desired, to model a large tokamak with realistic parameters is beyond the capability of existing sequential supercomputers. Parallel supercomputers might be a cost effect tool for performing such large scale 3D tokamak simulations. The goal of the work was to develop algorithms for performing PIC codes on coarse-grained message passing parallel computers and to evaluate the performance of such parallel computers on PIC codes. This algorithm would be used in a large scale PIC production code such as the UCLA 3D gyrokinetic code

One-dimensional hybrid-direct kinetic simulation of the discharge plasma in a Hall thruster

International Nuclear Information System (INIS)

Hara, Kentaro; Boyd, Iain D.; Kolobov, Vladimir I.

2012-01-01

In order to model the non-equilibrium plasma within the discharge region of a Hall thruster, the velocity distribution functions (VDFs) must be obtained accurately. A direct kinetic (DK) simulation method that directly solves the plasma Boltzmann equation can achieve better resolution of VDFs in comparison to particle simulations, such as the particle-in-cell (PIC) method that inherently include statistical noise. In this paper, a one-dimensional hybrid-DK simulation, which uses a DK simulation for heavy species and a fluid model for electrons, is developed and compared to a hybrid-PIC simulation. Time-averaged results obtained from the hybrid-DK simulation are in good agreement with hybrid-PIC results and experimental data. It is shown from a comparison of using a kinetic simulation and solving the continuity equation that modeling of the neutral atoms plays an important role for simulations of the Hall thruster discharge plasma. In addition, low and high frequency plasma oscillations are observed. Although the kinetic nature of electrons is not resolved due to the use of a fluid model, the hybrid-DK model provides spatially and temporally well-resolved plasma properties and an improved resolution of VDFs for heavy species with less statistical noise in comparison to the hybrid-PIC method.

End-to-end plasma bubble PIC simulations on GPUs

Science.gov (United States)

Germaschewski, Kai; Fox, William; Matteucci, Jackson; Bhattacharjee, Amitava

2017-10-01

Accelerator technologies play a crucial role in eventually achieving exascale computing capabilities. The current and upcoming leadership machines at ORNL (Titan and Summit) employ Nvidia GPUs, which provide vast computational power but also need specifically adapted computational kernels to fully exploit them. In this work, we will show end-to-end particle-in-cell simulations of the formation, evolution and coalescence of laser-generated plasma bubbles. This work showcases the GPU capabilities of the PSC particle-in-cell code, which has been adapted for this problem to support particle injection, a heating operator and a collision operator on GPUs.

Dynamic Load Balancing for PIC code using Eulerian/Lagrangian partitioning

OpenAIRE

Sauget, Marc; Latu, Guillaume

2017-01-01

This document presents an analysis of different load balance strategies for a Plasma physics code that models high energy particle beams with PIC method. A comparison of different load balancing algorithms is given: static or dynamic ones. Lagrangian and Eulerian partitioning techniques have been investigated.

Electromagnetic direct implicit PIC simulation

International Nuclear Information System (INIS)

Langdon, A.B.

1983-01-01

Interesting modelling of intense electron flow has been done with implicit particle-in-cell simulation codes. In this report, the direct implicit PIC simulation approach is applied to simulations that include full electromagnetic fields. The resulting algorithm offers advantages relative to moment implicit electromagnetic algorithms and may help in our quest for robust and simpler implicit codes

Spectral properties and associated plasma energization by magnetosonic waves in the Earth's magnetosphere: Particle-in-cell simulations

Science.gov (United States)

Sun, Jicheng; Gao, Xinliang; Lu, Quanming; Chen, Lunjin; Liu, Xu; Wang, Xueyi; Tao, Xin; Wang, Shui

2017-05-01

In this paper, we perform a 1-D particle-in-cell (PIC) simulation model consisting of three species, cold electrons, cold ions, and energetic ion ring, to investigate spectral structures of magnetosonic waves excited by ring distribution protons in the Earth's magnetosphere, and dynamics of charged particles during the excitation of magnetosonic waves. As the wave normal angle decreases, the spectral range of excited magnetosonic waves becomes broader with upper frequency limit extending beyond the lower hybrid resonant frequency, and the discrete spectra tends to merge into a continuous one. This dependence on wave normal angle is consistent with the linear theory. The effects of magnetosonic waves on the background cold plasma populations also vary with wave normal angle. For exactly perpendicular magnetosonic waves (parallel wave number k|| = 0), there is no energization in the parallel direction for both background cold protons and electrons due to the negligible fluctuating electric field component in the parallel direction. In contrast, the perpendicular energization of background plasmas is rather significant, where cold protons follow unmagnetized motion while cold electrons follow drift motion due to wave electric fields. For magnetosonic waves with a finite k||, there exists a nonnegligible parallel fluctuating electric field, leading to a significant and rapid energization in the parallel direction for cold electrons. These cold electrons can also be efficiently energized in the perpendicular direction due to the interaction with the magnetosonic wave fields in the perpendicular direction. However, cold protons can be only heated in the perpendicular direction, which is likely caused by the higher-order resonances with magnetosonic waves. The potential impacts of magnetosonic waves on the energization of the background cold plasmas in the Earth's inner magnetosphere are also discussed in this paper.

Nonlinear Plasma Waves Excitation by Intense Ion Beams in Background Plasma

International Nuclear Information System (INIS)

Kaganovich, Igor D.; Startsev, Edward A.; Davidson, Ronald C.

2004-01-01

Plasma neutralization of an intense ion pulse is of interest for many applications, including plasma lenses, heavy ion fusion, cosmic ray propagation, etc. An analytical electron fluid model has been developed to describe the plasma response to a propagating ion beam. The model predicts very good charge neutralization during quasi-steady-state propagation, provided the beam pulse duration τ b is much longer than the electron plasma period 2π/ω p , where ω p = (4πe 2 n p /m) 1/2 is the electron plasma frequency and n p is the background plasma density. In the opposite limit, the beam pulse excites large-amplitude plasma waves. If the beam density is larger than the background plasma density, the plasma waves break. Theoretical predictions are compared with the results of calculations utilizing a particle-in-cell (PIC) code. The cold electron fluid results agree well with the PIC simulations for ion beam propagation through a background plasma. The reduced fluid description derived in this paper can provide an important benchmark for numerical codes and yield scaling relations for different beam and plasma parameters. The visualization of numerical simulation data shows complex collective phenomena during beam entry and exit from the plasma

Nonlinear plasma waves excitation by intense ion beams in background plasma

International Nuclear Information System (INIS)

Kaganovich, Igor D.; Startsev, Edward A.; Davidson, Ronald C.

2004-01-01

Plasma neutralization of an intense ion pulse is of interest for many applications, including plasma lenses, heavy ion fusion, cosmic ray propagation, etc. An analytical electron fluid model has been developed to describe the plasma response to a propagating ion beam. The model predicts very good charge neutralization during quasi-steady-state propagation, provided the beam pulse duration τ b is much longer than the electron plasma period 2π/ω p , where ω p =(4πe 2 n p /m) 1/2 is the electron plasma frequency, and n p is the background plasma density. In the opposite limit, the beam pulse excites large-amplitude plasma waves. If the beam density is larger than the background plasma density, the plasma waves break. Theoretical predictions are compared with the results of calculations utilizing a particle-in-cell (PIC) code. The cold electron fluid results agree well with the PIC simulations for ion beam propagation through a background plasma. The reduced fluid description derived in this paper can provide an important benchmark for numerical codes and yield scaling relations for different beam and plasma parameters. The visualization of numerical simulation data shows complex collective phenomena during beam entry and exit from the plasma

Self-modulation instability of a long proton bunch in plasmas

CERN Document Server

Kumar, Naveen; Lotov, Konstantin

2010-01-01

An analytical model for the self-modulation instability of a long relativistic proton bunch propagating in uniform plasmas is developed. The self-modulated proton bunch resonantly excites a large amplitude plasma wave (wake field), which can be used for acceleration of plasma electrons. Analytical expressions for the linear growth rate and the number of exponentiations are given. We use the full three-dimensional particle-in-cell (PIC) simulations to study the beam self-modulation and the transition to the nonlinear stage. It is shown that the self-modulation of the proton bunch competes with the hosing instability which tends to destroy the plasma wave. A method is proposed and studied through PIC simulations to circumvent this problem which relies on the seeding of the self-modulation instability in the bunch.

Thrust calculation of electric solar wind sail by particle-in-cell simulation

Energy Technology Data Exchange (ETDEWEB)

Hoshi, Kento [Kyoto Univ. (Japan). Dept. of Electrical Engineering; Kojima, Hirotsugu; Yamakawa, Hiroshi [Kyoto Univ. (Japan). Research Inst. for Sustainable Humanosphere; Muranaka, Takanobu [Chukyo Univ., Nagoya (Japan). Dept. of Electrical Engineering

2016-07-01

In this study, thrust characteristics of an electric solar wind sail were numerically evaluated using full threedimensional particle-in-cell (PIC) simulation. The thrust obtained from the PIC simulation was lower than the thrust estimations obtained in previous studies. The PIC simulation indicated that ambient electrons strongly shield the electrostatic potential of the tether of the sail, and the strong shield effect causes a greater thrust reduction than has been obtained in previous studies. Additionally, previous expressions of the thrust estimation were modified by using the shielded potential structure derived from the present simulation results. The modified thrust estimation agreed very well with the thrust obtained from the PIC simulation.

Thrust calculation of electric solar wind sail by particle-in-cell simulation

International Nuclear Information System (INIS)

Hoshi, Kento; Kojima, Hirotsugu; Yamakawa, Hiroshi; Muranaka, Takanobu

2016-01-01

In this study, thrust characteristics of an electric solar wind sail were numerically evaluated using full threedimensional particle-in-cell (PIC) simulation. The thrust obtained from the PIC simulation was lower than the thrust estimations obtained in previous studies. The PIC simulation indicated that ambient electrons strongly shield the electrostatic potential of the tether of the sail, and the strong shield effect causes a greater thrust reduction than has been obtained in previous studies. Additionally, previous expressions of the thrust estimation were modified by using the shielded potential structure derived from the present simulation results. The modified thrust estimation agreed very well with the thrust obtained from the PIC simulation.

Particle-in-cell simulations of plasma interaction with shaped and unshaped gaps in TEXTOR

Czech Academy of Sciences Publication Activity Database

Komm, Michael; Dejarnac, Renaud; Gunn, J. P.; Kirschner, A.; Litnovsky, A.; Matveev, D.; Pekarek, Z.

2011-01-01

Roč. 53, č. 11 (2011), s. 115004-115004 ISSN 0741-3335 R&D Projects: GA ČR GA202/09/1467 Institutional research plan: CEZ:AV0Z20430508 Keywords : Tokamak * gap * plasma PIC * retention power flux Subject RIV: BL - Plasma and Gas Discharge Physics Impact factor: 2.425, year: 2011 http://iopscience.iop.org/0741-3335/53/11/115004/pdf/0741-3335_53_11_115004.pdf

Study of plasma meniscus formation and beam halo in negative ion source using the 3D3VPIC model

International Nuclear Information System (INIS)

Nishioka, S.; Goto, I.; Hatayama, A.; Miyamoto, K.; Fukano, A.

2015-01-01

In this paper, the effect of the electron confinement time on the plasma meniscus and the fraction of the beam halo is investigated by 3D3V-PIC (three dimension in real space and three dimension in velocity space) (Particle in Cell) simulation in the extraction region of negative ion source. The electron confinement time depends on the characteristic time of electron escape along the magnetic field as well as the characteristic time of diffusion across the magnetic field. Our 3D3V-PIC results support the previous result by 2D3V-PIC results i.e., it is confirmed that the penetration of the plasma meniscus becomes deep into the source plasma region when the effective confinement time is short

Study of plasma meniscus formation and beam halo in negative ion source using the 3D3VPIC model

Energy Technology Data Exchange (ETDEWEB)

Nishioka, S.; Goto, I.; Hatayama, A. [Graduate school of Science and Technology, Keio University, Hiyoshi, Kohoku-ku, Yokohama 223-8522 (Japan); Miyamoto, K. [Naruto University of Education, 748 Nakashima, Takashima, Naruto-cho, Naruto-shi, Tokushima 772-8502 (Japan); Fukano, A. [Tokyo Metropolitan Collage of Industrial Technology, Higashioi, Shinagawa, Tokyo 140-0011 (Japan)

2015-04-08

In this paper, the effect of the electron confinement time on the plasma meniscus and the fraction of the beam halo is investigated by 3D3V-PIC (three dimension in real space and three dimension in velocity space) (Particle in Cell) simulation in the extraction region of negative ion source. The electron confinement time depends on the characteristic time of electron escape along the magnetic field as well as the characteristic time of diffusion across the magnetic field. Our 3D3V-PIC results support the previous result by 2D3V-PIC results i.e., it is confirmed that the penetration of the plasma meniscus becomes deep into the source plasma region when the effective confinement time is short.

Effect of plasma distribution on propulsion performance in electrodeless plasma thrusters

Science.gov (United States)

Takao, Yoshinori; Takase, Kazuki; Takahashi, Kazunori

2016-09-01

A helicon plasma thruster consisting of a helicon plasma source and a magnetic nozzle is one of the candidates for long-lifetime thrusters because no electrodes are employed to generate or accelerate plasma. A recent experiment, however, detected the non-negligible axial momentum lost to the lateral wall boundary, which degrades thruster performance, when the source was operated with highly ionized gases. To investigate this mechanism, we have conducted two-dimensional axisymmetric particle-in-cell (PIC) simulations with the neutral distribution obtained by Direct Simulation Monte Carlo (DSMC) method. The numerical results have indicated that the axially asymmetric profiles of the plasma density and potential are obtained when the strong decay of neutrals occurs at the source downstream. This asymmetric potential profile leads to the accelerated ion towards the lateral wall, leading to the non-negligible net axial force in the opposite direction of the thrust. Hence, to reduce this asymmetric profile by increasing the neutral density at downstream and/or by confining plasma with external magnetic field would result in improvement of the propulsion performance. These effects are also analyzed by PIC/DSMC simulations.

Electromagnetic particle-in-cell simulations of the solar wind interaction with lunar magnetic anomalies.

Science.gov (United States)

Deca, J; Divin, A; Lapenta, G; Lembège, B; Markidis, S; Horányi, M

2014-04-18

We present the first three-dimensional fully kinetic and electromagnetic simulations of the solar wind interaction with lunar crustal magnetic anomalies (LMAs). Using the implicit particle-in-cell code iPic3D, we confirm that LMAs may indeed be strong enough to stand off the solar wind from directly impacting the lunar surface forming a mini-magnetosphere, as suggested by spacecraft observations and theory. In contrast to earlier magnetohydrodynamics and hybrid simulations, the fully kinetic nature of iPic3D allows us to investigate the space charge effects and in particular the electron dynamics dominating the near-surface lunar plasma environment. We describe for the first time the interaction of a dipole model centered just below the lunar surface under plasma conditions such that only the electron population is magnetized. The fully kinetic treatment identifies electromagnetic modes that alter the magnetic field at scales determined by the electron physics. Driven by strong pressure anisotropies, the mini-magnetosphere is unstable over time, leading to only temporal shielding of the surface underneath. Future human exploration as well as lunar science in general therefore hinges on a better understanding of LMAs.

Particle-in-Cell laser-plasma simulation on Xeon Phi coprocessors

Science.gov (United States)

Surmin, I. A.; Bastrakov, S. I.; Efimenko, E. S.; Gonoskov, A. A.; Korzhimanov, A. V.; Meyerov, I. B.

2016-05-01

This paper concerns the development of a high-performance implementation of the Particle-in-Cell method for plasma simulation on Intel Xeon Phi coprocessors. We discuss the suitability of the method for Xeon Phi architecture and present our experience in the porting and optimization of the existing parallel Particle-in-Cell code PICADOR. Direct porting without code modification gives performance on Xeon Phi close to that of an 8-core CPU on a benchmark problem with 50 particles per cell. We demonstrate step-by-step optimization techniques, such as improving data locality, enhancing parallelization efficiency and vectorization leading to an overall 4.2 × speedup on CPU and 7.5 × on Xeon Phi compared to the baseline version. The optimized version achieves 16.9 ns per particle update on an Intel Xeon E5-2660 CPU and 9.3 ns per particle update on an Intel Xeon Phi 5110P. For a real problem of laser ion acceleration in targets with surface grating, where a large number of macroparticles per cell is required, the speedup of Xeon Phi compared to CPU is 1.6 ×.

Modeling an emittance-dominated elliptical sheet beam with a 212-dimensional particle-in-cell code

International Nuclear Information System (INIS)

Carlsten, Bruce E.

2005-01-01

Modeling a 3-dimensional (3-D) elliptical beam with a 212-D particle-in-cell (PIC) code requires a reduction in the beam parameters. The 212-D PIC code can only model the center slice of the sheet beam, but that can still provide useful information about the beam transport and distribution evolution, even if the beam is emittance dominated. The reduction of beam parameters and resulting interpretation of the simulation is straightforward, but not trivial. In this paper, we describe the beam parameter reduction and emittance issues related to the initial beam distribution. As a numerical example, we use the case of a sheet beam designed for use with a planar traveling-wave amplifier for high power generator for RF ranging from 95 to 300GHz [Carlsten et al., IEEE Trans. Plasma Sci. 33 (2005) 85]. These numerical techniques also apply to modeling high-energy elliptical bunches in RF accelerators

Low-noise electromagnetic δf particle-in-cell simulation of electron Bernstein waves

International Nuclear Information System (INIS)

Xiang Nong; Cary, John R.; Barnes, Daniel C.; Carlsson, John

2006-01-01

The conversion of the extraordinary (X) mode to an electron Bernstein wave (EBW) is one way to get rf energy into an overdense plasma. Analysis of this is complex, as the EBW is a fully kinetic wave, and so its linear propagation is described by an intractable integro-differential equation. Nonlinear effects cannot be calculated within this rubric at all. Full particle-in-cell (PIC) simulations cannot be used for these analyses, as the noise levels for reasonable simulation parameters are much greater than the typical rf amplitudes. It is shown that the delta-f computations are effective for this analysis. In particular, the accuracy of those computations has been verified by comparison with full PIC, cold plasma theory, and small gyroradius theory. This computational method is then used to analyze mode conversion in different frequency regimes. In particular, reasonable agreement with the theoretical predictions of Ram and Schultz [Phys. Plasmas 7, 4084 (2000)] in the linear regime is found, where 100% X-B mode conversion has been obtained when the driving frequency is less than twice the electron gyrofrequency. The results show that cold-plasma theory well predicts the mode conversion efficiency, as is consistent with the phase-space picture of mode conversion. From this it can be shown that nearly 100% X-B mode conversion cannot be obtained when the frequency is higher than the electron second harmonic cyclotron frequency

Particle-in-cell simulations of electron transport from plasma turbulence: recent progress in gyrokinetic particle simulations of turbulent plasmas

International Nuclear Information System (INIS)

Lin, Z; Rewoldt, G; Ethier, S; Hahm, T S; Lee, W W; Lewandowski, J L V; Nishimura, Y; Wang, W X

2005-01-01

Recent progress in gyrokinetic particle-in-cell simulations of turbulent plasmas using the gyrokinetic toroidal code (GTC) is surveyed. In particular, recent results for electron temperature gradient (ETG) modes and their resulting transport are presented. Also, turbulence spreading, and the effects of the parallel nonlinearity, are described. The GTC code has also been generalized for non-circular plasma cross-section, and initial results are presented. In addition, two distinct methods of generalizing the GTC code to be electromagnetic are described, along with preliminary results. Finally, a related code, GTC-Neo, for calculating neoclassical fluxes, electric fields, and velocities, are described

Breathing oscillations in enlarged cylindrical-anode-layer Hall plasma accelerator

Energy Technology Data Exchange (ETDEWEB)

Geng, S. F.; Wang, C. X. [Southwestern Institute of Physics, Chengdu 610041 (China); Department of Physics and Materials Science, City University of Hong Kong, Tat Chee Avenue, Kowloon, Hong Kong (China); Tang, D. L.; Qiu, X. M. [Southwestern Institute of Physics, Chengdu 610041 (China); Fu, R. K. Y. [Plasma Technology Limited, Festival Walk Tower, Tat Chee Avenue, Kowloon, Hong Kong (China); Chu, Paul K. [Department of Physics and Materials Science, City University of Hong Kong, Tat Chee Avenue, Kowloon, Hong Kong (China)

2013-05-28

Breathing oscillations in the discharge of an enlarged cylindrical-anode-layer Hall plasma accelerator are investigated by three-dimensional particle-in-cell (PIC) simulation. Different from the traditional breathing mode in a circular Hall plasma accelerator, the bulk plasma oscillation here is trigged by the potential barrier generated by the concentrated ion beam and substantial enough to compete with the anode voltage. The electric field near the anode is suppressed by the potential barrier thereby decreasing the electron density by {approx}36%. The discharge is restored to the normal level after the concentrated beam explodes and then it completes one cycle of electro-driven breathing oscillation. The breathing mode identified by the PIC simulation has a frequency range of {approx}156 kHz-{approx}250 kHz and does not vary monotonically with the discharge voltage.

Particle-in-cell simulations of plasma accelerators and electron-neutral collisions

Directory of Open Access Journals (Sweden)

David L. Bruhwiler

2001-10-01

Full Text Available We present 2D simulations of both beam-driven and laser-driven plasma wakefield accelerators, using the object-oriented particle-in-cell code XOOPIC, which is time explicit, fully electromagnetic, and capable of running on massively parallel supercomputers. Simulations of laser-driven wakefields with low \\(∼10^{16} W/cm^{2}\\ and high \\(∼10^{18} W/cm^{2}\\ peak intensity laser pulses are conducted in slab geometry, showing agreement with theory and fluid simulations. Simulations of the E-157 beam wakefield experiment at the Stanford Linear Accelerator Center, in which a 30 GeV electron beam passes through 1 m of preionized lithium plasma, are conducted in cylindrical geometry, obtaining good agreement with previous work. We briefly describe some of the more significant modifications to XOOPIC required by this work, and summarize the issues relevant to modeling relativistic electron-neutral collisions in a particle-in-cell code.

Compact disposal of high-energy electron beams using passive or laser-driven plasma decelerating stage

Energy Technology Data Exchange (ETDEWEB)

Bonatto, A.; Schroeder, C. B.; Vay, J. -L.; Geddes, C. R.; Benedetti, C.; Esarey and, E.; Leemans, W. P.

2014-07-13

A plasma decelerating stage is investigated as a compact alternative for the disposal of high-energy beams (beam dumps). This could benefit the design of laser-driven plasma accelerator (LPA) applications that require transportability and or high-repetition-rate operation regimes. Passive and laser-driven (active) plasma-based beam dumps are studied analytically and with particle-in-cell (PIC) simulations in a 1D geometry. Analytical estimates for the beam energy loss are compared to and extended by the PIC simulations, showing that with the proposed schemes a beam can be efficiently decelerated in a centimeter-scale distance.

Young gamma-ray pulsar: from modeling the gamma-ray emission to the particle-in-cell simulations of the global magnetosphere

Science.gov (United States)

Brambilla, Gabriele; Kalapotharakos, Constantions; Timokhin, Andrey; Kust Harding, Alice; Kazanas, Demosthenes

2016-04-01

Accelerated charged particles flowing in the magnetosphere produce pulsar gamma-ray emission. Pair creation processes produce an electron-positron plasma that populates the magnetosphere, in which the plasma is very close to force-free. However, it is unknown how and where the plasma departs from the ideal force-free condition, which consequently inhibits the understanding of the emission generation. We found that a dissipative magnetosphere outside the light cylinder effectively reproduces many aspects of the young gamma-ray pulsar emission as seen by the Fermi Gamma-ray Space Telescope, and through particle-in-cell simulations (PIC), we started explaining this configuration self-consistently. These findings show that, together, a magnetic field structure close to force-free and the assumption of gamma-ray curvature radiation as the emission mechanism are strongly compatible with the observations. Two main issues from the previously used models that our work addresses are the inability to explain luminosity, spectra, and light curve features at the same time and the inconsistency of the electrodynamics. Moreover, using the PIC simulations, we explore the effects of different pair multiplicities on the magnetosphere configurations and the locations of the accelerating regions. Our work aims for a self-consistent modeling of the magnetosphere, connecting the microphysics of the pair-plasma to the global magnetosphere macroscopic quantities. This direction will lead to a greater understanding of pulsar emission at all wavelengths, as well as to concrete insights into the physics of the magnetosphere.

Experimental and numerical study of electromagnetically induced transparency in magnetized plasmas

International Nuclear Information System (INIS)

Kawamori, Eiichirou; Hsieh, Tung-Yuan; Nishida, Yasushi; Cheng, C-Z

2012-01-01

We present a demonstration of electromagnetically induced transparency (EIT) in magnetized plasmas by means of experiment and numerical simulation. EIT in magnetized plasmas is a phenomenon by which a plasma-absorbing electron cyclotron wave is rendered transparent by a pump wave, which is a classical analog to conventional quantum EIT although the plasma EIT is not a quantum-mechanics-based phenomenon. This paper describes an attempt to identify plasma oscillations excited by the mode coupling of a pump wave and a probe wave, which is a key mechanism for achieving magnetized plasma EIT, by an experiment and a particle-in-cell (PIC) simulation. A preliminary result of the longitudinal electric field measurement indicates an enhancement of the plasma oscillation in the vicinity of the beat frequency between the probe and pump waves. Also the PIC calculation, which simulated the real experiment, shows a plasma oscillation excited by the mode coupling between the probe and pump waves in the magnetized plasma EIT, showing agreement with theory and experiment. (paper)

Plasma-particle interaction effects in induction plasma modelling under dense loading conditions

International Nuclear Information System (INIS)

Proulx, P.; Mostaghimi, J.; Boulos, M.

1983-07-01

The injection of solid particles or aerosol droplets in the fire-ball of an inductively coupled plasma can substantially perturb the plasma and even quench it under high loading conditions. This can be mainly attributed to the local cooling of the plasma by the particles or their vapour cloud, combined with the possible change of the thermodynamic and transport properties of the plasma in the presence of the particle vapour. This paper reports the state-of-the-art in the mathematical modelling of the induction plasma. A particle-in-cell model is used in order to combine the continuum approach for the calculation of the flow, temperature and concentration fields in the plasma, with the stochastic single particle approach, for the calculation of the particle trajectories and temperature histories. Results are given for an argon induction plasma under atmospheric pressure in which fine copper particles are centrally injected in the coil region of the discharge

Dust removal in radio-frequency plasmas by a traveling potential modulation

International Nuclear Information System (INIS)

Li Yangfang; Jiang Ke; Thomas, Hubertus M.; Morfill, Gregor E.

2010-01-01

The dust contamination in plasma deposition processes plays a crucial role in the quality and the yield of the products. To improve the quality and the yield of plasma processing, a favorable way is to remove the dust particles actively from the plasma reactors.Our recent experiments in the striped electrode device show that a traveling plasma modulation allows for a systematic particle removal independent of the reactor size. Besides the rf powered electrode, the striped electrode device includes a segmented electrode that consists of 100 electrically insulated narrow stripes. A traveling potential profile is produced by the modulation of the voltage signals applied on the stripes. The dust particles are trapped in the potential wells and transported with the traveling of the potential profile.The particle-in-cell (PIC) simulation on the potential above the segmented electrode indicates that the traveling potential profile can be realized either by applying low-frequency (0.1-10 Hz) voltage signals with a fixed phase shift between adjacent stripes or high-frequency (10 kHz a circumflex AS 100 MHz) signals with the amplitudes modulated by a low-frequency envelope. The transportation of the dust particles is simulated with a two-dimensional molecular dynamics (MD) code with the potential profile obtained from the PIC simulation. The MD results reproduce the experimental observations successfully.This technology allows for an active removal of the contaminating particles in processing plasmas and it is independent of the reactor size. The removal velocity is controllable by adjusting the parameters for the modulation.

A Particle-In-Cell approach to particle flux shaping with a surface mask

Directory of Open Access Journals (Sweden)

G. Kawamura

2017-08-01

Full Text Available The Particle-In-Cell simulation code PICS has been developed to study plasma in front of a surface with two types of masks, step-type and roof-type. Parameter scans with regard to magnetic field angle, electron density, and mask height were carried out to understand their influence on ion particle flux distribution on a surface. A roof-type mask with a small mask height yields short decay length in the flux distribution which is consistent with that estimated experimentally. A roof-type mask with a large height yields very long decay length and the flux value does not depend on a mask height or an electron density, but rather on a mask length and a biasing voltage of the surface. Mask height also changes the flux distribution apart from the mask because of the shading effect of the mask. Electron density changes the distribution near the mask edge according to the Debye length. Dependence of distribution on parameters are complicated especially for a roof-type mask, and simulation study with various parameters are useful to understand the physical reasons of dependence and also is useful as a tool for experiment studies.

Three-dimensional particle-in-cell simulations of gap crossings in castellated plasma-facing components in tokamaks

Czech Academy of Sciences Publication Activity Database

Komm, Michael; Dejarnac, Renaud; Gunn, J. P.; Pekarek, Z.

2013-01-01

Roč. 55, č. 2 (2013), 025006-025006 ISSN 0741-3335 R&D Projects: GA ČR GA202/09/1467; GA ČR GAP205/11/2341; GA MŠk(CZ) LM2011021 Institutional support: RVO:61389021 Keywords : plasma * tokamak * PIC divertor * castellation gaps Subject RIV: BL - Plasma and Gas Discharge Physics Impact factor: 2.386, year: 2013 http://iopscience.iop.org/0741-3335/55/2/025006/pdf/0741-3335_55_2_025006.pdf

An electrostatic particle-in-cell model for a lower hybrid grill

International Nuclear Information System (INIS)

Rantamaeki, K.

1998-01-01

In recent lower hybrid (LH) current drive experiments, generation of hot spots and impurities in the grill region have been observed on Tore Supra and Tokamak de Varennes (TdeV). A possible explanation is the parasitic absorption of the LH power in front of the grill. In parasitic absorption, the short-wavelength part of the lower hybrid spectrum can resonantly interact with the cold edge electrons. In this work, the absorption of the LH waves and the generation of fast electrons near the waveguide mouth is investigated with a new tool in this context: particle-in-cell (PIC) simulations. The advantage of this new method is that the electric field is calculated self-consistently. The PIC simulations also provide the key parameters for the hot spot problem: the absorbed power, the radial deposition profiles and the absorption length. A grill model has been added to the 2d3v PIC code XPDP2. Two sets of simulations were made. The first simulations used a phenomenological grill model. Strong absorption in the edge plasma was obtained. About 5% of the coupled power was absorbed within 1.7 mm in the case with fairly large amount of power in the modes with large parallel refractive index. Consequently, a rapid generation of fast electrons took place in the same region. In order to model experiments with realistic wave spectra, the PIC code was coupled to the slow wave antenna coupling code SWAN. The absorption within 1.7 mm in front of the grill was found to be between 2 and 5%. In the short time of a few wave periods, part of the initially thermal electrons (T e = 100 eV) were accelerated to velocities corresponding to a few keV. (orig.)

An electrostatic particle-in-cell model for a lower hybrid grill

Energy Technology Data Exchange (ETDEWEB)

Rantamaeki, K

1998-07-01

In recent lower hybrid (LH) current drive experiments, generation of hot spots and impurities in the grill region have been observed on Tore Supra and Tokamak de Varennes (TdeV). A possible explanation is the parasitic absorption of the LH power in front of the grill. In parasitic absorption, the short-wavelength part of the lower hybrid spectrum can resonantly interact with the cold edge electrons. In this work, the absorption of the LH waves and the generation of fast electrons near the waveguide mouth is investigated with a new tool in this context: particle-in-cell (PIC) simulations. The advantage of this new method is that the electric field is calculated self-consistently. The PIC simulations also provide the key parameters for the hot spot problem: the absorbed power, the radial deposition profiles and the absorption length. A grill model has been added to the 2d3v PIC code XPDP2. Two sets of simulations were made. The first simulations used a phenomenological grill model. Strong absorption in the edge plasma was obtained. About 5% of the coupled power was absorbed within 1.7 mm in the case with fairly large amount of power in the modes with large parallel refractive index. Consequently, a rapid generation of fast electrons took place in the same region. In order to model experiments with realistic wave spectra, the PIC code was coupled to the slow wave antenna coupling code SWAN. The absorption within 1.7 mm in front of the grill was found to be between 2 and 5%. In the short time of a few wave periods, part of the initially thermal electrons (T{sub e} = 100 eV) were accelerated to velocities corresponding to a few keV. (orig.)

ERO and PIC simulations of gross and net erosion of tungsten in the outer strike-point region of ASDEX Upgrade

Directory of Open Access Journals (Sweden)

A. Hakola

2017-08-01

Full Text Available We have modelled net and gross erosion of W in low-density l-mode plasmas in the low-field side strike point region of ASDEX Upgrade by ERO and Particle-in-Cell (PIC simulations. The observed net-erosion peak at the strike point was mainly due to the light impurities present in the plasma while the noticeable net-deposition regions surrounding the erosion maximum could be attributed to the strong E ×B drift and the magnetic field bringing eroded particles from a distance of several meters towards the private flux region. Our results also imply that the role of cross-field diffusion is very small in the studied plasmas. The simulations indicate net/gross erosion ratio of 0.2–0.6, which is in line with the literature data and what was determined spectroscopically. The deviations from the estimates extracted from post-exposure ion-beam-analysis data (∼0.6–0.7 are most likely due to the measured re-deposition patterns showing the outcomes of multiple erosion-deposition cycles.

On thermionic emission from plasma-facing components in tokamak-relevant conditions.

Czech Academy of Sciences Publication Activity Database

Komm, Michael; Ratynskaia, S.; Tolias, P.; Cavalier, Jordan; Dejarnac, Renaud; Gunn, J. P.; Podolník, Aleš

2017-01-01

Roč. 59, č. 9 (2017), č. článku 094002. ISSN 0741-3335 R&D Projects: GA ČR(CZ) GA16-14228S; GA MŠk(CZ) 8D15001 EU Projects: European Commission(XE) 633053 - EUROfusion Institutional support: RVO:61389021 Keywords : thermionic * PIC * tungsten * tokamak * thermionic emission * plasma facing components * particle-in-cell Subject RIV: BL - Plasma and Gas Discharge Physics OBOR OECD: Fluids and plasma physics (including surface physics) Impact factor: 2.392, year: 2016 http://iopscience.iop.org/article/10.1088/1361-6587/aa78c4/pdf

Fusion PIC code performance analysis on the Cori KNL system

Energy Technology Data Exchange (ETDEWEB)

Koskela, Tuomas S. [Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States). National Energy Research Scientific Computing Center (NERSC); Deslippe, Jack [Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States). National Energy Research Scientific Computing Center (NERSC); Friesen, Brian [Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States). National Energy Research Scientific Computing Center (NERSC); Raman, Karthic [INTEL Corp. (United States)

2017-05-25

We study the attainable performance of Particle-In-Cell codes on the Cori KNL system by analyzing a miniature particle push application based on the fusion PIC code XGC1. We start from the most basic building blocks of a PIC code and build up the complexity to identify the kernels that cost the most in performance and focus optimization efforts there. Particle push kernels operate at high AI and are not likely to be memory bandwidth or even cache bandwidth bound on KNL. Therefore, we see only minor benefits from the high bandwidth memory available on KNL, and achieving good vectorization is shown to be the most beneficial optimization path with theoretical yield of up to 8x speedup on KNL. In practice we are able to obtain up to a 4x gain from vectorization due to limitations set by the data layout and memory latency.

A curvilinear, fully implicit, conservative electromagnetic PIC algorithm in multiple dimensions

Science.gov (United States)

Chacón, L.; Chen, G.

2016-07-01

We extend a recently proposed fully implicit PIC algorithm for the Vlasov-Darwin model in multiple dimensions (Chen and Chacón (2015) [1]) to curvilinear geometry. As in the Cartesian case, the approach is based on a potential formulation (ϕ, A), and overcomes many difficulties of traditional semi-implicit Darwin PIC algorithms. Conservation theorems for local charge and global energy are derived in curvilinear representation, and then enforced discretely by a careful choice of the discretization of field and particle equations. Additionally, the algorithm conserves canonical-momentum in any ignorable direction, and preserves the Coulomb gauge ∇ ṡ A = 0 exactly. An asymptotically well-posed fluid preconditioner allows efficient use of large cell sizes, which are determined by accuracy considerations, not stability, and can be orders of magnitude larger than required in a standard explicit electromagnetic PIC simulation. We demonstrate the accuracy and efficiency properties of the algorithm with numerical experiments in mapped meshes in 1D-3V and 2D-3V.

The three-dimensional particle-in-cell simulation analysis of cavity of high power subterahertz pulsed gyrotron

International Nuclear Information System (INIS)

Ito, Koyu; Jiang, Weihua

2013-01-01

High power sub-terahertz pulsed gyrotrons for Collective Thomson Scattering (CTS) diagnostics of fusion plasmas are being developed. The typical target parameters are: output power of 100-200 kW, operation frequency of 300 GHz, and pulsed length > 10 us. In order to support experimental development, numerical simulations were carried out by using Particle-In-Cell (PIC) code MAGIC. The oscillation mode of the electromagnetic radiation was selected as TE_1_5_,_2, for which the beam parameters and cavity dimensions were determined accordingly. The simulation results have showed maximum power of 144 kW at oscillation frequency of 292.80 GHz, with oscillation efficiency of 22.15%. (author)

Charged particle acceleration in nonuniform plasmas

International Nuclear Information System (INIS)

Bulanov, S.V.; Naumova, N.M.; Pegoraro, F.

1996-11-01

The high-gradient electron acceleration schemes that have been demonstrated using LWFA appear promising for the development of plasma-based laser accelerators into practical devices. However, a question still exists: how to avoid the wake field deterioration and the loss of the phase synchronism between the plasma wave and the electrons that prevent them from being accelerated up to the theoretical limit. In order to obtain the highest possible values of the wake electric field one must use as intense laser pulses as possible i.e., pulses with dimensionless amplitudes a much-gt 1. Pulses that have a dimensionless amplitude larger than one tend to be subject to a host of instabilities, such as relativistic self-focusing, self modulation and stimulated Raman scattering, that affect their propagation in the plasma. Such processes could be beneficial, in so far as they increase the pulse energy density, enhance the wake field generation, and provide the mechanism for transporting the laser radiation over several Rayleigh lengths without diffraction spreading. However, it is still far from certain that these processes can be exploited in a controlled form and can lead to regular, stationary wake fields. It is known that, in order to create good quality wake fields, it would be preferable to use laser pulses with steep fronts of order λ p . The present paper aims at analyzing the influence of the laser pulse shape and of the plasma nonuniformity on the charged particle acceleration. This study is based on the results obtained with one dimensional PIC simulations

Database structure for plasma modeling programs

International Nuclear Information System (INIS)

Dufresne, M.; Silvester, P.P.

1993-01-01

Continuum plasma models often use a finite element (FE) formulation. Another approach is simulation models based on particle-in-cell (PIC) formulation. The model equations generally include four nonlinear differential equations specifying the plasma parameters. In simulation a large number of equations must be integrated iteratively to determine the plasma evolution from an initial state. The complexity of the resulting programs is a combination of the physics involved and the numerical method used. The data structure requirements of plasma programs are stated by defining suitable abstract data types. These abstractions are then reduced to data structures and a group of associated algorithms. These are implemented in an object oriented language (C++) as object classes. Base classes encapsulate data management into a group of common functions such as input-output management, instance variable updating and selection of objects by Boolean operations on their instance variables. Operations are thereby isolated from specific element types and uniformity of treatment is guaranteed. Creation of the data structures and associated functions for a particular plasma model is reduced merely to defining the finite element matrices for each equation, or the equations of motion for PIC models. Changes in numerical method or equation alterations are readily accommodated through the mechanism of inheritance, without modification of the data management software. The central data type is an n-relation implemented as a tuple of variable internal structure. Any finite element program may be described in terms of five relational tables: nodes, boundary conditions, sources, material/particle descriptions, and elements. Equivalently, plasma simulation programs may be described using four relational tables: cells, boundary conditions, sources, and particle descriptions

Simulation of ionization effects for high-density positron drivers in future plasma wakefield experiments

International Nuclear Information System (INIS)

Bruhwiler, D.L.; Dimitrov, D.A.; Cary, J.R.; Esarey, E.; Leemans, W.P.

2003-01-01

The plasma wakefield accelerator (PWFA) concept has been proposed as a potential energy doubler for present or future electron-positron colliders. Recent particle-in-cell (PIC) simulations have shown that the self-fields of the required electron beam driver can tunnel ionize neutral Li, leading to plasma wake dynamics differing significantly from that of a preionized plasma. It has also been shown, for the case of a preionized plasma, that the plasma wake of a positron driver differs strongly from that of an electron driver. We will present new PIC simulations, using the OOPIC code, showing the effects of tunneling ionization on the plasma wake generated by high-density positron drivers. The results will be compared to previous work on electron drivers with tunneling ionization and positron drivers without ionization. Parameters relevant to the energy doubler and the upcoming E-164x experiment at the Stanford Linear Accelerator Center will be considered

Activation of AMP-Activated Protein Kinase α and Extracelluar Signal-Regulated Kinase Mediates CB-PIC-Induced Apoptosis in Hypoxic SW620 Colorectal Cancer Cells

Directory of Open Access Journals (Sweden)

Sung-Yun Cho

2013-01-01

Full Text Available Here, antitumor mechanism of cinnamaldehyde derivative CB-PIC was elucidated in human SW620 colon cancer cells. CB-PIC significantly exerted cytotoxicity, increased sub-G1 accumulation, and cleaved PARP with apoptotic features, while it enhanced the phosphorylation of AMPK alpha and ACC as well as activated the ERK in hypoxic SW620 cells. Furthermore, CB-PIC suppressed the expression of HIF1 alpha, Akt, and mTOR and activated the AMPK phosphorylation in hypoxic SW620 cells. Conversely, silencing of AMPKα blocked PARP cleavage and ERK activation induced by CB-PIC, while ERK inhibitor PD 98059 attenuated the phosphorylation of AMPKα in hypoxic SW620 cells, implying cross-talk between ERK and AMPKα. Furthermore, cotreatment of CB-PIC and metformin enhanced the inhibition of HIF1α and Akt/mTOR and the activation of AMPKα and pACC in hypoxic SW620 cells. In addition, CB-PIC suppressed the growth of SW620 cells inoculated in BALB/c athymic nude mice, and immunohistochemistry revealed that CB-PIC treatment attenuated the expression of Ki-67, CD34, and CAIX and increased the expression of pAMPKα in CB-PIC-treated group. Interestingly, CP-PIC showed better antitumor activity in SW620 colon cancer cells under hypoxia than under normoxia, since it may be applied to chemoresistance. Overall, our findings suggest that activation of AMPKα and ERK mediates CB-PIC-induced apoptosis in hypoxic SW620 colon cancer cells.

Large Scale Earth's Bow Shock with Northern IMF as Simulated by PIC Code in Parallel with MHD Model

Science.gov (United States)

Baraka, Suleiman

2016-06-01

In this paper, we propose a 3D kinetic model (particle-in-cell, PIC) for the description of the large scale Earth's bow shock. The proposed version is stable and does not require huge or extensive computer resources. Because PIC simulations work with scaled plasma and field parameters, we also propose to validate our code by comparing its results with the available MHD simulations under same scaled solar wind (SW) and (IMF) conditions. We report new results from the two models. In both codes the Earth's bow shock position is found to be ≈14.8 R E along the Sun-Earth line, and ≈29 R E on the dusk side. Those findings are consistent with past in situ observations. Both simulations reproduce the theoretical jump conditions at the shock. However, the PIC code density and temperature distributions are inflated and slightly shifted sunward when compared to the MHD results. Kinetic electron motions and reflected ions upstream may cause this sunward shift. Species distributions in the foreshock region are depicted within the transition of the shock (measured ≈2 c/ ω pi for Θ Bn = 90° and M MS = 4.7) and in the downstream. The size of the foot jump in the magnetic field at the shock is measured to be (1.7 c/ ω pi ). In the foreshocked region, the thermal velocity is found equal to 213 km s-1 at 15 R E and is equal to 63 km s -1 at 12 R E (magnetosheath region). Despite the large cell size of the current version of the PIC code, it is powerful to retain macrostructure of planets magnetospheres in very short time, thus it can be used for pedagogical test purposes. It is also likely complementary with MHD to deepen our understanding of the large scale magnetosphere.

Raman laser amplification in preformed and ionizing plasmas

International Nuclear Information System (INIS)

Clark, D S; Fisch, N J

2004-01-01

The recently proposed backward Raman laser amplification scheme utilizes the stimulated Raman backscattering in plasma of a long pumping laser pulse to amplify a short, frequency downshifted seed pulse. The output intensity for this scheme is limited by the development of forward Raman scattering (FRS) or modulational instabilities of the highly amplified seed. Theoretically, focused output intensities as high as 1025 W/cm 2 and pulse lengths of less than 100 fs could be accessible by this technique for 1 (micro)m lasers--an improvement of 10 4 -10 5 in focused intensity over current techniques. Simulations with the particle-in-cell (PIC) code Zohar are presented which investigate the effects of FRS and modulational instabilities and of Langmuir wave breaking on the output intensity for Raman amplification. Using the intense seed pulse to photoionize the plasma simultaneous with its amplification (and hence avoid plasmas-based instabilities of the pump) is also investigated by PIC simulations. It is shown that both approaches can access focused intensities in the 1025 W/cm 2 range

Relationship Between Particle and Plasma Properties and Coating Characteristics of Samaria-Doped Ceria Prepared by Atmospheric Plasma Spraying for Use in Solid Oxide Fuel Cells

Science.gov (United States)

Cuglietta, Mark; Kesler, Olivera

2012-06-01

Samaria-doped ceria (SDC) has become a promising material for the fabrication of high-performance, intermediate-temperature solid oxide fuel cells (SOFCs). In this study, the in-flight characteristics, such as particle velocity and surface temperature, of spray-dried SDC agglomerates were measured and correlated to the resulting microstructures of SDC coatings fabricated using atmospheric plasma spraying, a manufacturing technique with the capability of producing full cells in minutes. Plasmas containing argon, nitrogen and hydrogen led to particle surface temperatures higher than those in plasmas containing only argon and nitrogen. A threshold temperature for the successful deposition of SDC on porous stainless steel substrates was calculated to be 2570 °C. Coating porosity was found to be linked to average particle temperature, suggesting that plasma conditions leading to lower particle temperatures may be most suitable for fabricating porous SOFC electrode layers.

Particle-in-cell plasma simulations of the modified two-stream instability

Directory of Open Access Journals (Sweden)

K. Schlegel

1994-08-01

Full Text Available We model the modified two-stream plasma instability occurring in the ionospheric E-region using a 2.5-dimensional particle-in-cell code. Compared to previous similar work we concentrate on simulated quantities that can easily be measured in the real ionosphere by coherent radars or rockets, such as the Doppler velocity, the backscattered power, backscattered spectra, aspect angle behaviour and electron temperature enhancement. Despite using a relatively small simulation model, we obtain remarkably good agreement between actual observed and simulated plasma parameters. The advantage of such a small system is that we were able to perform (other than in previous related work many simulation runs with different sets of input parameters, thus studying the unstable plasma under various conditions.

Enhanced stopping of macro-particles in particle-in-cell simulations

International Nuclear Information System (INIS)

May, J.; Tonge, J.; Ellis, I.; Mori, W. B.; Fiuza, F.; Fonseca, R. A.; Silva, L. O.; Ren, C.

2014-01-01

We derive an equation for energy transfer from relativistic charged particles to a cold background plasma appropriate for finite-size particles that are used in particle-in-cell simulation codes. Expressions for one-, two-, and three-dimensional particles are presented, with special attention given to the two-dimensional case. This energy transfer is due to the electric field of the wake set up in the background plasma by the relativistic particle. The enhanced stopping is dependent on the q 2 /m, where q is the charge and m is the mass of the relativistic particle, and therefore simulation macro-particles with large charge but identical q/m will stop more rapidly. The stopping power also depends on the effective particle shape of the macro-particle. These conclusions are verified in particle-in-cell simulations. We present 2D simulations of test particles, relaxation of high-energy tails, and integrated fast ignition simulations showing that the enhanced drag on macro-particles may adversely affect the results of these simulations in a wide range of high-energy density plasma scenarios. We also describe a particle splitting algorithm which can potentially overcome this problem and show its effect in controlling the stopping of macro-particles

Particle-in-cell Simulations of Raman Laser Amplification in Ionizing Plasmas

International Nuclear Information System (INIS)

Clark, Daniel S.; Fisch, Nathaniel J.

2003-01-01

By using the amplifying laser pulse in a plasma-based backward Raman laser amplifier to generate the plasma by photo-ionization of a gas simultaneous with the amplification process, possible instabilities of the pumping laser pulse can be avoided. Particle-in-cell simulations are used to study this amplification mechanism, and earlier results using more elementary models of the Raman interaction are verified [D.S. Clark and N.J. Fisch, Phys. Plasmas, 9 (6): 2772-2780, 2002]. The effects (unique to amplification in ionizing plasmas and not included in previous simulations) of blue-shifting of the pump and seed laser pulses and the generation of a wake are observed not significantly to impact the amplification process. As expected theoretically, the peak output intensity is found to be limited to I ∼ 10 17 W/cm 2 by forward Raman scattering of the amplifying seed. The integrity of the ionization front of the seed pulse against the development of a possible transverse modulation instability is also demonstrated

Resonant and non-resonant whistlers-particle interaction in the radiation belts

NARCIS (Netherlands)

E. Camporeale (Enrico)

2015-01-01

htmlabstractWe study the wave-particle interactions between lower band chorus whistlers and an anisotropic tenuous population of relativistic electrons. We present the first direct comparison of first-principle Particle-in-Cell (PIC) simulations with a quasi-linear diffusion code. In the PIC

Deployment of the OSIRIS EM-PIC code on the Intel Knights Landing architecture

Science.gov (United States)

Fonseca, Ricardo

2017-10-01

Electromagnetic particle-in-cell (EM-PIC) codes such as OSIRIS have found widespread use in modelling the highly nonlinear and kinetic processes that occur in several relevant plasma physics scenarios, ranging from astrophysical settings to high-intensity laser plasma interaction. Being computationally intensive, these codes require large scale HPC systems, and a continuous effort in adapting the algorithm to new hardware and computing paradigms. In this work, we report on our efforts on deploying the OSIRIS code on the new Intel Knights Landing (KNL) architecture. Unlike the previous generation (Knights Corner), these boards are standalone systems, and introduce several new features, include the new AVX-512 instructions and on-package MCDRAM. We will focus on the parallelization and vectorization strategies followed, as well as memory management, and present a detailed performance evaluation of code performance in comparison with the CPU code. This work was partially supported by Fundaçã para a Ciência e Tecnologia (FCT), Portugal, through Grant No. PTDC/FIS-PLA/2940/2014.

Simulation of damage to tokamaks plasma facing components during intense abnormal power deposition

International Nuclear Information System (INIS)

Genco, F.; Hassanein, A.

2014-01-01

Highlights: • HEIGHTS-PIC a new technique based on particle in cell method to study disruptions events, ELMS and VDE is benchmarked in this paper with the use of the MK-200 experiments. • Disruptions simulations results for erosion and erosion rate are proposed showing good agreement with published experimental available data for such conditions. • Results are also compared with other published results produced by FOREV1/FOREV2 computer package and the original HEIGHTS computer package. • Accuracy of the simulations results is proposed with specific aim to address the use of number of super particles adopted versus computational time. - Abstract: Intense power deposition on plasma facing components (PFC) is expected in tokamaks during loss of confinement events such as disruptions, vertical displacement events (VDE), runaway electrons (RE), or during normal operating conditions such as edge-localized modes (ELM). These highly energetic events are damaging enough to hinder long term operation and may not be easily mitigated without loss of structural or functional performance of the PFC. Surface erosion, melted/ablated-vaporized material splashing, and material transport into the bulk plasma are reliability-threatening for the machine and system performance. A novel particle-in-cell (PIC) technique has been developed and integrated into the existing HEIGHTS package in order to obtain a global view of the plasma evolution upon energy impingement. This newly developed PIC technique is benchmarked against plasma gun experimental data, the original HEIGHTS computer package, and laser experiments. Benchmarking results are shown in this paper for various relevant reactor and experimental devices. The evolution of the plasma vapor cloud is followed temporally and results are explained and commented as a function of the computational time needed and the accuracy of the calculation

Kinetic instabilities in plasmas: from electromagnetic fluctuations to collisionless shocks

International Nuclear Information System (INIS)

Ruyer, Charles

2014-01-01

Collisionless shocks play a major role in powerful astrophysical objects (e.g., gamma-ray bursts, supernova remnants, pulsar winds, etc.), where they are thought to be responsible for non-thermal particle acceleration and radiation. Numerical simulations have shown that, in the absence of an external magnetic field, these self-organizing structures originate from electromagnetic instabilities triggered by high-velocity colliding flows. These Weibel-like instabilities are indeed capable of producing the magnetic turbulence required for both efficient scattering and Fermi-type acceleration. Along with rapid advances in their theoretical understanding, intense effort is now underway to generate collisionless shocks in the laboratory using energetic lasers. In a first part we study the (w,k)-resolved electromagnetic thermal spectrum sustained by a drifting relativistic plasma. In particular, we obtain analytical formulae for the fluctuation spectra, the latter serving as seeds for growing magnetic modes in counterstreaming plasmas. Distinguishing between sub-luminal and supra-luminal thermal fluctuations, we derived analytical formulae of their respective spectral contributions. Comparisons with particle-in-cell (PIC) simulations are made, showing close agreement in the sub-luminal regime along with some discrepancy in the supra-luminal regime. Our formulae are then used to estimate the saturation time of the Weibel instability of relativistic pair plasmas. Our predictions are shown to match 2-D particle-in-cell (PIC) simulations over a three-decade range in flow energy. We then develop a predictive kinetic model of the nonlinear phase of the Weibel instability induced by two counter-streaming, symmetric and non-relativistic ion beams. This self consistent, fully analytical model allows us to follow the evolution of the beams' properties up to a stage close to complete isotropization and thus to shock formation. Its predictions are supported by 2D and 3D particle-in-cell

Fast 2D Fluid-Analytical Simulation of IEDs and Plasma Uniformity in Multi-frequency CCPs

Science.gov (United States)

Kawamura, E.; Lieberman, M. A.; Graves, D. B.

2014-10-01

A fast 2D axisymmetric fluid-analytical model using the finite elements tool COMSOL is interfaced with a 1D particle-in-cell (PIC) code to study ion energy distributions (IEDs) in multi-frequency argon capacitively coupled plasmas (CCPs). A bulk fluid plasma model which solves the time-dependent plasma fluid equations is coupled with an analytical sheath model which solves for the sheath parameters. The fluid-analytical results are used as input to a PIC simulation of the sheath region of the discharge to obtain the IEDs at the wafer electrode. Each fluid-analytical-PIC simulation on a moderate 2.2 GHz CPU workstation with 8 GB of memory took about 15-20 minutes. The 2D multi-frequency fluid-analytical model was compared to 1D PIC simulations of a symmetric parallel plate discharge, showing good agreement. Fluid-analytical simulations of a 2/60/162 MHz argon CCP with a typical asymmetric reactor geometry were also conducted. The low 2 MHz frequency controlled the sheath width and voltage while the higher frequencies controlled the plasma production. A standing wave was observable at the highest frequency of 162 MHz. Adding 2 MHz power to a 60 MHz discharge or 162 MHz to a dual frequency 2 MHz/60 MHz discharge enhanced the plasma uniformity. This work was supported by the Department of Energy Office of Fusion Energy Science Contract DE-SC000193, and in part by gifts from Lam Research Corporation and Micron Corporation.

Enhanced quasi-static particle-in-cell simulation of electron cloud instabilities in circular accelerators

Science.gov (United States)

Feng, Bing

Electron cloud instabilities have been observed in many circular accelerators around the world and raised concerns of future accelerators and possible upgrades. In this thesis, the electron cloud instabilities are studied with the quasi-static particle-in-cell (PIC) code QuickPIC. Modeling in three-dimensions the long timescale propagation of beam in electron clouds in circular accelerators requires faster and more efficient simulation codes. Thousands of processors are easily available for parallel computations. However, it is not straightforward to increase the effective speed of the simulation by running the same problem size on an increasingly number of processors because there is a limit to domain size in the decomposition of the two-dimensional part of the code. A pipelining algorithm applied on the fully parallelized particle-in-cell code QuickPIC is implemented to overcome this limit. The pipelining algorithm uses multiple groups of processors and optimizes the job allocation on the processors in parallel computing. With this novel algorithm, it is possible to use on the order of 102 processors, and to expand the scale and the speed of the simulation with QuickPIC by a similar factor. In addition to the efficiency improvement with the pipelining algorithm, the fidelity of QuickPIC is enhanced by adding two physics models, the beam space charge effect and the dispersion effect. Simulation of two specific circular machines is performed with the enhanced QuickPIC. First, the proposed upgrade to the Fermilab Main Injector is studied with an eye upon guiding the design of the upgrade and code validation. Moderate emittance growth is observed for the upgrade of increasing the bunch population by 5 times. But the simulation also shows that increasing the beam energy from 8GeV to 20GeV or above can effectively limit the emittance growth. Then the enhanced QuickPIC is used to simulate the electron cloud effect on electron beam in the Cornell Energy Recovery Linac

Parallel Finite Element Particle-In-Cell Code for Simulations of Space-charge Dominated Beam-Cavity Interactions

International Nuclear Information System (INIS)

Candel, A.; Kabel, A.; Ko, K.; Lee, L.; Li, Z.; Limborg, C.; Ng, C.; Prudencio, E.; Schussman, G.; Uplenchwar, R.

2007-01-01

Over the past years, SLAC's Advanced Computations Department (ACD) has developed the parallel finite element (FE) particle-in-cell code Pic3P (Pic2P) for simulations of beam-cavity interactions dominated by space-charge effects. As opposed to standard space-charge dominated beam transport codes, which are based on the electrostatic approximation, Pic3P (Pic2P) includes space-charge, retardation and boundary effects as it self-consistently solves the complete set of Maxwell-Lorentz equations using higher-order FE methods on conformal meshes. Use of efficient, large-scale parallel processing allows for the modeling of photoinjectors with unprecedented accuracy, aiding the design and operation of the next-generation of accelerator facilities. Applications to the Linac Coherent Light Source (LCLS) RF gun are presented

Self-consistent kinetic simulations of lower hybrid drift instability resulting in electron current driven by fusion products in tokamak plasmas

International Nuclear Information System (INIS)

Cook, J W S; Chapman, S C; Dendy, R O; Brady, C S

2011-01-01

We present particle-in-cell (PIC) simulations of minority energetic protons in deuterium plasmas, which demonstrate a collective instability responsible for emission near the lower hybrid frequency and its harmonics. The simulations capture the lower hybrid drift instability in a parameter regime motivated by tokamak fusion plasma conditions, and show further that the excited electromagnetic fields collectively and collisionlessly couple free energy from the protons to directed electron motion. This results in an asymmetric tail antiparallel to the magnetic field. We focus on obliquely propagating modes excited by energetic ions, whose ring-beam distribution is motivated by population inversions related to ion cyclotron emission, in a background plasma with a temperature similar to that of the core of a large tokamak plasma. A fully self-consistent electromagnetic relativistic PIC code representing all vector field quantities and particle velocities in three dimensions as functions of a single spatial dimension is used to model this situation, by evolving the initial antiparallel travelling ring-beam distribution of 3 MeV protons in a background 10 keV Maxwellian deuterium plasma with realistic ion-electron mass ratio. These simulations provide a proof-of-principle for a key plasma physics process that may be exploited in future alpha channelling scenarios for magnetically confined burning plasmas.

Controlling the numerical Cerenkov instability in PIC simulations using a customized finite difference Maxwell solver and a local FFT based current correction

International Nuclear Information System (INIS)

Li, Fei; Yu, Peicheng; Xu, Xinlu; Fiuza, Frederico; Decyk, Viktor K.

2017-01-01

In this study we present a customized finite-difference-time-domain (FDTD) Maxwell solver for the particle-in-cell (PIC) algorithm. The solver is customized to effectively eliminate the numerical Cerenkov instability (NCI) which arises when a plasma (neutral or non-neutral) relativistically drifts on a grid when using the PIC algorithm. We control the EM dispersion curve in the direction of the plasma drift of a FDTD Maxwell solver by using a customized higher order finite difference operator for the spatial derivative along the direction of the drift (1^ direction). We show that this eliminates the main NCI modes with moderate |k_1|, while keeps additional main NCI modes well outside the range of physical interest with higher |k_1|. These main NCI modes can be easily filtered out along with first spatial aliasing NCI modes which are also at the edge of the fundamental Brillouin zone. The customized solver has the possible advantage of improved parallel scalability because it can be easily partitioned along 1^ which typically has many more cells than other directions for the problems of interest. We show that FFTs can be performed locally to current on each partition to filter out the main and first spatial aliasing NCI modes, and to correct the current so that it satisfies the continuity equation for the customized spatial derivative. This ensures that Gauss’ Law is satisfied. Lastly, we present simulation examples of one relativistically drifting plasma, of two colliding relativistically drifting plasmas, and of nonlinear laser wakefield acceleration (LWFA) in a Lorentz boosted frame that show no evidence of the NCI can be observed when using this customized Maxwell solver together with its NCI elimination scheme.

Time-Domain Modeling of RF Antennas and Plasma-Surface Interactions

Directory of Open Access Journals (Sweden)

Jenkins Thomas G.

2017-01-01

Full Text Available Recent advances in finite-difference time-domain (FDTD modeling techniques allow plasma-surface interactions such as sheath formation and sputtering to be modeled concurrently with the physics of antenna near- and far-field behavior and ICRF power flow. Although typical sheath length scales (micrometers are much smaller than the wavelengths of fast (tens of cm and slow (millimeter waves excited by the antenna, sheath behavior near plasma-facing antenna components can be represented by a sub-grid kinetic sheath boundary condition, from which RF-rectified sheath potential variation over the surface is computed as a function of current flow and local plasma parameters near the wall. These local time-varying sheath potentials can then be used, in tandem with particle-in-cell (PIC models of the edge plasma, to study sputtering effects. Particle strike energies at the wall can be computed more accurately, consistent with their passage through the known potential of the sheath, such that correspondingly increased accuracy of sputtering yields and heat/particle fluxes to antenna surfaces is obtained. The new simulation capabilities enable time-domain modeling of plasma-surface interactions and ICRF physics in realistic experimental configurations at unprecedented spatial resolution. We will present results/animations from high-performance (10k-100k core FDTD/PIC simulations of Alcator C-Mod antenna operation.

Wavelet-based Poisson Solver for use in Particle-In-Cell Simulations

International Nuclear Information System (INIS)

Terzic, B.; Mihalcea, D.; Bohn, C.L.; Pogorelov, I.V.

2005-01-01

We report on a successful implementation of a wavelet based Poisson solver for use in 3D particle-in-cell (PIC) simulations. One new aspect of our algorithm is its ability to treat the general(inhomogeneous) Dirichlet boundary conditions (BCs). The solver harnesses advantages afforded by the wavelet formulation, such as sparsity of operators and data sets, existence of effective preconditioners, and the ability simultaneously to remove numerical noise and further compress relevant data sets. Having tested our method as a stand-alone solver on two model problems, we merged it into IMPACT-T to obtain a fully functional serial PIC code. We present and discuss preliminary results of application of the new code to the modeling of the Fermilab/NICADD and AES/JLab photoinjectors

Three-dimensional electromagnetic solitary waves in an underdense plasma in PIC simulations

International Nuclear Information System (INIS)

Bulanov, Sergei; Esirkepov, Timur; Nishihara, Katsunobu; Pegoraro, Francesco

2002-01-01

A three-dimensional sub-cycle relativistic electromagnetic soliton has been observed for the first time in a 3D Particle-in-Cell simulation of the propagation of an intense short laser pulse in an underdense plasma. The structure of the 3D soliton is identified. It resembles an oscillating electric dipole and has a strong charge separation and toroidal magnetic field component. We call this structure a TM-soliton (transverse magnetic). The 3D TM-soliton resembles a 2D P-soliton in the plane of electric field polarization, and a 2D S-soliton in the perpendicular plane. The core of the soliton is positively charged on average in time, and this results in its Coulomb explosion and in ion heating. Then the soliton evolves into a post-soliton, which is a slowly expanding quasi-neutral cavity in the plasma

Fast ignition: Dependence of the ignition energy on source and target parameters for particle-in-cell-modelled energy and angular distributions of the fast electrons

Energy Technology Data Exchange (ETDEWEB)

Bellei, C.; Divol, L.; Kemp, A. J.; Key, M. H.; Larson, D. J.; Strozzi, D. J.; Marinak, M. M.; Tabak, M.; Patel, P. K. [Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, California 94550 (United States)

2013-05-15

The energy and angular distributions of the fast electrons predicted by particle-in-cell (PIC) simulations differ from those historically assumed in ignition designs of the fast ignition scheme. Using a particular 3D PIC calculation, we show how the ignition energy varies as a function of source-fuel distance, source size, and density of the pre-compressed fuel. The large divergence of the electron beam implies that the ignition energy scales with density more weakly than the ρ{sup −2} scaling for an idealized beam [S. Atzeni, Phys. Plasmas 6, 3316 (1999)], for any realistic source that is at some distance from the dense deuterium-tritium fuel. Due to the strong dependence of ignition energy with source-fuel distance, the use of magnetic or electric fields seems essential for the purpose of decreasing the ignition energy.

Theoretical and numerical studies on the transport of transverse beam quality in plasma-based accelerators

International Nuclear Information System (INIS)

Mehrling, Timon Johannes

2014-11-01

This work examines effects, which impact the transverse quality of electron-beams in plasma-based accelerators, by means of theoretical and numerical methods. Plasma-based acceleration is a promising candidate for future particle accelerator technologies. In plasma-based acceleration, highly intense laser beams or high-current relativistic particle beams are focused into a plasma to excite plasma-waves with extreme transverse and longitudinal electric fields. The amplitude of these fields exceed with 10-100 GV/m the ones in today's radio-frequency accelerators by several orders of magnitude, hence, in principle allowing for accordingly shorter and cheaper accelerators based on plasma. Despite the tremendous progress in the recent decade, beams from plasma accelerators are not yet achieving the quality as demanded for pivotal applications of relativistic electron-beams, e.g. free-electron lasers (FELs).Studies within this work examine how the quality can be optimized in the production of the beams and preserved during the acceleration and transport to the interaction region. Such studies cannot be approached purely analytical but necessitate numerical methods, such as the Particle-In-Cell (PIC) method, which can model kinetic, electrodynamic and relativistic plasma phenomena. However, this method is computationally too expensive for parameter-scans in three-dimensional geometries. Hence, a quasi-static PIC code was developed in connection with this work, which is significantly more effective than the full PIC method for a class of problems in plasma-based acceleration.The evolution of the emittance of beams which are injected into plasma modules was studied in this work by means of theoretical and the above numerical methods. It was shown that the beam parameters need to be matched accurately into the focusing plasma-channel in order to allow for beam-quality preservation. This suggested that new extraction and injection-techniques are required in staged plasma

Particle-in-cell modeling of laser Thomson scattering in low-density plasmas at elevated laser intensities

Science.gov (United States)

Powis, Andrew T.; Shneider, Mikhail N.

2018-05-01

Incoherent Thomson scattering is a non-intrusive technique commonly used for measuring local plasma density. Within low-density, low-temperature plasmas and for sufficient laser intensity, the laser may perturb the local electron density via the ponderomotive force, causing the diagnostic to become intrusive and leading to erroneous results. A theoretical model for this effect is validated numerically via kinetic simulations of a quasi-neutral plasma using the particle-in-cell technique.

A multidimensional theory for electron trapping by a plasma wake generated in the bubble regime

International Nuclear Information System (INIS)

Kostyukov, I; Nerush, E; Pukhov, A; Seredov, V

2010-01-01

We present a theory for electron self-injection in nonlinear, multidimensional plasma waves excited by a short laser pulse in the bubble regime or by a short electron beam in the blowout regime. In these regimes, which are typical for electron acceleration in the last impressive experiments, the laser radiation pressure or the electron beam charge pushes out plasma electrons from some region, forming a plasma cavity or a bubble with a huge ion charge. The plasma electrons can be trapped in the bubble and accelerated by the plasma wakefields up to a very high energy. We derive the condition of the electron trapping in the bubble. The developed theory predicts the trapping cross section in terms of the bubble radius and the bubble velocity. It is found that the dynamic bubble deformations observed in the three-dimensional (3D) particle-in-cell (PIC) simulations influence the trapping process significantly. The bubble elongation reduces the gamma-factor of the bubble, thereby strongly enhancing self-injection. The obtained analytical results are in good agreement with the 3D PIC simulations.

Characterizing Hypervelocity Impact Plasma Through Experiments and Simulations

Science.gov (United States)

Close, Sigrid; Lee, Nicolas; Fletcher, Alex; Nuttall, Andrew; Hew, Monica; Tarantino, Paul

2017-10-01

Hypervelocity micro particles, including meteoroids and space debris with masses produce a strong electromagnetic pulse (EMP) with a broad frequency spectrum. Subsequent plasma oscillations resulting from instabilities can also emit significant power and may be responsible for many reported satellite anomalies. We present theory and recent results from ground-based impact tests aimed at characterizing hypervelocity impact plasma. We also show results from particle-in-cell (PIC) and computational fluid dynamics (CFD) simulations that allow us to extend to regimes not currently possible with ground-based technology. We show that significant impact-produced radio frequency (RF) emissions occurred in frequencies ranging from VHF through L-band and that these emissions were highly correlated with fast (>20 km/s) impacts that produced a fully ionized plasma.

Stochastic plasma heating by electrostatic waves: a comparison between a particle-in-cell simulation and a laboratory experiment

International Nuclear Information System (INIS)

Fivaz, M.; Fasoli, A.; Appert, K.; Trans, T.M.; Tran, M.Q.; Skiff, F.

1993-08-01

Dynamical chaos is produced by the interaction between plasma particles and two electrostatic waves. Experiments performed in a linear magnetized plasma and a 1D particle-in-cell simulation agree qualitatively: above a threshold wave amplitude, ion stochastic diffusion and heating occur on a fast time scale. Self-consistency appears to limit the extent of the heating process. (author) 5 figs., 18 refs

Particle-in-cell vs straight line Gaussian calculations for an area of complex topography

International Nuclear Information System (INIS)

Lange, R.; Sherman, C.

1977-01-01

Two numerical models for the calculation of time integrated air concentraton and ground deposition of airborne radioactive effluent releases are compared. The time dependent Particle-in-Cell (PIC) model and the steady state Gaussian plume model were used for the simulation. The area selected for the comparison was the Hudson River Valley, New York. Input for the models was synthesized from meteorological data gathered in previous studies by various investigators. It was found that the PIC model more closely simulated the three-dimensional effects of the meteorology and topography. Overall, the Gaussian model calculated higher concentrations under stable conditions. In addition, because of its consideration of exposure from the returning plume after flow reversal, the PIC model calculated air concentrations over larger areas than did the Gaussian model

Global fully kinetic models of planetary magnetospheres with iPic3D

Science.gov (United States)

Gonzalez, D.; Sanna, L.; Amaya, J.; Zitz, A.; Lembege, B.; Markidis, S.; Schriver, D.; Walker, R. J.; Berchem, J.; Peng, I. B.; Travnicek, P. M.; Lapenta, G.

2016-12-01

We report on the latest developments of our approach to model planetary magnetospheres, mini magnetospheres and the Earth's magnetosphere with the fully kinetic, electromagnetic particle in cell code iPic3D. The code treats electrons and multiple species of ions as full kinetic particles. We review: 1) Why a fully kinetic model and in particular why kinetic electrons are needed for capturing some of the most important aspects of the physics processes of planetary magnetospheres. 2) Why the energy conserving implicit method (ECIM) in its newest implementation [1] is the right approach to reach this goal. We consider the different electron scales and study how the new IECIM can be tuned to resolve only the electron scales of interest while averaging over the unresolved scales preserving their contribution to the evolution. 3) How with modern computing planetary magnetospheres, mini magnetosphere and eventually Earth's magnetosphere can be modeled with fully kinetic electrons. The path from petascale to exascale for iPiC3D is outlined based on the DEEP-ER project [2], using dynamic allocation of different processor architectures (Xeon and Xeon Phi) and innovative I/O technologies.Specifically results from models of Mercury are presented and compared with MESSENGER observations and with previous hybrid (fluid electrons and kinetic ions) simulations. The plasma convection around the planets includes the development of hydrodynamic instabilities at the flanks, the presence of the collisionless shocks, the magnetosheath, the magnetopause, reconnection zones, the formation of the plasma sheet and the magnetotail, and the variation of ion/electron plasma flows when crossing these frontiers. Given the full kinetic nature of our approach we focus on detailed particle dynamics and distribution at locations that can be used for comparison with satellite data. [1] Lapenta, G. (2016). Exactly Energy Conserving Implicit Moment Particle in Cell Formulation. arXiv preprint ar

Novel methods in the Particle-In-Cell accelerator Code-Framework Warp

Energy Technology Data Exchange (ETDEWEB)

Vay, J-L [Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States); Grote, D. P. [Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States); Cohen, R. H. [Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States); Friedman, A. [Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)

2012-12-26

The Particle-In-Cell (PIC) Code-Framework Warp is being developed by the Heavy Ion Fusion Science Virtual National Laboratory (HIFS-VNL) to guide the development of accelerators that can deliver beams suitable for high-energy density experiments and implosion of inertial fusion capsules. It is also applied in various areas outside the Heavy Ion Fusion program to the study and design of existing and next-generation high-energy accelerators, including the study of electron cloud effects and laser wakefield acceleration for example. This study presents an overview of Warp's capabilities, summarizing recent original numerical methods that were developed by the HIFS-VNL (including PIC with adaptive mesh refinement, a large-timestep 'drift-Lorentz' mover for arbitrarily magnetized species, a relativistic Lorentz invariant leapfrog particle pusher, simulations in Lorentz-boosted frames, an electromagnetic solver with tunable numerical dispersion and efficient stride-based digital filtering), with special emphasis on the description of the mesh refinement capability. In addition, selected examples of the applications of the methods to the abovementioned fields are given.

Particle-in-cell Simulations of Raman Laser Amplification in Preformed Plasmas

International Nuclear Information System (INIS)

Clark, Daniel S.; Fisch, Nathaniel J.

2003-01-01

Two critical issues in the amplification of laser pulses by backward Raman scattering in plasma slabs are the saturation mechanism of the amplification effect (which determines the maximum attainable output intensity of a Raman amplifier) and the optimal plasma density for amplification. Previous investigations [V.M. Malkin, et al., Phys. Rev. Lett., 82 (22):4448-4451, 1999] identified forward Raman scattering and modulational instabilities of the amplifying seed as the likely saturation mechanisms and lead to an estimated unfocused output intensities of 10 17 W/cm 2 . The optimal density for amplification is determined by the competing constraints of minimizing the plasma density so as to minimize the growth rate of the instabilities leading to saturation but also maintaining the plasma sufficiently dense that the driven Langmuir wave responsible for backscattering does not break prematurely. Here, particle-in-cell code are simulations presented which verify that saturation of backward Raman amplification does occur at intensities of ∼10 17 W/cm 2 by forward Raman scattering and modulational instabilities. The optimal density for amplification in a plasma with the representative temperature of T(sub)e = 200 eV is also shown in these simulations to be intermediate between the cold plasma wave-breaking density and the density limit found by assuming a water bag electron distribution function

PIC Simulation of Laser Plasma Interactions with Temporal Bandwidths

Science.gov (United States)

Tsung, Frank; Weaver, J.; Lehmberg, R.

2015-11-01

We are performing particle-in-cell simulations using the code OSIRIS to study the effects of laser plasma interactions in the presence of temperal bandwidths under conditions relevant to current and future shock ignition experiments on the NIKE laser. Our simulations show that, for sufficiently large bandwidth, the saturation level, and the distribution of hot electrons, can be effected by the addition of temporal bandwidths (which can be accomplished in experiments using smoothing techniques such as SSD or ISI). We will show that temporal bandwidth along play an important role in the control of LPI's in these lasers and discuss future directions. This work is conducted under the auspices of NRL.

Dynamics of electron bunches at the laser–plasma interaction in the bubble regime

Energy Technology Data Exchange (ETDEWEB)

Maslov, V.I., E-mail: vmaslov@kipt.kharkov.ua; Svystun, O.M., E-mail: svistun_elena@mail.ru; Onishchenko, I.N.; Tkachenko, V.I.

2016-09-01

The multi-bunches self-injection, observed in laser–plasma accelerators in the bubble regime, affects the energy gain of electrons accelerated by laser wakefield. However, understanding of dynamics of the electron bunches formed at laser–plasma interaction may be challenging. We present here the results of fully relativistic electromagnetic particle-in-cell (PIC) simulation of laser wakefield acceleration driven by a short laser pulse in an underdense plasma. The trapping and acceleration of three witness electron bunches by the bubble-like structures were observed. It has been shown that with time the first two witness bunches turn into drivers and contribute to acceleration of the last witness bunch.

3-D electromagnetic plasma particle simulations on the Intel Delta parallel computer

International Nuclear Information System (INIS)

Wang, J.; Liewer, P.C.

1994-01-01

A three-dimensional electromagnetic PIC code has been developed on the 512 node Intel Touchstone Delta MIMD parallel computer. This code is based on the General Concurrent PIC algorithm which uses a domain decomposition to divide the computation among the processors. The 3D simulation domain can be partitioned into 1-, 2-, or 3-dimensional sub-domains. Particles must be exchanged between processors as they move among the subdomains. The Intel Delta allows one to use this code for very-large-scale simulations (i.e. over 10 8 particles and 10 6 grid cells). The parallel efficiency of this code is measured, and the overall code performance on the Delta is compared with that on Cray supercomputers. It is shown that their code runs with a high parallel efficiency of ≥ 95% for large size problems. The particle push time achieved is 115 nsecs/particle/time step for 162 million particles on 512 nodes. Comparing with the performance on a single processor Cray C90, this represents a factor of 58 speedup. The code uses a finite-difference leap frog method for field solve which is significantly more efficient than fast fourier transforms on parallel computers. The performance of this code on the 128 node Cray T3D will also be discussed

Modeling of Perpendicularly Driven Dual-Frequency Capacitively Coupled Plasma

International Nuclear Information System (INIS)

Wang Hongyu; Sun Peng; Zhao Shuangyun; Li Yang; Jiang Wei

2016-01-01

We analyzed perpendicularly configured dual-frequency (DF) capacitively coupled plasmas (CCP). In this configuration, two pairs of electrodes are arranged oppositely, and the discharging is perpendicularly driven by two radio frequency (RF) sources. Particle-in-cell/Monte Carlo (PIC/MC) simulation showed that the configuration had some advantages as this configuration eliminated some dual frequency coupling effects. Some variation and potential application of the discharging configuration is discussed briefly. (paper)

Development of 2D particle-in-cell code to simulate high current, low ...

Indian Academy of Sciences (India)

Abstract. A code for 2D space-charge dominated beam dynamics study in beam trans- port lines is developed. The code is used for particle-in-cell (PIC) simulation of z-uniform beam in a channel containing solenoids and drift space. It can also simulate a transport line where quadrupoles are used for focusing the beam.

Feasibility study of the plasma electron density measurement by electromagnetic radiation from the laser-driven plasma wave

International Nuclear Information System (INIS)

Jang, D G; Kim, J J; Suk, H; Hur, M S

2012-01-01

When an intense laser beam is focused in a plasma, a plasma wake wave is generated and the oscillatary motion of the plasma electrons produces a strong electromagnetic wave by a Cherenkov-like process. Spectrum of the genetated electromagnetic wave has dependence on the plasma density. In this paper, we propose to use the emitted electromagnetic radiation for plasma diagnostic, which may provide an accurate information for local electron densities of the plasma and will be very useful for three-dimensional plasma density profiles by changing the focal point location of the laser beam. Two-dimensional (2-D) particle-in-cell (PIC) simulation is used to study the correlation between the spectrum of the emitted radiation and plasma density, and the results demonstrate that this method is promising for the electron density measurement in the plasma.

Collimated fast electron beam generation in critical density plasma

Energy Technology Data Exchange (ETDEWEB)

Iwawaki, T., E-mail: iwawaki-t@eie.eng.osaka-u.ac.jp; Habara, H.; Morita, K.; Tanaka, K. A. [Graduate School of Engineering, Osaka University, 2-1, Yamada-oka, Suita, Osaka 565-0871 (Japan); Baton, S.; Fuchs, J.; Chen, S. [LULI, CNRS-Ecole Polytechnique-Université Pierre et Marie Curie-CEA, 91128 Palaiseau (France); Nakatsutsumi, M. [LULI, CNRS-Ecole Polytechnique-Université Pierre et Marie Curie-CEA, 91128 Palaiseau (France); European X-Ray Free-Electron Laser Facility (XFEL) GmbH (Germany); Rousseaux, C. [CEA, DAM, DIF, F-91297 Arpajon (France); Filippi, F. [La SAPIENZA, University of Rome, Dip. SBAI, 00161 Rome (Italy); Nazarov, W. [School of Chemistry, University of St. Andrews, North Haugh, St. Andrews, Fife KY16 9ST, Scotland (United Kingdom)

2014-11-15

Significantly collimated fast electron beam with a divergence angle 10° (FWHM) is observed when an ultra-intense laser pulse (I = 10{sup 14 }W/cm{sup 2}, 300 fs) irradiates a uniform critical density plasma. The uniform plasma is created through the ionization of an ultra-low density (5 mg/c.c.) plastic foam by X-ray burst from the interaction of intense laser (I = 10{sup 14 }W/cm{sup 2}, 600 ps) with a thin Cu foil. 2D Particle-In-Cell (PIC) simulation well reproduces the collimated electron beam with a strong magnetic field in the region of the laser pulse propagation. To understand the physical mechanism of the collimation, we calculate energetic electron motion in the magnetic field obtained from the 2D PIC simulation. As the results, the strong magnetic field (300 MG) collimates electrons with energy over a few MeV. This collimation mechanism may attract attention in many applications such as electron acceleration, electron microscope and fast ignition of laser fusion.

Memory-efficient optimization of Gyrokinetic particle-to-grid interpolation for multicore processors

Energy Technology Data Exchange (ETDEWEB)

Madduri, Kamesh [Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States); Williams, Samuel [Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States); Ethier, Stephane [Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States); Oliker, Leonid [Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States); Shalf, John [Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States); Strohmaier, Erich [Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States); Yelicky, Katherine [Univ. of California, Berkeley, CA (United States)

2009-01-01

We present multicore parallelization strategies for the particle-to-grid interpolation step in the Gyrokinetic Toroidal Code (GTC), a 3D particle-in-cell (PIC) application to study turbulent transport in magnetic-confinement fusion devices. Particle-grid interpolation is a known performance bottleneck in several PIC applications. In GTC, this step involves particles depositing charges to a 3D toroidal mesh, and multiple particles may contribute to the charge at a grid point. We design new parallel algorithms for the GTC charge deposition kernel, and analyze their performance on three leading multicore platforms. We implement thirteen different variants for this kernel and identify the best-performing ones given typical PIC parameters such as the grid size, number of particles per cell, and the GTC-specific particle Larmor radius variation. We find that our best strategies can be 2x faster than the reference optimized MPI implementation, and our analysis provides insight into desirable architectural features for high-performance PIC simulation codes.

Intense EM filamentation in relativistic hot plasmas

Energy Technology Data Exchange (ETDEWEB)

Hu, Qiang-Lin [Department of Physics, Jinggangshan University, Ji' an, Jiangxi 343009 (China); Chen, Zhong-Ping [Department of Physics and Institute for Fusion Studies, The University of Texas at Austin, Austin, TX 78712 (United States); Mahajan, Swadesh M., E-mail: mahajan@mail.utexas.edu [Department of Physics and Institute for Fusion Studies, The University of Texas at Austin, Austin, TX 78712 (United States); Department of Physics, School of Natural Sciences, Shiv Nadar University, Uttar Pradesh 201314 (India)

2017-03-03

Highlights: • Breaking up of an intense EM pulse into filaments is a spectacular demonstration of the nonlinear wave-plasma interaction. • Filaments are spectacularly sharper, highly extended and longer lived at relativistic temperatures. • EM energy concentration can trigger new nonlinear phenomena with absolute consequences for high energy density matter. - Abstract: Through 2D particle-in-cell (PIC) simulations, we demonstrate that the nature of filamentation of a high intensity electromagnetic (EM) pulse propagating in an underdense plasma, is profoundly affected at relativistically high temperatures. The “relativistic” filaments are sharper, are dramatically extended (along the direction of propagation), and live much longer than their lower temperature counterparts. The thermally boosted electron inertia is invoked to understand this very interesting and powerful phenomenon.

PIC Simulations of Velocity-space Instabilities in a Decreasing Magnetic Field: Viscosity and Thermal Conduction

Science.gov (United States)

Riquelme, Mario; Quataert, Eliot; Verscharen, Daniel

2018-02-01

We use particle-in-cell (PIC) simulations of a collisionless, electron–ion plasma with a decreasing background magnetic field, {\\boldsymbol{B}}, to study the effect of velocity-space instabilities on the viscous heating and thermal conduction of the plasma. If | {\\boldsymbol{B}}| decreases, the adiabatic invariance of the magnetic moment gives rise to pressure anisotropies with {p}| | ,j> {p}\\perp ,j ({p}| | ,j and {p}\\perp ,j represent the pressure of species j (electron or ion) parallel and perpendicular to B ). Linear theory indicates that, for sufficiently large anisotropies, different velocity-space instabilities can be triggered. These instabilities in principle have the ability to pitch-angle scatter the particles, limiting the growth of the anisotropies. Our simulations focus on the nonlinear, saturated regime of the instabilities. This is done through the permanent decrease of | {\\boldsymbol{B}}| by an imposed plasma shear. We show that, in the regime 2≲ {β }j≲ 20 ({β }j\\equiv 8π {p}j/| {\\boldsymbol{B}}{| }2), the saturated ion and electron pressure anisotropies are controlled by the combined effect of the oblique ion firehose and the fast magnetosonic/whistler instabilities. These instabilities grow preferentially on the scale of the ion Larmor radius, and make {{Δ }}{p}e/{p}| | ,e≈ {{Δ }}{p}i/{p}| | ,i (where {{Δ }}{p}j={p}\\perp ,j-{p}| | ,j). We also quantify the thermal conduction of the plasma by directly calculating the mean free path of electrons, {λ }e, along the mean magnetic field, finding that {λ }e depends strongly on whether | {\\boldsymbol{B}}| decreases or increases. Our results can be applied in studies of low-collisionality plasmas such as the solar wind, the intracluster medium, and some accretion disks around black holes.

BOA, Beam Optics Analyzer A Particle-In-Cell Code

International Nuclear Information System (INIS)

Bui, Thuc

2007-01-01

The program was tasked with implementing time dependent analysis of charges particles into an existing finite element code with adaptive meshing, called Beam Optics Analyzer (BOA). BOA was initially funded by a DOE Phase II program to use the finite element method with adaptive meshing to track particles in unstructured meshes. It uses modern programming techniques, state-of-the-art data structures, so that new methods, features and capabilities are easily added and maintained. This Phase II program was funded to implement plasma simulations in BOA and extend its capabilities to model thermal electrons, secondary emissions, self magnetic field and implement a more comprehensive post-processing and feature-rich GUI. The program was successful in implementing thermal electrons, secondary emissions, and self magnetic field calculations. The BOA GUI was also upgraded significantly, and CCR is receiving interest from the microwave tube and semiconductor equipment industry for the code. Implementation of PIC analysis was partially successful. Computational resource requirements for modeling more than 2000 particles begin to exceed the capability of most readily available computers. Modern plasma analysis typically requires modeling of approximately 2 million particles or more. The problem is that tracking many particles in an unstructured mesh that is adapting becomes inefficient. In particular memory requirements become excessive. This probably makes particle tracking in unstructured meshes currently unfeasible with commonly available computer resources. Consequently, Calabazas Creek Research, Inc. is exploring hybrid codes where the electromagnetic fields are solved on the unstructured, adaptive mesh while particles are tracked on a fixed mesh. Efficient interpolation routines should be able to transfer information between nodes of the two meshes. If successfully developed, this could provide high accuracy and reasonable computational efficiency.

Development and Benchmarking of a Hybrid PIC Code For Dense Plasmas and Fast Ignition

Energy Technology Data Exchange (ETDEWEB)

Witherspoon, F. Douglas [HyperV Technologies Corp.; Welch, Dale R. [Voss Scientific, LLC; Thompson, John R. [FAR-TECH, Inc.; MacFarlane, Joeseph J. [Prism Computational Sciences Inc.; Phillips, Michael W. [Advanced Energy Systems, Inc.; Bruner, Nicki [Voss Scientific, LLC; Mostrom, Chris [Voss Scientific, LLC; Thoma, Carsten [Voss Scientific, LLC; Clark, R. E. [Voss Scientific, LLC; Bogatu, Nick [FAR-TECH, Inc.; Kim, Jin-Soo [FAR-TECH, Inc.; Galkin, Sergei [FAR-TECH, Inc.; Golovkin, Igor E. [Prism Computational Sciences, Inc.; Woodruff, P. R. [Prism Computational Sciences, Inc.; Wu, Linchun [HyperV Technologies Corp.; Messer, Sarah J. [HyperV Technologies Corp.

2014-05-20

Radiation processes play an important role in the study of both fast ignition and other inertial confinement schemes, such as plasma jet driven magneto-inertial fusion, both in their effect on energy balance, and in generating diagnostic signals. In the latter case, warm and hot dense matter may be produced by the convergence of a plasma shell formed by the merging of an assembly of high Mach number plasma jets. This innovative approach has the potential advantage of creating matter of high energy densities in voluminous amount compared with high power lasers or particle beams. An important application of this technology is as a plasma liner for the flux compression of magnetized plasma to create ultra-high magnetic fields and burning plasmas. HyperV Technologies Corp. has been developing plasma jet accelerator technology in both coaxial and linear railgun geometries to produce plasma jets of sufficient mass, density, and velocity to create such imploding plasma liners. An enabling tool for the development of this technology is the ability to model the plasma dynamics, not only in the accelerators themselves, but also in the resulting magnetized target plasma and within the merging/interacting plasma jets during transport to the target. Welch pioneered numerical modeling of such plasmas (including for fast ignition) using the LSP simulation code. Lsp is an electromagnetic, parallelized, plasma simulation code under development since 1995. It has a number of innovative features making it uniquely suitable for modeling high energy density plasmas including a hybrid fluid model for electrons that allows electrons in dense plasmas to be modeled with a kinetic or fluid treatment as appropriate. In addition to in-house use at Voss Scientific, several groups carrying out research in Fast Ignition (LLNL, SNL, UCSD, AWE (UK), and Imperial College (UK)) also use LSP. A collaborative team consisting of HyperV Technologies Corp., Voss Scientific LLC, FAR-TECH, Inc., Prism

A One-Dimensional Particle-in-Cell Model of Plasma Build-Up in Vacuum Arcs

CERN Document Server

Timko, H; Kovermann, J; Taborelli, M; Nordlund, K; Descoeudres, A; Schneider, R; Calatroni, S; Matyash, K; Wuensch, W; Hansen, A; Grudiev, A

2011-01-01

Understanding the mechanism of plasma build-up in vacuum arcs is essential in many fields of physics. A one-dimensional particle-in-cell computer simulation model is presented, which models the plasma developing from a field emitter tip under electrical breakdown conditions, taking into account the relevant physical phenomena. As a starting point, only an external electric field and an initial enhancement factor of the tip are assumed. General requirements for plasma formation have been identified and formulated in terms of the initial local field and a critical neutral density. The dependence of plasma build-up on tip melting current, the evaporation rate of neutrals and external circuit time constant has been investigated for copper and simulations imply that arcing involves melting currents around 0.5-1 A/mu m(2),evaporation of neutrals to electron field emission ratios in the regime 0.01 - 0.05, plasma build-up timescales in the order of similar to 1 - 10 ns and two different regimes depending on initial ...

Analysis of the beam halo in negative ion sources by using 3D3V PIC code

Energy Technology Data Exchange (ETDEWEB)

Miyamoto, K., E-mail: kmiyamot@naruto-u.ac.jp [Naruto University of Education, 748 Nakashima, Takashima, Naruto-cho, Naruto-shi, Tokushima 772-8502 (Japan); Nishioka, S.; Goto, I.; Hatayama, A. [Faculty of Science and Technology, Keio University, 3-14-1 Hiyoshi, Kohoku-ku, Yokohama 223-8522 (Japan); Hanada, M.; Kojima, A.; Hiratsuka, J. [Japan Atomic Energy Agency, 801-1 Mukouyama, Naka 319-0913 (Japan)

2016-02-15

The physical mechanism of the formation of the negative ion beam halo and the heat loads of the multi-stage acceleration grids are investigated with the 3D PIC (particle in cell) simulation. The following physical mechanism of the beam halo formation is verified: The beam core and the halo consist of the negative ions extracted from the center and the periphery of the meniscus, respectively. This difference of negative ion extraction location results in a geometrical aberration. Furthermore, it is shown that the heat loads on the first acceleration grid and the second acceleration grid are quantitatively improved compared with those for the 2D PIC simulation result.

Two-dimensional PIC simulations of ion beam instabilities in Supernova-driven plasma flows

Energy Technology Data Exchange (ETDEWEB)

Dieckmann, M E; Shukla, P K [Institut fuer Theoretische Physik IV, Ruhr-Universitaet Bochum, D-44780 Bochum (Germany); Meli, A; Mastichiadis, A [Department of Physics, National University of Athens, Panepistimiopolis, Zografos 15783 (Greece); Drury, L O C [Dublin Institute for Advanced Studies, Dublin 2 (Ireland)], E-mail: markd@tp4.rub.de

2008-06-15

Supernova remnant blast shells can reach the flow speed v{sub s} = 0.1c and shocks form at its front. Instabilities driven by shock-reflected ion beams heat the plasma in the foreshock, which may inject particles into diffusive acceleration. The ion beams can have the speed v{sub b} {approx} v{sub s}. For v{sub b} << v{sub s} the Buneman or upper-hybrid instabilities dominate, while for v{sub b} >> v{sub s} the filamentation and mixed modes grow faster. Here the relevant waves for v{sub b} {approx} v{sub s} are examined and how they interact nonlinearly with the particles. The collision of two plasma clouds at the speed v{sub s} is modelled with particle-in-cell simulations, which convect with them magnetic fields oriented perpendicular to their flow velocity vector. One simulation models equally dense clouds and the other one uses a density ratio of 2. Both simulations show upper-hybrid waves that are planar over large spatial intervals and that accelerate electrons to {approx}10 keV. The symmetric collision yields only short oscillatory wave pulses, while the asymmetric collision also produces large-scale electric fields, probably through a magnetic pressure gradient. The large-scale fields destroy the electron phase space holes and they accelerate the ions, which facilitates the formation of a precursor shock.

Towards the petascale in electromagnetic modeling of plasma-based accelerators for high-energy physics

International Nuclear Information System (INIS)

Bruhwiler, D L; Antonsen, T; Cary, J R; Cooley, J; Decyk, V K; Esarey, E; Geddes, C G R; Huang, C; Hakim, A; Katsouleas, T; Messmer, P; Mori, W B; Tsung, F S; Vieira, J; Zhou, M

2006-01-01

Plasma-based lepton acceleration concepts are a key element of the long-term R and D portfolio for the U.S. Office of High Energy Physics. There are many such concepts, but we consider only the laser (LWFA) and plasma (PWFA) wakefield accelerators. We present a summary of electromagnetic particle-in-cell (PIC) simulations for recent LWFA and PWFA experiments. These simulations, including both time explicit algorithms and reduced models, have effectively used terascale computing resources to support and guide experiments in this rapidly developing field. We briefly discuss the challenges and opportunities posed by the near-term availability of petascale computing hardware

Low-temperature plasma simulations with the LSP PIC code

Science.gov (United States)

Carlsson, Johan; Khrabrov, Alex; Kaganovich, Igor; Keating, David; Selezneva, Svetlana; Sommerer, Timothy

2014-10-01

The LSP (Large-Scale Plasma) PIC-MCC code has been used to simulate several low-temperature plasma configurations, including a gas switch for high-power AC/DC conversion, a glow discharge and a Hall thruster. Simulation results will be presented with an emphasis on code comparison and validation against experiment. High-voltage, direct-current (HVDC) power transmission is becoming more common as it can reduce construction costs and power losses. Solid-state power-electronics devices are presently used, but it has been proposed that gas switches could become a compact, less costly, alternative. A gas-switch conversion device would be based on a glow discharge, with a magnetically insulated cold cathode. Its operation is similar to that of a sputtering magnetron, but with much higher pressure (0.1 to 0.3 Torr) in order to achieve high current density. We have performed 1D (axial) and 2D (axial/radial) simulations of such a gas switch using LSP. The 1D results were compared with results from the EDIPIC code. To test and compare the collision models used by the LSP and EDIPIC codes in more detail, a validation exercise was performed for the cathode fall of a glow discharge. We will also present some 2D (radial/azimuthal) LSP simulations of a Hall thruster. The information, data, or work presented herein was funded in part by the Advanced Research Projects Agency-Energy (ARPA-E), U.S. Department of Energy, under Award Number DE-AR0000298.

Two-stream instability for a light ion beam-plasma system with external magnetic field

International Nuclear Information System (INIS)

Okada, T.; Tazawa, H.

1992-12-01

For inertial confinement fusion, a focused light ion beam (LIB) is required to propagate stably through a chamber to a target. We have pointed out that the applied external magnetic field is important for LIB propagation. To investigate the influence of the external magnetic field on the LIB propagation, we analysed the electrostatic dispersion relation of magnetized light ion beam-plasma system. The particle in-cell (PIC) simulation results are presented for a light ion beam-plasma system with external magnetic field. (author)

Hybrid 3D model for the interaction of plasma thruster plumes with nearby objects

Science.gov (United States)

Cichocki, Filippo; Domínguez-Vázquez, Adrián; Merino, Mario; Ahedo, Eduardo

2017-12-01

This paper presents a hybrid particle-in-cell (PIC) fluid approach to model the interaction of a plasma plume with a spacecraft and/or any nearby object. Ions and neutrals are modeled with a PIC approach, while electrons are treated as a fluid. After a first iteration of the code, the domain is split into quasineutral and non-neutral regions, based on non-neutrality criteria, such as the relative charge density and the Debye length-to-cell size ratio. At the material boundaries of the former quasineutral region, a dedicated algorithm ensures that the Bohm condition is met. In the latter non-neutral regions, the electron density and electric potential are obtained by solving the coupled electron momentum balance and Poisson equations. Boundary conditions for both the electric current and potential are finally obtained with a plasma sheath sub-code and an equivalent circuit model. The hybrid code is validated by applying it to a typical plasma plume-spacecraft interaction scenario, and the physics and capabilities of the model are finally discussed.

Particle Control in Phase Space by Global K-Means Clustering

DEFF Research Database (Denmark)

Frederiksen, Jacob Trier; Lapenta, G.; Pessah, M. E.

2015-01-01

We devise and explore an iterative optimization procedure for controlling particle populations in particle-in-cell (PIC) codes via merging and splitting of computational macro-particles. Our approach, is to compute an optimal representation of the global particle phase space structure while decre...

Test particle method for incorporation of the kinetic effects into the envelope simulations of Raman backscattering

International Nuclear Information System (INIS)

Hur, Min Sup; Suk, Hyyong

2007-01-01

A new test particle method is presented for self-consistent incorporation of the kinetic effects into the fluid three-wave model. One of the most important kinetic effects is the electron trapping and it has been found that the trapping affects significantly the behavior of Raman backscatter and Raman backward laser amplification. The conventional fluid three-wave model cannot reproduce the kinetic simulations in the trapping regime. The test particle scheme utilizes the same equations for the laser evolution as in the three-wave model. However, the plasma wave is treated by the envelope-kinetic equation, which consists of envelope evolution and the kinetic term. The core of the new scheme is employing test particles to compute the kinetic term self-consistently. The benchmarking results against the averaged particle-in-cell (aPIC) code show excellent agreements, and the computation speed gain over the aPIC is from 2 to 20 depending on parameters

Dusty-Plasma Particle Accelerator

Science.gov (United States)

Foster, John E.

2005-01-01

A dusty-plasma apparatus is being investigated as means of accelerating nanometer- and micrometer-sized particles. Applications for the dusty-plasma particle accelerators fall into two classes: Simulation of a variety of rapidly moving dust particles and micrometeoroids in outer-space environments that include micrometeoroid streams, comet tails, planetary rings, and nebulae and Deposition or implantation of nanoparticles on substrates for diverse industrial purposes that could include hardening, increasing thermal insulation, altering optical properties, and/or increasing permittivities of substrate materials. Relative to prior apparatuses used for similar applications, dusty-plasma particle accelerators offer such potential advantages as smaller size, lower cost, less complexity, and increased particle flux densities. A dusty-plasma particle accelerator exploits the fact that an isolated particle immersed in plasma acquires a net electric charge that depends on the relative mobilities of electrons and ions. Typically, a particle that is immersed in a low-temperature, partially ionized gas, wherein the average kinetic energy of electrons exceeds that of ions, causes the particle to become negatively charged. The particle can then be accelerated by applying an appropriate electric field. A dusty-plasma particle accelerator (see figure) includes a plasma source such as a radio-frequency induction discharge apparatus containing (1) a shallow cup with a biasable electrode to hold the particles to be accelerated and (2) a holder for the substrate on which the particles are to impinge. Depending on the specific design, a pair of electrostatic-acceleration grids between the substrate and discharge plasma can be used to both collimate and further accelerate particles exiting the particle holder. Once exposed to the discharge plasma, the particles in the cup quickly acquire a negative charge. Application of a negative voltage pulse to the biasable electrode results in the

Recent progress in 3D EM/EM-PIC simulation with ARGUS and parallel ARGUS

International Nuclear Information System (INIS)

Mankofsky, A.; Petillo, J.; Krueger, W.; Mondelli, A.; McNamara, B.; Philp, R.

1994-01-01

ARGUS is an integrated, 3-D, volumetric simulation model for systems involving electric and magnetic fields and charged particles, including materials embedded in the simulation region. The code offers the capability to carry out time domain and frequency domain electromagnetic simulations of complex physical systems. ARGUS offers a boolean solid model structure input capability that can include essentially arbitrary structures on the computational domain, and a modular architecture that allows multiple physics packages to access the same data structure and to share common code utilities. Physics modules are in place to compute electrostatic and electromagnetic fields, the normal modes of RF structures, and self-consistent particle-in-cell (PIC) simulation in either a time dependent mode or a steady state mode. The PIC modules include multiple particle species, the Lorentz equations of motion, and algorithms for the creation of particles by emission from material surfaces, injection onto the grid, and ionization. In this paper, we present an updated overview of ARGUS, with particular emphasis given in recent algorithmic and computational advances. These include a completely rewritten frequency domain solver which efficiently treats lossy materials and periodic structures, a parallel version of ARGUS with support for both shared memory parallel vector (i.e. CRAY) machines and distributed memory massively parallel MIMD systems, and numerous new applications of the code

Particle-in-cell simulations of asymmetric guide-field reconnection: quadrupolar structure of Hall magnetic field

Science.gov (United States)

Schmitz, R. G.; Alves, M. V.; Barbosa, M. V. G.

2017-12-01

One of the most important processes that occurs in Earth's magnetosphere is known as magnetic reconnection (MR). This process can be symmetric or asymmetric, depending basically on the plasma density and magnetic field in both sides of the current sheet. A good example of symmetric reconnection in terrestrial magnetosphere occurs in the magnetotail, where these quantities are similar on the north and south lobes. In the dayside magnetopause MR is asymmetric, since the plasma regimes and magnetic fields of magnetosheath and magnetosphere are quite different. Symmetric reconnection has some unique signatures. For example, the formation of a quadrupolar structure of Hall magnetic field and a bipolar Hall electric field that points to the center of the current sheet. The different particle motions in the presence of asymmetries change these signatures, causing the quadrupolar pattern to be distorted and forming a bipolar structure. Also, the bipolar Hall electric field is modified and gives rise to a single peak pointing toward the magnetosheat, considering an example of magnetopause reconnection. The presence of a guide-field can also distort the quadrupolar pattern, by giving a shear angle across the current sheet and altering the symmetric patterns, according to previous simulations and observations. Recently, a quadrupolar structure was observed in an asymmetric guide-field MR event using MMS (Magnetospheric Multiscale) mission data [Peng et al., JGR, 2017]. This event shows clearly that the density asymmetry and the guide-field were not sufficient to form signatures of asymmetric reconnection. Using the particle-in-cell code iPIC3D [Markidis et al, Mathematics and Computers in Simulation, 2010] with the MMS data from this event used to define input parameters, we found a quadrupolar structure of Hall magnetic field and a bipolar pattern of Hall electric field in ion scales, showing that our results are in an excellent agreement with the MMS observations. To our

Relativistic extension of a charge-conservative finite element solver for time-dependent Maxwell-Vlasov equations

Science.gov (United States)

Na, D.-Y.; Moon, H.; Omelchenko, Y. A.; Teixeira, F. L.

2018-01-01

Accurate modeling of relativistic particle motion is essential for physical predictions in many problems involving vacuum electronic devices, particle accelerators, and relativistic plasmas. A local, explicit, and charge-conserving finite-element time-domain (FETD) particle-in-cell (PIC) algorithm for time-dependent (non-relativistic) Maxwell-Vlasov equations on irregular (unstructured) meshes was recently developed by Moon et al. [Comput. Phys. Commun. 194, 43 (2015); IEEE Trans. Plasma Sci. 44, 1353 (2016)]. Here, we extend this FETD-PIC algorithm to the relativistic regime by implementing and comparing three relativistic particle-pushers: (relativistic) Boris, Vay, and Higuera-Cary. We illustrate the application of the proposed relativistic FETD-PIC algorithm for the analysis of particle cyclotron motion at relativistic speeds, harmonic particle oscillation in the Lorentz-boosted frame, and relativistic Bernstein modes in magnetized charge-neutral (pair) plasmas.

Particle-In-Cell simulations of the Ball-pen probe

Czech Academy of Sciences Publication Activity Database

Komm, M.; Adámek, Jiří; Pekárek, Z.; Pánek, Radomír

2010-01-01

Roč. 50, č. 9 (2010), s. 814-818 ISSN 0863-1042 R&D Projects: GA AV ČR KJB100430901 Institutional research plan: CEZ:AV0Z20430508 Keywords : Ball- pen * tokamak * plasma * plasma potential * PIC * simulation * I-V characteristics Subject RIV: BL - Plasma and Gas Discharge Physics Impact factor: 1.006, year: 2010 http://onlinelibrary.wiley.com/doi/10.1002/ctpp.201010137/pdf

Particle-in-cell simulations of plasma opening switch with external magnetic field

International Nuclear Information System (INIS)

Chen Yulan; Zeng Zhengzhong; Sun Fengju

2003-01-01

Fully electromagnetic particle-in-cell simulations are performed to study the effects of an external magnetic field on coaxial plasma opening switch (POS). The simulation results show that POS opening performance can be significantly improved only when external longitudinal magnetic field coils are placed at the cathode side, and an additional azimuthal magnetic field is effective whether the central electrode is of positive or negative polarity. Voltage multiplication coefficient K rises with the additional magnetic field increasing till the electron current is completely magnetically insulated during the opening of POS

Particle-in-cell simulations of high energy electron production by intense laser pulses in underdense plasmas

International Nuclear Information System (INIS)

Susumu, Kato; Eisuke, Miura; Kazuyoshi, Koyama; Mitsumori, Tanimoto; Masahiro, Adachi

2004-01-01

The propagation of intense laser pulses and the generation of high energy electrons from underdense plasmas are investigated using two dimensional particle-in-cell simulations. When the ratio of the laser power to the critical power of relativistic self-focusing gets the optimal value, the laser pulse propagates in a steady way and electrons have maximum energies. (author)

Particle-in-cell simulations of high energy electron production by intense laser pulses in underdense plasmas

Energy Technology Data Exchange (ETDEWEB)

Susumu, Kato; Eisuke, Miura; Kazuyoshi, Koyama [National Institute of Advanced Industrial Science and Technology (AIST), Tsukuba, Ibaraki (Japan); Mitsumori, Tanimoto [Meisei Univ., Dept. of Electrical Engineering, Hino, Tokyo (Japan); Masahiro, Adachi [Hiroshima Univ., Graduate school of Advanced Science of Matter, Higashi-Hiroshima, Hiroshima (Japan)

2004-07-01

The propagation of intense laser pulses and the generation of high energy electrons from underdense plasmas are investigated using two dimensional particle-in-cell simulations. When the ratio of the laser power to the critical power of relativistic self-focusing gets the optimal value, the laser pulse propagates in a steady way and electrons have maximum energies. (author)

Plasma ignition and steady state simulations of the Linac4 H$^{-}$ ion source

CERN Document Server

Mattei, S; Yasumoto, M; Hatayama, A; Lettry, J; Grudiev, A

2014-01-01

The RF heating of the plasma in the Linac4 H- ion source has been simulated using an Particle-in-Cell Monte Carlo Collision method (PIC-MCC). This model is applied to investigate the plasma formation starting from an initial low electron density of 1012 m-3 and its stabilization at 1018 m-3. The plasma discharge at low electron density is driven by the capacitive coupling with the electric field generated by the antenna, and as the electron density increases the capacitive electric field is shielded by the plasma and induction drives the plasma heating process. Plasma properties such as e-/ion densities and energies, sheath formation and shielding effect are presented and provide insight to the plasma properties of the hydrogen plasma.

Simulation of plasma erosion opening switches

International Nuclear Information System (INIS)

Mason, R.J.; Jones, M.E.

1988-01-01

The plasma erosion opening switch (PEOS) has been studied with the ANTHEM and ISIS implicit simulation codes. The switch consists of plasma fill injected into a transmission line. The plasma initially shorts out the circuit, but eventually it is removed by self-electrical forces, allowing for the delivery of energy to a load. ANTHEM models the plasma by multiple fluids with electron inertia retained, or by the particle-in-cell (PIC) technique. ISIS is an optimized PIC code. Both codes determine electric and magnetic fields by the implicit moment method. This allows for the study of long time full-switch behavior with simulational zone sizes and time steps that are large compared to a Debye length and plasma period, respectively. Thus, the authors have modeled switch behavior at densities ranging from 5 x 10 11 to 5 x 10 14 electrons/cm -3 over drive pulses ranging from 5 to 250 ns. Here, the magnetic field rose linearly from zero to 0.8 or 3.0 Tesla. Switch gaps spanned from 1.0 to 8.0 cm, and inner radii ranged from 0.5 to 20.0 cm. Opening dynamics is shown to depend sensitively on the assumed electron emission thresholds at the cathode, and on the effective conductivity of the anode. The particle simulations predict broader current channels than the multi-fluid calculations - reasons for this are discussed. The effect of numerical diffusion in implicit simulations is examined. The response to realistic load impedances (10 Ohms for Sandia National Laboratory's PBFA II accelerator) of the opening characteristics is described. Advantages from plasma fill near the load are investigated. The action of preset initial magnetic fields aligned with the power flow, and of trigger magnetic fields for controlled removal of the plasma is discussed

PIC 16 F84

International Nuclear Information System (INIS)

Jung, Gi Cheol; Min, Han Sik

2001-11-01

The contents of this book are introduction of microprocessor, basic for microcomputer practice, introduction of one chip micro computer, basic command of PIC, instructions simulator and in circuit emulator ; what a simulator of PIC is, and MPLAB direction, making PIC rom writer and instructions, of La's PIC Micro Programmer, PIC programming ; learning Command with examples, and controlling hardware with C-language, practical task for PIC application, a line tracer automobile and making ultrasonic radar ; circuit, source program and monitor program.

PIC 16 F84

Energy Technology Data Exchange (ETDEWEB)

Jung, Gi Cheol; Min, Han Sik

2001-11-15

The contents of this book are introduction of microprocessor, basic for microcomputer practice, introduction of one chip micro computer, basic command of PIC, instructions simulator and in circuit emulator ; what a simulator of PIC is, and MPLAB direction, making PIC rom writer and instructions, of La's PIC Micro Programmer, PIC programming ; learning Command with examples, and controlling hardware with C-language, practical task for PIC application, a line tracer automobile and making ultrasonic radar ; circuit, source program and monitor program.

A parallel electrostatic Particle-in-Cell method on unstructured tetrahedral grids for large-scale bounded collisionless plasma simulations

Science.gov (United States)

Averkin, Sergey N.; Gatsonis, Nikolaos A.

2018-06-01

An unstructured electrostatic Particle-In-Cell (EUPIC) method is developed on arbitrary tetrahedral grids for simulation of plasmas bounded by arbitrary geometries. The electric potential in EUPIC is obtained on cell vertices from a finite volume Multi-Point Flux Approximation of Gauss' law using the indirect dual cell with Dirichlet, Neumann and external circuit boundary conditions. The resulting matrix equation for the nodal potential is solved with a restarted generalized minimal residual method (GMRES) and an ILU(0) preconditioner algorithm, parallelized using a combination of node coloring and level scheduling approaches. The electric field on vertices is obtained using the gradient theorem applied to the indirect dual cell. The algorithms for injection, particle loading, particle motion, and particle tracking are parallelized for unstructured tetrahedral grids. The algorithms for the potential solver, electric field evaluation, loading, scatter-gather algorithms are verified using analytic solutions for test cases subject to Laplace and Poisson equations. Grid sensitivity analysis examines the L2 and L∞ norms of the relative error in potential, field, and charge density as a function of edge-averaged and volume-averaged cell size. Analysis shows second order of convergence for the potential and first order of convergence for the electric field and charge density. Temporal sensitivity analysis is performed and the momentum and energy conservation properties of the particle integrators in EUPIC are examined. The effects of cell size and timestep on heating, slowing-down and the deflection times are quantified. The heating, slowing-down and the deflection times are found to be almost linearly dependent on number of particles per cell. EUPIC simulations of current collection by cylindrical Langmuir probes in collisionless plasmas show good comparison with previous experimentally validated numerical results. These simulations were also used in a parallelization

Self-focusing, self modulation and stability properties of laser beam propagating in plasma: A variational approach

International Nuclear Information System (INIS)

Kaur, Ravinder; Gill, Tarsem Singh; Mahajan, Ranju

2010-01-01

Laboratory as well as Particle in cell (PIC) simulation experiments reveal the strong flow of energetic electrons co-moving with laser beam in laser plasma interaction. Equation governing the evolution of complex envelope in slowly varying envelope approximation is nonlinear parabolic equation. A Lagrangian for the problem is set up and assuming a trial Gaussian profile, we solve the reduced Lagrangian problem for beam width and curvature. Besides self-focusing and self-modulation of laser beam, we observe that stability properties of such plasma system are studied about equilibrium values using this variational approach. We obtained an eigen value equation, which is cubic in nature and investigated the criterion for stability using Hurwitz conditions for laser beam plasma system.

A study of x-ray emission from the anode region in a plasma focus device

International Nuclear Information System (INIS)

Jia Wang; Tsinchi Yang

1988-01-01

The physical process of x-ray emission from the anode region in a plasma focus device due to the interaction of a powerful electron beam with the metal anode and with ionised metallic vapour from the anode is investigated. The influence of the magnetic field of the beam is taken into consideration. A MC-PIC model (Monte Carlo-particle in cell) is proposed for the process, in which an electron-photon collision cascade is simulated by the MC approach and the time-dependent state of metallic vapour is determined by PIC computation. The time-resolved energy spectra and angular distributions of x-ray emission from the extending anode region are calculated. The time-integrated characteristics of the x-ray emission can be compared with the results of experiments as far as they are available. (author)

Electromagnetic particle in cell modeling of the plasma focus: Current sheath formation and lift off

International Nuclear Information System (INIS)

Seng, Y. S.; Lee, P.; Rawat, R. S.

2014-01-01

The shaping and formation of the current sheath takes place in the breakdown phase of a plasma focus device and critically controls the device performance. Electrostatic particle in cell codes, with magnetic effects ignored, have been used to model the breakdown phase. This Letter reports the successful development and implementation of an electromagnetic particle in cell (EMPIC) code, including magnetic effects self-consistently, to simulate the breakdown phase; from the ionization, localization and gliding discharge along the insulator to the time instant of current sheath lift off. The magnetic field was found to be appreciable from the time the current sheath came into contact with the anode with increased local current, initiating the voltage breakdown of the device as a result

A parallel code named NEPTUNE for 3D fully electromagnetic and pic simulations

International Nuclear Information System (INIS)

Dong Ye; Yang Wenyuan; Chen Jun; Zhao Qiang; Xia Fang; Ma Yan; Xiao Li; Sun Huifang; Chen Hong; Zhou Haijing; Mao Zeyao; Dong Zhiwei

2010-01-01

A parallel code named NEPTUNE for 3D fully electromagnetic and particle-in-cell (PIC) simulations is introduced, which could run on the Linux system with hundreds to thousand CPUs. NEPTUNE is suitable to simulate entire 3D HPM devices; many HPM devices are simulated and designed by using it. In NEPTUNE code, the electromagnetic fields are updated by using the finite-difference in time domain (FDTD) method of solving Maxwell equations and the particles are moved by using Buneman-Boris advance method of solving relativistic Newton-Lorentz equation. Electromagnetic fields and particles are coupled by using liner weighing interpolation PIC method, and the electric filed components are corrected by using Boris method of solve Poisson equation in order to ensure charge-conservation. NEPTUNE code could construct many complicated geometric structures, such as arbitrary axial-symmetric structures, plane transforming structures, slow-wave-structures, coupling holes, foils, and so on. The boundary conditions used in NEPTUNE code are introduced in brief, including perfectly electric conductor boundary, external wave boundary, and particle boundary. Finally, some typical HPM devices are simulated and test by using NEPTUNE code, including MILO, RBWO, VCO, and RKA. The simulation results are with correct and credible physical images, and the parallel efficiencies are also given. (authors)

Programming 16-Bit PIC Microcontrollers in C Learning to Fly the PIC 24

CERN Document Server

Di Jasio, Lucio

2011-01-01

New in the second edition: * MPLAB X support and MPLAB C for the PIC24F v3 and later libraries * I2C™ interface * 100% assembly free solutions * Improved video, PAL/NTSC * Improved audio, RIFF files decoding * PIC24F GA1, GA2, GB1 and GB2 support   Most readers will associate Microchip's name with the ubiquitous 8-bit PIC microcontrollers but it is the new 16-bit PIC24F family that is truly stealing the scene. Orders of magnitude increases of performance, memory size and the rich peripheral set make programming these devices in C a must. This new guide by Microchip insid

Efficient Modeling of Laser-Plasma Accelerators with INF&RNO

Energy Technology Data Exchange (ETDEWEB)

Benedetti, C.; Schroeder, C. B.; Esarey, E.; Geddes, C. G. R.; Leemans, W. P.

2010-06-01

The numerical modeling code INF&RNO (INtegrated Fluid& paRticle simulatioN cOde, pronounced"inferno") is presented. INF&RNO is an efficient 2D cylindrical code to model the interaction of a short laser pulse with an underdense plasma. The code is based on an envelope model for the laser while either a PIC or a fluid description can be used for the plasma. The effect of the laser pulse on the plasma is modeled with the time-averaged poderomotive force. These and other features allow for a speedup of 2-4 orders of magnitude compared to standard full PIC simulations while still retaining physical fidelity. The code has been benchmarked against analytical solutions and 3D PIC simulations and here a set of validation tests together with a discussion of the performances are presented.

Particle-in-cell simulations of the plasma interaction with poloidal gaps in the ITER divertor outer vertical target.

Czech Academy of Sciences Publication Activity Database

Komm, Michael; Gunn, J. P.; Dejarnac, Renaud; Pánek, Radomír; Pitts, R.A.; Podolník, Aleš

2017-01-01

Roč. 57, č. 12 (2017), č. článku 126047. ISSN 0029-5515 R&D Projects: GA ČR(CZ) GA16-14228S; GA MŠk(CZ) 8D15001 Grant - others:Ga MŠk(CZ) LM2015070 Institutional support: RVO:61389021 Keywords : tokamak * plasma * ITER * particle-in-cell * heat loads * monoblock Subject RIV: BL - Plasma and Gas Discharge Physics OBOR OECD: Fluids and plasma physics (including surface physics) Impact factor: 3.307, year: 2016 http://iopscience.iop.org/article/10.1088/1741-4326/aa8a9a/meta

Electron heating in low pressure capacitive discharges revisited

Energy Technology Data Exchange (ETDEWEB)

Kawamura, E.; Lieberman, M. A.; Lichtenberg, A. J. [Department of Electrical Engineering and Computer Sciences University of California, Berkeley, California 94720 (United States)

2014-12-15

The electrons in capacitively coupled plasmas (CCPs) absorb energy via ohmic heating due to electron-neutral collisions and stochastic heating due to momentum transfer from high voltage moving sheaths. We use Particle-in-Cell (PIC) simulations to explore these heating mechanisms and to compare the PIC results with available theories on ohmic and stochastic heating. The PIC results for ohmic heating show good agreement with the ohmic heating calculation of Lafleur et al. [Phys. Plasmas 20, 124503 (2013)]. The PIC results for stochastic heating in low pressure CCPs with collisionless sheaths show good agreement with the stochastic heating model of Kaganovich et al. [IEEE Trans. Plasma Sci. 34, 696 (2006)], which revises the hard wall asymptotic model of Lieberman [IEEE Trans. Plasma Sci. 16, 638 (1988)] by taking current continuity and bulk oscillation into account.

Electron heating in low pressure capacitive discharges revisited

International Nuclear Information System (INIS)

Kawamura, E.; Lieberman, M. A.; Lichtenberg, A. J.

2014-01-01

The electrons in capacitively coupled plasmas (CCPs) absorb energy via ohmic heating due to electron-neutral collisions and stochastic heating due to momentum transfer from high voltage moving sheaths. We use Particle-in-Cell (PIC) simulations to explore these heating mechanisms and to compare the PIC results with available theories on ohmic and stochastic heating. The PIC results for ohmic heating show good agreement with the ohmic heating calculation of Lafleur et al. [Phys. Plasmas 20, 124503 (2013)]. The PIC results for stochastic heating in low pressure CCPs with collisionless sheaths show good agreement with the stochastic heating model of Kaganovich et al. [IEEE Trans. Plasma Sci. 34, 696 (2006)], which revises the hard wall asymptotic model of Lieberman [IEEE Trans. Plasma Sci. 16, 638 (1988)] by taking current continuity and bulk oscillation into account

Electron heating in low pressure capacitive discharges revisited

Science.gov (United States)

Kawamura, E.; Lieberman, M. A.; Lichtenberg, A. J.

2014-12-01

The electrons in capacitively coupled plasmas (CCPs) absorb energy via ohmic heating due to electron-neutral collisions and stochastic heating due to momentum transfer from high voltage moving sheaths. We use Particle-in-Cell (PIC) simulations to explore these heating mechanisms and to compare the PIC results with available theories on ohmic and stochastic heating. The PIC results for ohmic heating show good agreement with the ohmic heating calculation of Lafleur et al. [Phys. Plasmas 20, 124503 (2013)]. The PIC results for stochastic heating in low pressure CCPs with collisionless sheaths show good agreement with the stochastic heating model of Kaganovich et al. [IEEE Trans. Plasma Sci. 34, 696 (2006)], which revises the hard wall asymptotic model of Lieberman [IEEE Trans. Plasma Sci. 16, 638 (1988)] by taking current continuity and bulk oscillation into account.

Particle in cell simulation of peaking switch for breakdown evaluation

Energy Technology Data Exchange (ETDEWEB)

Umbarkar, Sachin B.; Bindu, S.; Mangalvedekar, H.A.; Saxena, A.; Singh, N.M., E-mail: sachin.b.umbarkar@gmail.com [Department of Electric Engineering, Veermata Jijabai Technological Institute, Mumbai (India); Sharma, Archana; Saroj, P.C.; Mittal, K.C. [Accelerator Pulse Power Division, Bhabha Atomic Research Centre, Mumbai (India)

2014-07-01

Marx generator connected to peaking capacitor and peaking switch can generate Ultra-Wideband (UWB) radiation. A new peaking switch is designed for converting the existing nanosecond Marx generator to a UWB source. The paper explains the particle in cell (PIC) simulation for this peaking switch, using MAGIC 3D software. This peaking switch electrode is made up of copper tungsten material and is fixed inside the hermitically sealed derlin material. The switch can withstand a gas pressure up to 13.5 kg/cm{sup 2}. The lower electrode of the switch is connected to the last stage of the Marx generator. Initially Marx generator (without peaking stage) in air; gives the output pulse with peak amplitude of 113.75 kV and pulse rise time of 25 ns. Thus, we design a new peaking switch to improve the rise time of output pulse and to pressurize this peaking switch separately (i.e. Marx and peaking switch is at different pressure). The PIC simulation gives the particle charge density, current density, E counter plot, emitted electron current, and particle energy along the axis of gap between electrodes. The charge injection and electric field dependence on ionic dissociation phenomenon are briefly analyzed using this simulation. The model is simulated with different gases (N{sub 2}, H{sub 2}, and Air) under different pressure (2 kg/cm{sup 2}, 5 kg/cm{sup 2}, 10 kg/cm{sup 2}). (author)

Laser propagation and soliton generation in strongly magnetized plasmas

Energy Technology Data Exchange (ETDEWEB)

Feng, W.; Li, J. Q.; Kishimoto, Y. [Graduate School of Energy Science, Kyoto University, Gokasho, Uji, Kyoto 611-0011 (Japan)

2016-03-15

The propagation characteristics of various laser modes with different polarization, as well as the soliton generation in strongly magnetized plasmas are studied numerically through one-dimensional (1D) particle-in-cell (PIC) simulations and analytically by solving the laser wave equation. PIC simulations show that the laser heating efficiency substantially depends on the magnetic field strength, the propagation modes of the laser pulse and their intensities. Generally, large amplitude laser can efficiently heat the plasma with strong magnetic field. Theoretical analyses on the linear propagation of the laser pulse in both under-dense and over-dense magnetized plasmas are well confirmed by the numerical observations. Most interestingly, it is found that a standing or moving soliton with frequency lower than the laser frequency is generated in certain magnetic field strength and laser intensity range, which can greatly enhance the laser heating efficiency. The range of magnetic field strength for the right-hand circularly polarized (RCP) soliton formation with high and low frequencies is identified by solving the soliton equations including the contribution of ion's motion and the finite temperature effects under the quasi-neutral approximation. In the limit of immobile ions, the RCP soliton tends to be peaked and stronger as the magnetic field increases, while the enhanced soliton becomes broader as the temperature increases. These findings in 1D model are well validated by 2D simulations.

Plasma particle accelerators

International Nuclear Information System (INIS)

Dawson, J.M.

1988-01-01

The Superconducting Supercollider (SSC) will require an 87-kilometer accelerator ring to boost particles to 40 TeV. The SSC's size is due in part to the fact that its operating principle is the same one that has dominated accelerator design for 50 years: it guides particles by means of magnetic fields and propels them by strong electric fields. If one were to build an equally powerful but smaller accelerator, one would need to increase the strength of the guiding and propelling fields. Actually, however, conventional technology may not be able to provide significant increases in field strength. There are two reasons. First, the forces from magnetic fields are becoming greater than the structural forces that hold a magnetic material together; the magnets that produce these fields would themselves be torn apart. Second, the energy from electric fields is reaching the energies that bind electrons to atoms; it would tear electrons from nuclei in the accelerator's support structures. It is the electric field problem that plasma accelerators can overcome. Plasma particle accelerators are based on the principle that particles can be accelerated by the electric fields generated within a plasma. Because the plasma has already been ionized, plasma particle accelerators are not susceptible to electron dissociation. They can in theory sustain accelerating fields thousands of times stronger that conventional technologies. So far two methods for creating plasma waves for accelerators have been proposed and tested: the wakefield and the beat wave. Although promising electric fields have been produced, more research is necessary to determine whether plasma particle accelerators can compete with the existing accelerators. 7 figs

A fully-implicit Particle-In-Cell Monte Carlo Collision code for the simulation of inductively coupled plasmas

Science.gov (United States)

Mattei, S.; Nishida, K.; Onai, M.; Lettry, J.; Tran, M. Q.; Hatayama, A.

2017-12-01

We present a fully-implicit electromagnetic Particle-In-Cell Monte Carlo collision code, called NINJA, written for the simulation of inductively coupled plasmas. NINJA employs a kinetic enslaved Jacobian-Free Newton Krylov method to solve self-consistently the interaction between the electromagnetic field generated by the radio-frequency coil and the plasma response. The simulated plasma includes a kinetic description of charged and neutral species as well as the collision processes between them. The algorithm allows simulations with cell sizes much larger than the Debye length and time steps in excess of the Courant-Friedrichs-Lewy condition whilst preserving the conservation of the total energy. The code is applied to the simulation of the plasma discharge of the Linac4 H- ion source at CERN. Simulation results of plasma density, temperature and EEDF are discussed and compared with optical emission spectroscopy measurements. A systematic study of the energy conservation as a function of the numerical parameters is presented.

Particle melting and particle/plasma interactions in DC and RF plasmas: a modeling study. (Volumes I and II)

International Nuclear Information System (INIS)

Wei, D.Y.C.

1987-01-01

Integral process models were developed to predict particle melting in both DC and RF plasmas. Specifically, a numerical model has been developed to predict the temperature history of particles injected in a low pressure DC plasma jet. The temperature and velocity fields of the plasma jet are predicted as a free jet by solving the parabolized Navier-Stokes equations using a spatial marching scheme. Correction factors were introduced to take into account non continuum effects encountered in the low pressure environment. The plasma jet profiles as well as the particle/plasma interactions under different jet pressure ratios (from underexpanded to overexpanded) were investigated. The flow and temperature fields in the RF plasma torch are calculated using the axisymmetric Navier-Stokes equations based on the primitive variables, along with pseudo two-dimensional electromagnetic field equations. Particle trajectories and heat transfer characteristics in both DC and RF plasmas are calculated using predicted plasma jet profiles. Particle melting efficiencies in both DC and RF plasmas are evaluated and compared using model alloy systems. Based on the theoretical considerations, an alternative route of plasma spraying process (hybrid plasma spraying process) is proposed. An evaluation of particle melting in hybrid plasma jets had indicated that further improvement in deposit properties could be made

Cell-centered particle weighting algorithm for PIC simulations in a non-uniform 2D axisymmetric mesh

Science.gov (United States)

Araki, Samuel J.; Wirz, Richard E.

2014-09-01

Standard area weighting methods for particle-in-cell simulations result in systematic errors on particle densities for a non-uniform mesh in cylindrical coordinates. These errors can be significantly reduced by using weighted cell volumes for density calculations. A detailed description on the corrected volume calculations and cell-centered weighting algorithm in a non-uniform mesh is provided. The simple formulas for the corrected volume can be used for any type of quadrilateral and/or triangular mesh in cylindrical coordinates. Density errors arising from the cell-centered weighting algorithm are computed for radial density profiles of uniform, linearly decreasing, and Bessel function in an adaptive Cartesian mesh and an unstructured mesh. For all the density profiles, it is shown that the weighting algorithm provides a significant improvement for density calculations. However, relatively large density errors may persist at outermost cells for monotonically decreasing density profiles. A further analysis has been performed to investigate the effect of the density errors in potential calculations, and it is shown that the error at the outermost cell does not propagate into the potential solution for the density profiles investigated.

Reduced 3d modeling on injection schemes for laser wakefield acceleration at plasma scale lengths

Science.gov (United States)

Helm, Anton; Vieira, Jorge; Silva, Luis; Fonseca, Ricardo

2017-10-01

Current modelling techniques for laser wakefield acceleration (LWFA) are based on particle-in-cell (PIC) codes which are computationally demanding. In PIC simulations the laser wavelength λ0, in μm-range, has to be resolved over the acceleration lengths in meter-range. A promising approach is the ponderomotive guiding center solver (PGC) by only considering the laser envelope for laser pulse propagation. Therefore only the plasma skin depth λp has to be resolved, leading to speedups of (λp /λ0) 2. This allows to perform a wide-range of parameter studies and use it for λ0 Tecnologia (FCT), Portugal, through Grant No. PTDC/FIS-PLA/2940/2014 and PD/BD/105882/2014.

Plasma wall particle balance in Tore Supra

International Nuclear Information System (INIS)

Grisolia, C.; Ghendrih, P.; Pegourie, B.; Grosman, A.

1992-01-01

A comprehensive study of the particle balance between the carbon wall and the plasma is presented. One finds that the effective particle content of the wall which governs the plasma equilibrium density departs from the deposited number of particles. This effect is dominant for the fully desaturated wall. A scaling law of the plasma density in terms of the wall effective particle content has been obtained. Moreover, the experimental data allows to estimate the plasma particle confinement time. Values ranging from 0.2 s to 0.5 s are found depending on the density. An analytical functional dependence of the particle confinement time is obtained

Localized Models of Charged Particle Motion in Martian Crustal Magnetic Cusps

Science.gov (United States)

Brain, D. A.; Poppe, A. R.; Jarvinen, R.; Dong, Y.; Egan, H. L.; Fang, X.

2017-12-01

The induced magnetosphere of Mars is punctuated by localized but strong crustal magnetic fields that are observed to play host to a variety of phenomena typically associated with global magnetic fields, such as auroral processes and particle precipitation, field-aligned current systems, and ion outflow. Each of these phenomena occur on the night side, in small-scale magnetic `cusp' regions of vertically aligned field. Cusp regions are not yet capable of being spatially resolved in global scale models that include the ion kinetics necessary for simulating charged particle transport along cusps. Local models are therefore necessary if we are to understand how cusp processes operate at Mars. Here we present the first results of an effort to model the kinetic particle motion and electric fields in Martian cusps. We are adapting both a 1.5D Particle-in-Cell (PIC) model for lunar magnetic cusps regions to the Martian case and a hybrid model framework (used previously for the global Martian plasma interaction and for lunar magnetic anomaly regions) to cusps in 2D. By comparing the models we can asses the importance of electron kinetics in particle transport along cusp field lines. In this first stage of our study we model a moderately strong nightside cusp, with incident hot hydrogen plasma from above, and cold planetary (oxygen) plasma entering the simulation from below. We report on the spatial and temporal distribution of plasma along cusp field lines for this initial case.

Plasma Wakefield Acceleration of an Intense Positron Beam

Energy Technology Data Exchange (ETDEWEB)

Blue, B

2004-04-21

The Plasma Wakefield Accelerator (PWFA) is an advanced accelerator concept which possess a high acceleration gradient and a long interaction length for accelerating both electrons and positrons. Although electron beam-plasma interactions have been extensively studied in connection with the PWFA, very little work has been done with respect to positron beam-plasma interactions. This dissertation addresses three issues relating to a positron beam driven plasma wakefield accelerator. These issues are (a) the suitability of employing a positron drive bunch to excite a wake; (b) the transverse stability of the drive bunch; and (c) the acceleration of positrons by the plasma wake that is driven by a positron bunch. These three issues are explored first through computer simulations and then through experiments. First, a theory is developed on the impulse response of plasma to a short drive beam which is valid for small perturbations to the plasma density. This is followed up with several particle-in-cell (PIC) simulations which study the experimental parameter (bunch length, charge, radius, and plasma density) range. Next, the experimental setup is described with an emphasis on the equipment used to measure the longitudinal energy variations of the positron beam. Then, the transverse dynamics of a positron beam in a plasma are described. Special attention is given to the way focusing, defocusing, and a tilted beam would appear to be energy variations as viewed on our diagnostics. Finally, the energy dynamics imparted on a 730 {micro}m long, 40 {micro}m radius, 28.5 GeV positron beam with 1.2 x 10{sup 10} particles in a 1.4 meter long 0-2 x 10{sup 14} e{sup -}/cm{sup 3} plasma is described. First the energy loss was measured as a function of plasma density and the measurements are compared to theory. Then, an energy gain of 79 {+-} 15 MeV is shown. This is the first demonstration of energy gain of a positron beam in a plasma and it is in good agreement with the predictions

Spectroscopic investigation of species separation in opening switch plasmas

Science.gov (United States)

Jackson, S. L.; Phipps, D. G.; Richardson, A. S.; Commisso, R. J.; Hinshelwood, D. D.; Murphy, D. P.; Schumer, J. W.; Weber, B. V.; Boyer, C. N.; Doron, R.; Biswas, S.; Maron, Y.

2015-11-01

Interactions between magnetic fields and current-carrying plasmas that lead to the separation of plasma species in multi-species plasmas are being studied in a plasma opening switch geometry. Several Marshall guns are used to inject single or multi-species plasmas between coaxial conductors connected to the output of the Naval Research Laboratory's Hawk pulsed-power generator. Following injection of the plasma, the generator is used at roughly half power to apply an electrical pulse with a peak current of 450 kA, a peak voltage of 400 kV, and a rise time of 1.2 μs. The resulting magnetic field interacts with the plasma through a combination of field penetration and magnetohydrodynamic (MHD) pushing that is not well understood but can lead to the separation of plasma species in multi-species plasmas. An ICCD-coupled spectrometer has been used in combination with magnetic probes, a ribbon-beam interferometer, and particle-in-cell (PIC) modeling to diagnose and understand conditions in the plasma from the time it is injected until the end of the conduction phase of the opening switch. This work supported by the Naval Research Laboratory Base Program and the Office of Naval Research.

Particle-in-cell simulations of multi-MeV pulsed X-ray induced air plasmas at low pressures

International Nuclear Information System (INIS)

Ribière, M.; D'Almeida, T.; Gaufridy de Dortan, F. de; Maulois, M.; Delbos, C.; Garrigues, A.; Cessenat, O.; Azaïs, B.

2016-01-01

A full kinetic modelling of the charge particles dynamics generated upon the irradiation of an air-filled cavity by a multi-MeV pulsed x-ray is performed. From the calculated radiative source generated by the ASTERIX generator, we calculated the electromagnetic fields generated by x-ray induced air plasmas in a metallic cavity at different pressures. Simulations are carried out based on a Particle-In-Cell interpolation method which uses 3D Maxwell-Vlasov calculations of the constitutive charged species densities of air plasmas at different pressures at equilibrium. The resulting electromagnetic fields within the cavity are calculated for different electron densities up to 4 × 10"1"0" cm"−"3. For each air pressure, we show electronic plasma waves formation followed by Landau damping. As electron density increases, the calculations exhibit space-charged neutralization and return current formation.

Progress in the study of mesh refinement for particle-in-cell plasma simulations and its application to heavy ion fusion

International Nuclear Information System (INIS)

Vay, J.-L.; Friedman, A.; Grote, D.P.

2002-01-01

The numerical simulation of the driving beams in a heavy ion fusion power plant is a challenging task, and, despite rapid progress in computer power, one must consider the use of the most advanced numerical techniques. One of the difficulties of these simulations resides in the disparity of scales in time and in space which must be resolved. When these disparities are in distinctive zones of the simulation region, a method which has proven to be effective in other areas (e.g. fluid dynamics simulations) is the Adaptive-Mesh-Refinement (AMR) technique. We follow in this article the progress accomplished in the last few months in the merging of the AMR technique with Particle-In-Cell (PIC) method. This includes a detailed modeling of the Lampel-Tiefenback solution for the one-dimensional diode using novel techniques to suppress undesirable numerical oscillations and an AMR patch to follow the head of the particle distribution. We also report new results concerning the modeling of ion sources using the axisymmetric WARPRZ-AMR prototype showing the utility of an AMR patch resolving the emitter vicinity and the beam edge

Comparison of multi-fluid moment models with particle-in-cell simulations of collisionless magnetic reconnection

International Nuclear Information System (INIS)

Wang, Liang; Germaschewski, K.; Hakim, Ammar H.; Bhattacharjee, A.

2015-01-01

We introduce an extensible multi-fluid moment model in the context of collisionless magnetic reconnection. This model evolves full Maxwell equations and simultaneously moments of the Vlasov-Maxwell equation for each species in the plasma. Effects like electron inertia and pressure gradient are self-consistently embedded in the resulting multi-fluid moment equations, without the need to explicitly solving a generalized Ohm's law. Two limits of the multi-fluid moment model are discussed, namely, the five-moment limit that evolves a scalar pressures for each species and the ten-moment limit that evolves the full anisotropic, non-gyrotropic pressure tensor for each species. We first demonstrate analytically and numerically that the five-moment model reduces to the widely used Hall magnetohydrodynamics (Hall MHD) model under the assumptions of vanishing electron inertia, infinite speed of light, and quasi-neutrality. Then, we compare ten-moment and fully kinetic particle-in-cell (PIC) simulations of a large scale Harris sheet reconnection problem, where the ten-moment equations are closed with a local linear collisionless approximation for the heat flux. The ten-moment simulation gives reasonable agreement with the PIC results regarding the structures and magnitudes of the electron flows, the polarities and magnitudes of elements of the electron pressure tensor, and the decomposition of the generalized Ohm's law. Possible ways to improve the simple local closure towards a nonlocal fully three-dimensional closure are also discussed

Surface roughness effects on plasma near a divertor plate and local impact angle

Directory of Open Access Journals (Sweden)

Wanpeng Hu

2017-08-01

Full Text Available The impact of rough surface topography on the electric potential and electric field is generally neglected due to the small scale of surface roughness compared to the width of the plasma sheath. However, the distributions of the electric potential and field on rough surfaces are expected to influence the characteristics of edge plasma and the local impact angle. The distributions of plasma sheath and local impact angle on rough surfaces are investigated by a two dimension-in-space and three dimension-in-velocity (2d3v Particle-In-Cell (PIC code. The influences of the plasma temperature andsurface morphology on the plasma sheath, local impact angle and resulting physical sputtering yield on rough surfaces are investigated.

Study of negative hydrogen ion beam optics using the 3D3V PIC model

International Nuclear Information System (INIS)

Miyamoto, K.; Nishioka, S.; Goto, I.; Hatayama, A.; Hanada, M.; Kojima, A.

2015-01-01

The mechanism of negative ion extraction under real conditions with the complex magnetic field is studied by using the 3D PIC simulation code. The extraction region of the negative ion source for the negative ion based neutral beam injection system in fusion reactors is modelled. It is shown that the E x B drift of electrons is caused by the magnetic filter and the electron suppression magnetic field, and the resultant asymmetry of the plasma meniscus. Furthermore, it is indicated that that the asymmetry of the plasma meniscus results in the asymmetry of negative ion beam profile including the beam halo. It could be demonstrated theoretically that the E x B drift is not significantly weakened by the elastic collisions of the electrons with neutral particles

Study of negative hydrogen ion beam optics using the 3D3V PIC model

Energy Technology Data Exchange (ETDEWEB)

Miyamoto, K., E-mail: kmiyamot@naruto-u.ac.jp [Naruto University of Education, 748 Nakashima, Takashima, Naruto-cho, Naruto-shi, Tokushima, 772-8502 (Japan); Nishioka, S.; Goto, I.; Hatayama, A. [Faculty of Science and Technology, Keio University, 3-14-1, Hiyoshi, Kohoku-ku, Yokohama, 223-8522 (Japan); Hanada, M.; Kojima, A. [Japan Atomic Energy Agency, 801-1,Mukoyama, Naka, 319-0913 (Japan)

2015-04-08

The mechanism of negative ion extraction under real conditions with the complex magnetic field is studied by using the 3D PIC simulation code. The extraction region of the negative ion source for the negative ion based neutral beam injection system in fusion reactors is modelled. It is shown that the E x B drift of electrons is caused by the magnetic filter and the electron suppression magnetic field, and the resultant asymmetry of the plasma meniscus. Furthermore, it is indicated that that the asymmetry of the plasma meniscus results in the asymmetry of negative ion beam profile including the beam halo. It could be demonstrated theoretically that the E x B drift is not significantly weakened by the elastic collisions of the electrons with neutral particles.

Massive parallel 3D PIC simulation of negative ion extraction

Science.gov (United States)

Revel, Adrien; Mochalskyy, Serhiy; Montellano, Ivar Mauricio; Wünderlich, Dirk; Fantz, Ursel; Minea, Tiberiu

2017-09-01

The 3D PIC-MCC code ONIX is dedicated to modeling Negative hydrogen/deuterium Ion (NI) extraction and co-extraction of electrons from radio-frequency driven, low pressure plasma sources. It provides valuable insight on the complex phenomena involved in the extraction process. In previous calculations, a mesh size larger than the Debye length was used, implying numerical electron heating. Important steps have been achieved in terms of computation performance and parallelization efficiency allowing successful massive parallel calculations (4096 cores), imperative to resolve the Debye length. In addition, the numerical algorithms have been improved in terms of grid treatment, i.e., the electric field near the complex geometry boundaries (plasma grid) is calculated more accurately. The revised model preserves the full 3D treatment, but can take advantage of a highly refined mesh. ONIX was used to investigate the role of the mesh size, the re-injection scheme for lost particles (extracted or wall absorbed), and the electron thermalization process on the calculated extracted current and plasma characteristics. It is demonstrated that all numerical schemes give the same NI current distribution for extracted ions. Concerning the electrons, the pair-injection technique is found well-adapted to simulate the sheath in front of the plasma grid.

Canonical momenta and numerical instabilities in particle codes

International Nuclear Information System (INIS)

Godfrey, B.B.

1975-01-01

A set of warm plasma dispersion relations appropriate to a large class of electromagnetic plasma simulation codes is derived. The numerical Cherenkov instability is shown by analytic and numerical analysis of these dispersion relations to be the most significant nonphysical effect involving transverse electromagnetic waves. The instability arises due to a spurious phase shift between resonant particles and light waves, caused by a basic incompatibility between the Lagrangian treatment of particle positions and the Eulerian treatment of particle velocities characteristic of most PIC--CIC algorithms. It is demonstrated that, through the use of canonical momentum, this mismatch is alleviated sufficiently to completely eliminate the Cherenkov instability. Collateral effects on simulation accuracy and on other numerical instabilities appear to be minor

Concave pulse shaping of a circularly polarized laser pulse from non-uniform overdense plasmas

Energy Technology Data Exchange (ETDEWEB)

Hur, Min Sup [School of Natural Science, UNIST, BanYeon-Ri 100, Ulju-gun, Ulsan, 689-798 (Korea, Republic of); Kulagin, Victor V. [Sternberg Astronomical Institute, Moscow State University, Universitetsky prosp. 13, Moscow, 119992 (Russian Federation); Suk, Hyyong, E-mail: hysuk@gist.ac.kr [Department of Physics and Photon Science, GIST, 123 Cheomdan-gwangiro, Buk-gu, Gwangju, 500-712 (Korea, Republic of)

2015-03-20

Pulse shaping of circularly polarized laser pulses in nonuniform overdense plasmas are investigated numerically. Specifically we show by two-dimensional particle-in-cell simulations the generation of a concave pulse front of a circularly polarized, a few tens of petawatt laser pulse from a density-tapered, overdense plasma slab. The concept used for the transverse-directional shaping is the differential transmittance depending on the plasma density, and the laser intensity. For suitable selection of the slab parameters for the concave pulse shaping, we studied numerically the pulse transmittance, which can be used for further parameter design of the pulse shaping. The concavely shaped circularly polarized pulse is expected to add more freedom in controlling the ion-beam characteristics in the RPDA regime. - Highlights: • Laser pulse shaping for a concave front by non-uniform overdense plasma was studied. • Particle-in-cell (PIC) simulations were used for the investigation. • A laser pulse can be shaped by a density-tapered overdense plasma. • The concave and sharp pulse front are useful in many laser–plasma applications. • They are important for ion acceleration, especially in the radiation pressure dominant regime.

Resonant and non-resonant whistlers-particle interaction in the radiation belts

NARCIS (Netherlands)

E. Camporeale (Enrico)

2014-01-01

htmlabstractWe study the wave-particle interactions between lower band chorus whistlers and an anisotropic tenuous population of relativistic electrons. We present the first direct comparison of first-principle Particle-in-Cell (PIC) simulations with a quasi-linear diffusion code, in this context.

Dynamics of the spectral behaviour of an ultrashort laser pulse in an argon-gas-filled capillary discharge-preformed plasma channel

Directory of Open Access Journals (Sweden)

Sakai S.

2013-11-01

Full Text Available We have reported the argon plasma waveguide produced in an alumina (Al2O3 capillary discharge and used to guide ultrashort laser pulses at intensities of the order of 1016  W/cm2. A one-dimensional magnetohydrodynamic (MHD code was used to evaluate the average degree of ionization of Ar in the preformed plasma channel. The spectrum of the propagated laser pulse in the Ar plasma waveguide was not modified and was well reproduced by a particle-in-cell (PIC simulation under initial ion charge state of Ar3+ in the preformed plasma waveguide. The optimum timing for the laser pulse injection was around 150 ns after initiation of a discharge with a peak current of 200 A.

Modeling of negative ion extraction from a magnetized plasma source: Derivation of scaling laws and description of the origins of aberrations in the ion beam

Science.gov (United States)

Fubiani, G.; Garrigues, L.; Boeuf, J. P.

2018-02-01

We model the extraction of negative ions from a high brightness high power magnetized negative ion source. The model is a Particle-In-Cell (PIC) algorithm with Monte-Carlo Collisions. The negative ions are generated only on the plasma grid surface (which separates the plasma from the electrostatic accelerator downstream). The scope of this work is to derive scaling laws for the negative ion beam properties versus the extraction voltage (potential of the first grid of the accelerator) and plasma density and investigate the origins of aberrations on the ion beam. We show that a given value of the negative ion beam perveance correlates rather well with the beam profile on the extraction grid independent of the simulated plasma density. Furthermore, the extracted beam current may be scaled to any value of the plasma density. The scaling factor must be derived numerically but the overall gain of computational cost compared to performing a PIC simulation at the real plasma density is significant. Aberrations appear for a meniscus curvature radius of the order of the radius of the grid aperture. These aberrations cannot be cancelled out by switching to a chamfered grid aperture (as in the case of positive ions).

Resolution of the Vlasov-Maxwell system by PIC discontinuous Galerkin method on GPU with OpenCL

Directory of Open Access Journals (Sweden)

Crestetto Anaïs

2013-01-01

Full Text Available We present an implementation of a Vlasov-Maxwell solver for multicore processors. The Vlasov equation describes the evolution of charged particles in an electromagnetic field, solution of the Maxwell equations. The Vlasov equation is solved by a Particle-In-Cell method (PIC, while the Maxwell system is computed by a Discontinuous Galerkin method. We use the OpenCL framework, which allows our code to run on multicore processors or recent Graphic Processing Units (GPU. We present several numerical applications to two-dimensional test cases.

SPECT3D - A multi-dimensional collisional-radiative code for generating diagnostic signatures based on hydrodynamics and PIC simulation output

Science.gov (United States)

MacFarlane, J. J.; Golovkin, I. E.; Wang, P.; Woodruff, P. R.; Pereyra, N. A.

2007-05-01

SPECT3D is a multi-dimensional collisional-radiative code used to post-process the output from radiation-hydrodynamics (RH) and particle-in-cell (PIC) codes to generate diagnostic signatures (e.g. images, spectra) that can be compared directly with experimental measurements. This ability to post-process simulation code output plays a pivotal role in assessing the reliability of RH and PIC simulation codes and their physics models. SPECT3D has the capability to operate on plasmas in 1D, 2D, and 3D geometries. It computes a variety of diagnostic signatures that can be compared with experimental measurements, including: time-resolved and time-integrated spectra, space-resolved spectra and streaked spectra; filtered and monochromatic images; and X-ray diode signals. Simulated images and spectra can include the effects of backlighters, as well as the effects of instrumental broadening and time-gating. SPECT3D also includes a drilldown capability that shows where frequency-dependent radiation is emitted and absorbed as it propagates through the plasma towards the detector, thereby providing insights on where the radiation seen by a detector originates within the plasma. SPECT3D has the capability to model a variety of complex atomic and radiative processes that affect the radiation seen by imaging and spectral detectors in high energy density physics (HEDP) experiments. LTE (local thermodynamic equilibrium) or non-LTE atomic level populations can be computed for plasmas. Photoabsorption rates can be computed using either escape probability models or, for selected 1D and 2D geometries, multi-angle radiative transfer models. The effects of non-thermal (i.e. non-Maxwellian) electron distributions can also be included. To study the influence of energetic particles on spectra and images recorded in intense short-pulse laser experiments, the effects of both relativistic electrons and energetic proton beams can be simulated. SPECT3D is a user-friendly software package that runs

AP-Cloud: Adaptive Particle-in-Cloud method for optimal solutions to Vlasov–Poisson equation

International Nuclear Information System (INIS)

Wang, Xingyu; Samulyak, Roman; Jiao, Xiangmin; Yu, Kwangmin

2016-01-01

We propose a new adaptive Particle-in-Cloud (AP-Cloud) method for obtaining optimal numerical solutions to the Vlasov–Poisson equation. Unlike the traditional particle-in-cell (PIC) method, which is commonly used for solving this problem, the AP-Cloud adaptively selects computational nodes or particles to deliver higher accuracy and efficiency when the particle distribution is highly non-uniform. Unlike other adaptive techniques for PIC, our method balances the errors in PDE discretization and Monte Carlo integration, and discretizes the differential operators using a generalized finite difference (GFD) method based on a weighted least square formulation. As a result, AP-Cloud is independent of the geometric shapes of computational domains and is free of artificial parameters. Efficient and robust implementation is achieved through an octree data structure with 2:1 balance. We analyze the accuracy and convergence order of AP-Cloud theoretically, and verify the method using an electrostatic problem of a particle beam with halo. Simulation results show that the AP-Cloud method is substantially more accurate and faster than the traditional PIC, and it is free of artificial forces that are typical for some adaptive PIC techniques.

Temporally asymmetric laser pulse for magnetic-field generation in plasmas

Energy Technology Data Exchange (ETDEWEB)

Singh, Mamta; Gopal, Krishna; Gupta, Devki Nandan, E-mail: dngupta@physics.du.ac.in

2016-04-01

Of particular interest in this article, the case study of an asymmetric laser pulse interaction with a plasma for magnetic field enhancement has been investigated. The strong ponderomotive force due to the short leading edge of the propagating laser pulse drives a large nonlinear current, producing a stronger quasistatic magnetic field. An analytical expression for the magnetic field is derived and the strength of the magnetic field is estimated for the current laser-plasma parameters. The theoretical results are validated through the particle-in-cell (PIC) simulations and are in very close agreement with the simulation based estimations. This kind of magnetic field can be useful in the plasma based accelerators as well as in the laser-fusion based experiments. - Highlights: • We employ an asymmetric laser pulse to enhance the magnetic field strength in a plasma. • Short leading front of the pulse drives a strong ponderomotive force. • An analytical expression for the magnetic field is derived. • The strength of the magnetic field is estimated for the current laser–plasma parameters.

Temporally asymmetric laser pulse for magnetic-field generation in plasmas

International Nuclear Information System (INIS)

Singh, Mamta; Gopal, Krishna; Gupta, Devki Nandan

2016-01-01

Of particular interest in this article, the case study of an asymmetric laser pulse interaction with a plasma for magnetic field enhancement has been investigated. The strong ponderomotive force due to the short leading edge of the propagating laser pulse drives a large nonlinear current, producing a stronger quasistatic magnetic field. An analytical expression for the magnetic field is derived and the strength of the magnetic field is estimated for the current laser-plasma parameters. The theoretical results are validated through the particle-in-cell (PIC) simulations and are in very close agreement with the simulation based estimations. This kind of magnetic field can be useful in the plasma based accelerators as well as in the laser-fusion based experiments. - Highlights: • We employ an asymmetric laser pulse to enhance the magnetic field strength in a plasma. • Short leading front of the pulse drives a strong ponderomotive force. • An analytical expression for the magnetic field is derived. • The strength of the magnetic field is estimated for the current laser–plasma parameters.

Effects of Temperature and Residence Time on the Emissions of PIC and Fine Particles during Fixed Bed Combustion of Conifer Stemwood Pellets

Energy Technology Data Exchange (ETDEWEB)

Boman, Christoffer; Lindmark, Fredrik; Oehman, Marcus; Nordin, Anders [Umeaa Univ. (Sweden). Energy Technology and Thermal Process Chemistry; Pettersson, Esbjoern [Energy Technology Centre, Piteaa (Sweden); Westerholm, Roger [Stockholm Univ., Arrhenius Laboratory (Sweden). Dept. of Analytical Chemistry

2006-07-15

The use of wood fuel Pellets has proved to be well suited for the small-scale market enabling controlled and efficient combustion with low emission of products of incomplete combustion (PIC). Still a potential for further emission reduction exists and a thorough understanding of the influence of combustion conditions on the emission characteristics of air pollutants like PAH and particulate matter (PM) is important. The objective was to determine the effects of temperature and residence time on the emission performance and characteristics with focus on hydrocarbons and PM during combustion of conifer stemwood Pellets in a laboratory fixed bed reactor (<5 kW). Temperature and residence time after the bed section were varied according to statistical experimental designs (650-970 deg C and 0.5-3.5 s) with the emission responses; CO, organic gaseous carbon, NO, 20 VOC compounds, 43 PAH compounds, PM{sub tot}, fine particle mass/count median diameter (MMD and CMD) and number concentration. Temperature was negatively correlated with the emissions of all studied PIC with limited effects of residence time. The PM{sub tot} emissions of 15-20 mg/MJ was in all cases dominated by fine (<1 {mu}m) particles of K, Na, S, Cl, C, O and Zn. Increased residence time resulted in increased fine particle sizes (i.e. MMD and CMD) and decreased number concentrations. The importance of high temperature (>850 deg C) in the bed zone with intensive, air rich and well mixed isothermal conditions for 0.5-1.0 s in the post combustion zone was illustrated for wood Pellets combustion with almost a total depletion of all studied PIC. The results emphasize the need for further verification studies and technology development work.

The plasma-wall transition layers in the presence of collisions with a magnetic field parallel to the wall

Science.gov (United States)

Moritz, J.; Faudot, E.; Devaux, S.; Heuraux, S.

2018-01-01

The plasma-wall transition is studied by means of a particle-in-cell (PIC) simulation in the configuration of a parallel to the wall magnetic field (B), with collisions between charged particles vs. neutral atoms taken into account. The investigated system consists of a plasma bounded by two absorbing walls separated by 200 electron Debye lengths (λd). The strength of the magnetic field is chosen such as the ratio λ d / r l , with rl being the electron Larmor radius, is smaller or larger than unity. Collisions are modelled with a simple operator that reorients randomly ion or electron velocity, keeping constant the total kinetic energy of both the neutral atom (target) and the incident charged particle. The PIC simulations show that the plasma-wall transition consists in a quasi-neutral region (pre-sheath), from the center of the plasma towards the walls, where the electric potential or electric field profiles are well described by an ambipolar diffusion model, and in a second region at the vicinity of the walls, called the sheath, where the quasi-neutrality breaks down. In this peculiar geometry of B and for a certain range of the mean-free-path, the sheath is found to be composed of two charged layers: the positive one, close to the walls, and the negative one, towards the plasma and before the neutral pre-sheath. Depending on the amplitude of B, the spatial variation of the electric potential can be non-monotonic and presents a maximum within the sheath region. More generally, the sheath extent as well as the potential drop within the sheath and the pre-sheath is studied with respect to B, the mean-free-path, and the ion and electron temperatures.

Two-dimensional particle-in-cell simulation of the expansion of a plasma into a rarefied medium

Energy Technology Data Exchange (ETDEWEB)

Sarri, G; Quinn, K; Kourakis, I; Borghesi, M [Centre for Plasma Physics, The Queens University of Belfast, Belfast BT7 1NN (United Kingdom); Murphy, G C; Drury, L O C [Dublin Institute for Advanced Studies, 31 Fitzwilliam Place, Dublin 2 (Ireland); Dieckmann, M E; Ynnerman, A [Department of Science and Technology (ITN), Linkoeping University, 60174 Norrkoping (Sweden); Bret, A, E-mail: gsarri01@qub.ac.uk [ETSI Industriales, Universidad de Castilla-La Mancha, 13071 Ciudad Real (Spain)

2011-07-15

The expansion of a dense plasma through a more rarefied ionized medium has been studied by means of two-dimensional particle-in-cell simulations. The initial conditions involve a density jump by a factor of 100, located in the middle of an otherwise equally dense electron-proton plasma with uniform proton and electron temperatures of 10 eV and 1 keV, respectively. Simulations show the creation of a purely electrostatic collisionless shock together with an ion-acoustic soliton tied to its downstream region. The shock front is seen to evolve in filamentary structures consistently with the onset of the ion-ion instability. Meanwhile, an un-magnetized drift instability is triggered in the core part of the dense plasma. Such results explain recent experimental laser-plasma experiments, carried out in similar conditions, and are of intrinsic relevance to non-relativistic shock scenarios in the solar and astrophysical systems.

Two-dimensional particle-in-cell simulation of the expansion of a plasma into a rarefied medium

International Nuclear Information System (INIS)

Sarri, G; Quinn, K; Kourakis, I; Borghesi, M; Murphy, G C; Drury, L O C; Dieckmann, M E; Ynnerman, A; Bret, A

2011-01-01

The expansion of a dense plasma through a more rarefied ionized medium has been studied by means of two-dimensional particle-in-cell simulations. The initial conditions involve a density jump by a factor of 100, located in the middle of an otherwise equally dense electron-proton plasma with uniform proton and electron temperatures of 10 eV and 1 keV, respectively. Simulations show the creation of a purely electrostatic collisionless shock together with an ion-acoustic soliton tied to its downstream region. The shock front is seen to evolve in filamentary structures consistently with the onset of the ion-ion instability. Meanwhile, an un-magnetized drift instability is triggered in the core part of the dense plasma. Such results explain recent experimental laser-plasma experiments, carried out in similar conditions, and are of intrinsic relevance to non-relativistic shock scenarios in the solar and astrophysical systems.

From particles to plasmas

International Nuclear Information System (INIS)

Van Dam, J.W.

1989-01-01

The title of this book, From Particles to Plasmas, has more than one meaning. First, it reflects how the scientific career of Marshall Rosenbluth has evolved, beginning in the field of elementary particle physics and extending into his major area of plasma physics. Secondly, it is meant to suggest the wide spectrum of subject matters addressed in the individual lectures, ranging from numerical simulation and space physics and accelerators to various subfields in the physics of plasmas. In the third place, the title is a reference to the way in which the theoretical description of plasmas is often constructed, namely starting from the motion of single particles and then incorporating collective effects. Most of the contributions in this book do concern various aspects of fusion plasma physics, which is the field in which most of Marshall Rosenbluth's scientific contributions have been and are being made. In this field his eminence and authority are indicated by the sobriquet pope of plasma physics that is often applied to him

Anomalous electron transport in Hall-effect thrusters: Comparison between quasi-linear kinetic theory and particle-in-cell simulations

Science.gov (United States)

Lafleur, T.; Martorelli, R.; Chabert, P.; Bourdon, A.

2018-06-01

Kinetic drift instabilities have been implicated as a possible mechanism leading to anomalous electron cross-field transport in E × B discharges, such as Hall-effect thrusters. Such instabilities, which are driven by the large disparity in electron and ion drift velocities, present a significant challenge to modelling efforts without resorting to time-consuming particle-in-cell (PIC) simulations. Here, we test aspects of quasi-linear kinetic theory with 2D PIC simulations with the aim of developing a self-consistent treatment of these instabilities. The specific quantities of interest are the instability growth rate (which determines the spatial and temporal evolution of the instability amplitude), and the instability-enhanced electron-ion friction force (which leads to "anomalous" electron transport). By using the self-consistently obtained electron distribution functions from the PIC simulations (which are in general non-Maxwellian), we find that the predictions of the quasi-linear kinetic theory are in good agreement with the simulation results. By contrast, the use of Maxwellian distributions leads to a growth rate and electron-ion friction force that is around 2-4 times higher, and consequently significantly overestimates the electron transport. A possible method for self-consistently modelling the distribution functions without requiring PIC simulations is discussed.

Plasma based charged-particle accelerators

International Nuclear Information System (INIS)

Bingham, R; Mendonca, J T; Shukla, P K

2004-01-01

Studies of charged-particle acceleration processes remain one of the most important areas of research in laboratory, space and astrophysical plasmas. In this paper, we present the underlying physics and the present status of high gradient and high energy plasma accelerators. We will focus on the acceleration of charged particles to relativistic energies by plasma waves that are created by intense laser and particle beams. The generation of relativistic plasma waves by intense lasers or electron beams in plasmas is important in the quest for producing ultra-high acceleration gradients for accelerators. With the development of compact short pulse high brightness lasers and electron positron beams, new areas of studies for laser/particle beam-matter interactions is opening up. A number of methods are being pursued vigorously to achieve ultra-high acceleration gradients. These include the plasma beat wave accelerator mechanism, which uses conventional long pulse (∼100 ps) modest intensity lasers (I ∼ 10 14 -10 16 W cm -2 ), the laser wakefield accelerator (LWFA), which uses the new breed of compact high brightness lasers ( 10 18 W cm -2 , the self-modulated LWFA concept, which combines elements of stimulated Raman forward scattering, and electron acceleration by nonlinear plasma waves excited by relativistic electron and positron bunches. In the ultra-high intensity regime, laser/particle beam-plasma interactions are highly nonlinear and relativistic, leading to new phenomena such as the plasma wakefield excitation for particle acceleration, relativistic self-focusing and guiding of laser beams, high-harmonic generation, acceleration of electrons, positrons, protons and photons. Fields greater than 1 GV cm -1 have been generated with particles being accelerated to 200 MeV over a distance of millimetre. Plasma wakefields driven by positron beams at the Stanford Linear Accelerator Center facility have accelerated the tail of the positron beam. In the near future

SELF-CONSISTENT LANGEVIN SIMULATION OF COULOMB COLLISIONS IN CHARGED-PARTICLE BEAMS

International Nuclear Information System (INIS)

QIANG, J.; RYNE, R.; HABIB, S.

2000-01-01

In many plasma physics and charged-particle beam dynamics problems, Coulomb collisions are modeled by a Fokker-Planck equation. In order to incorporate these collisions, we present a three-dimensional parallel Langevin simulation method using a Particle-In-Cell (PIC) approach implemented on high-performance parallel computers. We perform, for the first time, a fully self-consistent simulation, in which the FR-iction and diffusion coefficients are computed FR-om first principles. We employ a two-dimensional domain decomposition approach within a message passing programming paradigm along with dynamic load balancing. Object oriented programming is used to encapsulate details of the communication syntax as well as to enhance reusability and extensibility. Performance tests on the SGI Origin 2000 and the Cray T3E-900 have demonstrated good scalability. Work is in progress to apply our technique to intrabeam scattering in accelerators

Acceleration of particles in plasmas

CERN Multimedia

CERN. Geneva

2007-01-01

The accelerating fields in radio-frequency accelerators are limited to roughly 100 MV/m due to material breakdown which occurs on the walls of the structure. In contrast, a plasma, being already ionized, can support electric fields in excess of 100 GV/m. Such high accelerating gradients hold the promise of compact particle accelerators. Plasma acceleration has been an emerging and fast growing field of research in the past two decades. In this series of lectures, we will review the principles of plasma acceleration. We will see how relativistic plasma waves can be excited using an ultra-intense laser or using a particle beam. We will see how these plasma waves can be used to accelerate electrons to high energy in short distances. Throughout the lectures, we will also review recent experimental results. Current laser-plasma experiments throughout the world have shown that monoenergetic electron beams from 100 MeV to 1 GeV can be obtained in distances ranging from the millimetre to the centimetre. Experiments a...

Verification of high voltage rf capacitive sheath models with particle-in-cell simulations

Science.gov (United States)

Wang, Ying; Lieberman, Michael; Verboncoeur, John

2009-10-01

Collisionless and collisional high voltage rf capacitive sheath models were developed in the late 1980's [1]. Given the external parameters of a single-frequency capacitively coupled discharge, plasma parameters including sheath width, electron and ion temperature, plasma density, power, and ion bombarding energy can be estimated. One-dimensional electrostatic PIC codes XPDP1 [2] and OOPD1 [3] are used to investigate plasma behaviors within rf sheaths and bulk plasma. Electron-neutral collisions only are considered for collisionless sheaths, while ion-neutral collisions are taken into account for collisional sheaths. The collisionless sheath model is verified very well by PIC simulations for the rf current-driven and voltage-driven cases. Results will be reported for collisional sheaths also. [1] M. A. Lieberman, IEEE Trans. Plasma Sci. 16 (1988) 638; 17 (1989) 338 [2] J. P. Verboncoeur, M. V. Alves, V. Vahedi, and C. K. Birdsall, J. Comp. Phys. 104 (1993) 321 [3] J. P. Verboncoeur, A. B. Langdon and N. T. Gladd, Comp. Phys. Comm. 87 (1995) 199

Boltzmann electron PIC simulation of the E-sail effect

Directory of Open Access Journals (Sweden)

P. Janhunen

2015-12-01

Full Text Available The solar wind electric sail (E-sail is a planned in-space propulsion device that uses the natural solar wind momentum flux for spacecraft propulsion with the help of long, charged, centrifugally stretched tethers. The problem of accurately predicting the E-sail thrust is still somewhat open, however, due to a possible electron population trapped by the tether. Here we develop a new type of particle-in-cell (PIC simulation for predicting E-sail thrust. In the new simulation, electrons are modelled as a fluid, hence resembling hybrid simulation, but in contrast to normal hybrid simulation, the Poisson equation is used as in normal PIC to calculate the self-consistent electrostatic field. For electron-repulsive parts of the potential, the Boltzmann relation is used. For electron-attractive parts of the potential we employ a power law which contains a parameter that can be used to control the number of trapped electrons. We perform a set of runs varying the parameter and select the one with the smallest number of trapped electrons which still behaves in a physically meaningful way in the sense of producing not more than one solar wind ion deflection shock upstream of the tether. By this prescription we obtain thrust per tether length values that are in line with earlier estimates, although somewhat smaller. We conclude that the Boltzmann PIC simulation is a new tool for simulating the E-sail thrust. This tool enables us to calculate solutions rapidly and allows to easily study different scenarios for trapped electrons.

Advances in petascale kinetic plasma simulation with VPIC and Roadrunner

Energy Technology Data Exchange (ETDEWEB)

Bowers, Kevin J [Los Alamos National Laboratory; Albright, Brian J [Los Alamos National Laboratory; Yin, Lin [Los Alamos National Laboratory; Daughton, William S [Los Alamos National Laboratory; Roytershteyn, Vadim [Los Alamos National Laboratory; Kwan, Thomas J T [Los Alamos National Laboratory

2009-01-01

VPIC, a first-principles 3d electromagnetic charge-conserving relativistic kinetic particle-in-cell (PIC) code, was recently adapted to run on Los Alamos's Roadrunner, the first supercomputer to break a petaflop (10{sup 15} floating point operations per second) in the TOP500 supercomputer performance rankings. They give a brief overview of the modeling capabilities and optimization techniques used in VPIC and the computational characteristics of petascale supercomputers like Roadrunner. They then discuss three applications enabled by VPIC's unprecedented performance on Roadrunner: modeling laser plasma interaction in upcoming inertial confinement fusion experiments at the National Ignition Facility (NIF), modeling short pulse laser GeV ion acceleration and modeling reconnection in magnetic confinement fusion experiments.

Mixed-Language High-Performance Computing for Plasma Simulations

Directory of Open Access Journals (Sweden)

Quanming Lu

2003-01-01

Full Text Available Java is receiving increasing attention as the most popular platform for distributed computing. However, programmers are still reluctant to embrace Java as a tool for writing scientific and engineering applications due to its still noticeable performance drawbacks compared with other programming languages such as Fortran or C. In this paper, we present a hybrid Java/Fortran implementation of a parallel particle-in-cell (PIC algorithm for plasma simulations. In our approach, the time-consuming components of this application are designed and implemented as Fortran subroutines, while less calculation-intensive components usually involved in building the user interface are written in Java. The two types of software modules have been glued together using the Java native interface (JNI. Our mixed-language PIC code was tested and its performance compared with pure Java and Fortran versions of the same algorithm on a Sun E6500 SMP system and a Linux cluster of Pentium~III machines.

Convergence of lateral dynamic measurements in the plasma membrane of live cells from single particle tracking and STED-FCS

DEFF Research Database (Denmark)

Lagerholm, B. Christoffer; Andrade, Débora M.; Clausen, Mathias P.

2017-01-01

Fluorescence correlation spectroscopy (FCS) in combination with the super-resolution imaging method STED (STED-FCS), and single-particle tracking (SPT) are able to directly probe the lateral dynamics of lipids and proteins in the plasma membrane of live cells at spatial scales much below the diff...

AMITIS: A 3D GPU-Based Hybrid-PIC Model for Space and Plasma Physics

Science.gov (United States)

Fatemi, Shahab; Poppe, Andrew R.; Delory, Gregory T.; Farrell, William M.

2017-05-01

We have developed, for the first time, an advanced modeling infrastructure in space simulations (AMITIS) with an embedded three-dimensional self-consistent grid-based hybrid model of plasma (kinetic ions and fluid electrons) that runs entirely on graphics processing units (GPUs). The model uses NVIDIA GPUs and their associated parallel computing platform, CUDA, developed for general purpose processing on GPUs. The model uses a single CPU-GPU pair, where the CPU transfers data between the system and GPU memory, executes CUDA kernels, and writes simulation outputs on the disk. All computations, including moving particles, calculating macroscopic properties of particles on a grid, and solving hybrid model equations are processed on a single GPU. We explain various computing kernels within AMITIS and compare their performance with an already existing well-tested hybrid model of plasma that runs in parallel using multi-CPU platforms. We show that AMITIS runs ∼10 times faster than the parallel CPU-based hybrid model. We also introduce an implicit solver for computation of Faraday’s Equation, resulting in an explicit-implicit scheme for the hybrid model equation. We show that the proposed scheme is stable and accurate. We examine the AMITIS energy conservation and show that the energy is conserved with an error < 0.2% after 500,000 timesteps, even when a very low number of particles per cell is used.

Particle-in-cell simulation of Trichel pulses in pure oxygen

International Nuclear Information System (INIS)

Soria-Hoyo, C; Pontiga, F; Castellanos, A

2007-01-01

The development and propagation of Trichel pulses in oxygen have been numerically simulated using an improved fluid particle-in-cell (PIC) method. The numerical method has been optimized to accurately reproduce sequences of about 100 Trichel pulses (∼1 ms). A classical one-dimensional model of negative corona in sphere-to-plane geometry has been used to formulate the continuity equations for electrons and ions. The effects of ionization, attachment and secondary-electron emission from the cathode have all been considered. The electric field has been obtained from the solution of Poisson's equation in two dimensions. Using this model, the temporal and electrical characteristics of Trichel pulses have been investigated, in particular, the relation between applied voltage, pulse frequency and time-averaged current intensity and charge

A survey of elementary plasma instabilities and ECH wave noise properties relevant to plasma sounding by means of particle in cell simulations

International Nuclear Information System (INIS)

Dieckmann, M.E.

1999-01-01

In this work the emission of high amplitude wave packets into a plasma is examined. The plasma is modelled by an 1 1/2D electromagnetic and relativistic particle in cell code. The antenna is modelled by applying forced electrostatic field oscillations to a subset of the simulation grid cells. The emitted wave packets are followed in space and time. It is investigated how the wave packets are affected by instabilities. The detected instabilities affecting ECH waves have been identified as wave decay, nonlinear damping due to trapping and modulational instabilities. These instabilities have been discussed with hindsight to the plasma sounding experiment. A plasma sounder is an experiment emitting short wave packets into the ambient plasma and then it listens to the response. The assumption that the emitted waves are linear waves then allows to determine the plasma magnetic field strength, the electron density and possibly the electron thermal velocity from the response spectrum. The impact of the non-linear instabilities on the plasma wave response spectrum provided by a sounder have been predicted in this work and the predictions have been shown to match a wide range of experimental observations. A dependence of the instabilities on the simulation noise levels, for example the dependence of the wave interaction time in a wave decay on the noise electric field amplitudes, required it to investigate the simulation noise properties (spectral distribution) and to compare it to real plasma thermal noise. It has also been examined how a finite length antenna would filter the simulation noise. (author)

Research of transportation efficiency of low-energy high- current electron beam in plasma channel in external magnetic field

International Nuclear Information System (INIS)

Vagin, E S; Grigoriev, V P

2015-01-01

Effective high current (5-20 kA) and low energy (tens of keV) electrons beam transportation is possible only with almost complete charging neutralization. It is also necessary to use quite high current neutralization for elimination beam self-pinching effect. The research is based on the self-consistent mathematical model that takes into account beam and plasma particles dynamic, current and charge neutralization of electron beam and examines the transportation of electron beam into a chamber with low-pressure plasma in magnetic field. A numerical study was conducted using particle in cell (PIC) method. The study was performed with various system parameters: rise time and magnitude of the beam current, gas pressure and plasma density and geometry of the system. Regularities of local virtual cathode field generated by the beam in the plasma channel, as well as ranges of parameters that let transportation beam with minimal losses, depending on the external magnetic field were determined through a series of numerical studies. In addition, the assessment of the impact of the plasma ion mobility during the transition period and during steady beam was performed. (paper)

Particle balance in diverted plasmas in TEXT-U

International Nuclear Information System (INIS)

Rowan, W.L.; Bengtson, R.D.; Bonnin, X.; Edmonds, P.H.; Hurwitz, P.D.; Solano, E.R.; Tsui, H.Y.W.; Uglum, J.R.; Wootton, A.J.

1995-01-01

Tokamak plasmas in an open divertor configuration with the X-point at the inner equator were compared with plasmas limited by an inner toroidal belt. This paper describes and compares the particle balance in these two types of discharge. The plasma parameters n e (r), φ(r), T e (r), n e , and φ were measured in the scrape-off layer and in the plasma periphery at one poloidal location and mapped onto the rest of the plasma by assuming constancy on flux surfaces. Emission from neutral hydrogen was measured throughout the plasma. The particle source and then the global particle confinement were inferred from these measurements using a 3-D neutral transport simulation. The SOL profiles are significantly steeper in the diverted discharge. Both the source measurements and estimates from the SOL profiles indicate significantly better confinement in the diverted discharge. ((orig.))

Three-dimensional relativistic pair plasma reconnection with radiative feedback in the Crab Nebula

Energy Technology Data Exchange (ETDEWEB)

Cerutti, B. [Department of Astrophysical Sciences, Princeton University, Princeton, NJ 08544 (United States); Werner, G. R.; Uzdensky, D. A. [Center for Integrated Plasma Studies, Physics Department, University of Colorado, UCB 390, Boulder, CO 80309-0390 (United States); Begelman, M. C., E-mail: bcerutti@astro.princeton.edu, E-mail: greg.werner@colorado.edu, E-mail: uzdensky@colorado.edu, E-mail: mitch@jila.colorado.edu [JILA, University of Colorado and National Institute of Standards and Technology, UCB 440, Boulder, CO 80309-0440 (United States)

2014-02-20

The discovery of rapid synchrotron gamma-ray flares above 100 MeV from the Crab Nebula has attracted new interest in alternative particle acceleration mechanisms in pulsar wind nebulae. Diffuse shock-acceleration fails to explain the flares because particle acceleration and emission occur during a single or even sub-Larmor timescale. In this regime, the synchrotron energy losses induce a drag force on the particle motion that balances the electric acceleration and prevents the emission of synchrotron radiation above 160 MeV. Previous analytical studies and two-dimensional (2D) particle-in-cell (PIC) simulations indicate that relativistic reconnection is a viable mechanism to circumvent the above difficulties. The reconnection electric field localized at X-points linearly accelerates particles with little radiative energy losses. In this paper, we check whether this mechanism survives in three dimension (3D), using a set of large PIC simulations with radiation reaction force and with a guide field. In agreement with earlier works, we find that the relativistic drift kink instability deforms and then disrupts the layer, resulting in significant plasma heating but few non-thermal particles. A moderate guide field stabilizes the layer and enables particle acceleration. We report that 3D magnetic reconnection can accelerate particles above the standard radiation reaction limit, although the effect is less pronounced than in 2D with no guide field. We confirm that the highest-energy particles form compact bunches within magnetic flux ropes, and a beam tightly confined within the reconnection layer, which could result in the observed Crab flares when, by chance, the beam crosses our line of sight.

PIC simulation of a thermal anisotropy-driven Weibel instability in a circular rarefaction wave

International Nuclear Information System (INIS)

Dieckmann, M E; Sarri, G; Kourakis, I; Borghesi, M; Murphy, G C; O'C Drury, L; Bret, A; Romagnani, L; Ynnerman, A

2012-01-01

The expansion of an initially unmagnetized planar rarefaction wave has recently been shown to trigger a thermal anisotropy-driven Weibel instability (TAWI), which can generate magnetic fields from noise levels. It is examined here whether the TAWI can also grow in a curved rarefaction wave. The expansion of an initially unmagnetized circular plasma cloud, which consists of protons and hot electrons, into a vacuum is modelled for this purpose with a two-dimensional particle-in-cell (PIC) simulation. It is shown that the momentum transfer from the electrons to the radially accelerating protons can indeed trigger a TAWI. Radial current channels form and the aperiodic growth of a magnetowave is observed, which has a magnetic field that is oriented orthogonal to the simulation plane. The induced electric field implies that the electron density gradient is no longer parallel to the electric field. Evidence is presented here that this electric field modification triggers a second magnetic instability, which results in a rotational low-frequency magnetowave. The relevance of the TAWI is discussed for the growth of small-scale magnetic fields in astrophysical environments, which are needed to explain the electromagnetic emissions by astrophysical jets. It is outlined how this instability could be examined experimentally. (paper)

PIC simulation of a thermal anisotropy-driven Weibel instability in a circular rarefaction wave

Science.gov (United States)

Dieckmann, M. E.; Sarri, G.; Murphy, G. C.; Bret, A.; Romagnani, L.; Kourakis, I.; Borghesi, M.; Ynnerman, A.; O'C Drury, L.

2012-02-01

The expansion of an initially unmagnetized planar rarefaction wave has recently been shown to trigger a thermal anisotropy-driven Weibel instability (TAWI), which can generate magnetic fields from noise levels. It is examined here whether the TAWI can also grow in a curved rarefaction wave. The expansion of an initially unmagnetized circular plasma cloud, which consists of protons and hot electrons, into a vacuum is modelled for this purpose with a two-dimensional particle-in-cell (PIC) simulation. It is shown that the momentum transfer from the electrons to the radially accelerating protons can indeed trigger a TAWI. Radial current channels form and the aperiodic growth of a magnetowave is observed, which has a magnetic field that is oriented orthogonal to the simulation plane. The induced electric field implies that the electron density gradient is no longer parallel to the electric field. Evidence is presented here that this electric field modification triggers a second magnetic instability, which results in a rotational low-frequency magnetowave. The relevance of the TAWI is discussed for the growth of small-scale magnetic fields in astrophysical environments, which are needed to explain the electromagnetic emissions by astrophysical jets. It is outlined how this instability could be examined experimentally.

New relativistic particle-in-cell simulation studies of prompt and early afterglows from GRBs

International Nuclear Information System (INIS)

Ken-Ichi Nishikawa

2008-01-01

Nonthermal radiation observed from astrophysical systems containing relativistic jets and shocks, e.g., gamma-ray bursts (GRBs), active galactic nuclei (AGNs), and Galactic microquasar systems usually have power-law emission spectra. Recent PIC simulations of relativistic electron-ion (electro-positron) jets injected into a stationary medium show that particle acceleration occurs within the downstream jet. In the collisionless relativistic shock particle acceleration is due to plasma waves and their associated instabilities (e.g., the Buneman instability, other two-streaming instability, and the Weibel (filamentation) instability) created in the shocks are responsible for particle (electron, positron, and ion) acceleration. The simulation results show that the Weibel instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields. These magnetic fields contribute to the electrons' transverse deflection behind the jet head. The '' jitter '' radiation from deflected electrons has different properties than synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation may be important to understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants. (author)

Plasma physics via particle simulation

International Nuclear Information System (INIS)

Birdsall, C.K.

1981-01-01

Plasmas are studied by following the motion of many particles in applied and self fields, analytically, experimentally and computationally. Plasmas for magnetic fusion energy devices are very hot, nearly collisionless and magnetized, with scale lengths of many ion gyroradii and Debye lengths. The analytic studies of such plasmas are very difficult as the plasma is nonuniform, anisotropic and nonlinear. The experimental studies have become very expensive in time and money, as the size, density and temperature approach fusion reactor values. Computational studies using many particles and/or fluids have complemented both theories and experiments for many years and have progressed to fully three dimensional electromagnetic models, albeit with hours of running times on the fastest largest computers. Particle simulation methods are presented in some detail, showing particle advance from acceleration to velocity to position, followed by calculation of the fields from charge and current densities and then further particle advance, and so on. Limitations due to the time stepping and use of a spatial grid are given, to avoid inaccuracies and instabilities. Examples are given for an electrostatic program in one dimension of an orbit averaging program, and for a three dimensional electromagnetic program. Applications of particle simulations of plasmas in magnetic and inertial fusion devices continue to grow, as well as to plasmas and beams in peripheral devices, such as sources, accelerators, and converters. (orig.)

Introducing a distributed unstructured mesh into gyrokinetic particle-in-cell code, XGC

Science.gov (United States)

Yoon, Eisung; Shephard, Mark; Seol, E. Seegyoung; Kalyanaraman, Kaushik

2017-10-01

XGC has shown good scalability for large leadership supercomputers. The current production version uses a copy of the entire unstructured finite element mesh on every MPI rank. Although an obvious scalability issue if the mesh sizes are to be dramatically increased, the current approach is also not optimal with respect to data locality of particles and mesh information. To address these issues we have initiated the development of a distributed mesh PIC method. This approach directly addresses the base scalability issue with respect to mesh size and, through the use of a mesh entity centric view of the particle mesh relationship, provides opportunities to address data locality needs of many core and GPU supported heterogeneous systems. The parallel mesh PIC capabilities are being built on the Parallel Unstructured Mesh Infrastructure (PUMI). The presentation will first overview the form of mesh distribution used and indicate the structures and functions used to support the mesh, the particles and their interaction. Attention will then focus on the node-level optimizations being carried out to ensure performant operation of all PIC operations on the distributed mesh. Partnership for Edge Physics Simulation (EPSI) Grant No. DE-SC0008449 and Center for Extended Magnetohydrodynamic Modeling (CEMM) Grant No. DE-SC0006618.

Generation of electromagnetic waves and Alfven waves during coalescence of magnetic islands in pair plasmas

International Nuclear Information System (INIS)

Sakai, J.I.; Haruki, T.; Kazimura, Y.

2000-01-01

It is shown by using a 2-D fully relativistic electromagnetic particle-in-cell (PIC) code that the tearing instability in a current sheet of pair plasmas is caused by Landau resonances of both electrons and positrons. Strong magnetic flux can be generated during coalescence of magnetic islands in the nonlinear phase of the tearing instability. The magnetic flux produced in an O-type magnetic island is caused from the counter-streaming instability found by Kazimura et al. (1998). It is also shown that charge separation with a quadrupole-like structure is generated from the localized strong magnetic flux. During the decay of the quadrupole-like charge structure as well as the magnetic flux, there appear wave emission with high-frequency electromagnetic waves and Alfven waves as well as Langmuir waves. We also show by using a 3-D PIC code that current filaments associated with the O-type magnetic islands become unstable against the kink instability during the coalescence of current filaments. (orig.)

Computationally efficient description of relativistic electron beam transport in dense plasma

Science.gov (United States)

Polomarov, Oleg; Sefkov, Adam; Kaganovich, Igor; Shvets, Gennady

2006-10-01

A reduced model of the Weibel instability and electron beam transport in dense plasma is developed. Beam electrons are modeled by macro-particles and the background plasma is represented by electron fluid. Conservation of generalized vorticity and quasineutrality of the plasma-beam system are used to simplify the governing equations. Our approach is motivated by the conditions of the FI scenario, where the beam density is likely to be much smaller than the plasma density and the beam energy is likely to be very high. For this case the growth rate of the Weibel instability is small, making the modeling of it by conventional PICs exceedingly time consuming. The present approach does not require resolving the plasma period and only resolves a plasma collisionless skin depth and is suitable for modeling a long-time behavior of beam-plasma interaction. An efficient code based on this reduced description is developed and benchmarked against the LSP PIC code. The dynamics of low and high current electron beams in dense plasma is simulated. Special emphasis is on peculiarities of its non-linear stages, such as filament formation and merger, saturation and post-saturation field and energy oscillations. *Supported by DOE Fusion Science through grant DE-FG02-05ER54840.

INTERACTIONS OF FLOWING PLASMA WITH COLLECTING SPHERE/OBJECTS

Energy Technology Data Exchange (ETDEWEB)

HUTCHINSON, IAN H.

2013-10-31

The proposer’s computational flowing plasma group, supported in part by NSF/DOE since 2006, has been enormously productive. It has published 22 journal articles, of which 3 are in Physical Review Letters; received a total of 330 citations, and made 5 invited and numerous contributed presentations at international conferences. In addition, other groups have published at least 3 papers [20-22] based upon results they have obtained from the SCEPTIC code, developed by the proposer. The codes developed for this subject are Particle In Cell (PIC) codes, representing the entire ion distribution function and its variation in space in the vicinity of an absorbing object.

Load-balancing techniques for a parallel electromagnetic particle-in-cell code

Energy Technology Data Exchange (ETDEWEB)

PLIMPTON,STEVEN J.; SEIDEL,DAVID B.; PASIK,MICHAEL F.; COATS,REBECCA S.

2000-01-01

QUICKSILVER is a 3-d electromagnetic particle-in-cell simulation code developed and used at Sandia to model relativistic charged particle transport. It models the time-response of electromagnetic fields and low-density-plasmas in a self-consistent manner: the fields push the plasma particles and the plasma current modifies the fields. Through an LDRD project a new parallel version of QUICKSILVER was created to enable large-scale plasma simulations to be run on massively-parallel distributed-memory supercomputers with thousands of processors, such as the Intel Tflops and DEC CPlant machines at Sandia. The new parallel code implements nearly all the features of the original serial QUICKSILVER and can be run on any platform which supports the message-passing interface (MPI) standard as well as on single-processor workstations. This report describes basic strategies useful for parallelizing and load-balancing particle-in-cell codes, outlines the parallel algorithms used in this implementation, and provides a summary of the modifications made to QUICKSILVER. It also highlights a series of benchmark simulations which have been run with the new code that illustrate its performance and parallel efficiency. These calculations have up to a billion grid cells and particles and were run on thousands of processors. This report also serves as a user manual for people wishing to run parallel QUICKSILVER.

Load-balancing techniques for a parallel electromagnetic particle-in-cell code

International Nuclear Information System (INIS)

Plimpton, Steven J.; Seidel, David B.; Pasik, Michael F.; Coats, Rebecca S.

2000-01-01

QUICKSILVER is a 3-d electromagnetic particle-in-cell simulation code developed and used at Sandia to model relativistic charged particle transport. It models the time-response of electromagnetic fields and low-density-plasmas in a self-consistent manner: the fields push the plasma particles and the plasma current modifies the fields. Through an LDRD project a new parallel version of QUICKSILVER was created to enable large-scale plasma simulations to be run on massively-parallel distributed-memory supercomputers with thousands of processors, such as the Intel Tflops and DEC CPlant machines at Sandia. The new parallel code implements nearly all the features of the original serial QUICKSILVER and can be run on any platform which supports the message-passing interface (MPI) standard as well as on single-processor workstations. This report describes basic strategies useful for parallelizing and load-balancing particle-in-cell codes, outlines the parallel algorithms used in this implementation, and provides a summary of the modifications made to QUICKSILVER. It also highlights a series of benchmark simulations which have been run with the new code that illustrate its performance and parallel efficiency. These calculations have up to a billion grid cells and particles and were run on thousands of processors. This report also serves as a user manual for people wishing to run parallel QUICKSILVER

Simulation of QED effects in ultrahigh intensity laser-plasma interaction

International Nuclear Information System (INIS)

Kostyukov, I.; Nerush, E.

2010-01-01

Complete text of publication follows. Due to an impressive progress in laser technology, laser pulses with peak intensity of nearly 2 x 10 22 W/cm 2 are now available in laboratory. When the matter is irradiated by so intense laser pulses high energy density plasma is produced. Besides of fundamental interest such plasma is the efficient source of particles and radiation with extreme parameters that opens bright perspectives in developments of advanced particle accelerators, next generation of radiation sources, laboratory modelling of astrophysics phenomena etc. Even high laser intensity the radiation reaction and QED effects become important. One of the QED effects, which recently attracts much attention, is the electron-positron plasma creation in strong laser field. The plasma can be produced via electromagnetic cascades: the seeded charged particles is accelerated in the field of counter-propagating laser pulses, then they emit energetic photons, the photons by turn decay in the laser field and create electron-positron pairs. The pair particles accelerated in the laser field produce new generation of the photons and pairs. For self-consistent study of the electron-positron plasma dynamics in the laser field we develop 2D code based on particle-in-cell and Monte-Carlo methods. The electron, positron and photon dynamics as well as evolution of the plasma and laser fields are calculated by PIC technique while photon emission and pair production are calculated by Monte-Carlo method. We simulate pair production in the field of counter-propagating linearly polarized laser pulses. It is shown that for the laser intensity above threshold the plasma production becomes so intense that the laser pulse are strongly absorbed in the plasma. The laser intensity threshold and the rate of laser field absorption are calculated. Acknowledgements. This work has been supported by federal target 'The scientific and scientific-pedagogical personnel of innovation in Russia' and by

Simulation of high-energy particle production through sausage and kink instabilities in pinched plasma discharges

International Nuclear Information System (INIS)

Haruki, Takayuki; Yousefi, Hamid Reza; Masugata, Katsumi; Sakai, Jun-Ichi; Mizuguchi, Yusuke; Makino, Nao; Ito, Hiroaki

2006-01-01

In an experimental plasma, high-energy particles were observed by using a plasma focus device, to obtain energies of a few hundred keV for electrons, up to MeV for ions. In order to study the mechanism of high-energy particle production in pinched plasma discharges, a numerical simulation was introduced. By use of a three-dimensional relativistic and fully electromagnetic particle-in-cell code, the dynamics of a Z-pinch plasma, thought to be unstable against sausage and kink instabilities, are investigated. In this work, the development of sausage and kink instabilities and subsequent high-energy particle production are shown. In the model used here, cylindrically distributed electrons and ions are driven by an external electric field. The driven particles spontaneously produce a current, which begins to pinch by the Lorentz force. Initially the pinched current is unstable against a sausage instability, and then becomes unstable against a kink instability. As a result high-energy particles are observed

A PIC-MCC code RFdinity1d for simulation of discharge initiation by ICRF antenna

Science.gov (United States)

Tripský, M.; Wauters, T.; Lyssoivan, A.; Bobkov, V.; Schneider, P. A.; Stepanov, I.; Douai, D.; Van Eester, D.; Noterdaeme, J.-M.; Van Schoor, M.; ASDEX Upgrade Team; EUROfusion MST1 Team

2017-12-01

Discharges produced and sustained by ion cyclotron range of frequency (ICRF) waves in absence of plasma current will be used on ITER for (ion cyclotron-) wall conditioning (ICWC, Te = 3{-}5 eV, ne 18 m-3 ). In this paper, we present the 1D particle-in-cell Monte Carlo collision (PIC-MCC) RFdinity1d for the study the breakdown phase of ICRF discharges, and its dependency on the RF discharge parameters (i) antenna input power P i , (ii) RF frequency f, (iii) shape of the electric field and (iv) the neutral gas pressure pH_2 . The code traces the motion of both electrons and ions in a narrow bundle of magnetic field lines close to the antenna straps. The charged particles are accelerated in the parallel direction with respect to the magnetic field B T by two electric fields: (i) the vacuum RF field of the ICRF antenna E_z^RF and (ii) the electrostatic field E_zP determined by the solution of Poisson’s equation. The electron density evolution in simulations follows exponential increase, {\\dot{n_e} ∼ ν_ion t } . The ionization rate varies with increasing electron density as different mechanisms become important. The charged particles are affected solely by the antenna RF field E_z^RF at low electron density ({ne < 1011} m-3 , {≤ft \\vert E_z^RF \\right \\vert \\gg ≤ft \\vert E_zP \\right \\vert } ). At higher densities, when the electrostatic field E_zP is comparable to the antenna RF field E_z^RF , the ionization frequency reaches the maximum. Plasma oscillations propagating toroidally away from the antenna are observed. The simulated energy distributions of ions and electrons at {ne ∼ 1015} m-3 correspond a power-law Kappa energy distribution. This energy distribution was also observed in NPA measurements at ASDEX Upgrade in ICWC experiments.

Loading relativistic Maxwell distributions in particle simulations

Energy Technology Data Exchange (ETDEWEB)

Zenitani, Seiji, E-mail: seiji.zenitani@nao.ac.jp [National Astronomical Observatory of Japan, 2-21-1 Osawa, Mitaka, Tokyo 181-8588 (Japan)

2015-04-15

Numerical algorithms to load relativistic Maxwell distributions in particle-in-cell (PIC) and Monte-Carlo simulations are presented. For stationary relativistic Maxwellian, the inverse transform method and the Sobol algorithm are reviewed. To boost particles to obtain relativistic shifted-Maxwellian, two rejection methods are proposed in a physically transparent manner. Their acceptance efficiencies are ≈50% for generic cases and 100% for symmetric distributions. They can be combined with arbitrary base algorithms.

Loading relativistic Maxwell distributions in particle simulations

International Nuclear Information System (INIS)

Zenitani, Seiji

2015-01-01

Numerical algorithms to load relativistic Maxwell distributions in particle-in-cell (PIC) and Monte-Carlo simulations are presented. For stationary relativistic Maxwellian, the inverse transform method and the Sobol algorithm are reviewed. To boost particles to obtain relativistic shifted-Maxwellian, two rejection methods are proposed in a physically transparent manner. Their acceptance efficiencies are ≈50% for generic cases and 100% for symmetric distributions. They can be combined with arbitrary base algorithms

Particle-in-cell simulation of Trichel pulses in pure oxygen

Energy Technology Data Exchange (ETDEWEB)

Soria-Hoyo, C [Dpto. Electronica y Electromagnetismo, Universidad de Sevilla, Avda. Reina Mercedes s/n, Sevilla 41012 (Spain); Pontiga, F [Dpto. Fisica Aplicada II, Universidad de Sevilla, Avda. Reina Mercedes s/n, Sevilla 41012 (Spain); Castellanos, A [Dpto. Electronica y Electromagnetismo, Universidad de Sevilla, Avda. Reina Mercedes s/n, Sevilla 41012 (Spain)

2007-08-07

The development and propagation of Trichel pulses in oxygen have been numerically simulated using an improved fluid particle-in-cell (PIC) method. The numerical method has been optimized to accurately reproduce sequences of about 100 Trichel pulses ({approx}1 ms). A classical one-dimensional model of negative corona in sphere-to-plane geometry has been used to formulate the continuity equations for electrons and ions. The effects of ionization, attachment and secondary-electron emission from the cathode have all been considered. The electric field has been obtained from the solution of Poisson's equation in two dimensions. Using this model, the temporal and electrical characteristics of Trichel pulses have been investigated, in particular, the relation between applied voltage, pulse frequency and time-averaged current intensity and charge.

Fly ash particles spheroidization using low temperature plasma energy

Science.gov (United States)

Shekhovtsov, V. V.; Volokitin, O. G.; Kondratyuk, A. A.; Vitske, R. E.

2016-11-01

The paper presents the investigations on producing spherical particles 65-110 μm in size using the energy of low temperature plasma (LTP). These particles are based on flow ash produced by the thermal power plant in Seversk, Tomsk region, Russia. The obtained spherical particles have no defects and are characterized by a smooth exterior surface. The test bench is designed to produce these particles. With due regard for plasma temperature field distribution, it is shown that the transition of fly ash particles to a state of viscous flow occurs at 20 mm distance from the plasma jet. The X-ray phase analysis is carried out for the both original state of fly ash powders and the particles obtained. This analysis shows that fly ash contains 56.23 wt.% SiO2; 20.61 wt.% Al2O3 and 17.55 wt.% Fe2O3 phases that mostly contribute to the integral (experimental) intensity of the diffraction maximum. The LTP treatment results in a complex redistribution of the amorphous phase amount in the obtained spherical particles, including the reduction of O2Si, phase, increase of O22Al20 and Fe2O3 phases and change in Al, O density of O22Al20 chemical unit cell.

Wavelet-Based Poisson Solver for Use in Particle-in-Cell Simulations

CERN Document Server

Terzic, Balsa; Mihalcea, Daniel; Pogorelov, Ilya V

2005-01-01

We report on a successful implementation of a wavelet-based Poisson solver for use in 3D particle-in-cell simulations. One new aspect of our algorithm is its ability to treat the general (inhomogeneous) Dirichlet boundary conditions. The solver harnesses advantages afforded by the wavelet formulation, such as sparsity of operators and data sets, existence of effective preconditioners, and the ability simultaneously to remove numerical noise and further compress relevant data sets. Having tested our method as a stand-alone solver on two model problems, we merged it into IMPACT-T to obtain a fully functional serial PIC code. We present and discuss preliminary results of application of the new code to the modelling of the Fermilab/NICADD and AES/JLab photoinjectors.

Two dimensional numerical simulation of gas discharges: comparison between particle-in-cell and FCT techniques

Energy Technology Data Exchange (ETDEWEB)

Soria-Hoyo, C; Castellanos, A [Departamento de Electronica y Electromagnetismo, Facultad de Fisica, Universidad de Sevilla, Avda. Reina Mercedes s/n, 41012 Sevilla (Spain); Pontiga, F [Departamento de Fisica Aplicada II, EUAT, Universidad de Sevilla, Avda. Reina Mercedes s/n, 41012 Sevilla (Spain)], E-mail: cshoyo@us.es

2008-10-21

Two different numerical techniques have been applied to the numerical integration of equations modelling gas discharges: a finite-difference flux corrected transport (FD-FCT) technique and a particle-in-cell (PIC) technique. The PIC technique here implemented has been specifically designed for the simulation of 2D electrical discharges using cylindrical coordinates. The development and propagation of a streamer between two parallel electrodes has been used as a convenient test to compare the performance of both techniques. In particular, the phase velocity of the cathode directed streamer has been used to check the internal consistency of the numerical simulations. The results obtained from the two techniques are in reasonable agreement with each other, and both techniques have proved their ability to follow the high gradients of charge density and electric field present in this type of problems. Moreover, the streamer velocities predicted by the simulation are in accordance with the typical experimental values.

Two dimensional numerical simulation of gas discharges: comparison between particle-in-cell and FCT techniques

International Nuclear Information System (INIS)

Soria-Hoyo, C; Castellanos, A; Pontiga, F

2008-01-01

Two different numerical techniques have been applied to the numerical integration of equations modelling gas discharges: a finite-difference flux corrected transport (FD-FCT) technique and a particle-in-cell (PIC) technique. The PIC technique here implemented has been specifically designed for the simulation of 2D electrical discharges using cylindrical coordinates. The development and propagation of a streamer between two parallel electrodes has been used as a convenient test to compare the performance of both techniques. In particular, the phase velocity of the cathode directed streamer has been used to check the internal consistency of the numerical simulations. The results obtained from the two techniques are in reasonable agreement with each other, and both techniques have proved their ability to follow the high gradients of charge density and electric field present in this type of problems. Moreover, the streamer velocities predicted by the simulation are in accordance with the typical experimental values.

Charged particle and photon acceleration by wakefield plasma waves in non-uniform plasmas

International Nuclear Information System (INIS)

Bulanov, S.V.; Kirsanov, V.I.; Sakharov, A.S.; Pegoraro, F.

1993-01-01

We discuss the acceleration of charged particles and the upshift of the frequency of short wave packets of laser radiation. The acceleration and the upshift are caused by wake plasma waves excited by a strong laser pulse in a non-uniform plasma. We show that unlimited acceleration of charged particles is possible for specific spatial dependencies of the plasma density. In this unlimited acceleration regime, particles have a fixed phase relationship with respect to the plasma wave, while their energy increases with time. When the wave breaking limit is approached and surpassed, the efficiency of the acceleration of the charged particles and of the frequency upshift of the photons can be increased significantly. (author) 3 refs

Particle behavior in an ECR plasma etch tool

International Nuclear Information System (INIS)

Blain, M.G.; Tipton, G.D.; Holber, W.M.; Westerfield, P.L.; Maxwell, K.L.

1993-01-01

Sources of particles in a close-coupled electron cyclotron resonance (ECR) polysilicon plasma etch source include flaking of films deposited on chamber surfaces, and shedding of material from electrostatic wafer chucks. A large, episodic increase in the number of particles added to a wafer in a clean system is observed more frequently for a plasma-on than for a gas-only source condition. For polymer forming process conditions, particles were added to wafers by a polymer film which was observed to fracture and flake away from chamber surfaces. The presence of a plasma, especially when rf bias is applied to the wafer, caused more particles to be ejected from the walls and added to wafers than the gas-only condition; however, no significant influence was observed with different microwave powers. A study of effect of electrode temperatures on particles added showed that thermophoretic forces are not significant for this ECR configuration. Particles originating from the electrostatic chuck were observed to be deposited on wafers in much larger numbers in the presence of the plasma as compared to gas-only conditions

Rise time of proton cut-off energy in 2D and 3D PIC simulations

Science.gov (United States)

Babaei, J.; Gizzi, L. A.; Londrillo, P.; Mirzanejad, S.; Rovelli, T.; Sinigardi, S.; Turchetti, G.

2017-04-01

The Target Normal Sheath Acceleration regime for proton acceleration by laser pulses is experimentally consolidated and fairly well understood. However, uncertainties remain in the analysis of particle-in-cell simulation results. The energy spectrum is exponential with a cut-off, but the maximum energy depends on the simulation time, following different laws in two and three dimensional (2D, 3D) PIC simulations so that the determination of an asymptotic value has some arbitrariness. We propose two empirical laws for the rise time of the cut-off energy in 2D and 3D PIC simulations, suggested by a model in which the proton acceleration is due to a surface charge distribution on the target rear side. The kinetic energy of the protons that we obtain follows two distinct laws, which appear to be nicely satisfied by PIC simulations, for a model target given by a uniform foil plus a contaminant layer that is hydrogen-rich. The laws depend on two parameters: the scaling time, at which the energy starts to rise, and the asymptotic cut-off energy. The values of the cut-off energy, obtained by fitting 2D and 3D simulations for the same target and laser pulse configuration, are comparable. This suggests that parametric scans can be performed with 2D simulations since 3D ones are computationally very expensive, delegating their role only to a correspondence check. In this paper, the simulations are carried out with the PIC code ALaDyn by changing the target thickness L and the incidence angle α, with a fixed a0 = 3. A monotonic dependence, on L for normal incidence and on α for fixed L, is found, as in the experimental results for high temporal contrast pulses.

Investigation of the Electron Acceleration by a High-Power Laser and a Density-Tapered Mixed-Gas Cell

Science.gov (United States)

Kim, Jinju; Phung, Vanessa L. J.; Kim, Minseok; Hur, Min-Sup; Suk, Hyyong

2017-10-01

Plasma-based accelerators can generate about 1000 times stronger acceleration field compared with RF-based conventional accelerators, which can be done by high power laser and plasma. There are many issues in this research and one of them is development of a good plasma source for higher electron beam energy. For this purpose, we are investigating a special type of plasma source, which is a density-tapered gas cell with a mixed-gas for easy injection. By this type of special gas cell, we expect higher electron beam energies with easy injection in the wakefield. In this poster, some experimental results for electron beam generation with the density-tapered mixed-gas cell are presented. In addition to the experimental results, CFD (Computational-Fluid-Dynamics) and PIC (Particle-In-Cell) simulation results are also presented for comparison studies.

Modelling of neutral particle transport in divertor plasma

International Nuclear Information System (INIS)

Kakizuka, Tomonori; Shimizu, Katsuhiro

1995-01-01

An outline of the modelling of neutral particle transport in the diverter plasma was described in the paper. The characteristic properties of divertor plasma were largely affected by interaction between neutral particles and divertor plasma. Accordingly, the behavior of neutral particle should be investigated quantitatively. Moreover, plasma and neutral gas should be traced consistently in the plasma simulation. There are Monte Carlo modelling and the neutral gas fluid modelling as the transport modelling. The former need long calculation time, but it is able to make the physical process modelling. A ultra-large parallel computer is good for the former. In spite of proposing some kinds of models, the latter has not been established. At the view point of reducing calculation time, a work station is good for the simulation of the latter, although some physical problems have not been solved. On the Monte Carlo method particle modelling, reducing the calculation time and introducing the interaction of particles are important subjects to develop 'the evolutional Monte Carlo Method'. To reduce the calculation time, two new methods: 'Implicit Monte Carlo method' and 'Free-and Diffusive-Motion Hybrid Monte-Carlo method' have been developing. (S.Y.)

Plasmoid statistics in relativistic magnetic reconnection

Science.gov (United States)

Petropoulou, M.; Christie, I. M.; Sironi, L.; Giannios, D.

2018-04-01

Plasmoids, overdense blobs of plasma containing magnetic fields and high-energy particles, are a self-consistent outcome of the reconnection process in the relativistic regime. Recent two-dimensional particle-in-cell (PIC) simulations have shown that plasmoids can undergo a variety of processes (e.g. mergers, bulk acceleration, growth, and advection) within the reconnection layer. We developed a Monte Carlo code, benchmarked with the recent PIC simulations, to examine the effects of these processes on the steady-state size and momentum distributions of the plasmoid chain. The differential plasmoid size distribution is shown to be a power law, ranging from a few plasma skin depths to ˜0.1 of the reconnection layer's length. The power-law slope is shown to be linearly dependent upon the ratio of the plasmoid acceleration and growth rates, which slightly decreases with increasing plasma magnetization. We perform a detailed comparison of our results with those of recent PIC simulations and briefly discuss the astrophysical implications of our findings through the representative case of flaring events from blazar jets.

A high-order particle-in-cell method for low density plasma flow and the simulation of gyrotron resonator devices

International Nuclear Information System (INIS)

Stock, Andreas

2013-01-01

Within this thesis a parallelized, transient, three-dimensional, high-order discontinuous Galerkin Particle-in-Cell solver is developed and used to simulate the resonant cavity of a gyrotron. The high-order discontinuous Galerkin approach - a Finite-Element type method - provides a fast and efficient algorithm to numerically solve Maxwell's equations used within this thesis. Besides its outstanding dissipation and dispersion properties, the discontinuous Galerkin approach easily allows for using unstructured grids, as required to simulate complex-shaped engineering devices. The discontinuous Galerkin approach approximates a wavelength with significantly less degrees of freedom compared to other methods, e.g. Finite Difference methods. Furthermore, the parallelization capabilities of the discontinuous Galerkin framework are excellent due to the very local dependencies between the elements. These properties are essential for the efficient numerical treatment of the Vlasov-Maxwell system with the Particle-in-Cell method. This system describes the self-consistent interaction of charged particles and the electromagnetic field. As central application within this thesis gyrotron resonators are simulated with the discontinuous Galerkin Particle-in-Cell method on high-performance-computers. The gyrotron is a high-power millimeter wave source, used for the electron cyclotron resonance heating of magnetically confined fusion plasma, e.g. in the Wendelstein 7-X experimental fusion-reactor. Compared to state-of-the-art simulation tools used for the design of gyrotron resonators the Particle-in-Cell method does not use any significant physically simplifications w.r.t. the modelling of the particle-field-interaction, the geometry and the wave-spectrum. Hence, it is the method of choice for validation of current simulation tools being restricted by these simplifications. So far, the Particle-in-Cell method was restricted to be used for demonstration calculations only, because

A high-order particle-in-cell method for low density plasma flow and the simulation of gyrotron resonator devices

Energy Technology Data Exchange (ETDEWEB)

Stock, Andreas

2013-04-26

Within this thesis a parallelized, transient, three-dimensional, high-order discontinuous Galerkin Particle-in-Cell solver is developed and used to simulate the resonant cavity of a gyrotron. The high-order discontinuous Galerkin approach - a Finite-Element type method - provides a fast and efficient algorithm to numerically solve Maxwell's equations used within this thesis. Besides its outstanding dissipation and dispersion properties, the discontinuous Galerkin approach easily allows for using unstructured grids, as required to simulate complex-shaped engineering devices. The discontinuous Galerkin approach approximates a wavelength with significantly less degrees of freedom compared to other methods, e.g. Finite Difference methods. Furthermore, the parallelization capabilities of the discontinuous Galerkin framework are excellent due to the very local dependencies between the elements. These properties are essential for the efficient numerical treatment of the Vlasov-Maxwell system with the Particle-in-Cell method. This system describes the self-consistent interaction of charged particles and the electromagnetic field. As central application within this thesis gyrotron resonators are simulated with the discontinuous Galerkin Particle-in-Cell method on high-performance-computers. The gyrotron is a high-power millimeter wave source, used for the electron cyclotron resonance heating of magnetically confined fusion plasma, e.g. in the Wendelstein 7-X experimental fusion-reactor. Compared to state-of-the-art simulation tools used for the design of gyrotron resonators the Particle-in-Cell method does not use any significant physically simplifications w.r.t. the modelling of the particle-field-interaction, the geometry and the wave-spectrum. Hence, it is the method of choice for validation of current simulation tools being restricted by these simplifications. So far, the Particle-in-Cell method was restricted to be used for demonstration calculations only, because

Current-Sheet Formation and Reconnection at a Magnetic X Line in Particle-in-Cell Simulations

Science.gov (United States)

Black, C.; Antiochos, S. K.; Hesse, M.; Karpen, J. T.; Kuznetsova, M. M.; Zenitani, S.

2011-01-01

The integration of kinetic effects into macroscopic numerical models is currently of great interest to the heliophysics community, particularly in the context of magnetic reconnection. Reconnection governs the large-scale energy release and topological rearrangement of magnetic fields in a wide variety of laboratory, heliophysical, and astrophysical systems. We are examining the formation and reconnection of current sheets in a simple, two-dimensional X-line configuration using high-resolution particle-in-cell (PIC) simulations. The initial minimum-energy, potential magnetic field is perturbed by excess thermal pressure introduced into the particle distribution function far from the X line. Subsequently, the relaxation of this added stress leads self-consistently to the development of a current sheet that reconnects for imposed stress of sufficient strength. We compare the time-dependent evolution and final state of our PIC simulations with macroscopic magnetohydrodynamic simulations assuming both uniform and localized electrical resistivities (C. R. DeVore et al., this meeting), as well as with force-free magnetic-field equilibria in which the amount of reconnection across the X line can be constrained to be zero (ideal evolution) or optimal (minimum final magnetic energy). We will discuss implications of our results for understanding magnetic-reconnection onset and cessation at kinetic scales in dynamically formed current sheets, such as those occurring in the solar corona and terrestrial magnetotail.

Thermonuclear Tokamak plasmas in the presence of fusion alpha particles

International Nuclear Information System (INIS)

Anderson, D.; Hamnen, H.; Lisak, M.

1988-01-01

In this overview, we have focused on several results of the thermonuclear plasma research pertaining to the alpha particle physics and diagnostics in a fusion tokamak plasma. As regards the discussion of alpha particle effects, two distinct classes of phenomena have been distinguished: the simpler class containing phenomena exhibited by individual alpha particles under the influence of bulk plasma properties and, the more complex class including collective effects which become important for increasing alpha particle density. We have also discussed several possibilities to investigate alpha particle effects by simulation experiments using an equivalent population of highly energetic ions in the plasma. Generally, we find that the present theoretical knowledge on the role of fusion alpha particles in a fusion tokamak plasma is incomplete. There are still uncertainties and partial lack of quantitative results in this area. Consequently, further theoretical work and, as far a possible, simulation experiments are needed to improve the situation. Concerning the alpha particle diagnostics, the various diagnostic techniques and the status of their development have been discussed in two different contexts: the escaping alpha particles and the confined alpha particles in the fusion plasma. A general conclusion is that many of the different diagnostic methods for alpha particle measurements require further major development. (authors)

Thrust evaluation of magneto plasma sail that obtains an electromagnetic thrust from the solar wind

International Nuclear Information System (INIS)

Kajimura, Yoshihiro; Funaki, Ikkoh; Usui, Hideyuki; Yamakawa, Hiroshi

2011-01-01

Magneto Plasma Sail (MPS) is a propulsion system used in space, which generates its force by the interaction between the solar wind and an inflated magnetic field via a plasma injection. The quantitative evaluation of the thrust increment generated by injecting a plasma jet with a β in less than unity was conducted by three-dimensional hybrid particle-in-cell (PIC) simulations in an ion inertia scale. The injected plasma β in is 0.02 and the ratio of Larmor radius of injected ion to the representative length of the magnetic field is 0.5 at the injection point. In this situation, the obtained thrust of the MPS is 1.6 mN compared with the 0.2 mN of the thrust obtained by the pure magnetic sail since the induced current region on magnetosphere expanded by the magnetic inflation. (author)

Low-noise Collision Operators for Particle-in-cell Simulations

International Nuclear Information System (INIS)

Lewandowski, J.L.V.

2005-01-01

A new method to implement low-noise collision operators in particle-in-cell simulations is presented. The method is based on the fact that relevant collision operators can be included naturally in the Lagrangian formulation that exemplifies the particle-in-cell simulation method. Numerical simulations show that the momentum and energy conservation properties of the simulated plasma associated with the low-noise collision operator are improved as compared with standard collision algorithms based on random numbers

The TESS [Tandem Experiment Simulation Studies] computer code user's manual

International Nuclear Information System (INIS)

Procassini, R.J.

1990-01-01

TESS (Tandem Experiment Simulation Studies) is a one-dimensional, bounded particle-in-cell (PIC) simulation code designed to investigate the confinement and transport of plasma in a magnetic mirror device, including tandem mirror configurations. Mirror plasmas may be modeled in a system which includes an applied magnetic field and/or a self-consistent or applied electrostatic potential. The PIC code TESS is similar to the PIC code DIPSI (Direct Implicit Plasma Surface Interactions) which is designed to study plasma transport to and interaction with a solid surface. The codes TESS and DIPSI are direct descendants of the PIC code ES1 that was created by A. B. Langdon. This document provides the user with a brief description of the methods used in the code and a tutorial on the use of the code. 10 refs., 2 tabs

Overcoming Challenges in Kinetic Modeling of Magnetized Plasmas and Vacuum Electronic Devices

Science.gov (United States)

Omelchenko, Yuri; Na, Dong-Yeop; Teixeira, Fernando

2017-10-01

We transform the state-of-the art of plasma modeling by taking advantage of novel computational techniques for fast and robust integration of multiscale hybrid (full particle ions, fluid electrons, no displacement current) and full-PIC models. These models are implemented in 3D HYPERS and axisymmetric full-PIC CONPIC codes. HYPERS is a massively parallel, asynchronous code. The HYPERS solver does not step fields and particles synchronously in time but instead executes local variable updates (events) at their self-adaptive rates while preserving fundamental conservation laws. The charge-conserving CONPIC code has a matrix-free explicit finite-element (FE) solver based on a sparse-approximate inverse (SPAI) algorithm. This explicit solver approximates the inverse FE system matrix (``mass'' matrix) using successive sparsity pattern orders of the original matrix. It does not reduce the set of Maxwell's equations to a vector-wave (curl-curl) equation of second order but instead utilizes the standard coupled first-order Maxwell's system. We discuss the ability of our codes to accurately and efficiently account for multiscale physical phenomena in 3D magnetized space and laboratory plasmas and axisymmetric vacuum electronic devices.

Destabilizing effect of alpha particles in a Maxwellian plasma

International Nuclear Information System (INIS)

Wang, M.Y.

1976-01-01

Various plasma waves which are possibly excited by MeV alphas have been investigated. For a delta birth distribution it is found that: a) The right-circularly polarized Alfven wave can be excited. Its growth rate is linearly proportional to the α-particle density. b) The drift Alfven wave is stable against α-particles. c) For a uniform temperature, the plasma wave spectrum changes from three branches with n/sub α/ = 0 to four branches for n/sub α/ not equal to 0 case. d) α-particles can destabilize the ion drift acoustic wave even with uniform temperature. However, the ion acoustic wave appears to be stable against fusion products in a fusion grade plasma. e) If their effect on the background plasma spectrum is neglected, α-particles can excite the electromagnetic cyclotron wave in a range of harmonics (band structure). The growth rate is proportional to the square root of α-particle density. f) If the effect of α-particle on the plasma spectrum is included, we find that electromagnetic cyclotron wave is stable

Study of ultra-high gradient wakefield excitation by intense ultrashort laser pulses in plasma

International Nuclear Information System (INIS)

Kotaki, Hideyuki

2002-12-01

We investigate a mechanism of nonlinear phenomena in laser-plasma interaction, a laser wakefield excited by intense laser pulses, and the possibility of generating an intense bright electron source by an intense laser pulse. We need to understand and further employ some of these phenomena for our purposes. We measure self-focusing, filamentation, and the anomalous blueshift of the laser pulse. The ionization of gas with the self-focusing causes a broad continuous spectrum with blueshift. The normal blueshift depends on the laser intensity and the plasma density. We, however, have found different phenomenon. The laser spectrum shifts to fixed wavelength independent of the laser power and gas pressure above some critical power. We call the phenomenon 'anomalous blueshift'. The results are explained by the formation of filaments. An intense laser pulse can excite a laser wakefield in plasma. The coherent wakefield excited by 2 TW, 50 fs laser pulses in a gas-jet plasma around 10 18 cm -3 is measured with a time-resolved frequency domain interferometer (FDI). The density distribution of the helium gas is measured with a time-resolved Mach-Zehnder interferometer to search for the optimum laser focus position and timing in the gas-jet. The results show an accelerating wakefield excitation of 20 GeV/m with good coherency, which is useful for ultrahigh gradient particle acceleration in a compact system. This is the first time-resolved measurement of laser wakefield excitation in a gas-jet plasma. The experimental results are compared with a Particle-in-Cell (PIC) simulation. The pump-probe interferometer system of FDI and the anomalous blueshift will be modified to the optical injection system as a relativistic electron beam injector. In 1D PIC simulation we obtain the results of high quality intense electron beam acceleration. These results illuminate the possibility of a high energy and a high quality electron beam acceleration. (author)

Journal of Astrophysics and Astronomy | Indian Academy of Sciences

Indian Academy of Sciences (India)

In this paper, we propose a 3D kinetic model (particle-in-cell, PIC) for the description of the large scale Earth's bow shock. The proposed version is stable and does not require huge or extensive computer resources. Because PIC simulations work with scaled plasma and field parameters, we also propose to validate our ...

Two-dimensional particle-in-cell plasma source ion implantation of a prolate spheroid target

International Nuclear Information System (INIS)

Cheng-Sen, Liu; Hong-Ying, Han; Xiao-Qing, Peng; Ye, Chang; De-Zhen, Wang

2010-01-01

A two-dimensional particle-in-cell simulation is used to study the time-dependent evolution of the sheath surrounding a prolate spheroid target during a high voltage pulse in plasma source ion implantation. Our study shows that the potential contour lines pack more closely in the plasma sheath near the vertex of the major axis, i.e. where a thinner sheath is formed, and a non-uniform total ion dose distribution is incident along the surface of the prolate spheroid target due to the focusing of ions by the potential structure. Ion focusing takes place not only at the vertex of the major axis, where dense potential contour lines exist, but also at the vertex of the minor axis, where sparse contour lines exist. This results in two peaks of the received ion dose, locating at the vertices of the major and minor axes of the prolate spheroid target, and an ion dose valley, staying always between the vertices, rather than at the vertex of the minor axis

Two-dimensional particle-in-cell plasma source ion implantation of a prolate spheroid target

Science.gov (United States)

Liu, Cheng-Sen; Han, Hong-Ying; Peng, Xiao-Qing; Chang, Ye; Wang, De-Zhen

2010-03-01

A two-dimensional particle-in-cell simulation is used to study the time-dependent evolution of the sheath surrounding a prolate spheroid target during a high voltage pulse in plasma source ion implantation. Our study shows that the potential contour lines pack more closely in the plasma sheath near the vertex of the major axis, i.e. where a thinner sheath is formed, and a non-uniform total ion dose distribution is incident along the surface of the prolate spheroid target due to the focusing of ions by the potential structure. Ion focusing takes place not only at the vertex of the major axis, where dense potential contour lines exist, but also at the vertex of the minor axis, where sparse contour lines exist. This results in two peaks of the received ion dose, locating at the vertices of the major and minor axes of the prolate spheroid target, and an ion dose valley, staying always between the vertices, rather than at the vertex of the minor axis.

Particle-in-cell studies of fast-ion slowing-down rates in cool tenuous magnetized plasma

Science.gov (United States)

Evans, Eugene S.; Cohen, Samuel A.; Welch, Dale R.

2018-04-01

We report on 3D-3V particle-in-cell simulations of fast-ion energy-loss rates in a cold, weakly-magnetized, weakly-coupled plasma where the electron gyroradius, ρe, is comparable to or less than the Debye length, λDe, and the fast-ion velocity exceeds the electron thermal velocity, a regime in which the electron response may be impeded. These simulations use explicit algorithms, spatially resolve ρe and λDe, and temporally resolve the electron cyclotron and plasma frequencies. For mono-energetic dilute fast ions with isotropic velocity distributions, these scaling studies of the slowing-down time, τs, versus fast-ion charge are in agreement with unmagnetized slowing-down theory; with an applied magnetic field, no consistent anisotropy between τs in the cross-field and field-parallel directions could be resolved. Scaling the fast-ion charge is confirmed as a viable way to reduce the required computational time for each simulation. The implications of these slowing down processes are described for one magnetic-confinement fusion concept, the small, advanced-fuel, field-reversed configuration device.

The influence of laser-particle interaction in laser induced breakdown spectroscopy and laser ablation inductively coupled plasma spectrometry

International Nuclear Information System (INIS)

Lindner, Helmut; Loper, Kristofer H.; Hahn, David W.; Niemax, Kay

2011-01-01

Particles produced by previous laser shots may have significant influence on the analytical signal in laser-induced breakdown spectroscopy (LIBS) and laser ablation inductively coupled plasma (LA-ICP) spectrometry if they remain close to the position of laser sampling. The effects of these particles on the laser-induced breakdown event are demonstrated in several ways. LIBS-experiments were conducted in an ablation cell at atmospheric conditions in argon or air applying a dual-pulse arrangement with orthogonal pre-pulse, i.e., plasma breakdown in a gas generated by a focussed laser beam parallel and close to the sample surface followed by a delayed crossing laser pulse in orthogonal direction which actually ablates material from the sample and produces the LIBS plasma. The optical emission of the LIBS plasma as well as the absorption of the pre-pulse laser was measured. In the presence of particles in the focus of the pre-pulse laser, the plasma breakdown is affected and more energy of the pre-pulse laser is absorbed than without particles. As a result, the analyte line emission from the LIBS plasma of the second laser is enhanced. It is assumed that the enhancement is not only due to an increase of mass ablated by the second laser but also to better atomization and excitation conditions favored by a reduced gas density in the pre-pulse plasma. Higher laser pulse frequencies increase the probability of particle-laser interaction and, therefore, reduce the shot-to-shot line intensity variation as compared to lower particle loadings in the cell. Additional experiments using an aerosol chamber were performed to further quantify the laser absorption by the plasma in dependence on time both with and without the presence of particles. The overall implication of laser-particle interactions for LIBS and LA-ICP-MS/OES are discussed.

The influence of laser-particle interaction in laser induced breakdown spectroscopy and laser ablation inductively coupled plasma spectrometry

Science.gov (United States)

Lindner, Helmut; Loper, Kristofer H.; Hahn, David W.; Niemax, Kay

2011-02-01

Particles produced by previous laser shots may have significant influence on the analytical signal in laser-induced breakdown spectroscopy (LIBS) and laser ablation inductively coupled plasma (LA-ICP) spectrometry if they remain close to the position of laser sampling. The effects of these particles on the laser-induced breakdown event are demonstrated in several ways. LIBS-experiments were conducted in an ablation cell at atmospheric conditions in argon or air applying a dual-pulse arrangement with orthogonal pre-pulse, i.e., plasma breakdown in a gas generated by a focussed laser beam parallel and close to the sample surface followed by a delayed crossing laser pulse in orthogonal direction which actually ablates material from the sample and produces the LIBS plasma. The optical emission of the LIBS plasma as well as the absorption of the pre-pulse laser was measured. In the presence of particles in the focus of the pre-pulse laser, the plasma breakdown is affected and more energy of the pre-pulse laser is absorbed than without particles. As a result, the analyte line emission from the LIBS plasma of the second laser is enhanced. It is assumed that the enhancement is not only due to an increase of mass ablated by the second laser but also to better atomization and excitation conditions favored by a reduced gas density in the pre-pulse plasma. Higher laser pulse frequencies increase the probability of particle-laser interaction and, therefore, reduce the shot-to-shot line intensity variation as compared to lower particle loadings in the cell. Additional experiments using an aerosol chamber were performed to further quantify the laser absorption by the plasma in dependence on time both with and without the presence of particles. The overall implication of laser-particle interactions for LIBS and LA-ICP-MS/OES are discussed.

Wave-particle Interactions in Space and Laboratory Plasmas

Science.gov (United States)

An, Xin

This dissertation presents a study of wave-particle interactions in space and in the laboratory. To be concrete, the excitation of whistler-mode chorus waves in space and in the laboratory is studied in the first part. The relaxation of whistler anisotropy instability relevant to whistler-mode chorus waves in space is examined. Using a linear growth rate analysis and kinetic particle-in-cell simulations, the electron distributions are demonstrated to be well-constrained by the whistler anisotropy instability to a marginal-stability state, consistent with measurements by Van Allen Probes. The electron parallel beta beta ∥e separates the excited whistler waves into two groups: (i) quasi-parallel whistler waves for beta∥e > 0.02 and (ii) oblique whistler waves close to the resonance cone for beta∥e cell simulations. Motivated by the puzzles of chorus waves in space and by their recognized importance, the excitation of whistler-mode chorus waves is studied in the Large Plasma Device by the injection of a helical electron beam into a cold plasma. Incoherent broadband whistler waves similar to magnetospheric hiss are observed in the laboratory plasma. Their mode structures are identified by the phase-correlation technique. It is demonstrated that the waves are excited through a combination of Landau resonance, cyclotron resonance and anomalous cyclotron resonance. To account for the finite size effect of the electron beam, linear unstable eigenmodes of whistler waves are calculated by matching the eigenmode solution at the boundary. It is shown that the perpendicular wave number inside the beam is quantized due to the constraint imposed by the boundary condition. Darwin particle-in-cell simulations are carried out to study the simultaneous excitation of Langmuir and whistler waves in a beam-plasma system. The electron beam is first slowed down and relaxed by the rapidly growing Langmuir wave parallel to the background magnetic field. The tail of the core electrons

External circuit integration with electromagnetic particle in cell modeling of plasma focus devices

International Nuclear Information System (INIS)

Seng, Y. S.; Lee, P.; Rawat, R. S.

2015-01-01

The pinch performance of a plasma focus (PF) device is sensitive to the physical conditions of the breakdown phase. It is therefore essential to model and study the initial phase in order to optimize device performance. An external circuit is self consistently coupled to the electromagnetic particle in cell code to model the breakdown and initial lift phase of the United Nations University/International Centre for Theoretical Physics (UNU-ICTP) plasma focus device. Gas breakdown during the breakdown phase is simulated successfully, following a drop in the applied voltage across the device and a concurrent substantial rise in the circuit current. As a result, the plasma becomes magnetized, with the growing value of the magnetic field over time leading to the gradual lift off of the well formed current sheath into the axial acceleration phase. This lifting off, with simultaneous outward sheath motion along the anode and vertical cathode, and the strong magnetic fields in the current sheath region, was demonstrated in this work, and hence validates our method of coupling the external circuit to PF devices. Our method produces voltage waveforms that are qualitatively similar to the observed experimental voltage profiles of the UNU-ICTP device. Values of the mean electron energy before and after voltage breakdown turned out to be different, with the values after breakdown being much lower. In both cases, the electron energy density function turned out to be non-Maxwellian

A parallel 3D particle-in-cell code with dynamic load balancing

International Nuclear Information System (INIS)

Wolfheimer, Felix; Gjonaj, Erion; Weiland, Thomas

2006-01-01

A parallel 3D electrostatic Particle-In-Cell (PIC) code including an algorithm for modelling Space Charge Limited (SCL) emission [E. Gjonaj, T. Weiland, 3D-modeling of space-charge-limited electron emission. A charge conserving algorithm, Proceedings of the 11th Biennial IEEE Conference on Electromagnetic Field Computation, 2004] is presented. A domain decomposition technique based on orthogonal recursive bisection is used to parallelize the computation on a distributed memory environment of clustered workstations. For problems with a highly nonuniform and time dependent distribution of particles, e.g., bunch dynamics, a dynamic load balancing between the processes is needed to preserve the parallel performance. The algorithm for the detection of a load imbalance and the redistribution of the tasks among the processes is based on a weight function criterion, where the weight of a cell measures the computational load associated with it. The algorithm is studied with two examples. In the first example, multiple electron bunches as occurring in the S-DALINAC [A. Richter, Operational experience at the S-DALINAC, Proceedings of the Fifth European Particle Accelerator Conference, 1996] accelerator are simulated in the absence of space charge fields. In the second example, the SCL emission and electron trajectories in an electron gun are simulated

A parallel 3D particle-in-cell code with dynamic load balancing

Energy Technology Data Exchange (ETDEWEB)

Wolfheimer, Felix [Technische Universitaet Darmstadt, Institut fuer Theorie Elektromagnetischer Felder, Schlossgartenstr.8, 64283 Darmstadt (Germany)]. E-mail: wolfheimer@temf.de; Gjonaj, Erion [Technische Universitaet Darmstadt, Institut fuer Theorie Elektromagnetischer Felder, Schlossgartenstr.8, 64283 Darmstadt (Germany); Weiland, Thomas [Technische Universitaet Darmstadt, Institut fuer Theorie Elektromagnetischer Felder, Schlossgartenstr.8, 64283 Darmstadt (Germany)

2006-03-01

A parallel 3D electrostatic Particle-In-Cell (PIC) code including an algorithm for modelling Space Charge Limited (SCL) emission [E. Gjonaj, T. Weiland, 3D-modeling of space-charge-limited electron emission. A charge conserving algorithm, Proceedings of the 11th Biennial IEEE Conference on Electromagnetic Field Computation, 2004] is presented. A domain decomposition technique based on orthogonal recursive bisection is used to parallelize the computation on a distributed memory environment of clustered workstations. For problems with a highly nonuniform and time dependent distribution of particles, e.g., bunch dynamics, a dynamic load balancing between the processes is needed to preserve the parallel performance. The algorithm for the detection of a load imbalance and the redistribution of the tasks among the processes is based on a weight function criterion, where the weight of a cell measures the computational load associated with it. The algorithm is studied with two examples. In the first example, multiple electron bunches as occurring in the S-DALINAC [A. Richter, Operational experience at the S-DALINAC, Proceedings of the Fifth European Particle Accelerator Conference, 1996] accelerator are simulated in the absence of space charge fields. In the second example, the SCL emission and electron trajectories in an electron gun are simulated.

Modelling of plasma generation and expansion in a vacuum arc: application to the vacuum arc remelting process

International Nuclear Information System (INIS)

Chapelle, P.; Bellot, J.P.; Duval, H.; Jardy, A.; Ablitzer, D.

2002-01-01

As part of a complete theoretical description of the behaviour of the electric arc in the vacuum arc remelting process, a model has been developed for the column of plasma generated by a single cluster of cathode spots. The model combines a kinetic approach, taking into account the formation of the plasma in the cathodic region, and a hydrodynamic approach, describing the expansion of the plasma in the vacuum between the electrodes. The kinetic model is based on a system of Boltzmann-Vlasov-Poisson equations and uses a particle-type simulation procedure, combining the PIC (particle in cell) and FPM (finite point set method) methods. In the two-dimensional hydrodynamic model, the plasma is assimilated to a mixture of two continuous fluids (the electrons and the ions), each described by a system of coupled transport equations. Finally, a simplified method has been defined for calculating the electric current density and the energy flux density transmitted by the plasma to the anode. The results of the numerical simulation presented are consistent with a certain number of experimental data available in the literature. In particular, the model predicts a percentage of the electric power of the cluster transmitted to the anode (25%) in good agreement with the value indicated in the literature. (author)

Particle modeling of plasmas computational plasma physics

International Nuclear Information System (INIS)

Dawson, J.M.

1991-01-01

Recently, through the development of supercomputers, a powerful new method for exploring plasmas has emerged; it is computer modeling of plasmas. Such modeling can duplicate many of the complex processes that go on in a plasma and allow scientists to understand what the important processes are. It helps scientists gain an intuition about this complex state of matter. It allows scientists and engineers to explore new ideas on how to use plasma before building costly experiments; it allows them to determine if they are on the right track. It can duplicate the operation of devices and thus reduce the need to build complex and expensive devices for research and development. This is an exciting new endeavor that is in its infancy, but which can play an important role in the scientific and technological competitiveness of the US. There are a wide range of plasma models that are in use. There are particle models, fluid models, hybrid particle fluid models. These can come in many forms, such as explicit models, implicit models, reduced dimensional models, electrostatic models, magnetostatic models, electromagnetic models, and almost an endless variety of other models. Here the author will only discuss particle models. He will give a few examples of the use of such models; these will be taken from work done by the Plasma Modeling Group at UCLA because he is most familiar with work. However, it only gives a small view of the wide range of work being done around the US, or for that matter around the world

Electron vortex magnetic holes: A nonlinear coherent plasma structure

Energy Technology Data Exchange (ETDEWEB)

Haynes, Christopher T., E-mail: c.t.haynes@qmul.ac.uk; Burgess, David; Sundberg, Torbjorn [School of Physics and Astronomy, Queen Mary University of London, Mile End Road, London E1 4NS (United Kingdom); Camporeale, Enrico [Multiscale Dynamics, Centrum Wiskunde and Informatica (CWI), Amsterdam (Netherlands)

2015-01-15

We report the properties of a novel type of sub-proton scale magnetic hole found in two dimensional particle-in-cell simulations of decaying turbulence with a guide field. The simulations were performed with a realistic value for ion to electron mass ratio. These structures, electron vortex magnetic holes (EVMHs), have circular cross-section. The magnetic field depression is associated with a diamagnetic azimuthal current provided by a population of trapped electrons in petal-like orbits. The trapped electron population provides a mean azimuthal velocity and since trapping preferentially selects high pitch angles, a perpendicular temperature anisotropy. The structures arise out of initial perturbations in the course of the turbulent evolution of the plasma, and are stable over at least 100 electron gyroperiods. We have verified the model for the EVMH by carrying out test particle and PIC simulations of isolated structures in a uniform plasma. It is found that (quasi-)stable structures can be formed provided that there is some initial perpendicular temperature anisotropy at the structure location. The properties of these structures (scale size, trapped population, etc.) are able to explain the observed properties of magnetic holes in the terrestrial plasma sheet. EVMHs may also contribute to turbulence properties, such as intermittency, at short scale lengths in other astrophysical plasmas.

Propagation of an intense laser pulse in an under-dense plasma: channeling and stimulated Raman scattering

International Nuclear Information System (INIS)

Friou, A.

2012-01-01

This thesis is divided in two parts: i) the laser channeling in hundreds of microns long under-dense plasmas (0.1 nc ≤ n ≤ nc, nc being the critical density) of a laser pulse of intensity 10 18-20 W/cm 2 and duration 1-10 ps; ii) the saturation mechanisms of stimulated Raman back-scattering of a laser pulse of intensity 10 14 to 10 16 W/cm 2 and duration of about 1 ps. A parametric study was performed to study the channeling of a very intense laser pulse, using a 2D PIC (Particle In Cell) code. Various kinds of channels were obtained depending on the laser and plasma parameters, thereby reproducing and enlarging previous studies. Moreover, the channeling velocity was measured and scaling laws were established for homogeneous plasmas. They are then applied to inhomogeneous plasmas, similar to those encountered in inertial confinement fusion (ICF). It is then possible to estimate the energy necessary to channel to the critical density, an important step for the fast ignition scheme of ICF. Raman saturation was studied using numerical simulations, in order to determine if it is due to dephasing or to the growth of sidebands, using different approaches. The first is to study Raman simulations (electromagnetic) performed with kinetic PIC and Vlasov codes. The second, is to study the evolution of a plasma initialized with a distribution function after the adiabatic theory, using a Vlasov code (electrostatic). In this case, we observe the growth of a sideband, with dominant wave number and growth rate in good agreement with kinetic simulations. The saturation of the plasma wave can be caused by both saturation mechanisms. [fr

Computer simulation of phenomena in plasma via particles

International Nuclear Information System (INIS)

Alves, M.V.; Bittencourt, J.A.

1988-06-01

The method of plasma computer simulation via particles has become an efficient tool to investigate the time and spatial evolution of various physical phenomena in plasmas. This method is based on the study of the individual plasma particle motions interacting with one another and with the externally applied fields. Although fairly simple, it allows a non-linear analysis of complex plasma physical phenomena and to obtain diagnostics even for regions of the system where experimental measurements would be difficult to make. In this report, a general view of the electrostatic one-dimensional computer code ES1, originally developed by A. Bruce Langdon, is presented. The main mathematical artifice in this code is the use of a spatial grid in which various plasma particles are represented by ''superparticles'', using a given shape function. The principal characteristics of the model, the approximations made and the mathematical methods used to solve the equations involved, are described. The specification of the input parameters which characterize the system, the initial conditions and the graphic diagnostics which can be utilized, are also described. Results are presented illustrating graphically the behavior of the plasma oscillations, the two-stream instability and the beam-plasma instability. (author) [pt

Plasma simulation and fusion calculation

International Nuclear Information System (INIS)