Engineering PapersSearch

SEARCH · Engineering Papers

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.

700 records · Page 7

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

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

Analysis of an irradiated uranium sample for source attribution without chemical separation using microplasma ionization and ultrahigh resolution mass spectrometry

The use of element isotope ratios has great potential in not only determining the reactor type used to produce plutonium (Pu) but also in determining the burnup and the time since irradiation. While a powerful nuclear forensic technique, determining element isotope ratios is complicated by severe isobaric interferences when performed on typical inductively coupled plasma mass spectrometers. Such analyses require extensive chemical separations prior to analysis to alleviate the inter-elemental isobars. Ultrahigh mass resolution spectrometry provides a potential alternative, greatly reducing the complexity of sample preparation and turnaround times for these critical measurements. To demonstrate the power of the approach, a sample of irradiated, depleted uranium was analyzed with the liquid sampling—atmospheric pressure glow discharge ion source coupled to an Orbitrap mass spectrometer. The Orbitrap is augmented with an external data acquisition system, Spectroswiss’s FTMS-Booster X2T, allowing collection of extended ion transients, providing higher mass resolution. In using this approach, the 150 Sm/ 149 Sm and 152 Sm/ 149 Sm isotope ratios were found to be within 20% of predicted values without any chemical separations and without mass bias corrections. In addition, the 240 Pu/ 239 Pu isotope ratio was determined, free from the 238 UH + interferences common to the ICP-MS platforms, while at the same time allowing for the determination of U isotopic signatures. While these demonstrative results are from a single sample, the advantages of the microplasma/ultrahigh mass resolution approach to intra-element isotope ratio determinations are clear.

Fuel burnup

Decadal-Scale Hydrothermal Alteration of Basalt: Implications for Long-Term Reactivity in Subsurface Technologies

Global demands to secure critical mineral supply chains and managing water resources require innovative resource management strategies. Subsurface basalt reservoirs offer a dual-purpose application, where CO2 enhanced critical metal recovery can be coupled with water remediation strategies for the extraction of economically important and environmentally relevant elements. This study utilizes static batch experiments conducted on Columbia River Basalt Group samples to evaluate long-term fluid–rock interactions relevant to coupled subsurface applications. Samples were reacted under elevated temperature and pressure conditions, for 296 and 4,061 days, providing a unique dataset to assess geochemical evolution, element partitioning, and aqueous chemistry over decadal timescales. The alteration was dominated by silicate, aluminosilicate, and oxide phases. Progressive hydrothermal alteration increased the cation exchange capacity through the formation of secondary minerals (e.g. clays and zeolites), influencing solute retention, ion exchange behavior, and aqueous chemistry. Grain-scale and bulk analyses revealed selective cation release and retention patterns, establishing a mechanistic baseline for predicting element mobility. These geochemical constraints inform the prospective use of CO2-enhanced metal recovery for critical resources (e.g. Co, Ni, Mn, Mg, Ti, Fe) and associated water-quality impacts, providing a foundation for the design and management of engineered fluid–rock systems in basalt reservoirs.

Gibson, Marie N.

A novel digital lifecycle for Material‐Process‐Microstructure‐Performance relationships of thermoplastic olefins foams manufactured via supercritical fluid assisted foam injection molding

Abstract This research significantly enhances the applicability of thermoplastic olefins (TPOs) in the automotive industry using supercritical N 2 as a physical foaming agent, effectively addressing the limitations of traditional chemical agents. It merges experimental results with simulations to establish detailed material‐process‐microstructure‐performance (MP2) relationships, targeting 5–20% weight reductions. This innovative approach labeled digital lifecycle (DLC) helps accurately predict tensile, flexural, and impact properties based on the foam microstructure, along with experimentally demonstrating improved paintability. The study combines process simulations with finite element models to develop a comprehensive digital model for accurately predicting mechanical properties. Our findings demonstrate a strong correlation between simulated and experimental data, with about a 5% error across various weight reduction targets, marking significant improvements over existing analytical models. This research highlights the efficacy of physical foaming agents in TPO enhancement and emphasizes the importance of integrating experimental and simulation methods to capture the underlying foaming mechanism to establish material‐process‐microstructure‐performance (MP2) relationships. Highlights Establishes a material‐process‐microstructure‐performance (MP2) for TPO foams Sustainably produces TPO foams using supercritical (ScF) N 2 with 20% lightweighting Shows enhanced paintability for TPO foam improved surface aesthetics Digital lifecycle (DLC) that predicts both foam microstructure and properties DLC maps process effects & microstructure onto FEA mesh for precise prediction

