Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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.

At least 37 records · Page 2

CFD-DEM Onset of Motion Analysis for Application to Bed Scour Risk Assessment

This CFD study with DEM was done as a part of the Federal Highway Administration’s (FHWA’s) effort to improve scour design procedures. The Computational Fluid Dynamics-Discrete Element Method (CFD-DEM) model, available in CD-Adapco’s StarCCM+ software, was used to simulate multiphase systems, mainly those which combine fluids and solids. In this method the motion of discrete solids is accounted for by DEM, which applies Newton's laws of motion to every particle. The flow of the fluid is determined by the local averaged Navier–Stokes equations that can be solved using the traditional CFD approach. The interactions between the fluid phase and solids phase are modeled by use of Newton's third law. The inter-particle contact forces are included in the equations of motion. Soft-particle formulation is used, which allows particles to overlap. In this study DEM was used to model separate sediment grains and spherical particles laying on the bed with the aim to analyze their movement due to flow conditions. Critical shear stress causing the incipient movement of the sediment was established and compared to the available experimental data. An example of scour around a cylindrical pier is considered. Various depths of the scoured bed and flow conditions were taken into account to gain a better understanding of the erosion forces existing around bridge foundations. The decay of these forces with increasing scour depth was quantified with a ‘decay function’, which shows that particles become increasingly less likely to be set in motion by flow forces as a scour hole increases in depth. Computational and experimental examples of the scoured bed around a cylindrical pier are presented.

97 MATHEMATICS AND COMPUTING↗

BOREAS HYP-8 DEM Data Over The NSA-MSA and SSA-MSA in The AEAC Projection

These data were derived from the original Digital Elevation Models (DEMs) produced by the Boreal Ecosystem-Atmosphere Study (BOREAS) Hydrology (HYD)-8 team. The original DEMs were in the Universal Transverse Mercator (UTM) projection, while this product is projected in the Albers Equal-Area Conic (AEAC) projection. The pixel size of the data is 100 meters, which is appropriate for the 1:50,000-scale contours from which the DEMs were made. The original data were compiled from information available in the 1970s and 1980s. This data set covers the two Modeling Sub-Areas (MSAs) that are contained within the Southern Study Area (SSA) and the Northern Study Area (NSA). The data are stored in binary, image format files. The DEM data over the NSA-MSA and SSA-MSA in the AEAC projection are available from the Earth Observing System Data and Information System (EOSDIS) Oak Ridge National Laboratory (ORNL) Distributed Active Archive Center (DAAC). The data files are available on a CD-ROM (see document number 20010000884).

Knapp, David E.↗

Discrete-Element and Material-Point Method (DEM and MPM) Based Solvers for Sustainable Technologies

We present the use of discrete element method (DEM) and material point method (MPM) in three relevant green technology applications that include biomass feedstock handling, lithium-ion battery manufacturing, and high-pressure reverse osmosis. Our open-source DEM and MPM solvers are developed using performance portable grid and particle management library, AMReX, thus enabling superior performance on NVIDIA and AMD GPUs with > 100 million particles. Our DEM solver resolves the motion of individual particles in a granular system and includes a bonded sphere method for modeling non-spherical particles along with Hertzian and liquid bridge-based contact models. We simulate highly variable biomass feedstock flows in large-scale hoppers for biofuel production and electrode calendering in battery manufacturing using DEM. Our simulations predict flow blockage in large scale biomass hoppers and electrode microstructure variations, thus providing valuable information for biofuel and battery manufacturers, respectively. The second half of the talk will be on MPM and its application towards pore resolved simulations of reverse osmosis membranes under compressive loads. We present a validation study of our MPM simulations with membrane microscopy imaging thus providing useful insights on membrane stability under high pressure conditions. We also present a spectral stability analysis of using linear hat, quadratic and cubic spline basis in MPM indicating regions of numerical stability.

BIOMASS FUELS,MATHEMATICS AND COMPUTING↗

Utah FORGE: Area LiDAR Bare Earth DEM Mosaic

This is a mosaic of 154 bare earth DEM panels (0.5 m resolution) acquired during Phase II of the Utah FORGE project. The DEM covers the western Mineral Mountains and adjoining basin including the Roosevelt Hot Springs area. The coordinates are in UTM Zone 12, NAD 83 and elevation units are in meters. The DEM is in GeoTiff format.

15 GEOTHERMAL ENERGY↗

