Search NASA⌕ Search

SEARCH · Search NASA

Results for “variational method”

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.

618 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS2026

Evaluating Energy Absorption Methods for Integrated Composite Seat Designs

Composite materials have become ubiquitous in the aerospace industry due to their exceptionally light weight and high strength characteristics, as well as their unique ability to be engineered and tailored to meet specific loading conditions and performance requirements. These advanced materials offer superior strength-to-weight ratios compared to traditional metallic materials, making them particularly valuable in weight-critical aerospace applications where every pound saved translates to improved efficiency and performance. In currently operating fleets of commercial and military aircraft, composite materials have been successfully applied to critical structural components, including primary load-bearing elements such as the fuselage sections and flooring structures, which must withstand significant in-flight loads and provide passenger safety. Additionally, these materials have been specifically tailored and optimized for aerodynamic components such as wings and tail assemblies, where their ability to be molded into complex shapes while maintaining structural integrity is particularly advantageous. The application of composite materials extends beyond primary structural elements into the realm of internal cabin components, most notably in innovative seat designs where weight reduction and structural integration are paramount concerns. Modern composite seat structures can be designed to integrate multiple functions, including structural support, comfort features, and safety systems, all while maintaining the lightweight characteristics essential for aircraft performance.

Digital image correlation↗
DOE OSTI2026

Rapid Coal-Ash Characterization using Geophysical Methods & Machine Learning

Coal combustion products (CCP) are challenging to delineate in heterogeneous field settings. Conventional methods (test pits, coring, and laboratory analyses) are labor-intensive, slow, invasive, and provide sparse spatial coverage. This study evaluates whether rapid non-invasive geophysical screening methods—induced polarization (IP), magnetic susceptibility, and nuclear magnetic resonance (NMR) —combined with surface colorimetry (RGB_24), can discriminate CCP-soil mixtures and provide reliable estimates of CCP content. Laboratory measurements were collected on five CCP-soil mixtures (series) and modeled using (i) a linear baseline, (ii) a calibrated non-linear (power-mean) model, and (iii) a machine-learning (ML) Random Forest approach, with validation via leave-one-series-out and site-specific tests. Across the five series, individual signals—particularly IP and magnetic susceptibility—were strongly predictive of ash content but were consistently outperformed by combined models. The pooled calibrated non-linear and ML models captured the observed non-linearity and achieved high accuracy and precision, improving on linear fits. Colorimetry showed the weakest direct relationship with ash content for the tested samples but improved performance when included in multi-signal models. At pre-selected 3.5% decision threshold, calibrated and ML approaches yielded near-perfect classification (Matthews correlation coefficient ˜ 1), suggesting strong practical operability for field screening. Additionally, field-analog tests highlighted the role of endmembers—accuracy declined without access to end-member measurements but was largely recovered by collecting a minimal labeled pair for local recalibration. With end members, accuracy remained high. Globally trained models performed well on three operational unknowns; however, series-specific refits provided the most accurate predictions. Overall, these results highlight the potential of combining rapid geophysics and minimal local calibration for improved coal-ash delineation.

Peshtani, Klaudio↗
DOE OSTI2026

