Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase space methods”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

The crystal structure of Fe2S at 90 GPa based on single-crystal X-ray diffraction techniques

The Fe-S system was explored in a laser-heated diamond-anvil cell at 89(2) GPa and 2380(120) K to better understand the phase stability of Fe 2 S. Upon temperature quenching, crystallites of Fe 2 S were identified, and their structure was investigated using single-crystal X-ray diffraction techniques. At these conditions, Fe 2 S adopts the C23 structure (anti-PbCl 2 , Co 2 P) with space group Pnma (Z = 4). This structure consists of columns of corner-sharing, FeS 4 tetrahedra, and columns of edge-sharing FeS 5 square pyramids linked along edges in the b direction. Sulfur is in ninefold coordination with Fe. This work marks the first high-pressure structural solution and refinement of Fe 2 S synthesized in a multigrain Fe+FeS sample at 90 GPa and 2400 K and establishes the stability of C23 Fe 2 S at these conditions. A previous powder diffraction study reports an orthorhombic Fe 2 S phase with a C37, Co 2 Si-like unit cell above 190 GPa. A C23–C37 structural transition is inferred to explain the previously observed unit-cell parameters at higher pressures and temperatures. These results highlight the utility of applying single-crystal X-ray diffraction techniques to high P-T multigrain samples to explore the structural properties of iron-rich phases in Earth and planetary cores.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Performance of Windows in Walls With Continuous Insulation

Window openings in walls are a significant contributor to poor thermal performance because of thermal bridging through the framing members (e.g., studs, joists, plates, bracing) and because windows lack the thermal properties of insulation. Window installation guidance for walls with continuous insulation (CI) is critical for continued market growth of this energy efficiency technology. This research project offers window manufacturers a starting point and a potential path toward developing installation instructions for windows over CI. The objectives of the research include evaluating the common method for installing windows in walls with CI, as well as establishing acceptance criteria for evaluating the performance of windows installed in walls with and without CI. The research measures: 1. The performance characteristics (e.g., water management, structural integrity) of windows in walls without CI. 2. The performance of different thicknesses and types of CI used in walls. 3. The performance of different types of window assemblies (e.g., double-hung windows, mulled double-hung windows, mulled casement windows, and slider windows) installed over CI. 4. The performance of window flange types (e.g., rigid mounting and less robust flanges) installed over CI. 5. Installing windows over CI using baseline installation instructions versus window manufacturer installation instructions. The project’s sequential testing protocol consists of the following: • A water penetration resistance testing adapted from two ASTM standards: E331 (uniform static air pressure in four steps) and E547 (cyclic static air pressure) • A temperature cycling adapted from ASTM E2264 Method B (convective hot air) • A service condition wind loading test adapted from ASTM E330 • A six-month vertical displacement observation phase prior to the structural performance testing • A final water penetration resistance test after vertical displacement observation • A structural performance test adapted from ASTM E330. Key research findings include: • The criterion for passing a water penetration resistance test is that there is no water overflowing at the interior face of the studs. If there is any bubbling or slight pooling of water at the sill, then it must recede after the pressure is removed. Excessive leakage and/or water leaking to the interior face of the framing around the window constitutes a failure. • All single double-hung windows installed directly to lumber or over oriented strand board passed all test protocols. • For most wall specimens, the test results showed that the use of foam sheathing did not affect the performance of the window for water leakage. • All wall specimens underwent temperature cycling. The results indicated that temperature cycling had little to no effect on windows installed over foam sheathing. • For wall specimens that underwent six-month vertical displacement monitoring, the results showed that windows installed over foam sheathing do not sag over time. • The single-hung and double-hung windows installed using window manufacturer installation instructions passed the structural performance test, compared to failures observed in windows that were installed using generic installation instructions. The generic and manufacturer installation methods differed on the following construction details: type of fasteners, fastening patterns on the flanges, and window shimming details. • Additional testing is required to determine methods to improve structural pressure performance of slider windows. Potential solutions that would require additional testing may include fastener spacing, different types of fasteners, masonry window clips, construction adhesive, foam sealant, stronger window flange material, and/or straps.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Analytical gradient-based optimization of CALPHAD model parameters