Engineering

Spinel high-entropy oxides (FeNiCrMnZnX) 3 O 4 (X = Al, mg) as anode materials for high-performance lithium-ion batteries

To address the high cost, cobalt dependency, and resource constraints typical of conventional high-entropy oxide (HEO) anodes, this study reports the successful synthesis of two Co-free, six-component spinel-type HEOs(FeNiCrMnZnAl) 3 O 4 (HEO-Al) and (FeNiCrMnZnMg) 3 O 4 (HEO-Mg), via a sol-gel method. The distinct effects of Al 3+ and Mg 2+ incorporation on the electrochemical performance and lithium storage kinetics were systematically investigated. XRD, Raman, and TEM characterizations confirm that both materials possess a pure spinel phase, uniform particle size, and homogeneous elemental distribution. Notably, electrochemical evaluations reveal that HEO-Al delivers a superior reversible capacity of 480.7 mAh g −1 after 100 cycles at 0.1 A g −1 , and maintains 354.5 mAh g −1 after 1000 long-term cycles at 1 A g −1 , significantly outperforming HEO-Mg. Kinetic analysis indicates that HEO-Al exhibits lower charge transfer resistance, a higher Li + diffusion coefficient, and a pseudocapacitive contribution of up to 82%. Furthermore, these findings demonstrate that Al substitution effectively optimizes the structural stability and lithium storage kinetics of Co-free HEOs, providing a viable strategy for designing low-cost, highly stable HEO anode systems.

Anode materials

Battery case shear

Hand operated shear removes a battery case without disturbing the internal components which are to be tested. It consists of three tool-steel elements, the cutter blade, and a hand lever that provides the mechanical advantage required to cut steel.

Patro, S.

Aragonite saturation horizon variability along North Pacific seamounts and implications for deep-sea coral reefs

The 2014 discovery of living deep-sea coral reefs along the Northwest Hawaiian Islands (NWHI) and lower Emperor Seamount Chain (ESC), despite the North Pacific’s shallow aragonite saturation horizon (ASH) and high CaCO3 dissolution rates, underscores the need to understand the local seawater chemistry that allows these reefs to persist. We investigated seawater carbonate chemistry using discrete samples along the NWHI and ESC from two cruises ~1 year apart (08/26/21 – 09/26/21, 09/09/22 – 10/24/22). Across the two cruises, ASH depth difference ranged from 15 to 77 m. Since the Pacific ASH shoals by 1–2 m yr?¹, this long-term trend cannot explain the magnitude of ASH change observed between cruises. Temperature-salinity plots similarly indicate no major shifts in intermediate water masses that could account for these changes. Instead, ASH depth variability was primarily governed by localized biogeochemical processes, namely changes in intermediate water respiration and CaCO3 dissolution. Indicators for dissolution (TA*) and respiration (AOU) suggest changes in ASH depth were driven by changes in dissolution at the northern- and southern-most sites, whereas respiration exerted stronger control at central sites. Combining 2021 and 2022 data with data from 2014 – 2019 (excluding 2018) revealed high interannual ASH variability, by as much as >200 m at one site. Deep-sea coral reefs across the NWHI and ESC currently reside close to the ASH depth and likely experience seasonal and interannual shifts between under- and supersaturation. As ocean acidification progresses however, persistent exposure to undersaturated seawater could further threaten these critical deep-sea ecosystems.

coral reefs

Numerical Investigation of Fluid Flow and Space Charge in Liquid Argon Time Projection Chamber (LArTPC) Detectors

