Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “edge effects”

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 1,333 records · Page 74
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Predicting nutrition and environmental factors associated with female reproductive disorders using a knowledge graph and random forests

Female reproductive disorders (FRDs) are common health conditions that may present with significant symptoms. Diet and environment are potential areas for FRD interventions. We utilized a knowledge graph (KG) method to predict factors associated with common FRDs (for example, endometriosis, ovarian cyst, and uterine fibroids). We harmonized survey data from the Personalized Environment and Genes Study (PEGS) on internal and external environmental exposures and health conditions with biomedical ontology content. We merged the harmonized data and ontologies with supplemental nutrient and agricultural chemical data to create a KG. We analyzed the KG by embedding edges and applying a random forest for edge prediction to identify variables potentially associated with FRDs. We also conducted logistic regression analysis for comparison. Across 9765 PEGS respondents, the KG analysis resulted in 8535 significant or suggestive predicted links between FRDs and chemicals, phenotypes, and diseases. Amongst these links, 32 were exact matches when compared with the logistic regression results, including comorbidities, medications, foods, and occupational exposures. Mechanistic underpinnings of predicted links documented in the literature may support some of our findings. Our KG methods are useful for predicting possible associations in large, survey-based datasets with added information on directionality and magnitude of effect from logistic regression. These results should not be construed as causal but can support hypothesis generation. This investigation enabled the generation of hypotheses on a variety of potential links between FRDs and exposures. Future investigations should prospectively evaluate the variables hypothesized to impact FRDs.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Navigating Large Chemical Spaces Using Graph Theory and Integer Programming

Navigating and analyzing large chemical spaces are necessary to accelerate the design and discovery of new molecules and chemical processes. In this work, we introduce a computational framework that integrates graph theory and integer programming to enable the efficient navigation of large chemical spaces. Our framework represents the chemical space as a graph, wherein nodes represent molecules and edges represent the degree of similarity or connectivity based on domain-specific information. Using the graph representation, we identify representative molecules by computing the so-called minimum dominating set (MDS), which in our context is the minimum set of molecules that is connected to all other molecules. We present a suite of solution strategies for the MDS problem including heuristic and rigorous integer programming (IP) approaches. We show that these approaches allow us to capture physicochemical properties and domain-specific logic and constraints, facilitating the identification of molecules with the target properties. We demonstrate the effectiveness of the proposed approach by navigating the chemical space of per- and polyfluoroalkyl substances (PFAS); this comprises approximately 15,000 molecular structures. We compare our framework against traditional dimensionality reduction and clustering methods such as t-SNE and K-means clustering.

Chemical structure↗
DOE OSTI2021

All-Solid-State Lithium–Sulfur Batteries Enhanced by Redox Mediators

We report redox mediators (RMs) play a vital role in some liquid electrolyte-based electrochemical energy storage systems. However, the concept of redox mediator in solid-state batteries remains unexplored. Here, we selected a group of RM candidates and investigated their behaviors and roles in all-solid-state lithium-sulfur batteries (ASSLSBs). The soluble-type quinone-based RM (AQT) shows the most favorable redox potential and the best redox reversibility that functions well for lithium sulfide (Li 2 S) oxidation in solid polymer electrolytes. Accordingly, Li 2 S cathodes with AQT RMs present a significantly reduced energy barrier (average oxidation potential of 2.4 V) during initial charging at 0.1 C at 60 °C and the following discharge capacity of 1133 mAh g s -1 . Using operando sulfur K-edge X-ray absorption spectroscopy, we directly tracked the sulfur speciation in ASSLSBs and proved that the solid-polysulfide-solid reaction of Li 2 S cathodes with RMs facilitated Li 2 S oxidation. In contrast, for bare Li 2 S cathodes, the solid-solid Li 2 S-sulfur direct conversion in the first charge cycle results in a high energy barrier for activation (charge to ~4 V) and low sulfur utilization. The Li 2 S@AQT cell demonstrates superior cycling stability (average Coulombic efficiency 98.9% for 150 cycles) and rate capability owing to the effective AQT-enhanced Li-S reaction kinetics. This work reveals the evolution of sulfur species in ASSLSBs and realizes the fast Li-S reaction kinetics by designing an effective sulfur speciation pathway

25 ENERGY STORAGE↗
DOE OSTI2025

Au 76 (SC 6 H 4 - p -CH 3 ) 42 Square Quantum Platelet: One-Dimensional Growth of Quantum Rods Turns 90 Degrees

