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.

445 records

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

Simulations of Yarn Micro-Mechanics of Woven Heat Shield Materials

Carbon and phenolic fibers are commonly used in ablative thermal protection materials, such as 3-dimensional Mid-Density Carbon Phenolic (3MDCP), a 3D-woven composite comprised of mixed-fiber yarn bundles. Predicting the micro-mechanical response and fracture of twisted yarns composed of brittle and ductile fibers requires a modeling approach that captures per-fiber yielding, fiber fracture, and inter-fiber friction and contact. This work presents an extended bonded particle model (BPM) for discrete element method (DEM) simulation of fiber and yarn mechanics, implemented in LAMMPS. The model builds on the incremental bond formulation of Guo et al. and introduces a piecewise elasto-plastic constitutive law for axial extension, enabling representation of fibers that yield before failure. 3MDCP yarns were constructed using measured fiber radius distributions and helical twist geometry. Tensile simulations of single-ply 3MDCP yarns show good agreement with vender stress–strain results. Fiber breakage models also show details on yarn breakage propagration, centered radially in the yarn. Yarn breakage of multi-ply 3MDCP also matched experimental observations in per-ply breakage; however, predicted yarn breakage strength were found higher than experimental observations.

Discrete Element Method↗

NASA NTRSDate not supplied

Simulations of Yarn Micro-Mechanics of Woven Heat Shield Materials (3MDCP)

Woven↗

NASA NTRSDate not supplied

Thermomechanical Modeling of Woven Materials With Particle-Based, Explicit-Fiber Simulations

Fiber-based materials are extensively used to protect spacecraft during entry. Insulative fibers, often in a fiber network or woven, provide rigidity, strength, and control of material anisotropy and density. Woven thermal protection materials, such as ADEPT (Adaptable, Deployable Entry and Placement Technology), 3D-MAT (3-Dimensional Multifunctional Ablative Thermal Protection), and 3MDCP (3D Woven Mid-Density Carbon Phenolic), enable missions with stronger and denser materials for entry profiles with high shear and heat flux. Vulnerabilities to woven thermal protection materials include manufacturing-induced material property variation, and impact from micrometeoroids. Simulating woven materials under these conditions require models that can resolve hierarchal structures, thermomechanical behavior, and failure. To address this, we simulate weave thermal conduction and mechanical deformation. We simulate the full weave with a coarse-grained yarn model is presented. The model combines a validated, high-resolution single 3MDCP yarn model and phenolic resin model. Instead of modeling every fiber, each yarn ply with order 10, instead of order 1000, fibers. The discrete element bonded particle model (DEM-BPM) of fibers captures the thermal and mechanical behavior within and between fibers. We study the proportion of heat transfer and stress via the contact network, fiber bonds, and overall weave geometry.

bonded particle↗

NASA NTRSDate not supplied

Discrete Rotational Energy for Polyatomic Molecules in Direct Simulation Monte Carlo

Accurate prediction of aerothermodynamic loads in thermal non-equilibrium flows requires precise modeling of internal energy exchange. While previous direct simulation Monte Carlo frameworks have successfully implemented discrete rotational energy models for diatomic species, the treatment of polyatomic molecules has traditionally relied on continuous energy assumptions that break down at low temperatures and neglect critical high-temperature corrections. This study extends the discrete rotational energy models of Boyd and Gimelshein to fully encompass polyatomic molecules. The proposed framework implements quantized rotational energy level sampling for linear, spherical, and symmetric/asymmetric top rotors. Crucially, the model incorporates centrifugal distortion to address the limitations of the rigid-rotor assumption at hypersonic temperatures, and accounts for nuclear spin parity, which dictates the permissible rotational states and macroscopic specific heats at low temperatures. The model is verified through equilibrium sampling procedures, demonstrating agreement with theoretical quantum Boltzmann distributions and accurately reproducing thermophysical properties across a wide range of temperatures.

DSMC↗

NASA NTRSDate not supplied

Discrete Rotational Energy for Polyatomic Molecules in Direct Simulation Monte Carlo

