Search for: All records

Abstract Self-consistent 1D modeling of streamers in ammonia-oxygen-nitrogen-water mixtures has been performed in this work. A fluid model that includes species transport, electrostatic potential, and detailed chemistry was developed and verified. This model is then used to simulate the avalanche, streamer formation and propagation phases, driven by a nanosecond voltage pulse, at different thermochemical conditions derived from a 1D laminar premixed ammonia-air flame. The applicability of the Meek’s criterion in predicting the streamer inception location was successfully confirmed. Streamer formation and propagation duration were found to vary significantly with different thermochemical conditions, due to the difference in ionization rates. The thermochemical state also affected the breakdown characteristics which was tested by maintaining the background reduced electric field constant. Detailed kinetic analyses revealed the importance of $\mathrm{O}{(}^1\mathrm{D})$in the production of key radicals, such as O, OH, and NH₂. Furthermore, the contributions of the dissociative electronic excitation of NH₃towards the production of H and NH₂radicals have also been reported. Spatial and temporal evolution of the electron energy loss fractions for various inelastic collision processes at different thermochemical states uncovered the input plasma energy spent of fuel dissociation and the large variability in the dominant processes during the avalanche and streamer propagation phases. The methodology and analyses reported in this work are key towards developing effective strategies for controlled nanosecond-pulsed non-equilibrium plasma sources used for ammonia ignition and flame stabilization. 

Abstract The present work improves a phenomenological plasma-assisted combustion model by integrating the spatiotemporal distribution of plasma power density, thereby considering the evolution of plasma streamers in the modeling, and subsequently, better predicting the ignition kernel evolution. The improved phenomenological model is validated against experiments representing the plasma discharge and post-discharge ignition kernel evolution. Specifically, the new model demonstrates a more accurate prediction of ultrafast gas heating and O₂dissociation during the plasma discharge, compared to the original model. In addition, the new model is found to closely match the experimental pressure wave and heated channel profiles post-discharge without the need for tuning the energy deposition (unlike the original model), highlighting its accuracy of post-discharge ignition kernel dynamics. The improved phenomenological model is then employed to investigate ignition kernel evolution for a stoichiometric methane-air discharge across various discharge gap configurations. Simulations reveal a non-uniform temperature and streamer distribution progressing from the electrode tips toward the center, contrasting uniform cylindrical discharges previously described in the original model. Streamer propagation is observed to be faster for larger gaps when maintained at the same average electric field for different discharge gaps. The tendency of smaller gaps to produce detached toroidal ignition kernels is observed, while larger gaps promote cylindrical and attached ignition kernels. Interactions between successive ignition kernels from consecutive discharges varied significantly, with the smallest gap (1 mm) promoting the quenching of the preceding ignition kernel due to the initial kernel–kernel separation. The intermediate gap (2 mm) promotes detached kernel growth. In contrast, in the largest gap (4 mm), kernels consistently combine and expand attached to electrodes. The impact of homogeneous isotropic turbulence is also explored, showing the persistence of ignition kernels early on but eventually quenching due to enhanced radical and heat losses with pronounced turbulence intensity. 

The impacts of the pulse repetition frequency (PRF), number of pulses, and energy per pulse in a train of nanosecond discharge pulses on the ignition of a flowing lean premixed methane–air mixture are investigated using numerical simulations. A phenomenological plasma model coupled with a compressible reacting flow solver is used for these simulations. The simulation strategy has been well validated by comparing the experimental schlieren and OH planar laser induced fluorescence (PLIF) results with the numerical schlieren (i.e., density gradient) and OH density profiles, respectively. The characteristics of the ignition kernels produced by each discharge pulse and their interaction with each other as functions of the PRF are investigated. Three regimes were defined in the literature based on this interaction of the ignition kernels — fully coupled, partially coupled, and decoupled. This study uses numerical simulations to probe into the constructive and destructive effects, that ultimately determine ignition success, in these different regimes. The complete overlap of kernels and the complete lack of synergy between kernels produced by consecutive pulses are attributed to the success and failure of ignition and flame propagation in the fully coupled and decoupled regimes, respectively. In the partially coupled regime, the convection heat loss driven by the shock-turned-acoustic wave of the next discharge pulse, on the kernel produced by the previous discharge pulse, in addition to diffusion losses, contribute to ignition failure. However, the expansion of the next kernel in a region of higher average temperature and radical concentration created by the previous kernel could help to bridge the gap between the two kernels and result in successful ignition. The important parameters of energy per pulse, number of pulses, and equivalence ratio affect the competition between these constructive and destructive effects, which eventually determines the ignition success in this regime. Finally, the change in the nature of interaction between consecutive kernels from decoupled to partially coupled, at the same frequency but with different energies per pulse, is also shown. 

