# Investigation into the Mechanism of the Steel Particle Feeding Process: A Numerical and Experimental Study

Authors: Lei Li; Xueyang Xing; Fangxiang Wang; Xiaodong Dai
Journal: Geofluids
 (2021)

---

## Abstract

The particle feeding system is a prerequisite for the realization of particle impact drilling technology. Because of the high density, the storage and flow of the steel particle are different from those of the other nonviscous particles. The differential equation of the particle movement was built with the liquid bridge force model and the discrete element method. The dynamic movement process and the distribution state of particles in the high-pressure tank were analyzed. For 1 mm steel particles, the mass flow rate decreased with the increase in water content. For 2 mm and 3 mm steel particles, the water content of 15% and 20% was the dividing point of the mass flow rate from increasing to decreasing. When the water content was 10% and 20%, the mass flow rate increased with the steel particle size. But when the water content was 30% and 40%, the mass flow rate decreased with the steel particle size. The study of the control mechanism of the uniformity and stability of particles showed that the funnel flow was the major reason causing the instability and blocking of the feeding process. This research results can provide a basis for the further improvement of the differential pressure feeding system.


---

## Body

## 1. Introduction

Due to the poor drillability of the hard formation and the strong abrasive formation for the unconventional oil and gas reservoir, many problems are met during the drilling operation of the unconventional oil and gas reservoir, such as the low ROP and the complex downhole problems. The rock breaking efficiency of the steel particles with large velocity and high frequency is satisfactory, so the particle impact drilling technology can significantly improve the mechanical drilling speed of the hard formation and the highly abrasive formation. Therefore, the particle impact drilling technology is a good choice for exploration and development of unconventional reservoirs [1]. The particle feeding system is a guarantee of the particle impact drilling technology [2]. How to feed particles from the high-pressure drilling fluid to the mixer continuously and stably is an interesting research topic [3].The wet particle system is a multiphase structure composed of the dense solid particles and the interstitial fluid [4]. The liquid bridge forms between contacting particles, and it results in viscous force. The increase in viscosity has a negative effect on the flow of the particles in the high-pressure tank [5]. But on the other hand, this strong combination effect is beneficial to the other processes [6]. The gas-solid flow in a pseudo 2D fluidized bed of wet particles with the presence of immersed tubes was studied. The gas-solid flow characteristics around the immersed tubes and the effects of immersed tubes, adhesive force, and fluidization velocity on the fluidization of wet particles were analyzed. The various models and experimental results for liquids at small and moderate viscosities were summarized. Several discrepancies between experiments at small liquid viscosities and models were found, and reasons for these differences were discussed. A fast and accurate framework for simulating spheroids in periodic Stokes flow was built, which was based on the completed double-layer boundary integral formulation [7–10]. In view of the particle impact drilling technology, the steel particles in the high-pressure tank contain a certain amount of water, and it often leads to the phenomenon of particle blocking and arching. High density is a typical feature of the steel particle, but there is no similar research about the particle with high density. The problems of discontinuity, instability, and arching often appear in the particle hydraulic feeding system, and it is due to the lack of research on the particle motion law in the high-pressure tank. Consequently, the reasons for the instability and blockage of particle feeding should be studied.In order to analyze the behavior of the wet particles, it is necessary to determine the geometric properties and the physicochemical properties of the particle and the liquid phase [11–14]. Based on the flow of the wet particles in the high-pressure tank, the motion of the wet particles can be described and the behavior characteristics of the wet particle system can be studied. This research was aimed at researching the particle feeding mode under differential pressure [15–17]. The liquid bridge force model and the discrete particle method were combined to simulate the particle motion [18–20]. The dynamic change process and the particle distribution state of particle motion in the high-pressure tank were analyzed. The internal mechanism of particle feeding instability and blockage was clarified. This research can promote the popularization and application of differential pressure feeding.


## 2. Methodology

### 2.1. Computational Model

#### 2.1.1. Mechanical Analysis

Based on the feeding mode under differential pressure, the movement rule of the particles in the high-pressure tank was studied. The high-pressure fluid entered from the inlet, and a part of the fluid entered the high-pressure tank through the pressure balance pipeline (Figure1). This ensured the pressure in the high-pressure tank consistent with the inlet pressure, and it played a role in preboosting.Figure 1
The working principle of the differential pressure feeding system.


#### 2.1.2. Physical Model

