Search NASA⌕ Search

SEARCH · Search NASA

Results for “design of computer experiments”

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

Stability analysis of the Eulerian–Lagrangian finite volume methods for nonlinear hyperbolic equations in one space dimension

In this paper, we construct a novel Eulerian–Lagrangian finite volume (ELFV) method for nonlinear scalar hyperbolic equations in one space dimension. It is well known that the exact solutions to such problems may contain shocks though the initial conditions are smooth, and direct numerical methods may suffer from restricted time step sizes. To relieve the restriction, we propose an ELFV method, where the space-time domain was separated by the partition lines originated from the cell interfaces whose slopes are obtained following the Rakine–Hugoniot junmp condition. Unfortunately, to avoid the intersection of the partition lines, the time step sizes are still limited. To fix this gap, we detect effective troubled cells (ETCs) and carefully design the influence region of each ETC, within which the partitioned space-time regions are merged together to form a new one. Then with the new partition of the space-time domain, we theoretically prove that the proposed first-order scheme with Euler forward time discretization is total-variation-diminishing and maximum-principle-preserving with at least twice larger time step constraints than the classical first order Eulerian method for Burgers’ equation. Numerical experiments verify the optimality of the designed time step sizes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Thermal Reservoir Networks for Modularly Expandable Thermal Microgrids

The Department of Defense (DoD) faces the substantial challenge of cost-effectively retrofitting one to two installations per month, each comprising approximately 1,000 buildings, to improve resilience, reduce energy consumption, and enhance energy supply security. Achieving these objectives requires optimal system selection and effective risk mitigation during system integration. To address this need, we introduce Platform-Based Design (PBD), a structured, hierarchical methodology adapted from other industrial sectors to the domain of energy system retrofits. We demonstrate the effectiveness of PBD through a techno-economic feasibility study comparing geothermal-coupled thermal energy networks (TENs) with conventional energy systems for heating, cooling, and powering 17 buildings at Joint Base Andrews (JBA) in Maryland. Our analysis illustrates that the PBD approach enables rigorous, data-driven, sequential decision making, resulting in a family of Pareto-optimal systems, among which the TEN emerged as the most promising solution. The selected TEN design integrates geothermal borefields, heat recovery heat pumps, photovoltaic (PV) arrays, and battery storage. Compared to the baseline system – gas heating combined with air-source chillers – the proposed TEN reduces annual imported energy by 74% and peak electricity demand by 45%, achieves a levelized cost of energy of $\$0.210$/kWh, and substantially enhances resilience. Life-cycle costs increase by approximately 6%, and initial investment costs are about 2.5 times higher than the baseline. However, if central plant infrastructure, district loops, and utility-scale PV and battery systems are privately funded and operated, the initial investment would fall below the baseline system cost. Critical to achieving these significant performance improvements were detailed nonlinear dynamic simulations coupling geothermal heat transfer, energy system operation, and realistic feedback control logic. These simulations identified essential design modifications and control strategy refinements that substantially reduced energy use, peak demand, and compressor shortcycling, thereby improving durability and reliability—issues that would have been significantly more expensive to resolve during operation. Additionally, the verification step highlighted sensitivities to key design parameters that could reduce initial investment by approximately $\$2$ million and reduce annual life-cycle costs more than $\$300,000$. We recommend adopting the PBD methodology for future feasibility studies and TEN pilot projects to gain valuable operational experience. Furthermore, we recommend that DoD invest in transferring and scaling the PBD methodology to other installations. This entails developing standardized computational frameworks and component libraries as well as training industry in conducting PBD. Such investments would enable rapid, robust, reliable, and cost-effective retrofits, supporting DoD’s ambitious energy system modernization goals.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Validation of Numerical Tools for Calculating Reactivity Feedback in Sodium Fast Reactors Using SEFOR Experimental Data

