SEARCH · Engineering Papers

Results for “Reacting”

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.

SourceResearch connectionFrom yearThrough year

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

At least 937 records · Page 52

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Imaging the photochemistry of cyclobutanone using ultrafast electron diffraction: Experimental results

We investigated the ultrafast structural dynamics of cyclobutanone following photoexcitation at λ = 200 nm using gas-phase megaelectronvolt ultrafast electron diffraction. Our investigation complements the simulation studies of the same process within this special issue. It provides information about both electronic state population and structural dynamics through well-separable inelastic and elastic electron scattering signatures. We observe the depopulation of the photoexcited S 2 state of cyclobutanone with n3s Rydberg character through its inelastic electron scattering signature with a time constant of (0.29 ± 0.2) ps toward the S 1 state. The S 1 state population undergoes ring-opening via a Norrish Type-I reaction, likely while passing through a conical intersection with S 0 . The corresponding structural changes can be tracked by elastic electron scattering signatures. These changes appear with a delay of (0.14 ± 0.05) ps with respect to the initial photoexcitation, which is less than the S 2 depopulation time constant. This behavior provides evidence for the ballistic nature of the ring-opening once the S 1 state is reached. The resulting biradical species react further within (1.2 ± 0.2) ps via two rival fragmentation channels yielding ketene and ethylene, or propene and carbon monoxide. Furthermore, our study showcases the value of both gas-phase ultrafast diffraction studies as an experimental benchmark for nonadiabatic dynamics simulation methods and the limits in the interpretation of such experimental data without comparison with such simulations.

Carbon monoxide↗

DOE OSTI2025

Potentials of mean force fail to describe chemical bond-breaking in solution

Many liquid phase studies assume that the potential energy surfaces of reacting molecules are the same as in the gas phase, neglecting complex solvent dynamics that can completely alter the nature of chemical reactivity. Even studies that include solvent effects typically only consider them in an average, equilibrium way as part of a potential of mean force (PMF). In this work, we use mixed quantum/classical simulations to compare how equilibrium and non-equilibrium solvent motions affect the photodissociation of a simple diatomic molecule, NaK + , in liquid tetrahydrofuran. A PMF analysis shows that as the excited-state molecule dissociates with the solvent at equilibrium, the bonding electron remains associated with K + at short bond distances but eventually localizes on Na + at the end of dissociation. When we examine non-equilibrium dynamical photodissociation trajectories, however, we find that they fall into three distinct categories: about a quarter of them have the bonding electron mainly associated with Na + , another quarter stay mainly associated with K + , and about half have the bonding electron shared roughly equally between the two ions. The results show that equilibrium PMFs cannot accurately describe the dynamics of bond-breaking chemical reactions in solution because there is insufficient time for the solvent to reach equilibrium on the time scale over which bond dissociation occurs. Furthermore, our analysis shows that the solvent coupling between the electronic energy surfaces is similar at and away from equilibrium, suggesting that other factors, such as solute velocity-driven solvent memory effects, play a more important role in explaining the failure of the equilibrium PMF to predict the non-equilibrium dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Structural basis for Zika envelope domain III recognition by a germline version of a recurrent neutralizing antibody

Recent epidemics demonstrate the global threat of Zika virus (ZIKV), a flavivirus transmitted by mosquitoes. Although infection is usually asymptomatic or mild, newborns of infected mothers can display severe symptoms, including neurodevelopmental abnormalities and microcephaly. Given the large-scale spread, symptom severity, and lack of treatment or prophylaxis, a safe and effective ZIKV vaccine is urgently needed. However, vaccine design is complicated by concern that elicited antibodies (Abs) may cross-react with other flaviviruses that share a similar envelope protein, such as dengue virus, West Nile virus, and yellow fever virus. This cross-reactivity may worsen symptoms of a subsequent infection through Ab-dependent enhancement. To better understand the neutralizing Ab response and risk of Ab-dependent enhancement, further information on germline Ab binding to ZIKV and the maturation process that gives rise to potently neutralizing Abs is needed. In this work, we use binding and structural studies to compare mature and inferred-germline Ab binding to envelope protein domain III of ZIKV and other flaviviruses. We show that affinity maturation of the light-chain variable domain is important for strong binding of the recurrent VH3-23/VK1-5 neutralizing Abs to ZIKV envelope protein domain III, and identify interacting residues that contribute to weak, cross-reactive binding to West Nile virus. These findings provide insight into the affinity maturation process and potential cross-reactivity of VH3-23/VK1-5 neutralizing Abs, informing precautions for protein-based vaccines designed to elicit germline versions of neutralizing Abs.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2020