The growth pattern of atomically precise nanoclusters (NCs) is of fundamental interest, and the structural effect on their photoluminescence (PL) is important owing to their PL in the second near-infrared (NIR-II) range (900–1700 nm optoelectronically or 1000–1700 nm biologically) that holds great promise for optoelectronic and biomedical applications. Because of the small energy gaps required for NIR-II emission, the PL performance of NIR-II luminophores is largely limited by nonradiative processes. In this work, we discovered a Au 76 (p-MBT) 42 (Au 76 ) (p-MBTH = p-methylbenzenethiol) nanocluster featuring a face-centered cubic (fcc) core in a square shape (edge length: 1 nm). This square quantum platelet can be viewed as a side-facet (010) growth of the Au 52 (p-MBT) 32 (Au 52 ) rod, as opposed to the (001) facet growth. We found that Au 76 exhibits bright emission centered at 970 nm with a PL quantum yield (PLQY) of 30% in solution under ambient conditions, which can be further enhanced to 40% when the solution is deaerated. X-ray crystallography analysis coupled with time-resolved spectroscopy revealed that the nearly doubled PLQY compared to Au 52 (18.3%) was resulted from shorter Au–Au bond lengths in Au 76 (average 2.835 Å) than that in Au 52 (3.04 Å). This work provides important insights into the design of highly luminescent NCs, which are promising for photovoltaics, photocatalysis, and optoelectronic applications.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Impact of atomic substitution on core-hole relaxation dynamics: A study of Br 2 and IBr

Understanding inner-shell decay processes in heavy-element molecules is essential for unraveling x-ray-induced photodynamics and advancing molecular imaging techniques. Here, in this study, we investigate the influence of atomic substitution on core-hole relaxation dynamics and molecular fragmentation in Br 2 and IBr, initiated by x-ray ionization absorption at the Br K-edge. Using a combination of x-ray/ion coincidence measurements and Monte Carlo/molecular dynamics simulations, we track the charge distribution and the kinetic energy release (KER) of fragment ions with a total charge from 2+ to 8+. For both molecules, the simulated KER values show good agreement with experiment across different fragmentation channels. Our comparison reveals that substituting Br with the heavier I atom in IBr has a minimal impact on the inner-shell electronic decay process but significantly influences nuclear motion, leading to slower dissociation and thus a KER close to the Coulomb limit—an effect attributed to the atomic mass. These findings highlight the interplay between electronic and nuclear effects in molecular fragmentation, particularly in heavy-element species, and provide new insights into medical therapies, structural biology, and astrophysics.

