A bounding box search algorithm for DEM simulation

Effect of variations in atelectasis on tumor displacement during radiation therapy for locally advanced lung cancer

A pseudo-coupled DEM–non-linear PBM approach for simulating the evolution of particle size during dry milling

HighlightsThe relationship between the angle of repose and compound fertilizer particle size distribution and sphericity was revealed.The range of the rolling friction coefficient obtained by different modeling was calibrated by the golden section method.Scanning modeling required 80% of the rolling friction coefficient of spherical modeling to reflect actual motion.The relative error of scanning modeling (7.91%) was better than that of spherical modeling (30.84%).Abstract. The discrete element method (DEM) is widely used to simulate the behavior of granular materials. The accurate establishment of the particle model has a great influence on the accuracy of the simulation. Most of the existing discrete element simulations of granular compound fertilizers use spherical modeling without considering the shape or particle size distribution. To compare the simulation accuracy of spherical modeling (spherical particles with equivalent diameters) and scanning modeling (a three-dimensional model of fertilizers was obtained by a scanner and then automatically filled by EDEM2020 software) on compound fertilizer discharge, three kinds of granular compound fertilizers (Changqingshu, Xinshidai, and Munvhong) were selected as the research objects. The characterization parameters of fertilizer shape and fluidity were measured. Principal component analysis and Pearson correlation analysis were used to analyze the influence of the compound fertilizer characterization parameters on the angle of repose. The results showed that the particle shape and particle size distribution of compound fertilizer were the main factors affecting its angle of repose. Through the golden section method and a single-factor experiment, the rolling friction coefficients of the three granular fertilizers under the two granular modeling methods for spherical modeling and scanning modeling were determined. The results showed that the rolling friction coefficient when using scanning modeling to approximate the actual angle of repose was approximately 80% of that when using the spherical modeling method. The influence of different modeling methods on the total rotational kinetic energy and the formation of the angle of repose was discussed. For the same fertilizer, the total rotational kinetic energy was greater in spherical modeling than in scanning modeling, and the formation time of the angle of repose was longer in scanning modeling than in spherical modeling. The results of an external tank wheel fertilization simulation and bench test showed that considering the scanning modeling of fertilizer shape and particle size distribution, the fertilization variation coefficient obtained from the simulation test was closer to that of the bench test. This research helps to better understand the influence of particle shape and particle size distribution on DEM simulations and provides references for discrete element modeling of other granular fertilizers. Keywords: Compound fertilizer, Modeling method, Shape, Simulated flow characteristics, Size distribution.

Modelling and simulation of fruit drop tests by discrete element method

Unwanted grain commingling impedes new quality-based grain handling systems and has proven to be an expensive and time-consuming issue to study experimentally. Experimentally validated models may reduce the time and expense of studying grain commingling while providing additional insight into details of the particle flow. In this study, grain commingling in a pilot-scale bucket elevator boot was first modeled in three-dimensional (3-D) discrete element method (DEM) simulations. Experiments on the pilot-scale boot were performed using red-dyed and clear (uncolored) soybeans to validate the 3-D DEM model. Predicted results from the 3-D boot model generally followed the experimental data but tended to underpredict commingling early in the process. To reduce computational time, quasi-two-dimensional (quasi-2-D) DEM simulations were also evaluated. Comparison of predicted average commingling of five quasi-2-D boot models with reduced control volumes (i.e., with boot widths from four to seven times the mean particle diameter) led to the selection of the quasi-2-D model with a boot width of 5.6 times the mean particle diameter (i.e., 5.6d) to reduce computational time. In addition, the 3-D and quasi-2-D (5.6d) models were refined by accounting for the initial surge of particles at the beginning of each test and correcting for the effective dynamic gap between the bucket cups and the boot wall. The quasi-2-D (5.6d) models reduced simulation run time by approximately 70% compared to the 3-D model of the pilot-scale boot. Results of this study can be used to accurately predict commingling levels and improve grain handling, which can help farmers and grain handlers reduce costs and maintain grain purity during transport and export of grain.

Improvement of a pharmaceutical powder mixing process in a tote blender via DEM simulations

