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. . 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 of O2 are attributed to the electronic excitation of N2 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. .
|Original language||English (US)|
|Title of host publication||AIAA SciTech Forum 2022|
|Publisher||American Institute of Aeronautics and Astronautics Inc, AIAA|
|State||Published - 2022|
|Event||AIAA Science and Technology Forum and Exposition, AIAA SciTech Forum 2022 - San Diego, United States|
Duration: Jan 3 2022 → Jan 7 2022
|Name||AIAA Science and Technology Forum and Exposition, AIAA SciTech Forum 2022|
|Conference||AIAA Science and Technology Forum and Exposition, AIAA SciTech Forum 2022|
|Period||1/3/22 → 1/7/22|
Bibliographical noteFunding Information:
S. Yang acknowledges the grant support from NSFCBET 2002635. T.S. Taneja acknowledges the support from the UMII MnDRIVE Graduate Assistantship Award. The authors acknowledge Dr. Dezhi Zhou for his contributions towards developing the base solvers which are extended in this work. T.S. Taneja acknowledges the useful discussions with Dr. Maria Castela, Mr. Surya Narayan and Mr. Navneeth Srinivasan.
S. Yang acknowledges the grant support from NSF CBET 2002635. T.S. Taneja acknowledges the support from the UMII MnDRIVE Graduate Assistantship Award. The authors acknowledge Dr. Dezhi Zhou for his contributions towards developing the base solvers which are extended in this work. T.S. Taneja acknowledges the useful discussions with Dr. Maria Castela, Mr. Surya Narayan and Mr. Navneeth Srinivasan.
© 2022, American Institute of Aeronautics and Astronautics Inc.. All rights reserved.