Bhat, Nivedita [Argonne National Laboratory (ANL),↗
DOE OSTI2023

Feasibility of raised inner strike point equilibria scenario in ITER for detritiation from beryllium co-deposits

Abstract In ITER, tritium retention primarily occurs through co-deposition with beryllium. To avoid exceeding the strict tritium inventory limit, efficient tritium recovery techniques are essential. Baking is the ITER baseline for tritium recovery, but its effectiveness in removing tritium from thick beryllium layers is limited. A raised strike point scenario is considered an alternative method for removing tritium from the ITER inner vertical divertor target by heating components via plasma flux. This paper presents SOLPS-ITER code simulations conducted under various conditions, assessing the divertor performance and tritium outgassing of the raised strike point scenario. As the strike point is raised, recycled neutrals are not efficiently baffled by the dome and scrape-off layer, significantly changing the neutral trajectory and ionization source distribution. This improves detachment accessibility but worsens core-edge compatibility compared to the baseline scenario. However, in the partially detached condition, the impact of raising the strike point, perpendicular transport, and q 95 on target heat flux is not significant, as it primarily scales with the input power. Target heat flux is translated to target surface temperature using a simplified heat transfer model that considers the 3D target monoblock geometry and active cooling condition, excluding Be layer thermal properties. For partially detached divertor conditions, the bulk tungsten monoblock surface temperature remains below the baking temperature, which is insufficient for efficient tritium outgassing under the actively cooled ITER divertor condition. However, considering the potential thermal contact resistance between the beryllium and tungsten layers, which may significantly impact temperature distribution, the temperature of the beryllium layer can be raised to a level sufficient for efficient tritium outgassing. Therefore, the raised strike point scenario can be considered as an alternative in-vessel tritium removal technique.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A Hybrid Data-Driven and Model-Based Anomaly Detection Scheme for DER Operation: Preprint

This paper proposes a hybrid data and model-based anomaly detection for securing the operation of distributed energy resources (DERs) in distribution grids. Data-driven autoencoders (AE) are set up at the edge level by taking local DER data and detect anomalous operations by leveraging the reconstruction ability. In parallel, model-based state estimation (SE) is running at the system level by taking system models and measurements, the anomalies are identified by analyzing the measurements residual. The hybrid scheme preserves the benefits of both data-driven and model-based analysis and thus improves the robustness and accuracy of anomaly detection. It can be established by getting full use of the existing infrastructures in distribution grids. Numerical tests on a realistic distribution feeder in Southern California highlight the effectiveness as well as benefits of the proposed scheme.

anomaly detection↗
DOE OSTI2026

Analysis of RF Sheath-Driven Tungsten Erosion at RF Antenna in the WEST Tokamak

This study applies the newly developed STRIPE (Simulated Transport of RF Impurity Production and Emission) framework to analyze tungsten (W) erosion at RF antenna structures in the WEST tokamak. STRIPE integrates SolEdge3x for edge plasma backgrounds, COMSOL for 3D RF sheath potentials, RustBCA for sputtering yields, and GITR for impurity transport and ion energy–angle distributions. Building on prior work by Kumar et al. (2025) Nuclear Fusion, 65, 076039, which validated STRIPE for WEST ICRH discharge #57877, the present study provides a spatially resolved assessment of gross W erosion at both Q2 antenna limiters under ohmic and ICRH conditions. Simulations using 2D SolEdge3x profiles in COMSOL capture rectified sheath potentials exceeding 300 V, leading to strong upper-limiter localization. Both poloidal and toroidal asymmetries are observed and attributed to RF sheath effects, with modeled erosion patterns deviating from experiment—highlighting sensitivity to sheath geometry and plasma resolution. Erosion is driven primarily by high-charge-state oxygen ions (O6+–O8+), while D+ plays a negligible role. Assuming a plasma composition of 1% oxygen and 98% deuterium, STRIPE predicts a 30-fold increase in gross W erosion from ohmic to ICRH phases, consistent with a >25-fold rise in W-I (400.9 nm) brightness. Quantitative agreement is within 5% in the ohmic phase and 30% under ICRH, demonstrating predictive capability. Importantly, the study shows that the magnitude of ICRH-driven W erosion depends strongly on the concentration of light impurities (O, B, N, C), which drive sputtering through high charge states. Cleaner plasma conditions with reduced impurity content are therefore expected to substantially mitigate antenna W sources in WEST and other toroidal fusion devices. These findings establish STRIPE as a predictive framework for RF-induced plasma–material interactions and support its application to reactor-scale antenna design.

Kumar, Atul [ORNL] (ORCID:0000000277475210)↗
DOE OSTI2021

Efficient Sampling of Complex Interdependent and Multiplex Networks

Efficient sampling of interdependent and multiplex infrastructure networks is critical for effectively applying failure and recovery algorithms in real-world settings, as well as to generate property-preserving reduced-order graph-based ensembles that address topological uncertainties. In this paper, we first explore the performance, i.e. the success in preserving graph properties, of graph sampling algorithms for interdependent and multiplex networks with synthetic and real-world graphs. We simulate sampling algorithms under different parameter settings. These settings include probabilistic graph generators, coupling patterns, and various performance metrics. Our results show that while Random Node and Random Walk sampling algorithms perform best for interdependent networks, Random Edge and Forest Fire sampling algorithms perform best for multiplex networks. Second, we propose and implement a novel similarity-based sampling algorithm for multiplex networks that samples only log(N) number of layers of an N-layer multiplex network while yielding computational savings with performance guarantees. Experimental results show that similarity sampling outperforms complete sampling of all layers while decreasing performance costs from a linear scale to a logarithmic one. Our results also indicate that similarity-based sampling outperforms complete sampling and random selection in nearly all scenarios when tested with real-world data.

Subasi, Omer↗
DOE OSTI2021

Interpretation of hydrothermal conditions, production-injection induced effects, and evidence for enhanced geothermal system-type heat exchange in response to >30 years of production at Roosevelt Hot Springs, Utah, USA

The Roosevelt Hot Springs hydrothermal system is located at the base of the Mineral Mountains in southwestern Utah on the eastern side of the Basin and Range. Hydrothermal activity is related to relatively recent bimodal magmatism, and the system is hosted in coarsely crystalline rock made of Oligocene–Miocene granitoids and Precambrian gneiss. The hydrothermal plume covers ~5 km 2 , with a maximum temperature of 268 °C at ~750 m depth, and a vertically extensive fault-fracture mesh east of the Opal Mound fault controls the upflow of hydrothermal fluids. Power generation (currently 38 MWe gross) began in 1984, and up through 2016, four wells were used for fluid production, and three wells were used for edge-field injection. Chemical analyses of produced fluids show that modern reservoir fluid compositions are similar to but more concentrated than those at the start of production, having near-neutral pH, total dissolved solids of 7000–10,000 mg/kg, and ionic ratios of Cl/HCO 3 ~50–100, Cl/SO4 ~50–100, and Na/K ~4–5. Chemical geothermometers indicate equilibration temperatures that mainly range between 240° and 300 °C. Early production induced a steep drop in pressure (~3.0–3.5 MPa), which was accompanied by a 250–300 m lowering of piezometric levels in wells and development of a shallow steam zone across the system. Hydrothermal fluid compositions evolved continuously in response to production-related steam-loss and injection breakthrough, which is reflected by gradual increases in chloride of up to 35% and stable isotope ratios of up to ~2‰δ 18 O and ~10‰δD. Simple mixing model calculations suggest that there has been a significant amount, ~10–20 MW th , of sustained multi-decadal heat mining and enhanced geothermal system (EGS)–type heat transfer by the injectate as it returns to the production zone. Overall, the two factors that have sustained long-term power production (currently 38 MWe gross) are the increased upflow of deep chloride water and, to a lesser extent, the mining of heat at <1 km depth.

58 GEOSCIENCES↗
DOE OSTI2025

Recrystallization, cracking, and erosion of dispersoid-strengthened tungsten materials during exposure to divertor plasmas

In this study, we investigated the effects of combined intense particle and heat flux exposure on advanced tungsten plasma-facing materials within the DIII-D fusion facility. Our test matrix included two types of dispersoid-strengthened tungsten (containing either 100 nm diameter TiO 2 or Ni particles), along with high-purity polycrystalline tungsten as a reference. This experiment relied on a sample geometry angled at 15° relative to the divertor surface, thereby allowing the surfaces to intercept steady-state perpendicular heat fluxes (q ⟂ ) ranging from 10.1 to 19.6 MW/m 2 . During each shot, the samples were exposed to 42 Hz edge-localized modes (ELMs), allowing us to test the material response to transient heating. We correlated the exposure conditions with extensive post-test surface composition analysis and microscopy to determine how the plasma modified each surface. The angled specimens closest to the strike point received the highest combined heat and particle flux and melted midway through the experiment. EBSD analysis revealed they were completely recrystallized throughout, with an average grain size >100 µm. On the other hand, the specimens that received a lower steady state heat flux survived with more superficial surface damage. Whereas the high-purity polycrystalline tungsten exhibited a higher surface roughness, the dispersoid-strengthened material exhibited more extensive shallow inter-granular cracking. In addition, the surface was depleted of dispersoids following plasma exposure, possibly because of evaporation and/or sputtering. The results described here provide insights into the performance of these materials in a fusion environment which can guide further optimization for use in long-pulse devices.

Kolasinski, Robert D. [Sandia National Laboratorie↗
DOE OSTI2021

Unraveling the Formation Mechanism of a Hybrid Zr-Based Chemical Conversion Coating with Organic and Copper Compounds for Corrosion Inhibition

Environmental-friendly chromate-free, zirconium (Zr)-based conversion coating is a promising green technology for corrosion protection. Additives in the surface treatment provides critical functionalities and performance improvements; however, mechanistic understanding as of how the additives influence the coatings remains unclear. In this study, a new organic-inorganic hybrid Zrbased conversion coating which combines copper (Cu) compounds and polyamidoamine (PAMAM), taking advantages of a complementary nature of organic and inorganic additives. Here, a multimodal approach combining electron and X-ray characterization is applied to study the interaction of Cu 2+ and PAMAM and the resulting impacts on the coating formation. Adding PAMAM changed the surface morphology, thickness, distribution of Cu in the cluster and void formation of the coatings. High PAMAM (100-200 pm) leads to little conversion coating formation, and low PAMAM (0-25 ppm) shows voids formation under the coatings. Moreover, PAMAM incorporates in the coating in a form of PAMAM-Cu complex with a higher concentration towards the surface, provides an organic layer at the surface of the coating. X-ray absorption near-edge structure (XANES) spectroscopy shows a difference between the conventional and hybrid coating treatments in an alkaline solution to simulate the E-coat process, suggesting the contribution of PAMAM in the enhanced chemical stability in an alkaline environment. Therefore, an intermediate range of addition of PAMAM (50 ppm) is optimal to 1) avoid excessive voids formation, 2) promote some Cu cluster formation and thus enhancing the Zr-based coating formation and 3) incorporate organic components into the coating to improve the adhesion of the subsequent coatings. Overall, this work furthers our knowledge on the formation mechanism of an effective and environmentally friendly hybrid conversion coating for corrosion inhibition, demonstrating a critical processing-structureproperty relationship. This study will benefit future development of green and effective surface treatment technology.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2026

