Higher-order particle representation for particle-in-cell simulations

JOURNAL OF COMPUTATIONAL PHYSICS

卷:435

Article

Brown, Dominic A. S.¹  Bettencourt, Matthew T.²  Wright, Steven A.³  Maheswaran, Satheesh⁴  Jones, John P.⁵  Jarvis, Stephen A.⁶ 

[1] Univ Warwick, Dept Comp Sci, Coventry, W Midlands, England

[2] Sandia Natl Labs, POB 5800, Albuquerque, NM 87185 USA

[3] Univ York, Dept Comp Sci, York, N Yorkshire, England

[4] Diamond Light Source Ltd, Diamond House,Harwell Sci & Innovat Campus, Didcot, Oxon, England

[5] Atom Weap Estab, Aldermaston, England

[6] Univ Birmingham, Coll Engn & Phys Sci, Birmingham, W Midlands, England

关键词: Particle-in-cell;    High-order;    Unstructured;    Particle representation;    Shape function;   

DOI  :  10.1016/j.jcp.2021.110255

来源: Elsevier

【 摘 要 】

In this paper we present an alternative approach to the representation of simulation particles for unstructured electrostatic and electromagnetic PIC simulations. In our modified PIC algorithm we represent particles as having a smooth shape function limited by some specified finite radius, r(0). A unique feature of our approach is the representation of this shape by surrounding simulation particles with a set of virtual particles with delta shape, with fixed offsets and weights derived from Gaussian quadrature rules and the value of r(0). As the virtual particles are purely computational, they provide the additional benefit of increasing the arithmetic intensity of traditionally memory bound particle kernels. The modified algorithm is implemented within Sandia National Laboratories' unstructured EMPIRE-PIC code, for electrostatic and electromagnetic simulations, using periodic boundary conditions. We show results for a representative set of benchmark problems, including electron orbit, a transverse electromagnetic wave propagating through a plasma, numerical heating, and a plasma slab expansion. Good error reduction across all of the chosen problems is achieved as the particles are made progressively smoother, with the optimal particle radius appearing to be problem-dependent. Crown Copyright (c) 2021 Published by Elsevier Inc. All rights reserved.

【 授权许可】

Free