An electrostatic Particle-In-Cell code on multi-block structured meshes
[1];
[2];
[2];
[3];
[2]
- Los Alamos National Lab. (LANL), Los Alamos, NM (United States). Intelligence and Space Research Division
- Los Alamos National Lab. (LANL), Los Alamos, NM (United States). Theoretical Division
- Los Alamos National Lab. (LANL), Los Alamos, NM (United States). Computer, Computational, and Statistical Sciences Division
We present an electrostatic Particle-In-Cell (PIC) code on multi-block, locally structured, curvilinear meshes called Curvilinear PIC (CPIC). Multi-block meshes are essential to capture complex geometries accurately and with good mesh quality, something that would not be possible with single-block structured meshes that are often used in PIC and for which CPIC was initially developed. In spite of the structured nature of the individual blocks, multi-block meshes resemble unstructured meshes in a global sense and introduce several new challenges, such as the presence of discontinuities in the mesh properties and coordinate orientation changes across adjacent blocks, and polyjunction points where an arbitrary number of blocks meet. In CPIC, these challenges have been met by an approach that features: (1) a curvilinear formulation of the PIC method: each mesh block is mapped from the physical space, where the mesh is curvilinear and arbitrarily distorted, to the logical space, where the mesh is uniform and Cartesian on the unit cube; (2) a mimetic discretization of Poisson's equation suitable for multi-block meshes; and (3) a hybrid (logical-space position/physical-space velocity), asynchronous particle mover that mitigates the performance degradation created by the necessity to track particles as they move across blocks. The numerical accuracy of CPIC was verified using two standard plasma–material interaction tests, which demonstrate good agreement with the corresponding analytic solutions. And compared to PIC codes on unstructured meshes, which have also been used for their flexibility in handling complex geometries but whose performance suffers from issues associated with data locality and indirect data access patterns, PIC codes on multi-block structured meshes may offer the best compromise for capturing complex geometries while also maintaining solution accuracy and computational efficiency.
- Research Organization:
- Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
- Sponsoring Organization:
- USDOE National Nuclear Security Administration (NNSA)
- Grant/Contract Number:
- AC52-06NA25396
- OSTI ID:
- 1392872
- Alternate ID(s):
- OSTI ID: 1495566
- Report Number(s):
- LA-UR-16-28543
- Journal Information:
- Journal of Computational Physics, Vol. 350, Issue C; ISSN 0021-9991
- Publisher:
- ElsevierCopyright Statement
- Country of Publication:
- United States
- Language:
- English
Web of Science
Molecular dynamics calculation of the spectral densities of plasma fluctuations
|
journal | June 2018 |
Similar Records
EMAPS: An Efficient Multiscale Approach to Plasma Systems with Non-MHD Scale Effects
Electromagnetic PIC simulations using finite elements on unstructured grids. [Particle-in-cell (PIC)]