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

Title: Timesaving techniques for decision of electron-molecule collisions in Monte Carlo simulation of electrical discharges

Abstract

Techniques to reduce the computational load for determination of electron-molecule collisions in Monte Carlo simulations of electrical discharges have been presented. By enhancing the detection efficiency of the no-collision case in the decision scheme of the collisional events, we can decrease the frequency of access to time-consuming subroutines to calculate the electron collision cross sections of the gas molecules for obtaining the collision probability. A benchmark test and an estimation to evaluate the present techniques have shown a practical timesaving efficiency.

Authors:
 [1];  [2];  [2];  [2]
  1. Division of Electronics for Informatics, Graduate School of Information Science and Technology, Hokkaido University, North 14 West 9, Hokkaido, Sapporo 060-0814 (Japan). E-mail: sugawara@ist.hokudai.ac.jp
  2. Division of Electronics for Informatics, Graduate School of Information Science and Technology, Hokkaido University, North 14 West 9, Hokkaido, Sapporo 060-0814 (Japan)
Publication Date:
OSTI Identifier:
20991573
Resource Type:
Journal Article
Resource Relation:
Journal Name: Journal of Computational Physics; Journal Volume: 223; Journal Issue: 1; Other Information: DOI: 10.1016/j.jcp.2006.09.007; PII: S0021-9991(06)00437-2; Copyright (c) 2006 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved; Country of input: International Atomic Energy Agency (IAEA)
Country of Publication:
United States
Language:
English
Subject:
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; ALGORITHMS; BENCHMARKS; COMPUTERIZED SIMULATION; CROSS SECTIONS; ELECTRIC DISCHARGES; ELECTRON-MOLECULE COLLISIONS; MONTE CARLO METHOD; PLASMA; PROBABILITY

Citation Formats