The calibration of CALPHAD (CALculation of PHAse Diagrams) models involves the solution of a very challenging high-dimensional multiobjective optimization problem. Traditional approaches to parameter fitting predominantly rely on gradient-free methods, which while robust, are computationally inefficient and often scale poorly with model complexity. In this work, we introduce and demonstrate a generalizable framework for analytic gradient-based optimization of the parameters of the CALPHAD model enabled by the recently formalized Jansson derivative technique. This method allows for efficient evaluation of gradients of thermodynamic properties at equilibrium with respect to model parameters, even in the presence of arbitrarily complex internal degrees of freedom. Leveraging these semi-analytic gradients, we employ the conjugate gradient (CG) method to optimize thermodynamic model parameters for four binary alloy systems: Cu-Mg, Fe-Ni, Cr-Ni, and Cr-Fe. Across all systems, CG achieves comparable or superior optimality relative to Bayesian ensemble Markov Chain Monte Carlo (MCMC) with improvements in computational efficiency ranging from one to three orders of magnitude. Furthermore, our results establish a new paradigm for CALPHAD assessments in which high fidelity data-rich model calibration becomes tractable using deterministic gradient-informed algorithms.

CALPHAD↗
DOE OSTI2026

Benchmarking of massively parallel phase-field codes for directional solidification

We present a detailed benchmark comparing two state-of-the-art phase-field implementations for simulating alloy solidification under experimentally relevant conditions. The study investigates the directional solidification of Al-3wt%Cu under high-velocity solidification conditions and SCN-0.46wt% camphor under microgravity conditions from National Aeronautics and Space Administration (NASA) DECLIC-DSI-R experiments. Both codes, one employing finite-difference discretization with uniform mesh and GPU-acceleration (GPU-PF) and the other one employing finite-element discretization with adaptive-mesh and CPU-parallelization (PRISMS-PF), solve the same quantitative phase-field formulation that incorporates an anti-trapping current for the solidification of dilute alloys. We evaluate the predictions of each code for dendritic morphology, primary spacing, and tip dynamics in both 2D and 3D, as well as their numerical convergence and computational performance. While existing benchmark problems have primarily focused on simplified or small-scale simulations, they do not reflect the computational and modeling challenges posed by employing experimentally relevant time and length scales. Our results provide a practical framework for assessing phase-field code performance as well as validating and facilitating their application in integrated computational materials engineering (ICME) workflows that require integration with realistic experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Solidification cracking of refractory alloys: a computational and machine learning study to investigate composition-dependence for improved weldability and additive manufacturability

Large-batch numerical, CALculation of PHAse Diagrams (CALPHAD)-based solidification cracking calculations are performed and then analyzed with machine learning methods to generate models that relate chemistry of refractory alloys to cracking susceptibility. Kou’s solidification cracking index is used to study the refractory alloys including O, N, C binary mixtures with Mo, Ta, Nb, and W, the molybdenum-based TZM, Niobium-based C103, and Tantalum-based T111 and Ta-10 W, as well as hypothetical refractory ternary alloys. Findings strongly validate Kou’s Crack Susceptibility Index (CSI) against Varestraint test data for Nb- and Ta-based alloys, establishing CSI thresholds where refractory alloys with CSI < 15,000 K are likely weldable, CSI > 15,000 K are prone to cracking, and CSI > 25,000 K are likely unweldable (or unprintable). Furthermore, interstitial elements C, N, and O significantly increase crack susceptibility, with some existing material specifications coinciding with peak cracking susceptibility concentrations. Finally, machine learning-derived elemental potency factors enable rapid prediction of CSI from alloy chemistry for C103, TZM, Ta-10 W, and T-111 alloys. These results provide practical guidance for feedstock selection, powder reuse limits, and alloy specification amendments for welding and additive manufacturing applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Searching for Dwarf Galaxies in Gaia DR2 Phase-space Data Using Wavelet Transforms

