Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “pressure gradient”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

The Dependence of Electron-only Magnetic Reconnection on Plasma Beta and Magnetic Shear: Theory and Observations

We investigate the onset conditions for electron-only magnetic reconnection through an analysis of 86 reconnection events observed by the Magnetospheric Multiscale mission in magnetosheath and magnetospheric current sheets. Our study presents the first observational evidence that the occurrence of electron-only reconnection depends on both the magnetic-field shear angle (θ) across the current sheet and the cross-sheet electron beta jump, quantified either by the inflow electron beta difference (Δβ e ) or by a modified electron beta jump (Δβ$^*_e$) defined by the cross-sheet electron pressure difference normalized by the average of the inflow magnetic pressure in the two inflow regions. We extend the diamagnetic-drift suppression mechanism from standard reconnection to electron-only reconnection by deriving relations applicable to electron-scale current sheets. Using scaling estimates, the electron-only reconnection suppression relation is derived as Δβ$^*_e$ ≳ 2(L/d e ) tan (θ/2), while a natural extension of the ion-coupled reconnection suppression criterion is Δβ e ≳ 2(L/d e ) tan (θ/2). At low Δβ e and Δβ$^*_e$ values, reconnection is detected across a large range of magnetic shear (from low to high), whereas high Δβ e and Δβ$^*_e$ conditions restrict reconnection to large shear configurations: shear angles as small as θ ≈ 5° occur when Δβ e (or Δβ$^*_e$) ≲0.1. However, only events exceeding θ ≳ 70° are found when Δβ e ≳ 1.0; this apparent restriction weakens when employing Δβ$^*_e$, as the corresponding data points shift toward smaller values of the modified electron beta jump. These findings quantitatively support theoretical predictions of reconnection suppression in high-β e plasmas at low shears, attributed to X-line drift induced by pressure gradients across electron-scale current sheets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Toward Higher-order Accuracy in Self-gravitating Hydrodynamics

High-order algorithms have emerged in numerical astrophysics as a promising avenue to reduce truncation error (proportional to a power of the linear resolution Δ x ) with only a moderate increase to computational expense. Significant effort has been placed in the development of finite-volume algorithms for (magneto)hydrodynamics; however, state-of-the-art astrophysical simulations tightly couple a plenitude of physics, additionally including gravity, photon transport, cosmic-ray transport, chemistry, and/or diffusion, to name a few. Algorithms frequently operator-split this additional physics (often a first-order error in time) and/or adopt a model wherein their evaluation is limited to second-order accuracy in space. In this work, we present a fourth-order-accurate finite-volume scheme for self-gravitating hydrodynamics on a uniform Cartesian grid. The method supplies source terms for the gravitational acceleration ( ρ g ) and gravitational energy release ( ρ v · g ) associated with fourth-order-accurate solutions to the Poisson equation. Our scheme (1) guarantees the conservation of total linear momentum while (2) decreasing (in proportion to Δ x 4 ) the effects of spurious heating and/or cooling associated with truncation error in the gravity. We demonstrate expected convergence rates for the algorithm by measuring errors in test problems evolving self-gravity modified linear waves and 3D polytropic equilibria. We test robustness of the algorithm by integrating an induced “inside-out” adiabatic collapse. We also discuss a method to smoothly downgrade the solution to second-order spatial accuracy to avoid spurious overshoots near steep density and/or pressure gradients.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Planet-induced Vortices with Dust Coagulation in Protoplanetary Disks

In this work, we study how the dust coagulation/fragmentation will influence the evolution and observational appearances of vortices induced by a massive planet embedded in a low-viscosity disk by performing global 2D high-resolution hydrodynamical simulations. Within the vortex, due to its higher gas surface density and steeper pressure gradients, dust coagulation, fragmentation, and drift (to the vortex center) are all quite efficient, producing dust particles ranging from 1 μm to ~1.0 cm, as well as an overall high dust-to-gas ratio (above unity). In addition, the dust size distribution is quite nonuniform inside the vortex, with the mass-weighted average dust size at the vortex center (~4.0 mm) being a factor of ~10 larger than other vortex regions. Both large (~millimeter) and small (tens of microns) particles contribute strongly to affect the gas motion within the vortex. As such, we find that the inclusion of dust coagulation has a significant impact on the vortex lifetime and the typical vortex lifetime is about 1000 orbits. Furthermore, after the initial gaseous vortex is destroyed, the dust spreads into a ring with a few remaining smaller gaseous vortices with a high dust concentration and a large maximum size (~millimeter). At late time, the synthetic dust continuum images for the coagulation case show as a ring inlaid with several hot spots at the 1.33 mm band, while only distinct hot spots remain at 7.0 mm.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Numerical Treatment of Shock-induced Nuclear Burning in Double Detonation Type Ia Supernovae

