Simulation Of Particulate Flows Research Articles

Study on the surface texture-microparticles synergistic lubrication mechanism of artificial joint

Abstract After total joint replacement surgery, frequent injection of lubricating substances is required due to lubrication failure caused by the absorption of lubricants by the human body, which will cause serious physiological and psychological burdens on patients. In this study, with the help of surface texture aggregation of particles, gelatin microgel particles with a diameter of about 4 µm were used as lubricating substances to achieve aggregation in the texture, thus improving the friction environment under low-concentration lubricants. The study found that the rectangular cross-section texture demonstrated the most effective anti-friction properties, resulting in a 26% reduction in friction coefficient compared to non-textured surfaces. Additionally, utilizing COMSOL for particle flow simulation, researchers observed the motion behavior of particles within the texture and clarified the mechanism of particle aggregation and the improvement of lubrication on the surface. This article confirms the beneficial effects of combining surface texture and particles on lubrication, thus providing valuable insight for improving lubrication in dilute particle solutions.

Particle flow simulation of crack propagation and penetration in Brazilian disc under uniaxial compression

The presence of rock mass fractures has always been a subject of study for the prevention and control of related natural disasters. To understand the effects of different dip angles and horizontal distances on crack development, numerical simulation experiments on Brazilian disks under uniaxial compression were con-ducted using the PFC2D particle flow program. A function module was utilized to monitor the expansion and quantity of cracks. The numerical simulation results under 0° conditions were in good agreement with the experimental results, validating the rationality of the numerical simulation. The simulation results indicate that: under single fracture conditions with different dip angles, samples with angles between 30° and 60° produced typical wing-shaped cracks. At 0°, cracks propagat-ed through the center of the fracture, while at 90°, cracks initiated from the tip of the fracture and propagated through the sample. The peak stress and the number of cracks in the samples first decreased and then increased with the increase of the dip angle, reaching a maximum at 90°. For samples with double fractures and varying horizontal distances, all produced wing-shaped cracks. Their peak stress and the number of cracks increased monotonically with the increase in distance, reaching a maximum at a distance of 30mm. The experimental results confirmed that the PFC2D program can effectively simulate the process of crack initiation and development, and the research findings provide a reference for correctly understanding the fracture mechanics of fractured rock masses.

A Review of the Homogenized Lattice Boltzmann Method for Particulate Flow Simulations: From Fundamentals to Applications

The homogenized lattice Boltzmann method (HLBM) has emerged as a flexible computational framework for studying particulate flows, providing a monolithic approach to modeling pure fluid flows and flows through porous media, including moving solid and porous particles, within a unified framework. This paper presents a thorough review of HLBM, elucidating its underlying principles and highlighting its diverse applications to particle-laden flows in various fields as reported in literature. These include studies leading to new fundamental knowledge on the settling of single arbitrarily shaped particles as well as application-oriented research on wall-flow filters, hindered settling, and evaluation of the damage potential during particle transport. Among the strengths of HLBM are its monolithic approach, which allows seamless simulation of different fluid-solid interactions, and its ability to handle arbitrary particle shapes, including irregular and concave geometries, while resolving surface interactions to capture local forces. In addition, its parallel scheme based on the lattice Boltzmann method (LBM) results in high computational efficiency, making it suitable for large-scale simulations, even though LBM requires small time steps. Important future development needs are identified, including the addition of a lubrication force correction model, performance enhancements, such as support for hybrid parallelization and GPU, and the extension of compatible contact models to accommodate concave shapes. These advances promise expanded capabilities for HLBM and broader applicability for solving complex real-world problems.

The One-Dimensional Flow Pressure Loss Correction Model Based on the Particle Flow through Concrete Bend

The rheological behavior of concrete pumping in the bend section is complex. The existing conversion calculation for pumping pressure loss based on the straight pipe section often fails to meet engineering precision requirements. This paper develops a pressure loss correction model for a concrete pumping bend section based on a new one-dimensional flow model for the straight pipe section simulation of particle flow; the study analyzes the impact of parameters such as elbow shape and pumping flow rate on pressure loss. An equivalent conversion relationship between the bend section and the straight section is established. The comparison with the measured pressure values from the project shows that the relative error between the pressure loss values calculated by the correction model and the actual measurements is 15.8%.

