Search NASA⌕ Search

SEARCH · Search NASA

Results for “Discrete element method (DEM)”

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.

83 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

CFD-DEM Modeling of Autothermal Pyrolysis of Corn Stover with a Coupled Particle- and Reactor-Scale Framework

Autothermal operation of fast pyrolysis is an efficient process-intensification technique wherein exothermic oxidation reactions are used to overcome the heat-transfer bottleneck of conventional pyrolysis. The development of accurate, reliable modeling toolsets is imperative to generating a deeper understanding of biomass autothermal pyrolysis systems to support scale-up and industrial deployment. This modeling effort describes the development of single-particle and reactor models which incorporate detailed reaction schemes and simultaneous exothermic oxidation reactions. The particle-scale model was parameterized for corn stover feedstock with particle morphology, density, ash content, and biopolymer composition, all of which impact the emergent conversion characteristics during pyrolysis. Results were then used to parameterize a reactor-scale autothermal pyrolysis model, which was developed using a coarse-grained computational fluid dynamic-discrete element method. The simulation results compared well with experimental results, with the predicted bio-oil, light gas, and biochar yield within 3.0 wt% of the experimental yields. Further analyses were performed to test the influence of equivalence ratio, biomass injection position, and particle size distribution on autothermal pyrolysis. The analysis of the physio-chemical properties of the fluid and solid phase inside the reactor and at the reactor outlet help reveal important process interactions of autothermal pyrolysis.

BIOMASS FUELS,INORGANIC, ORGANIC, PHYSICAL, AND AN↗
DOE OSTI2024

Two-Fluid and Discrete Element Modeling of a Parallel Plate Fluidized Bed Heat Exchanger for Concentrating Solar Power

A novel high-temperature particle solar receiver is developed using a light trapping planar cavity configuration. As particles fall through the cavity, the concentrated solar radiation warms the boundaries of the receiver and in turn heats the particles. Particles flow through the system, forming a fluidized bed at the lower section, leaving the system from the bottom at a constant flowrate. Air is introduced to the system as the fluidizing medium to improve particle heat transfer and mixing. A laboratory scale cavity receiver is built by collaborators at the Colorado School of Mines and their data are used for model validation. In this experimental setup, near IR quartz lamp is used to provide flux to the vertical wall of the heat exchanger. The system is modeled using the discrete element method and a continuum two-fluid method. The computational model matches the experimental system size and the particle size distribution is assumed monodisperse. A new continuum conduction model that accounts for the effects of solid concentration is implemented, and the heat flux boundary condition matches the experimental setup. Radiative heat transfer is estimated using a widely used correlation during the post-processing step to determine an overall heat transfer coefficient. The model is validated against testing data and achieves less than 30% discrepancy and a heat transfer coefficient greater than 1000 W/m2 K.

concentrating solar power↗
DOE OSTI2025

A rheological model for loose sands with insights from DEM

A rheological model for loose granular media is developed to capture both solid-like and fluid-like responses during shearing. The proposed model is built by following the mathematical structure of an extended Kelvin–Voigt model, where an elastic spring and plastic slider act in parallel to a viscous damper. This arrangement requires the partition of the total stress into rate-independent and rate-dependent stress components. To model the solid-like behavior, a simple frictional plasticity model is adopted without modifications, thus contributing to the rate-independent stress. Instead, the fluid-like or rate-dependent stress is further decomposed into deviatoric and volumetric parts, by proposing a new formulation based on a combination of the μ(I) relation, originally developed under pressure-controlled shear, with a pressure-shear rate relation derived under volume-controlled shear. The proposed formulation allows the model to capture both the increase in the friction coefficient and the enhanced dilation at high shear rates. High-fidelity simulation data, obtained from discrete element method and multiscale modelling, are used to evaluate the performance of the proposed constitutive model. The model provides accurate results under both drained and undrained simple shear paths across a wide range of shear rates. Furthermore, it successfully reproduces at much lower computational cost the flowslide mobility computed through multiscale simulations, which is primarily regulated by the shear rate dependence of the material properties during the dynamic runout stage.