We present a wavelet-based algorithm to identify dwarf galaxies in the Milky Way in Gaia DR2 data. Our algorithm detects overdensities in 4D position–proper-motion space, making it the first search to explicitly use velocity information to search for dwarf galaxy candidates. We optimize our algorithm and quantify its performance by searching for mock dwarfs injected into Gaia DR2 data and for known Milky Way satellite galaxies. Comparing our results with previous photometric searches, we find that our search is sensitive to undiscovered systems at Galactic latitudes |b| > 20° and with half-light radii larger than the 50% detection efficiency threshold for Pan-STARRS1 (PS1) at (i) absolute magnitudes of –7 < M V < –3 and distances of 32 kpc < D < 64 kpc, and (ii) M V < –4 and 64 kpc < D < 128 kpc. Based on these results, we predict that our search is expected to discover 5 ± 2 new satellite galaxies: four in the PS1 footprint and one outside the Dark Energy Survey and PS1 footprints. We apply our algorithm to the Gaia DR2 data set and recover ~830 high-significance candidates, out of which we identify a "gold standard" list of ~200 candidates based on cross-matching with potential candidates identified in a preliminary search using Gaia EDR3 data. All of our candidate lists are publicly distributed for future follow-up studies. Here, we show that improvements in astrometric measurements provided by Gaia EDR3 increase the sensitivity of this technique; we plan to continue to refine our candidate list using future data releases.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Lifetime of Almost Strong Edge-Mode Operators in One-Dimensional, Interacting, Symmetry Protected Topological Phases

Almost strong edge-mode operators arising at the boundaries of certain interacting one-dimensional symmetry protected topological phases with Z 2 symmetry have infinite temperature lifetimes that are nonperturbatively long in the integrability breaking terms, making them promising as bits for quantum information processing. We extract the lifetime of these edge-mode operators for small system sizes as well as in the thermodynamic limit. For the latter, a Lanczos scheme is employed to map the operator dynamics to a one-dimensional tight-binding model of a single particle in Krylov space. We suggest this model to be that of a spatially inhomogeneous Su-Schrieffer-Heeger model with a hopping amplitude that increases away from the boundary, and a dimerization that decreases away from the boundary. We associate this dimerized or staggered structure with the existence of the almost strong mode. Thus, the short time dynamics of the almost strong mode is that of the edge mode of the Su-Schrieffer-Heeger model, while the long time dynamics involves decay due to tunneling out of that mode, followed by chaotic operator spreading. Furthermore, we show that competing scattering processes can lead to interference effects that can significantly enhance the lifetime.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Characterization of Cs 3 Sb photocathodes at cryogenic temperatures

Here, we report measurements of quantum efficiency (QE) and mean transverse energy (MTE) from Cs 3 Sb photocathodes in a wide range of photon energies at both room and cryogenic temperatures. Our measurements show a strong temperature dependence of MTE even at photon energies well above threshold, indicating the presence of strong inelastic scattering of excited electrons during transport before emission into vacuum. We also demonstrate a cathode cooling method that largely preserves the QE while reducing MTE, allowing us to achieve MTEs as low as 58 meV with 3% QE in green light from Cs 3 Sb photocathodes. Our results are crucial for producing brighter electron beams for various photoinjector applications like ultrafast electron diffraction and microscopy, x-ray free-electron lasers, and particle colliders.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Assessing the sources of reionization: a spectroscopic case study of a 30× lensed galaxy at z ~ 5 with Ly α , C IV, Mg ιι, and [Ne ιιι]