Research on the Water Inrush Mechanism and Grouting Reinforcement of a Weathered Trough in a Submarine Tunnel

Based on the structural and geological characteristics of the F1 weathering trough of a submarine tunnel and its spatial relationship with the cavern, a simplified calculation model of the weathering trough water inrush was established, and the formation, development process and influencing factors of the water inrush channel in the water-resistant rock layer were carried out by a numerical simulation of particle flow. It shows that the integrity and stability of the critical water-resistant rock mass is the key to preventing water inrush, and the identification and positioning of the water inrush channel is the basis for the grouting reinforcement design of the weathering groove of the submarine tunnel. Based on above research results, the F1 weathering trough was blocked and reinforced by the composite grouting method, and the engineering reinforcement effect was good.

GPU-powered CFD-DEM framework for modelling large-scale gas–solid reacting flows (GPU- rCFD-DEM) and an industry application

The coupling of CFD (computational fluid dynamics) and DEM (discrete element method) is extensively used for simulating gas–solid reacting flows in various industrial processes, while its high computational cost limits its industry applications, especially large-scale systems with a large number of particles and complex geometries. This paper reports a robust highly efficient GPU (graphics processing unit) − powered CFD-DEM coupling approach that is, for the first time, capable of simulating large-scale gas–solid reacting flow systems with complex geometries (GPU- rCFD-DEM). The fluid flow calculations are performed using CPU parallelization, while the particle flow simulations leverage GPU parallelization, and the coupling calculations between CFD and DEM are fully implemented on GPU. The model includes advanced coupling strategies to enhance numerical stability, especially when handling complex geometries with unstructured CFD meshes. The developed model is validated through experimental measurements and its computational performance is evaluated by comparison with previous GPU-based simulations. It shows good agreement with the experiments and superior performance compared to the traditional coupling method. The model is then applied to simulate raceway dynamics and coke combustion in an industrial-scale blast furnace, showcasing its effectiveness in handling complex geometries and a huge number of particles in gas–solid reacting flows. This work provides an efficient and robust solution for numerically simulating industrial applications of gas–solid reacting flow systems.

Investigation of macro-micro mechanical behaviors and failure mechanisms in cemented tailings backfill with varying proportions of fine

This study investigates the impact of tailings characteristics, particularly fine tailings, on the efficiency and effectiveness of mining engineering backfill operations. Contrary to the common underestimation, fine tailings significantly affect the mechanical properties and failure mechanisms of Cemented Tailings Backfill (CTB). Through a series of experiments examining various particle size distributions and employing advanced Scanning Electron Microscopy (SEM) and Mercury Intrusion Porosimetry (MIP) techniques, this research clarifies the role of fine tailings in modifying the mechanical strength and failure patterns of CTB. Adding 20 % fine tailings optimally enhances both the early and long-term strength of backfill materials. With aging, adding 30 % fine tailings shows a rapid strength increase from 7d to 28d, with the 28d strength slightly lower than that of specimens with 20 % addition, but the specimens demonstrate better ductility after reaching peak strength. As the content of fine tailings increases, crack propagation in the specimens shifts from radial to axial, simplifying the fracture mode from complex to singular. Consequently, the fracture mechanism changes from ductile to brittle, and then reverts to ductile. Moreover, using particle flow simulation and moment tensor analysis to study the failure processes indicates that while increased fine particle content boosts the material's initial strength, it ultimately leads to a simpler crack network and volumetric expansion as the primary failure mode. In terms of energy conversion mechanisms, as the content of fine particles increases, the energy first increases and then decreases. During the elastic stage and the early stage of the plastic stage, the energy is mainly converted into strain energy and PB bonding energy, which affects the elasticity of the sample. In the later stage of the plastic phase, the energy is mainly converted into dissipated energy, leading to the failure of the sample.

Fast and precise approximation of Minkowski sum of two rotational ellipsoids with a superellipsoid