Sugawara, Hirotake, Mori, Naoki, Sakai, Yosuke, and Suda, Yoshiyuki. Timesaving techniques for decision of electron-molecule collisions in Monte Carlo simulation of electrical discharges. United States: N. p., 2007. Web. doi:10.1016/j.jcp.2006.09.007.
Sugawara, Hirotake, Mori, Naoki, Sakai, Yosuke, & Suda, Yoshiyuki. Timesaving techniques for decision of electron-molecule collisions in Monte Carlo simulation of electrical discharges. United States. doi:10.1016/j.jcp.2006.09.007.
Sugawara, Hirotake, Mori, Naoki, Sakai, Yosuke, and Suda, Yoshiyuki. Tue . "Timesaving techniques for decision of electron-molecule collisions in Monte Carlo simulation of electrical discharges". United States. doi:10.1016/j.jcp.2006.09.007.
@article{osti_20991573,
title = {Timesaving techniques for decision of electron-molecule collisions in Monte Carlo simulation of electrical discharges},
author = {Sugawara, Hirotake and Mori, Naoki and Sakai, Yosuke and Suda, Yoshiyuki},
abstractNote = {Techniques to reduce the computational load for determination of electron-molecule collisions in Monte Carlo simulations of electrical discharges have been presented. By enhancing the detection efficiency of the no-collision case in the decision scheme of the collisional events, we can decrease the frequency of access to time-consuming subroutines to calculate the electron collision cross sections of the gas molecules for obtaining the collision probability. A benchmark test and an estimation to evaluate the present techniques have shown a practical timesaving efficiency.},
doi = {10.1016/j.jcp.2006.09.007},
journal = {Journal of Computational Physics},
number = 1,
volume = 223,
place = {United States},
year = {Tue Apr 10 00:00:00 EDT 2007},
month = {Tue Apr 10 00:00:00 EDT 2007}
}
  • Electron properties in a parallel plate capacitively coupled rf discharge are studied with results from a Monte-Carlo simulation. Time averaged, spatially dependent electron distributions are computed by integrating, in time, electron trajectories as a function of position while oscillating the applied electric field at rf frequencies. The dc component of the sheath potential is solved for in a self-consistent manner during the simulation. For conditions where the secondary emission coefficient for electrons from the electrodes is large, the electron distribution is spatially differentiated, being dominated by an e-beam component near the electrodes while being nearly in equilibrium with the appliedmore » electric field in the body of the plasma. The dc component of the sheath potential is found to be a function of the ratio lambda/d, where lambda is the electron mean free path and d is the electrode spacing.« less
  • A model of one dimensional in position and three dimensional in velocity space self-consistent particle in cell with Monte Carlo collision technique was employed to simulate the argon discharge between the needle and plane electrodes at high pressure, in which a nanosecond rectangular pulse was applied to the needle electrode. The work focused on the investigation of the spatiotemporal evolution of the discharge versus the needle tip size and working gas pressure. The simulation results showed that the discharge occurred mainly in the region near the needle tip at atmospheric pressure, and that the small radius of the needle tipmore » led to easy discharge. Reducing the gas pressure gave rise to a transition from a corona discharge to a glowlike discharge along the needle-to-plane direction. The microscopic mechanism for the transition can arguably be attributed to the peak of high-energy electrons occurring before the breakdown; the magnitude of the number of these electrons determined whether the breakdown can take place.« less
  • Although surface-wave plasma (SWP) sources have many industrial applications, the ionization process for SWP discharges is not yet well understood. The resonant excitation of surface plasmon polaritons (SPPs) has recently been proposed to produce SWP efficiently, and this work presents a numerical study of the mechanism to produce SWP sources. Specifically, SWP resonantly excited by SPPs at low pressure (0.25 Torr) are modeled using a two-dimensional in the working space and three-dimensional in the velocity space particle-in-cell with the Monte Carlo collision method. Simulation results are sampled at different time steps, in which the detailed information about the distribution ofmore » electrons and electromagnetic fields is obtained. Results show that the mode conversion between surface waves of SPPs and electron plasma waves (EPWs) occurs efficiently at the location where the plasma density is higher than 3.57 Multiplication-Sign 10{sup 17} m{sup -3}. Due to the effect of the locally enhanced electric field of SPPs, the mode conversion between the surface waves of SPPs and EPWs is very strong, which plays a significant role in efficiently heating SWP to the overdense state.« less
  • The radio-frequency (RF) plasma production technique in the ion cyclotron range of frequency (ICRF) attracts growing attention among fusion experts because of its high potential for solving several basic problems of reactor-oriented superconducting fusion machines, such as ICRF wall conditioning in tokamaks and stellarators (T{sub e} = 3−5eV, n{sub e}<10{sup 12}cm{sup −3}), ICRF-assisted tokamak start-up and target plasma production (n{sub e} = 10{sup 13}cm{sup −3}) in stellarators. Plasma initiation by ICRF has been studied intensively using single particle descriptions and basic analytic models. To further improve the present understanding on plasma production employing the vacuum RF field of ICRF antennas in toroidal devices inmore » presence of the toroidal magnetic field, and its parametric dependencies a Monte Carlo code has been developed. The 1D code RFdinity1D describes the motion of electrons, accelerated by the RF field in front of the ICRF antenna, along one toroidal magnetic field line. Dependent on their individual energies and the related electron collision cross sections (ionisation, excitation and dissociation) weighted by a Monte Carlo procedure, an electron avalanche may occur. Breakdown conditions are discussed as function of RF discharge parameters (i) RF vacuum electric field strength, (ii) RF frequency and (iii) neutral pressure (H2). The slope of the exponential density increase, taken as measure for the breakdown speed, shows qualitative agreement to experimental breakdown times as found in literature and experimental data of the ASDEX upgrade and TEXTOR tokamak, and is interpreted by studying the characteristic electron velocity distribution functions.« less
  • The problems connected with the simulation by the diffusion quantum Monte Carlo method of atomic and molecular systems in the two-dimensional (2D) space have been investigated on the 1 {sup 2}S and 2{sup 2}S states of the hydrogen atom and on the 1{sup 1}S and 1 {sup 3}S states of the helium atom, assuming {tau}{sup {minus}1} Coulomb interactions between the particles. The potential surface of the {sup 1}{Sigma}{sub g}{sup +} state of the 2D hydrogen molecule has been calculated in the range 0.1-1.5 bohr: A dissociation energy of 1.2703 Hartrees is found at the equilibrium distance of 0.3639 bohr. 33more » refs., 4 figs., 3 tabs.« less