The focus of this research was the movement and the feeding law of particles in the high-pressure tank, and the aim was to study the flow instability and the particle blockage reasons. The physical model of the high-pressure tank was established (Figure2). The calculation model parameters, high-pressure tank parameters, and particle parameters required in the study are shown in Table 1.Figure 2
The experimental equipment and the procedure.Table 1
Calculation parameters.
ItemsParametersUnitsValueSteel particlesDiameterm0.0008-0.003Densitykg/m37400Elastic modulusPa2700000Poisson’s ratio0.3High-pressure tankDiametermm200Diameter of the feeding holemm25Heightmm400Half of apex angle°45°OperationPressureMPa20Diameter of the orificemm9.5CalculationFriction coefficient between particles and tank wall0.3Friction coefficient between particles0.3FluidSurface tensionN/m0.06ViscosityPa·s0.002The soft sphere model of the DEM method was used to study the movement of particles in the high-pressure tank [21–23]. The model of indirect contact force between particles was represented by the spring damper friction device. The resultant forces were composed of the elastic force, the viscous damping force, and the friction force. The relative displacement between particles determined the magnitude of the elastic force. The normal elastic force was determined by the Hertz contact theory model. The tangential elastic force was determined by the Mindlin and Deresiewicz theory [24, 25]. The normal force f→cn,ij and the tangential force f→ct,ij were, respectively, described as the following formula.
(1)f→cn,ij=43E∗R∗0.5δn,ij1.5n→,f→ct,ij=μf→cn,ij1−1−δt,ijδmax1.5t→,where E∗ is the relative elastic model of two particles (MPa). When the elastic modulus of both particles is E, E∗=E/21−ν2. ν is Poisson’s ratio of the particle. R∗ is the relative radius of two particles (m). R∗=RiRj/Ri+Rj, where Ri and Rj are the radius of particle i and particle j, respectively. μ is the particle friction coefficient. δn,ij is the normal displacement (m). δt,ij is the tangential displacement (m). δmax is the maximum tangential displacement without sliding (m). δmax=δn,ijμ2−ν/21−ν.The normal componentf→dn,ij and the tangential component f→dt,ij of the viscous damping force were described by the following formula, respectively.
(2)f→dn,ij=cnV→ij⋅n→n→,f→dt,ij=ctV→ij×n→×n→,where cn is the normal viscous damping coefficient. ct is the tangential viscous damping coefficient. Vij is the relative velocity at the contact of particle i and particle j (m/s).When multiple particles contacted with each other, the contact force and moment of particles were∑F→c and ∑M→.
(3)∑F→c=∑i=1nf→cn,ij+f→dn,ij+f→ct,ij+f→dt,ij,∑M→=∑i=1nR→i×F→t,ij.Because there was a certain amount of fluid between particles, it was easy to form the liquid bridge and generate the liquid bridge force. The liquid bridge forceFl included the static capillary force fcap and the dynamic viscous force fν. The static capillary force was calculated by Fisher’s method [26].
(4)fcap=πρ22Δp0+2πρ2σs=πρ2σs1+ρ2ρ1,ρ1=R01−cosϕ+h/2cosϕ+θ,ρ2=R0sinϕ−1−sinϕ+θρ1,where θ is the contact angle (rad). ϕ is the half filling angle (rad). σs is the fluid surface tension (N/m). Δp0 is the capillary pressure (Pa). ρ1 is the radius of the liquid bridge (m). ρ2 is the neck radius of the liquid bridge (m). R0 is the particle radius (m). h is the particle spacing (m).The dynamic viscous force was analyzed by Washino et al. [27]. The normal component and tangential component of the viscous force of the liquid bridge were, respectively, described by the following formulas.
(5)fνn=6πηVnR024h,(6)fνt=2415lnR02h+2.8764πηR0Vt,where Vn is the normal relative velocity (m/s). Vt is the tangential relative velocity (m/s). η is the viscosity coefficient of fluid (Pa·s).


#### 2.1.3. Motion Equation and Solution

The particlei was affected by the ejecting force Fg, the contact force ∑Fc of surrounding particles to particle i, the liquid bridge force Fl between particles, and the resultant moment ∑M. All the resultant forces and the moments of the particles could be described by Newton’s second law at any time.
(7)∂2xit∂t2=∑Fimi,∂2θit∂t2=∑MiIi,where i=1,2,3⋯,n. xi is the translational displacement (m). θi is the angular displacement (rad). mi is the particle mass (kg). Ii is the moment of inertia (kg·m2).The acceleration at timet was discretized by the central difference scheme, and the displacement at time t could be obtained.
(8)xit+Δt=xit+∂xit+Δt/2∂tΔt,θit+Δt=θit+∂θit+Δt/2∂tΔt.During the process of solving the particle motion equation, the calculation of force and displacement was carried out in the same time step. The displacement of timet+Δt was derived from the acceleration of time t in the form of central difference. The particle surface bore new contact force and moment at the new position. The resultant force and the resultant moment of particle i were calculated, and the results were returned to formulas (5) and (6) for cyclic iteration. The motion form of each particle in the particle aggregate was obtained until the end of the required calculation process.


### 2.2. Indoor Experiment