In this paper, we propose a fast and accurate method for approximating the Minkowski sum of two rotational ellipsoids with a superellipsoid. The Minkowski sum is used in a variety of applications such as robot motion planning and particle flow simulation requiring collision detection. Many of them are computed based on Minkowski sum, whose accuracy and processing time depend on how the mesh is created. We approximate Minkowski sum with a superellipsoid function. The superellipsoid has various shapes with the value of the exponents, and the computation of the parameters including the exponents is an algebraic computation with the time complexity O(1)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\varvec{O}(1)$$\\end{document}. As a result, approximation error is small and computation is fast in all combinations of rotational ellipsoids. In addition, collision is quickly detected by solving the inequality of superellipsoid.

Particle Flow Simulation Study of Damage Evolution in Expansive Slurry–Fractured Rock Mass Composites under Direct Shear Conditions

Particle Flow Simulation Study of Damage Evolution in Expansive Slurry–Fractured Rock Mass Composites under Direct Shear Conditions

Investigation of the Shear Mechanism at Sand-Concrete Interface under the Influence of the Concave Groove Angle of the Contact Surface

A “random-type” sand–concrete interface shear test was developed based on the sand cone method, with a focus on the most commonly encountered triangular contact surface morphology. A “regular-type” triangular interface, matched in roughness to the “random-type”, was meticulously designed. This “regular-type” interface features five distinct triangular groove inclinations: 18°, 33°, 50°, 70°, and 90°. A series of sand–concrete interface direct shear tests were conducted under consistent compaction conditions to investigate the impact of varying compaction densities and triangular groove inclinations on the shear strength at the interface. Particle flow simulations were utilized to examine the morphology of the shear band and the characteristics of particle migration influenced by the triangular contact surface. This analysis is aimed at elucidating the influence of the inclination of the triangular groove on the shear failure mechanism at the sand–concrete interface. The findings indicate that: (1) The morphology of the interface significantly impacts the shear strength of the sand–concrete interface, while the shape of the stress-displacement curve experiences minimal alteration. (2) At smaller inclination angles, particle contact forces are arranged in a wave-like configuration around the sawtooth tip, resulting in a non-uniform stress distribution along the sawtooth slope. However, as the inclination angle grows, the stress concentration at the sawtooth tip diminishes, and the stress distribution across the sawtooth slope becomes more consistent. (3) Particle migration is significantly influenced by the sawtooth’s inclination angle. At lower angles, particles climb the structure’s tip through sliding and rolling. As the angle increases, particle motion shifts to shear, accompanied by a transition in friction from surface friction to internal shear friction. This leads to the formation of a wider shear band and an increase in shear strength.

Experimental study on failure mode and fracture evolution characteristics of red shale in Kaiyang Phosphorus mining area

In order to study the failure mode and fracture evolution characteristics of red shale in Kaiyang Phosphorus mining area, conventional triaxial compression mechanical tests of red shale with different bedding dip angles were carried out by using DSTD-1000 electro-hydraulic servo rock mechanics experiment system. Based on the laboratory test results, the conventional triaxial particle flow simulation of red shale samples with different bedding dip angles was carried out using discrete element PFC2D. The results show that: (1) the failure mode of red shale is controlled by bedrock when the bedding dip angle is 0° and 60° ~ 90°. When the bedding dip angle is 15° ~ 45°, the rock failure mode is controlled by bedding. The compressive strength of rock is the minimum when the bedding dip angle is 30°and the maximum at 0°, which is about 2 times of the minimum. (2) In the failure process of red shale, the cracks with different bedding dip angles show slow growth stage, accelerated growth stage and stable stage with axial strain. The whole failure process is dominated by tensile cracks, accompanied by a few shear cracks. (3) The type of displacement field varies with the bedding dip angle: tensile failure and shear failure are the main displacement field types at 15° ~ 45°, and mixed failure is often the main mode at 60° ~ 90°and 0°. The research results provide the basis and reference for the safety control of red shale roadway.

Migration trajectories and blocking effect of the fine particles in porous media based on particle flow simulation

