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.

600 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

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↗
DOE OSTI2024

Selective solid-state isolation of NMR circuit elements using back-to-back field effect transistors

Nuclear Magnetic Resonance (NMR) electronics that employ selective solid-state isolation of circuit elements can include solid-state switches, such as back-to-back Field Effect Transistor (FET) pairs, and isolated gate drive electronics adapted to operate the solid-state switches in order to selectively decouple induction coils from receive electronics. The solid-state switches can be placed in series to achieve higher standoff voltages, and can be configured for low on resistance and short switching times. The gate drive electronics can include electrical isolation components adapted to enhance standoff voltages and reduce electrical noise at the selectively isolated receive electronics.

Walsh, David O.↗
NASA NTRSDate not supplied

Power and Propulsion Element Steerable High Gain Antenna Lunar Transit Thermal Analysis Tracking Methodology

The Power and Propulsion Element (PPE) is an ion thruster propulsion spacecraft developed as an element of Space Reactor (SR-1) Freedom to provide propulsion, communications and power for the spacecraft. PPE was originally being developed for the use with the lunar orbiting space station Gateway as one of the first two planned elements. PPE was to be launched with the Habitation and Logistics Outpost (HALO) element in a configuration called the Co-Manifested Vehicle (CMV) that would arrive at a Near-Rectilinear Halo Orbit (NRHO) around the Moon via a lunar transit spiral trajectory phase. The PPE communication system is equipped with two Steerable High Gain Antennas (SHGA) each steered by a two-axis gimbal (TAG) mechanism. A thermal analysis was conducted during the near-Earth spiral phase of the mission using Thermal Desktop (TD). This analysis utilized multiple axis Earth tracking articulators in combination with SINDA system internal environmental heating symbols to produce accurate Earth ground station tracking communication system temperatures. This presentation provides an overview of the communication system thermal model and the analysis methodology.

Thermal Analysis↗
DOE OSTI2026

High-Q superconducting lumped-element resonators for low-mass axion searches

Low-frequency superconducting lumped-element resonators have recently attracted significant attention in the context of axion dark matter searches. Here, we present the design and implementation of a fixed-frequency superconducting resonator operating near 250 kHz, possessing an inductor volume of ∼1 liter and achieving an unloaded quality factor Q ≈ 2.1 × 10 6 . This resonator represents a significant improvement over the state of the art and informs the design of searches for low-mass axions.

Bosons↗
DOE OSTI2025

Grain boundary segregation and chemical ordering in CoCrFeMnNi multi-principal element alloy

Owing to their far-from-dilute compositions, multi-principal element alloys (MPEAs) can exhibit unique combinations of engineering properties. As nearly all MPEAs are polycrystalline aggregates, it is necessary to understand the interactions of various elemental species with grain boundaries (GBs). This is of particular importance in extreme environments, such as radiation and elevated temperatures, where such interactions have implications on the properties of MPEAs. Herein, we employ atomistic simulations to generate a series of [001] asymmetric tilt GBs in a model CoCrFeMnNi MPEA and quantify solute interactions and segregation to these boundaries. We employ the Warren-Cowley order parameters to investigate the interplay between GB segregation and chemical short-range order (SRO). At temperatures above 800 K, simulation results reveal the segregation of Cr and Mn to CoCrFeMnNi GBs and show weak dependence of boundary solute excess on GB geometry, at least for the boundaries explored in this work. At temperatures in the range of 673–800 K, formation of domains rich in Cr is observed at GBs in agreement with experimental observations. Quantitative analysis shows that solute excess of various alloying elements decreases rapidly with the increase in temperature in the range of 1000–1200 K. Furthermore, we show that GB regions exhibit SRO characteristics that are distinct from the bulk crystals, leading to spatial variations in SRO. In broad terms, our study highlights the need to account for GB interactions with alloying elements when designing advanced MPEAs with novel chemistries.