We present a benchmark problem to assess the treatment of shock-induced nuclear burning in the context of double detonation Type Ia supernovae. In a stratified white dwarf model, we implement a shock-detection criterion that suppresses burning in zones characterized by compression and significant pressure gradients, controlled by a tunable parameter, f shock . One-dimensional simulations, using the open-source Castro suite, were conducted across three treatments—burning fully enabled, and burning suppressed with f shock = 2/3 and f shock = 1—across three spatial resolutions (5.0, 2.5, and 0.3125 km). At the finest resolution, the burning-enabled and f shock = 1 models converge, while the f shock = 2/3 front continues to show slight offset behavior. Since most simulations are carried out at much lower resolutions, our tests support the idea that burning in shocks should always be disabled in practice. We also observe that the behavior of lower-resolution simulations remains extremely sensitive to the choice of f shock .

GPU computing↗
DOE OSTI2025

Dry and warm conditions in Australia exacerbated by aerosol reduction in China

A substantial decline in anthropogenic aerosols in China has been observed since the initiation of clean air actions in 2013. This study reveals a linkage between aerosol reductions in China and drier and warmer conditions in Australia. Aerosol decline in China trigger alterations in temperature and pressure gradients between the two hemispheres, leading to intensified outflow from Asia towards the South Indian Ocean, strengthening the Southern Indian Subtropical High and its related Southern Trade Winds. Consequently, this atmospheric pattern results in a moisture divergence over Australia. The reduction in surface moisture further results in more surface energy being converted into sensible heat instead of evaporating as latent heat, warming the near-surface air. The intensified dry and warm climate conditions further cause the increase in wildfire risks during fire seasons in Australia. Our study illuminates the potential impact of distant aerosols on precipitation and temperature variations in Australia, offering valuable insights for drought and wildfire risk mitigation in Australia.

aerosol-cloud interaction↗
DOE OSTI2020

A new form of the Saint-Venant equations for variable topography

The solution stability of river models using the one-dimensional (1D) Saint-Venant equations can be easily undermined when source terms in the discrete equations do not satisfy the Lipschitz smoothness condition for partial differential equations. Although instability issues have been previously noted, they are typically treated as model implementation issues rather than as underlying problems associated with the form of the governing equations. This study proposes a new reference slope form of the Saint-Venant equations to ensure smooth slope source terms and eliminate one source of potential numerical oscillations. It is shown that a simple algebraic transformation of channel geometry provides a smooth reference slope while preserving the correct cross-section flow area and the total Piezometric pressure gradient that drives the flow. The reference slope method ensures the slope source term in the governing equations is Lipschitz continuous while maintaining all the underlying complexity of the real-world geometry. The validity of the mathematical concept is demonstrated with the open-source Simulation Program for River Networks (SPRNT) model in a series of artificial test cases and a simulation of a small urban creek. Validation comparisons are made with analytical solutions and the Hydrologic Engineering Center's River Analysis System (HEC-RAS) model. The new method reduces numerical oscillations and instabilities without requiring ad hoc smoothing algorithms.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Effect of Initial Water Saturation on The Performance of Fracturing Fluids With and Without Polyallylamine under Simulated EGS Conditions

