Regular sensitivity computation avoiding chaotic effects in particle-in-cell plasma methods

Seung Whan Chung, Stephen D. Bond, Eric C. Cyr, Jonathan B. Freund

Research output: Contribution to journalArticle

Abstract

Particle-in-cell (PIC) simulation methods are attractive for representing species distribution functions in plasmas. However, as a model, they introduce uncertain parameters, and for quantifying their prediction uncertainty it is useful to be able to assess the sensitivity of a quantity-of-interest (QoI) to these parameters. Such sensitivity information is likewise useful for optimization. However, computing sensitivity for PIC methods is challenging due to the chaotic particle dynamics, and sensitivity techniques remain underdeveloped compared to those for Eulerian discretizations. This challenge is examined from a dual particle–continuum perspective that motivates a new sensitivity discretization. Two routes to sensitivity computation are presented and compared: a direct fully-Lagrangian particle-exact approach provides sensitivities of each particle trajectory, and a new particle-pdf discretization, which is formulated from a continuum perspective but discretized by particles to take the advantages of the same type of Lagrangian particle description leveraged by PIC methods. Since the sensitivity particles in this approach are only indirectly linked to the plasma-PIC particles, they can be positioned and weighted independently for efficiency and accuracy. The corresponding numerical algorithms are presented in mathematical detail. The advantage of the particle-pdf approach in avoiding the spurious chaotic sensitivity of the particle-exact approach is demonstrated for Debye shielding and sheath configurations. In essence, the continuum perspective makes implicit the distinctness of the particles, which circumvents the Lyapunov instability of the N-body PIC system. The cost of the particle-pdf approach is comparable to the baseline PIC simulation.

Original languageEnglish (US)
Article number108969
JournalJournal of Computational Physics
Volume400
DOIs
StatePublished - Jan 1 2020

Fingerprint

Plasma
Plasmas
sensitivity
Cell
cells
Shielding
Distribution functions
Trajectories
Costs
Discretization
Continuum
continuums
Uncertainty
Particle Trajectory
particle trajectories
Uncertain Parameters
sheaths
Numerical Algorithms
Simulation Methods
Lyapunov

Keywords

  • Chaos
  • Particle-in-cell methods
  • Plasma models
  • Sensitivity
  • Uncertainty quantification

ASJC Scopus subject areas

  • Numerical Analysis
  • Modeling and Simulation
  • Physics and Astronomy (miscellaneous)
  • Physics and Astronomy(all)
  • Computer Science Applications
  • Computational Mathematics
  • Applied Mathematics

Cite this

Regular sensitivity computation avoiding chaotic effects in particle-in-cell plasma methods. / Chung, Seung Whan; Bond, Stephen D.; Cyr, Eric C.; Freund, Jonathan B.

In: Journal of Computational Physics, Vol. 400, 108969, 01.01.2020.

Research output: Contribution to journalArticle

@article{859989dc68e64e008d4073902631d83b,
title = "Regular sensitivity computation avoiding chaotic effects in particle-in-cell plasma methods",
abstract = "Particle-in-cell (PIC) simulation methods are attractive for representing species distribution functions in plasmas. However, as a model, they introduce uncertain parameters, and for quantifying their prediction uncertainty it is useful to be able to assess the sensitivity of a quantity-of-interest (QoI) to these parameters. Such sensitivity information is likewise useful for optimization. However, computing sensitivity for PIC methods is challenging due to the chaotic particle dynamics, and sensitivity techniques remain underdeveloped compared to those for Eulerian discretizations. This challenge is examined from a dual particle–continuum perspective that motivates a new sensitivity discretization. Two routes to sensitivity computation are presented and compared: a direct fully-Lagrangian particle-exact approach provides sensitivities of each particle trajectory, and a new particle-pdf discretization, which is formulated from a continuum perspective but discretized by particles to take the advantages of the same type of Lagrangian particle description leveraged by PIC methods. Since the sensitivity particles in this approach are only indirectly linked to the plasma-PIC particles, they can be positioned and weighted independently for efficiency and accuracy. The corresponding numerical algorithms are presented in mathematical detail. The advantage of the particle-pdf approach in avoiding the spurious chaotic sensitivity of the particle-exact approach is demonstrated for Debye shielding and sheath configurations. In essence, the continuum perspective makes implicit the distinctness of the particles, which circumvents the Lyapunov instability of the N-body PIC system. The cost of the particle-pdf approach is comparable to the baseline PIC simulation.",
keywords = "Chaos, Particle-in-cell methods, Plasma models, Sensitivity, Uncertainty quantification",
author = "Chung, {Seung Whan} and Bond, {Stephen D.} and Cyr, {Eric C.} and Freund, {Jonathan B.}",
year = "2020",
month = "1",
day = "1",
doi = "10.1016/j.jcp.2019.108969",
language = "English (US)",
volume = "400",
journal = "Journal of Computational Physics",
issn = "0021-9991",
publisher = "Academic Press Inc.",

}

