Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Geometric Scale-up Experiments on Fluidization of Geldart B Glass Beads

The objective of this work is to provide a valuable database from controlled experiments for validating computational models. Recently, coarse-grained techniques such as particle-in-cell (PIC) or coarse-grained discrete element modeling (DEM) have gained popularity due to their computational efficiency while modeling large-scale systems; however, the influence of model parameters and their sensitivities at different geometric scales and flow conditions remain to be analyzed. These datasets are critical for the multiphase flow research community to assess predictive capability of modeling techniques as well as elucidate the hydrodynamic behavior in these systems. This study performed fluidization experiments using three different test sections with internal diameters of 2.5, 4, and 6 in. The operating conditions, bed material, and range of flow velocities at the inlet were constant in all the units, which were not hydrodynamically scaled. Glass beads having a Sauter Mean Diameter of 332 μm were used. Superficial velocity was varied from 2.97 to 5.35 times the minimum fluidization velocity and the initial static bed height was 0.1524 m. The order in which the experiments were performed was randomized and replicates were included to provide uncertainty in measurements. Statistics of differential pressure and bed height from these tests were reported. Future plans include validating PIC methodology in the open-source software, MFiX (Multiphase Flow with Interphase Exchanges) using results from this study. This could further be extended to determine optimal model parameters using inverse techniques such as deterministic calibration or Bayesian inference.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2020

Geometric Scale-up Experiments on Fluidization of Geldart B Glass Beads

The objective of this work is to provide a valuable database from controlled experiments for validating computational models. Recently, coarse-grained techniques such as particle-in-cell (PIC) or coarse-grained discrete element modeling (DEM) have gained popularity due to their computational efficiency while modeling large-scale systems; however, the influence of model parameters and their sensitivities at different geometric scales and flow conditions remain to be analyzed. These datasets are critical for the multiphase flow research community to assess predictive capability of modeling techniques as well as elucidate the hydrodynamic behavior in these systems. This study performed fluidization experiments using three different test sections with internal diameters of 2.5, 4, and 6 in. The operating conditions, bed material, and range of flow velocities at the inlet were constant in all the units, which were not hydrodynamically scaled. Glass beads having a Sauter Mean Diameter of 332 μm were used. Superficial velocity was varied from 2.97 to 5.35 times the minimum fluidization velocity and the initial static bed height was 0.1524 m. The order in which the experiments were performed was randomized and replicates were included to provide uncertainty in measurements. Statistics of differential pressure and bed height from these tests were reported. Future plans include validating PIC methodology in the open-source software, MFiX (Multiphase Flow with Interphase Exchanges) using results from this study. This could further be extended to determine optimal model parameters using inverse techniques such as deterministic calibration or Bayesian inference.

42 ENGINEERING↗
DOE OSTI2021

Assessment of model parameters in MFiX particle-in-cell approach

The limitations in numerical treatment of solids-phase in conventional methods like Discrete Element Model and Two-Fluid Model have facilitated the development of alternative techniques such as Particle-In-Cell (PIC). However, a number of parameters are involved in PIC due to its empiricism. In this work, global sensitivity analysis of PIC model parameters is performed under three distinct operating regimes common in chemical engineering applications, viz. settling bed, bubbling fluidized bed and circulating fluidized bed. Simulations were performed using the PIC method in Multiphase Flow with Interphase eXchanges (MFiX) developed by National Energy Technology Laboratory (NETL). A non-intrusive uncertainty quantification (UQ) based approach is applied using Nodeworks to first construct an adequate surrogate model and then identify the most influential parameters in each case. This knowledge will aid in developing an effective design of experiments and determine optimal parameters through techniques such as deterministic or statistical calibration.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2022

Correction to: Investigating the rheology of fluidized and non-fluidized gas-particle beds: implications for the dynamics of geophysical flows and substrate entrainment

We report natural geophysical mass flows are among the most complex granular systems and their dynamics are often modified by the presence of an interstitial fluid. Prediction of their runout requires the development of models estimating the solid stresses in these hazardous currents wherein excess pore-fluid pressure can develop. We use discrete element modelling (DEM-CFD) with a Coarse-Graining post-processing technique (CG) to investigate the rheology of unsteady gas-particle fluidized to non-fluidized granular beds placed on horizontal and inclined planes. Similar to fluidized beds immersed in viscous fluids, the effective friction coefficient of air-fluidized beds can be defined as a function of the classic μ(I)-rheology and the non-dimensional fluid or solid pressure to explain the failure and dynamics of granular flows with excess pore pressure on inclines. However, dilation imposed by fluid drag and particle collisions in gas-particle fluidized beds can drastically change its effective frictional properties. In contrast with the common assumption in water-particle flows that granular temperature is negligible, in our gas-particle simulations, the contribution of the velocity fluctuations to the stress tensor is significant. Hence, the shear stress is found to be non-zero even when the flow is fully fluidized in the inertial regime. These results suggest the need to better understand velocity fluctuations to predict the effective viscosity of sheared fluidized mixtures and are relevant for many applications. Notably, a unified approach is useful for many geophysical flows that encompass a range of fluidization conditions in a single flow such as pyroclastic density currents and snow avalanches.