Objectives/Scope: StimuFrac (US Patents 9,873,828 B2 and 9,447,315 B2), a CO2-reactive polymer aqueous solution [polyallylamine (PAA) 1wt% in water] combined with CO2, can be used as a less water-intensive fracturing fluid for enhanced geothermal systems (EGS). Our previous results show that in hot dry rock (HDR), PAA/CO2 fracturing fluids outperformed other fluids such as water, CO2, and CO2/water in generating large fractures with less fluid consumed. The objective of this work is to study the effect of initial water saturation of rock on the performance of StimuFrac fluid in ½ foot cubic rock samples and under representative EGS pressure/temperature conditions using cyclic and constant flow rate injection strategies. The fracturing results are compared with results using different fracturing fluids in terms of controlling fracture propagation rates, fracture hydraulic conductivity, breakdown pressures and volumes of fluids required. Methods/Procedures/Process: In all tests, water was initially injected into the rock to increase the water saturation before the fracturing processes to simulate actual geothermal reservoir conditions. For the cyclic injection, one complete cycle consisted of (1) a PAA slug (or water slug) injection followed by (2) CO2 injection to initiate the fracture. In the second step of the first cycle, the pressure of CO2 is increased until a maximum pressure is reached (fracture is initiated at this moment), and then the injection of CO2 is allowed for another 30 seconds to propagate the fracture. Then, the two-step cycle of PAA followed by CO2 injection (up to 2-4 mL/min) was continued. For the constant flow rate injection strategies, the initial water saturation was increased by injecting water at 1000 psi and 200°C for three days. After that, an initial slug of water, CO2, or PAA was injected and then fracturing was initiated and propagated by injecting CO2 at a constant flow rate. Applications/Significance/Novelty: The results of this study suggest that water saturation, especially near the wellbore region, will significantly affect the fracturing fluid transmission into the rock porous media by changing the relative permeability of CO2 or water, thus affecting the fracture initiation and propagation. In this study, fracturing with cyclic injection or constant flow rate injection methods were performed using three different kinds of fluids systems. These fluids are water, CO2, or CO2 with PAA. Splitting the rock samples in half after fracturing reveals that the fracture propagation is significantly limited under these high water saturated conditions compared to dry initial conditions: The fractures propagate less than 1/3 length of the distance from the wellbore to rock surface, and in some cases no fracture is generated. This may be caused by the fact that leak-off is dominating the fracturing process and the injected fluid flow rate is not high enough to overcome the leak-off even under high flow rate injection conditions. Additionally, CO2 could be leaking off into the wellbore annulus and this may be making it more difficult to generate pressure gradients away from the near-wellbore region.

Jian, Guoqing↗
DOE OSTI2021

Enhanced heavy oil recovery by low salinity polymer flood combined with microgel treatment

“Heavy oil resources account for a large portion of the total oil reserves around the world. The target heavy oil reservoir is located on Alaska’s North Slope (ANS). Advantages of low-salinity HPAM polymer (LSP) over high-salinity polymer (HSP) were demonstrated. LSP could recover more oil with 40% less polymer consumption. No additional oil was recovered by HSP after LSP flood. The first-ever polymer flood pilot on ANS showed remarkable success regarding water cut reduction, oil production increase, delayed breakthrough, and projected oil recovery improvement. Polymer alone was insufficient to achieve satisfactory recovery as the reservoirs were highly heterogeneous. Microgels could improve the effectiveness of polymer flood by reducing water cut and increasing oil recovery. Favorable working conditions were identified. Microgel transport behavior was studied using superpermeable sandpacks (27-221 darcies) with multiple pressure sensors. The particle-to-pore matching size ratio significantly impacted the effectiveness of the gels. A threshold differential pressure (ΔPth) and critical pressure gradient (∇Pcr) were required to push the gels to penetrate and propagate through the channels. The ΔPth and ∇Pcr revealed the underlying mechanisms of selective penetration/placement behavior of microgels in heterogeneous reservoirs. Diagrams were developed to estimate the maximum propagation distance of the gels in channels in conceptual field applications. Sandwich-like channel models and methodologies were developed to comprehensively evaluate the effectiveness of gel materials. Gel retention in the channels was quantified. Results also indicated that the retained gels were dehydrated. Fluid diversion and sweep improvement after gel treatments were evaluated by tracer tests”--Abstract, page iv.

Zhao, Yang↗
DOE OSTI2023

Sorbent system for removing ammonia and organic compounds from a gaseous environment

A first process and sorbent for removing ammonia from a gaseous environment, the sorbent comprised of graphene oxide having supported thereon at least one compound selected from metal salts, metal oxides and acids, each of which is capable of adsorbing ammonia. A second process and sorbent system for removing ammonia and a volatile organic compound from a gaseous environment; the sorbent system comprised of two graphene-based materials: (a) the aforementioned graphene oxide, and (b) a nitrogen and oxygen-functionalized graphene. The sorbents are regenerable under a pressure gradient with little or no application of heat. The processes are operable through multiple adsorption-desorption cycles and are applicable to purifying and revitalizing air contaminated with ammonia and organic compounds as may be found in spacesuits, aerospace cabins, underwater vehicles, and other confined-entry environments.

