SEARCH · Engineering Papers

Results for “DEM simulation”

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.

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

NASA NTRS2010

A Review of Discrete Element Method (DEM) Particle Shapes and Size Distributions for Lunar Soil

As part of ongoing efforts to develop models of lunar soil mechanics, this report reviews two topics that are important to discrete element method (DEM) modeling the behavior of soils (such as lunar soils): (1) methods of modeling particle shapes and (2) analytical representations of particle size distribution. The choice of particle shape complexity is driven primarily by opposing tradeoffs with total number of particles, computer memory, and total simulation computer processing time. The choice is also dependent on available DEM software capabilities. For example, PFC2D/PFC3D and EDEM support clustering of spheres; MIMES incorporates superquadric particle shapes; and BLOKS3D provides polyhedra shapes. Most commercial and custom DEM software supports some type of complex particle shape beyond the standard sphere. Convex polyhedra, clusters of spheres and single parametric particle shapes such as the ellipsoid, polyellipsoid, and superquadric, are all motivated by the desire to introduce asymmetry into the particle shape, as well as edges and corners, in order to better simulate actual granular particle shapes and behavior. An empirical particle size distribution (PSD) formula is shown to fit desert sand data from Bagnold. Particle size data of JSC-1a obtained from a fine particle analyzer at the NASA Kennedy Space Center is also fitted to a similar empirical PSD function.

Lane, John E.↗

NASA NTRS2013

PowderSim: Lagrangian Discrete and Mesh-Free Continuum Simulation Code for Cohesive Soils

PowderSim is a calculation tool that combines a discrete-element method (DEM) module, including calibrated interparticle-interaction relationships, with a mesh-free, continuum, SPH (smoothed-particle hydrodynamics) based module that utilizes enhanced, calibrated, constitutive models capable of mimicking both large deformations and the flow behavior of regolith simulants and lunar regolith under conditions anticipated during in situ resource utilization (ISRU) operations. The major innovation introduced in PowderSim is to use a mesh-free method (SPH-based) with a calibrated and slightly modified critical-state soil mechanics constitutive model to extend the ability of the simulation tool to also address full-scale engineering systems in the continuum sense. The PowderSim software maintains the ability to address particle-scale problems, like size segregation, in selected regions with a traditional DEM module, which has improved contact physics and electrostatic interaction models.

Johnson, Scott↗

NASA NTRS2014

On Characterizing Particle Shape

It is well known that particle shape affects flow characteristics of granular materials, as well as a variety of other solids processing issues such as compaction, rheology, filtration and other two-phase flow problems. The impact of shape crosses many diverse and commercially important applications, including pharmaceuticals, civil engineering, metallurgy, health, and food processing. Two applications studied here include the dry solids flow of lunar simulants (e.g. JSC-1, NU-LHT-2M, OB-1), and the flow properties of wet concrete, including final compressive strength. A multi-dimensional generalized, engineering method to quantitatively characterize particle shapes has been developed, applicable to both single particle orientation and multi-particle assemblies. The two-dimension, three dimension inversion problem is also treated, and the application of these methods to DEM model particles will be discussed. In the case of lunar simulants, flow properties of six lunar simulants have been measured, and the impact of particle shape on flowability - as characterized by the shape method developed here -- is discussed, especially in the context of three simulants of similar size range. In the context of concrete processing, concrete construction is a major contributor to greenhouse gas production, of which the major contributor is cement binding loading. Any optimization in concrete rheology and packing that can reduce cement loading and improve strength loading can also reduce currently required construction safety factors. The characterization approach here is also demonstrated for the impact of rock aggregate shape on concrete slump rheology and dry compressive strength.

Ennis, Bryan J.↗

NASA NTRS2019

Software Development and Testing Support for the Avionics Systems Telemetry Tool Suite

The Customer Avionics Interface Development and Analysis (CAIDA) team provides modeling and simulation software for the verification of the Launch Control System (LCS). In late 2015, CAIDA began development of the Data Exchange Message (DEM) Generator (DEMGen) and development of the DEM Modifier (DEMMod) in early 2017. DEMGen is able to mimic the telemetry streams normally sent by different simulators allowing increased user control and more telemetry instances. DEMMod takes in a telemetry stream and allows the modification of entire packets or simple measurement values. Together, these tools provide the capabilities to simulate and test complex launch scenarios to ensure LCS is fully prepared for any anomalous behavior during an actual launch. In early 2018, CAIDA began a new project called the CAIDA Advanced Telemetry Tool (CATT), which combined the code and documentation of DEMGen and DEMMod. CATT now allows the integration of new tools into a singular program providing compact and easy access. In the fall of 2018, the author worked with another intern, Antonio Negron, to develop several new CATT features along with proper documentation and unit testing. In the following spring, the author continued that work by integrating the features into CATT, solidifying the documentation, and expanding the testing to cover more cases. The outcome of both semesters includes an inspection engine for automated packet analysis, a packet-masking feature, and a multi-measurement sequencing feature.

