DEM Simulations Research Articles

Impact of Non-Linear Contact Models on DEM Simulations of Cemented Angular Limestone Behavior

Impact of Non-Linear Contact Models on DEM Simulations of Cemented Angular Limestone Behavior

GPU based discrete element modeling for convex polyhedral shape particles: Development and validation

Particle dynamics simulations face a significant challenge in understanding the intricate behaviors of convex polyhedral particles due to their complex geometries and interactions. DEM emerges as a key method, illuminating the concealed intricacies of these geometric entities. Traditional algorithms often require cumbersome processes to check each type of contact individually. However, Gilbert-Johnson-Keerthi’s (GJK) and the expanding polytope algorithm (EPA) provide efficient numerical solutions for polyhedral contact detection and contact resolution. These Minkowski difference-based methods streamline contact detection and overlap computation, paving the way for deeper exploration of three-dimensional contact theory within DEM simulations. By leveraging GPU computational power, this paper outlines key algorithmic steps and verifies the code’s accuracy through comparison with simulated and experimental data, with an average deviation of less than 5%. This study explores the impact of particle shape on the dynamics and mechanical behavior of densely packed systems, particularly in hoppers and tumblers. Spherical particles discharge faster but mix more slowly than polyhedral shapes, with icosahedrons achieving quicker full mixing. These results align with experimental findings, further validating the simulation approach.

Influence of steam curing on cyclic triaxial characteristics of recycled aggregate concrete: Experimental analysis and DEM simulation

The engineering application of steam-cured recycled aggregate concrete (often in a cyclic triaxial stress state) can not only improve the recycling efficiency of resources, but also accelerate the construction progress. In this paper, we adopt a combination of laboratory experiments, theoretical analysis and numerical simulation to study the cyclic triaxial characteristics of recycled aggregate concrete (RAC) under different curing regimes (20℃, 40℃, 60℃ and 80℃). The results indicate that as the steam curing temperature increases, the internal damage caused by high temperature continues to intensify, and the cyclic triaxial failure mode of RAC transitions from shear failure to compression failure, and its peak strength, dynamic elastic modulus, and dilatancy angle all show a downward trend. A linear prediction model is established based on the strong correlation between peak strength and steam curing temperature. As the loading cycles increase, the dynamic elastic modulus and dilatancy angle of RAC show exponential and linear downward trends respectively, and the decline rate increases with the increase of steam curing temperature, and the prediction models for dynamic elastic modulus and dilatancy angle are established based on quantitative relationships between variables. On the basis of experimental analysis results, a cyclic triaxial DEM model considering real recycled aggregates is established by introducing steam-cured damage indexes into the mesoscopic parameters, and its applicability in predicting cyclic triaxial mechanical properties of RAC under different curing regimes is verified. The research outcomes can save a lot of test cost and time consumption for developing steam-cured concrete with better performance, and have important theoretical guidance and practical significance for the wide application of steam-cured concrete.

DEM simulations of particle dissolution effects on the passive earth pressure of retaining walls

Assessing passive earth pressure is fundamental in geotechnical engineering practice. Mineral dissolution in the soil can reduce the soil strength, causing an overestimation of passive earth pressure in design. In this study, the effect of dissolution on passive earth pressure on retaining walls is investigated by using discrete element method, taking into account three modes of motion: translation (mode T), rotation around the wall bottom (mode RB), and rotation around the wall top (mode RT). The simulations show that the particle dissolution results in a significant reduction in the passive earth pressure on retaining walls. The largest reductions of resultant force are about 73.5 % and 78.5 % in modes T and RT when wall displacement is minor. Nevertheless, when the wall displacement is large, the largest reduction in resultant force is approximately 74.7 % in mode RB. A detailed analysis is provided to explain this phenomenon. Dissolution also leads to weak force chains and an increase in soil porosity, thereby weakening wall-soil interactions. Dissolution has a significant effect on passive earth pressure of lower part's soil. This study suggests that reinforcing the lower part's soil and preventing seepage in this area can help mitigate the effects of dissolution on retaining walls.

