A new method for analyzing and visualizing plasma simulations using a phase-space tessellation
Abstract
We apply a novel phase-space interpolation technique referred to as the simplex-in-cell (SIC) method to analyze two- and three-dimensional particle-in-cell (PIC) simulations of electromagnetic plasmas. SIC relies on a discretization of the initial phase-space distribution function into simplices, which allows an approximation to the full, continuously defined distribution function to be constructed at any later time in the simulation. This allows densities, currents, and even full momentum distribution functions to be measured at any point in the simulation domain without averaging over control volumes. The SIC approach applies to any PIC simulation for which a tessellation of the initial particle distribution can be constructed. In this study, we use outputs from standard PIC simulations of the Weibel instability and compare physical quantities such as charge and current densities calculated in postprocessing using SIC and standard particle deposits. Using 2D simulations with 1–65 536 particles-per-cell, we find that SIC eliminates discrete particle noise and in some cases can reach a given noise level using ~1000 times fewer simulation particles than with standard particle deposition schemes. In regions of low density, such as between current filaments, SIC is able to capture small amplitude features even with fewer particles than gridpoints due tomore »
- Authors:
-
- Stanford Univ., CA (United States). Dept. of Physics, and Kavli Inst. for Particle Astrophysics and Cosmology; SLAC National Accelerator Lab., Menlo Park, CA (United States). High Energy Density Science Division
- SLAC National Accelerator Lab., Menlo Park, CA (United States). High Energy Density Science Division
- Stanford Univ., CA (United States). Dept. of Physics, and Kavli Inst. for Particle Astrophysics and Cosmology; SLAC National Accelerator Lab., Menlo Park, CA (United States)
- Publication Date:
- Research Org.:
- SLAC National Accelerator Lab., Menlo Park, CA (United States)
- Sponsoring Org.:
- USDOE Office of Science (SC), Fusion Energy Sciences (FES); National Science Foundation (NSF)
- OSTI Identifier:
- 1468691
- Alternate Identifier(s):
- OSTI ID: 1459754
- Grant/Contract Number:
- AC11339893; AC02-76SF00515; 100237
- Resource Type:
- Accepted Manuscript
- Journal Name:
- Physics of Plasmas
- Additional Journal Information:
- Journal Volume: 25; Journal Issue: 7; Journal ID: ISSN 1070-664X
- Publisher:
- American Institute of Physics (AIP)
- Country of Publication:
- United States
- Language:
- English
- Subject:
- 70 PLASMA PHYSICS AND FUSION TECHNOLOGY
Citation Formats
Totorica, Samuel R., Fiuza, Frederico, and Abel, Tom. A new method for analyzing and visualizing plasma simulations using a phase-space tessellation. United States: N. p., 2018.
Web. doi:10.1063/1.5037348.
Totorica, Samuel R., Fiuza, Frederico, & Abel, Tom. A new method for analyzing and visualizing plasma simulations using a phase-space tessellation. United States. doi:https://doi.org/10.1063/1.5037348
Totorica, Samuel R., Fiuza, Frederico, and Abel, Tom. Wed .
"A new method for analyzing and visualizing plasma simulations using a phase-space tessellation". United States. doi:https://doi.org/10.1063/1.5037348. https://www.osti.gov/servlets/purl/1468691.
@article{osti_1468691,
title = {A new method for analyzing and visualizing plasma simulations using a phase-space tessellation},
author = {Totorica, Samuel R. and Fiuza, Frederico and Abel, Tom},
abstractNote = {We apply a novel phase-space interpolation technique referred to as the simplex-in-cell (SIC) method to analyze two- and three-dimensional particle-in-cell (PIC) simulations of electromagnetic plasmas. SIC relies on a discretization of the initial phase-space distribution function into simplices, which allows an approximation to the full, continuously defined distribution function to be constructed at any later time in the simulation. This allows densities, currents, and even full momentum distribution functions to be measured at any point in the simulation domain without averaging over control volumes. The SIC approach applies to any PIC simulation for which a tessellation of the initial particle distribution can be constructed. In this study, we use outputs from standard PIC simulations of the Weibel instability and compare physical quantities such as charge and current densities calculated in postprocessing using SIC and standard particle deposits. Using 2D simulations with 1–65 536 particles-per-cell, we find that SIC eliminates discrete particle noise and in some cases can reach a given noise level using ~1000 times fewer simulation particles than with standard particle deposition schemes. In regions of low density, such as between current filaments, SIC is able to capture small amplitude features even with fewer particles than gridpoints due to the deformable nature of the SIC volume elements. Here, by calculating momentum distributions, we show how SIC can capture low density tails in the spectrum using far fewer particles than are necessary for standard particle deposits. We calculate the charge density on spatial grids of increasing resolution to demonstrate the ability of SIC to reveal fine-scale details that are not accessible with standard particle deposits. Finally, we show how SIC can be extended to 3D and give an example of its use to calculate the charge density from 3D PIC simulations of the Weibel instability. These results motivate the future implementation of SIC directly in the simulation force calculation for a novel low-noise electromagnetic plasma simulation method.},
doi = {10.1063/1.5037348},
journal = {Physics of Plasmas},
number = 7,
volume = 25,
place = {United States},
year = {2018},
month = {7}
}
Works referenced in this record:
Rigorous charge conservation for local electromagnetic field solvers
journal, March 1992
- Villasenor, John; Buneman, Oscar
- Computer Physics Communications, Vol. 69, Issue 2-3
A Novel Approach to Visualizing Dark Matter Simulations
journal, December 2012
- Kaehler, R.; Hahn, O.; Abel, T.
- IEEE Transactions on Visualization and Computer Graphics, Vol. 18, Issue 12
One-to-one direct modeling of experiments and astrophysical scenarios: pushing the envelope on kinetic plasma simulations
journal, November 2008
- Fonseca, R. A.; Martins, S. F.; Silva, L. O.
- Plasma Physics and Controlled Fusion, Vol. 50, Issue 12
Dense Electron-Positron Plasmas and Ultraintense rays from Laser-Irradiated Solids
journal, April 2012
- Ridgers, C. P.; Brady, C. S.; Duclous, R.
- Physical Review Letters, Vol. 108, Issue 16
Elimination of the numerical Cerenkov instability for spectral EM-PIC codes
journal, July 2015
- Yu, Peicheng; Xu, Xinlu; Decyk, Viktor K.
- Computer Physics Communications, Vol. 192
Spontaneously Growing Transverse Waves in a Plasma Due to an Anisotropic Velocity Distribution
journal, February 1959
- Weibel, Erich S.
- Physical Review Letters, Vol. 2, Issue 3
Mechanism for Instability of Transverse Plasma Waves
journal, January 1959
- Fried, Burton D.
- Physics of Fluids, Vol. 2, Issue 3
Seeded QED cascades in counterpropagating laser pulses
journal, February 2017
- Grismayer, T.; Vranic, M.; Martins, J. L.
- Physical Review E, Vol. 95, Issue 2
The properties of cosmic velocity fields
journal, October 2015
- Hahn, Oliver; Angulo, Raul E.; Abel, Tom
- Monthly Notices of the Royal Astronomical Society, Vol. 454, Issue 4
Efficient modeling of laser–plasma interactions in high energy density scenarios
journal, May 2011
- Fiuza, F.; Marti, M.; Fonseca, R. A.
- Plasma Physics and Controlled Fusion, Vol. 53, Issue 7
An exact general remeshing scheme applied to physically conservative voxelization
journal, September 2015
- Powell, Devon; Abel, Tom
- Journal of Computational Physics, Vol. 297
A new approach to simulating collisionless dark matter fluids
journal, July 2013
- Hahn, Oliver; Abel, Tom; Kaehler, Ralf
- Monthly Notices of the Royal Astronomical Society, Vol. 434, Issue 2
Particle acceleration in laser-driven magnetic reconnection
journal, April 2017
- Totorica, S. R.; Abel, T.; Fiuza, F.
- Physics of Plasmas, Vol. 24, Issue 4
Relativistic Electron Streaming Instabilities Modulate Proton Beams Accelerated in Laser-Plasma Interactions
journal, May 2017
- Göde, S.; Rödel, C.; Zeil, K.
- Physical Review Letters, Vol. 118, Issue 19
A method for incorporating Gauss' law into electromagnetic PIC codes
journal, January 1987
- Marder, Barry
- Journal of Computational Physics, Vol. 68, Issue 1
QED cascades induced by circularly polarized laser fields
journal, May 2011
- Elkina, N. V.; Fedotov, A. M.; Kostyukov, I. Yu.
- Physical Review Special Topics - Accelerators and Beams, Vol. 14, Issue 5
A binary collision model for plasma simulation with a particle code
journal, November 1977
- Takizuka, Tomonor; Abe, Hirotada
- Journal of Computational Physics, Vol. 25, Issue 3
ColDICE: A parallel Vlasov–Poisson solver using moving adaptive simplicial tessellation
journal, September 2016
- Sousbie, Thierry; Colombi, Stéphane
- Journal of Computational Physics, Vol. 321
Statistical kinetic treatment of relativistic binary collisions
journal, February 2009
- Peano, F.; Marti, M.; Silva, L. O.
- Physical Review E, Vol. 79, Issue 2
Exploiting multi-scale parallelism for large scale numerical modelling of laser wakefield accelerators
journal, November 2013
- Fonseca, R. A.; Vieira, J.; Fiuza, F.
- Plasma Physics and Controlled Fusion, Vol. 55, Issue 12
Observation of magnetic field generation via the Weibel instability in interpenetrating plasma flows
journal, January 2015
- Huntington, C. M.; Fiuza, F.; Ross, J. S.
- Nature Physics, Vol. 11, Issue 2
Cosmic web, multistream flows, and tessellations
journal, April 2012
- Shandarin, Sergei; Habib, Salman; Heitmann, Katrin
- Physical Review D, Vol. 85, Issue 8
Contemporary particle-in-cell approach to laser-plasma modelling
journal, September 2015
- Arber, T. D.; Bennett, K.; Brady, C. S.
- Plasma Physics and Controlled Fusion, Vol. 57, Issue 11
Controlling the numerical Cerenkov instability in PIC simulations using a customized finite difference Maxwell solver and a local FFT based current correction
journal, May 2017
- Li, Fei; Yu, Peicheng; Xu, Xinlu
- Computer Physics Communications, Vol. 214
Particle Acceleration in Relativistic Collisionless Shocks: Fermi Process at Last?
journal, July 2008
- Spitkovsky, Anatoly
- The Astrophysical Journal, Vol. 682, Issue 1
Relativistic Reconnection: an Efficient Source of Non-Thermal Particles
journal, February 2014
- Sironi, Lorenzo; Spitkovsky, Anatoly
- The Astrophysical Journal, Vol. 783, Issue 1
Numerical stability of relativistic beam multidimensional PIC simulations employing the Esirkepov algorithm
journal, September 2013
- Godfrey, Brendan B.; Vay, Jean-Luc
- Journal of Computational Physics, Vol. 248
An adaptively refined phase–space element method for cosmological simulations and collisionless dynamics
journal, November 2015
- Hahn, Oliver; Angulo, Raul E.
- Monthly Notices of the Royal Astronomical Society, Vol. 455, Issue 1
Tracing the dark matter sheet in phase space: Tracing the dark matter sheet
journal, October 2012
- Abel, Tom; Hahn, Oliver; Kaehler, Ralf
- Monthly Notices of the Royal Astronomical Society, Vol. 427, Issue 1
Numerical instability due to relativistic plasma drift in EM-PIC simulations
journal, November 2013
- Xu, Xinlu; Yu, Peicheng; Martins, Samual F.
- Computer Physics Communications, Vol. 184, Issue 11
Weighted Particles in Coulomb Collision Simulations Based on the Theory of a Cumulative Scattering Angle
journal, September 1998
- Nanbu, K.; Yonemura, S.
- Journal of Computational Physics, Vol. 145, Issue 2
Generation of Magnetic Fields in the Relativistic Shock of Gamma‐Ray Burst Sources
journal, December 1999
- Medvedev, Mikhail V.; Loeb, Abraham
- The Astrophysical Journal, Vol. 526, Issue 2
Interpenetrating Plasma Shells: Near-Equipartition Magnetic Field Generation and Nonthermal Particle Acceleration
journal, September 2003
- Silva, L. O.; Fonseca, R. A.; Tonge, J. W.
- The Astrophysical Journal, Vol. 596, Issue 1
Numerical methods for particle simulations at extreme densities and temperatures: Weighted particles, relativistic collisions and reduced currents
journal, July 2008
- Sentoku, Y.; Kemp, A. J.
- Journal of Computational Physics, Vol. 227, Issue 14
Particle merging algorithm for PIC codes
journal, June 2015
- Vranic, M.; Grismayer, T.; Martins, J. L.
- Computer Physics Communications, Vol. 191
Improved modeling of relativistic collisions and collisional ionization in particle-in-cell codes
journal, August 2012
- Pérez, F.; Gremillet, L.; Decoster, A.
- Physics of Plasmas, Vol. 19, Issue 8
Nonthermal Electron Energization from Magnetic Reconnection in Laser-Driven Plasmas
journal, March 2016
- Totorica, Samuel R.; Abel, Tom; Fiuza, Frederico
- Physical Review Letters, Vol. 116, Issue 9
Low-noise electromagnetic and relativistic particle-in-cell plasma simulation models
journal, September 1999
- Sydora, R. D.
- Journal of Computational and Applied Mathematics, Vol. 109, Issue 1-2
Computational Science — ICCS 2002: International Conference Amsterdam, The Netherlands, April 21–24, 2002 Proceedings, Part III
book, January 2002
- Goos, Gerhard; Hartmanis, Juris; van Leeuwen, Jan
- Lecture Notes in Computer Science
Numerical Cherenkov instabilities in electromagnetic particle codes
journal, August 1974
- Godfrey, Brendan B.
- Journal of Computational Physics, Vol. 15, Issue 4
Exact charge conservation scheme for Particle-in-Cell simulation with an arbitrary form-factor
journal, April 2001
- Esirkepov, T. Zh.
- Computer Physics Communications, Vol. 135, Issue 2
Weibel-Instability-Mediated Collisionless Shocks in the Laboratory with Ultraintense Lasers
journal, June 2012
- Fiuza, F.; Fonseca, R. A.; Tonge, J.
- Physical Review Letters, Vol. 108, Issue 23
Particle simulation of plasmas
journal, April 1983
- Dawson, John M.
- Reviews of Modern Physics, Vol. 55, Issue 2
Simplex-in-cell technique for collisionless plasma simulations
journal, January 2016
- Kates-Harbeck, Julian; Totorica, Samuel; Zrake, Jonathan
- Journal of Computational Physics, Vol. 304