We present a detailed spectroscopic analysis of a galaxy at z ≃ 4.88 that is, by chance, magnified ~30× by gravitational lensing. Only three sources at z ≳ 5 are known with such high magnification. This particular source has been shown to exhibit widespread, high equivalent width ${{\rm C\, \small {IV}}}\ \lambda 1549\, \mathring{\rm A}$ emission, implying it is a unique example of a metal-poor galaxy with a hard radiation field, likely representing the galaxy population responsible for cosmic reionization. Using ultraviolet (UV) nebular line ratio diagnostics, Very Large Telescope (VLT)/X-shooter observations rule out strong active galactic nuclei (AGN) activity, indicating a stellar origin of the hard radiation field instead. We present a new detection of ${[{\rm Ne\, \small {III}}]}\ \lambda 3870\, \mathring{\rm A}$ and use the [${\rm Ne\, \small {III}}$]/[${\rm O\, \small {II}}$] line ratio to constrain the ionization parameter and gas-phase metallicity. Closely related to the commonly used [${\rm O\, \small {III}}$]/[${\rm O\, \small {II}}$] ratio, our [${\rm Ne\, \small {III}}$]/[${\rm O\, \small {II}}$] measurement shows this source is similar to local ‘Green Pea’ galaxies and Lyman-continuum leakers. It furthermore suggests this galaxy is more metal poor than expected from the fundamental metallicity relation, possibly as a consequence of excess gas accretion diluting the metallicity. Finally, we present the highest redshift detection of ${{\rm Mg\, \small {II}}}\ \lambda 2796\, \mathring{\rm A}$, observed at high equivalent width in emission, in contrast to more evolved systems predominantly exhibiting ${\rm Mg\, \small {II}}$ absorption. Strong ${\rm Mg\, \small {II}}$ emission has been observed in most z ~ 0 Lyman-continuum leakers known and has recently been proposed as an indirect tracer of escaping ionizing radiation. In conclusion, this strongly lensed galaxy, observed just $300\, \mathrm{Myr}$ after reionization ends, enables testing of observational diagnostics proposed to constrain the physical properties of distant galaxies in the James Webb Space Telescope (JWST)/Extremely Large Telescope (ELT) era.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Spin and orbital metallic magnetism in rhombohedral trilayer graphene

We provide a theoretical interpretation of the metallic broken spin-valley (flavor) symmetry states recently discovered in hole-doped rhombohedral trilayer (ABC) graphene in large electric displacement fields. Our conclusions about the phase diagram and phase transitions combine insights from ABC graphene electronic structure models and mean-field theory, and are guided by the precise magneto-oscillation Fermi-surface-area measurements of recent experiments. We find that the principle of momentum-space condensation plays a key role in determining Fermi-surface reconstructions enabled by broken flavor symmetries when the single-particle bands imply thin annular Fermi seas. The reconstructed Fermi sea consists of one large outer Fermi-surface-enclosed majority-flavor states in reciprocal-space area A maj and one or more small inner holelike Fermi-surface-enclosed minority-flavor states in A min that are primarily responsible for nematic order. The competing ground states (valley-Ising, valley-XY, and spin-polarized state) have different A maj /A min and exchange energy maximizes this ratio and selects valley-XY nematic metal as the lowest-energy state. Here, we discuss how the nematic pockets explain the observed fractionalization of quantum oscillation frequencies, and propose anisotropic transport and the nonlinear Hall effect as additional observables.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Renormalized classical theory of quantum magnets

Here, we derive a renormalized classical spin (RCS) theory for 𝑆 >1/2 quantum magnets by constraining a generalized classical theory that includes all multipolar fluctuations to a reduced CP 1 phase space of dipolar SU(2) coherent states. When the spin Hamiltonian $\hat{ℋ}$(𝑆) is linear in the spin operators $\hat{𝑺}$ 𝑗 for each lattice site 𝑗, the RCS Hamiltonian $\tilde{ℋ}$ cl coincides with the usual classical model ℋ cl = lim 𝑆→∞⁡ $\hat{ℋ}$(𝑆). In the presence of nonlinear terms, however, the RCS theory is more accurate than ℋ cl . For the many materials modeled by spin Hamiltonians with (nonlinear) single-ion anisotropy terms, the use of the RCS theory is essential to accurately model phase diagrams and to extract the correct Hamiltonian parameters from neutron-scattering data.

magnetic anisotropy↗
DOE OSTI2024

Theoretical studies of chemical reactions related to the formation and growth of polycyclic aromatic hydrocarbons (PAH) and molecular properties of their key intermediates (Final Progress Report)