Novel Method for Domestic Stable Isotope Production (CRADA #667) Final Report

This report presents the progress achieved on our chlorine isotope separation efforts in the CRADA #667 agreement. We used a microchannel distillation (MCD) packing structure to enrich HCl isotopes (H35Cl, H37Cl) and an isotachophoresis (ITP) separation to enrich ionic chloride species (35Cl- and 37Cl-). A 2 m tall MCD column was constructed and achieved 264 separation stages over several days, which is well above the goal of 50 stages outlined in the CRADA. The resulting height equivalent to a theoretical plate (HETP) was 0.76 cm. ITP tests were run ~ 24 hours and achieved a single-stage separation factor of 1.559.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2026

SAM Finite Volume Method Development Status Update: GCR Application, Restart, and MultiApp

The System Analysis Module (SAM) is being developed as a modern system analysis code for advanced non-light-water-reactor safety analysis under the U.S. DOE NEAMS program. Previous feasibility studies have demonstrated that a staggered-grid finite volume method (SG-FVM), implemented under the MOOSE framework, can deliver more than an order of magnitude speedup over the existing continuous Galerkin finite element method (CG-FEM) solver for liquid-cooled, incompressible but thermally expandable flow systems. This work extends the previous effort to compressible, gas-cooled reactor applications, where pressure couples directly into the mass equation adding additional nonlinearity into the equation system. New code capabilities are implemented for pebble bed high-temperature gas-cooled reactor (PB-HTGR) analysis, including a pebble bed CoreChannel component, built-in pebble bed effective thermal conductivity model and channel-to-channel crossflow model. The capabilities are tested, benchmarked, and demonstrated for problems with increased level of model and physical complexities, including the HTTU effective thermal conductivity test, the SANA passive cooling test, and a demonstration case using the GPBR200 reactor design covering steady-state operation, DLOFC and PLOFC transients. Across all cases, the SG-FVM solver demonstrated strong robustness and efficiency, and the solutions agree well with reference results and data. The finding of this work proves that SG-FVM is a viable and efficient solver pathway for compressible, gas-cooled reactor system analysis in SAM. In addition, work has been done to successfully support SAM-FVM recover/restart code feature that is essential to reactor safety analysis applications, and MultiApp code feature that is essential to multi-scale and multi-physics simulations. In summary, this work continued from previous feasibility studies, and further demonstrated that the SG-FVM will serve as a strong foundation for SAM’s advanced solver algorithm for future deployment.

Zou, Ling↗
NASA NTRS2026

Thermochemically-Closed Sonic-Flow Inversion for Enthalpy and Temperature in Multispecies Arc-Jet Flows

A thermochemically-closed sonic-flow inversion framework (TSIF) is developed to infer bulk enthalpy and total temperature upstream of a choked nozzle in arc-jet flows. The formulation recasts a pressure-rise total enthalpy quantification technique as an inverse problem in characteristic-velocity c * space using measured mass flow rate, upstream total pressure, gas composition, and nozzle throat geometry as inputs. Unlike calorimetric energy-balance approaches or optical diagnostics, the method relies primarily on routinely measured facility quantities combined with explicit thermochemical closure. Thermochemical states are obtained using NASA’s open-source Chemical Equilibrium with Applications (CEA) code, enabling construction of a chemistry-consistent relation between characteristic velocity, total enthalpy, and total temperature under equilibrium or frozen assumptions. A discharge coefficient is self-calibrated using cold-flow (arc-off) operation data and applied to hot-flow (arc-on) measurements, enabling upstream losses to be accounted for without empirical correlations. The framework is applied to air, N 2 , and CO 2 –N 2 arc-jet flows and demonstrates expected trends for the inferred thermochemical states as function of arc power, specific energy input, mass-flow, heater configuration, and test gas. In the air limit, under equilibrium assumptions, the method recovers the classical high-enthalpy asymptotic correlation of Winovich with a mean residual of 4.4%, demonstrating compatibility with established sonic-flow scaling, while extending applicability to arbitrary multi-species mixtures and non-equilibrium chemistry. The framework provides a mixture-flexible methodology for determining bulk thermochemical states in modern arc-jet environments using routine facility pressure, mass-flow, gas-composition, and nozzle-geometry information together with a cold-flow calibration.

stagnation heat flux↗
NASA NTRS2026

Thermochemically-Closed Sonic-Flow Inversion for Enthalpy and Temperature in Multispecies Arc-Jet Flows

A thermochemically-closed sonic-flow inversion framework (TSIF) is developed to infer bulk enthalpy and total temperature upstream of a choked nozzle in arc-jet flows. The formulation recasts a pressure-rise total enthalpy quantification technique as an inverse problem in characteristic-velocity c * space using measured mass flow rate, upstream total pressure, gas composition, and nozzle throat geometry as inputs. Unlike calorimetric energy-balance approaches or optical diagnostics, the method relies primarily on routinely measured facility quantities combined with explicit thermochemical closure. Thermochemical states are obtained using NASA’s open-source Chemical Equilibrium with Applications (CEA) code, enabling construction of a chemistry-consistent relation between characteristic velocity, total enthalpy, and total temperature under equilibrium or frozen assumptions. A discharge coefficient is self-calibrated using cold-flow (arc-off) operation data and applied to hot-flow (arc-on) measurements, enabling upstream losses to be accounted for without empirical correlations. The framework is applied to air, N 2 , and CO 2 –N 2 arc-jet flows and demonstrates expected trends for the inferred thermochemical states as function of arc power, specific energy input, mass-flow, heater configuration, and test gas. In the air limit, under equilibrium assumptions, the method recovers the classical high-enthalpy asymptotic correlation of Winovich with a mean residual of 4.4%, demonstrating compatibility with established sonic-flow scaling, while extending applicability to arbitrary multi-species mixtures and non-equilibrium chemistry. The framework provides a mixture-flexible methodology for determining bulk thermochemical states in modern arc-jet environments using routine facility pressure, mass-flow, gas-composition, and nozzle-geometry information together with a cold-flow calibration.

inviscid theory↗
DOE OSTI2025

Strong Kitaev Interaction in BaCo 2⁢ (AsO 4 ) 2

The inelastic neutron scattering results and their analysis unequivocally point to a dominant Kitaev interaction in the honeycomb-lattice cobaltate BaCo 2 ⁢(AsO 4 ) 2 . Our anisotropic-exchange model closely describes all available neutron scattering data in the material’s field-polarized phase. Furthermore, the density-matrix renormalization group results for our model are in close accord with the unusual double-zigzag magnetic order and the low in-plane saturation field of BaCo 2 ⁢(AsO 4 ) 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
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 OSTI2025

Thermo-Mechanical Phase-Field Modeling of Fracture in High-Burnup UO2 Fuels Under Transient Conditions

This study presents a novel multiphysics phase-field fracture model to analyze high-burnup uranium dioxide (UO2) fuel behavior under transient reactor conditions. Fracture is treated as a stochastic phase transition, which inherently accounts for the random microstructural effects that lead to variations in the value of fracture strength. Moreover, the model takes into consideration the effects of temperature and burnup on thermal conductivity. Therefore, the model is able to predict crack initiation, propagation, and complex morphologies in response to thermal gradients and stress distributions. Several simulations were conducted to investigate the effects of operational and transient conditions on fracture behavior and the resulting cracking patterns. High-burnup fuels exhibit reduced thermal conductivity, elevating temperature gradients and resulting in extensive radial and circumferential cracks. Transient heating rates and temperatures significantly affect fracture patterns, with higher heating rates generating steeper gradients and more irregular crack trajectories. This approach provides critical insights into fuel integrity during accident scenarios and supports the safety evaluation of extended burnup limits.

Chemistry↗
DOE OSTI2026

Flow synthesis of iron oxide nanoparticles: using multiple precursor additions to improve size control

Controlling the size of iron oxide nanoparticles while maintaining uniformity in flow-based synthesis systems has been a great challenge in nanoparticle synthesis. Using an extended LaMer mechanism, we improve both the size and shape uniformity as compared to a conventional flow synthesis. The key to this approach is the injection of additional Fe precursor during the flow reaction, providing extra precursor during the growth process leading to a size focusing step. This approach has also allowed for the systematic variation of the nanoparticle size produced in the flow reaction. This work represent a new synthetic concept for improving the size uniformity of iron oxide nanoparticles in flow-based synthesis and has potential for application across a wide range of nanoparticle systems.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

Pragmatic Stress Prediction on Additively Manufactured Coupons

Prediction of residual stresses from process parameters for additively manufactured large metal parts is computationally expensive. NASA is currently developing meter-scale parts with direct energy deposition. Practically, the predictive computational methods need to efficiently scale-up to meter-scale parts. Coupled thermal-mechanical multi-physics simulations have been developed with the pragmatic method using ABAQUS, COMSOL Multiphysics, ALE3D software. The residual stresses are a result of the manufacturing process which creates thermal cycling of the build layers. The pragmatic method uses lumped thermal layers for stress predictions to reduce computational costs. The stress predictions as well as deformations of the different codes are compared with each other and with ANSYS Additive using identical material models, boundary and initial conditions. The codes were used to simulate three different geometries: a thin wall, hollow cylinder and twin-cantilever part. The coupon parts were then manufactured with Inconel-625. The residual stresses in these parts were measured using X-ray diffraction as well as neutron beam diffraction at NIST. The stress measurements for the two technologies are compared. The pragmatic stress prediction method enabled predictions of the multi-centimeter scale parts using desktop computer workstations in only a few hours for each coupon. The results of the simulated stress predictions compared favorably with the measured stresses even though thermally lumped layers were employed. Finally, a two-meter scale nozzle was simulated using ANSYS Additive. The simulations were used to examine the build orientation trade-space with respect to resulting geometric deformation. The predicted deformations were compared to measurements of an actual subscale part manufactured with direct energy deposition.

pragmatic method↗
DOE OSTI2026

Elucidating the impact of stress states on grain boundary passivation in Ni-5Cr using the Rhines pack method

The oxidation behavior of a Ni–5Cr (at.%) alloy was evaluated at 420 ° C using the Rhines pack method and simultaneous tensile and compressive stress states via a miniature four-point bending fixture. At this moderate temperature, grain boundaries dominate mass transport and the resulting oxidation response. Oxidation produced approximately 1μ⁢m-wide protective Cr 2 O 3 films capping some grain boundaries and penetrative, intergranular Cr-rich oxides at other grain boundaries. Externally applied tensile and compressive stress during oxidation increased the prevalence of Cr 2 O 3 cap formation compared to no applied stress, with tensile stress resulting in more Cr 2 O 3 caps than compressive stress. The observed oxide cap morphology was similar across all test conditions. Regions under compressive stress showed an order of magnitude greater Cr depletion depth along the grain boundary. Chromium nitrides (CrN), likely from N contamination of the Rhines pack cell, were observed both at the oxide–metal interface and intergranularly within the alloy. Collectively these results demonstrate that applied stress promotes localized protective oxide cap formation over grain boundaries, with compressive stress additionally promoting deeper Cr depletion. Furthermore the experimental approach helped separate out the effects of stress on local oxide formation and grain boundary passivation.

4 point bend↗
DOE OSTI2026

Weak, shallow, dry convection over Angola increases offshore stratocumulus cloud droplet number concentrations

Boundary-layer cloud interactions involving shortwave-absorbing aerosols remain one of the least understood aerosol influences on climate. Here, we find the highest stratocumulus cloud droplet number concentrations over the southeast Atlantic occur when agricultural fires coincide with synoptically-weakened surface warming over Angola, occurring June-early August. Dry convection fills a shallow continental boundary layer with smoke, and a nighttime (local solar time 2-9) land breeze transports the aerosol into the marine boundary layer. Offshore aerosol transport is strengthened by low-level easterlies from a continental pressure high southeast of Angola. Simultaneously, the South Atlantic subtropical high is weaker, allowing extensive dispersal of aerosol offshore into the boundary layer, obscuring cloud brightening from shipping. Meteorological co-variation at synoptic scales compensates for cloud brightening by the smoke. Outgoing shortwave radiation increases by 15–20% of the monthly mean in June and July when offshore droplet numbers are less but the stratocumulus deck is more developed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Dopant adsorption as a function of bulk concentration at a near 42º (100) twist grain boundary in SrTiO3

The enrichment of grain boundaries with dopant atoms is of critical importance for the macroscopic physical properties of materials. In thermodynamic equilibrium the Gibbs adsorption isotherm relates grain boundary excess of dopant atoms, their chemical potential in the adjacent bulk, and the respective interface energy. This study has used bicrystals with a near 42º (100) twist grain boundary in SrTiO3 to demonstrate that different kinetic pathways of Fe dopant atom additions converge towards comparable grain boundary configurations. Despite differences in bulk chemical potentials remarkably similar grain boundary excess quantities were observed. The experimental results indicate the general experimental feasibility to establish quantitative relationships between variations of grain boundary energy, interfacial excess, and overall dopant concentration. Improved experimental counting statistics are needed to distinguish solute interface excess as a function of bulk chemical potential.

Hahn, William [University of California, Davis]↗
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↗
NASA NTRS1993

Formal Functional Test Designs with a Test Representation Language

The application of the category-partition method to the test design phase of hardware, software, or system test development is discussed. The method provides a formal framework for reducing the total number of possible test cases to a minimum logical subset for effective testing. An automatic tool and a formal language were developed to implement the method and produce the specification of test cases.

J M Hops↗