The Southwest Experimental Fast Oxide Reactor (SEFOR) was an experimental sodium-cooled fast breeder reactor operated from 1969 to 1972 with experiments designed to measure Doppler reactivity feedback in a wide temperature range from around 350 °F to temperatures approaching the melting point of mixed oxide fuel of around 5000 °F, providing valuable data for code validations. Co-supported by the Department of Energy (DOE) Fast Reactor Program (FRP) and the DOE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, the SEFOR benchmark project focused on using the experimental data to validate numerical tools that are used in industry and academia to design and license sodium-cooled fast reactors (SFRs). By the end of FY-25, substantial progress was achieved in the SEFOR benchmark study. A variety of numerical tools commonly used for modeling SFRs were applied to develop models for SEFOR core configurations I-D, I-E, I-I, and I-J. These included Monte Carlo codes such as MCNP, Serpent, and Shift; deterministic codes such as the legacy Argonne Reactor Computation (ARC) suite and the high-fidelity NEAMS code Griffin; and the system analysis code SAS4A/SASSYS-1 (SAS). Using these models, both SEFOR zero-power experiments and power-ascending tests were successfully simulated. Comparisons were performed against experimental measurements of core criticalities, reflector worth, kinetics parameters (Λ/βeff), isothermal reactivity feedback (from 350 °F to 760 °F at zero power), and power-ascending reactivity feedback (as power increased from 0.4 MW to 17 MW). In general, these comparisons demonstrated very good agreement between numerical results and experimental data. In Fiscal Year 26 (FY-26), the SEFOR benchmark project will continue to address the modeling issues identified in FY-25. Effort will focus on the simulation of reactivity insertion transients in SEFOR core II using the ARC/SAS model. Future work will also focus on incorporating BISON into the SEFOR core modeling process to enable the first Multiphysics simulations of the isothermal tests based on the MOOSE framework.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Spin-Controllable Dynamics in Defect-Engineered Carbon Nanotubes as Single Photon Emitters: Data-Driven Modeling and Computations

Quantum technologies, such as quantum computing and sensing, require efficient single-photon emission (SPE) sources that operate at room temperature in telecom wavelengths. While several materials can serve as SPE sources, no single platform meets all the criteria for efficiency, ambient operation, and scalability. Single-walled carbon nanotubes (SWCNTs) with covalently attached molecules offer a promising solution. Their SPE can be easily tuned via modifications of the SWCNT's diameter, chirality, and bonded molecules, enabling emission across near-IR to telecom wavelengths at ambient conditions. However, to fully realize the potential of SWCNTs and unlock their quantum capabilities, a deeper understanding of how structural defects from molecular adducts affect their emission and competing photoexcited processes is essential. To address this gap in our knowledge, this project combined quantum chemistry calculations with data-driven methods of cheminformatics (QSAR) and machine learning (ML). The developed computational approaches have provided several design strategies for covalent functionalization of SWCNTs to improve their optical response. The collaboration with Los Alamos National Lab (LANL) enabled direct comparison of computational and experimental data, facilitating method validation. This partnership was enhanced through access to LANL's Center for Integrated Nanotechnologies (CINT) utilizing User Facility Program and summer internships, which provided three NDSU graduate students with hands-on experience at LANL. The outcomes of this project included (1) Advancing the current stage of computational methods in accurate modeling of non-adiabatic spin-dependent photoexcited dynamics and its applicability to nanosystems consisting of thousands of atoms, realized as open-access codes linked to existing DFT-based software; (2) Establishing the relationship between the structure of adducts and SWCNTs and intrinsic excitonic and spin properties of defect states for guiding novel synthetic strategies and experimental probes of chemically functionalized SWCNTs as near-IR emitting materials; (3) Generating virtual libraries of hypothetical functionalized SWCNTs for virtual screening of their chemical structures and optical properties, leveraging new functionalities of SWCNTs; (4) Offering a unique experience for NDSU graduate students that prepared them for future scientific careers related to materials modeling and big data processing. These results were summarized in 12 published journal papers and 3 recently submitted papers. One of a key finding is that the position of defect sites on the SWCNT surface primarily drives the emission redshift (up to 100 meV), while the polarity of the defect-inducing molecules has a much smaller effect (~10 meV). However, the electron-donating or withdrawing properties of a molecule influence selecting reactivity of defect sites. These insights important for optimizing synthetic protocols for desired emissions in SWCNTs. We also revealed that the interaction between two defects at various positions on the SWCNT enhances the redshift and optical activity of states, favoring strong near-IR emission. This suggests that manipulations in defect concentrations is a promising strategy for controlling efficient emission. Mostly important, the defect position was found controllable by the spin states of photoexcited intermediates: Excited aromatic molecules form ortho defects with SWCNTs at their singlet states in the presence of oxygen, while oxygen-free conditions favor para defects via the triplet-state mechanism. Additionally, a heat-activated [2+2] cycloaddition reaction facilitates divalent defect formation with fewer bonding positions that narrows emission bands. These groundbreaking findings have been experimentally validated and significantly advance our understanding of defect chemistry in SWCNTs. Using a novel encoding technique and 3D-MoRSE descriptors, we developed highly accurate ML/QSAR models to predict both the 3D structure and optical properties of SWCNTs with chemical defects. This model enabled the creation of a virtual library of 125,556 structures, providing new insights into the relationship between SWCNT-defect structure and emission.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