Overview This project focused on developing a high-fidelity numerical framework to simulate the multiphysics environment within Liquid Argon Time Projection Chamber (LArTPC) detectors. The primary objective was to characterize the complex interplay between ion transport, background fluid dynamics, and electric field distortions—a critical factor for the calibration and sensitivity of next-generation High Energy Physics experiments, such as DUNE. Technical Achievements The research successfully yielded a hybrid numerical space-charge solver utilizing a Cell-Centered Finite Volume Method (FVM) for ion transport coupled with a Finite Element Method (FEM) for electric potential. Key accomplishments include: • Verification & Validation: The 3-D solver was rigorously verified against 1-D analytical solutions, demonstrating high numerical accuracy in predicting space-charge-induced field deviations. • Field Distortion Analysis: 3D simulations revealed that space charge effects introduce significant non-uniformities in the electric field. Critically, the research identified that background LAr flow velocities, when comparable to ion drift velocities, markedly exacerbate these distortions. • Technology Transfer: The resulting source code and comprehensive user manuals were successfully transferred to collaborators at Fermilab, providing a portable computational tool for the broader scientific community. Challenges and Future Directions While the space-charge solver achieved all performance metrics, the integrated fluid dynamics modeling encountered convergence challenges stemming from the extreme 200-fold disparity in length scales between the detector's 37 mm inlet pipes and the 8-meter global domain. To address this, the project has identified a clear technical pivot toward Hierarchical Geometric Adaptive Mesh Refinement (HG-AMR). By implementing an h-type refinement strategy with hanging nodes, future iterations of this solver will be capable of resolving localized high-gradient inlet flows without the prohibitive computational costs of regular grids. This advancement, combined with data-driven uncertainty quantification based on MicroBooNE-style calibration, will enable the precise modeling of detector responses in large-scale cryogenic environments where direct measurement remains difficult. Impact The computational tools developed under this award provide a foundation for enhancing the energy resolution and spatial reconstruction of noble liquid detectors. By bridging the gap between theoretical fluid dynamics and experimental field calibration, this work supports the DOE’s mission to advance the frontiers of neutrino physics and dark matter detection.

42 ENGINEERING

Multiscale Modeling of Fracture Strength in Fibrous Thermal Protection System Materials

This work presents a multiscale modeling approach to predict the fracture strength of fibrous Thermal Protection System (TPS) materials. The model assumes that system failure is initiated at the joints between individual fibers. We investigated three distinct TPS compositions: amorphous silica, alumina and aluminosilicate fibers. Molecular dynamics (MD) simulations were employed to determine the fracture strength values of these fiber joints for both material systems. These fracture strength values were then integrated into simulations of 3D randomly populated fiber structures, where tensile load transfer occurs through the fiber joints. These microscale properties are upscaled through a renormalization approach [1] to predict macroscale tensile strength of 3D random fiber networks, accounting for joint-dominated failure and effective load-bearing area. The study concludes by demonstrating the resulting strength variation as a function of material composition, fiber density, and morphology. We also show validation of results by comparing them against explicit fiber finite element (FE) modeling [2] where fiber joint fracture is represented by cohesive elements.

Jaehyun Cho

General Purpose Data-Driven System Monitoring for Space Operations

Modern space propulsion and exploration system designs are becoming increasingly sophisticated and complex. Determining the health state of these systems using traditional methods is becoming more difficult as the number of sensors and component interactions grows. Data-driven monitoring techniques have been developed to address these issues by analyzing system operations data to automatically characterize normal system behavior. The Inductive Monitoring System (IMS) is a data-driven system health monitoring software tool that has been successfully applied to several aerospace applications. IMS uses a data mining technique called clustering to analyze archived system data and characterize normal interactions between parameters. This characterization, or model, of nominal operation is stored in a knowledge base that can be used for real-time system monitoring or for analysis of archived events. Ongoing and developing IMS space operations applications include International Space Station flight control, spacecraft vehicle system health management, launch vehicle ground operations, and fleet supportability. As a common thread of discussion this paper will employ the evolution of the IMS data-driven technique as related to several Integrated Systems Health Management (ISHM) elements. Thematically, the projects listed will be used as case studies. The maturation of IMS via projects where it has been deployed or is currently being integrated to aid in fault detection will be described. The paper will also explain how IMS can be used to complement a suite of other ISHM tools, providing initial fault detection support for diagnosis and recovery