58 GEOSCIENCES↗
DOE OSTI2020

Investigating biomass composition and size effects on fast pyrolysis using global sensitivity analysis and CFD simulations

It is notoriously difficult to build an accurate universal model for biomass pyrolysis due to its sensitivity to a wide number of critical material attributes such as chemical species and physical sizes. In this work, a biomass pyrolysis kinetics with 32 heterogeneous reactions and 59 species was implemented in an open-source multiphase computational fluid dynamics (CFD) software MFiX and validated against two different experimental pyrolysis data sets that provided detailed data describing chemical component yields. The reaction scheme was then used to build a surrogate model and assess the sensitivity of pyrolysis yields to feedstock compositions. The sensitivity analysis determined that the yield of bio-char showed a strong positive sensitivity to the carbon-rich lignin and tannin pseudo-species in the reaction scheme while the bio-oil and bio-gas were correlated to oxygen-rich lignin pseudo-species. The reaction scheme was then integrated into a coarse-grained discrete element model to simulate fast pyrolysis in a bubbling fluidized bed over a range of feedstock particle sizes. The reactor simulations showed further sensitivity to particle size and hydrodynamics. Notably, particles under 0.5 mm have small heat transfer limitations but left the reactor before completely converting and thus reduced the bio-oil yield. Results from this study can be used to guide future development of highly accurate models for fast pyrolysis reactors with a variety of feedstock properties and operating conditions.

09 BIOMASS FUELS↗
DOE OSTI2023

Fabric controls on fracture surface roughness of an architected rock material

Fluid flow through fractures is intimately linked to the fracture surfaces that define the void geometry through which fluids flow. Thus, an understanding of what controls fracture surface roughness is essential to the development of models for predicting fluid transport through fractured rock. The difficulty in predicting surface roughness arises from the complexity of rock which is inherently heterogeneous and nonuniform in composition, fabric, and structural components, even when samples are acquired from the same rock mass. Here, a benchmarked-simulation approach motivated from geo-architected 3D printed synthetic gypsum rocks is used to provide insight into the competing contributions from fabric and layering on fracture roughness formation. Simulation results from a discrete element model (Particle Flow Code, Itasca Consulting Group, Inc.) clearly indicate that the relative orientation between mineral layers and in-layer mineral fabric, and the variability in mineral bonding strengths determine whether anisotropic corrugated surfaces or isotropic surfaces are formed. Weak mineral layers oriented perpendicular to the applied load resulted in strong roughness anisotropy. Peak failure loads were found to vary up to 30% depending on the strength of the mineral fabric at the location of fracture initiation, which provides insight into the observed high variability in strength values of natural rock. The uniqueness of induced fracture roughness and peak failure load is intimately linked to layering, mineral fabric, and their distribution in the rock. These findings have important implications for any architected material fabricated through serial printing of layers with local compositional heterogeneity.

3D printed rock↗
DOE OSTI2023

Particle size and shape effect of Crumbler® rotary shear-milled granular woody biomass on the performance of Acrison® screw feeder: A computational and experimental investigation

Physical experiments and discrete element model (DEM) simulations are conducted to evaluate particle characteristics and operation parameter effects on screw feeding performance for rotary shear-milled Douglas fir. Three performance metrics are used: mass flow rate, shaft driving torque, and specific energy consumption. The impact of particle size, particle size distribution (PSD), shaft rotational speed (rpm), and hopper dimensions on the performance are investigated. All employed performance metrics reveal the superior flowability of the 2-mm particles in contrast to the larger 6-mm counterpart. Remarkably, wider PSD results in poorer flowability than the two mono-sized particles, proving the flowability enhancement achieved by narrower PSD of the woody feedstock. More importantly, DEM simulations unveil PSD-induced degradation in flowability is attributed to mechanical interlocking and particle segregation effects. Furthermore, higher shaft rpm causes higher mass flow rate at the cost of higher specific energy consumption due to viscous dissipation and changes in flow pattern.

