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 4
All sourcesNASA NTRSDOE 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↗
NASA NTRS1977

Lunar Elemental Analysis Obtained From the Apollo Gamma-Ray and X-Ray Remote Sensing Experiment

Gamma-ray and X-ray spectrometers carried in the Service Modules of the Apollo 15 and Apollo 16 spacecraft were employed for compositional mapping of the lunar surface. The measurements involved the observation of the intensity and characteristic energy distribution of gamma rays and X-rays emitted from the lunar surface. A large-scale compositional map of over 10 percent of the lunar surface was obtained from an analysis of the observed spectra. The Apollo 15 flight was at a lunar orbital inclination of 29 degrees as compared with a 9-degree inclination of the Apollo 16 flight; thus, the projected ground track of the Apollo 15 flight covered a larger projected surface area than that of the Apollo 16 flight. The objective of the X-ray experiment was to measure the K spectral lines from Mg, Al, and Si. Spectra were obtained and the data were reduced to Al/Si and Mg/Si intensity ratios and ultimately to chemical ratios. Analyses of the results have indicated (1) that the Al/Si ratios are highest in the lunar highlands and considerably lower in the maria, and (2) that the Mg/Si concentrations generally show the opposite relationship. There is a tendency for the Al/Si values to increase from the western mare areas to the eastern limb highlands. There are distinct chemical contrasts between such features: as the small mare basins and the highland rims. The objective of the gamma-ray experiment was to measure the natural and cosmic-ray-induced activity emission spectrum. At this time, the elemental abundances for Th, U, K, Fe, Ti, Si, and O have been determined over a number of major lunar regions. Regions of relatively high natural radioactivity were found in the Mare Imbrium and Oceanus Procellarum regions. High spots of natural radioactivity were also found south of Fra Mauro; somewhat southwest of Archimedes, and south of Aristarchus. An enhanced region of natural radioactivity was found around Van de Graaff on the far side of tile Moon. In regions other than Mare Imbrium and Oceanus Procellarum, an anticorrelation between natural radioactivity and lunar elevation, as determined from the Apollo laser, altimeter, has been found. From the combined results of the gamma-ray and X-ray spectrometer experiments almost complete information concerning the major element composition of over 10 percent of the lunar surface has been obtained. Distributions have also been mapped for K, Th, and Ti. Interesting correlations between lunar topography and magnetic and gravitational properties have been found.

J I Trombka↗
NASA NTRSDate not supplied

Mission Design Considerations for a Low-Thrust Spacecraft

Developing an executable low-thrust trajectory for use in a spaceflight mission requires the design and optimization of a deterministic trajectory as well as the validation that the selected architecture is robust to some set of uncertainties, execution errors, and potential contingencies. Uncertainty in the ability of the spacecraft and its launch vehicle to execute a trajectory as well as potential deviations such as in-flight anomalies combine with design-to constraints and requirements to complicate the optimization problem. The approach by which robust mission design was accomplished for the initial capability of NASA’s Gateway is presented as well as associated results.

NRHO↗
NASA NTRSDate not supplied

Mission Design Considerations for a Low-Thrust Spacecraft

Developing an executable low-thrust trajectory for use in a spaceflight mission requires the design and optimization of a deterministic trajectory as well as the validation that the selected architecture is robust to some set of uncertainties, execution errors, and potential contingencies. Uncertainty in the ability of the spacecraft and its launch vehicle to execute a trajectory as well as potential deviations such as in-flight anomalies combine with design-to constraints and requirements to complicate the optimization problem. The approach by which robust mission design was accomplished for the initial capability of NASA’s Gateway is presented as well as associated results.

mission design↗
NASA NTRS2026

Global Model Estimates of Atmospheric Al, Ca, Fe, Si, and Ti from Dust and Non-Dust Aerosols Informed by EMIT Surface Mineralogy and Evaluated Against Observations

