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 5
All sourcesNASA NTRSDOE OSTI
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 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↗
NASA NTRSDate not supplied

Battery Pack Shape Optimization using Transient Heat Conduction Coupled with Cell-Discharge Analysis

Battery electric systems exhibit significant time-dependence, especially when evaluated in the context of an aircraft mission profile with continually changing power demands. Additionally, when evaluating battery-powered aircraft concepts, it is important to accurately compute the temperature of the batteries and properly characterize the thermal response of the system. The temperature of the batteries has a significant impact on cell performance, in addition to safety considerations of maintaining battery temperatures below their operating limit. Because of these considerations, battery models for preliminary design and optimization of aircraft should include the capability to accurately compute the temperature distribution within the battery pack. Furthermore, battery pack designs should be as light-weight as possible to maximize the pack energy density, while also considering battery temperature limits. Here, we demonstrate a simultaneous trajectory and shape optimization of a battery pack concept, using a transient heat transfer finite element model coupled with a time-varying cell-discharge battery model to provide this capability. The transient finite-element analysis is done using TACS, and the cell-discharge battery model uses OpenMDAO and dymos. Including the transient finite element problem in the loop enables accurate temperatures that can be passed back to the cell discharge model, while the cell discharge model can supply the finite element model with time-varying heat boundary conditions, further benefiting the fidelity of the thermal response of the batteries. We first demonstrate the coupling capability between the battery cell-discharge model and the transient finite-element heat transfer through an optimization which computes the optimal current profile for the battery pack while ensuring the battery temperatures remain below their operational limit. We then build on this optimization by adding shape optimization to the problem, which allows us to consider a composite objective function which also minimizes the mass of the battery pack, while also producing an optimal current discharge profile.

Optimization↗
NASA NTRSDate not supplied

High-Fidelity Modeling Approach for Coupling Thermal, Fluid, and Neutronic Analysis in an NTP Engine

Nuclear thermal propulsion (NTP), with its increased efficiency compared to traditional rocket engines, is a potentially enabling technology for space missions of higher complexity, such as delivering heavy machinery or crew to the moon or other planets. Thus, exploring modeling approaches that can effectively couple the thermal, fluid, and neutronic behaviors of NTP reactors is essential to mitigate cost and schedule challenges with hardware development. This paper details one such approach which utilizes Cardinal to unify reactor physics analysis performed using OpenMC, heat transfer analysis performed using the MOOSE framework, and thermal hydraulic analysis performed using the MOOSE Thermal Hydraulic Module. This approach is demonstrated with a proof of concept that focuses on a legacy reactor geometry from the historic NERVA program. This preliminary model found an overprediction in, when compared to the reference data, axial power peaking and the temperatures of both solids and fluids by roughly 20%. It is clear, however, that the model conserved a number of broad trends, and a few model refinements should reduce the observed error. This approach has definite merit and should enable the production of informative results for more modern NTP systems. This work is a first step in creating a transient reactor model that can be updated across a sequence of time steps based on continually fluctuating, time-dependent boundary conditions.

Nuclear Thermal Propulsion↗
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↗
NASA NTRS1993

Micromechanics of crenulated fibers in carbon/carbon composites

The influence of crenulated noncircular fibers on the micromechanical stress states due to a transverse strain and to a temperature change in carbon/carbon composites is examined using the finite element method. Stresses at the interface of both fully bonded and fully disbonded fibers having two crenulation amplitudes and with two fiber volume fractions are presented. In each case, these interface stresses are compared to stresses at the interface of circular fibers which have the same degree of disbond and fiber volume fraction and are under the same loading conditions. For the disbonded cases, deformed meshes showing locations of fiber/matrix contact are also included. In addition to the interface stress states, selected composite properties are also computed and compared in each case examined. Interest in studying noncircular fibers stems from a desire to increase the transverse properties of carbon/carbon by introducing a mechanical interlocking between the fiber and the matrix. Results presented here indicate that this interlocking does in fact occur. Evidence from the interface stress data suggests, however, that any possible advantage of this interlocking may be outweighed by the disadvantage of stress concentrations which arise at the interface due to the crenulated geometry of the fibers.

E. E. Carapella↗
NASA NTRSDate not supplied