Thomas, Taylor Walter↗

NASA NTRS2020

Analysis of Sample Acquisition Dynamics Using Discrete Element Method

The analysis presented in this paper is conducted in the framework of the Ocean Worlds Autonomy Testbed for Exploration Research and Simulation (OceanWATERS) project, currently under development at NASA Ames Research Center. OceanWATERS aims at designing a simulation environment which allows for testing autonomy of scientific lander missions to the icy moons of our solar system. Mainly focused on reproducing the end effector interaction with the inherent terrain, this paper introduces a novel discrete element method (DEM)-based approach to determine forces and torques acting on the lander’s scoop during the sample acquisition process. An accurate force feedback from the terrain on the scoop is required by fault-detection and autonomous decision-making algorithms to identify when the requested torque on the robotic arm’s joints exceeds the maximum available torque. Knowledge of the terrain force feedback significantly helps evaluating the arm’s links structural properties and properly selecting actuators for the joints. Models available in literature constitute a partial representation of the dynamics of the interaction. As an example, Balovnev derived an analytical expression of the vertical and horizontal force acting on a bucket while collecting a sample as a function of its geometry and velocity, soil parameters and reached depth. Although the model represents an adequate approximation of the two force components, it ignores the direction orthogonal to the scoop motion and neglects the torque. This work relies on DEM analysis to compensate for analytical models’ deficiencies and inaccuracies, i. e. provide force and torque 3D vectors, defined in the moving reference (body) frame attached to the scoop, at each instant of the sample collection process. Results from the first presented analysis relate to the specific OceanWATERS sampling strategy, which consists of collecting the sample through five consecutive passes with increasing depth, each pass following the same circularlinear- circular trajectory. Data is collected given a specific scoop design interacting with two types of bulk materials, which may characterize the surface of icy planetary bodies: snow and ice. Although specifically concerned with the OceanWATERS design, this first analysis provides the expected force trends for similar sampling strategies and allows to deduce phenomenological information about the general scooping process. In order to further instruct the community on the use of DEM tools as a solution to the sampling collection problem, two more analyses have been carried out, mainly focused on reducing the DEM computation time, which increases with a decrease in particle size. After running a set of identical simulations, where the only changing parameter is the size of the spherical particle, it is observed that the resulting force trajectories, starting from a given particle size, converge to the true trend. It is deducible that a further decrease in size yields negligible improvements in the accuracy, while it sensibly increases computation time. A final analysis aims at discussing limitations of approximating bulk material particles having a complex shape, e. g. ice fragments, with spheres, by comparing force trends resulting in the two cases for the same simulation scenario.

Catanoso, Damiana↗

NASA NTRS2009

Physics of Reconnection and MMS Mission

Reconnection is the most important process driving the Earth's magnetosphere. Key to the success of the MMS science plan is the coupling of theory and observation. Determining the kinetic processes occurring in the diffusion region and physical parameters that control the rate of magnetic reconnection are among primary objectives of the MMS mission. Analysis of the role played by particle inertial effects in the diffusion region where the plasma is unmagnetized will be presented. The reconnection electric field in he diffusion region is supported primarily by particle non-gyrotropic effects. At the quasi-steady stage the reconnection electric field serves to accelerate and heat the incoming plasma population to maintain the current flow in the diffusion region the pressure balance. The primary mechanism controlling the dissipation in the vicinity of the reconnection site is incorporated into the fluid description in terms of non-gyrotropic corrections to the. induction and energy equations. The results of kinetic and fluid simulations illustrating the physics of magnetic reconnection will be presented. We will dem~:tistrate that kinetic nongyrotropic effects can significantly alter the global magnetosphere evolution and location of reconnection sites.

Kuznetsova, M. M.↗

NASA NTRS2009

Terminal Descent Sensor Simulation

