# 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, three-dimensional 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 three-dimensional 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 particles-in-cell 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; THREE-DIMENSIONAL CALCULATIONS; TRAPPING

### Citation Formats

```
Davidson, A., E-mail: davidsoa@physics.ucla.edu, Tableman, A., E-mail: Tableman@physics.ucla.edu, An, W., E-mail: anweiming@ucla.edu, Tsung, F.S., E-mail: tsung@physics.ucla.edu, Lu, W., E-mail: luwei@ucla.edu, Department of Engineering Physics, Tsinghua University, Beijing 100084, Vieira, J., E-mail: jorge.vieira@ist.utl.pt, Fonseca, R.A., E-mail: ricardo.fonseca@iscte.pt, Departamento Ciências e Tecnologias da Informação, ISCTE – Instituto Universitário de Lisboa, 1649-026 Lisboa, Silva, L.O., E-mail: luis.silva@ist.utl.pt, Mori, W.B., E-mail: 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., E-mail: davidsoa@physics.ucla.edu, Tableman, A., E-mail: Tableman@physics.ucla.edu, An, W., E-mail: anweiming@ucla.edu, Tsung, F.S., E-mail: tsung@physics.ucla.edu, Lu, W., E-mail: luwei@ucla.edu, Department of Engineering Physics, Tsinghua University, Beijing 100084, Vieira, J., E-mail: jorge.vieira@ist.utl.pt, Fonseca, R.A., E-mail: ricardo.fonseca@iscte.pt, Departamento Ciências e Tecnologias da Informação, ISCTE – Instituto Universitário de Lisboa, 1649-026 Lisboa, Silva, L.O., E-mail: luis.silva@ist.utl.pt, Mori, W.B., E-mail: 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., E-mail: davidsoa@physics.ucla.edu, Tableman, A., E-mail: Tableman@physics.ucla.edu, An, W., E-mail: anweiming@ucla.edu, Tsung, F.S., E-mail: tsung@physics.ucla.edu, Lu, W., E-mail: luwei@ucla.edu, Department of Engineering Physics, Tsinghua University, Beijing 100084, Vieira, J., E-mail: jorge.vieira@ist.utl.pt, Fonseca, R.A., E-mail: ricardo.fonseca@iscte.pt, Departamento Ciências e Tecnologias da Informação, ISCTE – Instituto Universitário de Lisboa, 1649-026 Lisboa, Silva, L.O., E-mail: luis.silva@ist.utl.pt, Mori, W.B., E-mail: mori@physics.ucla.edu, and Department of Electrical Engineering, University of California, Los Angeles, CA 90095. Thu .
"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., E-mail: davidsoa@physics.ucla.edu and Tableman, A., E-mail: Tableman@physics.ucla.edu and An, W., E-mail: anweiming@ucla.edu and Tsung, F.S., E-mail: tsung@physics.ucla.edu and Lu, W., E-mail: luwei@ucla.edu and Department of Engineering Physics, Tsinghua University, Beijing 100084 and Vieira, J., E-mail: jorge.vieira@ist.utl.pt and Fonseca, R.A., E-mail: ricardo.fonseca@iscte.pt and Departamento Ciências e Tecnologias da Informação, ISCTE – Instituto Universitário de Lisboa, 1649-026 Lisboa and Silva, L.O., E-mail: luis.silva@ist.utl.pt and Mori, W.B., E-mail: mori@physics.ucla.edu and Department of Electrical Engineering, University of California, Los Angeles, CA 90095},

abstractNote = {For many plasma physics problems, three-dimensional 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 three-dimensional 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 particles-in-cell 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 = {Thu Jan 15 00:00:00 EST 2015},

month = {Thu Jan 15 00:00:00 EST 2015}

}