Elasticity↗
DOE OSTI2025

Dense granular flows with MFIX-Exa

This report extends the linear spring dashpot collision model the discrete element method available in MFIX-Exa to include static a static tangential friction force. Additionally, two rolling friction models frequently used in the literature are also implemented. The governing equations are provided with an emphasis on the new terms. The new model is validated by comparison to existing experimental data of single particle oblique collisions. The model is then tested on three dense granular flow problems: the formation of static piles, the discharge from a flat-bottom hopper and the self-induced granular Rayleigh-Taylor instability.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

MFIX-Exa: A path toward exascale CFD-DEM simulations

MFIX-Exa is a computational fluid dynamics–discrete element model (CFD-DEM) code designed to run efficiently on current and next-generation supercomputing architectures. MFIX-Exa combines the CFD-DEM expertise embodied in the MFIX code—which was developed at NETL and is used widely in academia and industry—with the modern software framework, AMReX, developed at LBNL. The fundamental physics models follow those of the original MFIX, but the combination of new algorithmic approaches and a new software infrastructure will enable MFIX-Exa to leverage future exascale machines to optimize the modeling and design of multiphase chemical reactors.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Critical analysis of velocimetry methods for particulate flows from synthetic data

Particle tracking methods that extract high-fidelity particle velocity data from high speed video of particle laden flows is a common experimental technique applied to chemical processes. These measurements are used to better understand the motion of particles and fluids in complex systems and create data against which computational models are validated. However, the methods, codes, and experimental setups all have limitations. It is imperative that practitioners verify the methods and their implementation as well as understand the limitations of experimental setups. This work focuses on quantifying the visible depth of field in a high particle concentration fluidized bed. Following a precedent set by the particle imaging velocimetry community, a particle velocity field is manufactured using a computational fluid dynamics and discrete element method simulation. Photo realistic high-speed videos are rendered based on the simulated data using the three-dimensional creation software Blender. Particle velocities are extracted from the synthetic high-speed videos using three variants of Particle Tracking Velocimetry and Optical Flow Velocimetry methodologies. Here, the tracked results are then compared to the known solution, quantifying the error associated with the assumed visible depth. The results indicate that at depth of one particle diameter, all three particle tracking codes give accurate measurements, largely within 5%. However, the error increases when the full bed video measurements are compared to the known solution at one particle diameter, i.e., mimicking a validation study. Finally, for some statistics the constant depth assumption only increases the error slightly, for others significantly.

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

Exploring the Regimes of Particle Behavior upon Impact via the Discrete Element Method

Discrete element method simulations are conducted to probe the various regimes of post-impact behavior of particles with solid surfaces. The impacting particles are described as spherical agglomerates consisting of smaller constituent (or primary) particles held together via surface adhesion. Under the influence of a wide range of impact velocities and particle surface energies, five distinct behavioral regimes—rebounding, vibration, fragmentation, pancaking, and shattering—are identified, and force transmission patterns are linked to post-impact behavior. In the rebounding regime, the coefficient of restitution decreases linearly as impact velocity increases and the particle agglomerate experiences compaction. In the fragmentation regime, rebound velocity generally decreases with increasing fragment size. The rebound velocity of fragments decreases with time except for the smallest fragments, which can increase in velocity due to collisions with other fragments of high velocity. Particle breakage in the pancaking regime does not follow common mechanistic models of breakage.

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

A nonlinear elasto-plastic bond model for the discrete element modeling of woody biomass particles

This work presents an experiment-informed, semi-empirical, elasto-plastic bond model for the discrete element modeling of woody biomass particles. The model renders nonlinear plastic deformation of materials when subjected to compression/tension, bending, and twisting, essential for accurately simulating the behavior of biomass in comminution. The model is implemented in an open-source DEM package LIGGGHTS and is assessed in a number of verification tests. The model is applied to simulate the fracturing test of notched loblolly pine blocks. Reduced bond strength is prescribed to the weak planes along the growth rings. It is found that the square-root formulation for describing the elasto-plastic normal bond contact delivers the best agreement with the experimental data over its linear and quadratic counterparts. Furthermore, the model predicts a strong anisotropic behavior of wood blocks with regard to the grain orientations, indicating that the cutting angle should be an important attribute for the optimization of biomass comminution.