The environments experienced by vehicles upon entry into planetary atmospheres generally result in severe aerothermodynamic loading. These flows are characterized by the formation of strong shock waves, behind which high temperatures and non-equilibrium environments are generated, and the accurate prediction of internal energy relaxation and chemical kinetics becomes important. This manuscript focuses on the internal relaxation of the rotational energy of polyatomic molecules. Previous work has described the implementation of a discrete rotational energy model for diatomic molecules (linear rotators). Later, Gimelshein et. al extended this model to include polyatomic molecular internal energies. However, they did not include details on how to generate the particles at equilibrium conditions and centrifugal distortion was not included in their model. The purpose of the present study is to extend the afore-mentioned models to the treatment of discrete rotational energy in polyatomic molecules with the addition of centrifugal distortion. An overview of the background theory needed for the implementation of the model will be given followed by the model itself and verification of the model. For example, sampled versus Boltzmann distributions are compared in Figure 1 for CH3 (oblate rotator – two quantum numbers) and examples of thermal relaxation behavior for linear and non-linear rotators are presented in Figure 2. In addition, comparisons of thermophysical properties will be presented, as well as examples of adiabatic relaxation and application to hypersonic flow.

Rarefied Gas Dynamics↗

DOE OSTI2026

Optimal binning of correlated measurements

Experimental measurements are commonly represented on a discrete grid, requiring a balance between granularity and statistical noise. Two strategies have traditionally been used to improve such representations: selecting an appropriate bin width to control discretization error and applying kernel-based smoothing to suppress fluctuations. Despite their shared goal, these approaches have largely developed independently, without a unified statistical description of how discretization and correlation jointly determine measurement precision. Here, we extend the discussion of optimal interval averaging to a correlation-aware setting by Gaussian process regression, which explicitly accounts for correlations among neighboring bins. Starting from first principles, we derive the mean-squared error of discretized measurements and obtain closed-form asymptotic expressions for the optimal bin width and correlation length. When recast in reduced variables, the theory reveals distinct universal scaling laws governing the error in the correlation-free and correlation-controlled regimes. Characterized by intrinsically smooth intensity profiles and counting-based statistics, neutron scattering measurements are well suited for demonstrating the enhanced error contraction enabled by inter-bin correlations. We show that such improvement is achievable over the experimentally accessible Q-range and across multiple instruments and material systems. These results show that explicitly accounting for correlations systematically reshapes the limits of precision in discretized, noise-limited measurements. More broadly, the framework provides a transferable statistical foundation for optimizing data representation, inference, and experimental design across the physical and data sciences.

Tung, Chi-Huan [ORNL] (ORCID:0000000221972074)↗

NASA NTRSDate not supplied

An Efficient Approach to Collect Inlet Characterization Data

This report considers a novel process to acquire performance data for supersonic mixed compression inlets that is faster, and therefore less costly, than the previous one. In the previous process, measurements were recorded in discrete increments. The backpressuring actuator was advanced, there was a pause to allow for airflow fluctuations to diminish, and then the measurements were recorded. While this process was methodical, reliable and accurate, it consumed significant wind-on time. In general, high-speed wind tunnels are expensive to operate considering maintenance, energy and human resource requirements. Furthermore, many facilities, such as blowdown wind tunnels, are limited in the duration that they can maintain the desired test condition. Consequently, a quick method to collect the measurements is desired. A favored approach has the actuator moving in a continuous fashion with measurements recorded as the actuator continuously progresses through the increment points. This approach, identified as the dynamic inlet characteristic data acquisition procedure, uses in situ dynamic high-speed pressure sensors to measure pressure signals. While faster, continuous movement adds dynamic effects to the data. To better understand these dynamic effects and other potential influences, a study was undertaken to compare wind tunnel measurements taken at discrete actuator position points with data gathered during continuous actuator movement. Specially, plots of inlet pressure recovery versus mass capture ratio, or total pressure characteristic curves, were examined. The data were measured during experiments in the NASA Glenn Research Center Abe Silverstein 10- by 10-Foot Supersonic Wind Tunnel (SWT) with an inlet propulsion test article. In this report, the wind tunnel experiment and the study objectives are described. The processes to reduce the data and for analysis are explained. A discussion of the results along with features noted in the data is given. Finally, conclusions from this study that may be used to guide future wind tunnel experiment planning are presented.