Interpreting AI for fusion: An application to plasma profile analysis for tearing mode stability

Artificial intelligence models have demonstrated strong predictive capabilities for various instabilities in fusion devices such as Tokamaks, including tearing modes (TM), edge localized modes, and disruptive events, but their opaque nature raises concerns about safety and trustworthiness when applied to fusion power plants. Here, we present a physics-based interpretation framework using a TM prediction model as a demonstration that is validated through a dedicated DIII-D TM avoidance experiment. By applying Shapley analysis, we identify how profiles such as rotation, temperature, and density contribute to the model's prediction of TM stability. Our analysis shows that in our experimental scenario, core electron temperature and rotation peaking play the primary role in TM stability, while density changes have smaller effects on stability. We show that off-axis ion temperature stabilizes TMs, suggesting that off-axis neutral beam heating can further stabilize this scenario. This work presents a generalizable ML-based event prediction methodology, from training to physics-driven interpretation, bridging the gap between physics understanding and opaque ML models.

Farre-Kaga, Hiro J. [Princeton Univ., NJ (United S↗
DOE OSTI2022

A Hybrid Data-Driven and Model-Based Anomaly Detection Scheme for DER Operation

This paper proposes a hybrid data and model-based anomaly detection scheme to secure the operation of distributed energy resources (DERs) in distribution grids. Data-driven autoencoders are set up at the edge device level and they use local DER operational data as inputs. The abnormal statuses are detected by analyzing reconstruction errors. In parallel, modelbased state estimation (SE) is set up at the central level and it uses system-wide models and measurements as data inputs. The anomalies are identified by analyzing measurement residuals. The hybrid scheme preserves the benefits of both data-driven and model-based analyses and thus improves the robustness and the accuracy of anomaly detection. Numerical tests based on the model of a real distribution feeder in Southern California highlight the proposed scheme's effectiveness and benefits.

anomaly detection↗
DOE OSTI2022

Toward holistic understanding of the ITER-like resonant magnetic perturbation (RMP) ELM control on KSTAR

Here, KSTAR has clarified a set of unresolved 3D physics issues utilizing the ITER-like in-vessel, three-row, resonant magnetic perturbation (RMP) configurations. Since RMP-driven, edge-localized-modes (ELMs)-crash control elevates the divertor heat flux peak through its impact on edge plasma parameters and transport, a series of intentionally misaligned RMP configurations (IMCs) have been explored to investigate the relationship between RMP ELM control and divertor heat fluxes, while searching for an ideal IMC that could be favorable in both aspects. First of all, the contrasting influence of kink vs anti-kink phasing on the ELM-crash suppression has been articulated, demonstrating the synergistic benefit of 'kink' phasing on ELM-crash-suppression. On the other hand, the three-row IMC in the anti-kink phasing becomes more insensitive to the ELM-crashes at the sub-marginal level of RMP, consistent with theory. Meanwhile, the divertor 'wetted' area of ELM-crash-suppression gets narrower than that of ELM-crash-mitigation, suggesting that ELM-crash-mitigation remains advantageous over ELM-crash-suppression in terms of time-averaged divertor thermal loading. In comparison, based on a set of two-row IMCs, no evidence of divertor heat flux broadening was found during ELM-crash-suppression, supporting a hypothesis that the dispersal of the divertor heat flux in three-row IMCs cannot be driven by helically structured two-row RMPs alone. Among ITER-like three-rows, lower two-row RMPs have been found to be much more effective in suppressing the ELM-crashes than upper two-row RMPs. Although it is quite preliminary, the up/down asymmetric dependence of RMP coupling may be generically attributed to lower-single-null plasmas. Such a holistic understanding of RMP-driven, ELM-crash-control in KSTAR is expected not only to elucidate various subtle points in the vicinity of ELM-crash-suppression, but also to clarify the relevant divertor thermal loading issues for ITER and beyond.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Dynamics of quantum Hall interfaces

A quantum Hall (QH) interface is different from an ordinary QH edge, as the latter has its location determined by the confining potential, while the former can be unpinned and behave like a free string. In this paper, we demonstrate this difference by studying three different interfaces formed by (i) the Laughlin state and the vacuum, (ii) the Pfaffian state and the vacuum, and (iii) the Pfaffian and the anti-Pfaffian states. We find that stringlike interfaces propagating freely in the QH system lead to very different dynamical properties from edges. This qualitative difference gives rise to fascinating physics and suggests a different direction for future research on QH physics. We also discuss briefly possible analogies between QH interfaces and concepts in string theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

The effect of copper additives on hybrid Zr-based chemical conversion coating morphology and chemical compositions

In the realm of corrosion protection, Zr-based conversion coatings offer an environmentally friendly, chromate-free alternative to conventional coating. This study uses advanced X-ray, electron microscopy and electrochemical testing techniques to better understand the impact of varying Cu²⁺ ion concentrations on the characteristics of Zr-based coatings on Fe substrates. Our findings demonstrate that within the tested conditions, higher Cu²⁺ additive concentrations at 40 ppm enhance surface characteristics, increasing stability towards anti-corrosion capability, particularly under NaCl treatment. Conversely, at a lower Cu²⁺ concentration of 20 ppm, coatings exhibited more significant dissolution of Cu clusters and increased vulnerability to chloride-induced degradation. X-ray photoelectron spectroscopy and synchrotron X-ray fluorescence (XRF) and X-ray absorption near edge structure (XANES) spectroscopy analyses revealed Cu²⁺ formation in Cu20, while Cu40 retained its metallic state (Cu⁰) with a slower reaction rate. Cu20 offers some protection but lacks durability, whereas in Cu40, lower I corr , and enhance corrosion resistance, making it ideal for protective coatings in salt-water environments. This study underscores the importance of balancing Cu²⁺ ions concentration in the coating solution to optimize performance, highlighting the role of Cu in enhancing both surface properties and long-term stability.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Modeling of the effects of non-equilibrium excitation and electrode geometry on H 2 /air ignition in a nanosecond plasma discharge

In this work, we present the results of two-dimensional modeling of the effects of non-equilibrium excitation and electrode geometry on H 2 /air ignition in a nanosecond plasma discharge. A multiscale adaptive reduced chemistry solver for plasma assisted combustion (MARCS-PAC) based on PASSKEy discharge modeling package and compressible multi-component reactive flow solver ASURF+ is developed and validated. This model is applied to simulate the impact of non-equilibrium plasma excitation and electrode geometry and heat loss on the dynamics of the discharge from streamer to spark and ignition kernel development in a H 2 /air mixture with a pair of cylindrical electrodes. The results show that the plasmagenerated species (N 2 (A), N 2 (B), N 2 (a'), N 2 (C), O( 1 D), O and H) in the spark and afterglow significantly accelerate the ignition kernel development. The increase of discharge voltage at the same total discharge energy promotes the non-equilibrium active species production. It is found that the production of electronically excited species at higher reduced electric field strength is more efficient in enhancing ignition in comparison to the vibrational excitation and heating. Moreover, the 2D simulation clearly reveals that the electric field and active species distribution are highly non-uniform. The streamers are initiated at the sharp outer edges of the negative and positive electrodes by a strong electric field while the electric field is much weaker at the centerline of the electrodes. Furthermore, the simulations reveal that the ignition enhancement is sensitive to the variation of electrode shape, diameter, and gap size due to the changes of electric field distribution and location of streamer formation. A cylindrical electrode produces a larger discharge volume and ignition kernel than the parabolic and spherical electrodes, when the discharge is localized near the axis of the gap. It is found that there is a non-monotonic dependence of ignition kernel size on the electrode diameter and inter-electrode distance. The increase of electrode diameter and gap size above the optimal conditions leads to the reduction of ignition kernel volume, due to the decrease of active species concentration and gas temperature. At a larger electrode surface area and electrode diameter as well as smaller electrode gap size, the heat loss to electrode plays a greater role in reducing the ignition kernel size and slowing ignition kernel development. This work provides insights and guidance to understand the kinetic enhancement of non-equilibrium plasma and the effects of electrode geometries on ignition for the optimization ignitors in advanced engines.

42 ENGINEERING↗