42 ENGINEERING↗
DOE OSTI2022

Memory access optimization for particle operations in computational fluid dynamics-discrete element method simulations

Computational Fluid Dynamics - Discrete Element Method is used to model gas-solid systems in several applications in energy, pharmaceutical and petrochemical industries. Computational performance bottlenecks often limit the problem sizes that can be simulated at industrial scale. The data structures used to store several millions of particles in such large-scale simulations have a large memory footprint that does not fit into the processor cache hierarchies on current high-performance-computing platforms, leading to reduced computational performance. This paper specifically addresses this aspect of memory access bottlenecks in industrial scale simulations. The use of space-filling curves to improve memory access patterns is described and their impact on computational performance is quantified in both shared and distributed memory parallelization paradigms. The Morton space filling curve applied to uniform grids and k-dimensional tree partitions are used to reorder the particle data-structure thus improving spatial and temporal locality in memory. The performance impact of these techniques when applied to two benchmark problems, namely the homogeneous-cooling-system and a fluidized-bed, are presented. We report these optimization techniques lead to approximately two-fold performance improvement in particle focused operations such as neighbor-list creation and data-exchange, with ~ 1.5 times overall improvement in a fluidization simulation with 1.27 million particles.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Investigation of Best-Practices and Computationally Inexpensive Radiative Exchange Models for Discrete Element Method Modeling of Aluminosilicate Particles in Concentrating Solar Power Environments

Chemically inert, aluminosilicate based particles have been investigated as both a thermal transport and sensible energy storage medium for concentrating 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 concentrating solar power environments. In this study, coupled experimental and numerical slump testing of aluminosilicate particles was performed and computationally efficient radiative exchange models were evaluated to establish best-practices for discrete element method models for concentrating solar power environments. The three particle types investigated included Carbobead HSP 30 /60, Carbobead CP 30/60, and Granusil 4030. Existing modeling limitations and computationally-efficient multi-modal heat transfer models were evaluated using Aspherix®, a commercial discrete element method software. High-temperature (< 1073 K) slump testing of aluminosilicate particles was performed to investigate the deviation between experimentally-observed and numerically-predicted angles of repose introduced by computation-time reduction practices including the relaxation of the particle elastic modulus and coarse-graining. Coarse-graining is used to use a single modeled particle that is representative of a collection of smaller particles, decreasing the computational cost at the expense of geometric accuracy. Additionally, relaxation of the elastic modulus is used to reduce computational time at the expense of an increased, modeled particle overlap. Prior studies have determined that aluminosilicate particles retain a high elastic modulus at high temperatures (< 1073 K), requiring small simulation timesteps to ensure resolved contact forces resemble appropriate solid mechanics. A parametric study was performed to evaluate the influence of computation time improvements on the deviation between experimental and modeled angle of repose across high temperatures < 1073 K. Additionally, numerical case studies were performed on candidate particle systems at varying porosities and temperatures. These studies were performed to investigate the influence of computationally-efficient radiative-exchange modeling methods coupled to Aspherix® on modeled accuracy and computation time. The recently-developed distance-based approximation was evaluated in estimating radiative exchange between particles and participating surfaces located in close proximity. The distance based approximation was developed to use tabulated estimates of the radiative distribution factor between individual particles and surfaces in close proximity (< 40 particle radii). These methods were expanded to the aluminosilicate particles of interest, including the influence of particle size distributions. To capture radiative exchange between particles and surfaces not in close proximity (> 40 particle radii) and to capture the absorption of directional irradiation from concentrating solar resources, a volumetrically-averaged radiative distribution factor was calculated between the modeled granular flow and surfaces using Monte Carlo ray-tracing for participating media. Volume-averaged absorption and scattering coefficients were predicted using a volumetric discretization of the modeled domain with monodisperse approximations based on geometric optics and experimentally-determined scattering phase functions for aluminosilicate particles.

14 SOLAR ENERGY↗