Wang, Yitao [Lehigh University, Bethlehem, PA (Uni↗
NASA NTRS1977

Calculations of the Moon's Thermal History at Different Concentrations of Radioactive Elements, Taking Into Account Differentiation on Melting

Calculations of the thermal history of the moon were carried out by solving the thermal conductivity equation for the case in which the heat sources are the long-lived radioactive elements Th, U, and K 40 . The concentrations of these elements were adjusted to give four variations of the heat flow: 1.35 x 10 -8 and 0.91 x 10 -6 cal cm -2 s -1 (I and I**), 0.61 x 10 -6 (variant II, the terrestrial mixture of Lyubimova and Starkova), and 0.236 x 10 -6 (variant III, the chondrite model of Urey and MacDonald). In the same calculations, we considered layering of the differentiated material with transport to the surface of the radioactive elements after the temperature of the layer rose to 200 K above the melting temperature, which is considered in five variants that differ in the amount of transported radioactive elements: 100 percent ( n = 1), 80 percent n = 0.8), 40 percent, 20 percent ( n = 0.6, 0.4, 0.2). During fusion the heat capacity and heat conductivity were changed. We considered two variants of an initially cold (273 K) and hot (900 K) Moon. Calculations show that the interior of the Moon was heated to melting during the first 0.7 to 2.3 x 10 9 years. The maximum fusion involved practically the entire Moon to a distance from 15 to 45 km beneath the surface, and started 3.5 to 4.0 x 10 9 years ago (I, I**), or 2.5 to 3.0 x 10 9 years ago (II, III) and continued for 1 to 2 x 10 9 years. Today the Moon is cooling. The current thickness of the solid crust is from 150 to 200 km and the heat flow exceeds the stationary value 1.5 fold. Apparently the most realistic variant is II (terrestrial mixture) for an initially hot Moon, and gives, regardless of the moderate concentration of radioactive elements, a heat flow of 0.9 to 0.95 x 10 -6 cal cm -2 s -1 , which agrees with the radioastronomical measurements of Troitsky and Krotikov and with the Apollo 15 data.

O I Ornatskaya↗
DOE OSTI2026

Heavy Element Spectroscopy in the Gas Phase

Actinides are inherently unstable and undergo nuclear decay processes with a concurrent release of energy. Consequently, they are used for nuclear power generation, nuclear weapons, and nuclear medicine. However, the radioactive decay processes also pose significant technological problems for the safe treatment and storage of spent nuclear materials. Cost-effective extraction of the actinides is the key first step in the remediation of nuclear waste, but the appropriate chemical means have yet to be determined. Our present understanding of the chemistry of actinides is limited, with the role of the 5f electrons posing a set of particularly challenging questions. The work reported here is focused on the use of electronic spectroscopy to probe the bonding of small molecules in the gas phase that contains thorium or uranium. Analyses of these data, carried out within the framework of ligand field theory, reveal clear evidence that the 5f electrons are spectators that retain their atomic metal ion character.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1977

The Effect of Temperature and Pressure on the Distribution of Iron Group Elements Between Metal and Olivine Phases in the Process of Differentiation of Protoplanetary Material

The distribution patterns of Ni, Co, Mn, and Cr were studied in olivines of various origins: from meteorites (chondrites, achondrites, pallasites), which are likely analogs of the protoplanetary material, to peridotite inclusions in kimberlite pipes, which are analogs of mantle material. According to X-ray microanalysis data, each genetic group of olivines is characterized by a specific concentration of these elements. Nickel is concentrated (up to 0.34 percent) in peridotite olivines, while manganese is concentrated in meteoritic olivines. The maximum chromium content (0.2 percent) was found in ureilites, which were formed under reducing conditions. Experiments at pressures of 20 to 70 kbar and temperatures of 1100 to 2000°C have shown that in a mixture of olivine and Ni metal or NiO nickel enters the silicate phase (up to 4 percent), displacing Fe into the metallic phase. Equilibrium temperatures were estimated from the Fe, Ni distribution coefficients between the metal and olivine: 1500 K for pallasites, 1600 K for olivine-bronzite H6 chondrites, 1200 K for olivine-hypersthene L6, 900 K for LL6, and 1900 K for ureilites (at P = 1 atm). The equilibrium conditions of peridotites are close to T = 1800 K and P over 100 kbar. The distribution patterns of the transition elements are explained on the basis of physical-chemical properties. It is concluded that there is a sharp difference between the conditions of differentiation of the protoplanetary material at the time meteorites were formed and the conditions of differentiation of the planets into concentric layers.

A P Vinogradov↗
DOE OSTI2026

Graphitization of banana peel – synergistic effect of Ca and other elements of the biomass on carbon structure transformation

Recently, growing interest has been observed in using calcium as an effective catalyst for graphitizing disordered carbon materials derived from biomass and bioprecursors, such as cellulose and lignin, at relatively low tem peratures (<2000 °C). Herein, it is demonstrated that in the presence of elements such as calcium, silicon, sulfur, and potassium, which are abundant in banana peel composed of cellulose and lignin base, graphitization occurs much more effectively. The addition of external calcium source, in the form of calcium carbonate, to the banana peel-derived carbon and heat-treatment of the mixture under protective atmosphere up to 1750 °C results in a significant increase in the degree of the graphitization order and local formation of graphite crystals, whereas the same preparation and heat-treatment procedure applied to carbon material derived from pristine cellulose does not lead to such effective graphitization. It is indicated that the graphitization of banana peel carbon occurs due to the synergistic effect of the external Ca-based catalyst and the internal elements present in the raw biomass, supporting the transformation of the disordered graphene-like layers into graphitic structures. These findings are important for the future production of green graphite from biomass

09 BIOMASS FUELS↗
DOE OSTI2026

Compositional Homogeneity in Rapidly Solidified 316L stainless steel produced by laser powder-bed fusion processing

The segregation of alloying elements in rapidly solidified 316L stainless steel produced using additive manufacturing technologies is widely reported throughout literature. The microstructure of thin walls manufactured from 316L stainless steel using laser powder bed fusion are examined for the presence of alloy segregation, i.e. the rejection of solutes to grain boundaries. Metallographic methods including mechanical polishing and chemical etching are used to prepare samples for microanalysis m

Jankowski, Alan F. [Sandia National Laboratories (↗
DOE OSTI2026

Comparing the fluorination of Pb+ and Po+: insights into p-block relativistic effects for superheavy elements

It is expected that the chemical properties of the heaviest of the superheavy elements (SHEs, = 113 – 118) do not align with what is suggested by their current positions on the periodic table. Specifically, the onset of significant relativistic effects, including increased spin–orbit splitting of the p-orbitals may lead to enhanced stability of low-oxidation states. Notably, it is predicted that flerovium ( = 114) may exhibit pseudo-noble gas behavior from a electron configuration that acts as an ‘inert pair’. Even though it is expected that the – splitting becomes pronounced toward the end of the p-block’s sixth row, there is currently limited experimental evidence to confirm its extent or impact. Here, the production of gas-phase lead ( = 82) and polonium ( = 84) fluoride cations (PbF and PoF ) were compared to elucidate differences in accessible oxidation states. The PoF and PbF species were produced and identified with the FIONA spectrometer at the Lawrence Berkeley National Laboratory 88-Inch Cylctron Facility. Polonium showed notably different fluorination chemistry than lead, producing PoF and PoF as primary products compared to lead’s PbF and PbF . The distribution of PoF products observed offer insights as to the role of spin–orbit splitting for polonium. Similar studies of superheavy elements would elucidate the accessibility of their low-oxidation states as well as to inform chemical predictions for eighth-row elements not yet discovered.

Pore, Jennifer L [Nuclear Science Division, Lawren↗
DOE OSTI2026

Influence of Powder Characteristics and Processing Methods on Creep Performance of Powder Metallurgy Hot Isostatic Pressed SS316

The U.S. nuclear energy expansion goals are driving the demand for manufacturing routes that can rapidly produce large, complex, near net shape components. Powder metallurgy hot isostatic pressing (PM HIP) is an advanced manufacturing technique that can be economically scaled-up, while alleviating the supply chain challenges that forging and casting face in terms of cost and lead time constraints. This makes PM-HIP a viable technology to aid and accelerate large-scale part manufacturing for nuclear applications. However, large-scale qualification and deployment of this technology require a thorough understanding of the influence of powder feedstock quality, powder handling history, hot isostatic pressing (HIP) parameters, and subsequent heat treatment on microstructural evolution and elevated temperature mechanical performance. The present work focusses on 316 austenitic stainless steel (SS316) which is most commonly used in high temperature environments for nuclear applications Results from this study show that PM HIPed SS316 meets ASME tensile requirements at room temperature and at elevated temperature. However, creep performance of PM-HIPed SS316 remains inferior to its wrought counterpart, demonstrating that tensile performance alone is not a reliable metric for long duration high temperature integrity. Further, this report delineates powder derived microstructural features that govern creep damage, with key evidences pointing to “microstructural inheritance” from gas atomized powder feedstocks. Commercial SS316 powders of varying chemical compositions and recycling histories were studies, and the results showed large differences in elemental segregation, oxide surface layers and secondary phase distributions. Multi-scale characterization revealed segregation of chromium, molybdenum, manganese and silicon at the boundaries and the precipitation of manganese-, silicon-, and molybdenum-oxides. During HIP consolidation, these surface oxides transform into decorated prior particle boundaries (PPBs) and grain boundary inclusions that persist through conventional post-HIP solution annealing treatment. The retained oxides in post-HIP microstructures were found to influence grain growth, precipitation behavior, and ultimately creep cavitation and fracture. Such post-HIP heat treatments are therefore limited by a complex trade-off between grain growth, and oxide coarsening which aggravate creep damage by acting as nucleation sites for cavities. The objective of this work is to establish an integrated processing–structure–property framework for PM-HIP 316 stainless steel by investigating the influence of powder feedstock characteristics in pre- and post-HIP processing as well as to understand the significance of post-HIP heat treatment on microstructural evolution and creep properties. The results from this report emphasize the significance of powder feedstock integrity in improving creep performance of PM-HIPed SS316, by highlight the effect of rapid solidification, elemental segregation, oxide formation and powder recycling on microstructural defect inheritance following HIP consolidation. Rather than considering HIP processing, solution annealing, and mechanical performance independently, this report treats powder production, HIP consolidation, post-HIP thermal processing, and creep deformation as interconnected stages within a continuous metallurgical process. The resulting framework will provide a scientific basis for developing feedstock engineering strategies capable of improving long-term reliability of PM-HIP stainless steels and accelerating their qualification for advanced nuclear applications.

Ajjarapu, Pavan [Oak Ridge National Laboratory (OR↗
NASA NTRS1977

Geochemical Zoning and Early Differentiation in the Moon

The volatile elements (e.g., Rb, Pb, Tl, Bi, Cs) seem to have been depleted at the time of lunar accretion. Accordingly, it may be assumed that the moon initially accreted from refractory material. The good correlation between volatile/involatile element ratios (e.g., Cs/U, K/La, K/Zr) in both highland and maria samples means that element distribution in lunar crustal rocks is not governed by volatility differences. This and other evidence encourages the view that the moon was accreted homogeneously. A consequence of homogeneous accretion theories is that very efficient large-scale element fractionation is required to account both for the high near-surface concentrations of refractory elements (e.g., Th, U, REE, Zr, Ba, etc.) and for the Ca-Al-rich crust.

S R Taylor↗
NASA NTRSDate not supplied

The Use of Immersion Rigs for High Temperature Hydrogen Exposure Testing within the Nuclear Thermal Rocket Element Environmental Simulator (NTREES): Thermal Soak Rig (TSR)

The Nuclear Thermal Rocket Element Environmental Simulator (NTREES) facility was purpose constructed to perform non-nuclear evaluations of nuclear thermal propulsion (NTP) system fuel materials and structures within prototypic thermochemical environments. This system has been utilized steadily in its ability to subject test specimens to thermochemical and thermohydraulic environments simulating that of an operating nuclear rocket engine. Fission heat is simulated by induction power and experiments are conducted within a ~1000 psi pressure vessel. Hydrogen is conventionally passed through the heated fuel surrogate test specimen while pressure, temperature, and gas species data are collected at various points along the experiment. In order to test fuel and material coupon samples, a class of test apparatus named “immersion rigs” are being developed and employed to more rapidly test these smaller and more technically challenging test specimen. One example of a promising potential fuel structure, Tristructural-isotropic (TRISO) particles, presents unique challenges for testing of this type. TRISO fuel micro-particles are spheroids typically on the order of 500 – 1000 μm in diameter, and exposing a batch sample to hot hydrogen requires purpose-built special test equipment. Thusly, an immersion rig was developed and successfully demonstrated to expose ~1 g of micro-particles to hydrogen gas at temperatures and pressures relevant to NTP systems for the purpose of fuel evaluation. The rig, comprised primarily of graphite and pure tungsten, houses in its core a batch of micro-particles between pucks of porous silicon carbide (SiC). This approach permits gas flow while simultaneously retaining the particles in place. Herein is a discussion of the design, analysis, fabrication, and testing of the NTREES Thermal Soak Rig (TSR).

Space Nuclear Propulsion↗
DOE OSTI2026

Collapse of magnetized white dwarfs as site of heavy-element formation and kilonova signal

We present the first end-to-end calculation connecting the accretion-induced collapse (AIC) of a magnetized, rapidly rotating white dwarf to observable kilonova signatures, combining two-dimensional (2D) general-relativistic neutrino-magnetohydrodynamic simulations, followed by radiation hydrodynamics with in-situ nuclear network and 2D Monte Carlo radiative transfer with spatially resolved heating rates. Unlike all previous unmagnetized AIC models – which predicted proton-rich, $^{56}$Ni-dominated ejecta – strong magnetic fields eject ${\approx }\, 0.2\, \mathrm{ M}_\odot$ of neutron-rich material ($\langle Y_e \rangle \sim 0.24$) on dynamical time-scales, before neutrino irradiation can raise the electron fraction, enabling strong r-process nucleosynthesis up to and beyond the third peak. The resulting kilonova is lanthanide-rich ($X_{\rm lan} \approx 8~{{\ \rm per\ cent}}$) and dominated by near-infrared emission. We compute synthetic light curves in the Large Synoptic Survey Telescope and J ames Webb Space Telescope bands and find striking agreement, without parameter tuning, between the observations of AT 2023vfi/GRB 230307A and our broadband light curves for polar viewing angles. These results establish magnetized AIC as a viable channel for heavy r-process element production and a compelling progenitor candidate for long-duration gamma-ray bursts with kilonova signatures.

MHD↗
DOE OSTI2026

Adsorption of hydroxamic acid ligands for improved extraction of rare earth elements from monazite ores

Efficient separation of rare earth element (REE) ores via froth flotation requires the development of novel ligands with enhanced adsorption capacity and selectivity. To realize these advances, understanding the mechanisms underlying interactions between the ligand and mineral surfaces is essential. This study systematically evaluates the adsorption behavior of alkyl- and aromatic alkyl-substituted hydroxamic acid ligands on monazite surfaces using complementary spectroscopic techniques, including UV–visible (UV–vis) spectroscopy, Raman spectroscopy, infrared spectroscopy, and vibrational sum frequency generation (SFG) spectroscopy, together with the ab initio molecular dynamics (AIMD) simulations. Among the studied ligands, octanohydroxamic acid (OHA) and 4-ethoxy-N,2-dihydroxybenzamide (EDHBA) exhibit high adsorption capacity under basic pH (8–10) by forming multilayers on the surface. OHA has a higher equilibrium adsorption capacity compared to EDHBA, but it forms a less stable multilayer susceptible to disruption in the presence of interfering ions. AIMD results show that OHA adopts a single stable chelating geometry, while EDHBA exhibits multiple binding modes involving distinct interactions with La surface atoms and phosphate-bound oxygens, resulting in more complex adsorption kinetics. The variations in surface binding and intermolecular interactions observed between alkyl and aromatic molecules influence the differences in adsorption kinetics, equilibrium adsorption capacities on the mineral surface, and their flotation performance. This work provides valuable insight into the adsorption mechanism of ligands at mineral interfaces, which is crucial for guiding the design of new ligands with enhanced separation performance.

Zhou, Muchu [ORNL] (ORCID:0000000182650215)↗
DOE OSTI2026

Pinning ångström-size solid ionic channels for rare-earth element separation

High-purity rare-earth elements are essential for modern technologies, yet current solvent extraction processes are energy-intensive and environmentally harmful because of inadequate selectivity and ligand toxicity. Although combining size exclusion and binding affinity can improve lanthanide separation, the role of long-range confinement remains underexplored. Here we report lanthanide separation in aqueous systems using extremely confined manganese oxide solid ionic channels with optimized layer spacing. Different lanthanides induce distinct solid-state phase transformations in manganese oxide, creating a strong driving force for separation. Two lanthanide groups, differing by ~1.4 Å in spacing, were identified and confirmed to be stable by density functional theory. The narrower confinement of heavier Group II lanthanides improves cross-group separation by increasing the dehydration barrier for lighter Group I lanthanides without inducing strong binding. Here, we further developed a strategy to pin the confinement dimensions and enhance same-group separation, increasing enrichment factors for La–Nd and La–Pr pairs from 1.6 ± 0.1 and 1.5 ± 0.1 to 5.4 ± 0.1 and 4.2 ± 0.1, respectively.

Chemical engineering↗