Search NASA⌕ Search

SEARCH · Search NASA

Results for “Discrete Element Modeling”

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.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Discrete element modeling of granular hopper flow of irregular-shaped deformable particles

Many natural and engineered granular materials have relatively deformable particles. Besides particle size and shape, particle deformability is another salient factor that significantly impacts the material’s flow behavior. Here, in this work, the flow of irregular-shaped deformable particles in a wedge-shaped hopper is investigated using discrete element simulations. A bonded-sphere model is developed to simultaneously capture irregular particle shapes and particle-wise deformations (e.g., compression, deflection, and distortion). Quantitative analysis of the effects of irregular shapes and particle deformations shows that the increase in particle stiffness tends to increase initial packing porosity and decrease the flow rate in the hopper. Rigid particles tend to have clogging issues, whereas deformable particles have less chance to, indicating particle deformation reduces the critical bridging width in the hopper flow. Detailed analysis of stress fields is also conducted to provide insights into the mechanism of particle flow and clogging. Stresses and discharge rates calculated from numerical simulations are compared and show good agreement with Walker’s theory and the extended Beverloo formula. Simulations with various particle shape combinations are also performed and show that the initial packing porosity decreases with an increasing percentage of fibers while the discharge rate has a complex dependency on particle shapes.

09 BIOMASS FUELS↗
DOE OSTI2024

Implementation of the Glued Sphere Discrete Element Model for Non-Spherical Particles in MFiX Software

To enhance solver capabilities, simulation flexibility and model validation within the MFiX software, the U.S. Department of Energy (DOE) is funding efforts to develop and integrate the glued-sphere discrete element method into the latest version of MFiX as a dedicated computational module. The glued-sphere discrete element method is a numerical technique to depict the behavior of non-spherical particles in granular flows or particulate systems by representing them as a collection of component spheres. These spheres are bonded together to approximate the shape and mechanical/chemical properties of a more complex particle. The method effectively reuses the existing sphere-sphere collision algorithm, interphase momentum and heat transfer calculations utilized in the traditional discrete element method, extending these capabilities to non-spherical particles. Additionally, this method explicitly resolves intra-particle temperature and species distributions. The MFiX glued-sphere computational module includes tools for generating glued sphere configurations, a dedicated solver, and visualization capabilities in post-processing. More specifically within the computational module, collision detection and calculations were first performed on component spheres and then mapped onto non-spherical particles. The linear spring-dashpot model was utilized to simulate the sphere-sphere interactions.

Ke, Renjie↗
DOE OSTI2025

Predicting the evolution of biomass bulk density through feedstock preprocessing: Discrete element modeling, regression analysis, and pilot-scale validation

Bulk density is an important material property of biomass feedstocks, influencing handling, storage, transport costs, and conversion efficiency. In this study, predictive regression models for loose and tapped bulk densities of Alamo and Cave-in-Rock switchgrass are developed using a comprehensive dataset generated via calibrated bonded-sphere discrete element method (DEM) simulations. Here, a key contribution of this study is the use of a DEM-based approach, which correlates density with moisture content and particle size distribution parameters and enables analysis across a continuous particle size range, overcoming limitations of purely experimental data. For comparison, regression models are also developed using only experimental data from pilot-scale runs at the Biomass Feedstock National User Facility at Idaho National Laboratory. Validation against pilot-scale data showed reasonable prediction accuracy for both model types, particularly for smaller particle sizes (post-secondary grinding). While the experimental model showed slightly better performance matching the validation data in some cases, the DEM-based model benefits from a much larger dataset, reduced predictor multicollinearity, and continuous parameter coverage, highlighting the utility of validated simulation models for developing robust predictive tools for biomass preprocessing applications.

09 - BIOMASS FUELS↗
DOE OSTI2022

Detailed biomass fast pyrolysis kinetics integrated to computational fluid dynamic (CFD) and discrete element modeling framework: Predicting product yields at the bench-scale

Fast pyrolysis is an intricate process due to the variability and anisotropy of lignocellulosic biomass and the complicated chemistry and physics during conversion in a bubbling fluidized bed reactor (BFBR). The complexity of biomass fast pyrolysis lends itself well to computational fluid dynamics (CFD) and discrete element (DEM) analysis, which promises to reduce experimental time and its associated cost. This work investigated switchgrass fast pyrolysis simulated by computational fluid dynamics coupled with a discrete element method to track individual reacting biomass particles throughout a bench-scale BFBR reactor. We accounted for the fast pyrolysis chemistry through a comprehensive reaction scheme with secondary cracking reactions. We performed a three-step reduction for secondary cracking reactions to convert the full cracking scheme into a reduced scheme easily incorporated into our model. We assessed the impact of operational conditions on the steady-state yields of liquid bio-oil, non-condensable gases (NCG), at 550 °C over a range of fluidization numbers (2 – 6 Umf), reported as a ratio to the minimum fluidization velocity (Umf). At steady-state, the volatile bio-oil yield had a range of 49.3–50.4 wt%. Levoglucosan was the primary volatile component present with 21 wt% of the bio-oil while water was the second largest with 20 wt%. The reduction of the secondary reaction schemes did not appreciably affect the overall yields of switchgrass pyrolysis compared to the full secondary scheme.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

