A highly parallelizable fluid plasma simulation tool based upon the first-order drift-diffusion equations is discussed. Atmospheric pressure plasmas have densities and gradients that require small element sizes in order to accurately simulate the plasm resulting in computational meshes on the order of millions to tens of millions of elements for realistic size plasma reactors. To enable simulations of this nature, parallel computing is required and must be optimized for the particular problem. Here, a finite-volume, electrostatic drift-diffusion implementation for low-temperature plasma is discussed. The implementation is built upon the Message Passing Interface (MPI) library in C++ using Object Oriented Programming. The underlying numerical method is outlined in detail and benchmarked against simple streamer formation from other streamer codes. Electron densities, electric field, and propagation speeds are compared with the reference case and show good agreement. Convergence studies are also performed showing a minimal space step of approximately 4 μm required to reduce relative error to below 1% during early streamer simulation times and even finer space steps are required for longer times. Additionally, strong and weak scaling of the implementation are studied and demonstrate the excellent performance behavior of the implementation up to 100 million elements on 1024 processors. Lastly, different advection schemes are compared for the simple streamer problem to analyze the influence of numerical diffusion on the resulting quantities of interest.
Molecular dynamics simulations are used to test when the particle-in-cell (PIC) method applies to atmospheric pressure plasmas. It is found that PIC applies only when the plasma density and macroparticle weight are sufficiently small because of two effects associated with correlation heating. The first is the physical effect of disorder-induced heating (DIH). This occurs if the plasma density is large enough that a species (typically ions) is strongly correlated in the sense that the Coulomb coupling parameter exceeds one. In this situation, DIH causes ions to rapidly heat following ionization. PIC is not well suited to capture DIH because doing so requires using a macroparticle weight of one and a grid that well resolves the physical interparticle spacing. These criteria render PIC intractable for macroscale domains. The second effect is a numerical error due to Artificial Correlation Heating (ACH). ACH is like DIH in that it is caused by the Coulomb repulsion between particles, but differs in that it is a numerical effect caused by a macroparticle weight larger than one. Like DIH, it is associated with strong correlations. However, here the macroparticle coupling strength is found to scale as Γ w2/3, where Γ is the physical coupling strength and w is the macroparticle weight. So even if the physical coupling strength of a species is small, as is expected for electrons in atmospheric pressure plasmas, a sufficiently large macroparticle weight can cause the macroparticles to be strongly coupled and therefore heat due to ACH. Furthermore, it is shown that simulations in reduced dimensions exacerbate these issues.
Recent findings suggest that ions are strongly correlated in atmospheric pressure plasmas if the ionization fraction is sufficiently high ( ≳ 10 − 5 ). A consequence is that ionization causes disorder-induced heating (DIH), which triggers a significant rise in ion temperature on a picosecond timescale. This is followed by a rise in the neutral gas temperature on a longer timescale of up to nanoseconds due to ion-neutral temperature relaxation. The sequence of DIH and ion-neutral temperature relaxation suggests a new mechanism for ultrafast neutral gas heating. Previous work considered only the case of an instantaneous ionization pulse, whereas the ionization pulse extends over nanoseconds in many experiments. Here, molecular dynamics simulations are used to analyze the evolution of ion and neutral gas temperatures for a gradual ionization over several nanoseconds. The results are compared with published experimental results from a nanosecond pulsed discharge, showing good agreement with a measurement of fast neutral gas heating.
Ion diffusion in atmospheric pressure plasmas is examined and particular attention is paid to the fact that ion-ion interactions can be influenced by strong Coulomb coupling. Three regimes are identified. At low ionization fractions ( x i ≲ 10 − 6 ), standard weakly correlated ion-neutral interactions set the diffusion rate. At moderate ionization fractions ( 10 − 6 ≲ x i ≲ 10 − 2 ) there is a transition from ion-neutral to ion-ion collisions setting the diffusion rate. In this regime, the effect of strong Coulomb coupling in ion-ion collisions is accounted for by applying the mean force kinetic theory. Since both ion-neutral and ion-ion interactions contribute a comparable amount to the total diffusion rate, models (such as particle-in-cell or fluid) must account for both contributions. At high ionization fractions ( x i ≳ 10 − 2 ), strongly correlated ion-ion collisions dominate and the plasma is heated substantially by a disorder-induced heating (DIH) process associated with strong correlations. The temperature increase due to DIH strongly influences the ion diffusion rate. This effect becomes even more important, and occurs at lower ionization fractions, as the pressure increases above atmospheric pressure. In addition to ion diffusion, DIH affects the neutral gas temperature, therefore influencing the neutral diffusion rate. Model predictions are tested using molecular dynamics simulations, which included a Monte Carlo collision routine to simulate the effect of ion-neutral collisions at the lowest ionization fractions. The model and simulations show good agreement over a broad range of ionization fractions. The results provide a model for ion diffusion, on a wide range of ionization fractions and pressures, solely considering the elastic contribution to the diffusion coefficient—as an illustration of how strong Coulomb coupling influences diffusion processes in general.
Large, pulsed power and high voltage systems often employ a stack of insulators to separate a vacuum section away from water or oil sections. The size of this insulator stack often drives overall costs and feasibility of these systems. An electric breakdown along the insulator surface is a primary failure mechanism and is especially impactful if it occurs while power is still being delivered downstream. This report describes a set of experimental and modeling investigations into the cause of these breakdowns, especially focusing on the much less well-understood anode-initiated breakdowns that occur during early parts of power delivery. Additionally, new diagnostics for assessing relevant material properties and behavior of insulators are described. These results describe breakdown behavior and evolution at new temporal and spatial fidelities and provide hypotheses and some answers as to how these breakdowns can occur. This new understanding of the roles of different physics phenomena guide modifications and trade-offs in generating newer insulator stack designs that are smaller and/or have higher electrical stress thresholds.
It is necessary to establish confidence in high-consequence codes containing an extensive suite of physics algorithms in the regimes of interest. Verification problems allow code developers to assess numerical accuracy and increase confidence that specific sets of model physics were implemented correctly in the code. The two main verification techniques are code verification and solution verification. In this work, we present verification problems that can be used in other codes to increase confidence in simulations of relativistic beam transport. Specifically, we use the general plasma code EMPIRE to model and compare with the analytical solution to the evolution of the outer radial envelope of a relativistic charged particle beam. We also outline a benchmark test of a relativistic beam propagating through a vacuum and pressurized gas cell, and present the results between EMPIRE and the hybrid code GAZEL. Further, we discuss the subtle errors that were caught with these problems and detail lessons learned.
This presentation describes a new effort to better understand insulator flashover in high current, high voltage pulsed power systems. Both experimental and modeling investigations are described. Particular emphasis is put upon understand flashover that initiate in the anode triple junction (anode-vacuum-dielectric).
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.
We propose a new scheme for simulation of collisions with multiple possible outcomes in variable-weight DSMC computations. The scheme is applied to a 0-D ionization rate coefficient computation, and 1-D electrical breakdown simulation. We show that the scheme offers a significant (up to an order of magnitude) improvement in the level of stochastic noise over the usual acceptance-rejection algorithm, even when controlling for the slight additional computational costs. Furthermore, the benefits and performance of the scheme are analyzed in detail, and possible extensions are proposed.
We introduce a robust verification tool for computational codes, which we call Stochastic Robust Extrapolation based Error Quantification (StREEQ). Unlike the prevalent Grid Convergence Index (GCI) [1] method, our approach is suitable for both stochastic and deterministic computational codes and is generalizable to any number of discretization variables. Building on ideas introduced in the Robust Verification [2] approach, we estimate the converged solution and orders of convergence with uncertainty using multiple fits of a discretization error model. In contrast to Robust Verification, we perform these fits to many bootstrap samples yielding a larger set of predictions with smoother statistics. Here, bootstrap resampling is performed on the lack-of-fit errors for deterministic code responses, and directly on the noisy data set for stochastic responses. This approach lends a degree of robustness to the overall results, capable of yielding precise verification results for sufficiently resolved data sets, and appropriately expanding the uncertainty when the data set does not support a precise result. For stochastic responses, a credibility assessment is also performed to give the analyst an indication of the trustworthiness of the results. This approach is suitable for both code and solution verification, and is particularly useful for solution verification of high-consequence simulations.
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 describes the verification and validation (V&V) framework developed for the stochastic Particle-in-Cell, Direct Simulation Monte Carlo code Aleph. An ideal framework for V&V from the viewpoint of the authors is described where a physics problem is defined, and relevant physics models and parameters to the defined problem are assessed and captured in a Phenomena Identification and Ranking Table (PIRT). Furthermore, numerous V&V examples guided by the PIRT for a simple gas discharge are shown to demonstrate the V&V process applied to a real-world simulation tool with the overall goal to demonstrably increase the confidence in the results for the simulation tool and its predictive capability. Although many examples are provided here to demonstrate elements of the framework, the primary goal of this work is to introduce this framework and not to provide a fully complete implementation, which would be a much longer document. Comparisons and contrasts are made to more usual approaches to V&V, and techniques new to the low-temperature plasma community are introduced. Specific challenges relating to the sufficiency of available data (e.g., cross sections), the limits of ad hoc validation approaches, the additional difficulty of utilizing a stochastic simulation tool, and the extreme cost of formal validation are discussed.
Helium is frequently used as a working medium for the generation of plasmas and is capable of energetic photon emissions. These energetic photon emissions are often attributed to the formation of helium excimer and subsequent photon emission. When the plasma device is exposed to another gas, such as nitrogen, this energetic photon emission can cause photoionization and further ionization wave penetration into the additional gas. Often ignored are the helium resonance emissions that are assumed to be radiation trapped and therefore not pertinent to photoionization. Here, experimental evidence for the presence of helium atomic emission in a pulsed discharge at ten's of Torr is shown. Simulations of a discharge in similar conditions agree with the experimental measurements. In this context, the role of atomic and molecular helium light emission on photoionization of molecular nitrogen in an ionization wave is studied using a kinetic modeling approach that accounts for radiation dynamics in a developing low-temperature plasma. Three different mixtures of helium at a total pressure of 250 Torr are studied in simulation. Photoionization of the nitrogen molecule by vacuum ultraviolet helium emission is used as the only seed source ahead of the ionization front. It is found that even though radiation trapped, the atomic helium emission lines are the significant source of photoionization of nitrogen. The significant effect of radiation trapped photon emission on ionization wave dynamics demonstrates the need to consider these radiation dynamics in plasma reactors where self-absorbed radiation is ignored.
A novel method based on combining Direct Simulation Monte Carlo (DSMC) and Discrete Velocity Method (DVM) representations of the velocity distribution functions in velocity space is applied to rarefied ionized gas flows in order to study its efficiency and accuracy. The objective is to improve the efficiency of modeling of flows where trace populations have a significant effect on the flow physics. Numerical results are obtained for a 0-dimensional flow of a Ar/Ar+ /e− mixture and compared with the BOLSIG+ solver.
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.
In the present research, a new method for simulation of rarefied gas flows is proposed, a velocity-space hybrid of both a DSMC representation of particles and a discrete velocity quasi-particle representation of the distribution function. The hybridization scheme is discussed in detail, and is numerically verified for two test-cases: the BKW relaxation problem and a stationary Maxwellian distribution. It is demonstrated that such a velocity-space hybridization can provide computational benefits when compared to a pure discrete velocity method or pure DSMC approach, while retaining some of the more attractive properties of discrete velocity methods. Further possible improvements to the velocity-space hybrid approach are discussed.
3D Particle-In-Cell Direct Simulation Monte Carlo (PIC-DSMC) simulations of cm-sized devices cannot resolve atomic-scale (nm) surface features and thus one must generate micron-scale models for an effective “local” work function, field enhancement factor, and emission area. Here we report on development of a stochastic effective model based on atomic-scale characterization of as-built electrode surfaces. Representative probability density distributions of the work function and geometric field enhancement factor (beta) for a sputter-deposited Pt surface are generated from atomic-scale surface characterization using Scanning Tunneling Microscopy (STM), Atomic Force Microscopy (AFM), and Photoemission Electron Microscopy (PEEM). In the micron-scale model every simulated PIC-DSMC surface element draws work functions and betas for many independent “atomic emitters”. During the simulation the field emitted current from an element is computed by summing each “atomic emitter's” current. This model has reasonable agreement with measured micron-scale emitted currents across a range of electric field values.
The purpose of this paper is to characterize the need for improved predictive capabilities in low-temperature plasma (LTP) science, and to identify possible means of accomplishing this. While these means may constitute an initiative of their own, we consider these ideas to have widespread importance to discovery plasma science. Therefore, it is our hope that these ideas are more generally incorporated in future work.
We report on the verification of elastic collisions in EMPIRE-PIC and EMPIRE-Fluid in support of the ATDM L2 V&V Milestone. The thermalization verification problem and the theory behind it is presented along with an analytic solution for the temperature of each species over time. The problem is run with both codes under multiple parameter regimes. The temperature over time is compared between the two codes and the theoretical results. A preliminary convergence analysis is performed on the results from EMPIRE-PIC and EMPIRE-Fluid showing the rate at which the codes converge to the analytic solution in time (EMPIRE-Fluid) and particles (EMPIRE-PIC).
Modern computational validation efforts rely on comparison of known experimental quantities such as current, voltage, particle densities, and other plasma properties with the same values determined through simulation. A discrete photon approach for radiation transport was recently incorporated into a particle-in-cell/direct simulation Monte Carlo code. As a result, spatially and temporally resolved synthetic spectra may be generated even for non-equilibrium plasmas. The generation of this synthetic spectra lends itself to potentially new validation opportunities. In this work, initial comparisons of synthetic spectra are made with experimentally gathered optical emission spectroscopy. A custom test apparatus was constructed that contains a 0.5 cm gap distance parallel plane discharge in ultra high purity helium gas (99.9999%) at a pressure of 75 Torr. Plasma generation is initiated with the application of a fast rise-time, 100 ns full-width half maximum, 2.0 kV voltage pulse. Transient electrical diagnostics are captured along with time-resolved emission spectra. A one-dimensional simulation is run under the same conditions and compared against the experiment to determine if sufficient physics are included to model the discharge. To sync the current measurements from experiment and simulation, significant effort was undertaken to understand the kinetic scheme required to reproduce the observed features. Additionally, the role of the helium molecule excimer emission and atomic helium resonance emission on photocurrent from the cathode are studied to understand which effect dominates photo-feedback processes. Results indicate that during discharge development, atomic helium resonance emission dominates the photo-flux at the cathode even though it is strongly self-absorbed. A comparison between the experiment and simulation demonstrates that the simulation reproduces observed features in the experimental discharge current waveform. Furthermore, the synthesized spectra from the kinetic method produces more favorable agreement with the experimental data than a simple local thermodynamic equilibrium calculation and is a first step towards using spectra generated from a kinetic method in validation procedures. The results of this study produced a detailed compilation of important helium plasma chemistry reactions for simulating transient helium plasma discharges.
A fully resolved kinetic model (particle-in-cell and direct simulation Monte Carlo for particle/photon collisions) of a near atmospheric pressure ionization wave is presented here. Fully resolving the required numerical spatial (sub-μm) and temporal scales (tens of fs) for atmospheric pressure discharges in three-dimensions is still a challenging task on modern super computers. To keep the overall problem tractable, the total number of elements are reduced by only simulating a 10° wedge rather than a full 360° geometry. The ionization wave is generated in a needle-plane configuration with a gap size of 250 μm and a background of nitrogen and helium gas. A voltage of 1500 V is applied to the anode and an initial electron and ion density of 109 cm-3 is seeded in a region near the anode electrode tip and extending towards the cathode. As these initial electrons are swept away, photoionization and photoemission create new electrons and allow the ionization front to propagate towards the cathode. Results from the 90% N2, 10% He discharge indicate that photoionization has minimal impact on plasma formation processes and cathode photoemission is the dominant mechanism for new electrons. In the 90% He, 10% N2 discharge case, however, photoionization likely has an impact as the observed locations of photoionization occur far enough away from the ionization front to allow for sufficient avalanche processes that contribute to the propagation of the ionization wave. Additionally, the electron energy distribution functions in the 90% He, 10% N2 case indicate that there is less energy loss to the low lying molecular N2 electronic states as well as the vibrational and rotational modes. This leads to higher electron energies and faster plasma development times of ∼0.4 ns for the 90% He, 10% N2 case, and ∼1.5 ns for the 90% N2, 10% He case. In addition to analysis of the ionization wave results, the overall challenges associated with simulations near atmospheric pressure discharges in three-dimensions are discussed, including the limitations of the 10° wedge that produces, at least qualitatively, minimal 3D effects.
This report describes a model for the time development of carrier distributions within a metallic or semiconductor target after the onset of an incident laser pulse. The dynamics of electron and hole populations in momentum-resolved conduction- and valence-band states are treated at the level of carrier-carrier and carrier-phonon scattering. These scattering events result in plasma and lattice heating, which in turn lead to electron thermionic emission and tunneling, and target material ablation. A fairly phenomenological approach is taken to mitigate numerical computation demands, in order to facilitate parametric studies. Two examples of application are presented. One involve s the incident of an intense near-infrared laser pulse on a solid aluminum target, where the goal is to connect excited species emission to physics at a band-structure level. The second involves modeling the trigger mechanism in laser-triggered high-voltage switches, where the results are used as input to highly intensive particle-in-cell (PIC) plasma simulations of switch operation.
A kinetic description for electronic excitation of helium for principal quantum number n 4 has been included into a particle-in-cell (PIC) simulation utilizing direct simulation Monte Carlo (DSMC) for electron-neutral interactions. The excited electronic levels radiate state-dependent photons with wavelengths from the extreme ultraviolet (EUV) to visible regimes. Photon wavelengths are chosen according to a Voigt distribution accounting for the natural, pressure, and Doppler broadened linewidths. This method allows for reconstruction of the emission spectrum for a non-thermalized electron energy distribution function (EEDF) and investigation of high energy photon effects on surfaces, specifically photoemission. A parallel plate discharge with a fixed field (i.e. space charge neglected) is used to investigate the effects of including photoemission for a Townsend discharge. When operating at a voltage near the self-sustaining discharge threshold, it is observed that the electron current into the anode is higher when including photoemission from the cathode than without even when accounting for self-absorption from ground state atoms. The photocurrent has been observed to account for as much as 20% of the total current from the cathode under steady-state conditions.
The goal of this report is to document the current status of Aleph with regards to electron collisions under an electric field. Aleph and the community-accepted BOLSIG+ code are both used to compute reactions rates for a set of 25 electron-nitrogen interactions. A reasonable comparison is found (see below) providing evidence that Aleph is successfully simulating or implementing: (1) Particle-particle collision cross-sections via DSMC methodology, (2) Energy balance for simple particle interactions, and (3) Electron energy distribution function (EEDF) evolution