Avalanche-to-streamer transition in particle simulations
Citation for published version (APA):Li, C., Ebert, U., & Brok, W. J. M. (2008). Avalanche-to-streamer transition in particle simulations. IEEE Transactions on Plasma Science, 36(4), 910-911. https://doi.org/10.1109/TPS.2008.922487
DOI:
10.1109/TPS.2008.922487 Document status and date: Published: 01/01/2008
Document Version:
Publisher’s PDF, also known as Version of Record (includes final page, issue and volume numbers)
Please check the document version of this publication:
• A submitted manuscript is the version of the article upon submission and before peer-review. There can be important differences between the submitted version and the official published version of record. People interested in the research are advised to contact the author for the final version of the publication, or visit the DOI to the publisher's website.
• The final author version and the galley proof are versions of the publication after peer review.
• The final published version features the final layout of the paper including the volume, issue and page numbers.
Link to publication
General rights
Copyright and moral rights for the publications made accessible in the public portal are retained by the authors and/or other copyright owners and it is a condition of accessing publications that users recognise and abide by the legal requirements associated with these rights. • Users may download and print one copy of any publication from the public portal for the purpose of private study or research. • You may not further distribute the material or use it for any profit-making activity or commercial gain
• You may freely distribute the URL identifying the publication in the public portal.
If the publication is distributed under the terms of Article 25fa of the Dutch Copyright Act, indicated by the “Taverne” license above, please follow below link for the End User Agreement:
www.tue.nl/taverne Take down policy
If you believe that this document breaches copyright please contact us at: openaccess@tue.nl
providing details and we will investigate your claim.
910 IEEE TRANSACTIONS ON PLASMA SCIENCE, VOL. 36, NO. 4, AUGUST 2008
Avalanche-to-Streamer Transition
in Particle Simulations
Chao Li, Ute Ebert, and W. J. M. Brok
Abstract—The avalanche-to-streamer transition is studied and
illustrated in a particle model. The results are similar to those of fluid models. However, when superparticles are introduced, numerical artifacts become visible. This underscores the need of models that are hybrid in space.
Index Terms—Particle simulation, streamers.
S
TREAMERS are growing filaments of weakly ionized nonstationary plasma produced by a sharp ionization front that propagates into nonionized matter. Streamers are used in industrial applications such as lighting [1], ozone generation, and gas and water purification, and they occur in natural processes such as lightning and transient luminous events in the upper atmosphere [2]. Streamers can emerge from ionization avalanches in free space when the self-induced field becomes comparable to the applied field; this was recently reinvestigated in the framework of a fluid model in [3]. Reference [4] dis-cussed that fluid models cannot be applied to avalanches in noble gas in a very low field; and in [5], striations were found in this case. In nitrogen and in intermediate fields, particle and fluid models give essentially the same results for planar ionization fronts, whereas there are growing deviations again at larger fields as discussed in [6].Here, we show how a streamer emerges in a particle model from few seed electrons and a consecutive avalanche. The system is pure nitrogen at standard temperature and pressure in a constant high field of 100 kV/cm (equivalent to 42 Td). We assume that the electrodes are far away or do not exist as in many natural discharges. Qualitatively similar results in a fluid model are shown in [3, Fig. 1].
The particle model includes the complete electron veloc-ity and energy distribution, as well as the discrete nature of particles. However, the required computation resources grow with the number of particles and eventually exceed the limits of any computer. This difficulty can be counteracted by us-ing superparticles carryus-ing the charge and the mass of many physical particles, but superparticles in turn create unphysical
Manuscript received December 1, 2007; revised February 14, 2008. This work was supported by the Dutch national program BSIK in the ICT project BRICKS (theme MSV1).
C. Li is with the Center for Mathematics and Computer Science (CWI), 1090 GB Amsterdam, The Netherlands.
U. Ebert is with the Center for Mathematics and Computer Science (CWI), 1090 GB Amsterdam, The Netherlands, and also with the Department of Applied Physics, Eindhoven University of Technology, 5600 MB Eindhoven, The Netherlands.
W. J. M. Brok is with the Department of Applied Physics, Eindhoven University of Technology, 5600 MB Eindhoven, The Netherlands.
Digital Object Identifier 10.1109/TPS.2008.922487
fluctuations and stochastic heating, as we will demonstrate hereafter.
The simulation starts from 100 pairs of electrons and ions at one point and follows the initial particles and their offspring up to 0.36-ns time. In Fig. 1, the simulated streamer is shown at two different times: t = 0.18 ns with real particles (first row),
t = 0.36 ns with real particles (second row), and once again t = 0.36 ns, now in a superparticle simulation (third row).
At time t = 0.18 ns (first row), we have approximately 105 electrons, and the maximal field enhancement is ∼1%. Therefore, the discharge is still in the avalanche phase. For the second row, de computation is continued with real particles, whereas for the third row, superparticles are introduced as follows. If a specified number of particles, in our case 105, is reached, a particle remapping is applied to reduce the number of computational particles; in this step, half of them are thrown away at random, and the weight of the remaining computational particles is doubled. At t = 0.36 ns, the total number of elec-trons is roughly 1.5× 107 within the real particle simulation. The space-charge layer has clearly formed. The maximal field enhancement is∼50%. At t = 0.36 ns within the superparticle simulation, the total number of superparticles is about 7× 104, each representing 256 real electrons. No charge layer has formed, but the charge density is noisy with a maximum inside the discharge. Clearly, superparticles here are not appropriate to deal computationally with large electron numbers. In [7], we describe how to circumvent the unwieldy runtimes of real particle models by using a hybrid model which couples particle and fluid models in different spatial regions.
REFERENCES
[1] A. Bhoj and M. Kushner, “Plasma dynamics during breakdown in an HID lamp,” IEEE Trans. Plasma Sci., vol. 33, no. 2, pp. 518–519, Apr. 2005.
[2] V. P. Pasko, “Red sprite discharges in the atmosphere at high altitude: The molecular physics and the similarity with laboratory discharges,” Plasma
Sources Sci. Technol., vol. 16, no. 1, pp. S13–S29, Feb. 2007.
[3] C. Montijn and U. Ebert, “Diffusion correction to the Raether–Meek cri-terion for the avalanche-to-streamer transition,” J. Phys. D, Appl. Phys., vol. 39, no. 14, pp. 2979–2992, Jul. 2006.
[4] B. J. P. Dowds, R. K. Barrett, and D. A. Diver, “Streamer initiation in atmospheric pressure gas discharges by direct particle simulation,” Phys.
Rev. E, Stat. Phys. Plasmas Fluids Relat. Interdiscip. Top., vol. 68, no. 2,
p. 026 412, Aug. 2003.
[5] W. J. M. Brok, “Modelling of transient phenomena in gas discharges,” Ph.D. dissertation, Eindhoven Univ. Tech., Eindhoven, The Netherlands, 2005. [Online]. Available: http://alexandria.tue.nl/extra2/200512799.pdf [6] C. Li, W. J. M. Brok, U. Ebert, and J. J. A. M. van der Mullen, “Deviations
from the local field approximation in negative streamer heads,” J. Appl.
Phys., vol. 101, no. 12, p. 123 305, Jun. 2007.
[7] C. Li, U. Ebert, W. J. M. Brok, and W. Hundsdorfer, “Spatial coupling of particle and fluid models for streamers: where nonlocality matters,” J. Phys.
D, Appl. Phys., vol. 41, no. 3, p. 032 005, Feb. 2008.
LI et al.: AVALANCHE-TO-STREAMER TRANSITION IN PARTICLE SIMULATIONS 911
Fig. 1. Avalanche to streamer transition in nitrogen at standard temperature and pressure. A constant background electric field of 100 kV/cm points downwards. First row: avalanche of real particles at t = 0.18 ns, second row: streamer formed by real particles at t = 0.36 ns, third row: the same system at the same time
t = 0.36 ns, but using super-particles with the mass of 256 real particles. The columns show from left to right: electron density, charge density, and electric field