Direct functionalization of C-H bonds by electrophilic anions

Bond formation between two of the most inert molecules, alkanes and [B12X12]2- (X = Cl, Br), is performed in a two step process. Fragmentation of [B12X12]2- in the gas phase generates highly reactive [B12X11]- ions which spontaneously react with alkanes. The reaction mechanism was investigated using tandem mass spectrometry and gas-phase vibrational spectroscopy combined with electronic structure calculations. This effort revealed the unprecedented and counterintuitive electrophilic substitution of a proton in an alkane by a negatively charged ion resulting in a B-C bond formation. The product is a dianionic [B12X11CnH2n+1]2- species, in which H+ is electrostatically bound to the dianion. High flux ion soft-landing was performed to co-deposit [B12X11]- and complex organic molecules (phthalates) in thin layers on surfaces. Molecular structure analysis of the product films revealed that C-H functionalization by [B12X11]- occured in the presence of other more reactive functional groups. This observation demonstrates the utility of highly reactive fragment ions for defined bond formation processes and may pave the way to the use of gas phase ion chemistry for chemical synthesis in the condensed phase.

Warneke, Jonas↗

DOE OSTI2022

Imaging active site chemistry and protonation states: NMR crystallography of the tryptophan synthase α-aminoacrylate intermediate

NMR-assisted crystallography-the integrated application of solid-state NMR, X-ray crystallography, and first-principles computational chemistry-holds significant promise for mechanistic enzymology: by providing atomic-resolution characterization of stable intermediates in enzyme active sites, including hydrogen atom locations and tautomeric equilibria, NMR crystallography offers insight into both structure and chemical dynamics. Here, this integrated approach is used to characterize the tryptophan synthase α-aminoacrylate intermediate, a defining species for pyridoxal-5'-phosphate-dependent enzymes that catalyze β-elimination and replacement reactions. For this intermediate, NMR-assisted crystallography is able to identify the protonation states of the ionizable sites on the cofactor, substrate, and catalytic side chains as well as the location and orientation of crystallographic waters within the active site. Most notable is the water molecule immediately adjacent to the substrate β-carbon, which serves as a hydrogen bond donor to the ε-amino group of the acid-base catalytic residue βLys87. From this analysis, a detailed three-dimensional picture of structure and reactivity emerges, highlighting the fate of the L-serine hydroxyl leaving group and the reaction pathway back to the preceding transition state. Reaction of the α-aminoacrylate intermediate with benzimidazole, an isostere of the natural substrate indole, shows benzimidazole bound in the active site and poised for, but unable to initiate, the subsequent bond formation step. When modeled into the benzimidazole position, indole is positioned with C3 in contact with the α-aminoacrylate C $\bar{β}$ and aligned for nucleophilic attack. Here, the chemically detailed, three-dimensional structure from NMR-assisted crystallography is key to understanding why benzimidazole does not react, while indole does.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Colloidal superionic conductors

Nanoparticles with highly asymmetric sizes and charges that self-assemble into crystals via electrostatics may exhibit behaviors reminiscent of those of metals or superionic materials. Here, we use coarse-grained molecular simulations with underdamped Langevin dynamics to explore how a binary charged colloidal crystal reacts to an external electric field. As the field strength increases, we find transitions from insulator (ionic state), to superionic (conductive state), to laning, to complete melting (liquid state). In the superionic state, the resistivity decreases with increasing temperature, which is contrary to metals, yet the increment decreases as the electric field becomes stronger. Additionally, we verify that the dissipation of the system and the fluctuation of charge currents obey recently developed thermodynamic uncertainty relation. Our results describe charge transport mechanisms in colloidal superionic conductors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Methanediol from cloud-processed formaldehyde is only a minor source of atmospheric formic acid

