Search NASA⌕ Search

SEARCH · Search NASA

Results for “Discrete Element 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.

700 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Flammability Test Method for Materials Intended for Use as Fire Barriers in Crewed Habitats

In some crewed-spacecraft applications, onboard materials are intentionally selected to function as hazard control barriers to prevent powered hardware and/or ignited flammable items from igniting nearby flammable materials by acting as non-breaching, non-flammable containment (e.g., cargo transfer bags, zero gravity stowage racks, and jettison stowage bags on the International Space Station). Historically, flammability self-extinguishment testing of broad acreage or outermost layer materials alone was believed to be sufficient to verify overall flammability compliance. Though required, testing on design features such as threaded seams and zippers was not performed due to the incorrect assumption that flammability performance would not be impacted. However, recent coupon level testing of layups and configurations consisting of multiple materials have revealed that these features may be more susceptible to ignition and fire propagation than the bulk material. Additionally, configurational coupon ignition and propagation testing does not explicitly determine a material’s ability to provide protection to vulnerable flammable hardware. The need to evaluate this gap in assessing material or layup effectiveness as fire barriers between ignition sources and flammable materials led to the development of a new configuration-based fire barrier test. Test setup, relevant configurations, and test validation are discussed. The new barrier test method establishes a configuration for the test apparatus, identifies a conservative ignitor configuration, and defines pass/fail criteria. Test setup, relevant configurations, and test validation are discussed.

Susana Harper↗
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

Innovations Driven by Advanced Characterization to Strategize Critical Mineral Production and Beneficial Reuse from Fossil Energy Waste

Critical minerals (CM), such as rare earth elements (REE), cobalt, nickel, and lithium, have important uses in modern electronics and advanced manufacturing, yet are vulnerable to potential supply chain disruptions. Relatively abundant and readily available fossil energy (FE) wastes, such as coal combustion ash, acid mine drainage (AMD) and treatment solids (AMD solids), and Oil and Gas (O&G) drilling wastes (drill cuttings and produced waters) are under consideration as CM feedstocks. The National Energy Technology Laboratory (NETL) has studied CM resources for various FE wastes as part of the U.S. Department of Energy’s mission of bolstering the domestic CM supply, and makes the data available to the public on EDX at sites such as the NEWTS group. Advanced characterization utilizing synchrotron x-ray techniques coupled with laboratory extractions has been performed to identify CM hosting phases in these FE wastes to inform CM recoverability mechanisms. Novel methods to selectively recover CMs while co-producing other valuable byproducts have been developed. Successful examples discussed here include: (1) The identification of REE/Co/Ni/Sc binding and hosting phases in select FE waste (coal combustion ash and AMD solids), resulting in the development of a patented CM step-extraction process, (2) coupled production of functional sorbents from these extraction wastes and for CM recovery. A pilot-scale testing to evaluate the patent’s technical feasibility for extracting REE from coal ash on a barrel scale has been successfully performed. Additionally, (3) evaluation and measurements of brine geochemistry from U.S. O&G produced waters has informed a high Li recovery potential from Marcellus Shale produced water. NETL researchers have been developing tailored pre-treatment processes, an innovative and highly durable lithium sorbent, and geochemical model guided precipitation to accelerate Li production from the Marcellus Shale produced waters. These innovations driven by characterization are integral for maximizing and advancing the potential for CM recovery while offsetting the cost and environmental footprint for FE waste management.

critical mineral processing↗
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

Compatibility of Type 304H Stainless Steel in Static Molten FLiBe for Inertial Confinement Fusion Reactors: Role of Impurities and Redox Control

Molten fluoride salts, such as FLiBe (LiF-BeF2), are promising candidates for tritium breeding and heat transfer in fusion reactors, but corrosion of structural materials remains a major challenge. This study investigates the corrosion behavior of austenitic stainless steel 304H in purified and NiF2-containing FLiBe at 500°C and 600°C, focusing on the effects of impurities and redox control. Exposure to purified FLiBe resulted in the concurrent depletion of Cr, Mn, and Fe, with corrosion at 500 °C dominated by the combined oxide formation and elemental dissolution, while at 600°C elemental depletion was predominant. The addition of a controlled NiF2 impurity significantly accelerated corrosion at both temperatures, demonstrating the sensitivity of 304H to the salt redox state. Beryllium additions were effective in mitigating corrosion for both baseline and NiF2-containing FLiBe; minimal depletion of Cr, Mn, and Fe occurred with Be additions as low as 2.5 mg (147 wppm), and no NiBe intermetallics formed at 5 mg (294 wppm), indicating that small Be inventories can provide substantial protection without deleterious phase formation. Thermodynamic equilibrium and coupled thermodynamic-kinetic analyses at the salt-alloy interface suggested low corrosion rates in systems with limited hydrogen fluoride (HF) generation, highlighting the importance of salt redox control. Estimates for a Be addition rate were calculated for the HYLIFE-II fusion reactor that can mitigate corrosion-induced degradation, assuming complete conversion of tritium to tritium fluoride (TF). Overall, 304H shows reasonable compatibility with FLiBe under optimized redox conditions. These results provide quantitative guidance for material selection and salt management in fusion blanket and heat exchanger systems and motivate validation under flowing, nonisothermal, and irradiated conditions.