Developing Digital Twin Visualizations: A Methodology and Case Study on Chemical Separation Processing

As advances in digital engineering continue to push the technological boundaries, digital twin (DT) visualizations for diagnostics and safeguards advancement become much more feasible and practical. DTs generate large and complex data streams that require effective user interfaces to provide monitoring and diagnostic capabilities. Unfortunately, while these frameworks exist, there is not much research on the systematic documentation of human–computer interaction (HCI) for DT visualization. This work presents a dual-mode visualization methodology (two dimensional [2D] graphical user interface dashboard and 3D mixed reality) designed to support diagnostic tasks in DT systems and building on a validated framework and applying established HCI principles. The methodology is demonstrated through a case study of aqueous processing at Idaho National Laboratory, using experimental data from the chemical solvent extraction runs. Our interfaces display real-time alerts and monitoring to inform users of safeguards anomalies. The interfaces use immersive 3D mixed-reality visualization for further system and experiment investigation. This work demonstrates how the systematic application of HCI principles can inform DT visualization design for diagnostic and safeguards applications. While formal user evaluation studies remain as future work, this paper documents the systematic design methodology and demonstrates a proof-of-concept implementation.

3D visualization↗
DOE OSTI2024

Exploring the Neutron Substructure with Advanced Polarized Helium-3 Targets (Or: How I Learned to Stop Worrying and Love Spectroscopy)

As we seek to understand the smallest, physical aspects of our universe, we cannot simply rely on our senses to probe the world around us as we did in the past. The smallest physical elements of our universe behave in strange, probabilistic ways and are completely invisible to the naked eye/ear/etc. So, we design clever experiments (such as scattering experiments) to probe these minute realms. Then, just as with the larger, observable world, we devise models and equations to describe what we think is happening. Due to the nature of the physical universe at the quantum scale and with the aid of symmetries such as Lorentz invariance, we can write down equations that describe the scattering, but the expressions contain functions, which we call ?form factors? and ?structure functions?, that we cannot compute from first principles. We can, however, formulate models that make predictions for these functions. By comparing our predictions with the observed data, we can gain insight into the validity of our models and thus a better physical understanding of what is happening at these minuscule scales. Studying the constituents inside of the nucleus of an atom adds another layer of difficulty if we can?t remove those components from the nucleus. This is the case with the neutron. When not bound in the nucleus with protons and other neutrons, the neutron will decay into a proton after about 15 minutes. So, we?re forced to study the neutron while it is still bound in the nucleus of an atom such as helium-3 (3He). For the last 1,000 years (rounding up), our group has developed high quality, polarized 3He targets made of an aluminosilicate glass. These targets are made in order to perform experiments at Jefferson Lab (JLab), experiments which let us determine the form factors and structure functions of the neutron by scattering polarized electrons from polarized neutrons (or rather polarized 3He). The specific experiments reported on in this thesis push the bounds of our understanding of the internal structure of the neutron. Good science is often about pushing experimental techniques to a new level. Toward that goal we study our polarized 3He targets both to advance the technology and to choose the best ones for our experiments. We do this using a process called nuclear magnetic resonance (NMR) to gauge the maximum polarization of a target and how fast the polarization decays with time. While these tests primarily provide us information that make analysis of our experimental scattering data possible, they also let us determine whether or not a target-cell is useful or even, dare I say, of spectacular quality. Our latest targets utilize a novel convection design allowing 3He to be polarized and quickly moved in front of the electron-beam, making it possible to use larger targets with higher electron-beam currents than ever before. This means more electrons scatter and we get more data. And by studying our targets in detail prior to using them in our experiments, we have found techniques to take effects which could have been detrimental to target quality and turn them to our advantage! It?s a real case of making lemonade out of lemons. We also use laser spectroscopy to study the absorption lines of alkali-metals in the target (potassium and rubidium, specifically). We add these alkali-metals to our target to facilitate polarizing the 3He. We can use the measurement of these pressure broadened absorption lines to determine the 3He density inside of the target with great precision. Historically, we understood the width of these lines would be dependent on the temperature of the target. Specifically, if I raise the temperature, the width should get bigger. I found that was not the case, which was very confusing at first, though very exciting now that I realize the data are self-consistent and suggestive of unexpected behavior. This thesis details the development of high quality, glass, polarized 3He targets for the 2020 An 1 /dn 2 and 2023 Gn E experiments, which utilized the first 3He convection targets and broke records in target quality. This thesis also covers the initial development of metal windows for the next-generation of 3He target-cells. Finally, this thesis documents the temperature dependence of the width of potassium (K) and rubidium (Rb) absorption lines as measured with laser spectroscopy.