Biomass has potential as a carbon-neutral alternative to petroleum for chemical and energy products. However, complete replacement of fossil fuel is contingent upon efficient processes to eliminate undesirable characteristics of biomass, e.g., low bulk density, variability, and storage-induced quality problems. Mechanical size reduction via comminution is a processing operation to engineer favorable biomass flowability in handling. Crumbler rotary shear mill has been empirically demonstrated to produce more uniformly shaped particles with higher flowability than hammermilled biomass. This study combines modeling and experimentation to unveil fundamental understandings of the relation between granular particle characteristics and biomass flow behavior, which elucidate underlying mechanisms and guide selection of critical processing parameters. For this purpose, the impact of critical material attributes, including particle size (2–6 mm), particle shape (briquette, chip, clumped-sphere, cube, etc.), and surface roughness, on the angle of repose (AOR) of milled pine chips were investigated using discrete element method (DEM) simulations. Forest Concepts Crumbler rotary shear system is used to produce milled pine particles within the same size range considered in DEM simulations. AOR of different sets of these particles were measured experimentally to benchmark DEM results against experimental data. Specific energy consumption for the comminution of biomass with different particle size and moisture content are measured for technoeconomic analysis. Our results show that the smaller size (2 mm) of pine particle achieves better followability (i.e., smaller AOR) while the energy cost of comminution is significantly higher and bulk density is almost the same as the 6-mm pine particles. For the 2-mm particle size, Crumbles from veneer have better flow properties than Crumbles from chips. Contrarily, no significant difference was observed between the AOR of the two materials for the 6-mm particle size. Furthermore, from DEM simulations, mechanical interlocking between particles was found as a dominant factor in determining AOR of complex-shaped particles such as milled pine, which cannot be accurately captured by using simple particle shapes (e.g., mono-sphere) with a rolling resistance model. Conversely, clumped-sphere model alleviates this limitation without increasing computational cost significantly and can be used for accurate representation of biomass granular particles when simulating free-flow behavior.

The presence of overlarge gravel particles poses significant challenges for laboratory testing on prototype gravelly soils due to sample size limitations. To address this issue, replacement techniques, such as substituting overlarge particles with finer materials, offer practical solutions. However, the impact of these techniques on the mechanical behavior of gravelly soils, particularly shear strength and stiffness, remains poorly understood. This study aims to bridge this knowledge gap by investigating the particle size effect on the shear behaviors of binary mixtures using a series of Discrete Element Method (DEM) simulations. Updated scaling relations, based on Iai's generalized scaling relations, were proposed to correct for particle size effects. DEM simulations, including drained triaxial tests and shear modulus measurements, were performed to validate the proposed law. The results indicate that the gravel replacement technique has a minor effect on peak shear strength but significantly reduces soil stiffness, especially at high gravel contents. The scaling relations effectively correct for the particle size effect, enabling the accurate prediction of shear behaviors of the prototype gravelly soils from those of the model gravelly soils. These validations demonstrate that for addressing the soil deformation problem instead of the stability problem in ultimate state, the developed scaling relations are highly effective for correcting the particle size effect. Based on the developed scaling relations, engineers can predict prototype-scale shear behaviors of gravelly soils with overlarge particles using scaled laboratory models, reducing reliance on costly large-scale equipment. Additionally, future studies, through both DEM simulations and laboratory experiments, are recommended to further validate and refine the proposed method across diverse soil conditions and loading scenarios, such as cyclic loadings.

In simple shear testing, the specimen boundaries play a pivotal role in the transmission of shear forces. Waffle-style porous stones, plates with ribs or similar types of projections are used in experiments to reduce slippage at the boundary and transmit shear throughout the specimen. Conventionally, Discrete Element Method (DEM) simulations often model the top and bottom caps as flat boundaries with artificially enlarged friction coefficients, or use geometrical configurations that are computationally efficient. This research compares flat high-friction boundaries with ribbed and more novel-designed boundaries incorporating large pyramid, and small pyramid projections, aiming to improve shear transmission capability while ensuring computational efficiency. DEM simulations were conducted on specimens of steel bearings, with experimental validation using identical setups and particle properties. Specimens featuring different boundaries showed different void ratios in the layers close to the boundaries, with boundary effects diminishing towards the central zones. DEM simulations with projection boundaries demonstrated good agreement with experiments in terms of the macroscopic response, validating the effectiveness of the projection boundaries. The conventional flat boundaries exhibited limited shear transmission capability, resulting in insufficient development of shear stress and inadequate particle engagement. Conversely, projections on boundaries significantly improved shear stress transmission, ensuring the simple shear condition throughout the entire specimens. Projection boundaries introduced manageable increases in computational cost despite increased mesh complexity. This study highlights the importance of boundary design and recommends the adoption of projection-based boundaries in both experimental and numerical simple shear tests to ensure effective shear transmission and reliable results.

Discrete element method simulations and experimental study of interactions in 3D granular bedding during low-velocity impact

