Implementation of a hybrid particle code with a PIC description in r–z and a gridless description in ϕ into OSIRIS
Abstract
For many plasma physics problems, threedimensional and kinetic effects are very important. However, such simulations are very computationally intensive. Fortunately, there is a class of problems for which there is nearly azimuthal symmetry and the dominant threedimensional physics is captured by the inclusion of only a few azimuthal harmonics. Recently, it was proposed [1] to model one such problem, laser wakefield acceleration, by expanding the fields and currents in azimuthal harmonics and truncating the expansion. The complex amplitudes of the fundamental and first harmonic for the fields were solved on an r–z grid and a procedure for calculating the complex current amplitudes for each particle based on its motion in Cartesian geometry was presented using a Marder's correction to maintain the validity of Gauss's law. In this paper, we describe an implementation of this algorithm into OSIRIS using a rigorous charge conserving current deposition method to maintain the validity of Gauss's law. We show that this algorithm is a hybrid method which uses a particlesincell description in r–z and a gridless description in ϕ. We include the ability to keep an arbitrary number of harmonics and higher order particle shapes. Examples for laser wakefield acceleration, plasma wakefield acceleration, and beammore »
 Authors:
 Department of Physics and Astronomy, University of California, Los Angeles, CA 90095 (United States)
 Department of Electrical Engineering, University of California, Los Angeles, CA 90095 (United States)
 (China)
 GoLP/Instituto de Plasmas e Fusão Nuclear, Instituto Superior Técnico, Universidade de Lisboa, Lisbon (Portugal)
 (Portugal)
 (United States)
 Publication Date:
 OSTI Identifier:
 22382177
 Resource Type:
 Journal Article
 Resource Relation:
 Journal Name: Journal of Computational Physics; Journal Volume: 281; Other Information: Copyright (c) 2014 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved.; Country of input: International Atomic Energy Agency (IAEA)
 Country of Publication:
 United States
 Language:
 English
 Subject:
 71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; 70 PLASMA PHYSICS AND FUSION TECHNOLOGY; ACCELERATION; ALGORITHMS; AMPLITUDES; COMPUTERIZED SIMULATION; CORRECTIONS; HARMONICS; IMPLEMENTATION; LOADING; PLASMA; PLASMA SIMULATION; SYMMETRY; THREEDIMENSIONAL CALCULATIONS; TRAPPING