Data Acquisition↗

NASA NTRS2026

A Modular Conjugate Heat Transfer Optimization Framework for Thermal Management of Electric Aircraft

Conjugate heat transfer (CHT) analysis and optimization is a powerful method for improving thermal management, as it simultaneously resolves the temperature distribution in both fluid and solid domains. This paper presents a modular, discrete adjoint-based CHT optimization capability integrated within the OpenMDAO/MPhys framework. A unique feature of the proposed framework is its flexibility to extend to multidisciplinary optimization, including aero-structural-thermal applications. The fluid domain is modeled using a finite-volume Computational Fluid Dynamics (CFD) solver, and the solid domain with a conduction heat transfer solver. A mixed Neumann-Dirichlet boundary condition is developed to enable full submersion of the solid geometry within the fluid domain, while ensuring consistent temperature and heat flux coupling at the CHT interface. Gradient-based optimization is performed; the gradients are efficiently computed using the discrete adjoint solvers implemented in DAFoam. To demonstrate the method, this paper considers two cases related to electric aircraft thermal management: a U-bend heat exchanger and an actively cooled battery pack. The U-bend case aims to minimize pressure loss while maximizing heat flux by changing the pipe geometry. The optimized design reduces pressure loss by 52.7% and increases total heat flux by 2.3%. In the battery pack case, a 3-by-3 cell configuration is cooled by ambient airflow, with constant heat generation prescribed in the cells. The battery casing shape serves as the design variable, and the objective function is a weighted sum of pressure loss and pack weight, subject to a maximum temperature constraint. The optimized design achieves a 44.6% reduction in pressure loss and a 1.5% reduction in weight, while satisfying the thermal constraint. To ensure the reliability of the optimized designs, this study validates coarse-mesh, steady-state predictions against fine-mesh unsteady simulations, demonstrating consistency within acceptable errors. This work demonstrates the potential of the developed framework to enable rapid, high-fidelity design of thermal management systems for electric aircraft.

heat transfer↗

NASA NTRSDate not supplied

Multidisciplinary Design Optimization and Analysis of an Open Rotor Stage: Part 1

Successful design of open rotor propulsors requires effective analysis across multiple disciplines, including aerodynamics, acoustics, and structures. A viable design must not only be efficient but must also produce an acceptable level of noise and meet all static and dynamic structural requirements. For design and optimization, this is especially challenging because running high fidelity analyses is resource-intensive, and optimizing a design may require many thousands of cases to be analyzed. For this reason, the NASA team has applied design methodology that utilizes low-cost aerodynamic methods, machine learning models, and high-fidelity analyses when necessary. This includes standard two-dimensional methods such as throughflow analysis and quasi-3D blade-to-blade CFD, as well as some newly developed methods. Optimization using 3D CFD is necessary to maximize performance, and this is considered as well. All optimizations are carried out subject to structural constraints evaluated using finite element analysis. Doing this accurately requires a robust trunnion design, capable of pitching the blade between cruise and takeoff conditions while maintaining acceptable factor of safety. Hot to cold analysis must also be applied in order to correctly determine the as-manufactured shape. For acoustics, the unsteady pressures on the blade surfaces must be predicted, and this can be done either through full-annulus unsteady CFD or through a nonlinear harmonic method (NLH). NLH can run much faster, allowing some acoustic considerations to be made earlier in the design process. The design process is ongoing, and this presentation will review the current status and planned next steps. This part of the talk will focus on aerodynamic performance and be followed by a talk on structures and acoustics.

Design↗

DOE OSTI2026

Bacterial Bioleaching and Biorecovery for Biomining Unconventional Rare Earth Element Feedstocks