The particle flow code method based on the discrete element method was used to establish the seepage migration model of fine particles [fine particles (FPs), i.e., suspended particles] in a porous medium. A series of numerical simulations were carried out by changing the particle size, seepage velocity, particle injection number, and wide particle size gradation. The research showed that large FPs play a major role in blocking porous medium channels when the injected FPs have a wide size gradation. Due to the blocking effect, small FPs that would not otherwise have deposited also deposit. Moreover, by increasing the number of large FPs in the mixed particles, the total number of particles deposited and the number of smaller FPs deposited will also increase. The distribution of FPs in porous mediums can be divided into three types: surface deposition, internal deposition, and non-deposition. When the seepage velocity increases and reaches a seepage threshold, which is the critical seepage velocity, the deposited FPs will once again be in a suspended state and undergo migration. On the contrary, the FPs will continue to maintain their sedimentary state, and the critical seepage velocity will also increase correspondingly with increasing particle size.

Lattice Boltzmann method for particulate multiphase flow system

This study proposes a numerical model for particulate three-phase flow in microchannels based on multiphase lattice Boltzmann method (LBM). The model combines the color-gradient method to track the immiscible fluid-fluid interface and the two-fluid model (TFM) to describe particle-particle and particle-fluid interactions, which can efficiently simulate transport and displacement processes involving large amounts of particles. A mixture-rheology TFM algorithm is proposed by introducing a mixture phase with rheology properties obtained from experiments instead of the conventional TFM particle phase with artificial viscosity models. Multi-relaxation-time (MRT) collision operator and GPU computing are adopted to enhance the numerical stability and efficiency. Various theoretical benchmarks for particle transport and two-phase flow are performed respectively to verify the accuracy of the proposed model. Exceptional consistency between results from particulate three-phase flow simulation and microfluidic experiments further confirms the reliability of our model, especially in capturing the inertial lagging and accumulation phenomena under multiphase and porous flow conditions. The proposed numerical framework will benefit our understanding of multiphase displacement with microgels in microchannels with complex geometries.

A novel particle decomposition scheme to improve parallel performance of fully resolved particulate flow simulations

This study addresses the challenge of simulating realistic particle systems by proposing a novel particle decomposition scheme that improves the parallel performance of surface resolved particle simulations. Realistic particle systems often involve large numbers of particles and complex particle shapes. The resulting need to account for shape factors requires the inclusion of even more particles to obtain statistically meaningful results. However, the computational cost increases with the number of particles, making efficient parallelization crucial. Therefore, the proposed scheme aims to improve the scalability by optimizing the communication and data management between processors. Through hindered settling experiments, the applicability and performance of the novel particle decomposition scheme are thoroughly investigated using the homogenized lattice Boltzmann method. The results show that the proposed method significantly improves the performance, especially in scenarios with a large number of particles, by reducing communication constraints and improving scalability. This research contributes to the advancement of computational methods for efficient study of complex particle systems and provides valuable insights for future developments in this field.

Step-path failure mechanism and stability analysis of water-bearing rock slopes based on particle flow simulation

Step-path failure is a typical damage mode of rock slope featuring intermittent joints. In order to investigate the effect of water-fracture interaction on the stepped failure, laboratory experiments of rock materials and numerical simulation of slopes with intermittent joints were conducted. At the laboratory scale, the influence of water-fracture interaction on the mechanical characteristics and acoustic emission response of step-path failure in rocks was examined, and the particle flow model parameters were calibrated for uniaxial testing. At the field scale, the step-path failure process of the rock slope near Longmenqiao reservoir was numerically analyzed using the particle flow method to assess the effects of water and joint distribution on failure. Laboratory findings suggest that compared with dry red sandstone, the average uniaxial compressive strength and average elastic modulus of saturated red sandstone decrease by 35.82% and 23.33%, respectively. The increase in fissure angle in sandstone results in enhanced strength and stiffness, alongside a transition from shear to tension failure mode, with a delayed occurrence of large-scale acoustic emission events. In particular, within the 30°–60° range, the uniaxial compressive strength and elastic modulus of saturated sandstone increase by 66.78% and 29.64%, respectively. The analysis of crack propagation and stress field during slope deformation shows that step-path failure is a progressive process from bottom to top, and revealed the control mechanism of joint and bridge angle on slope stability. The research results enhance the understanding of the step-path failure mechanism and provide important theoretical support for the stability evaluation and landslide monitoring of rock slopes with intermittent joints.

Particle Flow Simulation of Sand Loss through Cracks in Segmental Linings

Particle Flow Simulation of Sand Loss through Cracks in Segmental Linings

