An official website of the United States government
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.
more »
« less
Comparing Low-Mach and Fully-Compressible CFD Solvers for Phenomenological Modeling of Nanosecond Pulsed Plasma Discharges with and without Turbulence
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
- Award ID(s):
- 2002635
- PAR ID:
- 10387572
- Date Published:
- Journal Name:
- AIAA Scitech 2022 Forum
- Page Range / eLocation ID:
- AIAA 2022-0976
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Nanosecond pulsed plasmas have been demonstrated, both experimentally and numerically, to be beneficial for ignition, mainly through gas heating (at different timescales) and radical seeding. However, most studies focus on specific gas conditions, and little work has been done to understand how plasma performance is affected by fuel and oxygen content, at different gas temperatures and deposited energies. This is relevant to map the performance envelope of plasma-assisted combustion across different regimes, spanning from fuel-lean to fuel-rich operation, as well as oxygen-rich to oxygen-vitiated conditions, of interest to different industries. This work presents a computational effort to address a large parametric exploration of combustion environments and map out the actuation authority of plasmas under different conditions. The work uses a zero-dimensional plasma-combustion kinetics solver developed in-house to study the ignition of CH4/O2/N2 mixtures with plasma assistance. A main contribution of the study is the detailed tracking of the energy, from the electrical input all the way to the thermal and kinetic effects that result in combustion enhancement. This extends prior works that focus on the first step of the energy transfer: from the electrical input to the electron-impact processes. Independent of the composition, four pathways stand out: (i) vibrational-translational relaxation, (ii) fast gas heating, (iii) O2 dissociation, and (iv) CH4 dissociation. Results show that the activated energy pathways are highly dependent on gas state, composition, and pulse shape, and can explain the observed range in performance regarding ignition enhancement. The approach can be used to calculate the fractional energy deposition into the main pathways for any mixture or composition, including new fuels, and can be a valuable tool to construct phenomenological models of the plasma across combustion environments.more » « less
-
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 O2dissociation 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.more » « less
-
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.more » « less
-
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
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Conference Paper:
- https://doi.org/10.2514/6.2022-0976
