Search NASA⌕ Search

SEARCH · Search NASA

Results for “discrete element”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

GPU accelerated Discrete Element Method (DEM) molecular dynamics for conservative, faceted particle simulations

Faceted shapes, such as polyhedra, are commonly found in systems of nanoscale, colloidal, and granular particles. Many interesting physical phenomena, like crystal nucleation and growth, vacancy motion, and glassy dynamics are challenging to model in these systems because they require detailed dynamical information at the individual particle level. Within the granular materials community the Discrete Element Method has been used extensively to model systems of anisotropic particles under gravity, with friction. We provide an implementation of this method intended for simulation of hard, faceted nanoparticles, with a conservative Weeks–Chandler–Andersen (WCA) interparticle potential, coupled to a thermodynamic ensemble. This method is a natural extension of classical molecular dynamics and enables rigorous thermodynamic calculations for faceted particles.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

A discrete element method-based approach to predict the breakage of coal

Pulverization is an essential pre-combustion technique employed for solid fuels, such as coal, to reduce particle sizes. Smaller particles ensure rapid and complete combustion, leading to low carbon emissions. Traditionally, the resulting particle size distributions from pulverizers have been informed by empirical or semi-empirical approaches that rely on extensive data gathered over several decades during operations or experiments. However, the predictive capabilities for new coals and processes are limited. This work presents a Discrete Element Method based computational framework to predict particle size distribution resulting from the breakage of coal particles characterized by the coal’s physical properties. The effect of certain operating parameters on the breakage behavior of coal particles also is examined.

Brazilian Test.↗
DOE OSTI2018

Discrete element modeling of grain crushing and implications on reservoir compaction

Depletion-induced formation compaction is a concern for production in hydrocarbon reservoirs. In this study, numerical mechanical tests and digital samples are proposed as an alternative to uniaxial strain experiments for investigating the effects of compaction. We propose a grain crushing model with two-sized subgrains based on the discrete element method and bonded-particle model (DEM-BPM) and perform numerical single grain crushing tests to calibrate the model parameters against the experimental sand grain characteristic strength. Numerical uniaxial strain tests are used to predict porosity and permeability reduction caused by reservoir compaction. Grain crushing incurs a significant reduction of porosity under the uniaxial strain stress path, which becomes more evident for digital samples consisting of larger grains. The stress-dependent changes in porosity and permeability obtained from the DEM simulations can be used as inputs for a reservoir-scale correlation to predict the productivity loss due to reservoir compaction.

nuclear (including radiation effects), carbon sequ↗
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 OSTI2020

Three-Dimensional Evaluation of Sand Particle Fracture Using Discrete-Element Method and Synchrotron Microcomputed Tomography Images

Recent research showed that fracture of sand particles plays a significant role in determining the plastic bulk volumetric changes of granular materials under different loading conditions. One of the major tools used to better understand the influence of particle fracture on the behavior of granular materials is discrete-element modeling (DEM). This paper employed the bonded block model (BBM) to simulate the fracture behavior of sand. Each sand particle is modeled as an agglomerate of rigid blocks bonded at their contacts using the linear-parallel contact model, which can transmit both moment and force. DEM simulated particles closely matched the actual three-dimensional (3D) shape of sand particles acquired using high-resolution 3D synchrotron microcomputed tomography (SMT). Results from unconfined one-dimensional (1D) compression of a single synthetic silica cube were used to calibrate the model parameters. Particle fracture was investigated for specimens composed of three sand particles that were loaded under confined 1D compression. Breakage energy measured from DEM models matched well with that measured experimentally. The paper studied the effects of contact loading condition and particle interaction on the fracture mode of particles using BBM that can closely capture the 3D shape of real sand particles.

58 GEOSCIENCES↗
DOE OSTI2021

Assessment of a tomography-informed polyhedral discrete element modelling approach for complex-shaped granular woody biomass in stress consolidation