DEM, DSM, and Cleaned LiDAR Point Cloud Data from the NGEE Arctic UAS Campaigns at the Teller 27 Field Site from 2017 and 2018, Seward Peninsula, Alaska

A Digital Elevation Model (DEM) and Digital Surface Model (DSM) were derived from airborne Light Detection and Ranging (LiDAR) data collected from Los Alamos National Laboratory's (LANL) heavy-lift unoccupied aerial system (UAS) quadcopter and hexacopter platforms operated by Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic) scientists from the EES-14 group at LANL. These data were collected in August 2017 and July 2018 at the NGEE Arctic field site near mile marker 27 of the Bob Blodgett Nome-Teller Memorial Highway between Nome, Alaska and Teller, Alaska. A Vulcan Raven X8 Airframe (Mitcheldean, Gloucestershire, UK), DJI Matrice 600 Pro Airframe (Shenzhen, China), and Routescene UAV LiDARSystem (Edinburgh, Scotland, UK) were used to collect LiDAR data. Following pre-processing in Routescene LidarViewer Pro software, the LiDAR point clouds were cleaned and processed using CloudCompare software to separate ground and off-ground points. A high resolution DEM and DSM were then created using ArcGIS Pro software. This data package contains fully cleaned point clouds of ground and off-ground points (.las), a 25 cm DEM (.tif), and a 25 cm DSM (.tif) for the Teller 27 field site. Ancillary aircraft data, flight mission parameters, weather conditions, and raw lidar data and imagery can be found in the L0 datasets for these campaigns: NGA299 (2017) and NGA297 (2018). Minimally processed point clouds and auxiliary files can be found in the L1 dataset: NGA304 (2017 and 2018).The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a 15-year research effort (2012-2027) to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

A three-dimensional meso-scale approach to concrete fracture based on combined DEM with X-ray μCT images

Highlights: • We investigate fracture in concrete beams subjected to 3-point bending. • We use a discrete element method (DEM) to capture fracture. • Concrete is described as 4-phase material. • Evolution of micro-structural phenomena is investigated. • The 3D results are carefully compared with the 2D ones. - Abstract: The paper shows three-dimensional (3D) numerical mesoscopic results on progressive concrete fracture at aggregate level in a notched concrete beam under bending. A classical particle discrete element method was used to capture fracture in concrete. Heterogeneity was taken into account with a 4-phase concrete description: aggregate, cement matrix, macro-voids and interfacial transitional zones. Both the real shape and location of aggregate particles in concrete were assumed in DEM analyses based on X-ray μCT scans. A satisfactory agreement in terms of the vertical force versus crack mouth opening displacement evolution and crack geometry was achieved between DEM analyses and own laboratory tests. The evolution of broken contacts, contact normal forces, particle rotations, internal energies was studied at the aggregate level. The 3D results were also carefully compared with the 2D ones. The 3D model proved to be a promising simulation tool in predicting concrete fracture at the mesoscopic and macroscopic level.

36 MATERIALS SCIENCE↗

Effects of curvature-related DEM contact model on the macro- and micro-mechanical behaviours of granular soils

The presented paper conducted a comprehensive comparison between the Hertz-Mindlin model and the linear spring model in true triaxial shear simulations of granular soils using the discrete element method (DEM). The no-slip Hertz-Mindlin model for general elastic non-spherical particles with smooth surfaces was revisited, and implemented for superellipsoidal particles in an in-house DEM code. Three groups of specimens with a grain size distribution of Ottawa 20-30 sands, consisting of spheres, ellipsoids, and superellipsoids, respectively, were subjected to triaxial shear DEM simulations with the Hertz-Mindlin model and the linear spring model. The corresponding mechanical behaviours were examined in terms of a series of macro-and micro-parameters. It is found that the linear spring model is able to resemble the Hertz-Mindlin model on aspects of both microscopic and macroscopic mechanical behaviours of granular media with spherical and/or non-spherical particles. This finding suggests that the linear spring model can be used to investigate micro-mechanical behaviours of granular soils even with complex particle shapes.

42 ENGINEERING↗

BOREAS HYD-8 DEM Data Over the NSA-MSA and SSA-MSA in the UTM Projection

