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.

SourceResearch connectionFrom yearThrough year

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

417 records · Page 3

All sources NASA NTRS DOE OSTI

DOE OSTI2025

The Role of Asparagine as a Gatekeeper Residue in the Selective Binding of Rare Earth Elements by Lanthanide‐Binding Peptides

Abstract Lanthanide‐binding tag (LBT) peptides selectively complex lanthanide cations (Ln 3+ ) in their binding pockets and are promising for lanthanide separation. However, designing LBTs that selectively target specific Ln 3+ cations remains a challenge due to limited molecular‐level understanding and control of interactions within the lanthanide‐binding pocket. In this study, we reveal that the N5 asparagine residue acts as a gatekeeper in the binding pocket, resulting in a 100‐fold selectivity for smaller Lu 3+ over larger La 3+ cations. Nuclear magnetic resonance spectroscopy and molecular dynamics simulations show that the N5 residue weakly binds to the larger La 3+ cation, permitting H 2 O molecules inside the pocket. For the smaller Lu 3+ cations, the N5 residue forms an inter‐arm hydrogen bond with the E14 glutamic acid residue, locking the Lu 3+ cation in the pocket and preventing H 2 O infiltration. Mutating the N5 asparagine to a D5 aspartic acid prevents such a hydrogen bond, eliminating the gatekeeping mechanism and precipitously reducing selectivity. The resulting binding affinity to Ln 3+ cations is non‐monotonic but generally increases with cation size. These results suggest a molecular design paradigm: the reduced affinity for larger lanthanides is due to open pocket conformations, while the selectivity of smaller Ln 3+ cations over larger ones is due to the gatekeeping hydrogen bond.

Chemistry↗

DOE OSTI2026

SAVY-4000 Finite-Element Drop Test Analysis

PFE Auxiliary Systems conducted drop testing on SAVY-4000 containers to evaluate structural response under 12-foot drop conditions. In support of that effort, a finite-element modeling capability was developed to simulate drop response across multiple container sizes and impact orientations. The purpose of this work was to provide a consistent analysis framework that could support interpretation of testing, compare response trends across multiple configurations, and generate quantities of interest for later comparison with experimental data. More broadly, the analysis and testing were intended to assess whether the containers continued to perform their primary function after a 12-foot drop, namely maintaining structural integrity and containment of the contents. The modeling approach combined an implicit preload analysis with an explicit drop simulation so that each drop event began from a mechanically realistic assembled condition, including compression of the silicone O-ring. Separate models were developed for 2-quart, 5-quart, 12-quart, and 10-gallon containers. The results were evaluated in terms of strain-gauge response, collar-lid gap behavior, and accumulated plastic strain. In addition, parametric studies were performed on the 2-quart container to assess sensitivity to O-ring stiffness, friction, canister thickness, geometry tolerance, and mesh density. The simulations showed that predicted drop responses depended strongly on both container size and drop orientation. Gap metrics identified cases in which the predicted collar-lid opening exceeded the nominal O-ring cross-section threshold, while plastic strain metrics identified localized regions of elevated permanent deformation. Parametric studies showed that the predicted response was especially sensitive to the assumed O-ring stiffness and contact friction, while the geometry tolerance study produced smaller changes in the cases examined. The main value of this work was that it established a repeatable modeling and simulation workflow to support drop-test implementation, evaluate effects of future configuration changes, and understand modeling assumptions that most influenced predicted response. At the current stage, the results were viewed as preliminary model predictions rather than validated predictions. The next step would be to compare drop-test data to the model so that predictive values of the workflow could be refined and used with greater confidence to assess whether the containers maintained structural integrity and containment of the contents after a 12-foot drop.

42 ENGINEERING↗

DOE OSTI2023

Hybrid Data‐Driven Discovery of High‐Performance Silver Selenide‐Based Thermoelectric Composites