Atmospheric formic acid is severely underpredicted by models. A recent study proposed that this discrepancy can be resolved by abundant formic acid production from the reaction ( 1 ) between hydroxyl radical and methanediol derived from in-cloud formaldehyde processing and provided a chamber-experiment-derived rate constant, k 1 = 7.5 × 10 −12 cm 3 s −1 . High-level accuracy coupled cluster calculations in combination with E,J -resolved two-dimensional master equation analyses yield k 1 = (2.4 ± 0.5) × 10 −12 cm 3 s −1 for relevant atmospheric conditions ( T = 260–310 K and P = 0–1 atm). We attribute this significant discrepancy to HCOOH formation from other molecules in the chamber experiments. More importantly, we show that reversible aqueous processes result indirectly in the equilibration on a 10 min. time scale of the gas-phase reaction HCHO + H 2 O ⇌ HOCH 2 OH (2) with a HOCH 2 OH to HCHO ratio of only ca . 2%. Although HOCH 2 OH outgassing upon cloud evaporation typically increases this ratio by a factor of 1.5–5, as determined by numerical simulations, its in-cloud reprocessing is shown using a global model to strongly limit the gas-phase sink and the resulting production of formic acid. Based on the combined findings in this work, we derive a range of 1.2–8.5 Tg/y for the global HCOOH production from cloud-derived HOCH 2 OH reacting with OH. The best estimate, 3.3 Tg/y, is about 30 times less than recently reported. The theoretical equilibrium constant K eq (2) determined in this work also allows us to estimate the Henry’s law constant of methanediol (8.1 × 10 5 M atm −1 at 280 K).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Autonomous phototaxis of hydrogel swimmers

The design of synthetic soft matter capable of emulating the complex behaviors of living organisms, such as sensing and adapting to their environment, remains an important challenge in developing biomimetic materials. Functionalized hydrogels are ideal candidates for such materials since they are highly responsive to their environment and can be operated in water. In this work, we investigate a hybrid bonding hydrogel composed of peptide amphiphile supramolecular nanofibers covalently attached to a photoresponsive network, in which high-aspect-ratio ferromagnetic nanowires are aligned along the length of the sample, designed to swim under oscillating magnetic fields. This hybrid hydrogel swimmer can autonomously swim toward a light source by utilizing photoinduced interactions between supramolecular and covalent networks reminiscent of phototactic swimming in living systems. Using a combination of experimental techniques and a continuum model incorporating photochemistry, magnetoelasticity, and hydrodynamics, we explain the swimming mechanism and predict phototactic behavior. Our work highlights the potential role of hybrid bonding polymers, which leverage the interplay between supramolecular assemblies and covalent networks. We demonstrate how these polymers can be tailored to react dynamically to their environment, paving the way for developing intelligent and autonomous robotic systems.

Science & Technology - Other Topics↗

DOE OSTI2025

Allosteric inhibition of the IZUMO1–JUNO fertilization complex by the naturally occurring antisperm antibody OBF13

Sperm IZUMO1 binds to egg JUNO, and this interaction is essential for mammalian fertilization. Isolated from a female mouse immunized with syngeneic sperm, the antisperm antibody OBF13 recognizes IZUMO1 and inhibits murine fertilization. How OBF13 interferes with sperm–egg interactions was unknown. Here, we present the X-ray crystal structure of IZUMO1 in complex with OBF13. OBF13 binds to the apex of the four-helix domain of IZUMO1, distant from the JUNO-binding site. Our crystal structure of OBF13-bound IZUMO1 resembles apo-IZUMO1 and differs from the structure of IZUMO1 in complex with JUNO. We identify that OBF13 carries a low level of somatic hypermutation, and through deep mutational scanning, we engineer an affinity-enhanced OBF13 variant. This OBF13 variant single-chain fragment variable decreases the apparent affinity of IZUMO1 for membrane-bound murine JUNO and blocks the binding of acrosome-reacted sperm to eggs, thereby preventing fertilization. We propose allostery between the OBF13 epitope and the JUNO-binding site. OBF13 inhibits a conformational change in IZUMO1, preventing fusion-competent sperm from adhering to murine eggs during fertilization. Surprisingly, murine IZUMO1 binds to hamster JUNO with an affinity ~20-fold higher than to murine JUNO. The decreased affinity caused by OBF13 of murine IZUMO1 for hamster JUNO is sufficient for murine sperm to bind to and fuse with hamster eggs. Our studies provide a structural and mechanistic framework for species-specific, allosteric inhibition of IZUMO1 by a naturally occurring antisperm antibody and offer insights into the development of immunocontraceptives.

Science & Technology - Other Topics↗

DOE OSTI2020

Structure and function of a flavin-dependent S-monooxygenase from garlic ( Allium sativum )