This study investigates the physical parameters that affect the flow patterns of soybeans with various moisture content (12% to 60%) at varying orifice sizes (20, 40, and 60 mm) in a cylindrical silo. The flow conditions required to obtain a steady mass flow during discharge were evaluated via experiments and three-dimensional discrete element method (DEM) simulation. The discharged mass flow rates at different flow conditions provided the critical size of the orifice. If the reduced diameter (Dred) of an orifice is >5.59, the flow showed a steady state. Based on the mass flow index (MFI), the flow patterns at 40% and 60% moisture content at 40 and 60 mm orifice sizes, respectively, showed funnel flows. although these flow conditions were satisfied to maintain a steady flow. The maximum wall pressure for the funnel flow showed the location of the interlocking phenomenon where the stagnant zone began during discharging. DEM simulation was validated through the mass profiles using the parameters obtained by the experiments. This study demonstrates that the experimental and analytical results with DEM simulation predict the flow behaviors of soybeans well at various moisture contents. These results are useful for designing silos for continuous food processing.

Discrete element method (DEM) is an interesting alternative to classical approaches as the finite element method (FEM) to simulate homogeneous and heterogeneous materials. Indeed, although DEM was initially developed to simulate granular systems in motion, it also enables to model a continuous medium with the help of a dense granular packing in which the cohesion is introduced between each pair of particles in contact using beam or spring elements. However, among other issues, the local stress field obtained using DEM is heterogeneous, even in the case where it is theoretically homogeneous. In the present contribution, we investigate the stress field distribution in 3D DEM mechanical simulation, and propose an approach we named Halo, to evaluate the level of stress distribution. The idea of the method is to evaluate the stress at the scale of every discrete element (DE) taking into account the contribution of its neighbors, inside the Halo of the DE. The effect of the number of DEs per Halo on the stress distribution is discussed in homogeneous and heterogeneous media. The approach results for homogeneous media are compared to FEM calculations and to the theoretical Hertz solution via a Brazilian test. For the homogeneous media, we use the unidirectional flax/bio based epoxy composite, for which macroscopic longitudinal Young’s modulus is experimentally quantified. Furthermore, comparisons with FEM are performed in this context. Results highlight a good agreement between both approaches in terms of stress fields including extremal values.

This paper’s goal was to select methods and a calibration procedure which would lead to the determination of relevant parameters of a discrete element method (DEM) and virtual material creation. Seven particulates were selected with respect to their shape (spherical and non-spherical), size and density. The first calibration experiment involved “packing test” to determine the shape accuracy and bulk density of virtual packed particulates. The series of simulations were compared with real experiments, and the size, shape and density of virtual particles were optimized. Using three apparatuses, the input parameter values were experimentally determined for a contact model that defines the behavior of particulates in DEM simulations. The research part of the paper examines the influence of factors such as particle number; pile formation method; and the method of evaluation of the angle of repose on the process of the calibration of virtual material. The most reproducible results were achieved by the “pilling” method and by the rotating drum—both evaluated by the geometric method. However, it is always advisable to make an overall visual comparison of the slope shape between the calibration simulation and the experimental curves. The bowl’s diameter to particle size ratio should be greater than 25, and the calibration experiment should contain approximately 4000 particles to ensure representative results during angle of repose calibration experiment.

An impact energy erosion model with an energy allocation rule for the discrete element method

Particle sedimentation is a fundamental process encountered in various industrial applications. In this study, we used immersed boundary lattice Boltzmann method and discrete element method (IBM-LBM-DEM) to investigate two-particle sedimentation. A lattice Boltzmann method was used to simulate fluid flow, a discrete element method was used to simulate particle dynamics, and an immersed boundary method was used to handle particle–fluid interactions. Via the IBM-LBM-DEM, the particles collision process in fluid or between rigid walls can be calculated to capture the information of particles and the flow field more efficiently and accurately. The numerical method was verified by simulating settling of a single three-dimensional particle. Then, the effects of Reynolds number (Re), initial distance, and initial angle of particles on two-particle sedimentation were characterized. A specific focus was to reproduce, analyze, and define the well-known phenomenon of drafting-kissing-tumbling (DKT) interaction between two particles. Further kinematic analysis to define DKT is meaningful for two-particle sedimentation studies at different particle locations. Whether a pair of particles has experienced DKT can be viewed from time plots of the distance between the particles (for kissing), the second-order derivative of distance to time (for drafting), and angular velocities of particles (for tumbling). Simulation results show that DKT’s signatures, including attraction, (near) contact, rotation, and in the end, separation, is only completely demonstrated when particles have nearly vertically aligned initial positions. Hence, not all initial positions of particles and Reynolds numbers lead to DKT and not all particle–particle hydrodynamic interactions are DKT. Whether particle–particle interaction is attractive or repulsive depends on the relative positions of particles and Re. Collision occurs when Re is high and the initial angle is small (<20°), almost independent of the initial distance.