Optimizing material compositions often enhances thermoelectric performances. However, the large selection of possible base elements and dopants results in a vast composition design space that is too large to systematically search using solely domain knowledge. To address this challenge, a hybrid data-driven strategy that integrates Bayesian optimization (BO) and Gaussian process regression (GPR) is proposed to optimize the composition of five elements (Ag, Se, S, Cu, and Te) in AgSe-based thermoelectric materials. Data is collected from the literature to provide prior knowledge for the initial GPR model, which is updated by actively collected experimental data during the iteration between BO and experiments. Within seven iterations, the optimized AgSe-based materials prepared using a simple high-throughput ink mixing and blade coating method deliver a high power factor of 2100 µW m −1 K −2 , which is a 75% improvement from the baseline composite (nominal composition of Ag 2 Se 1 ). In conclusion, the success of this study provides opportunities to generalize the demonstrated active machine learning technique to accelerate the development and optimization of a wide range of material systems with reduced experimental trials.

36 MATERIALS SCIENCE↗

DOE OSTI2026

X-ray tomography of damage dynamics in advanced materials using a laser wakefield accelerator

Additively manufactured (AM) metals offer the potential for customizable, cost-effective components, but qualification and certification are crucial. Key to this process is understanding pore dynamics under stress, typically analyzed using micro-computed tomography. This study introduces laboratory-scale “betatron” x-rays from laser wakefield acceleration as a high-throughput alternative for x-ray tomography of advanced materials, such as AM AlSi10Mg alloys. Coupled with 3D finite element modeling, this method provides detailed insights into stress-porosity interactions. The approach delivers high-resolution scans, revealing that pore shape and local triaxiality significantly influence fracture dynamics, supporting advanced material characterization. This work also demonstrates the potential and versatility of laser-betatron x-ray μCT for generating large datasets to accelerate our understanding of the stochastic, process-specific nature of pore formation in AM alloys.

Senthilkumaran, Vigneshvar↗

DOE OSTI2026

Passive Confirmation of the Presence Of High-Explosive Material Via Neutron Transmission Spectroscopy

The goal of this project was to explore a novel approach for identifying presence and potential types of high explosives (HE) in treaty-controlled items with passive neutron sources by using passive neutron transmission spectroscopy. We predominantly look to measure relative elemental abundance of Carbon, Hydrogen, Oxygen, and Nitrogen (CHON) since most relevant materials are certain mix of these elements, as shown in Table 1. Previously, most studies exploring potential identification of CHON elemental content of targets in the vicinity of passive neutron source were focused solely on gamma spectroscopy using high resolution gamma detectors. Using neutron transmission spectroscopy with pulse-shape discrimination (PSD) capable organic scintillators is not only a novel idea in of itself, but arguably necessary for accurate identification of the CHON elemental content of an interrogated target. While high resolution gamma spectroscopy in principle can determine a presence of hydrogen and nitrogen, it is effectively blind to a quantitatively measuring areal densities of C, O. The initially proposed approach, described in detail in the project proposal, was to leverage previous Monte Carlo study on neutron transmission spectroscopy using slab shaped target interrogated with well-collimated (“pencil”) neutron beam. In this study, a simulated test CO 2 target was interrogated by a pencil neutron beam and transmitted neutron spectra were measured using PSD capable organic scintillator using MLEM based unfolding technique. To establish relative elemental content of carbon and oxygen, the unfolded neutron spectrum was fit with a parametrized combination of their respective elemental spectral templates. The individual spectral template was calculated by convoluting ENDF neutron crosssection with the known resolution of the PSD detector used in the study. This approach in the studied configuration successfully quantitatively established relative elemental composition of carbon and oxygen.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2026

The ABCs of phase retrieval: Connecting the acronyms of scanning transmission electron microscopy

