Search NASA⌕ Search

SEARCH · Search NASA

Results for “approximation methods”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

256 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

End-To-End Decentralized Transmission Line Protection in IBR-Dominated Weak Grids Using Interpretable Data-Driven Methods

Traditional transmission line protection relies on predictable synchronous-based fault signatures, which frequently fail under the non-standard, current-limited fault characteristics of Inverter-Based Resources (IBRs). This study investigates how to achieve secure, communication-free fault isolation in IBR-dominated weak grids without relying on opaque, computationally heavy "black-box" machine learning algorithms. To address this, we propose a novel, standalone, and inherently interpretable data-driven protection framework. Unlike centralized methods requiring multi-terminal communication, this decentralized approach relies solely on local measurements using a hierarchical linear-kernel Support Vector Machine (SVM). The methodology decomposes the protection task into four sequential stages that mimic traditional protection elements: fault detection and fault direction identification, fault type classification, zone classification, and location estimation. This multi-stage architecture allows for specialized feature engineering at each stage, combining high computational efficiency with logic traceability. The framework's end-to-end performance was validated via C-code and PSCAD/EMTDC co-simulation, utilizing a real-world utility network and an OEM black-box IBR model. The proposed relay achieves 97.2% overall accuracy and provides a reliable trip decision within a 2.5-cycle window. The results confirm 100% accuracy in fundamental fault detection, reliable zone selectivity across low to moderate fault resistances, and robust security against non-fault transients, proving its immediate viability for integration into commercial numerical relays.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

A Data-Driven Method for Modeling Creep-Fatigue Stress- Strain Behavior Using Neural ODEs

In this paper, we introduce a data-driven machine learning approach for modeling one-dimensional stress–strain behavior under cyclic loading, utilizing experimental data from the nickel-based Alloy 617. The study employs uniaxial creep–fatigue test data acquired under various loading histories and compares two distinct neural network-based ODE models. The first model, known as the black-box model, comprehensively describes the strain–stress relationship using a Neural ODE equation. To interpret this black-box model, we apply the Sparse Identification of Nonlinear Dynamical Systems (SINDy) technique, transforming the black-box model into an equation-based model using symbolic regression. The second model, the Neural flow rule model, incorporates Hooke’s Law for the linear elastic component, with the nonlinear part characterized by a Neural ODE. Both models are trained with experimental data to accurately reflect the observed stress–strain behavior. We conduct a detailed comparison with the standard Chaboche model, which includes three back stresses. Our results demonstrate that the neural network-based ODE models precisely capture the experimental creep–fatigue mechanical behavior, exceeding the standard Chaboche model’s accuracy. Furthermore, an interpretable model derived from the black-box neural ODE model through symbolic regression achieves accuracy comparable to the Chaboche model, enhancing its interpretability. The results highlight the potential of neural network-based ODE models to depict complex creep–fatigue behavior, eliminating the necessity for experts to define a specific, material-focused model form.

creep-fatigue↗
DOE OSTI2024

A First-Principles Study of the Structural and Thermo-Mechanical Properties of Tungsten-Based Plasma-Facing Materials

Tungsten (W) and tungsten alloys are being considered as leading candidates for structural and functional materials in future fusion energy devices. The most attractive properties of tungsten for the design of magnetic and inertial fusion energy reactors are its high melting point, high thermal conductivity, low sputtering yield, and low long-term disposal radioactive footprint. Despite these relevant features, there is a lack of understanding of how the structural and mechanical properties of W-based alloys are affected by the temperature in fusion power plants. In this work, we present a study on the thermo-mechanical properties of five W-based plasma-facing materials. First-principles density functional theory (DFT) calculations are combined with the quasi-harmonic approximation (QHA) theory to investigate the electronic, structural, mechanical, and thermal properties of these W-based alloys as a function of temperature. The coefficient of thermal expansion, temperature-dependent elastic constants, and several elastic parameters, including bulk and Young’s modulus, are calculated. Our work advances the understanding of the structural and thermo-mechanical behavior of W-based materials, thus providing insights into the design and selection of candidate plasma-facing materials in fusion energy devices.