Allicin is a component of the characteristic smell and flavor of garlic ( Allium sativum ). A flavin-containing monooxygenase (FMO) produced by A. sativum (AsFMO) was previously proposed to oxidize S-allyl- L -cysteine (SAC) to alliin, an allicin precursor. In this paper, we present a kinetic and structural characterization of AsFMO that suggests a possible contradiction to this proposal. Results of steady-state kinetic analyses revealed that AsFMO exhibited negligible activity with SAC; however, the enzyme was highly active with L -cysteine, N-acetyl- L -cysteine, and allyl mercaptan. We found that allyl mercaptan with NADPH was the preferred substrate-cofactor combination. Rapid-reaction kinetic analyses showed that NADPH binds tightly ($K_D$ of ~2 μM) to AsFMO and that the hydride transfer occurs with pro-R stereospecificity. We detected the formation of a long-wavelength band when AsFMO was reduced by NADPH, probably representing the formation of a charge-transfer complex. In the absence of substrate, the reduced enzyme, in complex with NADP + , reacted with oxygen and formed an intermediate with a spectrum characteristic of C4a-hydroperoxyflavin, which decays several orders of magnitude more slowly than the $k$ cat . The presence of substrate enhanced C4a-hydroperoxyflavin formation and, upon hydroxylation, oxidation occurred with a rate constant similar to the $k$ cat . The structure of AsFMO complexed with FAD at 2.08-Å resolution features two domains for binding of FAD and NADPH, representative of class B flavin monooxygenases. These biochemical and structural results are consistent with AsFMO being an S-monooxygenase involved in allicin biosynthesis through direct formation of sulfenic acid and not SAC oxidation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Determination of Synthetic Opioids Belonging to the Fentanyl Class in Silt Using Electron Ionization Gas Chromatography-Mass Spectrometry (GC-MS).

An extraction protocol from silt sediment of fentanyl and three analogs: acetylfentanyl, thiofentanyl and acetylthiofentanyl, spiked at two concentrations each and separately (at ~1 and ~10 µg/g), is described. In addition, the identity of the fentanyls preliminarily identified by electron ionization gas chromatography-mass spectrometry (EI-GC-MS) analysis, can be corroborated by reacting each opioid in the silt’s extract with 2,2,2-trichloroethoxycarbonyl chloride (Troc-Cl). Further, reaction between Troc-Cl and each opioid generates two unique products that can be used to retrospectively identify the original opioid therefore serving as a corroborating tool for known opioids as well as new, unknown fentanyl analogs.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2021

Halogenation of used aluminum matrix test reactor fuel – a bench-scale demonstration with surrogate materials

In this work, experiments with surrogate materials were performed at bench scale to demonstrate a halogenation technique applicable to treatment of used aluminum matrix test reactor fuel. The technique involves dissolution and separation of aluminum from used aluminum matrix test reactor fuel in molten-halide salt systems prior to treatment and disposition of the fuel’s uranium and fission products. Demonstration of the halogenation technique was performed with neodymium metal as a non-radiological surrogate for uranium metal. Experiments involved blending forms of aluminum and neodymium metal with ammonium and lithium chloride or ammonium and lithium bromide, which upon heating decomposed into ammonia gas and the respective hydrogen chloride or bromide gas. The latter reacted with the metals to form the respective aluminum and neodymium halides. At elevated temperatures, aluminum halides gasified away from the respective neodymium halides, which fused with their respective lithium halides. Samples of fused and distillate salts were collected and analyzed, yielding extents of aluminum removal that ranged from 94.5–98.2% for chlorination runs and 91.4–97.8% for bromination runs. No neodymium was detected in the distillate fractions. Some experiments were repeated with excess reactants, and a portion of aluminum chloride distillate was processed into a consolidated waste form.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

Conceptual Design of a Thermal–to-14-MeV Neutron Conversion Device for Use in the Advanced Test Reactor

