skip to main content

Title: 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 more » 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]. « less
Award ID(s):
Publication Date:
Journal Name:
AIAA Scitech 2022 Forum
Page Range or eLocation-ID:
AIAA 2022-0976
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The ignition of plasmas in liquids has applications from medical instrumentation to manipulation of liquid chemistry. Formation of plasmas directly in a liquid often requires prohibitively large voltages to initiate breakdown. Producing plasma streamers in bubbles submerged in a liquid with higher permittivity can significantly lower the voltage needed to initiate a discharge by reducing the electric field required to produce breakdown. The proximity of the bubble to the electrodes and the shape of the bubbles play critical roles in the manner in which the plasma is produced in, and propagates through, the bubble. In this paper, we discuss results from a three-dimensional direct numerical simulation (DNS) used to investigate the shapes of bubbles formed by injection of air into water. Comparisons are made to results from a companion experiment. A two-dimensional plasma hydrodynamics model was then used to capture the plasma streamer propagation in the bubble using a static bubble geometry generated by the DNS The simulations showed two different modes for streamer formation depending on the bubble shape. In an elliptical bubble, a short electron avalanche triggered a surface ionization wave (SIWs) resulting in plasma propagating along the surface of the bubble. In a circular bubble, anmore »electron avalanche first traveled through the middle of the bubble before two SIWs began to propagate from the point closest to the grounded electrode where a volumetric streamer intersected the surface. In an elliptical bubble approaching a powered electrode in a pin-to-pin configuration, we experimentally observed streamer behavior that qualitatively corresponds with computational results. Optical emission captured over the lifetime of the streamer curve along the path of deformed bubbles, suggesting propagation of the streamer along the liquid/gas boundary interface. Plasma generation supported by the local field enhancement of the deformed bubble surface boundaries is a mechanism that is likely responsible for initiating streamer formation.

    « less
  2. One of the key factors in simulating realistic wall-bounded flows at high Reynolds numbers is the selection of an appropriate turbulence model for the steady Reynolds Averaged Navier–Stokes equations (RANS) equations. In this investigation, the performance of several turbulence models was explored for the simulation of steady, compressible, turbulent flow on complex geometries (concave and convex surface curvatures) and unstructured grids. The turbulence models considered were the Spalart–Allmaras model, the Wilcox k- ω model and the Menter shear stress transport (SST) model. The FLITE3D flow solver was employed, which utilizes a stabilized finite volume method with discontinuity capturing. A numerical benchmarking of the different models was performed for classical Computational Fluid Dynamic (CFD) cases, such as supersonic flow over an isothermal flat plate, transonic flow over the RAE2822 airfoil, the ONERA M6 wing and a generic F15 aircraft configuration. Validation was performed by means of available experimental data from the literature as well as high spatial/temporal resolution Direct Numerical Simulation (DNS). For attached or mildly separated flows, the performance of all turbulence models was consistent. However, the contrary was observed in separated flows with recirculation zones. Particularly, the Menter SST model showed the best compromise between accurately describing the physicsmore »of the flow and numerical stability.« less
  3. Direct Numerical Simulation (DNS) of compressible spatially-developing turbulent boundary layers (SDTBL) is performed at a Mach number of 2.5 and low/high Reynolds numbers over isothermal Zero-Pressure Gradient (ZPG) flat plates. Turbulent inflow information is generated via a dynamic rescaling-recycling approach (J. Fluid Mech., 670, pp. 581-605, 2011), which avoids the use of empirical correlations in the computation of inlet turbulent scales. The range of the low Reynolds number case is approximately 400-800, based on the momentum thickness, freestream velocity and wall viscosity. DNS at higher Reynolds numbers (~3,000, about four-fold larger) is also carried out with the purpose of analyzing the effect of Reynolds number on the transport phenomena in the supersonic regime. Additionally, low/high order flow statistics are compared with DNS of an incompressible isothermal ZPG boundary layer at similar low Reynolds numbers and the temperature regarded as a passive scalar. Peaks of turbulence intensities move closer to the wall as the Reynolds number increases in the supersonic flat plate. Furthermore, Reynolds shear stresses depict a much larger "plateau" (constant shear layer) at the highest Reynolds number considered in present study.
  4. Direct Numerical Simulation (DNS) of turbulent spatially-developing boundary layers is performed over an isothermal flat plate at several flow regimes: incompressible, supersonic (Mach 2.5), and hypersonic (Mach 5). Similar low Reynolds numbers are considered in all cases with the purpose of assessing flow compressibility on low/high order flow statistics and on the dynamics of coherent structures of Zero Pressure Gradient (ZPG) flows. Turbulent inflow information is generated by following the concept of the rescaling-recycling approach introduced by Lund et al. (J. Comp. Phys. 140, 233-258, 1998); although, the proposed methodology is extended to high-speed flows. Furthermore, a dynamic approach is employed to connect the friction velocities at the inlet and recycle stations (i.e., there is no need of an empirical correlation as in Lund et al.). The Mach number effect has been mainly identified as significant changes in peak values of the streamwise velocity fluctuations. The vertical transport of Reynolds shear stresses is slightly away from the wall in the near wall region for the hypersonic case. Zones of low speed fluid exhibits a much more elongated shape in incompressible flow as compared with the compressible counterpart. Furthermore, low speed streaks exhibit a contorted, twisted and stretched form in incompressiblemore »flow while they are shorter and more isotropic in the supersonic flow.« less
  5. 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 memorymore »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.« less