42 ENGINEERING↗
DOE OSTI2026

Reduced-order CFD modeling of cryogenic hydrogen isotope extrusion for pellet fueling

This study presents a reduced-order model (ROM) for computational fluid dynamics (CFD) simulations of cryogenic hydrogen isotope extrusions, focusing on protium (H₂) and deuterium (D₂) piston extruders. Using a 2D axisymmetric ROM in ANSYS-Polyflow, significant computational savings were achieved (runtime reduced from 9∼24 h to 3∼5 min), with extrusion force discrepancies between the 2D ROM and 3D models being on the order of 1%. Parametric studies identified optimal cutoff shear rates in the viscosity model (0.01/s for H₂ and 0.001/s for D₂), providing recommendations for future simulations. Finally, a comprehensive comparison of ROM results with experimental data was performed across varying geometries, cryogenic materials, temperatures, extrusion lengths, and piston velocities. Predictions at low extrusion temperatures met the objective of providing quick and efficient solutions with an acceptable extrusion force error of approximately 10% or less, validating the effectiveness of the 2D ROM approach. However, at high temperatures closer to the triple point, extrusion force error grows, which necessitates developing an improved model that accounts for temperature effects, e.g. melting. Nevertheless, the findings still represent a significant improvement in efficiency of CFD modeling of cryogenic hydrogenic extrusion. The ROM framework can also be extended to tritium (T2) and screw extruders, which will ultimately provide a fast and effective tool for optimizing pellet injector design for ITER and future reactor systems.

Fan, Joy [ORNL] (ORCID:0000000229751735)↗
DOE OSTI2026

Collision Tracking in OpenMC: Methods and Applications in Neutron Noise, Neutron Imaging, Time-of-Flight, and Multiplicity Counting

We present the development and application of a collision tracking feature within the OpenMC Monte Carlo particle transport code, designed for diverse applications such as neutron spectroscopy, scatter camera system, neutron noise, and multiplicity counting simulations. This feature enables the tracking of individual particle collisions, with potential applications in nuclear nonproliferation, reactor physics, and nuclear security. Additionally, the feature holds potential for the calibration of neutron detectors, specifically in converting light output into energy deposited within the detectors. The implementation consists of a set of filters—such as reaction type, energy, cell, and material—that constrain the set of collisions that are tracked, extensions to the Python API to enable simple input specification, and support for writing either OpenMC’s native HDF5-based format or the Monte Carlo particle list format. This feature was added to the official OpenMC release in version 0.15.3. In this work, the feature will be applied to showcase scenarios such as time-of-flight simulations, scatter-camera imaging for neutron source localization, neutron-noise analysis to extract integral kinetic parameters such as the prompt decay constant α, and multiplicity counting to estimate the mass of special nuclear materials. Ultimately, this feature aims to expand the application scope of open-source Monte Carlo particle transport codes such as OpenMC.

Monte Carlo code↗
DOE OSTI2026

A parametric study of slow dynamic nonlinear elasticity with comparisons to models

Several phenomenological models that aspire to quantitative description of anomalous nonlinear mesoscopic elasticity are reviewed and compared with laboratory measurements. This class of nonlinearity, best known perhaps for slow dynamics and aging, is seen widely in imperfectly consolidated granular solids but is not well understood. Typical slow dynamic tests show that a modest conditioning oscillatory "pump" strain depresses material stiffness, which then recovers like the logarithm of time after conditioning ceases. Several phenomenological models based on physical arguments have been proposed that predict the material stiffness response to arbitrary pump strain histories during conditioning and recovery. Approximate closed form and numerical solutions to the models are presented that predict the quantitative influence of three key pump parameters: the pump's strain amplitude, the pump's strain rate, and the pump’s duration. Laboratory measurements on Berea sandstone, concrete and a confined single aluminum bead find that slow dynamic responses are linear in pump strain and independent of pump frequency. Measurements also show that, after pump-off, stiffness recovers over times far longer than the pump duration. These observations and others are compared to model predictions. One of the considered models, based on a picture of fast brittle damage and slow healing, successfully matches all these behaviors.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Next-Generation Energy Technologies for Connected and Automated On-Road Vehicles (NEXTCAR) - Predictive Data-Driven Vehicle Dynamics and Powertrain Control: from ECU to the Cloud (Final Scientific/Technical Report)