Junaedi, Christian↗
DOE OSTI2023

Scaling for Nuclear Reactor System: Overview

A scaling analysis is essential for predicting the behavior of a reactor system under various operating conditions, as well as for ensuring safe and reliable experimentation and adequate dataset generation for supporting the system code model development and assessment in accordance with USNRC guidelines 1.203. Testing and modeling can be carried out at various scales, from small test units to full-scale reactors, to validate the predictions produced by the analysis. However, due to the cost and challenges associated with full-scale testing, the nuclear industry uses scaled test facilities—integral effect test (IET) and separate effect test (SET)—to identify potential design issues and optimize the design to ensure regulatory compliance for safer operation and anticipated accident management. In general, IET facilities are non-nuclear thermal hydraulics test facilities. In contrast, SET facilities can be non-nuclear testing (e.g., steam generator, containment, pump system SETs) and nuclear testing (mostly fuel-material irradiation testing). IETs are used to study the system-level thermal-hydraulic behavior of the reactor coolant system (RCS) under different anticipated operational occurrences (AOOs) and design basis accidents (DBAs) such as loss-of-coolant-accident (LOCA) and loss-of-power (i.e., station blackout) and loss-of-flow (i.e., failure of RCS pump). Several LOCA scenarios could happen in reactor DBAs, such as main steam line break (MSLB), SG tube rapture (SGTR), break in direct vessel injection (DVI), and other small breaks LOCA. Scaling analysis determines the appropriate size and operating conditions for the test facility. This involves analyzing the system's thermal-hydraulic behavior at different scales and identifying the key parameters that affect the thermal-hydraulic behavior of the reactor. These key parameters can be identified so that the scaled facility can replicate the expected conditions in a prototypical setting. The scaled facility can simulate the same flow rates, pressure gradients, heat transfer coefficients, and fluid properties as the prototype. Scaling analysis is used to ensure that the results from these tests can be scaled up to predict the behavior of a full-scale prototype with an acceptable range of scaling distortions.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Experimental Study on Flow Condensation of Low Global Warming Potential Refrigerants in a Micro-fin Aluminum Tube

This study investigates the dynamic shifts in refrigerant technologies driven by environmental regulations, particularly emphasizing low global warming potential (GWP). Moreover, there is a rising trend in the adoption of aluminum tubes with internal axial micro-fin structures in heat exchangers to reduce costs. The research focuses on the condensation process within an expanded axial micro-fin aluminum tube with a 5.96 mm fin-tip diameter. Various refrigerants are analyzed, including both single compounds (R-32, R-1234yf, R-1234ze(E)) and zeotropic mixtures (R-454B, R-454C, R-455A). Experimental procedures cover a range of condensation temperatures (35~45 °C), reduced pressures (0.21~0.55), and mass fluxes (150~350 kg/(m² s)), providing crucial data on heat transfer coefficients (HTC) and frictional pressure gradient (FPG). This data is particularly significant for high-glide refrigerants and is instrumental in the design of advanced air conditioning and refrigeration systems aimed at mitigating global warming.

Hu, Yifeng↗
DOE OSTI2025

Sockeye Code Validation Against UNIST Heat Pipe - Poster - Internship 2025

Sockeye is an engineering level code being developed under the MOOSE framework for modeling heat pipes. A heat pipe is a sealed tube with a wick structure filled with a working fluid that transfers heat efficiently and passively via phase change of the working fluid. Heat applied to the evaporator end creates a pressure gradient which causes vapor to migrate to the condenser end where the vapor deposits its energy and condenses. Capillary force generated by the wick structure then draws the liquid back to the evaporator. A validation study was performed to compare the Sockeye code against data produced at the Ulsan National Institute of Standards and Technology for a slightly overfilled sodium heat pipe. The heat pipe was operated under natural convection cooling and a decreasing inactive length was observed. The experiment was modeled in Sockeye using the condenser pool model, the front non-condensable gas (NCG) model, and the mixture NCG model separately to reproduce the inactive length phenomenon. Also, a new capability was implemented in MOOSE to allow for conjugate heat transfer from the condenser based on the Churchill-Chu correlation for natural convection. The condenser pool model showed an increasing inactive length, demonstrating that the behavior was likely not caused by a pool of excess liquid. The front NCG model was able to show good agreement with the experiment, but it suffered convergence issues at the front. Finally, the mixture NCG model gave good results when the axial mesh was sufficiently refined. This work culminated in additions to the Sockeye documentation and a contribution to a journal article that will be published later. This poster is a summary of my work which I can take back to my university for presentation.