The effect of particle size distribution on the collapse of wet polydisperse granular materials

Low-saturation liquid-containing granular materials are commonly encountered in both natural and industrial settings, where interstitial liquids significantly affect the motion of particles, while particle size polydispersity plays a crucial role in determining the level of system cohesion. In this study, the collapse of wet polydisperse granular columns is numerically investigated based on the developed discrete element model, with corresponding dam-break experiments performed to validate our numerical model and methodology. The dependence of the dynamics and flow mobility on particle size distribution is primarily examined, and the underlying mechanisms are also explored by analyzing particle path lengths and average fidelity. Building upon the effective Bond number proposed using the mixing theory, a macroscopic cohesion parameter at the material scale is defined by considering the dependence of the collapse on the system size effect. The relevance of this cohesion parameter in describing different wet polydisperse granular collapses is further validated based on our designed experimental tests and DEM simulations. The approach of constructing the cohesion parameters at different scales can be extended to characterize cohesion effects in more complex wet polydisperse granular flows and describe their associated rheological behaviors.

Performance of Monotonic Pile Penetration in Sand: Model Test and DEM Simulation

By integrating laboratory tests and three-dimensional discrete element methods, this research extensively explores the macroscopic and microscopic mechanisms of static pile penetration in standard sand. Initially, the mesoscopic parameters of standard sand were established via flexible triaxial compression tests, and a three-dimensional discrete element model was created using the particle size magnification technique. The study results confirm the rationality of parameter selection and numerical modeling by comparing penetration resistance and displacement obtained from laboratory model tests and discrete element simulations. Initially, penetration resistance swiftly increases, then stabilizes progressively with increasing depth. The lateral friction resistance grows with penetration depth, especially peaking near the cone tip. Moreover, horizontal stress quickly rises during pile penetration, mainly caused by the pile foundation compressing the adjacent soil particles. Displacement of the foundation particles is primarily focused around the pile side and cone tip, affecting an area roughly twice the pile diameter. Soil particle displacement exhibits a pronounced vertical downward movement, primarily driven by lateral friction. The distribution of force chains among foundation particles indicates that the primary stressed areas are at the pile ends, highlighting stress concentration features. This research offers significant insights into the mechanical behaviors and soil responses during pile foundation penetration.

Shape effects in binary mixtures of PA12 powder in additive manufacturing

The quality of powder spread in additive manufacturing devices depends sensitively on particle shapes. Here, we study powder spreading for mixtures of spherical and irregularly shaped particles in Polyamide 12 powders. Using DEM simulations, including heat transfer, we find that spherical particles exhibit better flowability. Thus, the particles are deposited far ahead of the spreading blade. In contrast, a large fraction of non-spherical particles hinders the flow. Therefore, the cold particles are deposited near the front of the spreading blade. This results in a temperature drop of the deposited particles near the substrate, which cannot be seen with spherical particles. The particles of both shapes are homogeneously distributed in the deposited powder layer.

Study of triaxial loading of segregated granular assemblies through experiments and DEM simulations

Study of triaxial loading of segregated granular assemblies through experiments and DEM simulations

DEM simulation of subsoiling in tropical sugarcane fields: Effects of opposing subsoiler design and model parameters

During the subsoiling operation of clayey soil in sugarcane fields in tropical areas, there are key limitations such as high resistance caused by high adhesion and unsatisfactory subsoiling effects of machines and tools. To address these issues, this study first explored the machine's structural characteristics using a mechanism analysis method to identify key optimization parameters. Then, simulation tests were conducted based on the root-soil complex model. Single factor analysis was conducted to evaluate the effects of backsweep angle, plow spacing, and penetration angle on tillage resistance and soil disturbance, determining their optimal ranges. Finally, a multi-factor orthogonal rotation combination experiment was designed using Box-Behnken theory to optimize the subsoiler's structural parameters. The regression analysis results show that the optimized model significantly improves the reliability of simulation results. Field verification test results show that under the conditions of a sweep angle of 50°, a plow spacing of 60 cm, and a soil penetration angle of 45°, the average tillage resistance of subsoiling operations is significantly reduced to 76.2 kN, and the amount of soil disturbance is reduced to 1780cm². The optimization reduced tillage resistance and soil disturbance by 34.42% and 30.50% respectively, significantly improving the subsoiling performance. To our knowledge, this is the first work applying discrete element methods to subsoiling operations in red clay soils in sugarcane growing areas. The proposed machine structure has higher efficiency and better performance, providing an effective reference for the design and application of sugarcane field subsoiling machinery.