TY - JOUR

T1 - Regular sensitivity computation avoiding chaotic effects in particle-in-cell plasma methods

AU - Chung, Seung Whan

AU - Bond, Stephen D.

AU - Cyr, Eric C.

AU - Freund, Jonathan B.

PY - 2020/1/1

Y1 - 2020/1/1

N2 - Particle-in-cell (PIC) simulation methods are attractive for representing species distribution functions in plasmas. However, as a model, they introduce uncertain parameters, and for quantifying their prediction uncertainty it is useful to be able to assess the sensitivity of a quantity-of-interest (QoI) to these parameters. Such sensitivity information is likewise useful for optimization. However, computing sensitivity for PIC methods is challenging due to the chaotic particle dynamics, and sensitivity techniques remain underdeveloped compared to those for Eulerian discretizations. This challenge is examined from a dual particle–continuum perspective that motivates a new sensitivity discretization. Two routes to sensitivity computation are presented and compared: a direct fully-Lagrangian particle-exact approach provides sensitivities of each particle trajectory, and a new particle-pdf discretization, which is formulated from a continuum perspective but discretized by particles to take the advantages of the same type of Lagrangian particle description leveraged by PIC methods. Since the sensitivity particles in this approach are only indirectly linked to the plasma-PIC particles, they can be positioned and weighted independently for efficiency and accuracy. The corresponding numerical algorithms are presented in mathematical detail. The advantage of the particle-pdf approach in avoiding the spurious chaotic sensitivity of the particle-exact approach is demonstrated for Debye shielding and sheath configurations. In essence, the continuum perspective makes implicit the distinctness of the particles, which circumvents the Lyapunov instability of the N-body PIC system. The cost of the particle-pdf approach is comparable to the baseline PIC simulation.

AB - Particle-in-cell (PIC) simulation methods are attractive for representing species distribution functions in plasmas. However, as a model, they introduce uncertain parameters, and for quantifying their prediction uncertainty it is useful to be able to assess the sensitivity of a quantity-of-interest (QoI) to these parameters. Such sensitivity information is likewise useful for optimization. However, computing sensitivity for PIC methods is challenging due to the chaotic particle dynamics, and sensitivity techniques remain underdeveloped compared to those for Eulerian discretizations. This challenge is examined from a dual particle–continuum perspective that motivates a new sensitivity discretization. Two routes to sensitivity computation are presented and compared: a direct fully-Lagrangian particle-exact approach provides sensitivities of each particle trajectory, and a new particle-pdf discretization, which is formulated from a continuum perspective but discretized by particles to take the advantages of the same type of Lagrangian particle description leveraged by PIC methods. Since the sensitivity particles in this approach are only indirectly linked to the plasma-PIC particles, they can be positioned and weighted independently for efficiency and accuracy. The corresponding numerical algorithms are presented in mathematical detail. The advantage of the particle-pdf approach in avoiding the spurious chaotic sensitivity of the particle-exact approach is demonstrated for Debye shielding and sheath configurations. In essence, the continuum perspective makes implicit the distinctness of the particles, which circumvents the Lyapunov instability of the N-body PIC system. The cost of the particle-pdf approach is comparable to the baseline PIC simulation.

KW - Chaos

KW - Particle-in-cell methods

KW - Plasma models

KW - Sensitivity

KW - Uncertainty quantification

UR - http://www.scopus.com/inward/record.url?scp=85073379087&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85073379087&partnerID=8YFLogxK

U2 - 10.1016/j.jcp.2019.108969

DO - 10.1016/j.jcp.2019.108969

M3 - Article

AN - SCOPUS:85073379087

VL - 400

JO - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

M1 - 108969

ER -