Discrete element modeling of a particle heater for energy storage systems
Not Available
Schirck, Jason↗
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.
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↗
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↗
New electric force and charge exchange modules in discrete element model enables particle dynamics simulation in electric field
Discrete element modeling (DEM) is an important technique for particle dynamics simulation. The field of metal additive manufacturing often utilizes DEM to simulate the rheological behaviors of powder. Standard contact and short-range interactions are sufficient in most cases but insufficient to describe the particle dynamics with the influence of an electric field. Modeling such a system requires additional physics to describe the particle–field interactions. The relevant physics has been experimentally understood but is not yet available in DEM. Here, we develop a charge exchange and an electric force module. The electric force module governs particle response to the electric field, while the charge exchange module enables particles to acquire proper charge during contact with charged geometries. We validate the modules against analytical calculations and high-speed videos of electrostatic powder deposition experiments. Notably, the model struggles to capture the initial particle levitation. We later deploy a modified electric field, as supported by static electric field simulation, to better approximate the electric field penetration into the powder layer. This modification improves the model’s capability of simulating realistic particle levitation. The results highlight the challenges of modeling particle behaviors in the electric field while demonstrating the feasibility of obtaining quantitative results, which are difficult to measure experimentally.
charge exchange↗
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↗
A high-resolution pseudo-polygon discrete element model for regional sea ice
Here, this work presents a pseudo-polygon discrete element model for high-resolution sea ice simulations. A scale-invariant bonded particle contact model is proposed to model joints between sea ice floes based on the smeared fracture model and a lattice spring beam model, and the Mohr–Coulomb failure criterion is implemented to represent the shearing failure mechanism of sea ice packings under complex loadings. All mechanical parameters of the bond model can be directly determined from laboratory tests. Validations of the proposed model are made by investigations of mechanical response and failure criteria of field sea ice sheets. Compared with the field observations of sea ice from satellite radar and in situ stress sensors, the proposed model is capable of reproducing the typical constitutive behavior and the Coulomb friction envelope of field sea ice. Finally, the proposed discrete element sea ice model is used to study the effect of loading rates on mechanical behavior including failure strength of regional sea ice.
54 ENVIRONMENTAL SCIENCES↗
Machine learning accelerated discrete element modeling of granular flows
Granular flows are widely encountered in many industrial processes and natural phenomena. Discrete Element Modeling (DEM) is a useful tool for understanding and troubleshooting devices, which handle granular materials. However, its applicability is significantly limited by the huge computational cost associated with detecting and computing collisions. In this research, the computation speed of DEM was accelerated by orders of magnitude using a convolutional neural network to replace the direct calculation of particle-particle and particle-boundary collisions. The MFiX software was used to generate the training and testing dataset. Additionally, a GPU accelerated TensorFlow model was used to train the neural network and test the results. The model fluctuations caused by different training steps were reduced with a multi-scale loss function. The accuracy was improved with more frames within one training step. The modeling of a rotating drum and a hopper demonstrated the accuracy and efficiency of this machine learning accelerated DEM in the simulation of granular flows.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗
Spiral waves are stable in discrete element models of two-dimensional homogeneous excitable media
The spontaneous breakup of a single spiral wave of excitation into a turbulent wave pattern has been observed in both discrete element models and continuous reaction-diffusion models of spatially homogeneous 2D excitable media. These results have attracted considerable interest, since spiral breakup is thought to be an important mechanism of transition from the heart rhythm disturbance ventricular tachycardia to the fatal arrhythmia ventricular fibrillation. It is not known whether this process can occur in the absence of disease-induced spatial heterogeneity of the electrical properties of the ventricular tissue. Candidate mechanisms for spiral breakup in uniform 2D media have emerged, but the physical validity of the mechanisms and their applicability to myocardium require further scrutiny. In this letter, we examine the computer simulation results obtained in two discrete element models and show that the instability of each spiral is an artifact resulting from an unphysical dependence of wave speed on wave front curvature in the medium. We conclude that spiral breakup does not occur in these two models at the specified parameter values and that great care must be exercised in the representation of a continuous excitable medium via discrete elements.
NASA Discipline Cardiopulmonary↗
A hybrid calibration approach to Hertz-type contact parameters for discrete element models
This study aims at providing a hybrid calibration framework to estimate Hertz-type contact parameters (particle-scale shear modulus and Poisson ratio) for both two-dimensional and three-dimensional discrete element modelling (DEM). On the basis of statistically isotropic granular packings, a set of analytical formulae between macroscopic material parameters (Young modulus and Poisson ratio) and particle-scale Hertz-type contact parameters for granular systems are derived under small-strain isotropic stress conditions. However, the derived analytical solutions are only estimated values for general models. By viewing each DEM modelling as an implicit mathematical function taking the particle-level parameters as independent variables and employing the derived analytical solutions as the initial input parameters, an automatic iterative scheme is proposed to obtain the calibrated parameters with higher accuracies. Considering highly nonlinear features and discontinuities of the macro-micro relationship in Hertz-based discrete element models, the adaptive moment estimation algorithm is adopted in this study because of its capacity of dealing with noise gradients of cost functions. Here, the proposed method is validated with several numerical cases including randomly distributed monodisperse and polydisperse packings. Noticeable improvements in terms of calibration efficiency and accuracy have been made.
Constitutive law↗
Discrete Element Modeling (DEM) of Triboelectrically Charged Particles: Revised Experiments
In a previous work, the addition of basic screened Coulombic electrostatic forces to an existing commercial discrete element modeling (DEM) software was reported. Triboelectric experiments were performed to charge glass spheres rolling on inclined planes of various materials. Charge generation constants and the Q/m ratios for the test materials were calculated from the experimental data and compared to the simulation output of the DEM software. In this paper, we will discuss new values of the charge generation constants calculated from improved experimental procedures and data. Also, planned work to include dielectrophoretic, Van der Waals forces, and advanced mechanical forces into the software will be discussed.
Hogue, Michael D.↗
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↗
Polysiloxane coatings on microspheres based on Multiphase Flow with Interface Exchange-Discrete Element Modelling
Polysiloxane coatings on yttria stabilized zirconia (YSZ) microspheres of 500 μm were simulated in a spout fluidized bed coater using Multiphase Flow with Interface Exchange-Discrete Element Modelling (MFiX-DEM). In this work, two different coater configurations were developed to study the influence of gas velocity and its distribution on particle dynamics. The presence of the Wurster tube not only enhances the distribution but also increases the overall residence time of the particles. Investigations were also carried out with different Wurster tube positions (normal, 10% and 20% lowered from its initial position). Among these, 20% lowered Wurster tube position demonstrated the most effective coating process. The effects of gas inlet pressure on the average gas velocity and the distribution of particles were analyzed. More than 97% of the particles can be retained. The derived results, including average gas velocity, particle retention percentage, and distribution of particles with gas velocity, are being used to guide the experimental work in obtaining defect-free coatings for YSZ microspheres.
36 MATERIALS SCIENCE↗
Discrete element modeling of irregular-shaped soft pine particle flow in an FT4 powder rheometer
Pine residues are a commonly used biomass feedstock that consists of different anatomical fractions, each with distinct particle characteristics. Here, in this work, an experiment-informed discrete element model is developed to investigate the flowability of pine residues in an FT4 rheometer. Multi-sphere particles with distinct particle attributes are created to model each anatomical fraction type. A systematic analysis of specimens with varying particle characteristics (e.g., anatomical fraction type, particle shape, and size) is conducted to elucidate the relationship between particle attributes and flowability. The results show that stems recorded the highest axial force and torque and, correspondingly, the highest flow energy, which is attributed to their high stiffness and interlocking effect. Increasing their percentage in the mixture increases the flow energy while increasing needles tends to decrease flow energy. Knowledge gained in this study on the flow of anatomical fractions is important for the efficient and robust processing of pine residues.
09 BIOMASS FUELS↗
Determining Trajectory of Triboelectrically Charged Particles, Using Discrete Element Modeling
The Kennedy Space Center (KSC) Electrostatics and Surface Physics Laboratory is participating in an Innovative Partnership Program (IPP) project with an industry partner to modify a commercial off-the-shelf simulation software product to treat the electrodynamics of particulate systems. Discrete element modeling (DEM) is a numerical technique that can track the dynamics of particle systems. This technique, which was introduced in 1979 for analysis of rock mechanics, was recently refined to include the contact force interaction of particles with arbitrary surfaces and moving machinery. In our work, we endeavor to incorporate electrostatic forces into the DEM calculations to enhance the fidelity of the software and its applicability to (1) particle processes, such as electrophotography, that are greatly affected by electrostatic forces, (2) grain and dust transport, and (3) the study of lunar and Martian regoliths.
Source record↗
Discrete element modeling of switchgrass particles under compression and rotational shear
Switchgrass is a perennial herbaceous plant regarded as a biomass energy crop in the United States for its highadaptability and yield potential. Processing and handling of switchgrass particles are challenging due to the erratic mechanical and flow behavior originating from their intrinsic particulate properties. Here, we present a bonded-sphere discrete element model designed specifically for switchgrass particles. The model simultaneously captures three key particulate features, i.e., fibrous particle shapes, a wide range of particle sizes, and particle deformability. Realistic yet computationally efficient particle shape templates are created based on the image analysis data of switchgrass specimens. A fitting procedure is proposed to ensure both the particle width and length distributions are captured, a unique requirement for fibrous particles. Two full-scale numerical models, i.e., a uniaxial compression model and a Schulze ring shear model, are developed using information fromphysical experiments. The model is calibrated using experimental data of chopped-small switchgrass specimens, and then, is validated using data of chopped-large specimens in both compression and ring-shear tests. Numerical results show that the numerical models capture bulk densities accurately (with an error of 3%) while slightly underestimate the bulk friction angle. Furthermore, an extensive sensitivity analysis reveals that (1) switchgrass particles with rougher edges (due to different processing techniques) exhibit a higher shear strength and a lower flowability; (2) stiffer particles yield a lower bulk density (up to 21% lower) compared to more deformable particles, indicating particle deformability should be incorporated when modeling biomass flow in a preprocessing system.
09 BIOMASS FUELS↗