09 BIOMASS FUELS↗
DOE OSTI2023

Joint Geophysical and Numerical Insights of the Coupled Thermal‐Hydro‐Mechanical Processes During Heating in Salt

Abstract Salt offers an optimal medium for the permanent isolation of heat‐producing radioactive waste due to its impermeability, high thermal conductivity, and ability to close fractures through creep. A thorough understanding of the thermal‐hydrological‐mechanical (THM) processes, encompassing brine migration, is fundamental for secure radioactive waste disposal within salt formations. At the Waste Isolation Pilot Plant (WIPP), we conducted joint in situ geophysical monitoring experiments during active heating to investigate brine migration near excavations. This experiment incorporated electrical resistivity tomography (ERT) alongside high‐resolution fiber‐optic‐based distributed temperature sensing within a controlled heating experiment. Additionally, discrete element model (DEM) based numerical simulations were conducted to simulate THM processes during heating, providing a more mechanistic understanding of the coupled processes leading to the observed changes in geophysical measurements. During heating, resistivity shifts near the heater were reasonably explained by temperature effects. However, in more distant, cooler regions, the resistivity decrease exceeded predictions based solely on temperature. DEM simulations highlighted brine migration, propelled by pore pressure gradients, as the likely primary factor contributing to the additional resistivity decline beyond temperature effects. The comparison between the predicted ERT responses and observations was much improved when considering the effects of brine migration based on the DEM simulations. These geophysical and simulation findings shed light on brine migration in response to salt heating, enhancing our understanding of the coupled THM processes in salt for safe radioactive waste disposal.

58 GEOSCIENCES↗
DOE OSTI2021

StressNet - Deep learning to predict stress with fracture propagation in brittle materials

Abstract Catastrophic failure in brittle materials is often due to the rapid growth and coalescence of cracks aided by high internal stresses. Hence, accurate prediction of maximum internal stress is critical to predicting time to failure and improving the fracture resistance and reliability of materials. Existing high-fidelity methods, such as the Finite-Discrete Element Model (FDEM), are limited by their high computational cost. Therefore, to reduce computational cost while preserving accuracy, a deep learning model, StressNet, is proposed to predict the entire sequence of maximum internal stress based on fracture propagation and the initial stress data. More specifically, the Temporal Independent Convolutional Neural Network (TI-CNN) is designed to capture the spatial features of fractures like fracture path and spall regions, and the Bidirectional Long Short-term Memory (Bi-LSTM) Network is adapted to capture the temporal features. By fusing these features, the evolution in time of the maximum internal stress can be accurately predicted. Moreover, an adaptive loss function is designed by dynamically integrating the Mean Squared Error (MSE) and the Mean Absolute Percentage Error (MAPE), to reflect the fluctuations in maximum internal stress. After training, the proposed model is able to compute accurate multi-step predictions of maximum internal stress in approximately 20 seconds, as compared to the FDEM run time of 4 h, with an average MAPE of 2% relative to test data.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Granular packings with sliding, rolling, and twisting friction

Intuition tells us that a rolling or spinning sphere will eventually stop due to the presence of friction and other dissipative interactions. The resistance to rolling and spinning or twisting torque that stops a sphere also changes the microstructure of a granular packing of frictional spheres by increasing the number of constraints on the degrees of freedom of motion. We perform discrete element modeling simulations to construct sphere packings implementing a range of frictional constraints under a pressure-controlled protocol. Mechanically stable packings are achievable at volume fractions and average coordination numbers as low as 0.53 and 2.5, respectively, when the particles experience high resistance to sliding, rolling, and twisting. Only when the particle model includes rolling and twisting friction were experimental volume fractions reproduced.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Fluctuations and power-law scaling of dry, frictionless granular rheology near the hard-particle limit

Here the flow of frictionless granular particles is studied with stress-controlled discrete element modeling simulations for systems varying in size from 300 to 100,000 particles. The volume fraction and shear-stress ratio μ are relatively insensitive to system size for a wide range of inertial numbers I. Second-order effects in strain rate, such as normal stress differences, require large system sizes to accurately extract meaningful results, notably a nonmonotonic dependence in the first normal stress difference with strain rate. The rheological response represented by the μ(I) scalar model works well at describing the lower-order aspects of the rheology, except near the quasistatic limit of these stress-controlled flows. The pressure is varied over five decades, and a pressure dependence of the coordination number is observed, which is not captured by the inertial number. Large fluctuations observed for small systems N≤1000 near the quasistatic limit can lead to the arrest of flow resulting in challenges to fitting the data to rheological relationships. The inertial number is also insufficient for capturing the pressure-dependent behavior of property fluctuations. Fluctuations in the flow and microstructural properties are measured in both the quasistatic and inertial regimes, including shear stress, pressure, strain rate, normal stress differences, volume fraction, coordination number, and contact fabric anisotropy. The fluctuations in flow properties scale self-similarly with pressure and system size. A transition in the scaling of fluctuations of stress properties and contact fabric anisotropy are measured and proposed as a quantitative identification of the transition from inertial to quasistatic flow.