Bacterial metabolic interactions with rare earth elements (REEs) can be harnessed for biomining unconventional feedstocks like abandoned coal-mine drainage (AMD). Pennsylvania has ~500 AMD passive remediation systems that can precipitate REE rich solids. REEs include yttrium and the lanthanide series that are used in modern energy and technology. Bacteria that metabolically interact with REEs can be used for biomining in an affordable efficient process that does not require hazardous chemical additives. Currently, the microbial metal mechanisms that contribute to REE biorelease and biorecovery are poorly understood. Our work shows acidogenic bacterial isolates (Bacillus mycoides JR07 and Bacillus pseudomycoides KB7) successfully bioleach a mixed REE solution from AMD solids by their organic acid production and biofilm formation. Further, our work shows the potential for bacterial lanthanide-dependent enzymes to recover lanthanides from a mixed REE solution; here we have bacterial isolate Methylobacterium sp. B3 that can recover soluble lanthanum. Whole genome sequencing of Methylobacterium sp. B3 predict lanthanide-dependent methanol dehydrogenase XoxF. Understanding the microbial metabolism and genes involved in the REE release and recovery is crucial to optimize the biomining of AMD solids. Our work addresses the growing need to develop novel REE mining methods from unconventional feedstocks.

biogeochemistry↗

NASA NTRS1962

A Design Solution to the Spacecraft Look-Angle Problem

As the state of the art of the design of spacecraft has progressed, one of the difficult problems which has arisen is the positioning and articulating of spacecraft elements which must “see” in different directions without interference. Some of these elements are: solar panels or collectors, Sun sensors, communication antennas, star seekers, planet tracking and scanning devices, rocket motors, attitude-control jets, and scientific instruments. This Report presents a method of testing a spacecraft design to determine how well it satisfies these look-angle requirements. This method is applied to both simple and complex cases, and examples of its use are presented. A means of assessing spacecraft constraints on trajectories is discussed. The method was developed for the Mariner interplanetary spacecraft, but should be applicable to other cases.

MARINER PROGRAM↗

NASA NTRS2026

Evaluating Energy Absorption Methods for Integrated Composite Seat Designs

Composite materials have become ubiquitous in the aerospace industry due to their exceptionally light weight and high strength characteristics, as well as their unique ability to be engineered and tailored to meet specific loading conditions and performance requirements. These advanced materials offer superior strength-to-weight ratios compared to traditional metallic materials, making them particularly valuable in weight-critical aerospace applications where every pound saved translates to improved efficiency and performance. In currently operating fleets of commercial and military aircraft, composite materials have been successfully applied to critical structural components, including primary load-bearing elements such as the fuselage sections and flooring structures, which must withstand significant in-flight loads and provide passenger safety. Additionally, these materials have been specifically tailored and optimized for aerodynamic components such as wings and tail assemblies, where their ability to be molded into complex shapes while maintaining structural integrity is particularly advantageous. The application of composite materials extends beyond primary structural elements into the realm of internal cabin components, most notably in innovative seat designs where weight reduction and structural integration are paramount concerns. Modern composite seat structures can be designed to integrate multiple functions, including structural support, comfort features, and safety systems, all while maintaining the lightweight characteristics essential for aircraft performance.

Digital image correlation↗

DOE OSTI2026

First Wall Design of a Tokamak Pilot Plant Using a Monte Carlo Model for 3-D Heat Flux Deposition

We present a method for calculating the heat fluxes deposited on nonaxisymmetric tokamak first wall components, allowing for a first-of-its-kind model for power handling in the tokamak far scrape-off layer (SOL). The DIV3D Monte Carlo model features strict global power conservation and can calculate the finite cross-field plasma transport into magnetically-shadowed regions, which is significant when dealing with meter-scale shadows introduced by components such as poloidal limiters or antennas. As a case study, we apply the DIV3D model to inform the distribution of first wall poloidal limiters in an ARC-class reactor device. We demonstrate that discrete protection limiters can efficiently reduce peak heat fluxes on recessed breeder wall components in the presence of significant far-SOL plasma fluxes. By varying the toroidal periodicity and radial standoff depth of the limiters, we demonstrate one of the tradeoffs that must be considered in first wall design: more limiters provide greater protection, but at the cost of reduced breeding performance. We also present the impact that radial misalignments between limiters would have on first wall power loading.