An equivalent model of horizontal vibratory finishing process: Model construction and analysis based on similarity theory

The manufacturing cost of parts and experimental equipment seriously limits the development of experimental research on mass finishing. To tackle this issue, an equivalent model construction method was proposed based on similarity theory. First, the similarity criterion of relevant physical quantities in horizontal vibratory finishing was determined, and the calculation formula of the distortion coefficient was derived to correct the prediction results. Second, the effects of vibration parameters on the characteristics and similarities of particle velocity and normal force were analyzed by DEM simulation. Finally, the validity of the equivalent model was proved by PIV and force tests. The results show that the equivalent model can be used to reflect the variation of the actual model and anticipate the particle velocity and normal force. And the prediction accuracy can reach 99.56 % and 97.46 %, respectively. The research results provide a new efficient and low-cost method for researching the mass finishing process.

The influencing factors and mechanisms of granular flow dynamics

This study proposes an approach that combines experimental and numerical simulations to comprehensively elucidate the factors influencing the physical behaviour of granular flows. The rotating drum experiments are used to investigate the behaviour of mono-sized and poly-sized particle systems under increasing rotational speed of the drum. The experimental data is then integrated with GPU-based DEM simulations. This allows for the inverse calibration of key parameters like dynamic friction and damping ratios. These calibrated parameters are subsequently utilized throughout the study to explore the influence of various factors on granular flow dynamics. The increase in RPM leads to a rise in particle dispersion due to the interplay between centrifugal forces and particle collisions. The higher dynamic friction angles result in a steeper angle of repose and consequently shorter and longer runout distances for both mono-sized and poly-sized particles respectively. Conversely, increasing damping ratios lead to a decrease in the angle of repose and promote a well-organised flow patterns with increased frictional resistance, suggesting a significant impact on overall flow dynamics. The study paves the way for new insights into the behaviour of complex particulate systems, potentially benefiting various fields concerned with granular flow phenomena.

The influence of steep faults at various avoidance distances on the stability of hard rock caverns: Physical model experiments and DEM simulations

The influence of avoidance distances from steep faults on the stability of hard rock caverns is pivotal in deep engineering, yet research in this area remains scant. This study constructed 10 physical models of granite caverns, incorporating artificial cement-filled faults at various avoidance distances ranging from 0 to 2.4 cavern diameters. Biaxial compression experiments were conducted, utilizing acoustic emission and digital image correlation techniques for monitoring. Coupled with discrete element method simulations, the analysis explored the influence mechanisms of fault-cavern distance on the failure modes of hard rock caverns. The results reveal a positive correlation between the peak strength of specimens and fault-cavern distance. As the distance approaches 0.6 cavern diameters, the fault significantly influences the displacement field, crack number, distribution, and crack initiation location within the specimens. Furthermore, even at 1.2 cavern diameters, the fault still exhibits a notable influence on the stress field, external load transfer paths, cavern failure modes, and the size distribution of spalling debris. These findings can guide the selection of sites for deep hard rock caverns, ensuring a minimum fault-cavern distance of 1.2 cavern diameters to mitigate the risks posed by steep faults to such caverns.

DEM simulation on the AE characteristics and mechanical behaviour of soft-hard composite rock-like materials with parallel fissures subjected to uniaxial compression

DEM simulation on the AE characteristics and mechanical behaviour of soft-hard composite rock-like materials with parallel fissures subjected to uniaxial compression

Predicting cyclic liquefaction behavior of saturated granular materials using an updated state evolution model