Sulcata software simulates the operation of the Mars Science Laboratory (MSL) radar terminal descent sensor (TDS). The program models TDS radar antennas, RF hardware, and digital processing, as well as the physics of scattering from a coherent ground surface. This application is specific to this sensor and is flexible enough to handle end-to-end design validation. Sulcata is a high-fidelity simulation and is used for performance evaluation, anomaly resolution, and design validation. Within the trajectory frame, almost all internal vectors are represented in whatever coordinate system is used to represent platform position. The trajectory frame must be planet-fixed. The platform body frame is specified relative to arbitrary reference points relative to the platform (spacecraft or test vehicle). Its rotation is a function of time from the trajectory coordinate system specified via dynamics input (file for open loop, callback for closed loop). Orientation of the frame relative to the body is arbitrary, but constant over time. The TDS frame must have a constant rotation and translation from the platform body frame specified at run time. The DEM frame has an arbitrary, but time-constant, rotation and translation with respect to the simulation frame specified at run time. It has the same orientation as sigma0 frame, but is possibly translated. Surface sigma0 has the same arbitrary rotation and translation as DEM frame.

Chen, Curtis W.↗

NASA NTRS2008

Determining Trajectory of Triboelectrically Charged Particles, Using Discrete Element Modeling

The Kennedy Space Center (KSC) Electrostatics and Surface Physics Laboratory is participating in an Innovative Partnership Program (IPP) project with an industry partner to modify a commercial off-the-shelf simulation software product to treat the electrodynamics of particulate systems. Discrete element modeling (DEM) is a numerical technique that can track the dynamics of particle systems. This technique, which was introduced in 1979 for analysis of rock mechanics, was recently refined to include the contact force interaction of particles with arbitrary surfaces and moving machinery. In our work, we endeavor to incorporate electrostatic forces into the DEM calculations to enhance the fidelity of the software and its applicability to (1) particle processes, such as electrophotography, that are greatly affected by electrostatic forces, (2) grain and dust transport, and (3) the study of lunar and Martian regoliths.

Source record↗

NASA NTRS2010

Super Resolution Image Enhancement for a Flash Lidar: Back Projection Method

In this paper a new image processing technique for flash LIDAR data is presented as a potential tool to enable safe and precise spacecraft landings in future robotic or crewed lunar and planetary missions. Flash LIDARs can generate, in real-time, range data that can be interpreted as a 3-dimensional (3-D) image and transformed into a corresponding digital elevation map (DEM). The NASA Autonomous Landing and Hazard Avoidance (ALHAT) project is capitalizing on this new technology by developing, testing and analyzing flash LIDARs to detect hazardous terrain features such as craters, rocks, and slopes during the descent phase of spacecraft landings. Using a flash LIDAR for this application looks very promising, however through theoretical and simulation analysis the ALHAT team has determined that a single frame, or mosaic, of flash LIDAR data may not be sufficient to build a landing site DEM with acceptable spatial resolution, precision, size, or for a mosaic, in time, to meet current system requirements. One way to overcome this potential limitation is by enhancing the flash LIDAR output images. We propose a new super-resolution algorithm applicable to flash LIDAR range data that will create a DEM with sufficient accuracy, precision and size to meet current ALHAT requirements. The performance of our super-resolution algorithm is analyzed by processing data generated during a series of simulation runs by a high fidelity model of a flash LIDAR imaging a high resolution synthetic lunar elevation map. The flash LIDAR model is attached to a simulated spacecraft by a gimbal that points the LIDAR to a target landing site. For each simulation run, a sequence of flash LIDAR frames is recorded and processed as the spacecraft descends toward the landing site. Each run has a different trajectory profile with varying LIDAR look angles of the terrain. We process the output LIDAR frames using our SR algorithm and the results show that the achieved level of accuracy and precision of the SR generated landing site DEM is more than adequate for detecting hazardous terrain features and identifying safe areas.

Bulyshev, Alexander↗

DOE OSTI2025

Data-model files associated with the manuscript "Modeling the Effects of Wetland Restoration on Coastal Hydrology: A Case Study of Elkhorn Slough Watershed, California"