Based on the above analysis on the dynamic change process of the particles, the particle tank was machined, and the particle feeding experiment system was established. The accuracy of numerical simulation results was verified by the indoor experiment.The experimental equipment mainly included a high-pressure water supply system, a particle feeding experimental system, and an electronic balance (Figure2). The mass flow rate of steel particles represented the feeding rate of particles in this research. The steel particle mass flow rate was defined as the mass of steel particles flowing out of the high-pressure tank in unit time.Take 10 s as the unit time, and weight the total weight of particles and water after the bucket was connected. Remove the water and dry the particles, and weigh the particles. Then, convert it into the steel particle mass flow rate.


## 2.1. Computational Model

### 2.1.1. Mechanical Analysis

Based on the feeding mode under differential pressure, the movement rule of the particles in the high-pressure tank was studied. The high-pressure fluid entered from the inlet, and a part of the fluid entered the high-pressure tank through the pressure balance pipeline (Figure1). This ensured the pressure in the high-pressure tank consistent with the inlet pressure, and it played a role in preboosting.Figure 1
The working principle of the differential pressure feeding system.


### 2.1.2. Physical Model

The focus of this research was the movement and the feeding law of particles in the high-pressure tank, and the aim was to study the flow instability and the particle blockage reasons. The physical model of the high-pressure tank was established (Figure2). The calculation model parameters, high-pressure tank parameters, and particle parameters required in the study are shown in Table 1.Figure 2
The experimental equipment and the procedure.Table 1
Calculation parameters.
ItemsParametersUnitsValueSteel particlesDiameterm0.0008-0.003Densitykg/m37400Elastic modulusPa2700000Poisson’s ratio0.3High-pressure tankDiametermm200Diameter of the feeding holemm25Heightmm400Half of apex angle°45°OperationPressureMPa20Diameter of the orificemm9.5CalculationFriction coefficient between particles and tank wall0.3Friction coefficient between particles0.3FluidSurface tensionN/m0.06ViscosityPa·s0.002The soft sphere model of the DEM method was used to study the movement of particles in the high-pressure tank [21–23]. The model of indirect contact force between particles was represented by the spring damper friction device. The resultant forces were composed of the elastic force, the viscous damping force, and the friction force. The relative displacement between particles determined the magnitude of the elastic force. The normal elastic force was determined by the Hertz contact theory model. The tangential elastic force was determined by the Mindlin and Deresiewicz theory [24, 25]. The normal force f→cn,ij and the tangential force f→ct,ij were, respectively, described as the following formula.
(1)f→cn,ij=43E∗R∗0.5δn,ij1.5n→,f→ct,ij=μf→cn,ij1−1−δt,ijδmax1.5t→,where E∗ is the relative elastic model of two particles (MPa). When the elastic modulus of both particles is E, E∗=E/21−ν2. ν is Poisson’s ratio of the particle. R∗ is the relative radius of two particles (m). R∗=RiRj/Ri+Rj, where Ri and Rj are the radius of particle i and particle j, respectively. μ is the particle friction coefficient. δn,ij is the normal displacement (m). δt,ij is the tangential displacement (m). δmax is the maximum tangential displacement without sliding (m). δmax=δn,ijμ2−ν/21−ν.The normal componentf→dn,ij and the tangential component f→dt,ij of the viscous damping force were described by the following formula, respectively.
(2)f→dn,ij=cnV→ij⋅n→n→,f→dt,ij=ctV→ij×n→×n→,where cn is the normal viscous damping coefficient. ct is the tangential viscous damping coefficient. Vij is the relative velocity at the contact of particle i and particle j (m/s).When multiple particles contacted with each other, the contact force and moment of particles were∑F→c and ∑M→.
(3)∑F→c=∑i=1nf→cn,ij+f→dn,ij+f→ct,ij+f→dt,ij,∑M→=∑i=1nR→i×F→t,ij.Because there was a certain amount of fluid between particles, it was easy to form the liquid bridge and generate the liquid bridge force. The liquid bridge forceFl included the static capillary force fcap and the dynamic viscous force fν. The static capillary force was calculated by Fisher’s method [26].
(4)fcap=πρ22Δp0+2πρ2σs=πρ2σs1+ρ2ρ1,ρ1=R01−cosϕ+h/2cosϕ+θ,ρ2=R0sinϕ−1−sinϕ+θρ1,where θ is the contact angle (rad). ϕ is the half filling angle (rad). σs is the fluid surface tension (N/m). Δp0 is the capillary pressure (Pa). ρ1 is the radius of the liquid bridge (m). ρ2 is the neck radius of the liquid bridge (m). R0 is the particle radius (m). h is the particle spacing (m).The dynamic viscous force was analyzed by Washino et al. [27]. The normal component and tangential component of the viscous force of the liquid bridge were, respectively, described by the following formulas.
(5)fνn=6πηVnR024h,(6)fνt=2415lnR02h+2.8764πηR0Vt,where Vn is the normal relative velocity (m/s). Vt is the tangential relative velocity (m/s). η is the viscosity coefficient of fluid (Pa·s).


### 2.1.3. Motion Equation and Solution