High-resolution scanning transmission electron microscopy (S/TEM) is an indispensable tool for characterizing the structure and properties of materials down to the atomic scale. Conventional S/TEM imaging, however, is limited by the phase problem, whereby the phase of the electron exit wave is lost upon detection. Recent advances in diffractive imaging and 4D-STEM have enabled a range of phase-retrieval techniques that computationally reconstruct the missing information encoded in the phase of the transmission function. These approaches offer improved dose efficiency and enhanced sensitivity to weakly scattering signals, extending quantitative imaging to beam-sensitive materials composed of light elements. In this work, we introduce the phase problem in electron microscopy and survey the diverse landscape of phase-retrieval techniques used in the field. Despite their many acronyms and algorithmic variations, these techniques share a common physical and mathematical foundation. We present a unified framework that connects these seemingly distinct methods, from parallax imaging and tilt-corrected bright-field (tcBF-STEM), to aberration-corrected bright-field (acBF-STEM), optimum bright-field (OBF-STEM) and single-sideband (SSB) ptychography, as well as first-moment integrated center of mass techniques (iCOM) and iterative ptychographic algorithms. Based on these insights, we discuss the opportunities and practical limitations of applying these methods across different materials systems, detector designs, and microscope configurations.Graphical abstractRepresentative electron microscopy configurations used for phase retrieval and diffractive imaging in S/TEM: (a) Zernike phase-contrast transmission electron microscopy (TEM), (b) small-convergence-angle four-dimensional scanning transmission electron microscopy (4D-STEM) for nanobeam-based phase reconstruction methods, and (c) large-convergence-angle 4D-STEM for ptychographic and related diffractive imaging techniques reviewed in this work.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Systematic Construction of Time-Dependent Hamiltonians for Microwave-Driven Josephson Circuits

Time-dependent electromagnetic drives are fundamental for controlling complex quantum systems, including superconducting Josephson circuits. In these devices, accurate time-dependent Hamiltonian models are imperative for predicting their dynamics and designing high-fidelity quantum operations. Existing numerical methods, such as black-box quantization (BBQ) and energy-participation ratio (EPR), excel at modeling the static Hamiltonians of Josephson circuits. However, these techniques do not fully capture the behavior of driven circuits stimulated by external microwave drives, nor do they include a generalized approach to account for the inevitable noise and dissipation that enter through microwave ports. Here, we introduce numerical techniques that leverage classical microwave simulations, efficiently executable in finite-element solvers, to obtain the time-dependent Hamiltonian of microwave-driven superconducting circuits with arbitrary geometries under charge, flux, or mixed electromagnetic modulation. Importantly, our techniques do not rely on a lumped-element description of the superconducting circuit, in contrast to previous approaches to tackling this problem. We demonstrate the versatility of our approach by characterizing the driven properties of realistic circuit devices in complex electromagnetic environments, including coherent dynamics due to charge and flux modulation, as well as drive-induced relaxation and dephasing. Our techniques offer a powerful toolbox for optimizing circuit designs and advancing practical applications in superconducting quantum computing.