Atmospheric deposition of micro-nutrients like Fe has been shown to be important for ocean biogeochemistry. The largest source of atmospheric Fe and other elements (e.g., Ca, Al, Si, and Ti) is desert dust, although there are significant non-dust sources in some regions. However, past estimates of these elements have been substantially uncertain due to limited information about the composition of the desert source regions. Here we use elemental distributions estimated from new Earth Surface Mineral Dust Source Investigation (EMIT) observations, which provide mineralogical composition at the surface of the Earth based on imaging spectroscopy measurements from the International Space Station. We add in other sources of these elements (anthropogenic and natural) and compare to a compilation of available surface concentration data from stations over land and from shipborne observations. Our results suggest that the modeled distribution is similar to available observations, but discrepancies still exist in both natural desert dust regions as well as regions dominated by anthropogenic sources. Global budgets for the elements Ca, Al, Fe, Si, and Ti suggest that desert dust remains the dominant source for these elements but anthropogenic or volcanic sources are also important for these elements. Changes in elemental distributions since preindustrial times were also estimated.

aerosols↗
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↗
NASA NTRS2026

Parametric-Based Heat Rejection Trade Study for Lunar and Martian Surface Operations

Establishing and maintaining a sustained presence on the lunar and/or Martian surfaces will require a diverse portfolio of surface elements (e.g., habitation, mobility, power generation, etc.). Many of these systems generate excess heat that must be rejected across a wide range of magnitudes, temperatures, and duty cycles and under variable environmental conditions. To identify the most promising heat rejection approaches for this diverse portfolio, a heat rejection trade study was conducted to evaluate the performance of different technology approaches across a spectrum of surface environments and heat-load requirements. The trade study consisted of three stages: (1) development of a parametric-based modeling framework, (2) creation of a database of heat rejection technologies, surface elements, and environmental conditions for the Moon and Mars, and (3) execution of a quantitative analysis of various heat rejection technologies across different operating conditions and surface elements. The modeling framework is developed in Python and Excel to prioritize small model size and hence low computational cost to enable large parametric sweeps while avoiding the reliance on proprietary software. Individual heat rejection processes are represented as simple Excel models, and a centralized Python script interfaces with the models to coordinate the parametric study. These simple sizing models were developed to take heat load requirements and environmental parameters as inputs and compute mass, power, and volume as outputs. Rather than assess each heat rejection technology separately for each surface element, a unified parametric space was developed to evaluate all technologies across all elements. This parametric space includes factors related to heat load (e.g., magnitude or temperature) and environment (e.g., surface temperature, sky temperature, solar flux). This effort generated a database containing information on over 60 heat rejection technologies and 30 surface elements. For each surface element, the expected heat rejection requirements were documented and analyzed to determine the most common needs shared across all elements. Environmental conditions at various lunar and Martian latitudes were also established for worst-case hot and worst-case cold scenarios. High-fidelity heat rejection models are currently under development. Preliminary trades between heat rejection technologies including radiators, venting technologies, convective coolers, and more have been conducted to identify promising options. This presentation will summarize the preliminary trade results and provide an overview and discussion of the expected heat loads and thermal environments for sustained surface operations on the Moon and Mars.

Heat Rejection↗
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

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

A Missed Thrust Framework for Low-Thrust Spiral Trajectories to the NRHO

A framework is developed by which end-to-end optimization of many-revolution low-thrust spiral trajectories can be completed in the presence of missed thrust events. This framework is applied to the Lunar Transit trajectory by which the initial capability of NASA’s Gateway lunar space station will be delivered to a Near Rectilinear Halo Orbit. This low-thrust mission consists of three subphases, each designed according to the specific objectives and dynamical regimes encountered as the mission progresses from a medium Earth insertion orbit to cislunar space. The presented framework accounts for the unique considerations demanded by each mission phase and incorporates appropriate capabilities into a novel mission analysis tool. This methodology enables large scale and reliable analyses of missed thrust events across the end-to-end Lunar Transit to verify the robustness of flight trajectories across the full range of considered launch dates.

missed thrust↗
NASA NTRSDate not supplied

Softening the Gap between Wöhler and Paris – New Approaches for Fatigue Analysis –

Fatigue analysis tools can vary across industries. For example, automotive engineers often use the Wöhler (S-N) approach to design for safe-life, while aerospace engineers prioritize damage tolerance and inspection intervals, relying instead on crack growth models such as Paris’ law. Although both approaches may deal with the control of cracks in similar materials, their analysis tools and material characterizations are fundamentally distinct. This divide mirrors the classic split between stress-based strength analysis and linear elastic fracture mechanics. However, modern nonlinear models that incorporate material softening, such as cohesive laws, blur this boundary and capture fracture behaviors across scales. This presentation describes the CF23 fatigue model, which uses cohesive softening to link S-N crack initiation with crack propagation rates. CF23 spans the full fatigue spectrum, from initial propagation transients to steady-state growth and threshold conditions, offering a unified framework that bridges Wöhler and Paris-based methodologies. Example applications include fatigue crack propagation transients in adhesive interfaces and skin/stiffener separation.