42 - ENGINEERING↗
DOE OSTI2026

Validation of OpenPronghorn for Periodic Hill Flow Separation

OpenPronghorn is an open-source, MOOSE-based thermal-hydraulics simulation tool used for advanced reactor analysis. As an open-source code, it offers a transparent framework for validating governing equations, assumptions, and numerical methods against established benchmarks. This study evaluates OpenPronghorn's RANS turbulence model against the ERCOFTAC Case 81 periodic hill benchmark, a standard test case for separated turbulent flow featuring curved-wall separation, recirculation, shear-layer development, and reattachment. Streamwise velocity profiles predicted by OpenPronghorn were compared to reference LES data at multiple x/h locations. Results show that OpenPronghorn captures the overall trend of the velocity profiles, but the largest discrepancies occur in the separated-flow region, where turbulence is highly anisotropic and strongly affected by adverse pressure gradients and wall curvature—conditions that are inherently difficult for standard RANS models to resolve. Future work will test alternative k-e model variants and correction terms to improve prediction accuracy in this region.

42 - ENGINEERING↗
DOE OSTI2020

Grassy ELM regime at low pedestal collisionality in high-power tokamak plasma

Spontaneous mitigation of giant ELMs and appearance of grassy ELMs have been observed repeatedly at low pedestal collisionality (${\nu }_{\mathrm{e},\mathrm{p}\mathrm{e}\mathrm{d}}^{{\ast}}\sim 0.15$) in the high-power (P inj > 13 MW) hybrid scenario in DIII-D. Higher β p and higher q 95 appear to be beneficial to achieving the grassy ELM regime. The grassy ELM H-mode plasma shows high energy confinement performance (H 98y2 up to over 1.5) under the conditions of high neutral beam torque and high core rotation. The pedestal width appears to exceed the EPED1.0 model prediction by more than 50%. Pedestal stability analysis performed with the ELITE code indicates that the stability against low-n kink/peeling modes is improved with increased auxiliary heating power and the operational point in the grassy ELM regime is located near the ballooning boundary. The pedestal stability characteristics during the grassy ELM crashes have been investigated in comparison with the giant ELM crashes based on plasma profiles experimentally measured with high resolution and accuracy. It has been found that the underlying mechanism for the observed small-amplitude ELM crashes is mainly the expansion of the ballooning stability boundary induced by an initial radially localized collapse in the pedestal, which helps to stop the growth of instabilities and further collapse of the pedestal. The effect of electron density pedestal on mitigating edge localized instabilities has been analyzed by numerical simulation, suggesting that the electron density pedestal characterized by high n e,sep /n e,ped and low density gradient helps to stabilize peeling-ballooning modes because of a low pressure pedestal gradient and to lower the ballooning boundary mainly because of a low ion diamagnetic frequency in the pedestal region, thus triggering ballooning instabilities and producing the intrinsic grassy ELMs. Numerical simulation of the Chinese fusion engineering test reactor (CFETR) with the SOLPS code indicates that the separatrix density might be insensitive to the electron diffusivity in the pedestal region and increase with the power flowing from the core region to the edge region. Furthermore, pedestal stability analysis suggests that the flat density pedestal with high separatrix density obtained in the high-power plasma in CFETR would make the operational point close to the ballooning boundary, which is considered to help destabilize ballooning instabilities and facilitate the access to the grassy ELM regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Wind and Weather Variability within the Californian Offshore Wind Energy Areas