Pillai, Rishi [ORNL] (ORCID:0000000243688197)↗
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↗
NASA NTRS1990

Solar Array Automation Limitations

Significant progress in the automation of the spacecraft electrical power systems has been made within the past few years. This is especially important with the development of the space station and the increasing demand on the electrical power systems for future satellites. The key element of the spacecraft power system, the solar arrays which supply the power, will have to grow to supply many tens of kilowatts of power within the next twenty years. This growth will be accompanied by the problems associated with large distributed power systems. The growth of the arrays, the on-array management problems and potential solutions to array degradation or failure are discussed. Multilowatt arrays for unmanned spacecraft with comments on the implications of array degradation for manned spacecraft are discussed.

Terry M Trumble↗
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↗
NASA NTRS2026

Expanding Space with Inflatable Softgoods: Roadmap for In-Space Manufacturing of Resilient Space Structures

Overview of inflatable softgoods and their emerging role in enabling large, resilient space habitats and infrastructure. It highlights the advantages of inflatable systems, such as exceptional packing efficiency and scalable habitable volume, while also addressing challenges related to outfitting, complex material behavior, structural design, manufacturing precision, and testing limitations. Core architectural elements of crewed inflatable habitats are described, along with shell layer composition, structural interfaces, and examples of conceptual habitat configurations for transit, lunar, and surface applications. The presentation concludes by outlining key technology shortfalls, including structural health monitoring, ultra‑high‑strength materials, lifetime performance, and integration strategies, emphasizing the need for continued development to support future in‑space manufacturing and exploration missions.

Habitat↗
NASA NTRS1977

Radioactivity of the Moon, Planets, and Meteorites

This paper reviews and summarizes analytical data for the content of natural radioactive elements in meteorites, eruptive terrestrial rocks, and also, in lunar samples returned by Apollo missions and the Luna series of automatic stations. The K-U systematics of samples analyzed in the laboratory are combined with data for orbital gamma-ray measurements for Mars (Mars 5) and with the results of direct gamma-ray measurements of the surface of Venus by the Venera 8 lander. Using information about the radioactivity of solar system bodies and evaluations of the content of K, U, and Th in the terrestrial planets, we examine certain aspects of the evolution of material in the protoplanetary gas-dust cloud and then in the planets of the solar system.

Yu A Surkov↗
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 OSTI2025

Tree tensor network hierarchical equations of motion based on time-dependent variational principle for efficient open quantum dynamics in structured thermal environments

In this work, we introduce an efficient method, TTN-HEOM, for exactly calculating the open quantum dynamics for driven quantum systems interacting with highly structured bosonic baths by combining the tree tensor network (TTN) decomposition scheme with the bexcitonic generalization of the numerically exact hierarchical equations of motion (HEOM). The method yields a series of quantum master equations for all core tensors in the TTN that efficiently and accurately capture the open quantum dynamics for non-Markovian environments to all orders in the system–bath interaction. These master equations are constructed based on the time-dependent Dirac–Frenkel variational principle, which isolates the optimal dynamics for the core tensors given the TTN ansatz. The dynamics converges to the HEOM when increasing the rank of the core tensors, a limit in which the TTN ansatz becomes exact. We introduce TENSO, tensor equations for non-Markovian structured open systems, as a general-purpose Python code to propagate the TTN-HEOM dynamics. We implement three general propagators for the coupled master equations: two fixed-rank methods that require a constant memory footprint during the dynamics and one adaptive-rank method with a variable memory footprint controlled by the target level of computational error. We exemplify the utility of these methods by simulating a two-level system coupled to a structured bath containing one Drude–Lorentz component and eight Brownian oscillators, which is beyond what can presently be computed using the standard HEOM. Our results show that the TTN-HEOM is capable of simulating both dephasing and relaxation dynamics of driven quantum systems interacting with structured baths, even those of chemical complexity, with an affordable computational cost.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