The BOREAS HYD-8 team focused on describing the scaling behavior of water and carbon flux processes at local and regional scales. These DEMs were produced from digitized contours at a cell resolution of 100 meters. Vector contours of the area were used as input to a software package that interpolates between contours to create a DEM representing the terrain surface. The vector contours had a contour interval of 25 feet. The data cover the BOREAS MSAs of the SSA and NSA and are given in a UTM map projection. Most of the elevation data from which the DEM was produced were collected in the 1970s or 1980s. The data are stored in binary, image format files. The data files are available on a CD-ROM (see document number 20010000884) or from the Oak Ridge National Laboratory (ORNL) Distributed Active Archive Center (DAAC).

Wang, Xue-Wen↗

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.↗

Aster Global Dem (gdem) Version 3

The Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) is a 14-channel imaging instrument operating on NASA’s Terra satellite since 1999. ASTER’s visible–near infrared (VNIR) instrument, with three bands and a 15 m Instantaneous field of view (IFOV), is accompanied by an additional VNIR band using a second, backward-looking telescope. Collecting along-track stereo pairs, the geometry produces a base-to-height ratio of 0.6. In 2009, the ASTER Science Team released Version 1 of the global DEM (GDEM) based on stereo correlation of 1.2 million ASTER scenes. The DEM has 1 arc-second latitude and longitude postings (~30 m) and employed cloud masking to avoid cloud-contaminated pixels. The GDEM covers all of the Earth’s land surface from 83 degrees north to 83 degrees south latitude. Version 2 was released in 2011, with notable improvements in coverage and accuracy. In 2019, the final, Version 3, was released; again improving on coverage and removing almost all artifacts. Th GDEM is a unique, global high spatial resolution digital elevation dataset available to all users at no cost. In addition, a second unique dataset was produced and released. The raster-based ASTER Global Water Body Dataset (ASTWBD) identifies the presence of permanent water bodies, and marks them as ocean, lake, or river. An accompanying DEM file indicates the elevation for each water pixel. To date, over 110+ million 1x1 degree GDEM tiles have been distributed.

Crippen, R.↗

An Orthogonal Recursive Bisection (ORB) Based Time Advancement Algorithm for CFD-DEM Solvers

The time integration of the granular phase in coupled computational fluid dynamics (CFD) – discrete element method (DEM) simulations presents a unique computational challenge brought about by the large variations in particle collisional time scales. Particles in the dilute regions of the computational domain can be advanced with large time steps while dense regions require much smaller time increments. However, the time step size in most solvers is globally set as the limit for accuracy and stability imposed by the collisions and is typically orders of magnitude less than that required away from collisions. This work addresses this precise issue and provides a strategy to avoid the use of a global conservative small time step size for the entire set of particles.A novel time stepping algorithm for CFD-DEM solvers using a partitioning approach using orthogonal recursive bisection (ORB) that allows for variable time steps among particles is described and its computational performance is compared against baseline explicit methods, typically used in several CFD-DEM solvers. ORB has advantages of being relatively quick and easy to update incrementally and has the required heuristic behavior (i.e., it will split the region in half with a cluster on each side) when groups of particles are well separated (clustered). The algorithm presented in this work uses a local time stepping approach to resolve collisional time scales for subsets of particles that are present at the leaves of the ORB, thereby resulting in substantial reduction of computational cost. The parallel implementation of this method where a ``knapsack” algorithm is used in tandem with ORB for effective load-balancing is also presented, where a best possible partitioning is obtained based on number of particles and local time-stepping costs. The algorithm is tested against benchmark problems with varying particle distributions that include fluidized bed and riser flow scenarios. Preliminary results indicate that the approach is 2-3X faster than traditional explicit methods for problems that involve both dense and dilute regions, while maintaining the same level of accuracy.

adaptive timestepping↗

DEM and MELT: A Modular and Compact In-Can Melting Technology Dedicated to D and D and Remediation Waste - 20033

The DEM and MELT In-Can vitrification process is designed to process intermediate and high level waste coming for D and D and remediation operations. It is developed to treat liquid and solid waste, to produce a small amount of secondary waste and to minimize investment and operating costs. This In-Can vitrification process is also developed with a compact and modular design and can be adapted regarding nuclear operators' needs and requirements. Furthermore, this robust and simple process is flexible enough to accommodate most uncertainties on waste composition. The DEM and MELT process is developed through the French consortium coordinated by CEA which gathers Orano, ECM technologies and Andra. In the frame of this project, Orano has designed a compact and modular layout of the process allowing an efficient in situ implementation, close to the waste to be treated. This paper presents the compact and modular design of the DEM and MELT process and the way it can be used for a short duration, as a decommissioning tool aimed to be dismantled just after the treatment operation. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