Monte Carlo methods↗

NASA NTRS2026

Gateway Element and Payload Materials Outgassing Analyses: HALO, HERMES, and ERSA

Gateway was intended to be humanity’s first space station around the Moon, but its development has been paused as the National Aeronautics and Space Administration (NASA) shifts focus to achieving the United States’ National Space Policy goals. Instead of an orbiting lunar outpost, NASA will now pursue the development of a lunar surface base to support a sustained human presence on the Moon. Before the program’s pause, Gateway’s Induced Environments team worked to ensure payloads and elements (i.e., modules) complied with induced environment requirements. Methods developed and insights gained from this work will have applicability to NASA’s Moon Base and the potential repurposing of Gateway elements and payloads, as well as to induced environments modeling for future space stations. The Gateway program’s induced environment included molecular contamination, electric thruster plume sputter and redeposition, and lunar dust transfer from the Human Landing System (HLS). Primary sources of external molecular contamination included materials outgassing, chemical thruster plume contamination, and vacuum venting. The focus of this paper will be on element- and payload-level materials outgassing analyses performed for Gateway Configuration 1, extending the previously-developed framework for Gateway system-level external molecular contamination modeling. Gateway Configuration 1 consisted of the Power and Propulsion Element (PPE) and the Habitation and Logistics Outpost (HALO). It also included payloads like the European Radiation Sensor Array (ERSA) attached to PPE and the Heliophysics Environmental and Radiation Measurement Experiment Suite (HERMES) attached to HALO. The element- and payload-level analyses to be introduced in this paper for HALO, HERMES, and ERSA enabled high-fidelity descriptions of Gateway’s external molecular contamination environment. Approaches to geometric modeling, meshing, outgassing rate assignment, molecular transport modeling, and analysis methodology will be presented. Element and payload contaminant deposition onto sensitive Gateway receiver surfaces will be summarized and results compared to induced environment requirements. While these results incorporate refinements made over the course of the program, they were not intended to be final. Therefore, modeling assumptions and inputs, potential improvements, and lessons-learned will be documented to inform future work on Moon Base, repurposed elements and payloads, and other space stations.

Gateway↗

DOE OSTI2024

Determining Stellar Elemental Abundances from DESI Spectra with the Data-driven Payne

Abstract Stellar abundances for a large number of stars provide key information for the study of Galactic formation history. Large spectroscopic surveys such as the Dark Energy Spectroscopic Instrument (DESI) and LAMOST take median-to-low-resolution (R≲ 5000) spectra in the full optical wavelength range for millions of stars. However, the line-blending effect in these spectra causes great challenges for elemental abundance determination. Here we employDD-Payne, a data-driven method regularized by differential spectra from stellar physical models, to the DESI early data release spectra for stellar abundance determination. Our implementation delivers 15 labels, including effective temperatureT eff , surface gravity log g , microturbulence velocityv mic , and the abundances for 12 individual elements, namely C, N, O, Mg, Al, Si, Ca, Ti, Cr, Mn, Fe, and Ni. Given a spectral signal-to-noise ratio of 100 per pixel, the internal precisions of the label estimates are about 20 K forT eff , 0.05 dex for log g , and 0.05 dex for most elemental abundances. These results agree with the theoretical limits from the Crámer–Rao bound calculation within a factor of 2. The majority of the accreted halo stars contributed by the Gaia–Enceladus–Sausage are discernible from the disk and in situ halo populations in the resultant [Mg/Fe]–[Fe/H] and [Al/Fe]–[Fe/H] abundance spaces. We also provide distance and orbital parameters for the sample stars, which spread over a distance out to ∼100 kpc. The DESI sample has a significantly higher fraction of distant (or metal-poor) stars than the other existing spectroscopic surveys, making it a powerful data set for studying the Galactic outskirts. The catalog is publicly available.