Lu, Yao [Yale U.; Yale U. (main); Fermilab] (ORCID↗

DOE OSTI2026

Integrated Simulation of Weld Residual Stress Evolution and Crack Propagation Using XFEM

Nuclear power plant components operate in environments that promote multiple degradation mecha- nisms, several of which involve crack initiation and growth. An ongoing effort in the U.S. Department of Energy’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program is developing a general capability within the Multiphysics Object Oriented Simulation Environment (MOOSE) framework for simulating three-dimensional crack growth under a range of driving conditions, including fatigue, stress corrosion cracking (SCC), brittle fracture, and stress-relaxation cracking. This report demonstrates an end-to-end workflow that uses this capability to model weld-residual-stress-driven SCC in the J-groove weld of a pressurized-water reactor control rod drive mechanism penetration in the vessel head. The workflow consists of a thermomechanical welding simulation with temperature-dependent plasticity, followed by cooldown to ambient conditions, and a restart of the simulation using the MOOSE extended finite element method (XFEM) module to propagate a three-dimensional crack through the residual stress field. New welding capabilities were developed to properly initialize newly activated elements in the weld region, and robustness improvements were made to the mesh-based algorithm for defining cutting planes in the 3D XFEM algorithm, allowing it to handle complex crack fronts and stress fields. Together these advances allowed the simulated SCC crack to grow from an initial elliptical flaw in the weld, across the weld, through the tube wall, and almost to the triple point (where the weld, tube, and reactor pressure vessel head intersect) over roughly 36 years of simulated service. These results demonstrate a workflow that can be extended to fully three-dimensional welding simulations and more complex crack interaction problems.

42 - ENGINEERING↗

DOE OSTI2025

Time-Resolved Stochastic Dynamics of Quantum Thermal Machines

Steady-state quantum thermal machines are typically characterized by a continuous flow of heat between different reservoirs. However, at the level of discrete stochastic realizations, heat flow is unraveled as a series of abrupt quantum jumps, each representing an exchange of finite quanta with the environment. Here, in this work, we present a framework that resolves the dynamics of quantum thermal machines into cycles classified as enginelike, coolinglike, or idle. We analyze the statistics of individual cycle types and their durations, enabling us to determine both the fraction of cycles useful for thermodynamic tasks and the average waiting time between cycles of a given type. Central to our analysis is the notion of intermittency, which captures the operational consistency of the machine by assessing the frequency and distribution of idle cycles. Our framework offers a novel approach to characterizing thermal machines, with significant relevance to experiments involving mesoscopic transport through quantum dots.

full counting statistics↗

DOE OSTI2026

Selective depolymerization of polyolefin copolymer via high-entropy alloy catalyzed pyrolysis in a double-layered ultrafast heating reactor

Chemical recycling of polyolefin copolymers (PCPs) is essential for mitigating plastic pollution, yet achieving selective depolymerization into C 2 –C 3 monomers at high productivity remains challenging. Here, we introduce a double-layered ultrafast heating reactor (DUHR) that integrates Joule heating with a high-entropy alloy (HEA) catalyst composed of near-equimolar Co, Cu, Fe, Mn, and Ni. We achieve highly selective depolymerization of PCP into ethylene and propylene, yielding a 70% non-condensable fraction. Significant aromatic fractions form via radical-driven pathways—a behavior distinct from conventional externally heated pyrolyzers—highlighting DUHR's ability to generate reactive species and modulate reaction pathways. Furthermore, these findings demonstrate an effective strategy for chemical recycling of real-world plastic waste for polyolefin circularity and sustainable resource recovery.

Chemical recycling↗

DOE OSTI2026

Minimal dark 𝑆⁢𝑈⁡(2) origin of a massless Dirac neutrino

We propose a gauge-symmetry origin of a rank-two Dirac neutrino mass matrix that enforces one exactly massless neutrino, while being consistent with the oscillation data, as well as cosmological constraints. The mechanism relies on a minimal dark 𝑆⁢𝑈⁢(2) 𝐷 gauge symmetry under which one right-handed neutrinolike Weyl fermion is charged, thereby forbidding its Standard Model Yukawa coupling. Quantum consistency then fixes the minimal dark-sector completion: Cancellation of the Witten anomaly requires a second fermionic 𝑆⁢𝑈⁢(2) 𝐷 doublet, while a discrete 𝑍 4 symmetry that forbids Majorana masses allows the two dark doublets to form a vectorlike pair. This anomaly-free completion gives rise to a secluded, confining dark sector that can contain a potential dark matter candidate, linking the protected neutrino texture to dark infrared dynamics.

confinement↗

DOE OSTI2026

Multitask graph neural networks for elastoplastic response prediction in dual-phase polycrystals

Microstructure-sensitive prediction of elastoplastic response remains a recurring bottleneck in multiscale damage and fatigue modeling, where large ensembles of statistically distinct polycrystals are required to quantify variability and extreme-value behavior. In this work, we develop a multitask graph neural network (GNN) surrogate that maps dual-phase ferrite–martensite polycrystal microstructures to Statistical Volume Element (SVE)-level elastoplastic Quantities of Interest (QoIs). Each SVE is represented as a grain-adjacency graph, with node features encoding phase, geometry, and crystallographic orientation, and edge features encoding relative misorientation. A message-passing graph convolution generates node embeddings, which are pooled into a graph representation and passed to a multitask regression head that jointly predicts 10 scalar QoIs and vector-valued stress–strain responses in orthogonal loading directions across multiple martensite volume fractions and SVE sizes. Results show high accuracy for scalar QoIs and strong agreement for full stress–strain trajectories, with population envelopes reproducing both median behavior and finite-SVE variability across compositions and partition scales. A unified model trained on pooled volume-fraction data preserves most within-regime accuracy relative to regime-specific models while also capturing the broader cross-regime variation reflected in the pooled test set. Distributional comparisons further demonstrate that the surrogate preserves heterogeneity under SVE partitioning, enabling statistically consistent block-wise random-field construction for mesoscale analyses. Overall, the proposed grain-graph surrogate provides a practical pathway to accelerate ensemble-based studies of SVE-level constitutive variability in dual-phase polycrystals.

Crystal plasticity↗

DOE OSTI2024

Nanoscopic Imaging of Self-Propelled Ultrasmall Catalytic Nanomotors

Ultrasmall nanomotors (<100 nm) are highly desirable nanomachines for their size-specific advantages over their larger counterparts in applications spanning nanomedicine, directed assembly, active sensing, and environmental remediation. While there are extensive studies on motors larger than 100 nm, the design and understanding of ultrasmall nanomotors have been scant due to the lack of high-resolution imaging of their propelled motions with orientation and shape details resolved. Here, we report the imaging of the propelled motions of catalytically powered ultrasmall nanomotors─hundreds of them─at the nanometer resolution using liquid-phase transmission electron microscopy. These nanomotors are Pt nanoparticles of asymmetric shapes (“tadpoles” and “boomerangs”), which are colloidally synthesized and observed to be fueled by the catalyzed decomposition of NaBH4 in solution. Statistical analysis of the orientation and position trajectories of fueled and unfueled motors, coupled with finite element simulation, reveals that the shape asymmetry alone is sufficient to induce local chemical concentration gradient and self-diffusiophoresis to act against random Brownian motion. Our work elucidates the colloidal design and fundamental forces involved in the motions of ultrasmall nanomotors, which hold promise as active nanomachines to perform tasks in confined environments such as drug delivery and chemical sensing.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Electrochemical Behavior of Cerium at an Indium Tin-Doped Oxide Electrode in Acidic Media

The redox behavior and speciation of cerium at mesoporous thin films composed of nanoparticles of indium tin-doped oxide (nITO) electrodes were characterized in pH 4.8, 0.1 M acetate buffer and both 0.1 and 1 M HNO 3 using electrochemical techniques and X-ray photoelectron spectroscopy. Anodic deposition of ceria species from Ce(III) to the nITO electrode was achieved under all solvent conditions via spontaneous condensation of electrochemically generated ceric hydroxide species. In 1 M nitric acid, the rate of CeO 2 dissolution is on the same order as CeO 2 deposition, resulting in negligible amounts of CeO 2 electrodeposited at the nITO surface. The cathodic stripping of CeO 2 from the nITO substrate deposited in 0.1 M nitric acid or pH 4.8 acetate buffer follows a 2-step process where Ce(IV)-oxide is initially reduced to an unstable Ce(III)-oxide species that rapidly undergoes acid catalyzed dissolution to yield soluble Ce(III) (aq) . These findings provide a foundation for the pH and anodic potential controlled deposition of CeO 2 thin films to ITO substrates, which can aid in the development of materials composed of ceria. As a result, they can also be used to infer likely analogous actinide redox behavior and speciation at these electrodes.

Cerium↗

DOE OSTI2026

Final Technical Report for DE-SC0022206

This project developed foundational genetic, genomic, and epigenetic tools for anaerobic fungi (Neocallimastigomycota), a group of microorganisms with exceptional natural abilities to deconstruct lignocellulosic biomass. Efficient biomass deconstruction remains a major barrier to economical production of renewable fuels, chemicals, and materials from agricultural and forestry residues. The project sought to enable mechanistic studies and future engineering of anaerobic fungi by improving genomic resources, establishing methods for gene expression, and investigating epigenetic regulation of biomass-degrading pathways. Major accomplishments included generation of the first chromosome-scale genome assemblies for multiple anaerobic fungal species, providing publicly available genomic resources that support both engineering and fundamental biological research. The project established the first reproducible system for heterologous gene expression in anaerobic fungi and identified genomic features and mobile genetic elements that may support future development of stable transformation technologies. In parallel, the project demonstrated direct conversion of untreated lignocellulosic biomass into fuels and specialty chemicals through a fungal-yeast bioprocess and identified anaerobic fungal enzymes with utility for metabolic engineering. The research also revealed that epigenetic regulation plays an important role in controlling fungal gene expression and enzyme production, identifying potential strategies for enhancing biomass degradation. Collectively, this work established anaerobic fungi as a tractable emerging platform for bioenergy and biomanufacturing research, generated valuable public resources, trained the next generation of researchers, and advanced DOE-BER goals related to predictive biology, sustainable bioprocessing, and the circular bioeconomy.

Solomon, Kevin [University of Delaware] (ORCID:000↗

DOE OSTI1961

Laboratory Development of a Pot Calcination Process for Converting Liquid Wastes to Solids

Simulated Darex, Purex, Thorex, and TBP-25 wastes were batch-evaporated and calcined to solids. The weight of residue was decreased to ~ 8% of the original waste for Darex, to ~ 4% for Purex, to ~ 2% for Thorex, and ~ 7% for TBP-25 by calcination to 1200°C. Semicontinuous evaporation and calcination to solids at 900°C of Darex and Purex wastes gave volume reduction factors of ~ 8 and 7 to 10, respectively. The nitrate contents of the residues from batch calcination to 800°C were from ~ 0.07 to 0.5 wt % while the nitrate contents of the residues from semicontinuous calcination to 900°C were from ~ 0.1 to 0.7 wt %. Sodium, calcium, and magnesium additives to Purex waste decreased the percent of sulfate volatilized into the condensate from ~ 30% to <0.5%. Replacement of the atmosphere above the waste with nitric oxide decreased ruthenium volatility from 50 to 60% of that originally present in Purex waste to 0.5 to 3.5%. The thermal conductivities for the calcined wastes, measured in situ, were all >0.1 Btu/hr ft °F at ~ 400°F and increased almost linearly with increasing temperature to >0.3 Btu/hr ft °F at ~ 1600°F in all cases studied.

052001* -- Nuclear Fuels-- Waste Processing↗

DOE OSTI2026

Study of wrap mode impact on Pseudomonas aeruginosa motion in the chemotactic field of fungi

An agent-based discrete computational model biologically calibrated to Pseudomonas aeruginosa migration is used to explore the impacts of bacterial reversals and wrap mode on the efficiency of motion in different environments, both with and without chemotaxis. It is first shown that wrap mode increases the exploration of continuous multimodal chemotactic profiles such as those produced by biologically relevant fungal networks. For cells undergoing a run-reverse pattern, it is shown that the bacteria are likely to remain at the first local chemoattractant maximal production site on a hypha they find. However, with wrap mode, the bacteria can more easily escape these local sites to further explore their neighboring environment along the fungi, suggesting that wrap mode may be beneficial for migration along the fungi in liquid. In a different set of simulations of bacterial motion close to an isolated chemotactic source, wrap mode is shown to increase the ability of a bacterium to reorient toward the source while reducing the overall motion required for similar chemotactic efficiency as a run-reverse strategy, suggesting a potential metabolic benefit. In contrast, model simulations show that wrap mode can increase the rate of dispersal of P. aeruginosa in a nonchemotactic environment.

Hansen, Austin [University of California, Riversid↗

DOE OSTI2025

Large language model-driven database for thermoelectric materials

Thermoelectric materials have the ability to convert waste heat into electricity, offering a valuable solution for energy harvesting. However, their widespread use is hindered by low conversion efficiency, the reliance on expensive rare earth elements, and the environmental and regulatory concerns associated with lead-based materials. A fast and cost-effective way to identify highly efficient thermoelectric materials is through data-driven methods. These approaches rely on robust and comprehensive datasets to train models. Although there are several databases on thermoelectric materials, there is still a need to collect and integrate experimental data from peer-reviewed research articles to capture diverse compositions and properties of materials. Here, in this work, we developed a comprehensive database of 7,123 thermoelectric compounds, containing key information such as chemical composition, structural detail, seebeck coefficient, electrical and thermal conductivity, power factor, and figure of merit (ZT). We used the GPTArticleExtractor workflow, powered by large language models (LLM), to extract and curate data automatically from the scientific literature published in Elsevier journals. This process enabled the creation of a structured database that addresses the challenges of manual data collection. The open access database could stimulate data-driven research and advance thermoelectric material analysis and discovery.

Database↗