skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 11:00 PM ET on Friday, July 12 until 2:00 AM ET on Saturday, July 13 due to maintenance. We apologize for the inconvenience.

Title: Numerical Investigation of Ignition Kernel Development with Nanosecond Pulsed Plasma in Quiescent and Flowing Mixtures
Nanosecond Pulsed High Frequency Discharges (NPHFD) are gaining popularity over conventional spark and arc discharges as they have been shown to increase energy efficiency, enhance ignition probability and sustained kernel growth, and offer more flexibility and control for ignition applications under various conditions. Hence, it is important to determine the impact of different factors such as the optimal pulse energy, background flow conditions, inter-pulse time, mixture equivalence ratio, etc. on the success of ignition of premixed mixtures with NPHFD. This work presents a numerical investigation of the morphology of ignition kernel development with both single-pulse and multiple-pulse discharges. Nanosecond non-equilibrium plasma discharges are modeled between pin-pin electrodes in a subsonic ignition tunnel with quiescent and flowing premixed mixtures of methane and air. Large eddy simulations (LES) are conducted to investigate the reasons for successful and failed ignition in different scenarios. A single pulse discharge in the presence of electrodes, in a quiescent medium, elucidates the gas recirculation pattern caused by the plasma pulse which results in a separated toroidal kernel from the primary ignition kernel between the electrodes. Convection heat loss to the mean flow results in quenching of the high temperature, radical-rich hot-spots creeping on the electrode walls, and leaving only the semi-toroidal kernel to propagate downstream. Finally, simulations with multiple pulses with different inter-pulse times have been conducted to analyze the synergistic effect of overlapping kernels with high temperature and OH concentration, which has been attributed as the primary reason for higher ignition probabilities in the “fully coupled” regime reported in the experiments. Successful ignition kernel formation is reported with 3 pulses at a pulse repetition frequency of 300 kHz in the fully coupled regime. This kernel volume was almost 4 times, and develops in two-thirds the time, compared to the ignition kernel volume formed by the single pulse discharge with the same total energy. Ten pulses with twice as much total energy were deposited at a much lower frequency of 2 kHz, which resulted in disjoint hot-spots that fail to form an ignition kernel in the decoupled regime.  more » « less
Award ID(s):
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
American Institute of Aeronautics and Astronautics
Date Published:
Journal Name:
Medium: X
National Harbor, MD & Online
Sponsoring Org:
National Science Foundation
More Like this
  1. In this work, evolution of parameters of nanosecond repetitively pulsed (NRP) discharges in pin-to-pin configuration in air was studied during the transient stage of initial 20 discharge pulses. Gas and plasma parameters in the discharge gap were measured using coherent microwave scattering, optical emission spectroscopy, and laser Rayleigh scattering for NRP discharges at repetition frequencies of 1, 10, and 100 kHz. Memory effects (when perturbations induced by the previous discharge pulse would not decay fully until the subsequent pulse) were detected for the repetition frequencies of 10 and 100 kHz. For 10 kHz NRP discharge, the discharge parameters experienced significant change after the first pulse and continued to substantially fluctuate between subsequent pulses due to rapid evolution of gas density and temperature during the 100  μs inter-pulse time caused by intense redistribution of the flow field in the gap on that time scale. For 100 kHz NRP discharge, the discharge pulse parameters reached a new steady-state at about five pulses after initiation. This new steady-state was associated with well-reproducible parameters between the discharge pulses and substantial reduction in breakdown voltage, discharge pulse energy, and electron number density in comparison to the first discharge pulse. For repetition frequencies 1–100 kHz considered in this work, the memory effects can be likely attributed to the reduction in gas number density and increase in the gas temperature that cannot fully recover to ambient conditions before subsequent discharge pulses. 
    more » « less
  2. This work aims at comparing the accuracy and overall performance of a low-Mach CFD solver and a fully-compressible CFD solver for direct numerical simulation (DNS) of nonequilibrium plasma assisted ignition (PAI) using a phenomenological model described in Castela et al. [1]. The phenomenological model describes the impact of nanosecond pulsed plasma discharges by introducing source terms in the reacting flow equations, instead of solving the detailed plasma kinetics at every time step of the discharge. Ultra-fast gas heating and dissociation ofO2 are attributed to the electronic excitation ofN2 and the subsequent quenching to ground state. This process is highly exothermic, and is responsible for dissociation of O2 to form O radicals; both of which promote faster ignition. Another relatively slower process of gas heating associated with vibrational-to-translational relaxation is also accounted for, by solving an additional vibrational energy transport equation. A fully-compressible CFD solver for high Mach (M>0.2) reacting flows, developed by extending the default rhoCentralFoam solver in OpenFOAM, is used to perform DNS of PAI in a 2D domain representing a cross section of a pin-to-pin plasma discharge configuration. The same case is also simulated using a low-Mach, pressure-based CFD solver, built by extending the default reactingFoam solver. The lack of flow or wave dominated transport after the plasma-induced weak shock wave leaves the domain causes inaccurate computation of all the transport variables, with a rather small time step dictated by the CFL condition, with the fully-compressible solver. These issues are not encountered in the low-Mach solver. Finally, the low-Mach solver is used to perform DNS of PAI in lean, premixed, isotropic turbulent mixtures of CH4-air at two different Reynolds numbers of 44 and 395. Local convection of the radicals and vibrational energy from the discharge domain, and straining of the high temperature reaction zones resulted in slower ignition of the case with the higher Re. A cascade effect of temperature reduction in the more turbulent case also resulted in a five - six times smaller value of the vibrational to translational gas heating source term, which further inhibited ignition. Two pulses were sufficient for ignition of the Re = 44 case, whereas three pulses were required for the Re = 395 case; consistent with the results of Ref. [1]. 
    more » « less
  3. A robust and efficient framework for simulating NSPD in multiple dimensions is developed. The reactive Navier-Stokes equations are extended to include a drift-diffusion plasma-fluid model with a local field approximation (LFA) in a finite-volume solver, which uses an adaptive mesh refinement (AMR) strategy to address the wide separation of length scales in the problem. A two-way coupling strategy is used whereby the plasma-fluid model and reactive Navier-Stokes equations are integrated simultaneously. The oxidation of ethylene/air mixtures mediated by NSPD is simulated in a pin-to-pin configuration. All phases of the plasma discharge are simulated explicitly (including streamer ignition, propagation, and connection, as well as the subsequent spark phase), along with the evolution of the plasma during the inter-pulse period. Temporally and spatially-resolved results are presented, with an emphasis on the analysis of heating and energy deposition, as well as of the evolution of the concentration of active particles generated during the NSPD and their influence on ignition. 
    more » « less
  4. null (Ed.)
    The influence of pulse frequency (1−60 kHz) in a nanosecond filamentary discharge propagating along a flowing liquid water film was assessed with regards to the formation of chemical species with argon and helium carrier gasses. The production rate and energy yield for H₂O₂ and H₂ were measured for both carrier gases, and O₂ formation was determined for helium. The effect of pulse frequency on the energy dissipated per pulse as well as electron density was also investigated. The results indicate that the energy yield for H₂O₂ decreases with increasing pulse frequency while the energy yields of H₂ and O₂ remain relatively unaffected. It is proposed that the difference in the trends of the liquid versus the significantly longer residence time of the liquid phase allowing for more degradation of formed hydrogen peroxide before it is able exit the reactor. 
    more » « less
  5. Abstract

    Two dimensional arrays of streamer discharges were developed using electrical discharge machining of stainless steel sheets and stacking them together with spacers to allow gas flow between the sheets. A nanosecond pulsing circuit, capable of delivering 2–40 kV pulses with pulse widths of >20 ns by using two spark gaps as switches, was developed as a simple tunable pulsing power supply. High resolution imaging of the plasma for uniformity across the array tips was conducted. Optical emission spectroscopy was used to characterize the species created as well as probe the temperature of the discharge for various substrates, voltage pulse durations, voltage pulse magnitudes, and gas flows. The discharge properties were found to be substrate independent for a wide variety of conditions.

    more » « less