Liquefaction and dynamic response of granular materials under dynamic loading has been studied intensively in field and laboratory tests. However, theoretical modeling and analytical solutions on liquefaction are still lagging and investigations are mostly restricted to laboratory observations. To investigate undrained liquefaction shear deformation and fluidity of granular material, the updated state evolution model is proposed by introducing an excess pore water pressure ratio parameter. A series of undrained cyclic triaxial tests and DEM simulations are conducted to verify the proposed model. The result indicates that the liquefaction behavior of granular materials can be captured by the updated state evolution model both at constant and varying loading frequency. Furthermore, the state parameter based on the deviatoric strain and excess pore water pressure ratio is determined to quantify assess the fluidity of granular materials. It facilitates the refinement of the discriminative criteria for cyclic liquefaction of granular materials. This parameter increases slowly at the beginning of loading, followed by a rapid and fluctuating rise, and reaches the peak before the initial liquefaction. Another significant finding is that the turning point of the state parameter range from 0.89 to 0.95 in the θ−t/t0 plane and between 0.84 and 0.94 in the θ−ru plane, as affected by the cyclic loading conditions.

Optimizing structural anisotropy and thermal properties of hexagonal boron nitride ceramics through bimodal particle size distribution

Hexagonal boron nitride (h-BN) ceramics exhibit exceptional thermal conductivity, yet integrating their anisotropic thermal properties into bulk ceramics remains challenging. This study investigates the impact of bimodal particle size distribution on the structure, thermal and mechanical properties of h-BN ceramics. Through DEM simulations and hot-pressing techniques, we demonstrate that incorporating large particles into fine ones significantly enhances the packing density, structural anisotropy and thermal properties of h-BN ceramics. The resulting high-purity h-BN ceramic prepared with 10 vol% large particles (10LBN) exhibits an Index of Orientation Preference (IOP) of −2780.57, surpassing that of ceramics without large particles (0LBN) at −1168.61. By regulating the structure, 10LBN h-BN ceramic show a notable increase of the in-plane thermal conductivity from 81.78 for 0LBN ceramics to 153.20 W/(m·K), along with a flexural strength of 58.21 MPa. Additionally, structural characterization and performance testing show that, compared to adding sintering additives, incorporating large plate-like h-BN particles offers greater benefits in optimizing the orientation structure and thermal conductivity of h-BN ceramics. These findings offer new insights for powder compaction using dual-sized plate-like particles and into the sintering of h-BN ceramics with tailored structural and thermal properties.

Research on flow field characteristics of curved pipe in bulk grain cyclone conveying based on gas solid coupling

AbstractIn the transportation of grain particles, elbow wear is a pressing concern. A novel approach to mitigate this involves introducing a lateral air supplement rotation mechanism that enables swirling transportation. We assessed the integration of this device at angles of 45°, 55°, and 65° to the central axis of the main pipeline. While the multiphase flow velocity within the main conduit was held constant at 20 m/s, the lateral device's airflow velocities were tested at 20, 30, and 40 m/s. Utilizing the CFD–DEM simulation for grain particle transportation around bends, our findings indicated that an optimal arrangement was with the lateral device at 55° relative to the main pipeline, with an airflow speed of 30 m/s. This setup fostered a forward particle spiral, drastically diminishing wear on the pipe wall. Experimental evaluations corroborated the simulation's outcomes.

A sub-grid gas–solid interaction model for coarse-grained CFD–DEM simulations

Computational fluid dynamics has been coupled with coarse-grained discrete element method (CFD–CGDEM) to study industrial-scale gas–solid fluidized beds. However, solving the fluid flow at a grid size of several coarse-grained particles (CGPs) cannot accurately predict the gas–solid interaction, and thus reasonable correction model is required. This study develops a sub-grid model to correct the gas–solid interaction and thereby improve the accuracy of coarse-grid CFD–CGDEM simulation. First, the fluid grids with the size of several CGPs are divided into the auxiliary three-dimensional sub-grids. Then, the CGPs are mapped onto the sub-grids by the weight function to reconstruct the local solid velocities and solid fractions of the sub-grids. After that, the simplified mass and momentum conservation equations are solved to resolve the local gas velocities of sub-grids. Finally, the drag on each CGP is corrected based on the local gas-phase hydrodynamics of sub-grids. Results show that the sub-grid model resolving the local gas velocity in all the three directions works better than that just in the main stream direction. The coarse-grid CFD–CGDEM simulations with this sub-grid model agree well with the experiment for the bubbling, turbulent and rapid fluidized beds. Besides, by solving the simplified conservation equations with graphics processing units, the coarse-grid CFD–CGDEM coupled with the sub-grid model achieves 674 times speedup compared to the fine-grid simulation, demonstrating a much more efficient method to study industrial-scale gas–solid fluidized beds.