Space Propulsion

Circumventing data imbalance in magnetic ground state data for magnetic moment predictions

Abstract Magnetic materials play a crucial role in the transition to more sustainable forms of energy and electric vehicles. There is an anticipated shortage in magnetic materials in the future, and as a result there is an urgent need to discover and design new magnetic materials. Computational magnetic material design using density functional theory is daunting because of the challenge in identifying magnetic ground states from a combinatorially large set of possibilities. Machine learning offers a path forward by enabling efficient surrogate models that can more readily enumerate these states, but there is a dearth of training data available, and what is available tends to be imbalanced with too much non-magnetic data. In this work we show that the discrete and previously tackled data imbalance that exists at the level of the magnetic ordering leads to an imbalanced continuous distribution with many zeros when the data is unraveled at the atomic magnetic moment level, which subsequently leads to models with low accuracy for magnetic properties. We mitigate this by using a two-part model framework. Our scheme is able to classify atoms into magnetic and non-magnetic with an F1 score and Matthew’s correlation coefficient (MCC) of ~91% and then to provide an implicit embedding representation that maps directly onto the magnitude of the magnetic moment with a mean absolute error of 0.1 μ B . Beyond screening for new magnetic materials, we demonstrate an additional practical use case of our scheme: the provision of good initial guesses for magnetic moments in first-principles electronic relaxations. Such initialization is shown to lead to faster convergence to configurations that lie closer to the ground state.

Computer Science

Layout and Cabling Considerations for a Large Communications Antenna Array

Layout considerations for a large deep space communications antenna array are discussed. A novel fractal geometry for the antenna layout is described that provides optimal packing of antenna elements, efficient cable routing, and logical division of the array into identical sub-arrays.

R T Logan, Jr.

General Purpose Data-Driven System Monitoring for Space Operations

Modern space propulsion and exploration system designs are becoming increasingly sophisticated and complex. Determining the health state of these systems using traditional methods is becoming more difficult as the number of sensors and component interactions grows. Data-driven monitoring techniques have been developed to address these issues by analyzing system operations data to automatically characterize normal system behavior. The Inductive Monitoring System (IMS) is a data-driven system health monitoring software tool that has been successfully applied to several aerospace applications. IMS uses a data mining technique called clustering to analyze archived system data and characterize normal interactions between parameters. This characterization, or model, of nominal operation is stored in a knowledge base that can be used for real-time system monitoring or for analysis of archived events. Ongoing and developing IMS space operations applications include International Space Station flight control, satellite vehicle system health management, launch vehicle ground operations, and fleet supportability. As a common thread of discussion this paper will employ the evolution of the IMS data-driven technique as related to several Integrated Systems Health Management (ISHM) elements. Thematically, the projects listed will be used as case studies. The maturation of IMS via projects where it has been deployed, or is currently being integrated to aid in fault detection will be described. The paper will also explain how IMS can be used to complement a suite of other ISHM tools, providing initial fault detection support for diagnosis and recovery.

Satellites

A Survey of Lunar Rock Types and Comparison of the Crusts of Earth and Moon

The principal known types of lunar rocks are briefly reviewed, and their chemical relationships discussed. In the suite of low-KREEP highland rocks, Fe/(Fe + Mg) in the normative mafic minerals increases and the albite content of normative plagioclase decreases as the total amount of normative plagioclase increases, the opposite of the trend predicted by the Bowen reaction principle. Lunar highland samples analyzed are uniformly distributed in this sequence, in which normative plagioclase contents range from ~ 40 percent to ~ 100 percent. The distribution of compositions of rocks from terrestrial layered mafic intrusives is substantially different: here the analyses fall in several discrete clusters (anorthositic rocks, norites, granophyres and ferrogabbros, ultramafics), and the chemical trends noted above are not reproduced. It is suggested that the observed trends in lunar highland rocks could be produced by crystal fractionation in a deep global surface magma system if (1) plagioclase tended to float, upon crystallization, and (2) the magma was kept agitated and well mixed (probably by thermal convection) until crystallization was far advanced and relatively little residual liquid was left. When such a system was finally immobilized, the Fe-, Na-rich residual liquid would produce Fe-rich mafic minerals in the upper levels of the system, but could not much alter the composition of abundant calcic plagioclase. Conversely, the same liquid would produce sodic plagioclase deep in the sequence, but could not much alter the composition of abundant magnesian mafic minerals. After the crustal system solidified, but before extensive cooling had developed a thick, strong lithosphere, mantle convection was able to draw portions of the lunar anorthositic crust down into the mantle in a manner analogous to the present-day behavior of the terrestrial mantle and crust. At depth, the crustal material was heated; KREEP-rich norite was extracted by partial melting and erupted at the surface as a lava, analogous to terrestrial andesite eruptions.