36 MATERIALS SCIENCE↗
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 OSTI2022

Viewgraphs for w20_seismicsources [Slides]

Hybrid Optimization Software Suite (HOSS) is based on finite-discrete element modeling to capture the complex physics of explosions. Calibration of near-source modeling for underground explosions in geophysical models is presented, showing the effects of the amplitude of radial and vertical seismic components.

36 MATERIALS SCIENCE↗
DOE OSTI2021

MFIX DEM Enhancement for Industry-Relevant Flows (Final Report)

The overall goal of this two-phase project is to implement performance improvements of the Multiphase Flow with Interphase Exchanges (MFIX) Discrete Element Model (DEM) code that enable a transformative shift for industrial use. Prior to this effort, the largest simulations performed using MFIX are O(10 7 ) particles. This falls short of the O(10 9 ) particle simulations that must be completed on a timescale of days or weeks (vs. months or years) to enable simulations with physically-relevant domain sizes to be incorporated into industrial design cycles within five years. This was accomplished by tailoring best-in-class practices to bear on the unique challenges posed by the MFIX-DEM algorithm and code base. Scientific simulations (e.g., in cosmology, turbulent combustion) routinely use massively parallel computing to update far more particles in short wall clock times. Results from Phase 1 (1.5 years in duration) indicated significant gains in speed were possible for a wide range of benchmark cases. Moreover, a survey sent to >35 companies indicates that the timing is ideal for such an enhanced tool, with >80% of the respondents indicating that DEM is already value-added or will be within the next 5 years, and >70% of the respondents indicating that improved speed is the top computational priority. In Phase 2 (3.5 years in duration), the two major barriers that hinder industry from effectively using multiphase Computational Fluid Dynamics (CFD) to cut costs and improve performance, namely computational overhead and confidence in predictions, continued to be addressed. Regarding the former, the results from Phase 1 to guide the effort, with enhancements focused on an improved time-stepping algorithm and particle sorting. Four target problems of 1 billion particles each and increasing complexity were identified: homogeneous cooling, tumbler with continuous particle size distribution, discharge from a rectangular hopper and a cylindrical riser. Each of these were successfully simulated for relevant time scales (on order of seconds) using less than 24 hours of wall clock time. These represent the first 1-billion particle DEM simulations performed with MFIX, namely using the MFIX-Exa code. This code is currently under development at NETL in collaboration with Lawrence Berkeley National Laboratory. Regarding the second barrier on predictive uncertainty, experiments from Phase 1 (interacting nozzles - hydrodynamics only) and Phase 2 (very small-scale segregation experiments) were used to demonstrate the ability of two simplified approaches to uncertainty quantification (UQ). By limiting the number of particles, UQ based on the simplified treatment was compared to standard UQ, which was shown to have much higher computational demands. Experiments were also performed on a pilot-scale stripper unit to provide validation data for future CFD-DEM simulations and UQ.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2025

Modeling the formation of Sedan Crater using the FLAG and HOSS codes

Numerical modeling of explosion crater formation requires accounting for complex physical processes. Numerical validation of explosion cratering is an important step in modeling and requires experimental data for comparison. Models using discrete elements and continuum models have both benefits and drawbacks to their approaches. In this work, we consider both an arbitrary Lagrangian–Eulerian (ALE) hydrocode and a finite discrete element method (FDEM) approach to modeling the formation of the Sedan crater, the largest human-made crater in the United States. The Sedan crater formed from an underground nuclear detonation in the Nevada desert as part of Project Plowshare. Our models show that the continuum approach of the hydrocode matched well compared to early test time prior to the mound rupture and subsequent fireball venting, when most of the alluvium exhibited fluid behavior. Our FDEM approach matched the final crater dimensions well, after material had settled back into the crater, when material strength and solid mechanics play key roles. Our work shows how leveraging the benefits of multiple numerical approaches can lead to better understanding of complex physical problems, especially problems with limited experimental data. By using a continuum approach to early-time hydrodynamics and an FDEM approach to later-time solid mechanics, we can better understand the different physical regimes of explosion crater formation.

36 MATERIALS SCIENCE↗