Natural gas contains a significant fraction of methane, a strong greenhouse gas besides being a potent hydrogen carrier. Thus, reforming methane to a more reactive gas mixture could potentially abate the associated greenhouse heating by depleting methane and provide a pathway to generate hydrogen. The present study investigates the non-equilibrium plasma-assisted reforming of methane to produce hydrogen and reactive alkenes using repetitive nanosecond pulse discharges. A detailed gas-phase chemical kinetics mechanism along with plasma reforming kinetics derived from our previous work are used to perform 0D calculations to obtain the energy fractions for various plasma processes. A phenomenological model for the plasma-assisted reforming of methane/nitrogen mixtures is developed by considering the vibrational energy transport equations of both methane and nitrogen separately. The energy fractions involved in various plasma processes, such as ultra-fast gas heating and ultra-fast gas dissociation due to the electron excitation reactions, and slow gas heating due to the relaxation of vibrational excitation modes of methane and nitrogen, are accounted for in our new phenomenological model using energy fractions derived earlier. The newly developed phenomenological model is then used to perform 3D direct numerical simulation (DNS) of methane reforming diluted with 60% nitrogen in a pin-to-pin electrode configuration with a discharge gap of 1 mm. The effect of pulsing on the evolution of reformed mixture kernels is investigated by comparing two cases: a single-pulsed case with a pulse energy of 0.8 mJ, and another case using 4 pulses at 200 kHz, with a per pulse energy of 0.2 mJ. The single-pulsed case was observed to promote kernel separation and higher fractions of reformed products, while the multiple-pulsed case resulted in a more diffused kernel. 

The backward problem of plasma assisted combustion emphasizes evaluating the effect of the evolving thermochemical state on the plasma discharge. This paper investigates the dependence of avalanche to streamer to spark formation dynamics and kinetics on the gas composition and temperature at different points in an ammonia-air premixed laminar flame using a self-consistent multigrid-based 1D plasma solver. Different values of alpha, the coefficient for effective ionization events per unit length, have been reported for electron avalanches in air and stoichiometric NH3-air mixtures. The streamer inception has been shown to obey the Meek’s criterion. An exponential reduction in streamer and spark formation time has been observed from plasma simulations at different points in the unburnt, pre-heat zone, reaction zone and the fully burnt regions of the premixed flame. While the enhancement of the reduced electric field with increasing temperature affects effective ionization, there exists a minimum breakdown field for streamer formation, which does not vary proportionally with the changing number density of the gas. The change in the mixture from reactants (NH3, O2, N2) to products of complete combustion of ammonia in air (N2, H2O) has also been shown to affect the streamer and spark formation. Finally, the major pathways during the streamer and spark phases which are responsible for producing important radicals used in combustion of NH3 are also discussed. 