An Investigation into Solution Verification for CFD-DEM

This report presents the study of the convergence behavior of the computational fluid dynamicsdiscrete element method (CFD-DEM) method, specifically National Energy Technology Laboratory’s (NETL) open source MFiX code (MFiX-DEM) with a diffusion based particle-tocontinuum filtering scheme. In particular, this study focused on determining if the numerical method had a solution in the high-resolution limit where the grid size is smaller than the particle size. To address this uncertainty, fixed particle beds of two primary configurations were studied: i) fictitious beds where the particles are seeded with a random particle generator, and ii) instantaneous snapshots from a transient simulation of an experimentally relevant problem. Both problems considered a uniform inlet boundary and a pressure outflow. The CFD grid was refined from a few particle diameters down to 1/6 th of a particle diameter. The pressure drop between two vertical elevations, averaged across the bed cross-section was considered as the system response quantity of interest. A least-squares regression method was used to extrapolate the grid-dependent results to an approximate “grid-free” solution in the limit of infinite resolution. The results show that the diffusion based scheme does yield a converging solution. However, the convergence is more complicated than encountered in simpler, single-phase flow problems showing strong oscillations and, at times, oscillations superimposed on top of globally non-monotonic behavior. The challenging convergence behavior highlights the importance of using at least four grid resolutions in solution verification problems so that (over-determined) regression-based extrapolation methods may be applied to approximate the grid-free solution. The grid-free solution is very important in solution verification and VVUQ exercise in general as the difference between it and the reference solution largely determines the numerical uncertainty. By testing different randomized particle configurations of the same general problem (for the fictitious case) or different instances of freezing a transient simulation, the numerical uncertainties appeared to be on the same order of magnitude as ensemble or time averaging uncertainties. By testing different drag laws, almost all cases studied show that model form uncertainty in this one, very important closure relation was larger than the numerical uncertainty, at least with a reasonable CFD grid, roughly five particle diameters. In this study, the diffusion width (filtering length scale) was mostly set at a constant of six particle diameters. A few exploratory tests were performed to show that similar convergence behavior was observed for diffusion widths greater than approximately two particle diameters. However, this subject was not investigated in great detail because determining an appropriate filter size is really a validation question which must be determined by comparison to experimental or highly accurate numerical data. Future studies are being considered targeting solution verification of transient simulations as well as validation of the filter size with direct numerical simulation data.

97 MATHEMATICS AND COMPUTING↗

CFD-DEM study of bubble properties in a cylindrical fluidized bed of Geldart Group D particles and comparison with prior MRI data

Numerical simulations of a freely bubbling cylindrical fluidized bed are carried out using a coupled computational fluid dynamics and discrete element method (CFD-DEM) model and compared to recent experimental data. The experiments were conducted using high-resolution and high-frequency magnetic resonance imaging providing high-fidelity data of the bubbling within a central 10 mm slice of the bed. Qualitatively, we nd more regular (i.e., less chaotic) structures in the simulated beds than observed experimentally. Quantitatively, however, the bubble diameter and number of bubbles as a function of height within the bed is predicted well by the base model. Unfortunately, the regularity in the simulations manifests as a considerable discrepancy in the speed of the (dense) emulsion phase. The simulated velocity probability distribution functions show an accumulation of low-speed regions and deficiency of high-speed regions. A simple parametric study of the base model is also carried out considering many of the most common CFD-DEM modeling parameters. It is found that the fluid grid size, geometry resolution, transfer kernel and drag law did not have a significant effect on bubble or particle dynamics.

42 ENGINEERING↗

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↗

Small-scale validation tests for MFIX-Exa CFD-DEM

This report describes several bench-scale fluidization experiments that can be used to validate the CFD-DEM method as encapsulated in the MFIX-Exa code. The five cases considered are the cold-flow fluidized beds of Müller et al., Link et al. (spout-fluid), and Goldschmidt et al. (bi-disperse), the hot fluidized bed of Patil et al., and the adsorbing fluidized bed of Li et al. and Janssen. In most cases, MFIX-Exa with “standard” or “typical” CFD-DEM settings, the Gidaspow drag model, and the Gunn heat transfer provide a relatively good prediction of the quantities considered: mean void fraction profiles, mean velocity profiles, fluctuating velocity profiles, mean particle temperature and segregation index. These results, with other verification and validation tests reported elsewhere, contribute to a body of work providing confidence and credibility in CFD predictions from the MFIX-Exa code.

