Search NASA⌕ Search

SEARCH · Search NASA

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.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

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↗
DOE OSTI2021

Generation of excited species in a streamer discharge

At or near atmospheric pressure, most transient discharges, particularly in molecular gases or gas mixture containing molecular gases result in a space-charge dominated transport called a streamer discharge. The excited species generation in such discharges forms the basis for plasma chemistry in most technological applications. In this paper we simulate the propagation of streamers in atmospheric pressure N 2 to understand the energy partitioning in the formation of various excited species and compare the results to a uniform Townsend discharge. The model is fully two-dimensional with azimuthal symmetry. Results show a significantly larger fraction of the energy goes into vibrational excitation of the N 2 ground state in a streamer-type discharge in comparison to a Townsend discharge. For lower applied voltages, the anode-directed (negative) steamer is slightly more efficient in channeling energy for excited species production in comparison to a cathode-directed (positive) streamer. Near 70% overvoltage both types of streamers show very similar energy partitioning but quite different from a Townsend discharge.

42 ENGINEERING↗
DOE OSTI2024

Massively parallel axisymmetric fluid model for streamer discharges

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.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

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↗
DOE OSTI2024

Development of a Performance Portable Non-Equilibrium Plasma Fluid Solver on Adaptive Grids

This presentation will describe the numerical techniques, programming paradigms, verification, and performance of a non-equilibrium plasma fluid solver that can effectively utilize current and upcoming central processing and graphics processing unit (CPU+GPU) architectures. Our plasma fluid model solves the conservation equations for self-consistent electrostatic Poisson, electron and heavy species transport, and electron temperature on adaptive Cartesian grids. Our solver is written using performance portable adaptive mesh management library, AMReX (Zhang et al., JOSS, 4 (37) 1370, 2019), and can be built and run on widely available vendor specific GPU architectures (NVIDIA/AMD/Intel). We utilize a non-subcycled second order semi-implicit time-stepping method where all adaptive mesh refinement (AMR) levels are advanced with the same time step. The composite multi-level multigrid solver from within AMReX is used for each of the governing equations that are cast into a Helmholtz equation form. We have also developed a python based chemical mechanism parser framework that uses a similar format as CANTERA (Goodwin et al., Zenodo, 2018) yaml files as input. Our custom parser reads the yaml file and provides C++ files with transport and production rate functions that can be executed on both host (CPU) and device (GPU). 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 low-pressure capacitive and high-pressure streamer discharges. Our initial performance studies indicate 10X speed-up using 20 NVIDIA GPUs versus 200 CPUs for an atmospheric streamer discharge problem solved on a 512 x 1024 x 512 grid.

graphics processing units↗
DOE OSTI2021

Streamer self-focusing in an external longitudinal magnetic field

The numerical simulation of the development of a streamer discharge in a gap with an external longitudinal magnetic field was used to demonstrate the self-focusing of such discharges. Self-focusing is caused by a sharp deceleration of the radial ionization wave due to a change in the electron energy distribution function, a decrease in the average electron energy, the rate of gas ionization, and the electron mobility in crossed electric and magnetic fields as compared to the case of the discharge development without a magnetic field. The self-focusing effect of a streamer discharge in an external longitudinal magnetic field is observed for both positive and negative pulse polarities. The work proposes an estimate of the critical value of the magnetic field, which makes it possible to control the development of pulsed high-voltage discharges at various gas pressures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

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↗
DOE OSTI2023

Jacobian-free Newton–Krylov method for the simulation of non-thermal plasma discharges with high-order time integration and physics-based preconditioning

A preconditioning framework for the numerical simulation of non-thermal streamer discharges is developed using the Jacobian-free Newton-Krylov (JFNK) method. A reduced plasma fluid model is considered, consisting of electrons, one positive ion, one negative ion, and the electrostatic potential. Here, the plasma kinetics model includes ionization, electron-ion recombination, electron attachment, electron detachment, and ion-ion recombination. The governing equations are made dimensionless, discretized in space with finite differences, and integrated in time with a fully implicit method based on high-order backward differentiation formulas. The preconditioning framework is based on a linearized form of the governing equations and physics-based operator splitting. The efficiency of the preconditioning strategy is assessed through two test cases: streamer propagation between parallel plates and an axisymmetric pin-to-pin discharge. The fully implicit approach overcomes traditional restrictions in the time step size due to processes such as electron drift, electron diffusion, and dielectric relaxation. Excellent performance is observed through relevant statistics of the JFNK solver, although the number of linear iterations increases for the pin-to-pin discharge when nonlinear numerical boundary conditions are imposed at the electrodes. Performance studies show scalability with O(100-1000) processors for O(10M) unknowns with ample room for optimization.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Transient behavior of drift and ionization in atmospheric pressure nitrogen discharge