The Global Pathway Analysis (GPA) algorithm helps analyze the chemical kinetics of complex combustion systems by identifying important global reaction pathways connecting a source species to a sink species through various important intermediate species (i.e., hub species). The present work aims to extend GPA algorithm to plasma-assisted combustion and fuel reforming systems to identify the dominant global pathways in such systems at various conditions. In addition, the present study extends the ability of GPA algorithm to identify reaction cycles involving the excitation of high-concentration species (e.g., O2, N2, and fuel) to their vibrational and electronic states and the subsequent de-excitation to their ground state, based on their significance on the reactivity of plasma-assisted systems in terms of gas heating and radical production. Provisions are made in the GPA algorithm to evaluate the reactivity of identified reaction pathways and cycles based on the element-flux transfer (i.e., dominance), heat release, and radical production rate. The newly developed Plasma-based Global Pathway Analysis (PGPA) algorithm is then used to analyze the plasma-assisted combustion of ammonia and reforming of methane. The PGPA analyses elucidated the significance of vibrational-translational cycles on the reactivity of NH3/air mixtures. Further, analyses on the production of NO ascribed the early reforming of NH3 to N2 and H2 in impeding the production of NO during plasma-assisted NH3 ignition. Lastly, the enhanced reforming of CH4/N2 mixtures using plasma has been attributed to electron impact dissociation of CH4 when compared to thermal reforming. In contrast, conventional path-Flux analysis (PFA) was found to require significant manual effort and pre-analysis intuitions from expert knowledge, making it arduous to provide valuable insights into plasma chemistry. The user-friendly and automated nature of PGPA thus provides a valuable tool for assessing the kinetics of plasma-assisted systems helpful in analyzing and, further, a foundation in reducing plasma-assisted chemistry, without the needs of expert knowledge. 

Natural gas associated with oil wells and natural gas fields is a significant source of greenhouse gas emissions and airborne pollutants. Flaring of the associated gas removes greenhouse gases like methane and other hydrocarbons. The present study explores the possibility of enhancing the flaring of associated gas mixtures (C1 – C4 alkane mixture) using nanosecond pulsed non-equilibrium plasma discharges. Starting with a detailed chemistry for C0 – C4 hydrocarbons (Aramco mechanism 3.0 – 589 species), systematic reductions are performed to obtain a smaller reduced mechanism (156 species) yet retaining the relevant kinetics of C1 – C4 alkanes at atmospheric pressure and varying equivalence ratios (φ = 0.5 – 2.0). This conventional combustion chemistry for small alkanes is then coupled with the plasma kinetics of CH4, C2H6, C3H8, and N2, including electron-impact excitations, dissociations, and ionization reactions. The newly developed plasma-based flare gas chemistry is then utilized to investigate repetitively pulsed non-equilibrium plasma-assisted reforming and subsequent combustion of the flare gas mixture diluted with N2 at different conditions. The results indicate an enhanced production of hydrogen, ethylene and other species in the reformed gas mixture, owing to the electron-impact dissociation pathways and subsequent H-abstractions and recombination reactions, thereby resulting in a mixture of CH4, H2, C2H4, C2H2, and other unsaturated C3 species. The reformed mixture shows an enhanced reactivity as exhibited by their shorter ignition delays. The reformed mixture is also observed to undergo increased methane destruction and higher equilibrium temperatures compared to the original mixture as the gas temperature increases, thereby exhibiting a potential for reducing the unburnt emissions of methane and other hydrocarbons. 

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. 

Global Pathway Selection/Analysis (GPSA) algorithm helps in analyzing the chemical kinetics of complex combustion systems by identifying important global reaction pathways that connects a source and a sink species. The present work aims to extend the application of GPSA to plasma assisted combustion systems in order to identify the dominant global pathways that govern the plasma and combustion kinetics at various conditions. The reaction cycles involving the excitation of nitrogen to its vibrational and electronic states and the subsequent de-excitation to its ground state are found to control the reactivity of plasma assisted systems. Provisions are made in the GPSA algorithm to capture the dominant reaction pathways and cycles of plasma assisted combustion (i.e., p-GPSA). Further, the analysis of plasma assisted ammonia combustion are presented as an example, which includes the results obtained using both the traditional path flux analysis and p-GPSA. The dominant pathways for the plasma assisted combustion of ammonia are identified along with the dominant excitation--de-excitation loops and their importance are ascertained and verified using path flux analysis.