Cathode-directed streamer evolution in near atmospheric air is modeled in 3D pin-to-plane geometries using a 3D kinetic Particle-In-Cell (PIC) code that simulates particle-particle collisions via the Direct Simulation Monte Carlo (DSMC) method. Due to the computational challenges associated with a complete 360° volumetric domain, a practical alternative was achieved using a wedge domain and a range of azimuthal angles was explored (5°, 15°, 30°, and 45°) to study possible effects on the streamer growth and propagation due to the finite wedge angle. A DC voltage of 6 kV is administered to a hemispherical anode of radius 100 μm, with a planar cathode held at ground potential, generating an over-volted state with an electric field of 4 MV/m across a 1500 μm gap. The domain is seeded with an initial ion and electron density of 1018 m-3 at 1 eV temperature confined to a spherical region of radius 100 μm centered at the tip of the anode. The air chemistry model [1] includes standard Townsend breakdown mechanisms (electron-neutral elastic, excitation, ionization, attachment, and detachment collision chemistry and secondary electron emission) as well as streamer mechanisms (photoionization and ion-neutral collisions) via tracking excited state neutrals which can then either quench via collisions or spontaneously emit a photon based on specific Einstein-A coefficients [2, 3]. In this work, positive streamer dynamics are formally quantified for each wedge angle in terms of electron velocity and density as temporal functions of coordinates r, Φ, and z. Applying a random plasma seed for each simulation, particles of interest are tracked with near femtosecond temporal resolution out to 1.4 ns and spatially binned. This process is repeated six times and results are averaged. Prior 2D studies have shown that the reduced electric field, E/n, can significantly impact streamer evolution [4]. We extend the analysis to 3D wedge geometries, to limit computational costs, and examine the wedge angle’s effect on streamer branching, propagation, and velocity. Results indicate that the smallest wedge angle that produced an acceptably converged solution is 30°. The potential effects that a mesh, when under-resolved with respect to the Debye length, can impart on streamer dynamics and numerical heating were not investigated, and we explicitly state here that the smallest cell size was approximately 10 times the minimum λD in the streamer channel at late times. This constraint on cell size was the result of computational limitations on total mesh count.
This report documents the progress made in simulating the HERMES-III Magnetically Insulated Transmission Line (MITL) and courtyard with EMPIRE and ITS. This study focuses on the shots that were taken during the months of June and July of 2019 performed with the new MITL extension. There were a few shots where there was dose mapping of the courtyard, 11132, 11133, 11134, 11135, 11136, and 11146. This report focuses on these shots because there was full data return from the MITL electrical diagnostics and the radiation dose sensors in the courtyard. The comparison starts with improving the processing of the incoming voltage into the EMPIRE simulation from the experiment. The currents are then compared at several location along the MITL. The simulation results of the electrons impacting the anode are shown. The electron impact energy and angle is then handed off to ITS which calculates the dose on the faceplate and locations in the courtyard and they are compared to experimental measurements. ITS also calculates the photons and electrons that are injected into the courtyard, these quantities are then used by EMPIRE to calculated the photon and electron transport in the courtyard. The details for the algorithms used to perform the courtyard simulations are presented as well as qualitative comparisons of the electric field, magnetic field, and the conductivity in the courtyard. Because of the computational burden of these calculations the pressure was reduce in the courtyard to reduce the computational load. The computation performance is presented along with suggestion on how to improve both the computational performance as well as the algorithmic performance. Some of the algorithmic changed would reduce the accuracy of the models and detail comparison of these changes are left for a future study. As well as, list of code improvements there is also a list of suggested experimental improvements to improve the quality of the data return.
Poondla, Yasvanth; Goldstein, David; Varghese, Philip; Clarke, Peter; Moore, Christopher H.
The goal of this work is to build up the capability of quasi-particle simulation (QuiPS), a novel flow solver, such that it can adequately model the rarefied portion of an atmospheric reentry trajectory. Direct simulation Monte Carlo (DSMC) is the conventional solver for such conditions, but struggles to resolve transient flows, trace species, and high-level internal energy states due to stochastic noise. Quasi-particle simulation (QuiPS) is a novel Boltzmann solver that describes a system with a discretized, truncated velocity distribution function. The resulting fixed-velocity, variable weight quasi-particles enable smooth variation of macroscopic properties. The distribution function description enables the use of a variance-reduced collision model, greatly minimizing expense near equilibrium. This work presents the addition of a neutral air chemistry model to QuiPS and some demonstrative 0D simulations. The explicit representation of internal distributions in QuiPS reveals some of the flaws in existing physics models. Variance reduction, a key feature of QuiPS, can greatly reduce expense of multi-dimensional calculations, but is only cheaper when the gas composition is near chemical equilibrium.
Understanding the role of physical processes contributing to breakdown is critical for many applications in which breakdown is undesirable, such as capacitors, and applications in which controlled breakdown is intended, such as plasma medicine, lightning protection, and materials processing. The electron emission from the cathode is a critical source of electrons which then undergo impact ionization to produce electrical breakdown. In this study, the role of secondary electron yields due to photons (γ ph) and ions (γ i) in direct current breakdown is investigated using a particle-in-cell direct simulation Monte Carlo model. The plasma studied is a one-dimensional discharge in 50 Torr of pure helium with a platinum cathode, gap size of 1.15 cm, and voltages of 1.2-1.8 kV. The current traces are compared with experimental measurements. Larger values of γ ph generally result in a faster breakdown, while larger values of γ i result in a larger maximum current. The 58.4 nm photons emitted from He(21P) are the primary source of electrons at the cathode before the cathode fall is developed. Of the values of γ ph and γ i investigated, those which provide the best agreement with the experimental current measurements are γ ph = 0.005 and γ i = 0.01. These values are significantly lower than those in the literature for pristine platinum or for a graphitic carbon film which we speculate may cover the platinum. This difference is in part due to the limitations of a one-dimensional model but may also indicate surface conditions and exposure to a plasma can have a significant effect on the secondary electron yields. The effects of applied voltage and the current produced by a UV diode which was used to initiate the discharge, are also discussed.
The stability of low-index platinum surfaces and their electronic properties is investigated with density functional theory, toward the goal of understanding the surface structure and electron emission, and identifying precursors to electrical breakdown, on nonideal platinum surfaces. Propensity for electron emission can be related to a local work function, which, in turn, is intimately dependent on the local surface structure. The (1×N) missing row reconstruction of the Pt(110) surface is systematically examined. The (1×3) missing row reconstruction is found to be the lowest in energy, with the (1×2) and (1×4) slightly less stable. In the limit of large (1×N) with wider (111) nanoterraces, the energy accurately approaches the asymptotic limit of the infinite Pt(111) surface. This suggests a local energetic stability of narrow (111) nanoterraces on free Pt surfaces that could be a common structural feature in the complex surface morphologies, leading to work functions consistent with those on thermally grown Pt substrates.
This paper presents a new method for modeling rarefied gas flows based on hybridization of direct simulation Monte Carlo (DSMC) and discrete velocity method (DVM)-based quasi-particle representations of the velocity distribution function. It is aimed at improving the resolution of the tails of the distribution function (compared with DSMC) and computational efficiency (compared with DVM). Details of the method, such as the collision algorithm and the particle merging scheme, are discussed. The hybrid approach is applied to the study of noise in a Maxwellian distribution, computation of electron-impact ionization rate coefficient, as well as numerical simulation of a supersonic Couette flow. The hybrid-based solver is compared with pure DSMC and DVM approaches in terms of accuracy, computational speed, and memory use. It is shown that such a hybrid approach can provide a lower computational cost than a pure DVM approach, while being able to retain accuracy in modeling high-velocity tails of the distribution function. For problems where trace species have a significant impact on the flow physics, the proposed method is shown to be capable of providing better computational efficiency and accuracy compared with standard fixed-weight DSMC.
For high voltage electrical devices, prevention of high voltage breakdown is critical for device function. Use of polymeric encapsulation such as epoxies is common, but these may include air bubbles or other voids of varying size. The present work aimed to model and experimentally determine the size dependence of breakdown voltage for voids in an epoxy matrix, as a step toward establishing size criteria for void screening. Effects were investigated experimentally for both one-dimensional metal/epoxy/air/epoxy/metal gap sizes from 50 μm to 10 mm, as well as spherical voids of 250 μm, 500 μm, 1 mm and 2 mm sizes. These experimental results were compared to modified Paschen curve and particle-in-cell discharge models; minimum breakdown voltages of 6 - 8.5 kV appeared to be predicted by 1D models and experiments, with minimum breakdown voltage for void sizes of 0.2 - 1 mm. In a limited set of 3D experiments on 250 μm, 500 μm, 1 mm and 2 mm voids within epoxy, the minimum breakdown voltages observed were 18.5 - 20 kV, for 500 μm void sizes. These experiments and models are aimed at providing initial size and voltage criteria for tolerable void sizes and expected discharge voltages to support design of encapsulated high voltage components.
The influence of different quantum yields for photons and secondary emission yields for ions striking a surface is investigated. Using a one-dimensional particle-in-cell simulation, these secondary emission coefficients are varied to observe the impact on discharge current. The discharge is assumed to occur in pure helium gas at a pressure of 75 torr. To handle binary particle interactions, the Direct Simulation Monte Carlo (DSMC) method is utilized. The model includes electron-neutral interactions, neutral-neutral interactions, and photon-neutral interactions. It is observed that the discharge current in the early stages of discharge is heavily dependent upon the quantum yield due to photon impact. In the later stages of discharge, the current depends on both the quantum yield and secondary emission coefficient for ion impact.
In most models of vacuum breakdown, there is some initial emission of electrons from the cathodic surface, usually employing some form of Fowler-Nordheim emission. While this may be correct for 'textbook' surfaces, it is generally unreliable for real surfaces and fitted parameters are often used. For example, the beta employed is generally unphysical based on usual definitions (e.g., it incorporates more, but unexplained, physics than just a geometry-based field concentration effect). In this work, we describe experimental efforts to better characterize which surface structure parameters influence the vacuum field emission current.