Jantzi, Christopher↗
DOE OSTI2025

Exascale granular microstructure reconstruction in 3D volumes of arbitrary geometries with generative learning

Reconstructing 3D granular microstructures within volumes of arbitrary geometries from limited 2D image data is crucial for predicting the material properties, as well as performances of structural components accounting for material microstructural effects. We present a novel generative learning framework that enables exascale reconstruction of granular microstructures within complex 3D geometric volumes. Building upon existing transfer learning techniques using pre-trained convolutional neural networks (CNN), we introduce several key innovations to overcome the difficulties inherent in arbitrary geometries. Our framework incorporates periodic boundary conditions using circular padding techniques, ensuring continuity and representativeness of the reconstructed microstructures. We also introduce a novel seamless transition reconstruction (STR) method that creates statistically equivalent transition zones to integrate multiple pre-existing 3D microstructure volumes. Based on STR, we propose a cost-effective strategy for reconstructing microstructures within complex geometric volumes, minimizing computational waste. Validation through numerical experiments using kinetic Monte Carlo simulations demonstrates accurate reproduction of grain statistics, including grain size distributions and morphology. A case study involving the reconstruction of a 4-blade propeller microstructure illustrates the method’s capability to efficiently handle complex geometries. In conclusion, the proposed framework significantly reduces computational demands while maintaining high reconstruction quality, paving the way for scalable microstructure reconstruction in materials design and analysis.

36 MATERIALS SCIENCE↗
DOE OSTI2025

SAGIPS: a physics-inspired scalable asynchronous generative inverse-problem solver

Abstract Solving large-scale inverse problems using deep-learning algorithms have become an essential part of modern research and industrial applications. The complexity of the underlying inverse problem may require the utilization of high performance computing systems which poses a challenge on the algorithmic design of the inverse problem solver. Most deep learning algorithms require, due to their design, custom parallelization techniques in order to be resource efficient while showing a reasonable convergence. In this paper we introduce a S calable A synchronous G enerative I nverse P roblem S olver (SAGIPS) on high-performance computing systems. We present a workflow that utilizes an asynchronous ring-allreduce algorithm to transfer the gradients of the generator network across multiple GPUs. Experiments with a scientific proxy application demonstrate that SAGIPS shows near linear weak scaling, together with a convergence quality that is comparable to traditional methods. The approach presented here allows leveraging Generative Adverserial Network across multiple GPUs, promising advancements in solving complex inverse problems at scale.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

ATTNChecker: Highly-Optimized Fault Tolerant Attention for Large Language Model Training

Large Language Models (LLMs) have demonstrated remarkable performance in various natural language processing tasks. However, the training of these models is computationally intensive and susceptible to faults, particularly in the attention mechanism, which is a critical component of transformer-based LLMs. In this paper, we investigate the impact of faults on LLM training, focusing on INF, NaN, and near-INF values in the computation results with systematic fault injection experiments. We observe the propagation patterns of these errors, which can trigger non-trainable states in the model and disrupt training, forcing the procedure to load from checkpoints. To mitigate the impact of these faults, we propose ATTNChecker, the first Algorithm-Based Fault Tolerance (ABFT) technique tailored for the attention mechanism in LLMs. ATTNChecker is designed based on fault propagation patterns of LLM and incorporates performance optimization to adapt to both system reliability and model vulnerability while providing lightweight protection for fast LLM training. Evaluations on four LLMs show that ATTNChecker on average incurs on average 7% overhead on training while detecting and correcting all extreme errors. Compared with the state-of-the-art checkpoint/restore approach, ATTNChecker reduces recovery overhead by up to 49×.