The particlei was affected by the ejecting force Fg, the contact force ∑Fc of surrounding particles to particle i, the liquid bridge force Fl between particles, and the resultant moment ∑M. All the resultant forces and the moments of the particles could be described by Newton’s second law at any time.
(7)∂2xit∂t2=∑Fimi,∂2θit∂t2=∑MiIi,where i=1,2,3⋯,n. xi is the translational displacement (m). θi is the angular displacement (rad). mi is the particle mass (kg). Ii is the moment of inertia (kg·m2).The acceleration at timet was discretized by the central difference scheme, and the displacement at time t could be obtained.
(8)xit+Δt=xit+∂xit+Δt/2∂tΔt,θit+Δt=θit+∂θit+Δt/2∂tΔt.During the process of solving the particle motion equation, the calculation of force and displacement was carried out in the same time step. The displacement of timet+Δt was derived from the acceleration of time t in the form of central difference. The particle surface bore new contact force and moment at the new position. The resultant force and the resultant moment of particle i were calculated, and the results were returned to formulas (5) and (6) for cyclic iteration. The motion form of each particle in the particle aggregate was obtained until the end of the required calculation process.


## 2.1.1. Mechanical Analysis

Based on the feeding mode under differential pressure, the movement rule of the particles in the high-pressure tank was studied. The high-pressure fluid entered from the inlet, and a part of the fluid entered the high-pressure tank through the pressure balance pipeline (Figure1). This ensured the pressure in the high-pressure tank consistent with the inlet pressure, and it played a role in preboosting.Figure 1
The working principle of the differential pressure feeding system.


## 2.1.2. Physical Model

The focus of this research was the movement and the feeding law of particles in the high-pressure tank, and the aim was to study the flow instability and the particle blockage reasons. The physical model of the high-pressure tank was established (Figure2). The calculation model parameters, high-pressure tank parameters, and particle parameters required in the study are shown in Table 1.Figure 2
The experimental equipment and the procedure.Table 1
Calculation parameters.
ItemsParametersUnitsValueSteel particlesDiameterm0.0008-0.003Densitykg/m37400Elastic modulusPa2700000Poisson’s ratio0.3High-pressure tankDiametermm200Diameter of the feeding holemm25Heightmm400Half of apex angle°45°OperationPressureMPa20Diameter of the orificemm9.5CalculationFriction coefficient between particles and tank wall0.3Friction coefficient between particles0.3FluidSurface tensionN/m0.06ViscosityPa·s0.002The soft sphere model of the DEM method was used to study the movement of particles in the high-pressure tank [21–23]. The model of indirect contact force between particles was represented by the spring damper friction device. The resultant forces were composed of the elastic force, the viscous damping force, and the friction force. The relative displacement between particles determined the magnitude of the elastic force. The normal elastic force was determined by the Hertz contact theory model. The tangential elastic force was determined by the Mindlin and Deresiewicz theory [24, 25]. The normal force f→cn,ij and the tangential force f→ct,ij were, respectively, described as the following formula.
(1)f→cn,ij=43E∗R∗0.5δn,ij1.5n→,f→ct,ij=μf→cn,ij1−1−δt,ijδmax1.5t→,where E∗ is the relative elastic model of two particles (MPa). When the elastic modulus of both particles is E, E∗=E/21−ν2. ν is Poisson’s ratio of the particle. R∗ is the relative radius of two particles (m). R∗=RiRj/Ri+Rj, where Ri and Rj are the radius of particle i and particle j, respectively. μ is the particle friction coefficient. δn,ij is the normal displacement (m). δt,ij is the tangential displacement (m). δmax is the maximum tangential displacement without sliding (m). δmax=δn,ijμ2−ν/21−ν.The normal componentf→dn,ij and the tangential component f→dt,ij of the viscous damping force were described by the following formula, respectively.
(2)f→dn,ij=cnV→ij⋅n→n→,f→dt,ij=ctV→ij×n→×n→,where cn is the normal viscous damping coefficient. ct is the tangential viscous damping coefficient. Vij is the relative velocity at the contact of particle i and particle j (m/s).When multiple particles contacted with each other, the contact force and moment of particles were∑F→c and ∑M→.
(3)∑F→c=∑i=1nf→cn,ij+f→dn,ij+f→ct,ij+f→dt,ij,∑M→=∑i=1nR→i×F→t,ij.Because there was a certain amount of fluid between particles, it was easy to form the liquid bridge and generate the liquid bridge force. The liquid bridge forceFl included the static capillary force fcap and the dynamic viscous force fν. The static capillary force was calculated by Fisher’s method [26].
(4)fcap=πρ22Δp0+2πρ2σs=πρ2σs1+ρ2ρ1,ρ1=R01−cosϕ+h/2cosϕ+θ,ρ2=R0sinϕ−1−sinϕ+θρ1,where θ is the contact angle (rad). ϕ is the half filling angle (rad). σs is the fluid surface tension (N/m). Δp0 is the capillary pressure (Pa). ρ1 is the radius of the liquid bridge (m). ρ2 is the neck radius of the liquid bridge (m). R0 is the particle radius (m). h is the particle spacing (m).The dynamic viscous force was analyzed by Washino et al. [27]. The normal component and tangential component of the viscous force of the liquid bridge were, respectively, described by the following formulas.
(5)fνn=6πηVnR024h,(6)fνt=2415lnR02h+2.8764πηR0Vt,where Vn is the normal relative velocity (m/s). Vt is the tangential relative velocity (m/s). η is the viscosity coefficient of fluid (Pa·s).