The design of handling equipment for granular biomass primarily requires experiments. Discrete element models (DEM) can provide designers with detailed insight into the behaviour of granular materials. However, granular biomass comprising complex-shaped particles is difficult to model with DEM. A tomography-informed DEM approach and an exhaustive assessment of the ability of the approach to predict the bulk behaviour of milled pines using experimental data is presented. Nano-CT scan was conducted to obtain 3D particle surface geometries as the basis for particle shape approximation by a polyhedral model and a sphero-polyhedral model. These models were applied in the simulation of a compressibility test. Our parametric study showed that particle Young's modulus and restitution coefficient are the two main properties influencing the simulated bulk behaviour of the DEM particles and that the level of approximating DEM particle shape for real particles is critical for accurately replicating the bulk behaviour of milled pines. The polyhedral model demonstrated better suitability than its sphero-polyhedral counterpart for modelling pine particles. The polyhedral model calibrated in a compressibility test was then applied in the simulations of a friction test without additional parameter tuning. Interlocking was dominant in the shear of bulk pine particles. Remarkably, the polyhedral model predicted the frictional behaviour of the pine particles when compared to the experimental results. The limitations of the model, as well as possible ways for enhancement, are discussed. This work provided novel insights into the suitability of complex-shaped DEM models for granular woody biomass.

Biomass flow↗
DOE OSTI2023

Discrete element model for powder grain interactions under high compressive stress

A reduced order, nonlocal model is proposed for the contact force between initially spherical particles under compression. The model in effect provides the normal component of the interaction force between elements in the discrete element method (DEM). It is applicable to high relative density and large stress in powder compaction. It takes into account the mutual interaction between multiple points of contact, in contrast to the usual assumption in DEM of pair interactions. The mathematical form of the model is derived from a variational formulation that leads to the momentum balance for the forces on each grain. The model is calibrated mainly using detailed three dimensional peridynamic simulations of single grains under compressive loading by rigid plates that move radially with prescribed velocity. This calibration takes into account the large deformation and fracture of the grains. The interaction model also includes terms for the unloading behavior and adhesion. Finally, as validation, the model is applied to test data on the compaction of microcrystalline cellulose bulk powder.

36 MATERIALS SCIENCE↗
NASA NTRS1972

A design study for the addition of higher order parametric discrete elements to NASTRAN

The addition of discrete elements to NASTRAN poses significant interface problems with the level 15.1 assembly modules and geometry modules. Potential problems in designing new modules for higher-order parametric discrete elements are reviewed in both areas. An assembly procedure is suggested that separates grid point degrees of freedom on the basis of admissibility. New geometric input data are described that facilitate the definition of surfaces in parametric space.

Stanton, E. L.↗
DOE OSTI2020

A smooth contact algorithm for the combined finite discrete element method

From its inception, the combined finite discrete element method has used a distributed potential contact force algorithm to resolve interaction between finite elements. The contact interaction algorithm relies on evaluation of the contact force potential field. The problem with existing algorithms is that the potential field introduces artificial numerical non-smoothness in the contact force. This work introduces a smooth potential field based on the finite element topology, and a generalized contact interaction law is constructed on top of the smooth potential field. Further, a number of validation cases for the proposed algorithm, considering different shapes of discrete elements, are presented, and detailed aspects of the proposed contact interaction law are tested with numerical examples.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Simulation of crack induced nonlinear elasticity using the combined finite-discrete element method

Numerical simulation of nonlinear elastic wave propagation in solids with cracks is indispensable for decoding the complicated mechanisms associated with the nonlinear ultrasonic techniques in NonDestructive Testing (NDT). Here, we introduce a two-dimensional implementation of the combined finitediscrete element method (FDEM), which merges the finite element method (FEM) and the discrete element method (DEM), to explicitly simulate the crack induced nonlinear elasticity in solids with both horizontal and inclined cracks. In the FDEM model, the solid is discretized into finite elements to capture the wave propagation in the bulk material, and the finite elements along the two sides of the crack also behave as discrete elements to track the normal and tangential interactions between crack surfaces. The simulation results show that for cracked models, nonlinear elasticity is generated only when the excitation amplitude is large enough to trigger the contact between crack surfaces, and the nonlinear behavior is very sensitive to the crack surface contact. The simulations reveal the influence of normal and tangential contact on the nonlinear elasticity generation. Moreover, the results demonstrate the capabilities of FDEM for decoding the causality of nonlinear elasticity in cracked solid and its potential to assist in Non-Destructive Testing (NDT).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