John A Wood

Sierra/SD - Users Manual - 5.32

Sierra/SD provides a massively parallel implementation of structural dynamics finite element analysis, required for high-fidelity, validated models used in modal, vibration, static and shock analysis. This document provides a user's guide to the input for Sierra/SD. Details of input specifications for the different solution types, output options, element types and parameters are included.

Beale, Dagny Marie [Sandia National Laboratories (

Parabolic Deformable Mirror Project Development and Device Characterization

Existing state-of-the-art wavefront sensing and control setups for the imaging of exoplanets utilize two high-actuator count flat deformable mirrors (DMs) to control both phase and amplitude aberrations. One of the objectives of the Exoplanet Spectroscopy (ExoSpec) project at Goddard Space Flight Center (GSFC) is to improve upon this regime through the development of parabolic deformable mirrors (PDMs), which can address several of the key limitations in existing setups. Specifically, the need for flat DMs to be located in conjugate pupil planes drives unique packaging challenges that can be overcome by allowing the off-axis imaging elements to be controllable as well. In addition, simulations show that the use of deformable powered elements will increase the overall controllable bandwidth, up to 35% over 5-12 lambda/D. To pursue this effort, the ExoSpec group has begun testing a first-generation parabolic DM, which is already informing a second generation design. This paper will provide a comprehensive update on its characterization, including influence function measurements, linearity analysis, closed-loop wavefront control performance, and thermal stability assessment.

spectroscopy

Used Nuclear Fuel Management Using the Next Generation System Analysis Model

The U.S. Department of Energy (DOE) is leading the National effort to manage the back end of the nuclear fuel cycle, encompassing the safe transportation, storage/staging, and/or eventual disposal of used nuclear fuel (UNF) and high-level radioactive waste. The Next Generation System Analysis Model (NGSAM) is DOE’s discrete-event, agent-based simulation tool designed to model the full life cycle of UNF from reactor discharge to final disposal. NGSAM supports the DOE Office of Spent Fuel and High-Level Waste Disposition by enabling a detailed, scenario-based analysis of logistics, infrastructure, and shipping strategies. NGSAM replaces legacy models with a modern, flexible platform built on Repast Simphony and enhanced by the Process Analysis Tool. NGSAM simulates the movement and interaction of individual fuel assemblies with system components such as canisters, casks, railcars, and facilities. The model integrates with the Java Transportation Operations Model to plan and execute transportation scenarios, supporting both constrained and unconstrained resource allocation. Key features include customizable allocation and acceptance algorithms, detailed facility-level operations, and a Quick Edit tool for rapid scenario adjustments. NGSAM supports multimodal transportation modeling (e.g. rail, road, barge) and provides comprehensive cost, schedule, and infrastructure data. NGSAM utilizes data from sources such as DOE’s STANDARDS UNF database and DOE’s Stakeholder Tool for Assessing Radioactive Transportation, while also allowing user-defined inputs for scenario customization. NGSAM enables stakeholders to evaluate complex UNF management strategies, assess system performance under varying assumptions, and inform decision making for future infrastructure investments. Its modular architecture and integration with other Integrated Waste Management System tools make it a critical asset for planning the safe and efficient disposition of the Nation’s growing UNF inventory.

Craig, Brian [Argonne National Laboratory (ANL)]