## 2.1.3. Motion Equation and Solution

The particlei was affected by the ejecting force Fg, the contact force ∑Fc of surrounding particles to particle i, the liquid bridge force Fl between particles, and the resultant moment ∑M. All the resultant forces and the moments of the particles could be described by Newton’s second law at any time.
(7)∂2xit∂t2=∑Fimi,∂2θit∂t2=∑MiIi,where i=1,2,3⋯,n. xi is the translational displacement (m). θi is the angular displacement (rad). mi is the particle mass (kg). Ii is the moment of inertia (kg·m2).The acceleration at timet was discretized by the central difference scheme, and the displacement at time t could be obtained.
(8)xit+Δt=xit+∂xit+Δt/2∂tΔt,θit+Δt=θit+∂θit+Δt/2∂tΔt.During the process of solving the particle motion equation, the calculation of force and displacement was carried out in the same time step. The displacement of timet+Δt was derived from the acceleration of time t in the form of central difference. The particle surface bore new contact force and moment at the new position. The resultant force and the resultant moment of particle i were calculated, and the results were returned to formulas (5) and (6) for cyclic iteration. The motion form of each particle in the particle aggregate was obtained until the end of the required calculation process.


## 2.2. Indoor Experiment

Based on the above analysis on the dynamic change process of the particles, the particle tank was machined, and the particle feeding experiment system was established. The accuracy of numerical simulation results was verified by the indoor experiment.The experimental equipment mainly included a high-pressure water supply system, a particle feeding experimental system, and an electronic balance (Figure2). The mass flow rate of steel particles represented the feeding rate of particles in this research. The steel particle mass flow rate was defined as the mass of steel particles flowing out of the high-pressure tank in unit time.Take 10 s as the unit time, and weight the total weight of particles and water after the bucket was connected. Remove the water and dry the particles, and weigh the particles. Then, convert it into the steel particle mass flow rate.


## 3. Results and Discussions

There were two main flow patterns of the steel particles in the high-pressure tank during the feeding process [28–31]. One was the overall flow, and all the steel particles flowed at the same time during the feeding. It could be described as the “first in-first out” flow pattern, and there was relative movement between the steel particles and the wall. The other was the funnel flow, the center flow was faster, and the boundary flow was slower. The flow area was funnel-shaped during the feeding. In addition to the flow region, there was a “dead zone” in which the steel particles were stationary. Most of the steel particles on the top of the high-pressure tank flowed out of the high-pressure tank before the steel particles near the edge.The reason for the funnel flow was that the upper steel particles flowed from the inner surface to the active area in the middle, where they moved downward and flowed out of the high-pressure tank. As the particles flowed down the tube channel at the upper end of the feeding hole, the upper surface of the particles fell. The upper particles continuously flowed into the tube channel and flowed out of the tank; then, the “advanced backward” flow effect formed. For funnel flow, there was an extreme case that only the vertical center area has flow which is called plug flow.When the bottom of the high-pressure tank opened, the discharge mass flow rate began to increase with time. After a certain period of time, the mass flow rate began to fluctuate, but it basically changed around a constant value. The reason for the pulsation was the formation of dynamic arch during the process of steel particle flow. The arch was unstable, and it resulted in density wave during unloading. The magnitude of the pulsation depended on the contact force in the shear direction.The whole unloading process was divided into two stages, the initial acceleration stage and the flow stability stage until the end of unloading. When the bottom of the high-pressure tank was just opened, the force supporting the column disappeared, and the particles near the feeding hole fell freely. Then, the movement of other steel particles in the tank changed. The initial acceleration period was very short, and most of the unloading process was in stable flow.The discharge mass flow rate was an important index to study the influence of the geometric parameters and the steel particle characteristic parameters. Compared with the funnel flow, the mass flow rate of the whole flow was larger and the fluctuation of the flow rate was smaller. Consequently, its flow was relatively stable.


### 3.1. Analysis of the Particle Motion Process