The formation mechanisms of polycyclic aromatic hydrocarbons, (PAHs) – organic molecules carrying fused benzene rings – are of great interest to scientists and engineers due to their importance in combustion chemistry and astrochemistry. On Earth, PAHs are largely produced in incomplete combustion of fossil fuel and are considered as critical precursors to unwanted soot particles leading to combustion inefficiency and causing air pollution along with detrimental health effects. Simple PAH molecules initially formed in the gas phase, are further involved in a build-up process in combustion flames leading to larger PAH, bowl-shaped nanostructures, fullerenes, and solid-phase species including carbonaceous dust, graphene particles, and soot. In deep space, PAH and their derivatives are potential key intermediates and nucleation sites leading eventually to carbonaceous nanoparticles (“interstellar grains”). Therefore, the understanding of the key processes in the synthesis of PAHs along with their precursors and their degradation mechanisms in combustion systems and in interstellar, circumstellar, and planetary atmospheric environments will provide critical insights into how complex aromatic structures, carbonaceous nanoparticles, and fullerenes are formed and destroyed. Achieving this understanding is an important step in the development of the efficient combustion processes and of the ecofriendly devices with reduced environmental pollution as well as technological strategies for the production of hydrogen and solid carbon through thermal or plasma-assisted pyrolysis of natural gas and biomass. Also, the understanding of the key processes of PAH and soot growth will help in our comprehension of chemical evolution in the universe. Detailed information on the mechanisms and reliable rate constants of the key elementary chemical reactions involved in PAH formation and destruction processes and in inception of soot particles is often missing, with the main deficiencies being the absence of temperature- and pressure-dependent rate constants for the broad range of conditions occurring in various terrestrial and interstellar processes and the lack of data on the reaction products and their branching ratios. Complementary to experimental studies, these gaps in knowledge can be filled by using quantum chemical calculations of reaction potential energy surfaces providing us with accurate energies of reaction products, intermediates, and transition states, revealing the reaction mechanism, and giving the molecular properties required to compute rate constants for relevant reaction steps and product branching ratios using the RRKM-Master Equation (ME) method. Molecular dynamics (MD) simulations can be used in cases when a reaction rate cannot be properly described by statistical theories. During the terminal renewal project period we employed these ab initio/RRKM-ME and MD approaches to complete our studies on several key reactions relevant to the formation/growth of PAH and inception of soot particles including (1) the reaction mechanism and kinetics of the resonance stabilized fulvenallenyl radical with propargyl and C 3 H 4 isomers; (2) the reaction mechanism and kinetics for the C + indene and C 2 + styrene reactions producing naphthyl or azulenyl radicals in low-temperature environments; (3) the MD study of non-equilibrium dimerization of acepyrene and coronene and its radical. The information derived from our theoretical calculations contributed to a better fundamental understanding of the reaction mechanisms and provide missing critical kinetic data to improve combustion models of hydrocarbon fuels and astrochemical models of the growth of carbonaceous molecules and particles in cold molecular clouds, circumstellar envelopes, and planetary atmospheres.

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

Machine learning Frenkel Hamiltonian parameters to accelerate simulations of exciton dynamics

In this manuscript, we develop multiple machine learning (ML) models to accelerate a scheme for parameterizing site-based models of exciton dynamics from all-atom configurations of condensed phase sexithiophene systems. This scheme encodes the details of a system’s specific molecular morphology in the correlated distributions of model parameters through the analysis of many single-molecule excited-state electronic-structure calculations. These calculations yield excitation energies for each molecule in the system and the network of pair-wise intermolecular electronic couplings. In this work, we demonstrate that the excitation energies can be accurately predicted using a kernel ridge regression (KRR) model with Coulomb matrix featurization. We present two ML models for predicting intermolecular couplings. The first one utilizes a deep neural network and bi-molecular featurization to predict the coupling directly, which we find to perform poorly. The second one utilizes a KRR model to predict unimolecular transition densities, which can subsequently be analyzed to compute the coupling. We find that the latter approach performs excellently, indicating that an effective, generalizable strategy for predicting simple bimolecular properties is through the indirect application of ML to predict higher-order unimolecular properties. Such an approach necessitates a much smaller feature space and can incorporate the insight of well-established molecular physics.

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

Active learning emulators for nuclear two-body scattering in momentum space