Use of Assistive Technology to Augment API Capabilities

Application Programming Interfaces (APIs) allow for access to data and capabilities of computer applications by developers or users with experience in computer programming. Recent development with both Thermal Desktop and ESATAN-TMS have provided APIs to allow users to develop their own capabilities that interface with the Graphical User Interfaces (GUI) or manipulate the thermal model data. However, these APIs are only as good as the breadth of features in the native code accessible through the API; if a particular code’s feature is not accessible through the API, then users have very limited options besides waiting for updates to the API that expose the necessary functionality, particularly if model data access or user action, such as a button click, is required. However, Assistive Technology features that allow for differently-abled users to more fully experience a software’s capabilities may be creatively utilized to gain further access to data and capabilities not yet exposed by the API. This paper describes the process to augment the features of the OpenTD API via assistive technology and describes how to identify the application instance, navigate GUI elements, updates values on forms, and execute actions such as selecting a listbox item or clicking a button. It concludes with identifying some of the pitfalls to avoid and describes methods to best implement this approach.

Application Programming Interface↗
NASA NTRSDate not supplied

Use of Assistive Technology to Augment API Capabilities

Application Programming Interfaces (APIs) allow for access to data and capabilities of computer applications by developers or users with experience in computer programming. Recent development with both Thermal Desktop and ESATAN-TMS have provided APIs to allow users to develop their own capabilities that interface with the Graphical User Interfaces (GUI) or manipulate the thermal model data. However, these APIs are only as good as the breadth of features in the native code accessible through the API. If a particular code’s feature is not accessible through the API, then users have very limited options besides waiting for updates to the API that expose the necessary functionality, particularly if model data access or user action, such as a button click, is required. However, Assistive Technology features that allow for users with a disability to more fully experience a software’s capabilities may be creatively utilized to gain further access to data and capabilities not yet exposed by the API. This paper describes the process to augment the features of the OpenTD API via assistive technology and describes how to identify the application instance, navigate GUI elements, updates values on forms, and execute actions such as selecting a listbox item or clicking a button. It concludes with identifying some of the pitfalls to avoid and describes methods to best implement this approach.

Application Programming Interface↗
NASA NTRS2020

Evaluation of Fatigue Damage Accumulation Functions for Delamination Initiation and Propagation

The present report follows on the cohesive fatigue damage model methodology proposed in NASA-TP-2018-219838. In that report, an empirical function describing the incremental damage due to cyclic loading was used to calculate fatigue damage within a cohesive formulation. The form of the function was developed such that, when integrated at a constant stress amplitude from no damage to failure, it produces a life versus load response that is consistent with an S-N diagram. Therefore, the parameters of the damage model could be obtained by fitting the model predictions to an S-N diagram. The finite element analyses performed demonstrate that the cohesive fatigue accumulation function provides a link between the S-N diagram that describes crack initiation, and the Paris law that characterizes the rate of crack propagation. However, when the model was proposed, it was not known whether the form of the damage accumulation function associated with a desired S-N diagram is unique and, if not, if the link between S-N and the Paris law is unique and independent of the fatigue function selected. In the effort described herein, several alternative forms of the damage function that reproduce the desired features of S-N diagrams were found and evaluated. The effects of each of these functions on the predicted parameters of the Paris law and the propagation threshold are discussed. The results indicate that the predicted exponent m of the Paris law is indeed independent of the damage accumulation function. However, different functions predict different values for the pre-factor C of the Paris law. Therefore, the proper damage accumulation function must be selected by comparison with experiments. One of the new damage accumulation functions proposed herein was found to be particularly useful for analysis because of the ease with which the model parameters can be determined with a minimal amount of experimental information. The effectiveness of the proposed methodology and damage function was demonstrated by conducting analyses of a double cantilever beam test, a mixed-mode bending test, and a three-point bending test of a skin/doubler specimen. The results indicate that the same set of model parameters can provide accurate predictions of the rate of fatigue crack propagation for a variety of material interfaces, mode mixities, load levels, and stress ratios.

DKIN/Stiffener Debonding↗
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↗
NASA NTRS2026

An Investigation Regarding the Effects of Volume Dynamics on Gas Turbine Engine Performance and Operability

