Simulations of atmospheric pressure streamer discharges: Mechanisms of streamer propagation and related physics
Explore the source record for details and available documents.
Nikic, Dejan [Sandia National Laboratories (SNL-NM↗
SEARCH · Engineering Papers
Results for “streamer discharge”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Numerical-heating effects in atmospheric pressure streamer discharges simulated with a PIC code
Artificial heating in plasma simulations is a well-known phenomenon which occurs when, among other things, the Debye length is poorly resolved by the simulation mesh. Here, in this work, the degree to which numerical-heating occurs during a simulation of a nanosecond atmospheric pressure streamer discharge is examined. The streamer is simulated using a two-dimensional finite-element, particle-in-cell code Empire, which uses direct simulation Monte Carlo for binary particle interactions. Initially, an estimate of the numerical-heating rate applied to Empire is performed using a simple plasma model. Second, a positive atmospheric pressure streamer discharge simulation is performed to study the effects of numerical heating on plasma density, electron temperature, and streamer velocity. The nominal Debye length is approximately 1 μm and the amount of numerical heating introduced in the simulation is varied by using mesh sizes ranging from 2 μm to 20 μm. A measurable numerical heating quantity is proposed that can be used to estimate the appropriate element size and quantify the numerical-heating that can be expected over the simulation time for an atmospheric pressure streamer. In conclusion while Δx/λ D violations can be an issue it is not likely to be an issue with streamer discharges that are temporally short and occur in environments where collision frequencies are high. This result validates the rationale of grid size choices for a large amount of previously published works where Δx/λ D violation was not clearly addressed. Primary finding of this work is that numerical heating is of minor concern for plasma simulations where electron–neutral collisions are numerous such that multiple collisions can occur within a single plasma period.
Nikic, Dejan [University of New Mexico, Albuquerqu↗
Influence of photoionization on propagation of streamer discharges at atmospheric pressure
Explore the source record for details and available documents.
Nikic, Dejan [UNM]↗
Breakdown front dynamics of streamer-like discharge in electron-charged polymethyl methacrylate
We present the results of an experimental study of fast breakdown in electron-charged polymethyl methacrylate. We irradiate bulk polymethyl methacrylate disks with diameters up to one meter at different implanted charge densities and measure the discharge current during the forced electrical breakdown of the material. We infer the breakdown dynamics from these current waveforms, including the velocity and time dependence of the electric field driving the breakdown, and compare these results with the physical electrical tree patterns left behind in the material. We find that the dynamics and physical characteristics of the breakdown channels in electron-irradiated solids depart from typical expectations of electrical treeing behavior in solid materials. Thus, we interpret these results as an expression of streamer discharges in dense gases motivated by the existence of trapped gases in the solid due to radiation damage. We show that the dynamics of the breakdown channels in the electron-charged solid dielectric material is well described by applying standard streamer physics to this physical system. The results show that the dynamics of breakdown channels in solid dielectric material is a promising avenue to further understand streamer discharges in different media under extreme conditions.
36 MATERIALS SCIENCE↗
Enhancing photoionization rate calculations in low-temperature plasmas using spectral methods
Photoionization plays a central role in the development of streamer discharges and other non-equilibrium plasma phenomena. It creates seed electrons, which are essential for positive streamer propagation, allowing the ionization front to move forward. Because of this, accurate modeling of photoionization is very important for predicting streamer behavior and plasma evolution. The photoionization process in air (N 2 – O 2 mixture) is often described by the Zheleznyak model (1982). This model is usually solved through Helmholtz-type equations that approximate the Zheleznyak photoionization model (Zheleznyak et al. 1982) as Partial Differential Equations (PDEs). Conventional numerical methods, such as the Finite Difference Method (FDM) or Finite Volume Method (FVM), are widely used to solve these equations. Although they are prevalent, the computational cost of these methods due to their need for matrix operations and iterative solver is demanding. To address this challenge, this work develops a spectral solver based on the Fast Fourier Transform (FFT) combined with Discrete Cosine Transform (DCT) and Discrete Sine Transform (DST) to calculate the photoionization rate efficiently in an axisymmetric cylindrical domain. This method naturally satisfies the boundary conditions used in the model and converts the PDE into algebraic ones in spectral space. Thus, avoids the need for iterative matrix solvers. When compared with FDM results, it is demonstrated that the new solver not only maintains accuracy, but also reduces the computational cost, showing a performance increase of approximately 100 compared to FDM over a wide range of problem sizes. The method is parallelized using Message Passing Interface (MPI) and has been integrated into a fluid plasma model for streamer simulation. Here, this FFT-based approach provides a fast and reliable alternative for calculating photoionization in fluid models, helping large-scale plasma simulations run faster and efficiently, and allows higher-resolution simulation without extra computational cost.
Axisymmetric system↗
Experimental study of high-energy single-pulse nanosecond discharges in pin-to-pin configuration
Single-pulse pin-to-pin nanosecond pulsed discharges in ambient air (gap distances of 3–7 mm) were investigated at high pulse energies (∼20–30 mJ deposited energy per 11 ns full width at half maximum pulse). A Michelson interferometry setup (532 nm continuous-wave laser) was employed to record time-resolved interferograms of the discharge, enabling spatially resolved calculation of the electron number density. An intensified charge-coupled device camera was used to capture the spatiotemporal evolution of the discharge (streamer formation and spark channel development), and Coherent Anti-Stokes Raman Spectroscopy (CARS) was used to measure post-discharge N 2 vibrational temperatures. Discharge current and voltage were monitored with a back-current shunt. The discharge was initiated with simultaneous cathode-directed and anode-directed streamers that bridge the gap within ∼1 ns and form a luminous, filamentary plasma channel immediately after breakdown. Laser interferometry measurements showed peak electron number densities of the order of 10 1 7–10 18 cm −3 , occurring about 15–20 ns after pulse arrival at the discharge gap. Shorter gap discharges yielded higher peak electron densities, consistent with the higher energy density in the smaller gaps. Spatially, the electron density was highest near the electrodes and decreased toward the midgap region, with 5 and 7 mm gaps exhibiting a pronounced drop in the central channel. CARS indicated initial vibrational temperatures of approximately 4000–6000 K in the spark core of ∼50 ns after the discharge onset, decaying on a ∼500 ns timescale as the plasma cooled and recombined.
Coherent anti-Stokes Raman spectroscopy↗
Intense cavitation-assisted electric discharge as a promising tool for water treatment
This study investigates interrelations between one-electrode Cavitation-Assisted Electric Discharge (CAED), two-electrode CAED, and recently discovered Intense CAED (I-CAED). The one-electrode CAED is a self-triggered nanosecond discharge with pulse energy in the micro-Joule range, which can be generated even by a DC voltage. I-CAED consists of a non-equilibrium part within a low-pressure cavitating region and a micro-spark traversing a liquid film. We hypothesize that CAED propagates from the high-voltage electrode as an ionization wave through bubbles of saturated vapor. Subsequently, the streamer-like discharges in the bubbles may form a continuous plasma channel. Inside the cavitating region, the plasma is strongly non-equilibrium, providing an ideal environment for generating chemically unstable species such as hydrogen peroxide (H 2 O 2 ). Plasma of I-CAED spark is characterized by high electron density and near-thermal equilibrium, emitting a continuous ultraviolet spectrum. The combination of these different discharge parts makes I-CAED in water a highly effective tool for the Advanced Oxidation Process, particularly in water disinfection. Experimentally demonstrated Electric Energy per Order value for disinfection of E. coli-contaminated water is as low as 0.135 ± 0.035 kWh/m 3 /order. Estimates show that the implementation of “dry electrodes” configuration reduces the erosion rate of the electrode material by at least one order of magnitude. Spectral analysis reveals that the continuum emission generated by I-CAED in proximity to metal electrodes deviates from the spectra of discharges spatially decoupled from the electrodes. We assume that this spectral divergence is attributable to blackbody-like emission originating from metallic nanoparticles form during the electrode's erosion process.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
The Nature of Standing Plasma Waves Following Fast Breakdown in Thunderstorms
Narrow bipolar events (NBEs) are produced by large-scale systems of streamer corona discharges that occur inside thunderstorms and radiate intensely at very-high frequencies. Here we examine 14 high-power NBEs recorded in Florida which are followed by a lasting oscillatory signal. We isolate the dominant frequencies of the oscillations using empirical mode decomposition methods, and cross-reference this information with the radiator size inferred from broadband interferometry. We conclude that the oscillations are evidence of a standing electromagnetic wave in the plasma cavity created by the preceding streamer system, that is, the cavity's fundamental mode with wavelength equal to twice the radiator length. Contrasting the dispersive and attenuation properties of these waves across several different events allows us to conclude that the cavity is filled with ions resulting in an average conductivity of 0.01–0.1 mS/m, which are left behind after the streamers decay via electron attachment to oxygen molecules.
58 GEOSCIENCES↗
Hydrogen from low-density polyethylene via nonthermal plasma: Effects of energy density and process parameters
Nonthermal plasma processes are promising for the modular valorization of plastic waste, especially into hydrogen and carbon materials, due to their high intensity, lack of reliance on catalysts or consumables, and suitability to be directly powered by electricity. We investigate the production of hydrogen from low-density polyethylene (LDPE) as a plastic waste model using streamer Dielectric Barrier Discharge (sDBD) plasma in nitrogen at atmospheric pressure. Here, we examine the effects of process energy density (energy input per unit of feedstock mass), feedstock mass, and plasma intensity (electric voltage) on plasma properties, hydrogen yield and energy efficiency via gas chromatography, optical emission spectroscopy, and electrical diagnostics, together with reactor-scale and nonlinear electric circuit modeling. The characteristic temperature of free electrons in the sDBD plasma is approximately 15000 K (1.3 eV), and that of gas species 10 times lower, demonstrating strong thermal non-equilibrium that can lead to molecular bond scission via charged species impact rather than direct heating. Experimental results show that higher energy density leads to greater hydrogen production and diminishing energy efficiency, and that higher plasma intensity and larger feedstock mass lead to greater hydrogen yield due to higher plasma temperatures and enhanced energy fluxes to the feedstock.
08 HYDROGEN↗
High-rate heavy-ion tracking using MWPCs and PPACs with an Anti-Discharge Unit
This work reports on the development of two robust, heavy-ion beam tracking concepts operating at low pressure (< 15 Torr) for high rate applications (> 200 kHz). The first concept consists of a Multi-Wire Proportional Counter (MWPC) with a central anode consisting of 12 μm Au-plated Tungsten wires spaced 1 mm from each other. The anode grid is sandwiched between two segmented cathodes aligned orthogonally in the two dimensions for (x,y) particle localization. The second detector is a Parallel-Plate Avalanche Counter (PPAC) that uses the same readout geometry as the MWPC, but replaces the wire anode with a 150 nm silver layer deposited on both sides of a thin (< 1 mg/cm 2 ) polypropylene foil. Additionally, the bias circuitry for the PPAC anode central foil is equipped with an Anti-Discharge Unit (ADU) to prevent transitions from proportional operation to streamer formation, thereby avoiding damaging discharges. The localization capability of both detectors was tested with a low-rate alpha-particle source (241-Am). A position resolution of < 1 mm (FWHM) was achieved under stable, high-gas-gain (> 1000) operating conditions. Their performance in terms of detection efficiency as a function of the isotope charge (Z) was determined by irradiating the detectors with a cocktail beam (Z ≤ 15) with energy of ∼ 100 MeV/u. Full detection efficiency is maintained for all available fragments under optimal operational conditions (i.e., voltage bias). Full detection efficiency was achieved at rates above 200 kHz by irradiating the detectors with a 1 cm diameter 238 U beam at an energy of 143 MeV/u.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Vidyut3d: A GPU accelerated fluid solver for non-equilibrium plasmas on adaptive grids
We present the numerical methods, programming methodology, verification, and performance assessment of a non-equilibrium plasma fluid solver that can effectively utilize current and upcoming central processing and graphics processing unit (CPU+GPU) architectures, in this work. Our plasma fluid model solves the coupled conservation equations for species transport, electrostatic Poisson and electron temperature on adaptive Cartesian grids. Our solver is written using performance portable adaptive-grid/particle management library, AMReX, and is portable over widely available vendor specific GPU architectures. We present verification of our solver using method of manufactured solutions that indicate formal second order accuracy with central diffusion and fifth-order weighted-essentially-non-oscillatory (WENO) advection scheme. We also verify our solver with published literature on capacitive discharges and atmospheric pressure streamer propagation. We demonstrate the use of our solver on two 3D simulation cases: an atmospheric streamer propagation in Ar-H2 mixtures and a low pressure three-electrode radio frequency reactor. Our performance studies on three different CPU+GPU architectures indicate ~ 150-400X speed-up using AMD and NVIDIA GPUs per time step compared to a single CPU core for a 4 million cell simulation with 15 species.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Vidyut3d: A Gpu Accelerated Fluid Solver for Non-Equilibrium Plasmas on Adaptive Grids
We present the numerical methods, programming methodology, verification, and performance assessment of a non-equilibrium plasma fluid solver that can effectively utilize current and upcoming central processing and graphics processing unit (CPU+GPU) architectures, in this work. Our plasma fluid model solves the coupled conservation equations for species transport, electrostatic Poisson and electron temperature on adaptive Cartesian grids. Our solver is written using performance portable adaptive-grid/particle management library, AMReX, and is portable over widely available vendor specific GPU architectures. We present verification of our solver using method of manufactured solutions that indicate formal second order accuracy with central diffusion and fifth-order weighted-essentially-non-oscillatory (WENO) advection scheme. We also verify our solver with published literature on capacitive discharges and atmospheric pressure streamer propagation. We demonstrate the use of our solver on two 3D simulation cases: an atmospheric streamer propagation in Ar-H2 mixtures and a low pressure twin electrode radio frequency reactor. Our performance studies on three different CPU+GPU architectures indicate approximately 150-400X speed-up using AMD and NVIDIA GPUs per time step compared to a single CPU core for a 4 million cell simulation with 15 species.
Sitaraman, Hariswaran↗
Dynamics of the water-plasma interface in various discharge modes of atmospheric-pressure plasmas
This paper presents an experimental study of the dynamics of the water-plasma interface during the interaction high-voltage (HV) atmospheric pressure discharge with the water surface. Under the influence of the voltage applied to the pin type electrode, discharges are formed at the air layer between the sharp-tip HV electrode and deionized water. Three discharge regimes were identified from synchronized electrical waveforms and imaging: (i) a weak linear regime at 3–10.6 kV, (ii) a branching streamer regime at the maximum applied voltage of 12.6 kV, and (iii) a continuous arc-like regime occurring during the transition to lower voltage and higher current (e.g., U ≈ 3.2 kV, I ≈ 4.07 mA). High-speed shadowgraph images showed a symmetric interfacial cavity whose depth increased nonlinearly with voltage from h 0 ≈ 0.1 mm at 3 kV to h 0 ≈ 2.7 mm at 10.6 kV, reaching a maximum depth of h 0 ≈ 5.9 mm at 12.6 kV, while the cavity disappeared in the continuous-channel regime and was replaced by outward-propagating wave-like motion. The deformation of the water surface during the discharge-water interaction is governed by the balance between electric field forces, surface tension, gravitational forces, and electrohydrodynamic forces, whose relative contributions vary with the applied voltage. This effect was explained by quantitatively estimating the magnitudes of the acting forces and establishing the force balance. Optical emission spectroscopy in the continuous-channel regime indicated air plasma signatures and yielded a gas temperature of approximately 400 K, while prolonged operation (≈ 20 min) increased the water temperature to ~ 70 °C, reduced the water-layer thickness from 6 mm to 3 mm, and decreased pH from 7 to 4. These regime-resolved, quantitative results clarify how the dominant interfacial forcing shifts with discharge mode and provide a mechanistic basis for controlling plasma-water interactions in atmospheric-pressure applications.
atmospheric pressure plasma↗
Breakdown and Discharge Characteristics of High Repetition Frequency Nanosecond Pulsed Air Dielectric Barrier Discharge
This study explores the breakdown characteristics and discharge modes in a parallel plate air dielectric barrier discharge system using repetitive nanosecond pulses, with pulse repetition frequencies (PRFs) from 0.1 to 100 kHz. It examines how pulse parameters—PRF, pulse number and gas pressure—affect the memory effect, leveraging current and voltage measurements, fast imaging and optical emission spectroscopy. The findings show that higher PRFs lead to a reduction in breakdown voltage well below the streamer breakdown threshold. We argue that this effect may be attributed to the cumulative buildup of metastable species and negative ions in the discharge gap which could sustain free electrons in‐between the voltage pulses. Despite observing increased energy deposition with higher PRFs, the impact on filament formation was minimal, highlighting a strong dependence of discharge morphology on the accumulation of plasma‐produced species. This research studies provides valuable insights for controlling discharge regimes in applications such as plasma‐assisted combustion, surface treatment and air treatment by clarifying the interactions between discharge mechanisms at different PRFs and pressures.
42 ENGINEERING↗