The MFiX Particle-in-Cell Method (MFiX-PIC) Theory Guide

MFiX (Multiphase Flow with Interphase eXchanges) is an open-source multiphase flow solver developed at the National Energy Technology Laboratory. Within the code, users have access to a single phase or interpenetrating continua-based multiphase two-fluid model (TFM), a discrete element model (DEM), and a particle-in-cell model (PIC). TFM, DEM, and PIC can all be used to create multiphase simulations that include hydrodynamics, chemical reactions, and heat transfer.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2023

Coupled Experimental and High-Temperature Discrete-Element Method Modeling Studies of Aluminosilicate Particle Handling in Concentrated Solar Power Environments

Chemically inert, aluminosilicate based particles have been investigated as both a thermal transport and sensible energy storage medium for concentrated solar power facilities. These particles will experience a wide range of operating temperatures (300-1000 K) and handling conditions (dense to dilute falling particle curtains, dense granular flows, or dense structures), requiring specially-designed and optimized infrastructures. The relative influence of collisional and frictional interactions between particles varies based on temperature-dependent particulate properties and greatly impacts the bulk, granular flow behavior. These underlying physics are captured using discrete element method modeling tools. However, this modeling method is computationally expensive as each particle position and interaction is tracked during the simulation. These modeling methods are further complicated by introducing temperature-dependent particle properties, high-temperature radiative exchange, and directional irradiation sources experienced by granular flows in concentrated solar power environments. Coupled experimental and numerical studies of aluminosilicate particles in rotary kilns and dense particle curtains were performed for bulk temperatures up to 1073 K. The three particle types investigated included Carbobead HSP 30 /60, Carbobead CP 30/60, and Granusil 4030. Temperature, spatial, and velocity profile data were extracted from experimental runs using embedded K-type thermocouple probes and particle image velocimetry techniques. Experimental and numerical studies were compared using spatial temperature profiles, velocity fields, and shape profiles of the bulk, granular flows. Numerical models were developed using commercially available discrete element method modeling software, Aspherix®. Existing Aspherix® functionality was expanded by introducing coupled radiative exchange modeling tools. The laboratory-scale rotary kiln was developed to investigate the steady-state heat and mass transfer performance of aluminosilicate particles based on particle type, bulk handling temperature, and wall roughness. The rotational speed of the rotary kiln was varied to control the relative impact of collisional and frictional effects upon the granular flow behavior. Heat and mass transfer performance was categorized based on the Froude number and the observed flow regimes of slipping, rolling, cascading, and centrifuging. Coupled discrete element method modeling studies were used to evaluate the effects of temperature-dependent, particulate mechanical properties upon bulk flow behavior and upon the relative effects of radiative, advective, and/or conductive heat transfer. A high-temperature (< 1073 K) falling particle curtain was similarly fabricated to investigate the heat and mass transfer performance of aluminosilicate particles in particle handling situations dominated by inter-particle collisions. The impact of particle type, flow preheat temperatures (< 1073K), and bulk mass flow rates were investigated upon the particle curtain shape, temperature, and velocity profiles. Coupled discrete element method modeling studies were performed to evaluate the varying impact of temperature-dependent, particulate mechanical properties on the bulk flow behavior and the temperature profile of the particle curtain.

14 SOLAR ENERGY↗
DOE OSTI2021

Effects of finite element discretization and model simplification on calculations of ductile failure initiation

The finite element method is a scheme to discretize the infinite number of degrees of freedom in continuum-level problems down to a finite number of degrees of freedom. This discretization is done in conjunction with methods that also reduce the field differential equations to sets of algebraic ones that can be solved by arithmetical operations. Therefore, solutions attained by finite element models are approximations to the exact solutions of the field equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Deep Learning Predicts Stress–Strain Relations of Granular Materials Based on Triaxial Testing Data

This study presents an AI-based constitutive modelling framework wherein the prediction model directly learns from triaxial testing data by combining discrete element modelling (DEM) and deep learning. A constitutive learning strategy is proposed based on the generally accepted frame-indifference assumption in constructing material constitutive models. The low-dimensional principal stress-strain sequence pairs, measured from discrete element modelling of triaxial testing, are used to train recurrent neural networks, and then the predicted principal stress sequence is augmented to other high-dimensional or general stress tensor via coordinate transformation. Through detailed hyperparameter investigations, it is found that long short-term memory (LSTM) and gated recurrent unit (GRU) networks have similar prediction performance in constitutive modelling problems, and both satisfactorily predict the stress responses of granular materials subjected to a given unseen strain path. Furthermore, the unique merits and ongoing challenges of data-driven constitutive models for granular materials are discussed.

42 ENGINEERING↗
DOE OSTI2022

Precursory Off-Fault Deformation in Restraining and Releasing Step Overs: Insights From Discrete Element Method Models