Figure3 describes the flow state and the vertical velocity distribution of particles in the high-pressure tank at different times. In order to track the movement and the position of particles in the high-pressure tank, the tracer layer particles (blue particles) were set at four different heights. In the velocity distribution diagram of particles, the length of arrow represented the size of velocity, and the density of arrow represented the number of particles participating in the flow. In the process of the particle feeding, the particles in the middle part of the tank flowed faster than the particles near the tank wall. The tracer layer showed a parabola shape.Figure 3
Flow state and vertical velocity distribution of steel particles when the half of apex angle was 45° at different times. The upper part was the flow state of steel particles, and the lower part was the vertical velocity distribution of steel particles. (a) Time was 10 s. (b) Time was 25 s. (c) Time was 40 s.The flow in the middle of the high-pressure tank was the vertical flow (Figure3). The particles in the upper area on both sides of the vertical flow area were easy to enter the hopper under the sliding action, and it could be described as the fast flow area. Below the fast flow area, the particle flow area near the tank wall was the slow flow area. However, due to the friction between the particles and the viscous force of the liquid bridge, the velocity of the particles near the tank wall in the lower part of the slow flow region was almost zero, which could be described as the stagnation area. During the process of funnel flow particle feeding, the tangential friction caused by the relative motion of the particles increased. The friction between the particles and the viscous force of the liquid bridge caused the dissipation of injected kinetic energy of the particles, and this hindered the flow of particles.


### 3.2. Effect of the Water Content on the Mass Flow Rate

The water content of the wet steel particles in the high-pressure tank directly affected the liquid bridge force between the steel particles and the flow characteristics of the steel particles. The mass fraction of 1 mm steel particles with different water contents was calculated (Figure4). For 1 mm steel particles, the effect of the liquid bridge force became more obvious with the increase in the water content. When the moisture content was 30%, the mass flow rate of the steel particles decreased significantly. For example, the steel particle mass flow rate reduced from 39.2% to 15.8% when the unloading was in half. When the moisture content increased to 40%, the steel particles were agglomerated at the beginning of the feeding hole, and the unloading mass fraction remained at 11.2%. Therefore, the mass flow rate of steel particles decreased with the increase in water content (Figure 5).Figure 4
The mass fraction of 1 mm steel particles with different water contents.Figure 5
Comparison of simulation results and experimental results of 2 mm and 5 mm particles with different water contents.For steel particles with relatively small particle size (0.8 mm, 1.2 mm, and 1.4 mm), the above rules were also applicable. However, for steel particles with relatively large particle size (2 mm, 3 mm), the influence of water content on the mass flow rate was different. When the water content increased from 5% to 15%, the mass flow rate of steel particles increased. As the moisture content continued to increase, the mass flow rate of steel particles began to decrease. When the water content increased to 15%, the decrease was very obvious. In order to prove the correctness of the simulation results, the effect of the water content of 2 mm and 3 mm steel particles on the mass flow rate was tested. The experimental results were compared with the simulation results (Figure5).The results showed that the simulation results were in good agreement with the experimental results. The water content of 15% and 20% was the dividing point of the mass flow rate from increasing to decreasing. For 2 mm and 3 mm steel particles, the amount of water needed to build the liquid bridge was relatively large. When the water content was less than 15% and 20%, it was not enough to build a liquid bridge. The water infiltrated the surface of the steel particles. On the contrary, it played a certain role in lubricating the flow of the steel particles, so the mass flow rate increased. With the further increase in water content, the liquid bridges formed between the steel particles, and the viscous force between the steel particles increased. The kinetic energy of the steel particles dissipated, and the movement of steel particles was hindered. So the mass flow rate decreased. For 2 mm steel particles, the inflection point of water content was 15%. For larger steel particles, the inflection point should be different from that for 2 mm steel particles. With the increase in water content, the mass flow rate of steel particles must not have a single downward trend. There must be a process of rising first and then falling, and this was related to the amount of water required for the completion of the bridge.


### 3.3. Effect of the Steel Particle Size on the Mass Flow Rate

When the water content was 10% and 20%, the mass flow rate of the steel particles increased with the increase in steel particle size, and the relationship between particle size and mass flow rate was approximately linear (Figure6). When the water content was 30% and 40%, the mass flow rate of the steel particles decreased with the increase in steel particle size. This was mainly because the liquid bridge force hindered the flow of steel particles.Figure 6
Comparison of the mass flow rate with the different particle sizes.The smaller the particle size was, the more obvious the influence of liquid bridge force was. For 1 mm steel particles, the ratio of liquid bridge force to mass force changed with the steel particle spacing. The ratio was between 0.36 and 1.06. The liquid bridge force and mass force were of the same order of magnitude. The effect of liquid bridge force was significant. For 2 mm steel particles, the ratio of mass force to liquid bridge force was between 0.01 and 0.042, and the effect of the liquid bridge force on the flow of the steel particles was obviously less than that of 1 mm steel particles. The energy dissipation caused by the liquid bridge force hindered the movement of the steel particles and slowed down the movement.


### 3.4. Internal Mechanism of Particle Blockage