This package contains the data, simulation setups, notebooks and figures used in “Modeling the Effects of Wetland Restoration on Coastal Hydrology: A Case Study of Elkhorn Slough Watershed, California” (Xu et al., 2025). In this study, we selected Elkhorn Slough, a tidal estuary, in California, to investigate the impact of wetland restoration and sea level rise on coastal hydrology using the process-based coastal hydrologic model, Advanced Terrestrial Simulator (ATS), informed by site-specific data. We designed a novel modeling workflow for incorporating wetland restoration features into land cover and soil properties for the model parameterization. The validation results demonstrate a strong agreement between modeled and observed data. We studied the characteristics of coastal watershed hydrology, then focused on the surface water dynamics at two wetland sites within Elkhorn Slough, a reference site and a restored site. Our simulation results indicate that the restored site successfully maintains surface elevation, resulting in reduced surface inundation. We also examined the impact of wetland restoration under expected sea level rise over the next few decades. The low-lying Yampah Marsh, the reference site, is likely to be inundated due to future sea level rise when highest tides arrive; while a higher percentage of Hester Marsh, the restored site, would retain marsh vegetation in coming decades, regardless of tidal conditions. Our study provides important information for examining the outcome of restoration practices that include surface elevation in tidal wetlands under climate changes.Several files can be found from this data package.1. README.md: This file describes the title, journal, co-authors, abstract, repository structure and model version.2. Simulation_Setups.zip: The file contains the model configuration files (XML format) for ATS. 3. Notebooks.zip: The file contains the Jupyter notebooks for generating the pre- and post-restoration meshes and the meshes of future scenarios. 4. Figures.zip: The file contains the figures used in the manuscript.5. Data.zip: The file contains the data used to drive the model simulations, including watershed and wetlands boundaries, mesh files and references to additional datasets (e.g., meteorological forcing, tidal dataset, DEMs, land cover, soil properties). Also, it contains water level observations at the restored wetland.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

A High-Performance Discrete-Element Framework for Simulating Flow and Jamming of Moisture Bearing Biomass Feedstocks

We developed and verified a high-performance open-source discrete element method (DEM) solver with simultaneously-supported feedstock-specific interaction models, including bonded-sphere, liquid bridge, cohesion, and non-linear contact models. Our solver uses parallel data structures on hybrid central and graphics processing unit (CPU/GPU) architectures, with favorable strong scaling performance observed for large problem sizes comprised of (100 M particles), and 4X single-node GPU speedup. The particles for corn stover feedstock were conceptualized and calibrated based on experimental measurements and results. Sensitivity analyses demonstrate that the mass flow rate from a wedge hopper is governed primarily by moisture content, friction coefficient, and cohesion energy density. The model is used to reproduce experimentally observed hopper jamming results, highlighting that the experimental no-flow trends can only be achieved by using non-spherical particles, liquid bridge and cohesion models, highlighting the importance of using concurrent feedstock specialized models for the effective representation of biomass material handling problems.

bioenergy↗

NASA NTRS2024

Experimental Validation of Charged Lunar Dust Dynamics Simulants

The goal of this work is to experimentally verify the electrostatic interaction physics that have been incorporated into the discrete element method (DEM) modeling software package during the FY20 CIF project titled “Charged Particle Dynamics in the Lunar Environment.” Current state-of-the-art (SOA) granular gas dynamics models used to explain rocket plume impingement physics have not taken into account the natural or induced electrostatic environment of the lunar surface, nor the effect of charged regolith grains being present in the plume plasma. This work improved upon the SOA by adding long-range and contact inter-particle electrostatic interactions to a granular mechanics DEM modeling software package. Tribocharging mechanics between spherical particles has been experimentally examined along with the natural electrostatic plasma environments produced via emission from an electron gun and ultraviolet lamp.

Jay Phillips↗

NASA NTRS1993

The thermal structure of the magnetized solar transition region

The detailed thermal structure of the magnetized solar transition region, as measured by its differential emission measure DEM(T), is unknown. Proposals have been made that envision a significant lower-temperature contribution to the energy balance from cross-field (ion) heat flux. In this paper, we describe a self-consistent 2D MHD simulation (including the full effects of anisotropic thermal conduction) of a conceptual model due to Athay (1990). We display the detailed irregular thermal and magnetic structure of the transition region and demonstrate that the predicted DEM agrees with observations, particularly in the T less than 10 exp 5 K regime where previous theories had difficulty.

Mok, Y.↗

NASA NTRSDate not supplied

Particle-Based Fiber Models of Woven Materials for Earth Entry Thermal Protection

