Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “chemical modeling”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Imaging solvated oxygen atoms with a femtosecond laser

As a powerful oxidant, atomic oxygen (O) holds considerable promise for a variety of biomedical and industrial applications. However, the inability to quantify solvated oxygen atoms has prevented the determination of the fundamental parameters governing its behaviour in relevant aqueous environments. Here, we directly image ground-state oxygen atoms in water using femtosecond two-photon absorption laser-induced fluorescence. Measurements show that oxygen atoms persist for tens of microseconds in water, penetrating hundreds of micrometres into the liquid. This observed longevity has significant implications, suggesting a need to re-evaluate existing models of solvated atomic oxygen reactivity and transport. Beyond atomic oxygen, this technique is broadly applicable to other solvated atomic species of interest, including nitrogen (N) and hydrogen (H). This work establishes that radical atomic species can be quantified in liquid with ultrafast laser spectroscopy, providing the basis for the determination of key properties including reaction rates, chemical lifetimes, and Henry’s law constants.

Fluorescence spectroscopy↗
DOE OSTI2025

Atmospheric-pressure ammonia synthesis on AuRu catalysts enabled by plasmon-controlled hydrogenation and nitrogen-species desorption

The Haber–Bosch process for ammonia synthesis contributes up to ~3% of global greenhouse gas emissions. Plasmonic catalysts strongly concentrate light and can alter the reaction intermediates via out-of-equilibrium processes, providing the potential for an alternative, less-energy-intensive pathway to synthesize ammonia. Here, in this study, we show that gold-ruthenium (AuRu) bimetallic nanoparticles can synthesize ammonia at room temperature and pressure using visible light. We create AuRu alloys with varying compositions and achieve ammonia production rates of ~60 μmol per gram of catalyst bed per hour. In situ infrared spectroscopy reveals that light accelerates the hydrogenation of nitrogen intermediates compared to conventional thermal catalysis. Through computational modelling, we demonstrate that photo-excited electrons enable associative hydrogenation pathways for nitrogen activation rather than direct nitrogen–nitrogen bond breaking. This light-assisted mechanism requires both hydrogen and light working together to overcome the nitrogen activation barrier, mimicking how biological enzymes produce ammonia naturally and providing fundamental insights for developing sustainable, energy-efficient chemical synthesis.

