SEARCH · Engineering Papers
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.
Reactor noise is a nondestructive approach to measure neutron/gamma fluctuations. Extract kinetic parameters such the prompt neutron decay constant (a). Used in applications such as safeguards and reactor physics. This work address the impact on precision from frequency domain analysis methodology.
Department of Energy National Nuclear Security Administration Office of Material Management and Minimization’s mission includes the conversion of civilian research reactors from highly enriched uranium (HEU) to low-enriched uranium (LEU) fuel. Analyses have shown that the Oak Ridge National Laboratory High Flux Isotope Reactor (HFIR) will need to operate at 95 MW for the LEU silicide dispersion fuel designs to match key performance metrics obtained with HEU fuel at 85 MW. To prove safe operation of HFIR after installation of plant modifications to increase power, a high power HEU test cycle was proposed. Neutronics model updates and reactor physics analyses are performed to support the development of safety design reports for the high power (HP) HEU test cycle. Reactor physics metrics evaluated herein include fuel depletion, actinide production, cycle length, fission rate density distributions, axial power peaking factors, and reactor kinetic parameters. These reactor physics analyses support the development of future LEU safety design reports by providing key input for future HP HEU HFIR thermal hydraulics and reactor transient safety analyses.
Li-air batteries are considered a potential alternative to Li-ion batteries for transportation applications due to their high theoretical specific energy. Most works in this area focus on use of O 2 as the reactant. However, newer concepts for using gaseous reactants (such as CO 2 , which has a theoretical specific energy density of 1,876 Wh/kg) provide opportunities for further exploration. The main objective of this project was the development of a novel strategy that enables operation of Li-CO 2 batteries at high-capacity and high-rate, with a long-cycle-life. The team was able to: (1) Synthesize two novel transition metal chalcogenide (TMC) catalysts that work in synergy with ionic liquid-based electrolytes to enhance the efficiency of reactions during discharge and charge processes; (2) Fabricate high-porosity cathode electrodes with 3D printing to increase electrode surface area and gas permeability; (3) Develop a multiscale modeling framework that integrates Density Functional Theory (DFT), Ab-Initio Molecular Dynamics (AIMD), classical Molecular Dynamics (MD), and Finite Element Analysis (FEA) to investigate atomic and cell-level properties of Li-CO 2 batteries; (4) Assemble a stackable Li-CO 2 pouch-cell able to deliver a capacity of >200 mAh. These achievements were realized through an integrated approach based on materials synthesis, testing, characterization, analysis, and computation. This project produced a thorough understanding of key chemical, electronic, and kinetic parameters that govern the operation of Li- CO 2 batteries in realistic conditions. The methodologies employed, and the insight generated, will be valuable beyond advancing the field of Li-CO 2 batteries
Reactor noise is a nondestructive approach to measure neutron/gamma fluctuations. Extract kinetic parameters such the prompt neutron decay constant (a). Used in applications such as safeguards and reactor physics. This work address the impact on precision from frequency domain analysis methodology.
This is the final report for this project. The long-term objective of our project is to develop a general molecular-level picture of the surface chemistry associated with ALD processes. Our central hypothesis is that the chemistry of ALD precursors can differ significantly from that seen in solution. Awareness of such differences should provide general guidelines on what to watch for when synthesizing new ALD precursors and designing and optimizing new ALD processes, especially in terms of minimizing the incorporation of impurities in the growing films and of carrying out depositions under mild pressure and temperature conditions. The main objective of our research project has been to advance the fundamental knowledge of the surface chemistry of ALD precursors needed for the design and optimization of film deposition processes. A modern surface-science approach has been implemented to both elucidate the mechanism of the reactions of the precursors on the surface and characterize the composition and morphology of the growing films. In general, emphasis is being placed on: • Identifying the primary reactions that may lead to the deposition of clean films; •Identifying the secondary reactions that may help the ALD process, by, for instance, helping with the reduction (or oxidation) of the metal atom; • Identifying the secondary reactions that may lead to the deposition of undesirable impurities in the growing films; • Determining the kinetic parameters of the relevant surface reactions in order to define the optimum conditions for film deposition and to minimize impurity deposition. • Characterizing the nature of the resulting films, with focus on their stoichiometry and on the final oxidation states of the constituent elements; and • Using the information obtained to propose better precursors for given ALD processes.
Effective electrochemical CO 2 reduction to liquid fuels requires that the local catalytic environment facilitates the desired reactivity, yet a microscopic understanding of this environment is difficult to achieve from experiment alone. In this work, a 3D reaction-diffusion model was developed to explore the effects of electrode surface area and local geometry on the performance of a heterogeneous catalyst that performs a two-step CO 2 reduction cascade reaction to CO and then CH 3 OH under aqueous conditions. Kinetic parameters for the model were inspired by experimental results using a cobalt phthalocyanine (CoPc) catalyst. Three-dimensional architectures composed of arrays of square pillars with varying dimensions and either smooth or periodically modulated surfaces were tested, revealing the extent to which geometry modulates the performance of the cascade reactions. Although structural variations modulate local concentration gradients, we find that electrochemically active surface area predominantly governs the overall cascade reaction. Moreover, the results suggest that supersaturation of CO, with concentrations up to ten-fold higher than the equilibrium solubility limit, might be critical for more efficient conversion to CH 3 OH. For any given geometry, the spatially averaged ratio of [CO] to [CO 2 ] is dictated by the electrochemically active surface area and determines the yield of CH 3 OH. For a fixed surface area, geometries that spatially confine the electrolyte yield moderate local [CO] to [CO 2 ] ratios within small volumes. In contrast, less confining geometries result in a broader distribution of local ratios spread over larger volumes, with both configurations yielding the same spatially averaged [CO] to [CO 2 ] ratio. These insights provide valuable design principles—highlighting the critical importance of surface area and possibly CO supersaturation—for engineering advanced electrode architectures that leverage intermediate trapping and CO supersaturation to enhance overall performance in tandem CO 2 reduction systems.
This milestone report aims to identify and assess the cluster dynamics (CD) modeling requirements within the Department of Energy's Office of Nuclear Energy (DOE-NE) Advanced Materials and Manufacturing Technologies (AMMT) program and to communicate these needs to the DOE-NE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. The goal is to ensure NEAMS is well-informed about the CD modeling requirements to support AMMT's mission of accelerating the development, qualification, demonstration, and deployment of advanced structural materials and manufacturing for nuclear energy applications. CD modeling is an essential tool for predicting the degradation of structural materials under irradiation, which is a key component of AMMT's accelerated qualification process. The AMMT program focuses on both additively manufactured and wrought structural alloys, such as laser powder-bed fusion 316H austenitic stainless steel, alloy 709, Haynes 244, and alloy 617. These materials require a generalized CD modeling framework to facilitate rapid model development and computational simulation. A flexible, generalized CD software, similar to the Multiphysics Object-Oriented Simulation Environment (MOOSE) finite element framework, would enable modeling of various cluster types, including defect clusters, defect-solute clusters, and multicomponent clusters, incorporating thermodynamics and kinetics parameters. Radiation effects, microstructural feature evolution, and multi-dimensional modeling are critical considerations for the CD model. The usability of the CD code should allow for easy modification and coupling with MOOSE-based simulations. Additionally, the software should adhere to Nuclear Quality Assurance-1 standards, include a testing suite for verification and validation, and be version-controlled within a national laboratory-managed Git repository. Benchmark problems are needed to assess code predictions and performance.
Tortuosity-weighted interfacial flux for lithium (TWIF-Li) predicts through-thickness Li gradients in thick composite all-solid-state cathodes without fitted parameters. Image-derived microstructures, GITT-derived concentration-dependent solid diffusion, and tortuosity-weighted interfacial kinetics reproduce operando neutron radiography across practical rates, delivering transferable design rules to suppress transport-limited reaction fronts.
Titanium oxide site requirements for the anchoring of catalytically active forms of dispersed vanadia species are investigated. Using a series of oxygen-modified titanium nitride materials, we systematically modify the available vanadia-anchoring moieties on the titanium-bearing support surface. This strategy in turn results in a very narrow size distribution of well-dispersed vanadia clusters anchored on the support. This synthesis technique ensures that only highly dispersed vanadia is present on the catalyst surface even at nominal monolayer coverages. A deeper insight into the catalytic behavior for ethanol partial oxidation of the active vanadia species is obtained by comparing intrinsic kinetic and thermodynamic parameters of kinetically relevant reaction steps (heat of ethanol adsorption, hydrogen abstraction activation energy, and the energetics of catalyst reoxidation). In conclusion, it is found that these parameters are dependent on the relative amount of oxygen in the titanium-bearing support and the resulting distribution of vanadia species, which validates the premise of a potential participation of lattice oxygen atoms of the titania support into the catalytic activity of active vanadia species.
Microwave (MW)-assisted catalytic pyrolysis offers a promising pathway for efficient plastic upcycling. This work develops an integrated modeling framework combining dynamic data reconciliation, a temperature-dependent rate model, and a yield model to represent the time-varying production rate of components in MW-assisted LDPE pyrolysis conducted in a batch reactor. An Arrhenius-type rate model with a temperature-dependent reaction order is developed. A biexponential correlation is proposed for the yield of gaseous products that enables to capture the evolving product formation behavior during conversion. In the yield correlation, one term is used to represent the initial increase in yield, reflecting the rapid formation of intermediate or primary products at the early stages of the reaction when a larger fraction of the reactant remains available. As conversion progresses, the influence of this term gradually diminishes. The other term accounts for the subsequent decrease in the predicted yield, representing secondary reactions such as further cracking or coke formation that reduce the concentration of certain products at higher conversion. The model is found to accurately represent reconciled experimental flow rate profiles from an in-house MW-assisted catalytic batch reactor for major products, including ethylene, ethane, 1-butene, and benzene, across 250−350 °C. Ethylene remains the dominant product but decreases from about 41.95% at 250 °C to 30.14% at 350 °C, while heavier products increase significantly, with 1-butene rising to nearly 8.37% and benzene reaching 2.17% at intermediate temperatures. The model shows that the ethylene production rate can be maximized at around 270 °C. The models developed in this work can be utilized for process optimization, reactor design and scale-up of microwave-assisted plastic conversion technologies, and economic analysis.
Electrocatalytic activation of alkanes can further decarbonize chemical manufacturing by leveraging affordable renewable electricity and readily available shale gas reserves in the United States. Earlier works have identified the unique role of Pt in adsorbing and activating alkanes, like propane, at room temperature in acidic, aqueous electrolytes, revealing spontaneous formation of deeply dehydrogenated propane-derived surface species with an intact C 3 - backbone. Although an adsorption mechanism was hypothesized, it has not been explicitly investigated to date, preventing the quantification of kinetic rate parameters. A robust methodology to investigate and benchmark propane adsorption kinetics on Pt is critical for the rational design of electrocatalysts that exhibit higher selectivity toward desired partially oxidized products. Herein, we analyze an oxidative current transience that appears during the adsorption of propane on Pt in aqueous electrochemical conditions and develop a methodology that elucidates the adsorption mechanism and enables quantification of rate parameters such as order dependences and apparent activation barriers. This method yields an expected first-order dependence with respect to propane concentration at low coverage and reveals a second-order dependence with respect to the concentration of surface active sites. Additionally, the apparent activation barrier for propane adsorption was calculated using an Arrhenius analysis of the current transience under temperature control. The experimentally measured activation barrier of 35 kJ mol –1 is in excellent agreement with the theoretical barrier calculated by density functional theory (DFT). The kinetic analysis was extended, via the use of transition state theory, to extract entropy and enthalpy of activation, yielding consistent results with the proposed two-step adsorption mechanism and DFT calculations. These results demonstrate reliable quantification of kinetic parameters for electrocatalytic activation of C–H bonds in alkanes that can be employed for rational catalyst development for a versatile range of electrocatalytic conditions.
Plasma pyrolysis of methane is an effective route for zero-carbon hydrogen production. Yet, soot generated from pyrolysis of hydrocarbons is detrimental to the climate and human health. There is ample experimental and theoretical evidence that suggests polycyclic aromatic hydrocarbons (PAHs) are the molecular precursors to soot particles. The reaction pathways of PAH formation are intricately dependent on a multitude of process parameters, whose kinetic mechanisms are not well-understood in plasma pyrolysis. This project aims to leverage advances in the kinetic modeling of soot formation in combustion, as well as in surrogate modeling and active learning, to systematically investigate the effects of process parameter on the kinetics of PAH formation in plasma pyrolysis of methane. To this end, we propose to use the PAH formation kinetics model developed by the PPPL/PU group based on the well-established ABF and HACA mechanisms, coupled with low-temperature plasma models. We will develop an active learning (AL) framework based on Bayesian optimization to systematically and data-efficiently explore the complex and multivariable parameter space of plasma pyrolysis in order to quantify the effects of plasma and feed parameters on the ABF and HACA kinetic pathways. AL is the branch of machine learning concerned with systematically querying samples from a system (experimental or computational) to train a data-driven model that maps design parameters to a performance criterion. We will use the data generated via AL to perform global sensitivity analysis, combined with uncertainty quantification, to elucidate the impact of different reaction pathways on minimizing formation of soot precursors. This study will result in an improved understanding of kinetics of PAH formation in plasma pyrolysis and can pave the way for more advanced mechanistic studies (e.g., soot nucleation mechanisms). Additionally, the findings will be useful for establishing practical strategies for increasing the pyrolysis efficiency and producing high-grade carbon for synthesis of nanomaterials.
Inosine 5′-monophosphate dehydrogenase (IMPDH) is a promising antibiotic target. This enzyme catalyzes the NAD-dependent oxidation of inosine 5′-monophosphate (IMP) to xanthosine 5′-monophosphate (XMP), which is the rate-limiting step in guanine nucleotide biosynthesis. Bacterial IMPDH-specific inhibitors have been developed that bind to the NAD + site. These inhibitors display varied affinities to different bacterial IMPDHs that are not easily rationalized by X-ray crystal structures of enzyme–inhibitor complexes. Inspection of X-ray crystal structures of 25 enzyme–inhibitor complexes, including 10 newly described, suggested that a mobile active site flap may be a structural determinant of inhibitor potency. Saturation transfer difference NMR experiments also suggested that the flap may contact the inhibitors to varying extents in different IMPDHs. Flap residue Leu413 contacted some inhibitors but was not structured in the crystal structures of other inhibitor complexes. The substitution of Leu413 with Phe or Ala in Bacillus anthracis IMPDH had inhibitor-selective effects, suggesting residue 413 could be a structural determinant of affinity. Curiously, the Ala substitution increased the potency of most inhibitors, even those that contacted Leu413 in the crystal structures. Presteady-state and steady-state kinetics experiments showed that the Leu413Ala substitution had comparable effects on inhibitor binding to the noncovalent E·IMP complex and the covalent intermediate E-XMP*, suggesting that the flap had similar interactions in both complexes. These results demonstrate that contacts do not necessarily indicate favorable interactions, and poorly structured mobile regions should not be discounted when assessing binding determinants.
The microbial enzyme nitrogenase catalyzes the MgATP-dependent reduction of N 2 to 2NH 3 , a transformation central to the global nitrogen cycle. While the canonical Thorneley−Lowe (TL) kinetic model has long served as a mechanistic framework, it does not incorporate several recent insights. Here, we present an updated kinetic model for Monitrogenase that incorporates these new findings. A significant insight is that electron transfer (ET) from the reduced Fe protein to the FeMo-cofactor is gated by MgATP-dependent conformational transitions and can be described as a probabilistic event that is dependent on the ligand bound to the active-site metallocofactor. The updated kinetic model quantitatively reproduces steady-state product formation rates across a broad range of experimental conditions, yielding revised estimates for key rate constants. It is demonstrated that under N 2 turnover, the probability of productive ET to the active site decreases by ∼60%, resulting in a significant fraction of Fe protein cycles that are unproductive for electron delivery. This mechanistic feature explains the observed rate limitation in N 2 reduction and implies a revised minimum energetic cost of approximately 25 MgATP per N 2 reduced. Integrating these new features into the revised kinetic model provides a more complete and usable foundation for understanding nitrogenase catalysis.
We propose using highly excited cyclotron states of a trapped electron to detect meV axion and dark-photon dark matter, marking a significant improvement over our previous proposal and demonstration [One-electron quantum cyclotron as a milli-ev dark-photon detector, .]. When the axion mass matches the cyclotron frequency ω c , the cyclotron state is resonantly excited, with a transition probability proportional to its initial quantum number, n c . The sensitivity is enhanced by taking n c ∼ 10 6 ( 0.1 meV ω c ) 2 . By optimizing key experimental parameters, we minimize the required averaging time for cyclotron detection to t ave ∼ 10 − 6 s, permitting detection of such a highly excited state before its decay. An open–end-cap trap design enables the external photon signal to be directed into the trap, rendering our background-free detector compatible with large focusing cavities, such as the BREAD proposal, while capitalizing on their strong magnetic fields. Furthermore, the axion conversion rate can be coherently enhanced by incorporating layers of dielectrics with alternating refractive indices within the cavity. Collectively, these optimizations enable us to probe the QCD axion parameter space from 0.1 to 2.3 meV (25–560 GHz), covering a substantial portion of the predicted postinflationary QCD axion mass range. This sensitivity corresponds to probing the kinetic mixing parameter of the dark photon down to ε ≈ 2 × 10 − 16 . Published by the American Physical Society 2025
The ultralight dark photon is a well-motivated, hypothetical dark matter candidate. In a dilute plasma, they can resonantly convert into photons, and heat up the intergalactic medium between galaxies. In this work, we explore the dark photon dark matter parameter space by comparing synthetic Lyman- α forest data from cosmological hydrodynamical simulations to observational data from VLT/UVES of the quasar HE0940-1050 ( z em = 3.09 ). We use a novel flux normalization technique that targets underdense gas, reshaping the flux probability distribution. Not only do we place robust constraints on the kinetic mixing parameter of dark photon dark matter, but notably our findings suggest that this model can still reconcile simulated and observed Doppler parameter distributions of z ∼ 0 Lyman- α lines, as seen by HST/COS. This work opens new pathways for the use of the Lyman- α forest to explore new physics, and can be extended to other scenarios such as primordial black hole evaporation, dark matter decay, and annihilation.
The Jovian magnetic field, being the strongest and largest planetary one in the solar system, could offer us new insights into possible microscopic scale new physics, such as a non-zero mass of the Standard Model (SM) photon or a light dark photon kinetically mixing with the SM photon. We employ the immense data set from the latest Juno mission, which provides us unprecedented information about the magnetic field of the gas giant, together with a more rigorous statistical approach compared to the literature, to set strong constraints on the dark photon mass and kinetic mixing parameter, as well as the SM photon mass. The constraint on the dark photon parameters is independent of whether dark photon is (part of) dark matter or not, and serves as the most stringent one in a certain regime of the parameter space.