The fluid models are frequently used to describe a non-thermal plasma such as a streamer discharge. The required electron transport data and rate coefficients for the fluid model are parametrized using the local field approximation (LFA) in first order models and the local-mean-energy approximation (LMEA) in second order models. We performed Monte Carlo simulations in Nitrogen gas with step changes in the E/N (reduced electric field) to study the behavior of the transport properties in the transient phase. During the transient phase of the simulation, we extract the instantaneous electron mean energy, which is different from the steady state mean electron energy, and the corresponding transport parameters and rate coefficients. Our results indicate that the mean electron energy is not a suitable parameter for mobility/drift of electrons due to big difference in momentum relaxation and energy relaxation. Furthermore, the high energy threshold rates such as ionization show a strong correlation to mean electron energy. In second order models where the energy-balance equation is solved, we suggest that it would rather be appropriate to use the local electric field to find electron drift velocity in gases such as Nitrogen and the local mean electron energy to determine the ionization and excitation rates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

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↗
DOE OSTI2024

Enhancements of electric field and afterglow of non-equilibrium plasma by Pb(ZrxTi1−x)O3 ferroelectric electrode

Abstract Manipulating surface charge, electric field, and plasma afterglow in a non-equilibrium plasma is critical to control plasma-surface interaction for plasma catalysis and manufacturing. Here, we show enhancements of surface charge, electric field during breakdown, and afterglow by ferroelectric barrier discharge. The results show that the ferroelectrics manifest spontaneous electric polarization to increase the surface charge by two orders of magnitude compared to discharge with an alumina barrier. Time-resolved in-situ electric field measurements reveal that the fast polarization of ferroelectrics enhances the electric field during the breakdown in streamer discharge and doubles the electric field compared to the dielectric barrier discharge. Moreover, due to the existence of surface charge, the ferroelectric electrode extends the afterglow time and makes discharge sustained longer when alternating the external electric field polarity. The present results show that ferroelectric barrier discharge offers a promising technique to tune plasma properties for efficient plasma catalysis and electrified manufacturing.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Upward propagation of gigantic jets revealed by 3D radio and optical mapping

Occasionally, lightning will exit the top of a thunderstorm and connect to the lower edge of space, forming a gigantic jet. Here, we report on observations of a negative gigantic jet that transferred an extraordinary amount of charge between the troposphere and ionosphere (~300 C). It occurred in unusual circumstances, emerging from an area of weak convection. As the discharge ascended from the cloud top, tens of very high frequency (VHF) radio sources were detected from 22 to 45 km altitude, while simultaneous optical emissions (777.4 nm OI emitted from lightning leaders) remained near cloud top (15 to 20 km altitude). This implies that the high-altitude VHF sources were produced by streamers and the streamer discharge activity can extend all the way from near cloud top to the ionosphere. The simultaneous three-dimensional radio and optical data indicate that VHF lightning networks detect emissions from streamer corona rather than the leader channel, which has broad implications to lightning physics beyond that of gigantic jets.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Unraveling the Physics Associated with the Production of Extremely Dense Plasma States of Microscale Nanosecond-pulsed Discharges

The research carried out on this grant aimed to study the properties of anomalously dense plasmas formed in so-called ”fast” discharges, such as high pressure narrow-gap nanosecond pulsed discharges and short pulse (nanosecond-scale) laser driven discharges. These fast discharges are receiving much attention as they now enable new applications in aerodynamics, combustion, biology, and medicine. The mechanisms that generate these high electron densities (ne > 10 19 cm −3 ) that encroach on warm dense matter regimes are not well understood. Such densities are orders of magnitude higher than those (≈ 10 15 cm −3 ) routinely seen in classical high-pressures (streamer) discharges. Furthermore, the dense plasmas are reported to be highly nonequilibrium, with T gas <<, T i << T e . These properties would place these discharges in a regime of near full ionization (Z ≥ 1) where there is inadequate electron screening and therefore strong ion coupling - physics that is not generally accounted for in simulations. In this project work, we seek to study these fast discharges, validate simulations with experiments at pressures as high as 10-50 bar, and augment plasma production to even higher levels of n e (> 10 20 cm −3 ) with second-stage picosecond laser heating. In parallel, we have studied the use of these fast discharges in reactive gases, unraveling mechanisms related to the fast discharge dissociation of CO 2 .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