cohesive elements↗
NASA NTRS2011

Thermal Cycling and Isothermal Deformation Response of Polycrystalline NiTi: Simulations vs. Experiment

A recent microstructure-based FEM model that couples crystal-based plasticity, the B2<-> MB190 phase transformation and anisotropic elasticity at the grain scale is calibrated to recent data for polycrystalline NiTi (49.9 at.% Ni). Inputs include anisotropic elastic properties, texture and differential scanning calorimetry data, as well as a subset of recent isothermal deformation and load-biased thermal cycling data. The model is assessed against additional experimental data. Several experimental trends are captured - in particular, the transformation strain during thermal cycling monotonically increases and reaches a peak with increasing bias stress. This is achieved, in part, by modifying the martensite hardening matrix proposed by Patoor et al. [Patoor E, Eberhardt A, Berveiller M. J Phys IV 1996;6:277]. Some experimental trends are underestimated - in particular, the ratcheting of macrostrain during thermal cycling. This may reflect a model limitation that transformation-plasticity coupling is captured on a coarse (grain) scale but not on a fine (martensitic plate) scale.

Phase Transformations↗
NASA NTRS1992

On the isotopic signature of recent solar-wind nitrogen

One of the most intriguing discoveries yielded by the Apollo samples was evidence pointing towards a significant long-term change in the composition of the sun. Such a change, of the size inferred from the lunar sample data, is inconsistent with present theories of solar evolution. Consequently, there is much interest in exploring this phenomenon as closely as possible, to determine exactly what compositional changes have taken place and whether those changes really did take place in the sun, or whether the cause lies elsewhere. The reason why we can use the moon to analyze the elements in the sun is that the sun emits a stream of ions, known as the solar wind, whose composition, on average, is believed to be the same as that in the surface regions of the sun. When the solar-wind ions hit the surface of the moon, many of them penetrate a short distance into the dust grains lying on the lunar surface. Thus, after a grain has sat on the lunar surface for a while, it has a rim of material that is partly lunar and partly solar in composition. For most chemical elements, the difference between lunar and solar composition is so sufficiently small that the solar elements cannot be detected, but for a handful of elements that are missing from the moon, their solar 'signature' can be observed in samples of lunar soil brought back by the astronauts. Among those elements is nitrogen, the most common element in the air we breathe, but very rare indeed on the moon. Our analytical techniques are not sophisticated enough yet to enable us to analyze individual lunar soil grains for nitrogen, much less to zero in on just the nitrogen in the surface of such a grain. Consequently we are forced to analyze samples consisting of many different grains, each of which could have experienced its own individual history. This makes it difficult to identify the nitrogen implanted in grain surfaces, and also to define the age of a sample.

Kim, Y.↗
NASA NTRSDate not supplied

Construction of an Exact Pressure-Equilibrium Scheme for the Five-Equation Two-Phase Flow Model With Thermal Relaxation

Numerical simulation of compressible multiphase flows based on the four-equation (homogeneous relaxation) model is known to suffer from two fundamental difficulties with (a) wave propagation and (b) pressure equilibrium preservation. First, the mixture sound speed exhibits non-monotonic dependency with respect to the volume fraction, which leads to robustness issues in the resolution of shocks and acoustic wave propagation across two-phase regions. This difficulty can be mitigated by solving Allaire’s five-equation model augmented with infinitely fast phasic temperature equilibrium, from which solutions of the four-equation model can be recovered. However, when temperature is non-uniform, this augmented five-equation formulation still fails to preserve pressure equilibrium across material interfaces. In this work, we propose a fully conservative numerical scheme that exactly preserves pressure equilibrium at the discrete level for the augmented five-equation model, for arbitrary initial distributions of temperature and volume fraction. Combined with the monotonic sound speed property of the five-equation formulation, the proposed pressure-equilibrium preserving scheme significantly improves robustness in the presence of strong multiphase interactions, including shock–interface interactions and advection of material interfaces.

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

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

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↗
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. 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 to the finite element problem, 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↗