Liang, Yuhang [University of Alabama - Birmingham]↗
DOE OSTI2025

Optimizing the design and operation of water networks: Two decomposition approaches

We consider the design and operation of water networks simultaneously. Water network problems can be divided into two categories: the design problem and the operation problem. The design problem involves determining the appropriate pipe sizing and placements of pump stations, while the operation problem involves scheduling pump stations over multiple time periods to account for changes in supply and demand. Our focus is on networks that involve water co-produced with oil and gas. While solving the optimization formulation for such networks, we found that obtaining a primal (feasible) solution is more challenging than obtaining dual bounds using off-the-shelf mixed-integer nonlinear programming solvers. Therefore, we propose two methods to obtain good primal solutions. One method involves a decomposition framework that utilizes a convex reformulation, while the other is based on time decomposition. To test our proposed methods, we conduct computational experiments on a network derived from the PARETO case study.

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

Maximized Information Gain of Next Generation Pulsed Power Using Optimized Design of Z-Machine Experiments

This project develops a Bayesian optimization approach to extracting insights from Z Machine experimental data to determine if and how these insights can be used to extrapolate to a larger facility. The primary goal is to address the scientific challenge of informing how confidently experimental conditions can be predicted on a next generation facility, the design of which requires the reliable extrapolation of current high energy density technologies to regimes yet unobserved, except by costly high-fidelity computational models. Maximizing the use of presently available data and understanding how it informs future endeavors is critically important to enable transformative pulsed power and the science of extreme conditions. We explore a Bayesian optimization approach to experimental design which combines information theory, experimental data, and computational modeling to explore how information gain can be maximized.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Screening and qualification methodology for SiC end plug processing methods

Deployment of SiC-ceramic-based fuel cladding for light water reactors requires a hermetic end plug–to–cladding joint that can withstand neutron irradiation during normal operation and maintain integrity during design-basis accidents. Reactor experiments have shown that some SiC composite tubes with SiC end plugs can retain hermeticity after irradiation. However, achieving consistent joint performance under irradiation remains a key challenge. Resolving this issue is essential to enable integral irradiation testing and to demonstrate fuel integrity under commercial-reactor irradiation conditions. This report aims to: (1) provide guidance for designing radiation-tolerant end plug joints for SiC cladding; (2) demonstrate experimental methods to detect processing defects that are unstable under neutron irradiation at light-water-reactor-relevant temperatures and doses; and (3) outline a step-by-step approach for designing and conducting reactor experiments to screen joining methods. The resulting data will be used to improve joint processing and to define critical defect types and sizes that must be detected and eliminated through non-destructive evaluation for quality assurance. Based on prior irradiation experiments at the High Flux Isotope Reactor, differential swelling among the cladding, bonding layer, and end plug was identified as an underlying mechanism for irradiation-induced joint degradation. Accordingly, this effect must be considered in the design of radiation-tolerant joining techniques. In this work, miniature SiC end plug joint specimens irradiated during the previous project were analyzed using X-ray computed tomography to characterize the joint microstructure. Digital volume correlation of the tomography data quantified radiation-induced microstructural changes and enabled evaluation of defect-related risks. Finally, ongoing neutron irradiation efforts using larger specimen volumes are presented. These efforts aim to statistically assess joint performance and to build a microstructure–performance (e.g., leak-tightness) dataset to inform processing improvements and quality control.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Preventing Loss of Selectivity during the Oxidative Dehydrogenation of Propane over Supported Vanadium Catalysts

