Validation and benchmarking of two particle-in-cell codes for a glow discharge

Carlsson, J.; Khrabrov, A.; Kaganovich, I.; Sommerer, T.; Keating, D.
Issue date: 2017
Cite as:
Carlsson, J., Khrabrov, A., Kaganovich, I., Sommerer, T., & Keating, D. (2016). Validation and benchmarking of two particle-in-cell codes for a glow discharge [Data set]. Princeton Plasma Physics Laboratory, Princeton University.
@electronic{carlsson_j_2016,
  author      = {Carlsson, J. and
                Khrabrov, A. and
                Kaganovich, I. and
                Sommerer, T. and
                Keating, D.},
  title       = {{Validation and benchmarking of two parti
                cle-in-cell codes for a glow discharge}},
  publisher   = {{Princeton Plasma Physics Laboratory, Pri
                nceton University}},
  year        = 2016
}
Abstract:

The two particle-in-cell codes EDIPIC and LSP are benchmarked and validated for a parallel-plate glow discharge in helium, in which the axial electric field had been carefully measured, primarily to investigate and improve the fidelity of their collision models. The scattering anisotropy of electron-impact ionization, as well as the value of the secondary-electron emission yield, are not well known in this case. The experimental uncertainty for the emission yield corresponds to a factor of two variation in the cathode current. If the emission yield is tuned to make the cathode current computed by each code match the experiment, the computed electric fields are in excellent agreement with each other, and within about 10% of the experimental value. The non-monotonic variation of the width of the cathode fa ll with the applied voltage seen in the experiment is reproduced by both codes. The electron temperature in the negative glow is within experimental error bars for both codes, but the density of slow trapped electrons is underestimated. A more detailed code comparison don e for several synthetic cases of electron-beam injection into helium gas shows that the codes are in excellent agreement for ionization rate, as well as for elastic and excitation collisions with isotropic scattering pattern. The remaining significant discrepancies between the two codes are due to differences in their electron binary-collision models, and for anisotropic scattering due to elastic and excitation collisions.

Show More
# Filename Description Filesize
1 DenHartogGlowDischargeElectricField.txt 464 Bytes
2 fig05data.txt 20 KB
3 fig05plot.py 1.23 KB
4 fig06data.txt 20.1 KB
5 fig06plot.py 1.23 KB
6 fig07data.txt 19.3 KB
7 fig07plot.py 1.21 KB
8 fig09data.txt 27.6 KB
9 fig09plot.py 906 Bytes
10 fig10data.txt 29 KB
11 fig10plot.py 906 Bytes
12 fig11data.txt 29 KB
13 fig11plot.py 906 Bytes
14 intereHe4-consolidated.tab 25.6 KB
15 intereHe4.tab 66.8 KB
16 readme.txt 2.33 KB