In this work we extend the active learning emulators for two-body scattering in coordinate space with error estimation, recently developed by Maldonado et al. [Phys. Rev. C 112, 024002], to coupled-channel scattering in momentum space. Our full-order model (FOM) solver is based on the Lippmann-Schwinger integral equation for the scattering t-matrix as opposed to the radial Schrödinger equation. We use (Petrov-)Galerkin projections and high-fidelity calculations at a few snapshots across the parameter space of the interaction to construct efficient reduced-order models (ROMs), trained by a greedy algorithm for locally optimal snapshot selection. Both the FOM solver and the corresponding ROMs are implemented efficiently in Python using Google's JAX library. We present results for emulating scattering phase shifts in coupled and uncoupled channels and cross sections, and assess the accuracy of the developed ROMs and their computational speedup factors. We also develop emulator error estimation for both the t-matrix and the total cross section. The software framework for reproducing and extending our results is publicly available. Together with our recent advances in developing active-learning emulators for three-body scattering, these emulator frameworks set the stage for full Bayesian calibrations of chiral nuclear interactions and optical models against scattering data with quantified emulator errors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Thermal Property Modeling and Assessment of the Physical Properties of FLiNaK

Here, the thermodynamic properties of the LiF–NaF–KF system and its subsystems were re-evaluated using the CALPHAD method, focusing on the deviations from ideality in the heat capacity [C p (T)] of mixtures, temperature-dependence of mixing enthalpies (Δ mix H), vapor pressures of mixtures, and enthalpies of fusion (ΔH fus ) using reported values and measurements. The results of this work confirmed a pseudoternary eutectic at 732.9 K and 46.5LiF–11.5NaF–42KF mol % (FLiNaK) with determined Cp(T) = 69.822 + 0.000679T + 2,137,152T –2 J mol –1 K –1 (732.9 < T < 1200 K) and ΔH fus = 18,051 J mol –1 in good agreement to the experimental values. The resulting accurate thermodynamic representation can be used to find relevant excess thermochemical properties in the whole pseudoternary space and provide insights into their composition and temperature dependence. A detailed evaluation of select thermophysical properties was also conducted to determine recommended density, thermal diffusivity, thermal conductivity, vapor pressures, and boiling point with calculated uncertainties over the temperature range of interest for FLiNaK–salt energy applications.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

COMPUTATIONAL MODELING OF IGNITION AND PREMIXED FLAME PROPAGATION INITIATED BY A PRE-CHAMBER TURBULENT JET

Addressing the pressing need for reduced carbon emissions, Turbulent Jet Ignition (TJI) emerges as a promising technology for ultra-lean combustion, offering enhanced thermal efficiencies and minimized cyclic variability in spark-ignited engines. To facilitate rapid testing and integration of this technology, a robust computational modeling framework is crucial. This study delves into the predictive capabilities of computational models for main-chamber ignition and premixed flame propagation using a single-cycle TJI rig measured by Biswas et al. (Applied Thermal Engineering, vol 106, 2016). Employing an open-source compressible flow simulation solver with Large Eddy Simulation (LES) for turbulence modeling, the investigation integrates the conventional Laminar Finite Rate Chemistry (LFRC) model alongside the transported Probability Density Method (PDF) for turbulence-chemistry interaction. A fully-consistent Eulerian Monte-Carlo Fields (EMCF) method is utilized to approximate the transported PDF, while Interaction by Exchange with Mean is employed to close micro-mixing terms in stochastic differential equations. A reduced chemical reaction mechanism with 21 species and 84 reactions (DRM-19) is used for solving chemical kinetics, and a double Gaussian energy deposition model is used to approximate the spark ignition in the pre-chamber. An unstructured O-grid mesh with 0.3 million cells in the prechamber and 1 million cells in the main chamber is employed. Results are divided into two phases: pre-chamber initialization and full TJI simulations. Validation of the predicted pre-chamber flame propagation and the lean ignition in the main-chamber is carried out by using available experimental data. Under quiescent conditions, both the LFRC and transported PDF methods largely underestimate the flame speed and subsequent pressure growth in the pre-chamber. A linear momentum forcing technique is applied to investigate the impact of initial turbulence in the pre-chamber, demonstrating a notable influence on flame propagation. Fine-tuning of the forcing coefficient reproduces the sudden pressure growth observed in the experiment. The experimentally validated pre-chamber simulation serves as the initial condition for the full TJI simulations. It is found that the LFRC model fails to predict lean-ignition in the main-chamber, resulting in a misfiring event. Incorporation of turbulence-chemistry interaction using the transported PDF method substantially improves the prediction of the ignition event in the main-chamber, achieving fair qualitative agreement and quantitative validation of combustion parameters within ±10% of the reported experimental data. The rich simulation results consisting of a full set of statistical description of the thermo-chemical states enable us to gain deep insights into the ignition mechanisms in the main chamber, which is limited when done experimentally. A novel dual ignition phenomenon is revealed in the TJI rig for the first time. Initially, a primary ignition kernel is formed at a downstream location which eventually detaches from the main jet. As the jet momentum decreases, a secondary ignition event follows, this time at a more upstream location which eventually combines with the primary ignition kernel to form a single connected flame front. Investigation of these ignition sequences in chemical composition space reveal distinct differences between the two. The primary ignition event in the main-chamber is followed by a large concentration of active radicals from the pre-chamber jet, accelerating the chain-branching steps, characterizing what has been referred to as flame ignition. In contrast, the secondary ignition occurs in the absence of active radicals in the pre-chamber jet, hence characterized as jet ignition. Further analysis of the effect of pre-chamber jet characteristics on lean ignition in the main-chamber is conducted by setting up cases with different initial pressure ratios (por) between the two chambers, a non-dimensional parameter, ranging from 1.2 to 3.2. As the initial pressure ratio increases, jet momentum increases, with dual ignition observed in cases above por= 2.2. Case with por= 3.2 lead to misfiring. The effect of ignition sequence on global combustion characteristics of TJI is analyzed. Dual ignition events lead to non-monotonicity in combustion characteristics such as global reaction progress variable, flame penetration, and global heat release rate. In dual ignition events, although the rate of fuel consumption and global heat release rate is initially lower, the secondary ignition leads to a sudden increase in flame surface area, resulting in a sudden jump and promoting the overall performance of the TJI system.