Here, this paper details the conceptual design of a thermal-to-14-MeV neutron converter consisting of a mixture of lithium and deuterium in a blanket material. Such a device operates a two-step reaction, first generating tritons via thermal neutron absorption in the tritium breeding material, and in the second step, high-energy neutrons are produced either via deuterium-tritium fusion reaction or with tritium reacting with lithium. A thermal-to-14-MeV neutron converter significantly hardens the neutron spectrum by virtually removing thermal neutrons and adding a high-energy 14 MeV component to the neutron spectrum. While similar concepts have been previously proposed and tested in other reactors, the unique characteristics of the Advanced Test Reactor (ATR), namely its important thermal flux (up to $10^{15}$ n.cm$^{-2}$.s$^{-1}$) makes it markedly attractive for obtaining a very large fast neutron flux, usable for irradiation studies under neutron flux conditions prototypical of fusion reactors. The paper provides a description of a new computational scheme developed for handling the coupled neutron-triton transport mechanism using the Geant4 toolkit. Resulting neutron spectra and high-energy neutron yields are summarized for different irradiation positions and potential neutron breeder materials. Maximum predicted thermal-to-14-MeV neutron yields are on the order of $2.10^{-4}$, which is consistent with previous studies found in the literature. Thus, when placed inside the ATR, such a neutron converter will be providing the largest high-energy neutron source available for activation and irradiation studies of materials foreseen for use in fusion reactors. Future steps will involve qualifying the computational scheme using the ATR critical facility using activation foil measurements.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

Loss-of-Coolant Accident Analysis of a High-Burnup Pressurized Water Reactor Core Design Using Gadolinia-Doped UO 2

High-burnup and extended enrichment fuels are of interest for extending the cycle lengths of pressurized water reactors from 18 to 24 months. Changes to the fuel design and core loading scheme have potentially significant safety implications due to power distribution effects, reduced thermal conductivity, and/or increased plenum pressures due to additional burnable poison loadings. Additionally, higher burnups result in increased material degradation and risk of fuel fragmentation, relocation, and dispersal (FFRD). A representative 24-month core design using gadolinia-doped UO 2 was analyzed for performance under large-break (LB) loss-of-coolant accident (LOCA) conditions using PARCS, RELAP5-3D, and BISON. Furthermore, all considered acceptance criteria were met, with no cases exceeding the 1477 K maximum cladding temperature or the post-quench ductility oxidation limit of 17% equivalent cladding reacted. Full-core FFRD susceptibility was estimated to be approximately 455 kg, though high uncertainties exist with current approaches for computing susceptibility. Undoped fuel rods are more likely to be limiting due to higher linear heat rates. Relatively high burnup and linear heat rate rods located in second batch assemblies are of greatest safety significance during LB LOCA for this high-burnup core design.

High burnup↗

DOE OSTI2026

Review of ECCS Acceptance Criteria and Experimental Basis Evolution Toward Fuel Fragmentation, Relocation, and Dispersal Studies

The U.S. nuclear industry is pursuing extensions of light water reactor (LWR) fuel burnup and enrichment limits to approximately 75 GWd/t and 10 wt.% 235 U to achieve economic and operational benefits. A central safety consideration in this effort is the behavior of high burnup (HBu) fuel during loss-of-coolant accidents (LOCAs), particularly fuel fragmentation, relocation, and dispersal (FFRD). Here, this work provides a historical and technical review of U.S. LOCA regulation and experimentation, clarifying how the evolution of Emergency Core Cooling System (ECCS) acceptance criteria in 10 CFR 50.46 has shaped both testing approaches and interpretations of fuel safety. The study revisits the original intent of the ECCS criteria, showing that the peak cladding temperature and equivalent cladding reacted limits were developed as surrogates to preserve a coolable geometry. The explicit inclusion of the coolable geometry criterion in the regulation was intended to emphasize the underlying safety philosophy and as a safeguard against unforeseen failure modes, an intent that remains directly relevant to modern concerns regarding FFRD. The review traces the lineage of HBu LOCA experiments to the Argonne National Laboratory furnace tests, from which subsequent programs at Studsvik, Halden, and Oak Ridge National Laboratory were derived. These tests employed a 5 °C/s heating rate inherited from early embrittlement studies, a stylized temperature history that does not represent actual LWR LOCA thermal-hydraulics. Comparison of these test conditions to pressurized water reactor large break LOCAs and separate effects data indicates that the existing HBu LOCA database may not be fully applicable to all LWR LOCA scenarios, from which a qualitative framework for applicability is proposed.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2020

Criegee intermediates: production, detection and reactivity