97 MATHEMATICS AND COMPUTING↗

Implementation of Charged Particle Behavior in Discrete Element Method (DEM) Simulations

Lunar landers will agitate the surface of the Moon with an exhaust plume during descent which will, due to the particulate nature of the lunar regolith, loosen and eject grains from the surface. This ejection is not only coupled with the charged plume gas, but also results in significant particle-particle interactions. Settling of these grains after plume effects have subsided takes much longer than expected in a ballistic trajectory. The prevailing hypothesis attributes this behavior to the accumulated charge on the particles. We are thus developing a discrete element method (DEM) approach to explore these charged particle interactions on the lunar surface. The Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) Improved for General Granular and Granular Heat Transfer Simulations (LIGGGHTS) software package provides a DEM modeling framework for granular interactions. It includes many complexities such as non-spherical particle shapes, cohesion and frictional forces, and heat transfer, but has no provision for inter-particle electrostatic forces and charge transfer that are important to examine in the lunar environment. In this work, a standard Coulomb potential and a Yukawa potential are integrated into the LIGGGHTS framework to provide a basis for particle-particle electrostatic interactions, as well as a gravitational potential to enable inter-grain gravitational interactions. A preliminary approach to charge transfer between grains incorporating properties such as work function and electrical conductivity to the library of available material characteristics will be presented. Several scenarios have been simulated that include charged particle interactions within a diffuse granular gas, settling of charged grains into a regolith bed, sliding of granular material along an incline, and vibration of settled grains to produce a behavior similar to fluidization. There are numerous challenges to incorporate realistic interactions between complex lunar particles. Currently, grains are modeled to behave as if the entirety of the charge acts at the center of mass, such as conductors with spherical symmetry and insulators with homogeneously distributed charge. We are developing improvements that will include the use of non-spherical particle geometries, as well as reasonable approximations of insulating/dielectric materials that have non-uniform charge distributions. The cases simulated thus far will be examined in a relevant environment within a vacuum chamber to validate the simulations. These simulations will be bounded by experiments utilizing high-speed camera observations of the motion for validation. The grains in the experiment will exchange charge during their motion and this can be quantified by collection within a charge measurement device such as a Faraday cup. Such a device may be modeled within the software by defining an integration region and computing the contained charge as a function of simulation time, allowing for side-by-side comparison of simulated and measured bulk charging results. Any differences will be reconciled by updating the mathematical mechanisms described within the simulation suite. Successfully combining results from experiments within a relevant environment into the LIGGGHTS framework will improve modeling of the charged grain dynamics experienced on the Moon to provide insights into dust behavior for future lunar exploration missions.

Electrostatics↗

PRODEM: An Annual Series of Summer DEMs (2019 through 2022) of the Marginal Areas of the Greenland Ice Sheet

Surface topography across the marginal zone of the Greenland Ice Sheet is constantly evolving in response to changing weather, season, climate, and ice dynamics. However, current digital elevation models (DEMs) for the ice sheet are usually based on data from a multi-year period, thus obscuring these changes over time. Here we present four 500 m resolution summer DEMs (PRODEMs) of the Greenland Ice Sheet marginal zone for 2019 through 2022. The PRODEMs cover the marginal zone from the ice edge to 50 km inland, hence capturing all Greenland outlet glaciers. Each PRODEM is based on data fusion of CryoSat-2 radar altimetry and ICESat-2 laser altimetry using regionally varying kriging of elevation anomalies relative to ArcticDEM. The PRODEMs are validated using leave-one-out cross-validation, and PRODEM19 is further validated against an external data set, showcasing their ability to correctly represent surface elevations within the associated spatially varying prediction uncertainties. We observe a general lowering of surface elevations during the 4-year PRODEM period, but the spatial pattern of change is highly complex and with annual changes superimposed. The PRODEMs enable detailed studies of the marginal ice sheet elevation changes. With their high spatio-temporal resolution, the PRODEMs will be of value to a wide range of researchers and users studying ice sheet dynamics and monitoring how the ice sheet responds to changing environmental conditions. PRODEMs from summer 2019 through 2022 are available at https://doi.org/10.22008/FK2/52WWHG (Winstrup, 2024), and we plan to annually update the product henceforth.

Digital Elevation Models↗