Supported vanadium materials are promising catalysts for the oxidative dehydrogenation of propane to propylene (ODHP), but a lack of mechanistic understanding limits the rational design of catalysts with improved propylene selectivity. Adding Ta to V/SiO 2 increases the propylene selectivity, as well as the activity, leading to superior performance compared to state-of-the-art boron-based systems. In this contribution, we utilize this surprising promotional effect of Ta to elucidate key elements of the mechanistic cycle. Through a combination of characterization techniques, computational modeling, and kinetic experiments, we show that the catalytic cycle over V/SiO 2 likely involves the formation of an isopropyl alcohol intermediate, the fate of which is in kinetic competition between subsequent dehydration to propylene or further oxidation. Furthermore, we show that the relatively facile propylene overoxidation observed for these materials occurs via the epoxidation of propylene by a proposed peroxovanadium intermediate, rather than the abstraction of propylene’s allylic C–H bond as previously assumed. Using these key mechanistic features, we rationalize the enhanced selectivity and activity of Ta promotion. In conclusion, our mechanistic framework offers avenues for future catalyst development to improve supported vanadium materials for ODHP.

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

Electromagnetic Energy-Assisted Thermal Conversion of Fossil-Based Hydrocarbons to Low-Cost Hydrogen

The goal of this project was to develop and optimize catalysts for methane decomposition, particularly focusing on regeneration via an electromagnetic energy-assisted mechanism, to produce hydrogen more cost-effectively compared to electrolysis routes. To achieve this goal, the project pursued several key objectives. The project began with the preparation and testing of various catalysts. A nickel-silica based catalyst was identified as the most promising material for the pyrolysis of methane into carbon and hydrogen. Kinetic parameters for methane decomposition were determined, aiding in computational modeling efforts. Structured catalysts were investigated, highlighting the need for frequent cleaning or regeneration to maintain performance, with methane conversion rates exceeding 70% in tube furnace tests. Computational fluid dynamics modeling was employed to optimize reactor designs and electrode angles, leading the project team to propose a multi-compartment thermal conversion system for larger setups. This modeling work was important in understanding reaction characteristics, carbon deposition rates, and temperature profiles under various conditions. A bench-scale reactor system was assembled to evaluate catalyst regeneration using electromagnetic energy-assisted mechanisms. Experiments demonstrated the potential for carbon removal, though further optimization is needed. The carbon produced from the methane conversion process was evaluated for potential use in lithium-ion battery electrodes. The carbon exhibited properties similar to commercially available high-purity multi-walled carbon nanotubes and nanofibers, with a carbon content greater than 95%. Coin cell batteries assembled with this carbon showed that lower replacement levels (10% to 33%) outperformed the control group, improving specific capacity density and stability. However, higher replacement levels (100%) demonstrated poorer performance, suggesting that excessive carbon substitution negatively impacts battery performance. These findings indicate the potential marketability of the produced carbon as a component in lithium-ion batteries, though further testing is necessary to confirm long-term advantages and disadvantages associated with the use of the carbon product. These results however justified further technoeconomic assessments to determine if the process can provide low-cost hydrogen. The economic feasibility and technical performance of methane decomposition for hydrogen production were assessed, focusing on three plant configurations: 100E (electrically heated), 100C (combustion heated using produced hydrogen), and PE-Hybrid (a combination of pyrolysis (indicating decomposition) and electrolysis). The Levelized Cost of Hydrogen (LCOH) for the pyrolysis configurations was found to be approximately 25% lower than that of electrolysis. The 100E configuration had the lowest LCOH at $\$$3.12/kg. Including carbon product sales significantly improved the economics, with the 100C configuration achieving a negative LCOH of -$\$$0.35/kg. The PE-Hybrid configuration was not economically advantageous compared to pure pyrolysis plants due to its complexity and additional equipment requirements. Ultimately, methane pyrolysis presents a viable method for near carbon dioxide-free hydrogen production, with significant economic advantages over electrolysis, especially when considering the sale of carbon byproducts. The 100E and 100C configurations showed the most promise, with the choice between them ultimately depending on the prices of power and natural gas. In conclusion, this technology has the potential to lower hydrogen production costs by leveraging the methane decomposition process with the sale of valuable carbon byproducts. By optimizing catalyst performance and integrating electromagnetic energy-assisted regeneration, the process can achieve higher efficiency and economic viability, making it a competitive alternative to traditional hydrogen production methods.

08 HYDROGEN↗
DOE OSTI2024

Study of the heat pipe vapor flow using the MOOSE application Sockeye and the CFD code Nek5000