Astronomy & Astrophysics↗

DOE OSTI2026

A transport-based framework for solidification cracking in Ni-based superalloys processed by laser powder bed fusion

Solidification cracking remains a persistent barrier to laser powder bed fusion (LPBF) processing of solid-solution-strengthened (SSS) Ni-based superalloys and is commonly attributed to carbide formation, liquation-type failure, or elemental segregation. In this work, the cracking behavior of three SSS Ni-based superalloys—Inconel 625 (625), Inconel 617 (617), and Haynes 230 (230)—is examined under comparable LPBF conditions to evaluate the origins of their cracking susceptibility. Carbides were present in both 625 and 230, yet cracking behavior differed significantly. MC-type carbides in 625 remain discrete and preserve liquid connectivity, whereas Cr-rich M23C6 carbides in 230 form at cellular triple junctions that bottleneck the interdendritic liquid network, influencing cracking through liquid connectivity rather than as an independent cause. Alloy 617 exhibited anomalous segregation without resolvable secondary phases yet cracked mildly, indicating that extensive carbide formation is not a prerequisite for cracking. To rationalize these observations, a transport-limited framework is proposed in which cracking occurs when terminal-liquid feeding cannot accommodate solidification-induced strain; among the mechanisms evaluated, this framework is most consistent with the physical constraints imposed by LPBF solidification. Transport-based crack-susceptibility indexes incorporating permeability, viscosity, and solidification kinetics distinguish the alloy hierarchy; permeability ratios Krat increased from approximately 7 (625) to 11.5 (617) and 14.5 (230), with corresponding increases in the μ-weighted mushy zone risk index (CSIMZRM). These results are consistent with transport-limited liquid feeding as a rate-limiting contributor to solidification cracking in LPBF-processed SSS Ni-based superalloys, providing a physically grounded basis for alloy and process design.

Hyer, Holden [ORNL] (ORCID:0000000343915561)↗

DOE OSTI2025

Polaron catastrophe within quantum acoustics

The quantum acoustic framework has recently emerged as a nonperturbative, coherent approach to electron–lattice interactions, uncovering rich physics often obscured by perturbative methods with incoherent scattering events. Here, we model the strongly coupled dynamics of electrons and acoustic lattice vibrations within this framework, representing lattice vibrations as coherent states and electrons as quantum wave packets, in a manner distinctively different from tight-binding or discrete hopping-based approaches. We derive and numerically implement electron backaction on the lattice, providing both visual and quantitative insights into electron wave packet evolution and the formation of acoustic polarons. We investigate polaron binding energies across varying material parameters and compute key observables—including mean square displacement, kinetic energy, potential energy, and vibrational energy—over time. Our findings reveal the conditions that favor polaron formation, which is enhanced by low temperatures, high deformation potential constants, slow sound velocities, and high effective masses. Additionally, we explore the impact of external electric and magnetic fields, showing that while polaron formation remains robust under moderate fields, it is weakly suppressed at higher field strengths. These results deepen our understanding of polaron dynamics and pave the way for future studies into nontrivial transport behavior in quantum materials.

Science & Technology - Other Topics↗

DOE OSTI2026

A hybrid Monte Carlo-deterministic second moment method with efficient variance reduction

In this work, we present a hybrid method that combines Monte Carlo with deterministic finite element methods to solve a linear Boltzmann transport equation. Our hybrid method runs orders of magnitude faster than Monte Carlo, without sacrificing accuracy, for a proxy problem from radiative transfer that contains both optically-thick and optically-thin material. We believe that this is the first demonstration of a hybrid Second Moment Method in more than one spatial dimension, the first to consider more than one material, and the first to use variance reduction. Our variance reduction approach arises from an asymptotic analysis in which we show that the magnitude of the scattering source grows without bound. We transform the problem to compute the deviation of the radiation intensity from isotropy. The magnitude of the source in the transformed problem is bounded, and the quality of the hybrid method solution is dramatically improved by a substantial reduction in the variance.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