Implementation of a hybrid particle code with a PIC description in r-z and a gridless description in φ into OSIRIS

JOURNAL OF COMPUTATIONAL PHYSICS

卷:281

Article

Davidson, A.¹  Tableman, A.¹  An, W.¹  Tsung, F. S.¹  Lu, W.^4,5  Vieira, J.²  Fonseca, R. A.^2,3  Silva, L. O.²  Mori, W. B.^1,4 

[1] Univ Calif Los Angeles, Dept Phys & Astron, Los Angeles, CA 90095 USA

[2] Univ Nova Lisboa, Inst Super Tecn, GoLP Inst Plasmas & Fusao Nucl, P-1200 Lisbon, Portugal

[3] ISCTE Inst Univ Lisboa, Dept Ciencias & Tecnol Informacao, P-1649026 Lisbon, Portugal

[4] Univ Calif Los Angeles, Dept Elect Engn, Los Angeles, CA 90095 USA

[5] Tsinghua Univ, Dept Engn Phys, Beijing 100084, Peoples R China

关键词: LWFA;    PWFA;    PIC;    OSIRIS;    Self-trapping;    Hosing;   

DOI  :  10.1016/j.jcp.2014.10.064

来源: Elsevier

【 摘 要 】

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 phi. 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. (C) 2014 Elsevier Inc. All rights reserved.

【 授权许可】

Free