Heat pipes are efficient heat transfer devices used in various applications, including nuclear microreactors. Heat pipe-cooled microreactors offer advantages in size and cost, which can significantly accelerate their deployment and adoption. Developing accurate models of heat pipes is crucial for a heat pipe-cooled microreactor design and operation. In this work a verification and validation of the heat pipe application Sockeye is presented. Results are compared against RANS simulations, performed with the code Nek5000, and the SPHERE experiment, conducted at INL.

42 ENGINEERING↗
DOE OSTI2025

The design space of E(3)-equivariant atom-centred interatomic potentials

Abstract Molecular dynamics simulation is an important tool in computational materials science and chemistry, and in the past decade it has been revolutionized by machine learning. This rapid progress in machine learning interatomic potentials has produced a number of new architectures in just the past few years. Particularly notable among these are the atomic cluster expansion, which unified many of the earlier ideas around atom-density-based descriptors, and Neural Equivariant Interatomic Potentials (NequIP), a message-passing neural network with equivariant features that exhibited state-of-the-art accuracy at the time. Here we construct a mathematical framework that unifies these models: atomic cluster expansion is extended and recast as one layer of a multi-layer architecture, while the linearized version of NequIP is understood as a particular sparsification of a much larger polynomial model. Our framework also provides a practical tool for systematically probing different choices in this unified design space. An ablation study of NequIP, via a set of experiments looking at in- and out-of-domain accuracy and smooth extrapolation very far from the training data, sheds some light on which design choices are critical to achieving high accuracy. A much-simplified version of NequIP, which we call BOTnet (for body-ordered tensor network), has an interpretable architecture and maintains its accuracy on benchmark datasets.

Computer Science↗
DOE OSTI2025

Computational materials reliability assessment of hydrogen fueled gas turbine power generation engines

The use of blended fuel sources in land based gas turbine engines drives variations in the resulting operational profile (temperatures and pressures) which can impact engine reliability. Furthermore, variability in the manufacture of components affects the resulting microstructure which directly impacts material performance and reliability. Currently, data-driven models are typically used for maintaining and inspecting fleets of engines. Without explicitly capturing material and operational sources of variability conservatism must be used in developing component-level reliability models. Therefore, there exists an opportunity to use information from materials-scale physics models to better inform reliability modeling and reduce conservatism; the impact is more cost-efficient operation and maintenance of current and future fleets. Specifically, this work establishes a computational framework for evaluating the probabilistic high temperature creep performance of hot-section Ni-based superalloys where uncertainty comes from both microstructural and operational variability. A novel high-fidelity physics model which phenomenologically captures grain-boundary sensitive phenomena has been established. A probabilistic calibration procedure was used to calibrate the model and capture uncertainty in the parameterized model coefficients. A design of experiments methodology was established for identifying informative microstructural digital representations for suitable for forward model evaluation. Results show that training a machine-learning surrogate using this design criteria outperforms random selection of microstructural representations. Finally, two surrogate models were developed: (1) a deterministic surrogate model which predicts the local field response given microstructure, constitutive model parameters, and operating conditions (stress, temperature) and (2) a probabilistic model, where uncertainty comes from constitutive law uncertainty, built using denoising diffusion probabilistic models which samples responses given (1) microstructure and (2) operating conditions. These surrogate models enable partner Siemens Energy to rapidly perform UQ analysis specific to creep deformation across a range of microstructures and operating conditions. The impact is that these ML and physics codes can be used to establish more advanced reliability models for the inspection, servicing, and maintenance of land based gas turbine engines.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Identifying Bayesian optimal experiments for uncertain biochemical pathway models

Abstract Pharmacodynamic (PD) models are mathematical models of cellular reaction networks that include drug mechanisms of action. These models are useful for studying predictive therapeutic outcomes of novel drug therapies in silico. However, PD models are known to possess significant uncertainty with respect to constituent parameter data, leading to uncertainty in the model predictions. Furthermore, experimental data to calibrate these models is often limited or unavailable for novel pathways. In this study, we present a Bayesian optimal experimental design approach for improving PD model prediction accuracy. We then apply our method using simulated experimental data to account for uncertainty in hypothetical laboratory measurements. This leads to a probabilistic prediction of drug performance and a quantitative measure of which prospective laboratory experiment will optimally reduce prediction uncertainty in the PD model. The methods proposed here provide a way forward for uncertainty quantification and guided experimental design for models of novel biological pathways.

97 MATHEMATICS AND COMPUTING↗