In the context of tropospheric chemistry, Criegee intermediates denote carbonyl oxides with biradical/zwitterionic character (R 1 R 2 COO) that form during the ozonolysis of alkenes. First discovered almost 70 years ago, stabilised versions of Criegee intermediates formed via collisional removal of excess energy have interesting kinetic and mechanistic properties. The direct production and detection of these intermediates were not reported in the literature until 2008. However, recent advances in their generation through the ultraviolet irradiation of the corresponding diiodoalkanes in excess O 2 and detection by various spectroscopic techniques (photoionisation, ultraviolet, infrared, microwave and mass spectrometry) have shown that these species can react rapidly with closed-shell molecules, in many cases at or exceeding the classical gas-kinetic limit, via multiple reaction pathways. These reactions can be complex, and laboratory measurements of products and the temperature and pressure dependence of the reaction kinetics have also revealed unusual behaviour. The possible role of these intermediates in atmospheric chemistry is significant, altering models of the oxidising capacity of the Earth's atmosphere and the rate of generation of secondary organic aerosol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

A physics-constrained neural ordinary differential equations approach for robust learning of stiff chemical kinetics

The high computational cost associated with solving for detailed chemistry poses a significant challenge for predictive computational fluid dynamics (CFD) simulations of turbulent reacting flows. While deep learning techniques have been explored to develop faster surrogate models, they often fail to integrate reliably with CFD solvers. This instability arises because traditional deep learning approaches optimize for training error without ensuring compatibility with ordinary differential equation (ODE) solvers, resulting in accumulation of errors over time. Recently, neuralODE (NODE) based approaches have been shown to be a promising technique to emulate and accelerate detailed chemistry computations. Here, in the present work, we extend this NODE framework for stiff chemical kinetics by incorporating mass conservation constraints directly into the loss function during training. This ensures that the total mass as well as the individual elemental species masses are conserved in an a-posteriori manner. Proof-of-concept studies are performed with the novel physics-constrained NODE (PC-NODE) approach for homogeneous autoignition of hydrogen-air mixture over a range of composition and thermodynamic conditions. It is demonstrated that the PC-NODE framework not only improves the physical consistency of the resulting data-driven model with respect to mass conservation criteria, but also improves training efficiency. PC-NODE is shown to achieve 2–100× speedup relative to the hydrogen-air detailed chemical mechanism depending on the type of the ODE solver (implicit or explicit) used during autoregressive inference tests. Lastly, a-posteriori studies are performed wherein the trained PC-NODE model is coupled with a CFD solver. It is shown that higher accuracy is achieved with PC-NODE relative to the purely data-driven NODE approach. Moreover, PC-NODE also exhibits robustness and generalizability to unseen initial conditions from within (interpolative capability) as well as outside (extrapolative capability) the training regime.

computational combustion↗

DOE OSTI2026

Data-based filtered dissipation rate modelling for multi-modal turbulent combustion: evaluating a priori model generalizability

Manifold-based models offer a computationally efficient alternative to directly transporting the thermochemical state in computational simulations of turbulent reacting flows, projecting the high-dimensional thermochemical state-space onto a low-dimensional manifold. Recent efforts have yielded a manifold-based model applicable to multi-modal combustion, enabling reconstruction of the thermochemical state from solutions to two-dimensional manifold equations in mixture fraction and generalized progress variable that are parameterised by three scalar dissipation rates. In coarse-grained simulations such as Large Eddy Simulation (LES), closure of the multi-modal manifold equations and subfilter variances/covariance requires closure of three filtered scalar dissipation rates. Here, the present work adopts a data-based approach, providing closure for the three filtered scalar dissipation rates via deep neural networks (DNNs). High-fidelity datasets corresponding to an autoigniting n-dodecane jet flame and a bluff body swirl-stabilized confined lifted spray flame of two aviation fuels (Jet-A and C1) with different ignition propensities are leveraged to generate training data that spans a diverse range of thermodynamic conditions and combustion modes, including low- and high-temperature ignition regimes in addition to premixed and nonpremixed behaviour. A final DNN model is trained to enforce inherent physical constraints by learning nonlinear functional transformations of the three filtered scalar dissipation rates. The generalizability of this constrained DNN model is demonstrated a priori via conditional statistics evaluated on the lifted spray flame with C1–a configuration that had not been included in the training data. Excellent DNN agreement with conditional DNS statistics is observed, and integrated gradients are computed to identify the most sensitive input variables. The similarity of the marginal PDFs of the most informative input variables and outputs across configurations are quantified via the Wasserstein metric, demonstrating that data-based models may successfully generalize to unseen parametric conditions so long as the most informative input variables share similar distributions across training and testing datasets.

Data-based modelling↗