Applications requiring materials with layer-to-layer strength, from basketry to thermal protection systems, use interlaced, three-dimensional woven materials. NASA is developing and deploying woven material heat shields for missions, including Artemis I (3D-MAT for compression pads) and potentially Mars Sample Return - Earth Entry System. These materials are complex, hierarchical and must protect from extreme environments and phenomena, such as deformation, impact and high-enthalpy heating. Woven material performance depends on microstructure, damage and weave geometry. Therefore, fiber-specific models are needed to simulate fiber contacts within the weave hierarchical geometry (fiber to tow to yarn to weave) and the inherent directionality of fibers. Explicit fiber models can simulate how weave microstructure evolution affects thermal and mechanical properties. We parameterize a discrete element bonded particle models (DEM-BPM) of fibers to capture thermal and mechanical behavior within and between fibers, with bonded and contact forces, respectively. We study the proportion of heat transfer and stress via the contact network, fiber bonds and the weave geometry, for example, with respect to yarn warp-weft identity (whether it interlaces weave layers). Our results demonstrate the importance of explicit fiber modeling for connecting microstructure with thermal and mechanical properties.

woven material↗

NASA NTRSDate not supplied

Particle-Based Fiber Models of Woven Materials for Earth Entry Thermal Protection

woven material↗

NASA NTRS1989

Mars rover concept development

A structured study effort to develop an extensive, innovative set of mobility and navigation concepts for a planetary exploration vehicle along with the concomitant value system and evaluation tools is presented. A further objective is to submit these concepts to a rigorous, structured evaluation process to derive the most promising candidate systems. To support the evaluation process, a three-layer computer model of the Martian surface was developed, based on the 1/64 deg Digital Elevation Model (DEM) of Mars. Local surface roughness based on measured Martian slope distribution and power spectral density was superimposed on the DEM, and rocks based on Moore's distribution model were added. To assess performance, selected concepts were modeled using DADS, and simulations were run with the vehicle traversing the Martian surface model, including one-meter-high vertical steps and one-meter-wide crevasses. The design details of three promising candidate systems are presented, along with the discussion of their evolution with some recommendations.

Mctamaney, Louis S.↗

NASA NTRS1993

The differential emission measure of nested hot and cool magnetic loops

The detailed thermal structure of the magnetized solar transition region, as measured by its differential emission measure (DEM(T)), is poorly known. Building on the fact that the solar surface is strongly magnetized and thereby structured, proposals have been made that envision a significant lower-temperature contribution to the energy balance from (ion) heat flux across an arcade of different temperature loops. In this paper, we describe a self-consistent 2D MHD simulation, which includes the full thermal effects of parallel stability and anisotropic conduction, of a nested-loop model of the thermal and magnetic structure of the transition region. We then demonstrate that the predicted DEM agrees with observations in the conceptually elusive T less than 10 exp 5 K regime.

Van Hoven, G.↗

NASA NTRS2006

Forecasting of Storm-Surge Floods Using ADCIRC and Optimized DEMs

Increasing the accuracy of storm-surge flood forecasts is essential for improving preparedness for hurricanes and other severe storms and, in particular, for optimizing evacuation scenarios. An interactive database, developed by WorldWinds, Inc., contains atlases of storm-surge flood levels for the Louisiana/Mississippi gulf coast region. These atlases were developed to improve forecasting of flooding along the coastline and estuaries and in adjacent inland areas. Storm-surge heights depend on a complex interaction of several factors, including: storm size, central minimum pressure, forward speed of motion, bottom topography near the point of landfall, astronomical tides, and, most importantly, maximum wind speed. The information in the atlases was generated in over 100 computational simulations, partly by use of a parallel-processing version of the ADvanced CIRCulation (ADCIRC) model. ADCIRC is a nonlinear computational model of hydrodynamics, developed by the U.S. Army Corps of Engineers and the US Navy, as a family of two- and three-dimensional finite-element-based codes. It affords a capability for simulating tidal circulation and storm-surge propagation over very large computational domains, while simultaneously providing high-resolution output in areas of complex shoreline and bathymetry. The ADCIRC finite-element grid for this project covered the Gulf of Mexico and contiguous basins, extending into the deep Atlantic Ocean with progressively higher resolution approaching the study area. The advantage of using ADCIRC over other storm-surge models, such as SLOSH, is that input conditions can include all or part of wind stress, tides, wave stress, and river discharge, which serve to make the model output more accurate. To keep the computational load manageable, this work was conducted using only the wind stress, calculated by using historical data from Hurricane Camille, as the input condition for the model. Hurricane storm-surge simulations were performed on an eight-node Linux computer cluster. Each node contained dual 2-GHz processors, 2GB of memory, and a 40GB hard drive. The digital elevation model (DEM) for this region was specified using a combination of Navy data (over water), NOAA data (for the coastline), and optimized Interferometric Synthetic Aperture Radar data (over land). This high-resolution topographical data of the Mississippi coastal region provided the ADCIRC model with improved input with which to calculate improved storm-surge forecasts.

Valenti, Elizabeth↗