Yuan, Lin [Stanford Univ., CA (United States)] (OR↗
DOE OSTI2025

Temperature-based reactive flow model for triaminotrinitrobenzene (TATB) plastic bonded explosives

A new reactive flow model is presented for triaminotrinitrobenzene (TATB)-based plastic bonded explosives, applicable to shock initiation and steady detonation problems of differing initial temperature. Temperature disequilibrium is assumed between unreacted explosive, material in the vicinity of compressed defects (called hot spots), and reaction products. The model incorporates temperature-dependent decomposition reaction rates. Particularly, Arrhenius model parameters were derived from quantum-based molecular dynamics simulations of TATB decomposition. Further, a model of detonation carbon aggregation is incorporated, describing the slow release of energy inherent to detonation in TATB-based materials. Model parameters were calibrated against gas gun shock initiation experiments and steady detonation rate stick tests. The predictive ability of the model in the shock initiation regime is tested against recent thin pulse experiments. The model is found to perform equally well in predicting the size-effect curve of ambient, cold, and hot rate sticks. The present work demonstrates the viability of incorporating results from subscale simulations into a continuum-scale reactive flow model.

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

VZr_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys vanadium-zirconium (V-Zr). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach \cite{PAW}. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements V and Zr. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

TaTi_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys tantalum-titanium (Ta-Ti). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Hf and Nb. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on Air Force HPC11 cluster using the VASP 6.5.1. Additional methodology and file structure information is available in the dataset README.txt file.

36 MATERIALS SCIENCE↗
DOE OSTI2025

HfNb_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys hafnium-niobium (Hf-Nb). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Hf and Nb. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect BCC lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,029 randomized atomic structures over 31 chemical compositions. Further methodological and structural information is contained in the dataset README.txt file.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A Novel Integrated Fermentation Process with Engineered Microbial Consortia for Butanol Production from Lignocellulose Sugars without CO 2 Emission

The goal of this project was to develop a synthetic microbial consortium consisting of a lactic acid bacterium, a carboxydotrophic acetogen, and a solventogenic clostridia for the production of n butanol, an advanced biofuel and industrial chemical, from lignocellulose sugars (mainly glucose and xylose) and formate (produced from CO 2 by electrochemical reduction) in an integrated bioprocess (biorefinery), which can provide an effective solution to the technical challenges in developing energy and carbon optimized synthesis for the bioeconomy and achieve the program objectives of ARPA-E. The project focused on the design, modeling and construction of synthetic microbial consortia consisting of three bacterial species to maximize carbon conversion and butanol production with a 100% theoretical product yield from glucose and zero or negative CO 2 emission.

09 BIOMASS FUELS↗
DOE OSTI2025

Preliminary techno-economic assessment of gas switching reforming (GSR) of natural gas for pure hydrogen production and power generation with integrated CO2 capture

The increasing demand for hydrogen and the CO2 intensity of natural gas (NG) reforming motivate the development of low-carbon-emission hydrogen production technologies. Gas Switching Reforming (GSR) with integrated CO2 capture, a technology based on Chemical Looping Reforming (CLR), has been experimentally proven and shows potential for scale-up. In this study, select oxygen carriers (OC) (NiO/Al2O3, Fe2O3-CeO2/Al2O3, and magnetite) were tested in methane steam reforming in a fixed bed reactor to determine their relative reactivities under relevant conditions for GSR (800 °C, 7 bar total pressure). Process models were then developed to perform techno-economic analysis (TEA) of GSR for hydrogen production (GSR-H2) and a combined cycle (GSR-CC) in which high-purity H2 is fired in a gas turbine to produce electricity. Operating at 10 bar and 1100 °C and with the additional recovery steps implemented increased H2 production by ∼ 30% and improved efficiency relative to prior studies. For GSR-H2, the levelized cost of hydrogen (LCOH) is 1.61–1.64 $/kg-H2, competitive with a reference SMR case, though operating and maintenance costs are higher due to increased electricity demand. GSR-CC has a significantly higher levelized cost of electricity (LCOE) than its reference NGCC (natural gas combined cycle) plant, suggesting it is less competitive; however, increasing production scale could make it more attractive. Life-cycle results for GSR-H2 indicate NG consumption drives ∼ 75% of total global warming impacts (∼2.3 kg CO2 eq/kg H2). An environmental, health, and safety screening suggests iron-based carriers are comparatively safer, whereas NiO may pose greater risks. Overall, GSR-H2 is a scalable, competitive option for hydrogen production using nickel and non-nickel OC.

03 NATURAL GAS↗
DOE OSTI2025

Mesoscale Modeling Approach for Quantifying Microstructure-Aware Micromechanical Responses in Metal Hydrides

Metal hydrides can undergo significant volume changes upon hydrogen uptake and release, which induce a mechanical response that depends not only on the evolving hydrogen composition but also on the microstructure. We present a comprehensive mesoscale modeling framework based on microelasticity theory to quantify the micromechanical responses of metal hydrides, specifically focusing on a hydrogenating polycrystalline MgH 2x particle within a host material as a model micromechanical system. Utilizing digitally generated realistic microstructures and density-functional-theory-derived parameters, we analyzed highly nonuniform local stress profiles in the polycrystalline hydrides under the clamping force exerted by the host during hydrogenation. Our framework also allows us to predict the corresponding strain energy accumulation and mechanical hot spots formation in the hydrides, highlighting their roles in thermodynamic destabilization and mechanical failure, respectively. Through extensive parametric simulations, we further quantified the influence of interface type, crystallinity, grain size, loading ratio, and host stiffness, providing practical guidance for optimizing microstructural design and host material selection. This proposed approach is broadly applicable to micromechanical systems with complex microstructural features involving chemical reaction- and/or phase-transformation-induced deformation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Upcycling waste polystyrene to adipic acid through a hybrid chemical and biological process

Oxidative catalytic depolymerization of polystyrene (PS) can produce benzoic acid, but the annual consumption of benzoic acid is ~40 times lower than PS. For this catalytic oxidation method to be a viable means to manage PS waste, benzoic acid should be converted to higher-volume chemicals. We demonstrate a hybrid chemical and biological process that uses PS as feedstock for production of adipic acid, a high-volume co-monomer for nylon 6,6 via benzoic acid. Mn/Br co-catalyzed autoxidation of PS to benzoic acid proceeds with a yield of up to 94% in a solvent mixture of benzoic acid and water. The PS-derived benzoic acid undergoes bioconversion at near-quantitative yield to muconic acid, which is readily converted to adipic acid through catalytic hydrogenation. Process modeling, techno-economic analysis, and life cycle assessment estimate an adipic acid minimum selling price of $3.18/kg, with a 61% decrease in greenhouse gas emissions relative to production from fossil fuels.

09 BIOMASS FUELS↗
DOE OSTI2025

Expediting field-effect transistor chemical sensor design with neuromorphic spiking graph neural networks

Improving the sensitive and selective detection of analytes in a variety of applications requires accelerating the rational design of field-effect transistor (FET) chemical sensors. Achieving high-performance detection relies on identifying optimal probe materials that can effectively interact with target analytes, a process traditionally driven by chemical intuition and time-consuming trial-and-error methods. To address the difficulties in probe screening for FET sensor development, this work presents a methodology that combines neuromorphic machine learning (ML) architectures, specifically a hybrid spiking graph neural network (SGNN), with an enriched dataset of physicochemical properties through semi-automated data extraction using large language models. Achieving a classification accuracy of 0.89 in predicting sensor sensitivity categories, the SGNN model outperformed traditional ML techniques by leveraging its ability to capture both global physicochemical properties and sparse topological features through a hybrid modeling framework. Next-generation sensor design was informed by the actionable insights into the connections between material properties and sensing performance offered by the SGNN framework. Through virtual screening for the detection of per- and polyfluoroalkyl substances (PFAS) as a use case, the effectiveness of the SGNN model was further validated. Density functional theory simulations confirmed graphene as a promising active material for PFAS detection as suggested by the SGNN framework. By bridging gaps in predictive modeling and data availability, this integrated approach provides a strong foundation for accelerating advancements in FET sensor design and innovation.

Ferreira, Rodrigo Pires [Univ. of Chicago, IL (Uni↗
DOE OSTI2026

A Coupled Theoretical and Experimental Approach to Elucidating the Mechanisms of Methyl Esters

The goal of this project was to probe fundamental kinetics questions regarding the gas phase reactive behavior of oxygenates. In order to succeed, our program necessitated both the development of new experimental tools and the development of theory based kinetic models to explain the chemistry of oxygenate fuels. By using advanced hybrid additive manufacturing/ traditional manufacturing techniques to create specialized ceramic micro-reactors, the team was able to simulate extreme high-temperature environments with much greater precision and durability than was previously possible. These tiny, high-tech tubes allow researchers to capture and identify "fleeting" chemical species—molecules that exist for only a fraction of a millisecond—using sophisticated light sources and mass spectrometry. Additionally, the team supported the development of a tabletop VUV laser system for isomer detection without a synchrotron. Through this work, we have successfully mapped out the specific chemical pathways of various oxygen-rich fuels, solving long-standing mysteries about how these substances break down kinetically. Ultimately, this research provides the fundamental knowledge needed to design next-generation engines and fuels that are better for the environment.

09 BIOMASS FUELS↗
DOE OSTI2025

Environmental photosensitizers can exhibit enhanced actinic absorption in microhydrated clusters compared to solution

Brown carbon chromophores at environmental air-water interfaces often act as photosensitizers that absorb sunlight and subsequently transfer energy to nearby molecules, initiating a wide variety of chemical reactions. Despite their importance to understanding daytime chemistry at these airwater interfaces, little is known about the role of the solvation environment on the photophysical properties of these photosensitizers. In this work, we present a joint experimental-theoretical study of the vibrational and photophysical properties of microhydrated protonated and deprotonated 4-benzoylbenzoic acid (4-BBA), a key model system for environmental photosensitizers. We find that for protonated 4-BBAH + ∙(H 2 O) 0-1 , representing photosensitizers in very acidic conditions, a single bright state dominates the UV-Vis spectrum between 280 and 400 nm. Comparing the experimental UV-Vis spectra and quantum chemistry-predicted spectra of 4-BBA + ∙(H 2 O) 0-2 , we find that the degree of microhydration has little effect on the UV-Vis spectra or the orbitals of the dominant feature. For deprotonated 4-BBA ‒ , representing photosensitizers in basic conditions, quantum chemistry calculations predict that the UV-Vis spectra are ~3x weaker in intensity than the brightest 4-BBAH + ∙(H 2 O) 0-1 features and were not observed experimentally. Quantum chemistry calculations predict one dominant UV-Vis feature is present in the spectra of 4-BBA ‒ ∙(H 2 O) 0-2 , which exhibit minor shifts with degree of microhydration. While 4-BBA in bulk solution over a range of pH values has relatively weak absorption within the solar actinic region, we show that microhydrated 4-BBA has bright transitions within the actinic region. This indicates that the complex structure of environmental air-water interfaces can shift the absorption maximum of photosensitizers into the actinic region for enhanced absorption of sunlight and subsequent enhancement of photosensitizer-driven reactions.

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

Mechanically Accelerated Depolymerization of Entangled Linear Polymer Melts

Mechanical forces can enhance the chemical depolymerization of synthetic polymers when shear flow accelerates chain scission. To quantify the extent of mechanically-accelerated scission, the effect of simple shear flow (duration and strength) with low Weissenberg and Deborah numbers was investigated by considering the impact of applied work in both simple shear and shear dominated mixed flows. Hydrogenated polyisoprene was chosen as a model linear, entangled system. The conditions (strain amplitude, frequency, and shearing time) necessary to increase chain scission were assessed in the rubbery melt. Shear flow accelerated chain scission at higher temperatures, suggesting an activated process. Isothermal scission versus work curves were superposed by applying shift factors a T,S , whose Arrhenius-like temperature dependence gave an apparent activation energy for chain scission of ~ 110 kJ/mol, which is likely a combination of the activation energy of viscosity and bond energy. This work provides a base for quantifying the impact of shear on depolymerization of polymer melts and highlight the connection between viscous dissipation and scission chemistry.

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

Emergent topological polarization textures in relaxor ferroelectrics

Relaxor ferroelectrics underpin high-performance actuators and sensors, yet the nature of polar heterogeneities driving their broadband dielectric response remains debated. Using a unified, multimodal structural refinement framework— simultaneously fitting complementary X-ray and neutron total scattering, X-ray absorption spectra, and diffuse scattering—we reconstruct 3D mesoscale polarization maps in the classic relaxor system PbMg 1/3 Nb 2/3 O 3 –PbTiO 3 . We uncover self-organized swirling polarization textures with half-skyrmion (meron) vortices, challenging models of independent polar nanoregions. These textures, characterized by smooth changes in the polarization direction, originate from overlapping volumes in which the projections of locally correlated polarization vectors onto each volume’s long axis share the same sign. Vortex cores correlate strongly with local charge and strain gradients imposed by compositional heterogeneities. In this work, our results suggest that chemical disorder, acting via depolarizing and strain fields, stabilizes topological vortex textures of the polarization field, offering a route for engineering new dielectric and ferroelectric functionalities.

Eremenko, Maksim [National Institute of Standards ↗
DOE OSTI2026

Integrated System Planning: Emerging Software Requirements in the Power Industry

Power system planning software remains fragmented across organizational boundaries, with specialized tools for capacity expansion, production cost modeling, power flow, and dynamic analysis operating on incompatible data models and assumptions. This article argues that the fragmentation is not merely a technical problem but a predictable consequence of Conway's law: software architectures mirror the departmental structures within which they are developed. Regulatory milestones like Federal Energy Regulatory Commission (FERC) Order 888 formalized these divisions, but the roots trace back to the distinct engineering disciplines-mechanical, chemical, and electrical-that staffed generation and transmission planning departments in vertically integrated utilities. As the industry moves toward integrated system planning (ISP) that coordinates generation, transmission, and distribution investment decisions, the software ecosystem must evolve accordingly. We identify five categories of software requirements to enable this transition: coherent data inputs decoupled from individual applications, unified and extensible data schemas, modular component representations that support multiple abstraction levels, lifecycle management of planning datasets, and well-defined application programming interface (API) contracts that separate data exchange from algorithmic control. We examine how these requirements interact with three common workflow patterns-serial gate clearing, sequential multiapplication, and convergence oriented-and discuss the interface design principles each demands. We then outline a vision for platform-based planning architectures where specialized analytical services compose through standardized interfaces and where artificial intelligence (AI)/machine learning (ML) tools augment decision support within a disciplined software infrastructure. The practices proposed here offer a path from today's siloed tool collections toward collaborative planning ecosystems capable of handling the complexity of modern power system transformation.

24 POWER TRANSMISSION AND DISTRIBUTION↗