This project developed and demonstrated a predictive, data-driven vehicle control system designed to improve energy efficiency and driving performance. The team created intelligent self-driving car technology that optimizes fuel and electricity use by proactively planning vehicle actions. By combining Level 4 autonomous driving capabilities with vehicle-to-everything (V2X) connectivity, the system enables vehicles to adjust speed and change lanes in response to traffic signals, surrounding vehicles, and road conditions, reducing unnecessary stops and delays. In testing, the system improved vehicle fuel economy by more than 30% and reduced travel time by approximately 10%, compared to a conventional adaptive cruise control baseline. These results demonstrate the technical effectiveness of using predictive, V2X-enabled strategies, such as traffic light timing and surrounding traffic awareness, to inform real-time vehicle powertrain control and driving behavior. Additionally, a supporting cloud platform was developed to provide dispatch and route recommendations as well as to log vehicle data, demonstrating the economic feasibility of this approach at the fleet level. By optimizing dispatching and routing operations, this technology enables electric fleet operators to use their vehicles more efficiently and reduce reliance on diesel backups, lowering both operating costs and energy consumption. Overall, this project’s technology advances the future of clean, energy-efficient transportation, enabling vehicles and fleets to reduce energy waste, cut costs, and lower emissions through intelligent automation and connectivity.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

Variational Quantum Circuits to Prepare Low Energy Symmetry States

We explore how to build quantum circuits that compute the lowest energy state corresponding to a given Hamiltonian within a symmetry subspace by explicitly encoding it into the circuit. We create an explicit unitary and a variationally trained unitary that maps any vector output by ansatz A(α → ) from a defined subspace to a vector in the symmetry space. The parameters are trained varitionally to minimize the energy, thus keeping the output within the labelled symmetry value. The method was tested for a spin XXZ Hamiltonian using rotation and reflection symmetry and H 2 Hamiltonian within S z = 0 subspace using S 2 symmetry. We have found the variationally trained unitary gives good results with very low depth circuits and can thus be used to prepare symmetry states within near term quantum computers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Gear Test Assembly - Experimental Testing and Analysis of Gears and Bearings - FY2026

The Gear Test Assembly (GTA) is an experimental test apparatus built and installed in the Mechanisms and Engineering Test Loop (METL) at Argonne National Laboratory (ANL). While designed to accommodate a variety of components intended for use in advanced compact in-vessel transfer machine, the focus of GTA been the testing of large radial spur gears. The performance of the components used in GTA motivates the choices made in the design, components, and material choices for the forthcoming METL test articles; the Gripper Test Assembly (GrTA), and the Bearing Test Assembly (BTA). To date, GTA has completed nine experimental testing campaigns and achieved over 21 million revolutions under prototypic sodium fast reactor loads and operating conditions. This report provides an overview and results for the two most recent experiments Campaign #8 and Campaign #9. `Campaign #8 was the first campaign to make use of a set of nickel alloy radial spur gears, made from Hastelloy C-22HS, and combined with heat-treated tapered roller bearings. The bearings and torque profile used allowed Campaign #8 to be directly compared with Campaign #1 where Inconel 718 gears were used instead. Campaign #8 attained 1,058,880 shaft revolutions within 23% of the 1,314,855 revolutions achieved in Campaign #1. These results demonstrate that the heat-treated tapered roller bearings show consistent increased lifetimes, approximately three times longer than for non heat-treated taper rollers, and that the particular material of the gear does not play a substantial role in performance. Similar to Campaign #7 where bearing failure resulted in fragmentation and transport of bearing material into the gear teeth led to tooth damage in the Inconel gears, a similar process led to the damage of the Hastelloy gears. Subsequent NDE analysis showed that the damage was only limited to the tooth surface and does not extend into the gear interiors. Based on the NDE reports, the Hastelloy gears can be reconditioned and put back into service for future campaigns as was done for the Inconel gears. Reconditioned Inconel 718 gears damaged in Campaign #7 were returned to service in Campaign #9 and paired with non heat-treated tapered roller bearings. In Campaign #9 a novel set of operating conditions were implemented whereby the gears and bearings were cycled between hot, submerged refueling conditions where fuel handling operations occurred and cold, dry, and inert conditions where the sodium was drained and the vessel was allowed to cool simulating long, ex-vessel storage of the fuel handling machine. A total of eight of these operating cycles were achieved before a bearing failure occurred after the accumulation of 475,692 shaft revolutions were achieved. During Campaign #9, one of the drive motors was damaged and over the course of a 240 day standby period, the entire GTA drive system was upgraded from a 480VAC system to a 240VAC system that provided the opportunity to consolidate wall-space in the METL facility and provide more space for additional test articles being installed in the coming year. The over 475,000 shaft revolutions achieved in Campaign #9 eclipses the previous record for lifetime of non heat-treated tapered roller bearings held by Campaign #3 of 392,000 shaft revolutions. Further analysis and future campaigns hope to shed light on the cause of this approximately 17% increase in lifetime.