This work explores the impact of adding volume dynamics to an electrified turbofan engine model. The zero-dimensional single-input single-output volume elements capture the physics of fluid compressibility of flow in the engine core gas path turbomachinery during transient operation with greater fidelity. Some background regarding the motivation for this effort is provided. A model of the volume element is developed from first principles. The baseline engine model is described and an explanation is given of the modifications made to incorporate the volume elements into the new model. Details of the simulated inputs, designed to elicit exaggerated dynamic responses and to facilitate comparison between the baseline and modified engine models, are explicated. The results of these simulations are presented and discussed, demonstrating that the transient performance and operability of the modified engine model, measured by parameters such as shaft speeds, thrust, and high-pressure compressor (HPC) stall margin, are impacted by the additional dynamics. Commentary is provided regarding potential explanations for the differences in behavior exhibited by the two models. Finally, suggestions are offered for future work that might improve the volume models developed in this effort and leverage the improved fidelity enabled by the modeling of turbomachinery volume dynamics.

Volume dynamics↗
NASA NTRSDate not supplied

Thermostructural Testing of PICA-D for NASA Planetary Science Missions

Phenolic Impregnated Carbon Ablator–Domestic (PICA‑D) has been selected as the heatshield thermal protection system (TPS) material for two upcoming NASA planetary science missions: the Dragonfly mission to Titan and the Mars Sample Retrieval Lander (SRL). Early testing revealed differences between PICA‑D and heritage PICA, particularly in in‑plane (IP) tensile stiffness and thermal expansion. Thermostructural analyses using Finite Element Method (FEM) tools subsequently predicted the potential for IP compressive failure in the near‑surface layers of PICA‑D under both Dragonfly and SRL flight environments. Over the past three years, the Dragonfly and SRL teams have carried out an extensive thermostructural qualification campaign to address these concerns and validate PICA‑D for flight. This effort began with a comprehensive mechanical property characterization program at Kratos test laboratories, which significantly improved understanding of PICA‑D mechanical behavior and increased the fidelity of FEM predictions. The teams also conducted six large‑scale test entries at the National Solar Thermal Test Facility (NSTTF) solar tower, exposing PICA‑D articles—including gap fillers and representative design features or flaws—to combined thermal and mechanical loads. The talk will summarize key findings from these mechanical and thermostructural test campaigns and present the current status of PICA‑D qualification for NASA’s planetary science missions.

TPS↗
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 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

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↗
NASA NTRS1995

A Probabilistic Model of a Porous Heat Exchanger

This paper presents a probabilistic one-dimensional finite element model for heat transfer processes in porous heat exchangers. The Galerkin approach is used to develop the finite element matrices. Some of the submatrices are asymmetric due to the presence of the flow term. The Neumann expansion is used to write the temperature distribution as a series of random variables, and the expectation operator is applied to obtain the mean and deviation statistics. To demonstrate the feasibility of the formulation, a one-dimensional model of heat transfer phenomenon in superfluid flow through a porous media is considered. Results of this formulation agree well with the Monte-Carlo simulations and the analytical solutions. Although the numerical experiments are confined to parametric random variables, a formulation is presented to account for the random spatial variations.

O P Agrawal↗
DOE OSTI2025

Surface deposition of nanometer scale carbon structures on Bi 2 Sr 2 CaCu 2 O 8+δ using a low-cost scanning electron microscope

In this work, we investigate optically active nanostructures on Bi 2 Sr 2 CaCu 2 O 8+δ (BSCCO). The nanostructures are constructed from carbon nanocrystals formed using the electron beam of a scanning electron microscope to locally decompose residual hydrocarbons at the BSCCO surface. Atomic force microscopy is used to characterize the dimensions of these structures as a function of beam exposure time. This technique has been used to induce localized intercalation to form carbon nanoparticles up to tens of nanometers into a variety of 2D materials. Cross-sectional scanning transmission electron microscopy in BSCCO shows the presence of these carbon deposits on the surface, but there is no evidence they penetrate into the BSCCO substrate. Therefore, this technique is not suitable for localized intercalation in BSCCO, and it does have promise as a cost-effective approach to pattern nanometer scale etch stops.

2D materials↗