A non-locking composite tetrahedron element for the combined finite discrete element method

Here, a number of validation cases considering different loading and boundary conditions and different materials for the proposed element are presented. A practical application of the use of the composite tetrahedral finite element is presented by quantitative comparisons of numerical model results against simple theoretical estimates and results from acrylic fracturing experiments. All of these examples clearly show the capability of the composite element in eliminating volumetric locking.

42 ENGINEERING↗
DOE OSTI2022

Geometric remapping of particle distributions in the Discrete Element Model for Sea Ice (DEMSI v0.0)

Abstract. A new sea ice dynamical core, the Discrete Element Model for Sea Ice (DEMSI), is under development for use in coupled Earth system models. DEMSI is based on the discrete element method, which models collections of ice floes as interacting Lagrangian particles. In basin-scale sea ice simulations the Lagrangian motion results in significant convergence and ridging, which requires periodic remapping of sea ice variables from a deformed particle configuration back to an undeformed initial distribution. At the resolution required for Earth system models we cannot resolve individual sea ice floes, so we adopt the sub-grid-scale thickness distribution used in continuum sea ice models. This choice leads to a series of hierarchical tracers depending on ice fractional area or concentration that must be remapped consistently. The circular discrete elements employed in DEMSI help improve the computational efficiency at the cost of increased complexity in the effective element area definitions for sea ice cover that are required for the accurate enforcement of conservation. An additional challenge is the accurate remapping of element values along the ice edge, the location of which varies due to the Lagrangian motion of the particles. In this paper we describe a particle-to-particle remapping approach based on well-established geometric remapping ideas that enforces conservation, bounds preservation, and compatibility between associated tracer quantities, while also robustly managing remapping at the ice edge. One element of the remapping algorithm is a novel optimization-based flux correction that enforces concentration bounds in the case of nonuniform motion. We demonstrate the accuracy and utility of the algorithm in a series of numerical test cases.

58 GEOSCIENCES↗
DOE OSTI2018

A discrete element method representation of an anisotropic elastic continuum

A method for modeling cubically anisotropic elasticity within the discrete element method is presented. The discrete element method (DEM) is an approach originally intended for modeling granular materials (sand, soil, and powders); however, recent developments have usefully extended it to model stochastic mechanical processes in monolithic solids which, to date, have been assumed to be elastically isotropic. The method presented here for efficiently capturing cubic elasticity in DEM is an important prerequisite for further extending DEM to capture the influence of elastic anisotropy on the mechanical response of polycrystals, composites, etc. The system demonstrated here uses a directionally assigned stiffness in the bonds between adjacent elements and includes separate schemes for achieving anisotropy with Zener ratios greater and smaller than one. Finally, the model framework is presented along with an analysis of the accessible space of elastic properties that can be modeled and an artificial neural network interpolation scheme for mapping input parameters to model elastic behavior.

42 ENGINEERING↗
NASA NTRS1971

An Evaluation of Current NASTRAN Discrete Element Models for Monocoque and Semimonocoque Structures

The current discrete element analysis models in NASTRAN, Levels 8. 1.0 and 11. 1. 0, are evaluated for semi-rnonocoque and monocoque structures. Comparative results are presented with the FORMAT-I11 general-purpose code for a swept wing statics analysis and a Skylab substructure vibration analysis. Good agreement between Level i I. 1.0 and FORMAT-111 is demonstrated; however, a major error in Level 8. 1.0 was discovered on the Skylab vibration problem. The computational performance of NASTRAN was superior to FORMAT, but the vibration rigid format was found to have several inefficiencies. Comparative results are also presented between Level 11. 1. 0 and STACUSS, SLADE, BOND, UNIVALVE, and another special-purpose code. Results for two of the structures analyzed, a cylindrical shell and an elliptic plate, indicate curved discrete elements are significantly more accurate than the plate elements in NASTRAN for the same number of degrees of freedom.

Edward L. Stanton↗