Citation Formats
Davidson, A., Email: davidsoa@physics.ucla.edu, Tableman, A., Email: Tableman@physics.ucla.edu, An, W., Email: anweiming@ucla.edu, Tsung, F.S., Email: tsung@physics.ucla.edu, Lu, W., Email: luwei@ucla.edu, Department of Engineering Physics, Tsinghua University, Beijing 100084, Vieira, J., Email: jorge.vieira@ist.utl.pt, Fonseca, R.A., Email: ricardo.fonseca@iscte.pt, Departamento Ciências e Tecnologias da Informação, ISCTE – Instituto Universitário de Lisboa, 1649026 Lisboa, Silva, L.O., Email: luis.silva@ist.utl.pt, Mori, W.B., Email: mori@physics.ucla.edu, and Department of Electrical Engineering, University of California, Los Angeles, CA 90095. Implementation of a hybrid particle code with a PIC description in r–z and a gridless description in ϕ into OSIRIS. United States: N. p., 2015.
Web. doi:10.1016/J.JCP.2014.10.064.
Davidson, A., Email: davidsoa@physics.ucla.edu, Tableman, A., Email: Tableman@physics.ucla.edu, An, W., Email: anweiming@ucla.edu, Tsung, F.S., Email: tsung@physics.ucla.edu, Lu, W., Email: luwei@ucla.edu, Department of Engineering Physics, Tsinghua University, Beijing 100084, Vieira, J., Email: jorge.vieira@ist.utl.pt, Fonseca, R.A., Email: ricardo.fonseca@iscte.pt, Departamento Ciências e Tecnologias da Informação, ISCTE – Instituto Universitário de Lisboa, 1649026 Lisboa, Silva, L.O., Email: luis.silva@ist.utl.pt, Mori, W.B., Email: mori@physics.ucla.edu, & Department of Electrical Engineering, University of California, Los Angeles, CA 90095. Implementation of a hybrid particle code with a PIC description in r–z and a gridless description in ϕ into OSIRIS. United States. doi:10.1016/J.JCP.2014.10.064.
Davidson, A., Email: davidsoa@physics.ucla.edu, Tableman, A., Email: Tableman@physics.ucla.edu, An, W., Email: anweiming@ucla.edu, Tsung, F.S., Email: tsung@physics.ucla.edu, Lu, W., Email: luwei@ucla.edu, Department of Engineering Physics, Tsinghua University, Beijing 100084, Vieira, J., Email: jorge.vieira@ist.utl.pt, Fonseca, R.A., Email: ricardo.fonseca@iscte.pt, Departamento Ciências e Tecnologias da Informação, ISCTE – Instituto Universitário de Lisboa, 1649026 Lisboa, Silva, L.O., Email: luis.silva@ist.utl.pt, Mori, W.B., Email: mori@physics.ucla.edu, and Department of Electrical Engineering, University of California, Los Angeles, CA 90095. 2015.
"Implementation of a hybrid particle code with a PIC description in r–z and a gridless description in ϕ into OSIRIS". United States.
doi:10.1016/J.JCP.2014.10.064.
@article{osti_22382177,
title = {Implementation of a hybrid particle code with a PIC description in r–z and a gridless description in ϕ into OSIRIS},
author = {Davidson, A., Email: davidsoa@physics.ucla.edu and Tableman, A., Email: Tableman@physics.ucla.edu and An, W., Email: anweiming@ucla.edu and Tsung, F.S., Email: tsung@physics.ucla.edu and Lu, W., Email: luwei@ucla.edu and Department of Engineering Physics, Tsinghua University, Beijing 100084 and Vieira, J., Email: jorge.vieira@ist.utl.pt and Fonseca, R.A., Email: ricardo.fonseca@iscte.pt and Departamento Ciências e Tecnologias da Informação, ISCTE – Instituto Universitário de Lisboa, 1649026 Lisboa and Silva, L.O., Email: luis.silva@ist.utl.pt and Mori, W.B., Email: mori@physics.ucla.edu and Department of Electrical Engineering, University of California, Los Angeles, CA 90095},
abstractNote = {For many plasma physics problems, threedimensional and kinetic effects are very important. However, such simulations are very computationally intensive. Fortunately, there is a class of problems for which there is nearly azimuthal symmetry and the dominant threedimensional physics is captured by the inclusion of only a few azimuthal harmonics. Recently, it was proposed [1] to model one such problem, laser wakefield acceleration, by expanding the fields and currents in azimuthal harmonics and truncating the expansion. The complex amplitudes of the fundamental and first harmonic for the fields were solved on an r–z grid and a procedure for calculating the complex current amplitudes for each particle based on its motion in Cartesian geometry was presented using a Marder's correction to maintain the validity of Gauss's law. In this paper, we describe an implementation of this algorithm into OSIRIS using a rigorous charge conserving current deposition method to maintain the validity of Gauss's law. We show that this algorithm is a hybrid method which uses a particlesincell description in r–z and a gridless description in ϕ. We include the ability to keep an arbitrary number of harmonics and higher order particle shapes. Examples for laser wakefield acceleration, plasma wakefield acceleration, and beam loading are also presented and directions for future work are discussed.},
doi = {10.1016/J.JCP.2014.10.064},
journal = {Journal of Computational Physics},
number = ,
volume = 281,
place = {United States},
year = 2015,
month = 1
}

Gridless retarding potential analyzer for use in verylowenergy charged particle detection
The theory of the hyperbolic retarding potential analyzer in the electrostatic mode is developed in detail and verified in the laboratory. A monoenergetic electron beam is used for the laboratory investigation. The analyzer (acronym HARP) has advantages over other conventional electrostatic analyzers; among them are less contact potential influence and high throughput because of the symmetry shape of the analyzer. The most useful application of the HARP is in detecting lowenergy charged particles. A sample of lowenergy particle data obtained in the earth's ionosphere is given. (AIP) 
An efficient mixedprecision, hybrid CPUGPU implementation of a nonlinearly implicit onedimensional particleincell algorithm
Recently, a fully implicit, energy and chargeconserving particleincell method has been developed for multiscale, fullf kinetic simulations [G. Chen, et al., J. Comput. Phys. 230, 18 (2011)]. The method employs a Jacobianfree NewtonKrylov (JFNK) solver and is capable of using very large timesteps without loss of numerical stability or accuracy. A fundamental feature of the method is the segregation of particle orbit integrations from the field solver, while remaining fully selfconsistent. This provides great flexibility, and dramatically improves the solver efficiency by reducing the degrees of freedom of the associated nonlinear system. However, it requires a particle push per nonlinearmore »