skip to main content
DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: 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 » 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.« less

Authors:
ORCiD logo [1]; ORCiD logo [2];  [3]
  1. 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
  2. SLAC National Accelerator Lab., Menlo Park, CA (United States). High Energy Density Science Division
  3. 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) (SC-24); 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: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: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}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Save / Share:

Works referenced in this record:

Rigorous charge conservation for local electromagnetic field solvers
journal, March 1992


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
  • DOI: 10.1109/TVCG.2012.187

One-to-one direct modeling of experiments and astrophysical scenarios: pushing the envelope on kinetic plasma simulations
journal, November 2008


Dense Electron-Positron Plasmas and Ultraintense γ rays from Laser-Irradiated Solids
journal, April 2012


Elimination of the numerical Cerenkov instability for spectral EM-PIC codes
journal, July 2015


Spontaneously Growing Transverse Waves in a Plasma Due to an Anisotropic Velocity Distribution
journal, February 1959


Mechanism for Instability of Transverse Plasma Waves
journal, January 1959


Seeded QED cascades in counterpropagating laser pulses
journal, February 2017


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
  • DOI: 10.1093/mnras/stv2179

Efficient modeling of laser–plasma interactions in high energy density scenarios
journal, May 2011


An exact general remeshing scheme applied to physically conservative voxelization
journal, September 2015


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
  • DOI: 10.1093/mnras/stt1061

Particle acceleration in laser-driven magnetic reconnection
journal, April 2017

  • Totorica, S. R.; Abel, T.; Fiuza, F.
  • Physics of Plasmas, Vol. 24, Issue 4
  • DOI: 10.1063/1.4978627

A method for incorporating Gauss' law into electromagnetic PIC codes
journal, January 1987


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
  • DOI: 10.1103/PhysRevSTAB.14.054401

A binary collision model for plasma simulation with a particle code
journal, November 1977


ColDICE: A parallel Vlasov–Poisson solver using moving adaptive simplicial tessellation
journal, September 2016


Statistical kinetic treatment of relativistic binary collisions
journal, February 2009


Exploiting multi-scale parallelism for large scale numerical modelling of laser wakefield accelerators
journal, November 2013


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
  • DOI: 10.1038/nphys3178

Cosmic web, multistream flows, and tessellations
journal, April 2012


Contemporary particle-in-cell approach to laser-plasma modelling
journal, September 2015


Particle Acceleration in Relativistic Collisionless Shocks: Fermi Process at Last?
journal, July 2008

  • Spitkovsky, Anatoly
  • The Astrophysical Journal, Vol. 682, Issue 1
  • DOI: 10.1086/590248

Relativistic Reconnection: an Efficient Source of Non-Thermal Particles
journal, February 2014


Numerical stability of relativistic beam multidimensional PIC simulations employing the Esirkepov algorithm
journal, September 2013


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
  • DOI: 10.1093/mnras/stv2304

Tracing the dark matter sheet in phase space: Tracing the dark matter sheet
journal, October 2012


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
  • DOI: 10.1016/j.cpc.2013.07.003

Weighted Particles in Coulomb Collision Simulations Based on the Theory of a Cumulative Scattering Angle
journal, September 1998


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
  • DOI: 10.1086/308038

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
  • DOI: 10.1086/379156

Particle merging algorithm for PIC codes
journal, June 2015


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
  • DOI: 10.1063/1.4742167

Nonthermal Electron Energization from Magnetic Reconnection in Laser-Driven Plasmas
journal, March 2016


Low-noise electromagnetic and relativistic particle-in-cell plasma simulation models
journal, September 1999


Numerical Cherenkov instabilities in electromagnetic particle codes
journal, August 1974


Exact charge conservation scheme for Particle-in-Cell simulation with an arbitrary form-factor
journal, April 2001


Weibel-Instability-Mediated Collisionless Shocks in the Laboratory with Ultraintense Lasers
journal, June 2012


Particle simulation of plasmas
journal, April 1983


Simplex-in-cell technique for collisionless plasma simulations
journal, January 2016

  • Kates-Harbeck, Julian; Totorica, Samuel; Zrake, Jonathan
  • Journal of Computational Physics, Vol. 304
  • DOI: 10.1016/j.jcp.2015.10.017