The results of the numerical simulation showed that three kinds of effect resulted in the blockage of the feeding hole and the instability of the feeding process, including the viscous effect of the liquid bridge force, the arching effect of the asymmetric flow, and the collapse of the particle aggregate.(1)
The viscous action of the bridge force: the centrifugal pump was used to transport the fluid, and the pressure balance pipeline sent the fluid to the high-pressure tank, so the particles in the high-pressure tank were surrounded by the fluid. Due to the funnel flow, the particles in the high-pressure tank had a slow flow area and a stagnant area. The liquid bridge was easy to form between the particles, and the viscous force between the particles increased. Many particles were connected with each other to form bridges, and the particles aggregated into a group (Figure7)(2)
The arching effect of the asymmetric flow: the particle velocity distribution showed that the flow of the particles in the high-pressure tank was asymmetric under the effect of pressure difference and gravity. Sometimes the flow was fast on the left side, and sometimes the flow was fast on the right side. In the alternating moment of left and right flow, the flow was symmetrical and presented the form of converging flow. The particles rushed to the middle of the feeding hole, and the feeding speed of particles was fastest. It was easy to form an unstable dynamic arch. During the process of forming and destroying the dynamic arch, the feeding hole produced the particle flow density wave which caused the instability of particle feeding. Because the feeding speed of particles was too fast, it was easy to cause particles to squeeze each other to achieve a force balance and form a mechanical arch. Thus, a blocking phenomenon would occur(3)
The collapse of particle aggregate: during the process of particle feeding, the middle part of the high-pressure tank presented a flow state similar to that of the funnel. When the cone of the funnel was steep, the particles in the rapid flow area slid into the funnel with a large speed. It had a compaction effect on the particles in the lower part, or even the whole particle aggregate collapsed. Then, a large number of particles in the feeding hole turned compact, and this would result in the blocking accident of the feeding holeFigure 7
Particle aggregate with different sizes.


## 3.1. Analysis of the Particle Motion Process

Figure3 describes the flow state and the vertical velocity distribution of particles in the high-pressure tank at different times. In order to track the movement and the position of particles in the high-pressure tank, the tracer layer particles (blue particles) were set at four different heights. In the velocity distribution diagram of particles, the length of arrow represented the size of velocity, and the density of arrow represented the number of particles participating in the flow. In the process of the particle feeding, the particles in the middle part of the tank flowed faster than the particles near the tank wall. The tracer layer showed a parabola shape.Figure 3
Flow state and vertical velocity distribution of steel particles when the half of apex angle was 45° at different times. The upper part was the flow state of steel particles, and the lower part was the vertical velocity distribution of steel particles. (a) Time was 10 s. (b) Time was 25 s. (c) Time was 40 s.The flow in the middle of the high-pressure tank was the vertical flow (Figure3). The particles in the upper area on both sides of the vertical flow area were easy to enter the hopper under the sliding action, and it could be described as the fast flow area. Below the fast flow area, the particle flow area near the tank wall was the slow flow area. However, due to the friction between the particles and the viscous force of the liquid bridge, the velocity of the particles near the tank wall in the lower part of the slow flow region was almost zero, which could be described as the stagnation area. During the process of funnel flow particle feeding, the tangential friction caused by the relative motion of the particles increased. The friction between the particles and the viscous force of the liquid bridge caused the dissipation of injected kinetic energy of the particles, and this hindered the flow of particles.


## 3.2. Effect of the Water Content on the Mass Flow Rate

The water content of the wet steel particles in the high-pressure tank directly affected the liquid bridge force between the steel particles and the flow characteristics of the steel particles. The mass fraction of 1 mm steel particles with different water contents was calculated (Figure4). For 1 mm steel particles, the effect of the liquid bridge force became more obvious with the increase in the water content. When the moisture content was 30%, the mass flow rate of the steel particles decreased significantly. For example, the steel particle mass flow rate reduced from 39.2% to 15.8% when the unloading was in half. When the moisture content increased to 40%, the steel particles were agglomerated at the beginning of the feeding hole, and the unloading mass fraction remained at 11.2%. Therefore, the mass flow rate of steel particles decreased with the increase in water content (Figure 5).Figure 4
The mass fraction of 1 mm steel particles with different water contents.Figure 5
Comparison of simulation results and experimental results of 2 mm and 5 mm particles with different water contents.For steel particles with relatively small particle size (0.8 mm, 1.2 mm, and 1.4 mm), the above rules were also applicable. However, for steel particles with relatively large particle size (2 mm, 3 mm), the influence of water content on the mass flow rate was different. When the water content increased from 5% to 15%, the mass flow rate of steel particles increased. As the moisture content continued to increase, the mass flow rate of steel particles began to decrease. When the water content increased to 15%, the decrease was very obvious. In order to prove the correctness of the simulation results, the effect of the water content of 2 mm and 3 mm steel particles on the mass flow rate was tested. The experimental results were compared with the simulation results (Figure5).The results showed that the simulation results were in good agreement with the experimental results. The water content of 15% and 20% was the dividing point of the mass flow rate from increasing to decreasing. For 2 mm and 3 mm steel particles, the amount of water needed to build the liquid bridge was relatively large. When the water content was less than 15% and 20%, it was not enough to build a liquid bridge. The water infiltrated the surface of the steel particles. On the contrary, it played a certain role in lubricating the flow of the steel particles, so the mass flow rate increased. With the further increase in water content, the liquid bridges formed between the steel particles, and the viscous force between the steel particles increased. The kinetic energy of the steel particles dissipated, and the movement of steel particles was hindered. So the mass flow rate decreased. For 2 mm steel particles, the inflection point of water content was 15%. For larger steel particles, the inflection point should be different from that for 2 mm steel particles. With the increase in water content, the mass flow rate of steel particles must not have a single downward trend. There must be a process of rising first and then falling, and this was related to the amount of water required for the completion of the bridge.