42 ENGINEERING↗
DOE OSTI2024

Deep Learning-Based Dynamic Modeling of Three-Phase Voltage Source Inverters

Inverter-based resource (IBR) models are necessary to analyze modern power system stability and create effective control strategies. Modeling IBRs in converter-rich power systems is crucial, yet challenging due to the lack of commercial information on converter topologies and control parameters. This paper proposes novel convolutional neural network (CNN)–based data-driven techniques for modeling IBRs, addressing adaptability and proprietary concerns without requiring internal system physics knowledge. The proposed method is tested using real grid-tied commercial IBR transient data and demonstrates effectiveness and accuracy. Furthermore, the developed modeling approach is integrated and implemented in the open-source power distribution simulation and analysis tool, GridLAB-D, to illustrate the potentiality of dynamic analysis of large-scale power systems with high IBRs.

deep learning, artificial intelligence↗
DOE OSTI2022

Time-series forecasting using manifold learning, radial basis function interpolation, and geometric harmonics

We address a three-tier numerical framework based on nonlinear manifold learning for the forecasting of high-dimensional time series, relaxing the “curse of dimensionality” related to the training phase of surrogate/machine learning models. At the first step, we embed the high-dimensional time series into a reduced low-dimensional space using nonlinear manifold learning (local linear embedding and parsimonious diffusion maps). Then, we construct reduced-order surrogate models on the manifold (here, for our illustrations, we used multivariate autoregressive and Gaussian process regression models) to forecast the embedded dynamics. Finally, we solve the pre-image problem, thus lifting the embedded time series back to the original high-dimensional space using radial basis function interpolation and geometric harmonics. The proposed numerical data-driven scheme can also be applied as a reduced-order model procedure for the numerical solution/propagation of the (transient) dynamics of partial differential equations (PDEs). In conclusion, we assess the performance of the proposed scheme via three different families of problems: (a) the forecasting of synthetic time series generated by three simplistic linear and weakly nonlinear stochastic models resembling electroencephalography signals, (b) the prediction/propagation of the solution profiles of a linear parabolic PDE and the Brusselator model (a set of two nonlinear parabolic PDEs), and (c) the forecasting of a real-world data set containing daily time series of ten key foreign exchange rates spanning the time period 3 September 2001–29 October 2020.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