Rock block removal is the prevalent physical mechanism for rock erosion and could affect the stability of dam foundations and spillways.Despite this,understanding of block removal is still inadequate because of the co...Rock block removal is the prevalent physical mechanism for rock erosion and could affect the stability of dam foundations and spillways.Despite this,understanding of block removal is still inadequate because of the complex interactions among block characteristics,hydraulic forces,and erosive processes acting on the block.Herein,based on a previously conducted physical experiment of erosion of a single rock block,the removal processes of two different protruding blocks are represented by a coupled computational fluid dynamics-discrete element model(CFD-DEM)approach under varied flow conditions.Additionally,the blocks could be rotated with respect to the flow direction to consider the effect of the discontinuity orientation on the block removal process.Simulation results visualize the entire block removal process.The simulations reproduce the effects of the discontinuity orientation on the critical flow velocity inducing block incipient motion and the trajectory of the block motion observed in the physical experiments.The numerical results present a similar tendency of the critical velocities at different discontinuity orientations but have slightly lower values.The trajectory of the block in the simulations fits well with the experimental measurements.The relationship between the dimensionless critical shear stress and discontinuity orientation observed from the simulations shows that the effect of block protrusion becomes more dominant on the block incipient motion with the increase of relative protrusion height.To our knowledge,this present study is the first attempt to use the coupled finite volume method(FVM)-DEM approach for modelling the interaction behavior between the block and the flowing water so that the block removal process can be reproduced and analyzed.展开更多
To mitigate and alleviate low wheel-rail adhesion,a train-borne system is utilised to deposit sand particles into the wheel-rail interface via a jet of compressed air in a process called rail-sanding.Britain Rail Safe...To mitigate and alleviate low wheel-rail adhesion,a train-borne system is utilised to deposit sand particles into the wheel-rail interface via a jet of compressed air in a process called rail-sanding.Britain Rail Safety and Standards Board introduced guidelines on the sand particles’shape,size,and uniformity which needs to be adhered to for rail-sanding.To further inves-tigate these guidelines and help improve them,this research presents a parametric study on the particle characteristics that affect the rail-sanding process including density,size and size distribution,coefficient of uniformity,and shape,utilising a coupled computational fluid dynamics-discrete element method(CFD-DEM)model.The efficiency of rail-sanding is esti-mated for each case study and compared to the benchmark to optimise the sand characteristics for rail-sanding.It is concluded that particle size distribution(within the accepted range)has an insignificant effect on the efficiency while increasing particle size or the coefficient of uniformity decreases the efficiency.Particle shape is shown to highly affect the efficiency for flat,compact and elongated particles compared to the spherical shape.The current numerical model is capable of accurately predicting the trends in the efficiency compared to the actual values obtained from full-scale experiments.展开更多
Accurately predicting heat flux in coarse-grained CFD-DEM simulations is a significant challenge.Specifically,the rates of fluid-particle heat exchange,the effective thermal conductivity of a bed of particles,as well ...Accurately predicting heat flux in coarse-grained CFD-DEM simulations is a significant challenge.Specifically,the rates of fluid-particle heat exchange,the effective thermal conductivity of a bed of particles,as well as radiative heat transfer rates are difficult to predict.By using a novel algorithm,we significantly improve the accuracy and stability of such simulations by using a heat exchange limiter.This limiter enables realistic predictions even at time steps that are three orders of magnitude larger than a typical fluid heat relaxation time.Additionally,view-factor-based corrections for radiative heat exchange computations are developed.These corrections ensure an effective thermal bed conductivity with less than 3%error for a coarse-graining ratio of 10.The applicability of the P1 radiation model in coarse-grained settings is also examined,leading to recommendations for the CFD grid resolution to ensure accurate predictions.Our methods significantly enhance stability,accuracy,and computational efficiency,making coarse-grained CFD-DEM simulations more viable for industrial applications.These advancements enable more reliable modeling of high-temperature processes,accelerate optimization studies,and enable virtual equipment design of such processes.展开更多
Computational Fluid Dynamics coupled with Discrete Element Method (CFD-DEM) is a commonly used numerical method to model gas-solid flow in fluidised beds and other multiphase systems. A significant limitation of CFD-D...Computational Fluid Dynamics coupled with Discrete Element Method (CFD-DEM) is a commonly used numerical method to model gas-solid flow in fluidised beds and other multiphase systems. A significant limitation of CFD-DEM is the feasibility of the realistic simulation of large numbers of particles. Coarse-graining (CG) approaches, through which groups of multiple individual particles are represented by single, larger particles, can substantially reduce the total number of particles while maintaining similar system dynamics. As these three CG models have not previously been compared, there remains some debate, however, about the best practice in the application of CG in CFD-DEM simulations. In this paper, we evaluate the performance of three typical CG methods based on simulations of a bubbling fluidised bed. This is achieved through the use of a numerical validation framework, which makes full use of the high-resolution 3D positron emission particle tracking (PEPT) measurements to rigorously validate the outputs of CFD-DEM simulations conducted using various different coarse-graining models, and various different degrees of coarse-graining. The particle flow behaviours in terms of the particle occupancy field, velocity field, circulation time, and bubble size and velocity, are comprehensively analysed. It is shown that the CG simulation starts to fail when the size ratio between the bed chamber and the particles decreases to approximately 20. It is also observed, somewhat surprisingly, that the specific CG approach applied to interparticle contact parameters does not have a substantial effect on the simulation results for the bubbling bed simulations across a wide range of CG factors.展开更多
To fundamentally alleviate the excavation chamber clogging during slurry tunnel boring machine(TBM)advancing in hard rock,large-diameter short screw conveyor was adopted to slurry TBM of Qingdao Jiaozhou Bay Second Un...To fundamentally alleviate the excavation chamber clogging during slurry tunnel boring machine(TBM)advancing in hard rock,large-diameter short screw conveyor was adopted to slurry TBM of Qingdao Jiaozhou Bay Second Undersea Tunnel.To evaluate the discharging performance of short screw conveyor in different cases,the full-scale transient slurry-rock two-phase model for a short screw conveyor actively discharging rocks was established using computational fluid dynamics-discrete element method(CFD-DEM)coupling approach.In the fluid domain of coupling model,the sliding mesh technology was utilized to describe the rotations of the atmospheric composite cutterhead and the short screw conveyor.In the particle domain of coupling model,the dynamic particle factories were established to produce rock particles with the rotation of the cutterhead.And the accuracy and reliability of the CFD-DEM simulation results were validated via the field test and model test.Furthermore,a comprehensive parameter analysis was conducted to examine the effects of TBM operating parameters,the geometric design of screw conveyor and the size of rocks on the discharging performance of short screw conveyor.Accordingly,a reasonable rotational speed of screw conveyor was suggested and applied to Jiaozhou Bay Second Undersea Tunnel project.The findings in this paper could provide valuable references for addressing the excavation chamber clogging during ultra-large-diameter slurry TBM tunneling in hard rock for similar future.展开更多
The rapid discharge of cuttings from the air cushion chamber is crucial for the construction safety and excavation efficiency of the slurry shield.Previous studies have mainly focused on the transport characteristics ...The rapid discharge of cuttings from the air cushion chamber is crucial for the construction safety and excavation efficiency of the slurry shield.Previous studies have mainly focused on the transport characteristics of cuttings in the slurry discharge pipe,while the complete process of cuttings entering the slurry discharge pipe from the air cushion chamber until they are discharged was often overlooked.Based on the CFD-DEM coupling method and combined with actual engineering,this paper established a numerical model that can more completely reflect the cuttings discharge process during slurry shield tunneling,and the effects of the slurry velocity at slurry gate and scouring pipes,inclination angle of slurry discharge pipe,cuttings diameter and shape coefficient were investigated by analyzing the variation in the mass flow rate,mass of discharged cuttings,and discharged ratio.The results revealed that increasing the slurry velocity can promote the discharge of cuttings.To keep the discharged rate at a high level,it is recommended that the slurry velocity at the slurry gate should be greater than 0.15 m/s.Reducing the inclination angle of the discharge pipe is conducive to the rapid discharge of cuttings.Cuttings with large diameter or small shape coefficient are more prone to accumulate in the air cushion chamber and cause clogging risk.The research results not only help to improve engineers'understanding of cuttings discharge in slurry shield,but also provide practical guidance for formulating relevant construction measures.展开更多
A large amount of loose debris materials were deposited on the slope of mountainous areas after the 2008 Ms 8.0 Wenchuan earthquake. During and after the earthquake, these loose debris deposits collapsed and slide int...A large amount of loose debris materials were deposited on the slope of mountainous areas after the 2008 Ms 8.0 Wenchuan earthquake. During and after the earthquake, these loose debris deposits collapsed and slide into valleys or rivers, changing river sediment supply condition and channel morphology. To investigate the mechanisms of granular flow and deposition, the dynamics of slope failure and sediment transportation in typical mountainous rivers of different intersection angles were analyzed with a coupling model of Computational Fluid Dynamics and Discrete Element Method(CFD-DEM). The numerical results show that the change of intersection angle between the granular flow flume and the river channel can affect the deposit geometry and the fluid flow field significantly. As the intersection angle increases, the granular velocity perpendicular to the river channel increases, while the granular velocity parallel to the river channel decreases gradually. Compared to the test of dry granular flow, the CFD-DEM coupling tests show much higher granular velocity and larger volume of sediments entrained in the river. Due to the river flow, particles located at the edge of the deposition will move downstream gradually and the main section of sediments deposition moves from the center to the edge of the river channel. As a result, sediment supply in the downstream river will distribute unevenly. Under the erosion of fluid flow, the proportion of fine particles increases, while the proportion of coarse particles decreases gradually in the sediment deposition. The build-up of accumulated sediment mass will cause a significant increase in water level in the river channel, thus creating serious flooding hazard in mountainous rivers.展开更多
Hydraulic fracturing is a key technology for the development of unconventional hydrocarbon resources.The proppant placement morphology determines the fracture conductivity,thus affecting the reservoir stimulation effe...Hydraulic fracturing is a key technology for the development of unconventional hydrocarbon resources.The proppant placement morphology determines the fracture conductivity,thus affecting the reservoir stimulation effect.In this paper,the proppant migration and placement within complex fractures was studied by considering the fracture wall roughness through computational fluid mechanics-discrete element method(CFD-DEM)in numerical simulation,which is a key approach to study the proppant migration and placement.The results show that the proppant placement non-uniformity,proppant migration capacity,and proppant volume filled in the far-end and the secondary branched fracture are enhanced within the rough fracture compared with those within smooth fractures.The proppant migration capacity is increased within the fracture at low inclination angles(<60°)and low approach angles(<90°),and the proppant placement area is larger in the inclined fracture than that in the vertical fracture.The rise of injection rate and fracturing fluid viscosity causes more proppants migrate to far-end or secondary fractures,resulting in a non-proppant area within the near-wellbore fracture.An increase by 1.3 times in the injection rate and 3 times in the fracturing fluid viscosity leads to a decrease by 26.6%and 27%,respectively,in the proppant placement area within the near-wellbore fracture.The staged injection with small size proppants followed by large size proppants increases the proppant placement area in the primary fracture by 13%-26%,and that with large size proppants followed by small size proppants increases the proppant placement area by 19%-25%,which is due to that the latter method facilitates filling of the secondary branched fracture.The injection location mainly affects the proppant filling degree within the near-wellbore fractures.Compared with the upper injection,the middle and lower injection is not beneficial to filling of proppants within the near-wellbore fracture.展开更多
The solid particle impurities generated by pipe wall corrosion might deposit at the elbow of hilly pipelines during the production shutdown of oil pipelines.These solid particle impurities will seriously affect the sa...The solid particle impurities generated by pipe wall corrosion might deposit at the elbow of hilly pipelines during the production shutdown of oil pipelines.These solid particle impurities will seriously affect the safety of the pipeline operation and the quality of the petroleum products.Thus,it is necessary to study the methods of removing these trapped particles from pipelines.At present,the most common way to remove these solid particle impurities is pigging oil pipelines periodically by utilizing the mechanical pigging method,while the frequent pigging operation will increase the cost and risk of pipeline operation.It is very convenient and economical to remove the accumulated particles out from the pipeline by oil stream,which can be named Hydraulic Pigging Method(HPM).However,the behavior mechanism of particle in flowing oil is still unclear.This motivates the present research on the particles flushed out by the flowing oil.A numerical model governing the trapped particles displacement from the elbow of an inclined oil pipeline is established in the Euler-Lagrangian framework.The simulation is achieved via CFD coupling with DEM.The CFD method is employed to solving the continuous phase flow,while the discrete particle phase is tracked by the DEM.The numerical model is first validated by comparison with results taken from the published literature.From the simulation results,it is observed that the oil stream,carrier phase,can only flush out the solid particles in a certain diameter range under the given operation conditions,and the particles whose diameter beyond that diameter range will cannot be removed out from the pipeline.The influence of the pipe inclined angle,the oil bulk velocity and the particle diameter on the particle migration characteristics is examined in detail.Furthermore,in order to enhance the efficiency of HPM,an Enhanced Hydraulic Pigging Method based on Multi-Physical Field Collaboration(EHPM-MPFC)is also proposed in the present work.The EHPM-MPFC is validated for having high pigging efficiency via the comparison of the migration characteristics of particles during the EHPM-MPFC and HPM process.The present results can provide the guidance to the HPM operation of products pipelines.展开更多
Soil erosion is a critical process that is being studied in soil science, hydraulic engineering, and geotech- nical engineering. Among many societal and environmental impacts, soil erosion is a major cause for the fai...Soil erosion is a critical process that is being studied in soil science, hydraulic engineering, and geotech- nical engineering. Among many societal and environmental impacts, soil erosion is a major cause for the failures of bridges. The erodibility of soil is determined by its physical and geochemical properties and is also affected by surrounding biological activities. In most of the current models for soil erosion, erodibility of non-cohesive soil is characterized by its median grain size （Dso）, density, and porosity. The contribution to erodibility of the irregular shape of soil grains, which plays an important role in the mechanical and hydraulic properties of coarse-grained soils, is generally ignored. In this paper, a coupled computational fluid dynamics and discrete element method model is developed to analyze the influence of the shape of sand grain on soil erodibility. A numerical model for the drag force on spherical and non-spherical particles is verified by using the results from physical free settling experiments. Erosion of sand grains of different shapes is simulated in a virtual erosion function apparatus, a laboratory device used to mea- sure soil erodibility. The simulation results indicate that the grain shape has major effects on erodibility. Spherical particles do not show a critical velocity because of their low rolling resistance, but a critical velocity does exist for angular particles owing to grain interlocking. The erosion rate is proportional to the flow velocity for both spherical and non-spherical particles. The simulation result for angular particle erosion is fairly consistent with the experimental observations, implying that grain shape is an important factor affecting the erodibility of non-cohesive soils.展开更多
An understanding of the particle transport characteristics in a branched network helps to predict the particle distribution and prevent undesired plugging in various engineering systems.Quantitative analysis of partic...An understanding of the particle transport characteristics in a branched network helps to predict the particle distribution and prevent undesired plugging in various engineering systems.Quantitative analysis of particle flow characteristics is challenging in that experiments are expensive and particle flow is difficult to detect without disturbing the flow.To overcome this difficulty,man-made fractal tree-like branched networks were built,and a coupled computational fluid dynamic and discrete element method model was applied.A series of numerical simulations was carried out to analyze the influence of fractal structure parameters of networks on the particle flow characteristics.The joint influence of inertial,shunt capacity and superposition from upstream branches on particle flow was investigated.The injection position at the inlet determined the particle velocity and its future flow path.The particle density ratio,particle size and bifurcation angle had a greater influence on the shunting of K2 branches than that in the K1 level and N_(k22)/N_(k21) reached a maximum at 60°.Compared with a network with an even number of branches,there was a preferential branch when the branch number was odd.The preferential branch effect or asymmetry degree of the level(K2)branches had a more significant impact on particle shunting than that from the upstream branches(K1).展开更多
Fluid-particle interaction underpins important behavior of granular media. Particle-scale simulation may help to provide key microscopic information governing the interaction and offer better understanding of granular...Fluid-particle interaction underpins important behavior of granular media. Particle-scale simulation may help to provide key microscopic information governing the interaction and offer better understanding of granular media as a whole. This paper presents a coupled computational fluid dynamics and discrete element method （CFD-DEM） approach for this purpose. The granular particle system is modeled by DEM, while the fluid flow is simulated by solving the locally averaged Navier-Stokes equation with CFD. The coupling is considered by exchanging such interaction forces as drag force and buoyancy force between the DEM and CFD. The approach is benchmarked by two classic geomechanics problems for which analytical solutions are available, and is further applied to the prediction of sand heap formation in water through hopper flow. It is demonstrated that the key characteristic of granular materials interacting with pore water can be successfully captured by the proposed method.展开更多
A method to generate fractures with rough surfaces was proposed according to the fractal interpolation theory.Considering the particle-particle,particle-wall and particle-fluid interactions,a proppant-fracturing fluid...A method to generate fractures with rough surfaces was proposed according to the fractal interpolation theory.Considering the particle-particle,particle-wall and particle-fluid interactions,a proppant-fracturing fluid two-phase flow model based on computational fluid dynamics(CFD)-discrete element method(DEM)coupling was established.The simulation results were verified with relevant experimental data.It was proved that the model can match transport and accumulation of proppants in rough fractures well.Several cases of numerical simulations were carried out.Compared with proppant transport in smooth flat fractures,bulge on the rough fracture wall affects transport and settlement of proppants significantly in proppant transportation in rough fractures.The higher the roughness of fracture,the faster the settlement of proppant particles near the fracture inlet,the shorter the horizontal transport distance,and the more likely to accumulate near the fracture inlet to form a sand plugging in a short time.Fracture wall roughness could control the migration path of fracturing fluid to a certain degree and change the path of proppant filling in the fracture.On the one hand,the rough wall bulge raises the proppant transport path and the proppants flow out of the fracture,reducing the proppant sweep area.On the other hand,the sand-carrying fluid is prone to change flow direction near the contact point of bulge,thus expanding the proppant sweep area.展开更多
We investigate the effect of particle shape on the transportation mechanism in well-drilling using a three-dimensional model that couples computational fluid dynamics （CFD） with the discrete element method （DEM）. ...We investigate the effect of particle shape on the transportation mechanism in well-drilling using a three-dimensional model that couples computational fluid dynamics （CFD） with the discrete element method （DEM）. This numerical method allows us to incorporate the fluid-particle interactions （drag force, contact force, Saffman lift force, Magnus lift force, buoyancy force） using momentum exchange and the non-Newtonian behavior of the fluid. The interactions of particle-particle, particle-wall, and particle-drill pipe are taken into account with the Hertz-Mindlin model. We compare the transport of spheres with non-spherical particles （non-smooth sphere, disc, and cubic） constructed via the multi- sphere method for a range of fluid inlet velocities and drill pipe inclination angles. The simulations are carried out for laboratory-scale drilling configurations. Our results demonstrate good agreement with published experimental data. We evaluate the fluid-particle flow patterns, the particle velocities, and the particle concentration profiles. The results reveal that particle sphericity plays a major role in the fluid-solid interaction. The traditional assumption of an ideal spherical particle may cause inaccurate results.展开更多
The reservoir landslide is typically characterized by high-speed movement of a particle-fluid mixture,and its flow and deposit mechanisms are complex.This paper presents the mechanism of submerged granular column coll...The reservoir landslide is typically characterized by high-speed movement of a particle-fluid mixture,and its flow and deposit mechanisms are complex.This paper presents the mechanism of submerged granular column collapse under different densities ambient fluids based on coupled computational fluid dynamics and discrete element method(CFD-DEM)analysis.Important fluid-particle interaction forces,such as the drag force and the buoyancy,are considered by exchanging interaction forces between the CFD and DEM computations.We focus on the flow and deposit characteristics of submerged granular column collapse,namely the runout distance,the tail end height,the particle velocity,the energy,and deposit morphology,which are analyzed qualitatively and quantitatively.The change in fluid field caused by submerged granular column collapse and the formation of eddies are also discussed.A relatively dense fluid can significantly hinder the motion of granular flow,but can improve the conversion efficiency of kinetic energy from the vertical to the horizontal direction.Moreover,the eddies caused by fluid turbulence erode the surface of the granular pile,which is especially marked in a high-density fluid.The findings can provide vital theoretical support for the flow and deposit characteristics of granular flow under fluid and offer insights for the study of reservoir landslides.展开更多
The P1 approximation is a computationally efficient model for thermal radiation.Here,we present a P1 formulation in the context of the combined computational fluid dynamics and discrete element method(CFD-DEM),includi...The P1 approximation is a computationally efficient model for thermal radiation.Here,we present a P1 formulation in the context of the combined computational fluid dynamics and discrete element method(CFD-DEM),including closures for dependent scattering and coarse-graining.Using available analytical and semi-analytical solutions,we find agreement for steady-state and transient quantities in sizedisperse systems.Heat flux is identified as the most sensitive quantity to predict,displaying unphysical spatial oscillations.These oscillations are due to a temperature slip at the locations of abrupt change in solid fraction.We propose two techniques that mitigate this effect:smoothing of the radiative properties,and pseudo-scattering.Furthermore,using up to a million times enlarged particles,we demonstrate practically limitless compatibility with coarse-graining.Finally,we compare predictions made with our code to experimental data for a pebble bed under vacuum conditions,and in presence of nitrogen.We find that a carefully calibrated simulation can replicate trends observed in experiments,with relative temperature error of less than 10%.展开更多
Uniform crushed straw throwing and seed-sowing machines can achieve the processes of straw chopping,straw transport,sowing,fertilization,and straw mulching at the same time,which is widely used in many areas of China....Uniform crushed straw throwing and seed-sowing machines can achieve the processes of straw chopping,straw transport,sowing,fertilization,and straw mulching at the same time,which is widely used in many areas of China.Conveying device is one of the important components used to convey,elevate and throw straw.However,the problems of high power consumption and congestion affect the promotion of the machine.Therefore,the conveying device of uniform crushed straw throwing and seed-sowing machine was analyzed in order to determine its device operation mechanism.Kinematic and dynamic analyses of particles of crushed rice straw during lifting and dispersion are used to develop a flexible-body model of rod-shaped and agglomerate-shaped crushed straw and a coupling model including the mechanical structure of the device.By integrating computational fluid dynamics and the discrete element method,the gas-solid coupling theory in numerical simulations and motion analysis of crushed straw particles is used to determine how the flow field and motion characteristics affect the conveying performance.Besides,regression equations to describe the relationships between the factors and each assessment index were established by using the regression analysis and response surface analysis with the software Design-Expert.The effect of throwing blade speed X_(1),conveying volume of crushed straw X_(2),and pipeline diameter X_(3) on the throwing speed of crushed straw Y_(1) and specific power consumption Y2 were investigated.The highest throwing speed of crushed straw and lowest specific power consumption are the optimization goal.The results of optimization showed that the predict the best optimal parameters were 2000 r/min throwing blade rotational speed,1.4 kg/s conveying volume,and 220 mm pipeline diameter,the planter achieved a throwing speed of 12.2 m/s and specific power consumption of 9179 m^(2)/s^(2).And then a field test verification was conducted.The planter achieved a throwing speed 12.4 m/s and specific power consumption 9070 m^(2)/s^(2) while selecting the best optimal parameters.Thus,the optimal parameters can provide a high-performance operation and satisfy the actual operation requirements The results provide a theoretical basis and data support for seeding technology innovation and equipment optimization to ensure uniform crushed straw throwing in dense rice stubble fields.展开更多
The swelling particle grouting material has demonstrated remarkable plugging effectiveness in high-pressure and large-flow burst water within karst pipelines.Currently,current research on the rheolog-ical model,flow c...The swelling particle grouting material has demonstrated remarkable plugging effectiveness in high-pressure and large-flow burst water within karst pipelines.Currently,current research on the rheolog-ical model,flow computation theory,and plugging mechanism of this material is lacking.The conven-tional grouting slurry diffusion process,using the liquid-liquid two-phase flow method,fails to accurately simulate high solubility slurry and particle swelling.To address these limitations,this study established a precise constitutive model to describe the swelling particle slurry diffusion process in dynamic water.Additionally,a coupling calculation method was proposed to analyze the spatiotemporal heterogeneity of viscosity during slurry diffusion by considering the migration of slurry and the changes in viscosity.To investigate the interaction between particle swelling and flow field changes,a Compu-tational Fluid Dynamics-Discrete Element Method(CFD-DEM)coupling model was developed for the diffusion of swelling particle slurry.It is demonstrated that slurry viscosity increases exponentially within the diffusion front as the particle swelling rate rises,and the drag force exhibits an intriguing behavior of initially increasing and then decreasing as the slurry flows through the pipeline.Further-more,the CFD-DEM coupling model proved to be more accurate in describing viscosity distribution and diffusion distance compared to the finite element solution.The primary objective of this paper is to reveal the plugging mechanism and provide theoretical support for the engineering application of the swelling particle grouting material.展开更多
文摘Rock block removal is the prevalent physical mechanism for rock erosion and could affect the stability of dam foundations and spillways.Despite this,understanding of block removal is still inadequate because of the complex interactions among block characteristics,hydraulic forces,and erosive processes acting on the block.Herein,based on a previously conducted physical experiment of erosion of a single rock block,the removal processes of two different protruding blocks are represented by a coupled computational fluid dynamics-discrete element model(CFD-DEM)approach under varied flow conditions.Additionally,the blocks could be rotated with respect to the flow direction to consider the effect of the discontinuity orientation on the block removal process.Simulation results visualize the entire block removal process.The simulations reproduce the effects of the discontinuity orientation on the critical flow velocity inducing block incipient motion and the trajectory of the block motion observed in the physical experiments.The numerical results present a similar tendency of the critical velocities at different discontinuity orientations but have slightly lower values.The trajectory of the block in the simulations fits well with the experimental measurements.The relationship between the dimensionless critical shear stress and discontinuity orientation observed from the simulations shows that the effect of block protrusion becomes more dominant on the block incipient motion with the increase of relative protrusion height.To our knowledge,this present study is the first attempt to use the coupled finite volume method(FVM)-DEM approach for modelling the interaction behavior between the block and the flowing water so that the block removal process can be reproduced and analyzed.
基金funded by the UK Engineering and Physical Sciences Research Council (EPSRC) grant No. EP/V053655/1 RAILSANDING-Modelling Particle Behaviour in the Wheel–Rail Interface the support received from AltairEngineering and Physical Sciences Research CouncilEP/V053655/1, Sadegh Nadimi
文摘To mitigate and alleviate low wheel-rail adhesion,a train-borne system is utilised to deposit sand particles into the wheel-rail interface via a jet of compressed air in a process called rail-sanding.Britain Rail Safety and Standards Board introduced guidelines on the sand particles’shape,size,and uniformity which needs to be adhered to for rail-sanding.To further inves-tigate these guidelines and help improve them,this research presents a parametric study on the particle characteristics that affect the rail-sanding process including density,size and size distribution,coefficient of uniformity,and shape,utilising a coupled computational fluid dynamics-discrete element method(CFD-DEM)model.The efficiency of rail-sanding is esti-mated for each case study and compared to the benchmark to optimise the sand characteristics for rail-sanding.It is concluded that particle size distribution(within the accepted range)has an insignificant effect on the efficiency while increasing particle size or the coefficient of uniformity decreases the efficiency.Particle shape is shown to highly affect the efficiency for flat,compact and elongated particles compared to the spherical shape.The current numerical model is capable of accurately predicting the trends in the efficiency compared to the actual values obtained from full-scale experiments.
基金Rouven Weiler and Dominik Weis for their valuable insights,discussions,and feedback,which contributed to this work.We also gratefully acknowledge the financial support provided by BASF SE.
文摘Accurately predicting heat flux in coarse-grained CFD-DEM simulations is a significant challenge.Specifically,the rates of fluid-particle heat exchange,the effective thermal conductivity of a bed of particles,as well as radiative heat transfer rates are difficult to predict.By using a novel algorithm,we significantly improve the accuracy and stability of such simulations by using a heat exchange limiter.This limiter enables realistic predictions even at time steps that are three orders of magnitude larger than a typical fluid heat relaxation time.Additionally,view-factor-based corrections for radiative heat exchange computations are developed.These corrections ensure an effective thermal bed conductivity with less than 3%error for a coarse-graining ratio of 10.The applicability of the P1 radiation model in coarse-grained settings is also examined,leading to recommendations for the CFD grid resolution to ensure accurate predictions.Our methods significantly enhance stability,accuracy,and computational efficiency,making coarse-grained CFD-DEM simulations more viable for industrial applications.These advancements enable more reliable modeling of high-temperature processes,accelerate optimization studies,and enable virtual equipment design of such processes.
文摘Computational Fluid Dynamics coupled with Discrete Element Method (CFD-DEM) is a commonly used numerical method to model gas-solid flow in fluidised beds and other multiphase systems. A significant limitation of CFD-DEM is the feasibility of the realistic simulation of large numbers of particles. Coarse-graining (CG) approaches, through which groups of multiple individual particles are represented by single, larger particles, can substantially reduce the total number of particles while maintaining similar system dynamics. As these three CG models have not previously been compared, there remains some debate, however, about the best practice in the application of CG in CFD-DEM simulations. In this paper, we evaluate the performance of three typical CG methods based on simulations of a bubbling fluidised bed. This is achieved through the use of a numerical validation framework, which makes full use of the high-resolution 3D positron emission particle tracking (PEPT) measurements to rigorously validate the outputs of CFD-DEM simulations conducted using various different coarse-graining models, and various different degrees of coarse-graining. The particle flow behaviours in terms of the particle occupancy field, velocity field, circulation time, and bubble size and velocity, are comprehensively analysed. It is shown that the CG simulation starts to fail when the size ratio between the bed chamber and the particles decreases to approximately 20. It is also observed, somewhat surprisingly, that the specific CG approach applied to interparticle contact parameters does not have a substantial effect on the simulation results for the bubbling bed simulations across a wide range of CG factors.
基金supported by the Fundamental Research Funds for the Central Universities(Grant No.2023YJS053)the National Natural Science Foundation of China(Grant No.52278386).
文摘To fundamentally alleviate the excavation chamber clogging during slurry tunnel boring machine(TBM)advancing in hard rock,large-diameter short screw conveyor was adopted to slurry TBM of Qingdao Jiaozhou Bay Second Undersea Tunnel.To evaluate the discharging performance of short screw conveyor in different cases,the full-scale transient slurry-rock two-phase model for a short screw conveyor actively discharging rocks was established using computational fluid dynamics-discrete element method(CFD-DEM)coupling approach.In the fluid domain of coupling model,the sliding mesh technology was utilized to describe the rotations of the atmospheric composite cutterhead and the short screw conveyor.In the particle domain of coupling model,the dynamic particle factories were established to produce rock particles with the rotation of the cutterhead.And the accuracy and reliability of the CFD-DEM simulation results were validated via the field test and model test.Furthermore,a comprehensive parameter analysis was conducted to examine the effects of TBM operating parameters,the geometric design of screw conveyor and the size of rocks on the discharging performance of short screw conveyor.Accordingly,a reasonable rotational speed of screw conveyor was suggested and applied to Jiaozhou Bay Second Undersea Tunnel project.The findings in this paper could provide valuable references for addressing the excavation chamber clogging during ultra-large-diameter slurry TBM tunneling in hard rock for similar future.
基金support from the National Natural Science Foundation of China(grant number 52278403)the Natural Science Foundation of Shandong Province(grant number ZR2021ME135).
文摘The rapid discharge of cuttings from the air cushion chamber is crucial for the construction safety and excavation efficiency of the slurry shield.Previous studies have mainly focused on the transport characteristics of cuttings in the slurry discharge pipe,while the complete process of cuttings entering the slurry discharge pipe from the air cushion chamber until they are discharged was often overlooked.Based on the CFD-DEM coupling method and combined with actual engineering,this paper established a numerical model that can more completely reflect the cuttings discharge process during slurry shield tunneling,and the effects of the slurry velocity at slurry gate and scouring pipes,inclination angle of slurry discharge pipe,cuttings diameter and shape coefficient were investigated by analyzing the variation in the mass flow rate,mass of discharged cuttings,and discharged ratio.The results revealed that increasing the slurry velocity can promote the discharge of cuttings.To keep the discharged rate at a high level,it is recommended that the slurry velocity at the slurry gate should be greater than 0.15 m/s.Reducing the inclination angle of the discharge pipe is conducive to the rapid discharge of cuttings.Cuttings with large diameter or small shape coefficient are more prone to accumulate in the air cushion chamber and cause clogging risk.The research results not only help to improve engineers'understanding of cuttings discharge in slurry shield,but also provide practical guidance for formulating relevant construction measures.
基金supported by the National Natural Science Foundation of China (51579163 and 51639007)the National Key R&D Program of China (2017YFC1502504 and 2016YFC0402304)
文摘A large amount of loose debris materials were deposited on the slope of mountainous areas after the 2008 Ms 8.0 Wenchuan earthquake. During and after the earthquake, these loose debris deposits collapsed and slide into valleys or rivers, changing river sediment supply condition and channel morphology. To investigate the mechanisms of granular flow and deposition, the dynamics of slope failure and sediment transportation in typical mountainous rivers of different intersection angles were analyzed with a coupling model of Computational Fluid Dynamics and Discrete Element Method(CFD-DEM). The numerical results show that the change of intersection angle between the granular flow flume and the river channel can affect the deposit geometry and the fluid flow field significantly. As the intersection angle increases, the granular velocity perpendicular to the river channel increases, while the granular velocity parallel to the river channel decreases gradually. Compared to the test of dry granular flow, the CFD-DEM coupling tests show much higher granular velocity and larger volume of sediments entrained in the river. Due to the river flow, particles located at the edge of the deposition will move downstream gradually and the main section of sediments deposition moves from the center to the edge of the river channel. As a result, sediment supply in the downstream river will distribute unevenly. Under the erosion of fluid flow, the proportion of fine particles increases, while the proportion of coarse particles decreases gradually in the sediment deposition. The build-up of accumulated sediment mass will cause a significant increase in water level in the river channel, thus creating serious flooding hazard in mountainous rivers.
基金financial support of the National Natural Science Foundation of China(Grant No.52074332)express their gratitude to project ZR2020YQ36 supported by Shandong Provincial Science Fund for Excellent Young Scholars。
文摘Hydraulic fracturing is a key technology for the development of unconventional hydrocarbon resources.The proppant placement morphology determines the fracture conductivity,thus affecting the reservoir stimulation effect.In this paper,the proppant migration and placement within complex fractures was studied by considering the fracture wall roughness through computational fluid mechanics-discrete element method(CFD-DEM)in numerical simulation,which is a key approach to study the proppant migration and placement.The results show that the proppant placement non-uniformity,proppant migration capacity,and proppant volume filled in the far-end and the secondary branched fracture are enhanced within the rough fracture compared with those within smooth fractures.The proppant migration capacity is increased within the fracture at low inclination angles(<60°)and low approach angles(<90°),and the proppant placement area is larger in the inclined fracture than that in the vertical fracture.The rise of injection rate and fracturing fluid viscosity causes more proppants migrate to far-end or secondary fractures,resulting in a non-proppant area within the near-wellbore fracture.An increase by 1.3 times in the injection rate and 3 times in the fracturing fluid viscosity leads to a decrease by 26.6%and 27%,respectively,in the proppant placement area within the near-wellbore fracture.The staged injection with small size proppants followed by large size proppants increases the proppant placement area in the primary fracture by 13%-26%,and that with large size proppants followed by small size proppants increases the proppant placement area by 19%-25%,which is due to that the latter method facilitates filling of the secondary branched fracture.The injection location mainly affects the proppant filling degree within the near-wellbore fractures.Compared with the upper injection,the middle and lower injection is not beneficial to filling of proppants within the near-wellbore fracture.
基金This work is part of the program of“The research on the optimization and supply-side reliability of oil product logistics system(No.51874325)”,which is funded by the National Natural Science Foundation of China.
文摘The solid particle impurities generated by pipe wall corrosion might deposit at the elbow of hilly pipelines during the production shutdown of oil pipelines.These solid particle impurities will seriously affect the safety of the pipeline operation and the quality of the petroleum products.Thus,it is necessary to study the methods of removing these trapped particles from pipelines.At present,the most common way to remove these solid particle impurities is pigging oil pipelines periodically by utilizing the mechanical pigging method,while the frequent pigging operation will increase the cost and risk of pipeline operation.It is very convenient and economical to remove the accumulated particles out from the pipeline by oil stream,which can be named Hydraulic Pigging Method(HPM).However,the behavior mechanism of particle in flowing oil is still unclear.This motivates the present research on the particles flushed out by the flowing oil.A numerical model governing the trapped particles displacement from the elbow of an inclined oil pipeline is established in the Euler-Lagrangian framework.The simulation is achieved via CFD coupling with DEM.The CFD method is employed to solving the continuous phase flow,while the discrete particle phase is tracked by the DEM.The numerical model is first validated by comparison with results taken from the published literature.From the simulation results,it is observed that the oil stream,carrier phase,can only flush out the solid particles in a certain diameter range under the given operation conditions,and the particles whose diameter beyond that diameter range will cannot be removed out from the pipeline.The influence of the pipe inclined angle,the oil bulk velocity and the particle diameter on the particle migration characteristics is examined in detail.Furthermore,in order to enhance the efficiency of HPM,an Enhanced Hydraulic Pigging Method based on Multi-Physical Field Collaboration(EHPM-MPFC)is also proposed in the present work.The EHPM-MPFC is validated for having high pigging efficiency via the comparison of the migration characteristics of particles during the EHPM-MPFC and HPM process.The present results can provide the guidance to the HPM operation of products pipelines.
文摘Soil erosion is a critical process that is being studied in soil science, hydraulic engineering, and geotech- nical engineering. Among many societal and environmental impacts, soil erosion is a major cause for the failures of bridges. The erodibility of soil is determined by its physical and geochemical properties and is also affected by surrounding biological activities. In most of the current models for soil erosion, erodibility of non-cohesive soil is characterized by its median grain size （Dso）, density, and porosity. The contribution to erodibility of the irregular shape of soil grains, which plays an important role in the mechanical and hydraulic properties of coarse-grained soils, is generally ignored. In this paper, a coupled computational fluid dynamics and discrete element method model is developed to analyze the influence of the shape of sand grain on soil erodibility. A numerical model for the drag force on spherical and non-spherical particles is verified by using the results from physical free settling experiments. Erosion of sand grains of different shapes is simulated in a virtual erosion function apparatus, a laboratory device used to mea- sure soil erodibility. The simulation results indicate that the grain shape has major effects on erodibility. Spherical particles do not show a critical velocity because of their low rolling resistance, but a critical velocity does exist for angular particles owing to grain interlocking. The erosion rate is proportional to the flow velocity for both spherical and non-spherical particles. The simulation result for angular particle erosion is fairly consistent with the experimental observations, implying that grain shape is an important factor affecting the erodibility of non-cohesive soils.
基金Thig work was supportcd by the National Scicnce and Tech nology Major Project of the Ministry of Science and Technology of China(20172X05009-001)the National Natural Science Foun dation of China(No.J1930001,Nu.J1074208,Nu.J1304270,No.51504277.No.51774308 and No.51904321)+2 种基金the Shan dong Provincial Natural Science Foundation(ZR2019JQ21)the ulnllleltdl Resedltl Fulids fU1 the Celldl Uliveisities(Nu.17CX02008A,No.17CX05003,No.18CX02031A,No.18CX07012A and No.19CX05002A)Key Research and Development Plan of Shandong PToVince(2018GSF116009).
文摘An understanding of the particle transport characteristics in a branched network helps to predict the particle distribution and prevent undesired plugging in various engineering systems.Quantitative analysis of particle flow characteristics is challenging in that experiments are expensive and particle flow is difficult to detect without disturbing the flow.To overcome this difficulty,man-made fractal tree-like branched networks were built,and a coupled computational fluid dynamic and discrete element method model was applied.A series of numerical simulations was carried out to analyze the influence of fractal structure parameters of networks on the particle flow characteristics.The joint influence of inertial,shunt capacity and superposition from upstream branches on particle flow was investigated.The injection position at the inlet determined the particle velocity and its future flow path.The particle density ratio,particle size and bifurcation angle had a greater influence on the shunting of K2 branches than that in the K1 level and N_(k22)/N_(k21) reached a maximum at 60°.Compared with a network with an even number of branches,there was a preferential branch when the branch number was odd.The preferential branch effect or asymmetry degree of the level(K2)branches had a more significant impact on particle shunting than that from the upstream branches(K1).
基金supported by the Research Grants Council of Hong Kong (622910)
文摘Fluid-particle interaction underpins important behavior of granular media. Particle-scale simulation may help to provide key microscopic information governing the interaction and offer better understanding of granular media as a whole. This paper presents a coupled computational fluid dynamics and discrete element method （CFD-DEM） approach for this purpose. The granular particle system is modeled by DEM, while the fluid flow is simulated by solving the locally averaged Navier-Stokes equation with CFD. The coupling is considered by exchanging such interaction forces as drag force and buoyancy force between the DEM and CFD. The approach is benchmarked by two classic geomechanics problems for which analytical solutions are available, and is further applied to the prediction of sand heap formation in water through hopper flow. It is demonstrated that the key characteristic of granular materials interacting with pore water can be successfully captured by the proposed method.
基金Supported by National Natural Science Foundation of China(52274020,U21B2069,52288101)General Program of the Shandong Natural Science Foundation(ZR2020ME095)National Key Research and Development Program(2021YFC2800803).
文摘A method to generate fractures with rough surfaces was proposed according to the fractal interpolation theory.Considering the particle-particle,particle-wall and particle-fluid interactions,a proppant-fracturing fluid two-phase flow model based on computational fluid dynamics(CFD)-discrete element method(DEM)coupling was established.The simulation results were verified with relevant experimental data.It was proved that the model can match transport and accumulation of proppants in rough fractures well.Several cases of numerical simulations were carried out.Compared with proppant transport in smooth flat fractures,bulge on the rough fracture wall affects transport and settlement of proppants significantly in proppant transportation in rough fractures.The higher the roughness of fracture,the faster the settlement of proppant particles near the fracture inlet,the shorter the horizontal transport distance,and the more likely to accumulate near the fracture inlet to form a sand plugging in a short time.Fracture wall roughness could control the migration path of fracturing fluid to a certain degree and change the path of proppant filling in the fracture.On the one hand,the rough wall bulge raises the proppant transport path and the proppants flow out of the fracture,reducing the proppant sweep area.On the other hand,the sand-carrying fluid is prone to change flow direction near the contact point of bulge,thus expanding the proppant sweep area.
文摘We investigate the effect of particle shape on the transportation mechanism in well-drilling using a three-dimensional model that couples computational fluid dynamics （CFD） with the discrete element method （DEM）. This numerical method allows us to incorporate the fluid-particle interactions （drag force, contact force, Saffman lift force, Magnus lift force, buoyancy force） using momentum exchange and the non-Newtonian behavior of the fluid. The interactions of particle-particle, particle-wall, and particle-drill pipe are taken into account with the Hertz-Mindlin model. We compare the transport of spheres with non-spherical particles （non-smooth sphere, disc, and cubic） constructed via the multi- sphere method for a range of fluid inlet velocities and drill pipe inclination angles. The simulations are carried out for laboratory-scale drilling configurations. Our results demonstrate good agreement with published experimental data. We evaluate the fluid-particle flow patterns, the particle velocities, and the particle concentration profiles. The results reveal that particle sphericity plays a major role in the fluid-solid interaction. The traditional assumption of an ideal spherical particle may cause inaccurate results.
基金supported by the National Natural Science Foundation of China(No.51825905).
文摘The reservoir landslide is typically characterized by high-speed movement of a particle-fluid mixture,and its flow and deposit mechanisms are complex.This paper presents the mechanism of submerged granular column collapse under different densities ambient fluids based on coupled computational fluid dynamics and discrete element method(CFD-DEM)analysis.Important fluid-particle interaction forces,such as the drag force and the buoyancy,are considered by exchanging interaction forces between the CFD and DEM computations.We focus on the flow and deposit characteristics of submerged granular column collapse,namely the runout distance,the tail end height,the particle velocity,the energy,and deposit morphology,which are analyzed qualitatively and quantitatively.The change in fluid field caused by submerged granular column collapse and the formation of eddies are also discussed.A relatively dense fluid can significantly hinder the motion of granular flow,but can improve the conversion efficiency of kinetic energy from the vertical to the horizontal direction.Moreover,the eddies caused by fluid turbulence erode the surface of the granular pile,which is especially marked in a high-density fluid.The findings can provide vital theoretical support for the flow and deposit characteristics of granular flow under fluid and offer insights for the study of reservoir landslides.
基金funded through Marie SKEODOWSKA-CURIE Innovative Training Network MATHEGRAM,the People Programme(Marie SKLODOWSKA-CURIE Actions)of the European Union's Horizon 2020 Programme H2020 under REA grant agreement No.813202.
文摘The P1 approximation is a computationally efficient model for thermal radiation.Here,we present a P1 formulation in the context of the combined computational fluid dynamics and discrete element method(CFD-DEM),including closures for dependent scattering and coarse-graining.Using available analytical and semi-analytical solutions,we find agreement for steady-state and transient quantities in sizedisperse systems.Heat flux is identified as the most sensitive quantity to predict,displaying unphysical spatial oscillations.These oscillations are due to a temperature slip at the locations of abrupt change in solid fraction.We propose two techniques that mitigate this effect:smoothing of the radiative properties,and pseudo-scattering.Furthermore,using up to a million times enlarged particles,we demonstrate practically limitless compatibility with coarse-graining.Finally,we compare predictions made with our code to experimental data for a pebble bed under vacuum conditions,and in presence of nitrogen.We find that a carefully calibrated simulation can replicate trends observed in experiments,with relative temperature error of less than 10%.
基金supported by the earmarked fund for CARS-13Natural Science Foundation of Jiangsu Province (Grant No.BK20221187).
文摘Uniform crushed straw throwing and seed-sowing machines can achieve the processes of straw chopping,straw transport,sowing,fertilization,and straw mulching at the same time,which is widely used in many areas of China.Conveying device is one of the important components used to convey,elevate and throw straw.However,the problems of high power consumption and congestion affect the promotion of the machine.Therefore,the conveying device of uniform crushed straw throwing and seed-sowing machine was analyzed in order to determine its device operation mechanism.Kinematic and dynamic analyses of particles of crushed rice straw during lifting and dispersion are used to develop a flexible-body model of rod-shaped and agglomerate-shaped crushed straw and a coupling model including the mechanical structure of the device.By integrating computational fluid dynamics and the discrete element method,the gas-solid coupling theory in numerical simulations and motion analysis of crushed straw particles is used to determine how the flow field and motion characteristics affect the conveying performance.Besides,regression equations to describe the relationships between the factors and each assessment index were established by using the regression analysis and response surface analysis with the software Design-Expert.The effect of throwing blade speed X_(1),conveying volume of crushed straw X_(2),and pipeline diameter X_(3) on the throwing speed of crushed straw Y_(1) and specific power consumption Y2 were investigated.The highest throwing speed of crushed straw and lowest specific power consumption are the optimization goal.The results of optimization showed that the predict the best optimal parameters were 2000 r/min throwing blade rotational speed,1.4 kg/s conveying volume,and 220 mm pipeline diameter,the planter achieved a throwing speed of 12.2 m/s and specific power consumption of 9179 m^(2)/s^(2).And then a field test verification was conducted.The planter achieved a throwing speed 12.4 m/s and specific power consumption 9070 m^(2)/s^(2) while selecting the best optimal parameters.Thus,the optimal parameters can provide a high-performance operation and satisfy the actual operation requirements The results provide a theoretical basis and data support for seeding technology innovation and equipment optimization to ensure uniform crushed straw throwing in dense rice stubble fields.
基金supported by the National Natural Science Foundation of China for the Young Scientists Fund,China (grant No.52109126)the National Key Research and Development Program of China,China (grant No.202103AA080016).
文摘The swelling particle grouting material has demonstrated remarkable plugging effectiveness in high-pressure and large-flow burst water within karst pipelines.Currently,current research on the rheolog-ical model,flow computation theory,and plugging mechanism of this material is lacking.The conven-tional grouting slurry diffusion process,using the liquid-liquid two-phase flow method,fails to accurately simulate high solubility slurry and particle swelling.To address these limitations,this study established a precise constitutive model to describe the swelling particle slurry diffusion process in dynamic water.Additionally,a coupling calculation method was proposed to analyze the spatiotemporal heterogeneity of viscosity during slurry diffusion by considering the migration of slurry and the changes in viscosity.To investigate the interaction between particle swelling and flow field changes,a Compu-tational Fluid Dynamics-Discrete Element Method(CFD-DEM)coupling model was developed for the diffusion of swelling particle slurry.It is demonstrated that slurry viscosity increases exponentially within the diffusion front as the particle swelling rate rises,and the drag force exhibits an intriguing behavior of initially increasing and then decreasing as the slurry flows through the pipeline.Further-more,the CFD-DEM coupling model proved to be more accurate in describing viscosity distribution and diffusion distance compared to the finite element solution.The primary objective of this paper is to reveal the plugging mechanism and provide theoretical support for the engineering application of the swelling particle grouting material.