42 ENGINEERING↗
DOE OSTI2025

230Th/234U radiochronometry for uranium materials by alpha spectrometry for nuclear forensics analysis

We have developed an alpha spectrometry method for 230Th/234U radiochronometry to determine the model separation age for uranium materials. This method offers an alternative to the more commonly used, but significantly more resource intensive multi-collector inductively coupled plasma mass spectrometry (MC-ICP-MS) while providing similar accuracy and precision. The development included radiochemical separation of uranium and thorium, calibration of the 232U tracer and the alpha spectrometry measurement. The method was validated by analyzing a certified reference material, CRM 125-A, which is a uranium oxide pellet assay, isotopic, and radiochronometric standard. The results were compared to those obtained by our routine MC-ICP-MS 230Th/234U radiochronometry analysis method to assess the performance parameters for alpha spectrometry compared to the state-of-the-art technique.

Macsik, Zsuzsanna [Los Alamos National Laboratory ↗
DOE OSTI2026

Polynomial-time preparation of low-temperature Gibbs states for two-dimensional toric code

In this work, we propose a polynomial-time algorithm for preparing the Gibbs state of the two-dimensional toric code Hamiltonian at any temperature, starting from any initial state, significantly improving upon prior estimates that suggested exponential scaling with inverse temperature. We prove that fast mixing at low temperature for the two-dimensional toric code can be achieved by augmenting local jump operators with simple global jump operators, which enable efficient transitions between logical sectors. To establish tight lower bounds on the spectral gap, we introduce a new reduction method that eventually maps the problem to estimating the spectral gap of a perturbed graph Laplacian on a stair graph. Our proof also shows that the Lindblad dynamics with a digitally implemented low-temperature local Davies generator is able to efficiently drive the quantum state toward the ground state manifold.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Design Trade-Offs in Composite Fuel Cell Membranes: Effects of Reinforcement and Chemical Additives

Perfluorosulfonic acid (PFSA) membranes are critical components in proton exchange membrane fuel cells, where performance depends on balancing ionic conductivity, mechanical durability, and chemical stability. This study characterizes a composite membrane (NC700) featuring PFSA-impregnated expanded polytetrafluoroethylene (ePTFE) reinforcement and cerium-based radical scavengers, benchmarked against unreinforced NR211. Complementary techniques, including electron microscopy, X-ray scattering, infrared spectroscopy, thermogravimetric analysis, and dynamic mechanical analysis, identify the structural and compositional strategies employed in NC700. Water sorption isotherms reveal lower water uptake for NC700 across all conditions, attributed to reinforcement and cerium incorporation. Reinforcement reduces in-plane swelling from 11% to 2.1% at 90% RH, confirming strong swelling anisotropy, while maintaining mechanical properties at elevated temperatures. While the ionic conductivity of NC700 is approximately 10% lower than that of NR211, the reduced thickness yields a 40% decrease in calculated area-specific resistance, suggesting the composite architecture can favorably shift the conductivity-stability trade-off. The composite structure also reduces gas permeability, indicating potential for improved separator function alongside favorable transport properties. Systematic deconvolution of reinforcement and additive contributions shows that conductivity losses from cerium incorporation are largely offset by gains from the lower equivalent-weight polymer, providing quantitative relationships that may guide composite membrane design for fuel cells and other electrochemical applications.

25 ENERGY STORAGE↗
DOE OSTI2026

A transport-based framework for solidification cracking in Ni-based superalloys processed by laser powder bed fusion

Solidification cracking remains a persistent barrier to laser powder bed fusion (LPBF) processing of solid-solution-strengthened (SSS) Ni-based superalloys and is commonly attributed to carbide formation, liquation-type failure, or elemental segregation. In this work, the cracking behavior of three SSS Ni-based superalloys—Inconel 625 (625), Inconel 617 (617), and Haynes 230 (230)—is examined under comparable LPBF conditions to evaluate the origins of their cracking susceptibility. Carbides were present in both 625 and 230, yet cracking behavior differed significantly. MC-type carbides in 625 remain discrete and preserve liquid connectivity, whereas Cr-rich M23C6 carbides in 230 form at cellular triple junctions that bottleneck the interdendritic liquid network, influencing cracking through liquid connectivity rather than as an independent cause. Alloy 617 exhibited anomalous segregation without resolvable secondary phases yet cracked mildly, indicating that extensive carbide formation is not a prerequisite for cracking. To rationalize these observations, a transport-limited framework is proposed in which cracking occurs when terminal-liquid feeding cannot accommodate solidification-induced strain; among the mechanisms evaluated, this framework is most consistent with the physical constraints imposed by LPBF solidification. Transport-based crack-susceptibility indexes incorporating permeability, viscosity, and solidification kinetics distinguish the alloy hierarchy; permeability ratios Krat increased from approximately 7 (625) to 11.5 (617) and 14.5 (230), with corresponding increases in the μ-weighted mushy zone risk index (CSIMZRM). These results are consistent with transport-limited liquid feeding as a rate-limiting contributor to solidification cracking in LPBF-processed SSS Ni-based superalloys, providing a physically grounded basis for alloy and process design.

Hyer, Holden [ORNL] (ORCID:0000000343915561)↗
DOE OSTI2026

A Python Tool for Aqueous Plutonium Nitrate Density Law Input Preprocessing in MCNP6

Here, this work develops a predictive density tool in Python, named Plutonium Nitrate Solutions (PuNS), to reduce bias and uncertainty in nuclear criticality safety calculations for plutonium nitrate systems. The Pitzer method and an empirical method were implemented into the PuNS tool to generate atom densities for use in MCNP6 material cards. These material cards are directly prepared into an MCNP6 input text file and are calculated based on customizable user inputs of plutonium content, nitric acid content, temperature, and plutonium isotope weight percentages. The PuNS tool is validated and verified against the International Criticality Safety Benchmark Evaluation Project Handbook experiments and is observed to predict densities within a root mean square error of 0.89% for the Pitzer method and 1.82% for the empirical method. These errors in density lead to up to 1569 pcm difference in MCNP6 calculated k eff for the Pitzer method and up to a 1751 pcm difference for the empirical method when compared to experimental benchmarks. Simultaneous work is also being performed at Los Alamos National Laboratory and the University of New Mexico to create a similar tool for plutonium chloride solutions, named Plutonium Chloride Solution, which aims to provide the accreditation of the chlorine absorption. These capabilities will not only provide more accurate models but also facilitate an improved understanding of solution systems and a potential relaxation in the conservatism of current aqueous plutonium processing criticality safety limits.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2022

Conformer-Specific Dissociation Dynamics in Dimethyl Methylphosphonate Radical Cation

The dynamics of the dimethyl methylphosphonate (DMMP) radical cation after production by strong field adiabatic ionization have been investigated. Pump-probe experiments using strong field 1300 nm pulses to adiabatically ionize DMMP and a 800 nm non-ionizing probe induce coherent oscillations of the parent ion yield with a period of about 45 fs. The yields of two fragments, PO 2 C 2 H 7 + and PO 2 CH 4 + , oscillate approximately out of phase with the parent ion, but with a slight phase shift relative to each other. We use electronic structure theory and nonadiabatic surface hopping dynamics to understand the underlying dynamics. The results show that while the cation oscillates on the ground state along the P=O bond stretch coordinate, the probe excites population to higher electronic states that can lead to fragments PO 2 C 2 H 7 + and PO 2 CH 4 + . The computational results combined with the experimental observations indicate that the two conformers of DMMP that are populated under experimental conditions exhibit different dynamics after being excited to the higher electronic states of the cation leading to different dissociation products. These results highlight the potential usefulness of these pump-probe measurements as a tool to study conformer-specific dynamics in molecules of biological interest.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Lagrange-Remap strategy for multi-material fluid-solid simulations using compressive limiters

In the present work, the Lagrange-Remap strategy proposed in [1] is extended to multi-material fluid-solid simulations. Both hypo-elastic and hyper-elastic material models are considered to describe the mechanical behavior of the solids. In practice, the deviatoric stress tensor (for hypo-elastic materials) and the left Cauchy-Green tensor (for isotropic hyper-elastic materials) are remapped, while the use of compressive limiters effectively reduces numerical diffusion during the remapping step. The simplicity of this diffuse interface approach is emphasized in the context of multi-material fluid-solid simulations. A series of Lagrange-Remap test cases, involving both solids and fluids, are conducted and compared with reference Lagrangian simulations, demonstrating the robustness and accuracy of the overall numerical strategy.

Compressive limiters↗
DOE OSTI2026

High-Temperature SiC Cladding End Plug Irradiation Design and HFIR Readiness

This report documents the design of a High Flux Isotope Reactor (HFIR) irradiation experiment intended to evaluate irradiation effects on the hermeticity of silicon carbide (SiC) end plug specimens under a radial fast neutron flux gradient at representative light-water reactor (LWR) temperatures of approximately 300 °C. The overarching goal of this work is to statistically evaluate SiC end plug hermeticity and mechanical properties following irradiation using the high-throughput irradiation capability discussed here. Each specimen consists of a short section of SiC fiber–reinforced SiC (SiC/SiC) tube with a single monolithic SiC end plug joined to one end. The experiment allows up to 66 specimens to be irradiated in six different stacks within a dry, sealed irradiation capsule derived from the previously developed high-temperature SiC/SiC cladding bowing experiment. The neutronics basis, thermal analysis, and HFIR readiness of the experiment are discussed in this report for two possible design cases. The first design case is based on existing approval documentation and components that are on hand and approved for use, so the experiment insertion would require only specimen receipt, specimen pre-irradiation characterization, experiment assembly, and final fabrication package approval. The second design case provides improved thermal robustness and the preferred end plug geometry but requires fabrication of a modified holder and revisions to the HFIR approval documentation, in addition to the other activities required for the first design case, before insertion.

Hott, Daniel [Oak Ridge National Laboratory (ORNL)↗
DOE OSTI2026

First Wall Design of a Tokamak Pilot Plant Using a Monte Carlo Model for 3-D Heat Flux Deposition

We present a method for calculating the heat fluxes deposited on nonaxisymmetric tokamak first wall components, allowing for a first-of-its-kind model for power handling in the tokamak far scrape-off layer (SOL). The DIV3D Monte Carlo model features strict global power conservation and can calculate the finite cross-field plasma transport into magnetically-shadowed regions, which is significant when dealing with meter-scale shadows introduced by components such as poloidal limiters or antennas. As a case study, we apply the DIV3D model to inform the distribution of first wall poloidal limiters in an ARC-class reactor device. We demonstrate that discrete protection limiters can efficiently reduce peak heat fluxes on recessed breeder wall components in the presence of significant far-SOL plasma fluxes. By varying the toroidal periodicity and radial standoff depth of the limiters, we demonstrate one of the tradeoffs that must be considered in first wall design: more limiters provide greater protection, but at the cost of reduced breeding performance. We also present the impact that radial misalignments between limiters would have on first wall power loading.

Monte Carlo methods↗