Analysis of Mechanical Properties and Crack Failure Characteristics of Rock with Pre-Existing Fissures Using Particle Flow Simulation

Analysis of Mechanical Properties and Crack Failure Characteristics of Rock with Pre-Existing Fissures Using Particle Flow Simulation

Numerical simulation of particle flow around the impulsive monopile

In order to explore the particle flow field characteristics around the impulsive monopile, the particle flow field and particle distribution around impulsive monopile were simulated based on the classical impulsive flow around a monopile at medium and low Reynolds number (Re = 550). The Large Eddy Simulation (LES) technique was harnessed to capture the dynamic flow field around the monopile and to facilitate in-depth analysis of vortex evolution. Concurrently, the Lagrangian approach was employed to track and characterize particle trajectories, shedding light on the patterns of particle diffusion and congregation in the vicinity of the monopile under varying Reynolds numbers. Additionally, this investigation unveils the distinctive behavior of particle flow during the impulsive starting and stopping phases adjacent to the monopile. Leveraging these insights, the study elucidates the intricate influence of the unsteady flow field dynamics around the monopile on the progression of particle flow.

Discrete element simulation of particle flow and separation in a vibrating screen with variable rectangular hole screen

The screening operations aim to achieve large processing capacity and high screening accuracy. However, there is a certain inherent contradiction between these two goals. This paper proposes a variable rectangular hole screen for the negative pressure vibrating screen. Through numerical simulation of variable rectangular hole screens with different length-side ratios, the optimal ratio of adjacent rectangular holes is determined to be 60:37. This innovative screen addresses the limitations of traditional screening operations, effectively improving the processing capacity of drilling fluids. It achieves a high screening accuracy while reducing the adverse effects of screen paste and screen clogging on drilling fluid sieve penetration. The study also investigates the effects of amplitude, vibration frequency, feed rate, and screen transportation speed on particle behavior through the screen. These findings serve as a theoretical basis for optimizing the vibration parameters and working condition parameters of the negative pressure vibrating screen.

Study on the mechanical properties of sandstone-shale composite continental shale gas based on the discrete element method

China has abundant continental shale oil resources, and multiple shale oil exploration and development demonstration zones have been established, making it a key field for oil and gas exploration and development. However, drilling and fracturing treatment often encounter vertically interlaced and superimposed sandstone and shale reservoirs, which significantly affect the stability of the reservoir and the propagation of fractures. Therefore, it is of great significance to study the rock mechanics characteristics of sandstone-shale composite. Using the orthogonal experimental method for microscale parameter calibration, this study established the relationship between the macroscopic and microscopic parameters of sandstone-shale in the parallel bond model. The results of the orthogonal experimental scheme were analyzed using multi-factor variance analysis and linear regression analysis to determine the influence and ranking of each microscale parameter on the macroscopic parameters. Based on the results of uniaxial and triaxial compression experiments, the particle flow simulation of sandstone-shale composite rock samples was studied using the discrete element method. The study investigated the influence of the thickness ratio, combination form, and confining pressure of sandstone-shale on the mechanical properties of composite rock samples. The results showed that: (1) Under uniaxial loading conditions, for the combination form of sandstone on top and shale on the bottom (M-S type), the peak stress of the composite rock sample increased first and then decreased and then increased again with the increase of the thickness of sandstone; for the combination form of shale on top and sandstone on the bottom (S-M type), the peak stress of the composite sample decreased first and then increased and then decreased again with the increase of the thickness of shale. (2) After applying confining pressure, the compressive strength of the M-S type composite sample decreased with the increase of the thickness ratio of sandstone, while the compressive strength of the S-M type composite sample increased with the increase of the thickness ratio of shale as well as the compressive strength of the composite sample was lower than that of any complete rock sample. (3) Under uniaxial loading, the fracture mode of the composite sample was related to the combination form and the thickness ratio of different lithologies. For the M-S type sample, when the thickness ratio of sandstone was <30%, the sample was mainly fractured by tension, and when the thickness ratio of sandstone was >30%, it was mainly fractured by shear mode. For the S-M type sample, when the thickness ratio of shale was <70%, the sample had the most internal shear cracks, and when the thickness ratio of shale was >70%, it had the most tension cracks.