Origin of granular axial segregation bands in a rotating tumbler: An interface-mixing driven Rayleigh-Taylor instability

The origin of large and small particle axial bands in long rotating tumblers is a long-standing question. Using DEM simulations, we show that this axial segregation is due to a Rayleigh-Taylor instability which is characterized by the fact that the density of a granular medium increases with mixing and decreases with segregation. For initially mixed particles, segregation and collisional diffusion in the flowing layer balance and lead to a three-layer system, with a layer of large particles over a layer of small particles, and, interposed between these layers, a layer of more densely packed mixed particles. The higher density mixed particle layer over the lower density small particle layer induces a Rayleigh-Taylor instability, evident as waviness in the interface between the layers. The waviness destabilizes into ascending plumes of small particles and descending plumes of mixed particles with large particles enriched near the surface, which become evident as small and large particle bands visible at the free surface. Rolls driven by segregation at the tilted interface between plumes maintain the pattern of frozen plumes. Published by the American Physical Society 2024

Simulated Slidequakes: Insights From DEM Simulations Into the High‐Frequency Seismic Signal Generated by Geophysical Granular Flows

AbstractGeophysical granular flows generate seismic signals known as “slidequakes” or “landquakes”, with low‐frequency components whose generation by mean forces is widely used to infer hazard‐relevant flow properties. Many more such properties could be inferred by understanding the fluctuating forces that generate slidequakes' higher frequency components and, to do so, we conducted discrete‐element simulations that examined the fluctuating forces exerted by steady, downslope‐periodic granular flows on fixed, rough bases. Unlike our previous laboratory experiments, our simulations precluded basal slip. We show that, in its absence, simulated basal forces' power spectra have high‐frequency components more accurately predicted using mean shear rates than using depth‐averaged flow velocities, and can have intermediate‐frequency components which we relate to chains of prolonged interparticle contacts. We develop a “minimal model”, which uses a flow's collisional properties to even more accurately predict the high‐frequency components, and empirically parametrize this model in terms of mean flow properties, for practical application. Finally, we demonstrate that the bulk inertial number determines not only the magnitude ratio of rapidly fluctuating and mean forces on a unit basal area, consistent with previous experimental results, but also the relative magnitudes of the high and intermediate‐frequency force components.

Development of Models Relating Screw Conveying Capacity of Concrete to Operating Parameters and Their Use in Conveyor Operating Strategies to Consider Batch Production

The screw conveyor is the key equipment used to realize the casting and forming of concrete in prefabricated components (PC), and its performance affects the PC shape, quality, and cost. In batch production, there is a process variable, the residence time. It is affected by the quality of the downstream vibration process. This also results in operating parameters that are difficult to match to the time scales. Eventually, it can lead to problems such as low casting efficiency or poor molded quality. In this paper, the DEM simulation method is used to explain and quantify the relationship between the screw conveying capacity and three important operating parameters: the screw’s outer diameter, residence time, and screw speed. The axial and radial velocity vectors are used as features to analyze the changing rule of particle motion trajectory and mass flow rate. Based on the simulation data, the operating parameters and the mass flow rate are forward-fitted to establish the prediction model of the screw conveying capacity. In addition, the residence time is backward fitted from the screw speed and mass flow rate. It is used to estimate the concrete workability. Furthermore, the fitted forward and backward models explore how to propose feasible operational strategies to achieve automatic discharge during batch production.