Weather variability over the Northeast Pacific (NEP) region and its influence on wind resources within the Californian offshore wind energy areas (WEAs) at Humboldt and Morro Bay are characterized using 20-years reanalysis model and satellite data. The hub-height (180 m) winds at both locations are predominantly northwesterly driven by the NEP high pressure system, with strong coastal gradients in surface pressure, fluxes, planetary boundary layer (PBL) depths and cloudiness. These sharp coastal gradients and strong annual cycles of temperature and moisture advections pose potential challenges in accurately modeling the local wind resource. Hub-height wind speeds and power capacity factors significantly vary for different regimes of PBL depths, surface fluxes and rain area fractions. This highlights the importance of studying the physical mechanisms driving these weather regimes, hence our analysis of how large-scale NEP weather variability drives the local meteorology at the WEAs. Furthermore, at both WEAs, PBL tops and cloud boundaries intersect the rotor layer (80-280 m) more than 30% and 20% of the time, respectively. While PBL depths significantly modulates hub-height winds and power, cloud boundaries do not have a similar impact, likely due to reanalysis errors in simulating cloud boundaries accurately. These findings underscore the challenges in deploying tall wind turbines in shallow cloudy boundary layers, where the interaction between clouds, precipitation, and atmospheric layers can impact turbine efficiency. As turbines grow taller and are deployed in more complex meteorological conditions, understanding these interactions is crucial for improving wind power forecasting and optimizing energy production in coastal regions.

17 WIND ENERGY↗
DOE OSTI2021

Mitigating Injection‐Induced Seismicity Along Basement Faults by Extraction: Application to 2016–2018 Pohang Earthquakes

Abstract The risk of induced earthquake is considered as the major obstacle for sustainable subsurface energy activities. As a mitigation approach, mechanistic risk assessment can be essential to minimize the potential of seismic hazard. Coupled hydro‐mechanical modeling approach is integrated with geomechanical analysis to better assess the occurrence of seismicity induced by injection and/or extraction, which will provide a basis to update operational conditions (e.g., number and location of mitigating wells, rate and duration of operations) for mitigating earthquake nucleation. The Coulomb stability analysis and seismicity rate estimates from this generic study indicate that simultaneous extraction can relieve injection‐induced pore‐pressure buildup, but also transfer and accumulate elastic energy along the faults. Field‐scale simulation with variation in well design for the Pohang EGS site, South Korea, show that multiple wells aligned across a preexisting fault form gradients of pressure and stress fields normal to the fault, which can enhance the direct impact of diffusion process and poroelastic stressing on the fault stability perturbation. Our mechanistic study suggests multiphysics‐oriented mitigation strategies that consider (1) possible mechanical deformation, caused by hydro‐mechanical coupling effects, and (2) proper configuration of well operations based on site‐specific geological characteristics.

Chang, K. W.↗
DOE OSTI2024

Numerical modeling of plasma assisted deflagration to detonation transition in a microscale channel

Here, this work numerically studies the plasma assisted deflagration to detonation transition (DDT) of H 2 /O 2 mixtures in a microscale channel with detailed chemistry and transport. The results show that the DDT onset time is non-monotonically dependent on the discharge pulse number. The DDT is accelerated with small pulse numbers, whereas retarded with large ones. Two different DDT regimes, respectively at a small and large plasma discharge number, via acoustic choking of the burned gas and plasma-enhanced reactivity gradient without acoustic choking, are observed. Without plasma discharge, pronounced pressure and temperature gradients in front of the flame are generated by acoustic compression after the choking of the burned gas, triggering DDT via autoignition. With small plasma pulse numbers, the plasma-generated species enhance the ignition kinetics and lead to an increased reactivity in the boundary layer. After the choking of the burned gas, the plasma-enhanced reactivity advances the sequence of autoignition near the wall, strengthens ignition-shock wave coupling, and accelerates DDT. However, with a large discharge pulse number, a direct autoignition initiating DDT can occur without the acoustic choking of the burned gas due to the strongly accelerated reactivity and elevated temperature. In this case, DDT onset is retarded because the elevated temperature increases sonic velocity and the increased reactivity accelerates fuel oxidation in front of the flame, decelerating the formation of a leading shock and subsequent pressure buildup ahead of the flame. The present modeling reveals that no matter with or without plasma discharge, DDT is initiated by autoignition in thermal, pressure, and reactivity gradient fields via the Zel'dovich gradient mechanism. The acoustic choking of the burned gas may not be the necessary condition of DDT with strong plasma-enhanced reactivity gradient. This work provides an answer to the experimentally observed non-monotonic DDT onset time by plasma, which provides guidance to control DDT in advanced detonation engines and fire safety of hydrogen-fueled catalytic reactors in microchannels by non-equilibrium plasma discharge.

33 ADVANCED PROPULSION SYSTEMS↗