Accelerating geophysical activity is detected preceding some, but not all, large earthquakes. This observation may indicate that no precursors occur before some earthquakes, or that the instrumentation lacks the required sensitivity. Here, to aid crustal monitoring efforts, we use discrete element method models to identify the locations and styles of deformation that may provide useful information about approaching fault reactivation. We model the reactivation of two healed rough faults in a variety of step over configurations, embedded in a host rock with varying amounts of damage subject to shear velocity loading parallel to the faults. Both the fault geometry and ratio of fault to host rock strength control the amount of off-fault deformation. Consistent with field observations, models with larger steps and more preexisting host rock damage produce higher amounts of off-fault deformation. We assess the size of the continuous regions of high velocities and strains to compare the value of the precursory information of each velocity and strain component. Comparing the three components of the velocity vector suggests that the fault-parallel velocity produces the largest and most temporally continuous regions of elevated velocity. The size of these regions increases toward failure, indicating the usefulness of tracking this component. Comparing the volumetric and shear components of the three-dimensional strain tensor suggests that during most of the interseismic period, the shear strain provides more information about approaching fault slip than the volumetric strain. However, in the days and months preceding fault reactivation, both the shear and volumetric strains provide similarly valuable information.

58 GEOSCIENCES↗
DOE OSTI2026

LANL Institutional Computing Report: Towards a digital twin of Arctic sea ice

Our goal is to develop a digital twin of Arctic sea ice that combines high-resolution predictive modeling with observational products. This effort will provide an optimized model for use in seasonal to sub-seasonal forecasting and a tool for policymakers to better anticipate, plan for, and mitigate the national security impacts of rapidly changing Arctic conditions. The modeling component uses the Discrete Element Model for Sea Ice (DEMSI), which uses discrete elements to represent the sea ice with an explicit representation of forces between them.

58 GEOSCIENCES↗
DOE OSTI2020

Theoretical Review of the MFIX Fluid and Two-Fluid Models

MFIX (Multiphase Flow with Interphase eXchanges) is a general purpose code that can be used for describing the hydrodynamics, chemical reactions, and heat transfer of dense or dilute multiphase fluid flows. MFIX calculations give detailed information on pressure, temperature, composition, and velocity distributions in the defined system. This report provides an overview of MFIX’s expressions for single phase flow and continuous dispersed multiphase flow (i.e., the two-fluid model). In addition to the two fluid model (MFIX-TFM), MFIX contains a discrete element model (MFIX-DEM) and a multi-phase particle in cell model (MFIX-PIC) that are not reviewed within this document.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Protocol-dependent frictional granular jamming simulations: cyclical, compression, and expansion

Granular matter takes many paths to pack in natural and industrial processes. The path influences the packing microstructure, particularly for frictional grains. We perform discrete element modeling simulations of different paths to construct packings of frictional spheres. Specifically, we explore four stress-controlled protocols implementing packing expansions and compressions in various combinations thereof. We characterize the eventual packed states through their dependence of the packing fraction and coordination number on packing pressure, identifying non-monotonicities with pressure that correlate with the fraction of frictional contacts. These stress-controlled, bulk-like particle simulations access very low-pressure packings, namely, the marginally stable limit, and demonstrate the strong protocol dependence of frictional granular matter.

36 MATERIALS SCIENCE↗
DOE OSTI2021

DEM analysis of the effect of interface transition zone on dynamic splitting tensile behavior of high-strength concrete based on multi-phase model

Highlights: • Crushable aggregates with the real geometry shape are generated by using clump-cluster method. • Flat-Joint Model is adopted as the constitutive law. • The tensile strength is almost not affected by ITZ strength when its strength exceeds 0.8 times that of mortar. • A modified formula that can reflect the functional relationship between the ITZ strength and the TDIF is proposed. The tensile strength is one of the main parameters in the design of concrete structures and has an important influence on the cracking performance and durability of concrete materials. As the weak region in concrete, the strength of the interface transition zone (ITZ) directly affects the dynamic mechanical response of concrete. Meanwhile, considering that the tensile performance of concrete is far weaker than its compressive performance, the study of the strength of ITZ is of great significance to the dynamic tensile fracture characteristics and mechanical properties of concrete. In this paper, a method combining 3D scanning technology and Clump-Cluster was proposed to generate a real geometry shape model of crushable aggregate. On this basis, the 3D discrete element model (DEM) of high-strength concrete was constructed using the Flat-Joint Model (FJM). The effects of the ITZ strength on the failure mode, number of microcracks, tensile strength and strain rate effect of high-strength concrete under strain rates from 0.01/s to 100/s were studied by splitting tensile test. The simulation results indicate that the number of microcracks generated decreases with the increase of ITZ strength, and the proportion of ITZ cracks decreases obviously. When the ITZ strength exceeds the critical value of 0.8 times that of mortar, the tensile strength is almost not affected by ITZ strength. In addition, based on the numerical simulation results, a modified formula that can reflect the functional relationship between the ITZ strength and the dynamic increasing factor of the concrete tensile strength (TDIF) is proposed, and the rationality of the TDIF modified formula is verified.

36 MATERIALS SCIENCE↗