## 3.3. Effect of the Steel Particle Size on the Mass Flow Rate

When the water content was 10% and 20%, the mass flow rate of the steel particles increased with the increase in steel particle size, and the relationship between particle size and mass flow rate was approximately linear (Figure6). When the water content was 30% and 40%, the mass flow rate of the steel particles decreased with the increase in steel particle size. This was mainly because the liquid bridge force hindered the flow of steel particles.Figure 6
Comparison of the mass flow rate with the different particle sizes.The smaller the particle size was, the more obvious the influence of liquid bridge force was. For 1 mm steel particles, the ratio of liquid bridge force to mass force changed with the steel particle spacing. The ratio was between 0.36 and 1.06. The liquid bridge force and mass force were of the same order of magnitude. The effect of liquid bridge force was significant. For 2 mm steel particles, the ratio of mass force to liquid bridge force was between 0.01 and 0.042, and the effect of the liquid bridge force on the flow of the steel particles was obviously less than that of 1 mm steel particles. The energy dissipation caused by the liquid bridge force hindered the movement of the steel particles and slowed down the movement.


## 3.4. Internal Mechanism of Particle Blockage

The results of the numerical simulation showed that three kinds of effect resulted in the blockage of the feeding hole and the instability of the feeding process, including the viscous effect of the liquid bridge force, the arching effect of the asymmetric flow, and the collapse of the particle aggregate.(1)
The viscous action of the bridge force: the centrifugal pump was used to transport the fluid, and the pressure balance pipeline sent the fluid to the high-pressure tank, so the particles in the high-pressure tank were surrounded by the fluid. Due to the funnel flow, the particles in the high-pressure tank had a slow flow area and a stagnant area. The liquid bridge was easy to form between the particles, and the viscous force between the particles increased. Many particles were connected with each other to form bridges, and the particles aggregated into a group (Figure7)(2)
The arching effect of the asymmetric flow: the particle velocity distribution showed that the flow of the particles in the high-pressure tank was asymmetric under the effect of pressure difference and gravity. Sometimes the flow was fast on the left side, and sometimes the flow was fast on the right side. In the alternating moment of left and right flow, the flow was symmetrical and presented the form of converging flow. The particles rushed to the middle of the feeding hole, and the feeding speed of particles was fastest. It was easy to form an unstable dynamic arch. During the process of forming and destroying the dynamic arch, the feeding hole produced the particle flow density wave which caused the instability of particle feeding. Because the feeding speed of particles was too fast, it was easy to cause particles to squeeze each other to achieve a force balance and form a mechanical arch. Thus, a blocking phenomenon would occur(3)
The collapse of particle aggregate: during the process of particle feeding, the middle part of the high-pressure tank presented a flow state similar to that of the funnel. When the cone of the funnel was steep, the particles in the rapid flow area slid into the funnel with a large speed. It had a compaction effect on the particles in the lower part, or even the whole particle aggregate collapsed. Then, a large number of particles in the feeding hole turned compact, and this would result in the blocking accident of the feeding holeFigure 7
Particle aggregate with different sizes.


## 4. Conclusion

In this research, the liquid bridge force model and the discrete particle method were combined to simulate the particle motion. The dynamic change process and the particle distribution state of particle motion in the high-pressure tank were analyzed. This research can provide bases for the design of the particle feeding system.(1)
The numerical simulation results and the indoor experiment showed that the water content and the size of steel particles affect the mass flow rate. With the increase in water content, the mass flow rate of steel particles must not have a single downward trend. The smaller the particle size was, the more obvious the influence of liquid bridge force on the flow of steel particles was(2)
There are three kinds of effect that resulted in the blockage of the feeding hole and the instability of particle feeding, including the viscous effect of the liquid bridge force, the arching effect of the asymmetric flow, and the collapse of the particle aggregate(3)
The future research will focus on the structure design of the feeding tank, and the effect of different transport methods on the uniformity of particle injection should be investigated



---

*Source: 6691899-2021-01-22.xml*