Prediction of Horizontal Pneumatic Conveying of Large Coal Particles Using Discrete Phase Model

Yang, Daolong; Li, Ge; Wang, Yanxiang; Wang, Qingkai; Li, Jianping; Huang, Qianqian; Xia, Youtao; Li, Qian

doi:https://doi.org/10.1155/2020/1967052

Advances in Materials Science and Engineering

On this page

Abstract Introduction Results and Analysis Discussion Conclusions Data Availability Disclosure Conflicts of Interest Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2020 | Article ID 1967052 | https://doi.org/10.1155/2020/1967052

Prediction of Horizontal Pneumatic Conveying of Large Coal Particles Using Discrete Phase Model

Daolong Yang,^1,2Ge Li,³Yanxiang Wang,^1,2Qingkai Wang,²Jianping Li,³Qianqian Huang,¹Youtao Xia,¹and Qian Li¹

Academic Editor: Alicia E. Ares

Received07 Nov 2019

Revised12 Mar 2020

Accepted06 Apr 2020

Published12 May 2020

Abstract

The pneumatic conveying focusing on gas-solid two-phase flow plays an important role in a conveying system. Previous work has been conducted in the fields of small particles, where the size was less than 5 mm; however, there are few studies regarding large sizes (>5 mm). In order to predict the horizontal pneumatic conveying of large coal particles, the coupling methods based on the Euler–Lagrange approach and discrete phase model (DPM) have been used for the simulated research. Compared with the experimental results under the same working condition, the particle trajectory obtained by simulation is similar to the particle distribution at the same position in the experiment, and it turns out that the simulation method is feasible for the horizontal pneumatic conveying of large particles. Multifactor simulations are also carried out to analyse the effects of particle size, flow field velocity, solid-gas rate, and pipe diameter on the wall abrasion during horizontal pneumatic conveying, which provides simulation reference and design guide for pneumatic conveying of large particles.

1. Introduction

For a long time, research on fluid flow containing particles has become important in engineering [1–4]. The pneumatic conveying focusing on gas-solid two-phase flow plays an important role in a conveying system. Pneumatic conveying is widely used in energy [5, 6], the chemical industry [7], metallurgy [8, 9], food processing [10], and other fields. A two-phase flow system is more complex than a single-phase flow system due to the interface effects, relative velocity between the gas and solid, and the randomness of the phase interface. In recent years, scholars have conducted many studies on pneumatic conveying.

Advances in computational fluid mechanics have provided a basis for further insight into the dynamics of pneumatic conveying. Currently, there are two main approaches for pneumatic conveying [11–13]: the Euler-Euler approach and the Euler–Lagrange approach. In the Euler-Euler approach, the different phases are treated mathematically as interpenetrating continua [14, 15]. In the Euler–Lagrange approach, the fluid phase is treated as a continuum by solving the Navier–Stokes equations, while the solid phase is solved by tracking a large number of particles through the calculated flow field [16, 17]. The solid phase can exchange momentum, mass, and energy with the fluid phase. The particle trajectories are computed individually at specified intervals during the fluid phase calculation. For applications such as these, particle-particle interactions can be included using the discrete phase model (DPM). The DPM is suitable for simulating granular matter [18, 19]. Such simulations are characterized by a high volume fraction of particles, where particle-particle interaction is important. Previous work has been conducted in the fields of fine particle [20, 21], powder [22–24], and seed [25, 26], in which the size was less than 5 mm and in the range of Geldart A to Geldart C.

This paper mainly studies the pneumatic conveying behavior of large-size coal particles. At present, some scholars have carried out simulation and experimental research [27]. Zhou [28–30] used CFD-DEM numerical simulation to study the influence of particle shape and flow field flow pattern on the fragmentation of massive coal particles in pneumatic conveying. At the same time, in order to reduce the air flow speed in pneumatic conveying of coal particles, the influence of swirl strength on the capture speed of bulk coal particles was studied through experiments, and the particle breaking process and energy change law in the large coal particle pneumatic conveying with different pipeline structure and swirl strength were studied. Yang [8, 9] first applied the pneumatic conveying to the working face of the coal auger in the extremely thin coal seam and proposed the coal and gangue pneumatic conveying and filling system including the underground separation, the underground transportation, and the pneumatic filling; at the same time, the author designed the precise pressurization system, the high-efficiency dust removal system, and the antiblocking system for the pneumatic conveying of coal particles, which could improve the safety and reliability of pneumatic conveying system for coal with large particle. Yang [31, 32] carried out simulation and experimental research on the pneumatic suspension behavior of large irregular coal particles, obtained the suspension speed of coal particles under different particle sizes, and studied the influence of structural parameters of 5–30 mm coal particle gas-solid injection feeder on the pure flow field injection performance and particle injection performance through multifactor orthogonal experiment. Li [33] used CFD-DEM coupling method to carry out numerical simulation of swirling pneumatic conveying of large particle size coal, which shows that swirling pneumatic conveying can improve the integrity of particle pneumatic conveying and extend the service life of pipeline. The differences between the above studies and this paper are shown in Table 1.

As shown in Table 1, most simulation studies can only be based on a small number of particles because simulation using CFD-DEM method takes a lot of time. In order to predict the horizontal pneumatic conveying for large coal particles and reduce the simulation time, this paper uses the coupling method based on the Euler–Lagrange approach, DPM, and the particle trajectory equations. The experiment and the simulation of horizontal pneumatic conveying for large coal particles are carried out to verify the feasibility of the simulation method. Multifactor simulations are carried out to analyse the effects of the particle size, flow field velocity, solid-gas rate, and pipe diameter on the pressure drop and wall abrasion, which provides simulation reference and design guide for pneumatic conveying of large particles.

2. Theory

The flow field provides the energy required by the particles' motion, and the exchange of momentum and energy between the flow field and particles occurs in the pneumatic conveying flow field. The gas phase is a continuous medium, and considering the influence of the solid phase on the flow field, the continuity equation adds the volume fraction term to exclude the gas volume occupied by the solid phase. The solid phase is a discrete phase, and the motion law of solid particles obeys Newton's second law.

2.1. Gas Phase Equations

It is assumed that the temperature of both the gas and solid phases in pneumatic conveying is the same as that of the atmosphere, and no exothermic reaction occurs between the two phases. Therefore, the energy equation of the gas and solid phases can be ignored.

2.1.1. Gas Phase Continuity Equation

According to the law of mass conservation, the gas phase continuity equation can be obtained and is shown in

2.1.2. Gas Phase Momentum Equation

The momentum equation of the gas phase can be obtained from the law of momentum conservation, which is similar to the continuity equation.

The momentum transfer S_m refers to the sum of the fluid drag in the fluid unit.

2.1.3. Turbulence Transmission Equations

The realizable k-ε model [34] has the advantage of more accurate prediction for the divergence ratio of flat and cylindrical jets, and its transmission equations arewith

2.2. Motion Equations of the Solid Phase

2.2.1. Particle Motion Equation

The particle motion equation obtained from Newton's second law is

The drag force F_d of a particle in the uniform flow field is as follows:

The drag coefficient parameter c_D is assumed as a constant, that is, c_D = 0.44, and the Saffman lift force F_S of a particle is as follows:

Equation (10) generated from the velocity gradient is obtained by Li and Ahmadi [35] and was based on the analytical result of Saffman [36].

2.2.2. Particle Trajectory Equations

The force of a single particle in the flow field is shown in Figure 1 [37, 38]. The particle force in the x direction is mainly the fluid drag force; those in the y direction are the particle gravity, fluid drag force, and Saffman lift force; and those in the z direction are the fluid drag force and Saffman lift force. The force equations of a particle acquired are shown in equation (11)–equation (13):

(a)

(b)

For x direction,

For y direction,

For z direction,

There is an assumption that the initial positions of the particles are x_p0, y_p0, and z_p0 at the beginning and that c_D is constant when the particle moves on as time goes by. The particle trajectory equations in the three directions are obtained by integration twice.

For x direction,

For y direction,

For z direction,

The parameters a, b, , f, and n are given by

The particles will be impacted with the wall and other particles in the pipe, and the collision recovery factor is obtained by Forder’s recovery factor equations [39].

2.3. Pipe Abrasion

The pipe abrasion in pneumatic conveying belongs to the contact between the particles and the wall. The abrasion rate of the pipe wall is defined as

The variables d_p, m_p, and are obtained from the boundary conditions and pipe meshing. The variables α and are obtained in the solving process of the particle trajectory. The piecewise linear function of the impact angle f(α) used in the impact between the coal particles and the pipe wall is the abrasion parameter of sand and carbon steel [40].

3. Experiments

The experiment system of horizontal pneumatic conveying is shown in Figure 2. The total length of the conveying pipelines is 10 m, and the diameter is 70 mm. The test system included two pressure transducers, a signal amplifier, a data acquisition instrument, and a computer. The positions of the two pressure transducers and one transparent pipe are shown in Figure 2.

The air flow is pressurized by the screw air compressor and stabilized by the air receiver, and it then enters into the gas-solid injector through the flow valve. The coal particles gain kinetic energy from the air flow and are conveyed into the dust collection box through the conveying pipes. There is some back pressure in the dust collection box, but it is far less than the pressure of the air flow. Therefore, the dust collector can be considered as under atmospheric pressure. The first pressure transducer is installed 3 m from the outlet of the gas-solid injector, and the second is 4 m downstream of the first pressure transducer. The transparent pipe is installed 1.5 m downstream of the first pressure transducer to monitor the particle motion state.

When the output pressure of the air compressor is certain, the opening of the flow valve determines the air flow rate in the pipe. It is called pure flow field when no particles enter into the flow field. However, when the particles enter into the flow field, some of the flow field dynamic pressure turns to static pressure to transfer momentum and energy to the particles. At this point, it amounts to adding back pressure in the pipe. As a result, the air flow rate will be reduced by this back pressure. Therefore, the air flow rate of the pure flow field is regarded as the reference standard in the experiments. When the particles enter into the flow field, the opening of the flow valve will be appropriately increased to complement the reduction.

The size of the experiment coal particles is 5–10 mm, and the feeding rate is controlled by the frequency converter. The density of the experiment coal particles is 2100 kg/m³. The experiment scheme is shown in Table 2.

4. Simulations

4.1. Mesh Independence Simulations

The Euler–Lagrange approach and DPM are used in the simulations. The hexahedral meshes to mesh the 3 m pipe model and the simulation parameters are shown in Figure 3. The mesh parameters of mesh independence simulation are shown in Table 3. The inner diameter of the pipe is 150 mm, the distance between the particle injection surface and the pipe inlet is 400 mm, the distance between the measurement surface and the pipe inlet is 2.5 m, and the simulation time is 2 s. Taking the average velocity of the particles and the average velocity of the flow field of the measurement section as the index, 25 groups of simulation are carried out with 5 particle sizes.

The simulation results of mesh independence show that when the particle size is 2 times of the minimum grid size, the grid size has little effect on the particle velocity, which is also consistent with Aitaliyahia's research results [41]. Therefore, the minimum grid size used in the simulation of different particle sizes in the text is shown in Table 4.

4.2. Horizontal Pneumatic Conveying Simulations

In the horizontal pneumatic conveying simulations, the mesh parameters in Table 4 are used for simulations under different particle sizes and the boundary conditions and injection parameters are shown in Figure 4. The simulation pipe diameters are 70 mm, 100 mm, and 150 mm, and the length is 6 m. The wall roughness of seamless steel pipes used in the experiments is 0.05 mm. The particle density is 2100 kg/m³. The particles are uniformly generated on the particle face and enter the flow field at the speed of V_x = 1 m/s. The transmission medium is air, which is considered as an incompressible gas. The gas density is 1.225 kg/m³, and the dynamic viscosity is 1.8 × 10⁻⁵ kg/m³ (20°C, 1 atm). The simulation considered the interaction between the gas and particle phases; however, the shape characteristics of the particles are ignored.

5. Results and Analysis

5.1. Particle Trajectories in Experiment

A high-speed camera is used to record the particles' movement. The instantaneous images of the experiment videos in the stable conveying part are shown in Figure 5.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

Figure 5

Instantaneous images of experiment videos. (a) No.1, flow field velocity = 21 m/s, feeding rate f_p = 1.12 kg/s. (b) No.2, flow field velocity = 29 m/s, feeding rate f_p = 1.12 kg/s. (c) No.3, flow field velocity = 37 m/s, feeding rate f_p = 1.12 kg/s. (d) No.4, flow field velocity = 44 m/s, feeding rate f_p = 1.12 kg/s. (e) No.1, flow field velocity = 21 m/s, feeding rate f_p = 2.25 kg/s. (f) No.1, flow field velocity = 29 m/s, feeding rate f_p = 2.25 kg/s. (g) No.1, flow field velocity = 37 m/s, feeding rate f_p = 2.25 kg/s. (h) No.1, flow field velocity = 44 m/s, feeding rate f_p = 2.25 kg/s.

When the flow field velocity is small (Figures 5(a) and 5(e)), the images of the stable conveying part under the two feeding rates show an obvious stratification phenomenon. The flow field is divided into three layers: the fluid layer, the conveying layer, and the deposition layer. The fluid layer area is mainly high-velocity fluid. The conveying layer area mixes the particles and the fluid where the flow field promotes the particles' movement. The small particles are mainly in this layer. The deposition layer area is mainly the large particles with a lower velocity. The particle fluidity in this layer is poor, and it is difficult to promote the particles' movement. When the flow field velocity is gradually increased (Figures 5(b) and 5(f)), the images show that the layer boundary appears to be blurred, and the fluid layer and the deposition layer are reduced. The particle movement range increases, and the static and low velocity particles at the pipe bottom decrease. As the flow field velocity continues to increase (Figures 5(c) and 5(g)), the images show that the stratification has basically disappeared. The particle conveying status is relatively dense. The particles at the bottom of the pipe also move with the flow field. While the flow velocity increases further (Figures 5(d) and 5(h)), the images show that the spacings between the particles and the particle velocity are both increasing. There are some small interruptions in the conveying area under a large feed rate (shown in Figure 5(h)), which means that the particle conveying status is changing to the dilute phase state.

5.2. Comparison of the Particle Trajectories between Simulation and Experiment

The simulation and experiment results of particle motion are shown in Table 5. The experiment results are the particle positions from the instantaneous images in Figure 5. The simulation results are the particle trajectories in the range of 5–5.5 m from the particle factory in the simulation pipe. Comparing the experiment and simulation results, the particle motion state in the simulation results is basically similar to that in the experiment, which means that the simulation method is feasible for the horizontal pneumatic conveying of large particles. The feeding rate is 1.12 kg/s in groups No.1 to No.4, and the feeding rate is 2.25 kg/s in groups No.5 to No.8.

It can be found that the particle concentration in the results of simulations and experiments decreases with the flow field velocity under the same feeding rate. In No.1 and No.5, the conveying status is divided into three layers in both the simulation and experiment results due to the small flow field velocity. In the simulation results, the particle trajectories are dense in the pipe bottom, the middle layer has some particle trajectories, and the top layer has little particle trajectories. The stratification phenomenon gradually disappears with the increase in the flow field velocity. It also can be seen from the simulation results that the particle-particle collisions and particle-wall collisions become more severe with the increase in the flow field velocity.

5.3. Simulation Results

5.3.1. Particle Trajectories

The part of the particle trajectories in the horizontal pipe is shown in Figure 6. The solid-gas rate in the simulation is 10. The particle size is different in Figures 6(a) and 6(b), the flow field velocity is different in Figures 6(a) and 6(c), and the pipe diameter is different in Figures 6(a) and 6(d).

(a)

(b)

(c)

(d)

In the four cases, the variation trends of the particle velocity present an acceleration first followed by a deceleration and then an acceleration again with the simulation time. In the position of the particle factory, the particle trajectory is toward the negative direction of the y-axis due to gravity, and the particles are in the accelerated state during the falling process. The particle velocity suddenly decreases upon the first collision of the particles with the pipe wall, and the particles are bounced off the wall and reaccelerated by the flow field. Repeatedly, the particle velocity presents an acceleration first followed by a deceleration and then an acceleration again and finally becomes stabilized. The particle velocity shows an obvious stratification in the y-axis direction. The upper layer particle velocity is larger, and the lower layer particle velocity is smaller. This is due to the small number of particles in the upper layer and fewer particle collisions, so the particles accelerate significantly. However, the number of particles in the lower layer is large, and there are many particle collisions, and the particles are difficult to accelerate, resulting in an obvious stratification of the particle velocity in the y-axis direction.

Comparing Figures 6(a) and 6(b), under the same other conditions, the particle trajectories are messy and dispersed due to the increase in collision probability when the particle size increases. Comparing Figures 6(a) and 6(c), when the flow field velocity increases, the positions of the first particle collisions with the pipe wall are lagging behind. The pure flow field region appears in the upper layer, and the particle trajectories are mainly concentrated in the lower layer. Comparing Figures 6(a) and 6(d), when the pipe diameter increases, the positions of the first particle collisions with the pipe wall are also lagging behind and the pure flow field region above increases.

5.3.2. Pipe Abrasion

The maximum abrasion rate in the pipe is related to the working time. The service life of the pipe decreases with the increase in the maximum abrasion rate. The abrasion simulations are carried out using the abrasion model of equation (20) and considering the effects of the flow field velocity, particle size, solid-gas rate, and pipe diameter. Comparing the particle trajectories in Figure 6, the corresponding pipe abrasion is shown in Figure 7.

(a)

(b)

(c)

(d)

The main abrasion positions are concentrated in the first particle collision area and the pipe downstream. This is because the particles just enter the flow field and accelerate in the flow field during the falling process under the action of gravity. So the particle velocity is large when the particles first impact the pipe wall, which results in the large abrasion of the pipe wall. This is in line with the actual abrasion situation: the actual maximum abrasion areas in the pneumatic conveying system are often in the feeder mixing section where the particles first mix with the flow field. In the pipe downstream, some particles are impacted with the pipe wall again after they rebound in the first collision, which may lead to the high abrasion of the pipe wall.

The maximum abrasion areas of Figures 7(a) and 7(b) appear downstream in the pipe. The maximum abrasion rate decreases with the particle size, and the maximum abrasion position also changes. This is mainly because the increase in particle size reduces the number of particles in the pipe, and the particle collision probability also decreases under the same solid-gas rate.

Comparing Figures 7(a) and 7(c), the maximum abrasion rate increases greatly with the flow velocity. The maximum abrasion areas in Figure 7(c) appear in the positions of the first particle collision and the pipe downstream. The abrasion rate of the middle part of the pipe is less than the other parts. The particles extract more energy from the flow field due to the increase in the flow velocity, which leads to a faster particle velocity and smoother particle trajectories.

Comparing Figures 7(a) and 7(d), the maximum abrasion rate increases with the pipe diameter, and the maximum abrasion areas and abrasion situation are the same as those in Figure 7(c). As the pipe diameter increases, there is more space for the particles to accelerate. The particle velocity is larger upon the first collision, which results in the increase in the abrasion rate. In the middle of the pipe, the particle trajectories are concentrated and smooth, and the particle collisions are also moderated. However, in the pipe downstream, the particles reaccelerate and the abrasion rate increases again.

The abrasion simulations are carried out by using the simulation parameters in Table 5 to obtain the maximum abrasion rate for different particle sizes and flow field velocities. The simulation results are shown in Figure 8.

(a)

(b)

(c)

The maximum abrasion rate increases with the flow field velocity and decreases with the particle size. The maximum value occurs for a particle size of 5 mm and flow field velocity of 60 m/s; the minimum value occurs for a particle size of 25 mm and flow velocity of 20 m/s. When the flow field velocity is constant and the particle size is small, the particle size change has a great influence on the abrasion rate. However, the particle size change has less influence on the abrasion rate when the particle size is larger.

When the particle size is small (d_p = 5 mm and 10 mm), the maximum abrasion rate increases greatly with the flow field velocity. This is because the smaller the particle size is, the greater the particle number at the same solid-gas rate. When the flow velocity increases, the growth rate of the particle velocity is faster and the collision rate between the particles and the pipe wall is greater, which results in a large change in the maximum abrasion rate with the increase in the flow field velocity.

When the particle size is large (d_p = 15 mm, 20 mm, and 25 mm), the maximum abrasion rate increases slightly with the flow field velocity, and the change trend is relatively stable. At this time, the larger the particle size is, the smaller the particle number is. The collision between the particles and the pipe wall is more severe due to the increase in the flow field velocity, but the decrease in the particle number does not make the maximum abrasion rate increase significantly. The effect of particle size and flow field velocity on the maximum abrasion rate indicates that the increase in flow field velocity directly increases the maximum abrasion rate, while the particle size and the particle number have a great influence on the maximum abrasion rate under the same solid-gas rate.

The abrasion simulations are carried out using the simulation parameters in Table 5 to obtain the relations between the maximum abrasion rate, solid-gas rate, and pipe diameter, which are shown in Figure 9.

(a)

(b)

(c)

The maximum abrasion rate increases with the solid-gas rate and the pipe diameter. The maximum value occurs for a solid-gas rate of 20 and pipe diameter of 150 mm; the minimum value occurs for a solid-gas rate of 10 and pipe diameter of 70 mm.

Under the same pipe diameter, the maximum abrasion rate increases with the solid-gas rate. Since the flow velocity is constant, the increase in the solid-gas rate can increase the particle number, resulting in an increase in the collision rate and the abrasion rate.

Under the same solid-gas rate, the maximum abrasion rate increases linearly with the pipe diameter. Since the flow field velocity and the solid-gas rate are fixed, the mass flow rate of the particles and the particle number increase with the pipe diameter, resulting in an increase in the collision probability and the abrasion rate. However, a large pipe diameter leads to enough of an acceleration space for the particles, which also leads to an increase in the collision probability and abrasion rate.

6. Discussion

6.1. Accuracy of the Simulation Results

Comparing the simulation and experiment results of the pressure drop, the pressure drop obtained by simulation and experiment varies greatly with the flow field velocity. This is because the DPM used in this paper does not consider the particle shape characteristics. The particle shape characteristics play a key role in the particle-particle and particle-wall collisions. When the flow velocity is low, the stratification of the flow field leads to more frequent particle-particle and particle-wall interactions, and even some of the particles accumulate at the pipe bottom. Therefore, the accuracy of the DPM simulation results may be lower for low flow field velocity. However, when the flow velocity increases, the particles obtain more energy from the flow field, and some particles are even suspended in the flow field, which leads to a decrease in the collision probability. Therefore, the DPM simulation results are close to the experimental results for high flow field velocity.

6.2. Discussion of Pipe Abrasion

The abrasion experiments take a long time, and it is difficult to observe obvious experiment results. The analysis of pipe abrasion is based on the simulation results carried out in this paper. However, compared with the actual pneumatic conveying industrial system, the simulation results agree with the actual industrial situation in that the first collision position is the maximum abrasion position. Therefore, the DPM abrasion simulation results can be used as a reference for the design of pneumatic conveying systems.

During the pneumatic conveying process of coal particles, the particles will collide with the pipe wall under the action of gravity, which will cause the pipe wall to abrade. Under the same conditions, the flow field velocity, the pipe diameter, and solid-gas rate will increase the maximum abrasion rate. In the case of the pneumatic conveying industry, the abrasion of the pipe wall is an inevitable result of the interaction amount of the system operating conditions and particle collisions. For current pneumatic conveying systems, the service life is often extended by increasing the pipe thickness. However, based on the simulation results of this paper, reasonable configuration of the flow field velocity, particle diameter, pipe diameter, and solid-gas rate parameters help to extend the pipe life and increase the reliability of the pneumatic conveying piping system, but plenty of experiments are needed to confirm these results.

6.3. Limitation

The DPM simulation method used in this paper is only based on ANSYS Fluent software [42, 43]. Compared with the coupling simulation of EDEM software and ANSYS Fluent software [44, 45], it can predict the particle motion and has the advantage of a short calculation time, but the reliability needs to be improved due to ignoring the particle shape.

The shape of particles, such as coal and rock, has diverse characteristics [46], and the shape of particles has a large influence on particle collision and motion in the flow field. In the simulation of this paper, the particle shape characteristic is ignored due to simplification of the equation of particle motion in the flow field. The fluid drag and Saffman lift used in this paper are both based on the present theories, which came from small particles (size less than 5 mm). Therefore, there is a certain gap between the simulation results and the actual experiment. In order to improve the accuracy of the simulation, the next step will consider the influence of particle shape and particle size on the particle motion equation and particle collision behavior.

7. Conclusions

(1)The horizontal pneumatic conveying experiments are carried out with large coal particles. The results show that the flow field exhibits a stratification phenomenon and is divided into three layers when the flow velocity is small. The upper layer is the fluid layer and is mainly high-velocity fluid, the middle is the conveying layer and is a mixture of particles and fluid, and the bottom is the deposition layer and is mainly coal particles with low velocity or that are static. When the flow field velocity increases, the stratification phenomenon disappears.(2)Comparative simulations are carried out using simulation conditions similar to the experiment, and the results show that the DPM model is feasible for predicting the horizontal pneumatic conveying of large coal particles from the particle trajectories.(3)Under the condition of constant solid-gas rate and pipe diameter, the influences of the particle size and flow field velocity on the pipe abrasion are obtained. The maximum abrasion rate increases with the flow field velocity and decreases with the particle size. In the case of constant particle size and flow field velocity, the influences of the solid-gas rate and pipe diameter on the pipe abrasion are obtained. The maximum abrasion rate increases with the pipe diameter and the solid-gas rate.(4)Comparing the simulation and experiment results, the accuracy and limitations of the DPM simulation method used in this paper are discussed. The DPM simulation method has the advantage of a short calculation time; however, there is a certain gap between the simulation results and the actual experiment because the particle shape is not considered and the particle forces are both based on the present theories, which came from small particles. Therefore, the next step will be to consider the influence of the particle shape and particle size on the particle motion equation and particle collision behavior to improve the simulation accuracy and reliability.

Nomenclature

a:	The computing coefficient of drag force
b:	The computing coefficient of Saffman lift force
C(d_p):	The function of particle diameter, C(d_p) = 1.8 × 10⁻⁹
d_lk:	The deformation tenser
e_n:	The normal recovery factor
f:	The coefficient generated from the solution of elliptic integral
F_d:	The drag force vector of particle
F_n:	The normal contact forces
F_o:	The other forces vector of particle
k_n:	The normal contact coefficient
K:	The constant of Saffman lift force, K = 2.594
m_p:	The particle mass
n_k:	The number density of particle phase
:	The velocity component of gas phase on the j direction
:	The relative particle velocity
x_j:	The j direction of coordinate
α:	The normal contact displacement vector
θ:	The impact angle
ν_pki:	The velocity component of solid phase on the i direction
ρ:	The density of gas phase
τ_rk:	The particle relaxation time
:	The area of the cell face at the wall
b():	The function of relative particle velocity, b() = 2.6
d_kl:	The deformation tenser
d_p:	The particle diameter
e_r:	The tangential recovery factor
f(α):	The function of impact angle
:	The gravitational force vector of particle
F_t:	The tangential contact forces
F_S:	The Saffman lift force vector of particle
k_t:	The tangential contact coefficient
m_k:	The quality of single particle
n:	The coefficient generated from the solution of elliptic integral
S:	The volume fraction of solid phase in gas-solid mixture
v_p:	The particle velocity vector
:	The difference value between the gas and solid
α:	The impact angle of the particle path with the wall face
δ:	The tangential contact displacement vector
μ:	The gas dynamic viscosity
ξ:	The volume fraction term
ρ_p:	The density of solid phase.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Disclosure

Daolong Yang and Ge Li are co-first authors of the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The work was funded by the National Natural Science Foundation of China (51705222 and 51675521), the National Natural Science Foundation of Jiangsu Province (BK20170241), and the State Key Laboratory of Process Automation in Mining and Metallurgy, and the Beijing Key Laboratory of Process Automation in Mining and Metallurgy (BGRIMM-KZSKL-2019-07).

References

N. Huber and M. Sommerfeld, “Modelling and numerical calculation of dilute-phase pneumatic conveying in pipe systems,” Powder Technology, vol. 99, no. 1, pp. 90–101, 1998.
View at: Publisher Site | Google Scholar
K. Li, S. B. Kuang, R. H. Pan, and A. B. Yu, “Numerical study of horizontal pneumatic conveying: effect of material properties,” Powder Technology, vol. 251, pp. 15–24, 2014.
View at: Publisher Site | Google Scholar
R. Dhurandhar, J. P. Sarkar, and B. Das, “The recent progress in momentum, heat and mass transfer studies on pneumatic conveying: a review,” Heat and Mass Transfer, vol. 54, no. 9, pp. 2617–2634, 2018.
View at: Publisher Site | Google Scholar
K. Takabatake, Y. Mori, J. G. Khinast, and M. Sakai, “Numerical investigation of a coarse-grain discrete element method in solid mixing in a spouted bed,” Chemical Engineering Journal, vol. 346, pp. 416–426, 2018.
View at: Publisher Site | Google Scholar
B. J. Azzopardi, K. Jackson, J. P. Robinson, R. Kaji, M. Byars, and A. Hunt, “Fluctuations in dense phase pneumatic conveying of pulverised coal measured using electrical capacitance tomography,” Chemical Engineering Science, vol. 63, no. 9, pp. 2548–2558, 2008.
View at: Publisher Site | Google Scholar
J.-W. Zhou, L.-G. Xu, and C.-L. Du, “Prediction of lump coal particle pickup velocity in pneumatic conveying,” Powder Technology, vol. 343, pp. 599–606, 2019.
View at: Publisher Site | Google Scholar
C. Xu, C. Liang, B. Zhou, and S. Wang, “HHT analysis of electrostatic fluctuation signals in dense-phase pneumatic conveying of pulverized coal at high pressure,” Chemical Engineering Science, vol. 65, no. 4, pp. 1334–1344, 2010.
View at: Publisher Site | Google Scholar
D. Yang, J. Li, B. Xing, and Y. Wang, “Recent patents on gangue pneumatic filling for coal auger mining method,” Recent Patents on Mechanical Engineering, vol. 11, no. 1, pp. 31–40, 2018.
View at: Publisher Site | Google Scholar
D. Yang, B. Xing, J. Li, and Y. Wang, “Recent patents on pressurization and dedusting for pneumatic conveying,” Recent Patents on Mechanical Engineering, vol. 11, no. 3, pp. 180–189, 2018.
View at: Publisher Site | Google Scholar
J. J. Fitzpatrick, “Particle properties and the design of solid food particle processing operations,” Food and Bioproducts Processing, vol. 85, no. 4, pp. 308–314, 2007.
View at: Publisher Site | Google Scholar
A. Sokolichin, G. Eigenberger, A. Lapin, and A. Lübert, “Dynamic numerical simulation of gas-liquid two-phase flows Euler/Euler versus Euler/Lagrange,” Chemical Engineering Science, vol. 52, no. 4, pp. 611–626, 1997.
View at: Publisher Site | Google Scholar
M. Sanjosé, J. M. Senoner, F. Jaegle, B. Cuenot, S. Moreau, and T. Poinsot, “Fuel injection model for Euler-Euler and Euler-Lagrange large-eddy simulations of an evaporating spray inside an aeronautical combustor,” International Journal of Multiphase Flow, vol. 37, no. 5, pp. 514–529, 2011.
View at: Publisher Site | Google Scholar
W. P. Adamczyk, A. Klimanek, R. A. Białecki, G. Węcel, P. Kozołub, and T. Czakiert, “Comparison of the standard Euler-Euler and hybrid Euler-Lagrange approaches for modeling particle transport in a pilot-scale circulating fluidized bed,” Particuology, vol. 15, pp. 129–137, 2014.
View at: Publisher Site | Google Scholar
P. Mellin, Q. Zhang, E. Kantarelis, and W. Yang, “An Euler-Euler approach to modeling biomass fast pyrolysis in fluidized-bed reactors - focusing on the gas phase,” Applied Thermal Engineering, vol. 58, no. 1–2, pp. 344–353, 2013.
View at: Publisher Site | Google Scholar
W. Lou and M. Zhu, “Numerical simulation of gas and liquid two-phase flow in gas-stirred systems based on euler-euler approach,” Metallurgical and Materials Transactions B, vol. 44, no. 5, pp. 1251–1263, 2013.
View at: Publisher Site | Google Scholar
S. Laín and M. Sommerfeld, “Characterisation of pneumatic conveying systems using the Euler/Lagrange approach,” Powder Technology, vol. 235, pp. 764–782, 2013.
View at: Publisher Site | Google Scholar
P. Pepiot and O. Desjardins, “Numerical analysis of the dynamics of two- and three-dimensional fluidized bed reactors using an Euler-Lagrange approach,” Powder Technology, vol. 220, pp. 104–121, 2012.
View at: Publisher Site | Google Scholar
Y. Li, J. Zhang, and L.-S. Fan, “Numerical simulation of gas-liquid-solid fluidization systems using a combined CFD-VOF-DPM method: bubble wake behavior,” Chemical Engineering Science, vol. 54, no. 21, pp. 5101–5107, 1999.
View at: Publisher Site | Google Scholar
B. H. Xu and A. B. Yu, “Numerical simulation of the gas-solid flow in a fluidized bed by combining discrete particle method with computational fluid dynamics,” Chemical Engineering Science, vol. 52, no. 16, pp. 2785–2809, 1997.
View at: Publisher Site | Google Scholar
L. Sanchez, N. Vasquez, G. E. Klinzing, and S. Dhodapkar, “Characterization of bulk solids to assess dense phase pneumatic conveying,” Powder Technology, vol. 138, no. 2–3, pp. 93–117, 2003.
View at: Publisher Site | Google Scholar
F. Zhou, S. Hu, Y. Liu, C. Liu, and T. Xia, “CFD-DEM simulation of the pneumatic conveying of fine particles through a horizontal slit,” Particuology, vol. 16, pp. 196–205, 2014.
View at: Publisher Site | Google Scholar
S. Watano, S. Saito, and T. Suzuki, “Numerical simulation of electrostatic charge in powder pneumatic conveying process,” Powder Technology, vol. 135-136, pp. 112–117, 2003.
View at: Publisher Site | Google Scholar
S. Watano, “Mechanism and control of electrification in pneumatic conveying of powders,” Chemical Engineering Science, vol. 61, no. 7, pp. 2271–2278, 2006.
View at: Publisher Site | Google Scholar
K. Choi, Y. Endo, and T. Suzuki, “Experimental study on electrostatic charges and discharges inside storage silo during loading of polypropylene powders,” Powder Technology, vol. 331, pp. 68–73, 2018.
View at: Publisher Site | Google Scholar
M. Güner, “Pneumatic conveying characteristics of some agricultural seeds,” Journal of Food Engineering, vol. 80, no. 3, pp. 904–913, 2007.
View at: Publisher Site | Google Scholar
H. Ghafori and M. Sharifi, “Numerical and experimental study of an innovative design of elbow in the pipe line of a pneumatic conveying system,” Powder Technology, vol. 331, pp. 171–178, 2018.
View at: Publisher Site | Google Scholar
D. Yang, J. Li, K. Zheng, C. Du, and S. Liu, “Impact-crush separation characteristics of coal and gangue,” International Journal of Coal Preparation and Utilization, vol. 38, no. 3, pp. 127–134, 2018.
View at: Publisher Site | Google Scholar
J.-W. Zhou, Y. Liu, C.-l. Du, and S.-y. Liu, “Effect of the particle shape and swirling intensity on the breakage of lump coal particle in pneumatic conveying,” Powder Technology, vol. 317, pp. 438–448, 2017.
View at: Publisher Site | Google Scholar
J.-w. Zhou, C.-l. Du, and Z.-l. Ma, “Influence of swirling intensity on lump coal particle pickup velocity in pneumatic conveying,” Powder Technology, vol. 339, pp. 470–478, 2018.
View at: Publisher Site | Google Scholar
J. Zhou, Y. Liu, C. Du, S. Liu, and J. Li, “Numerical study of coarse coal particle breakage in pneumatic conveying,” Particuology, vol. 38, pp. 204–214, 2018.
View at: Publisher Site | Google Scholar
D. Yang, B. Xing, J. Li, Y. Wang, N. Hu, and S. Jiang, “Experiment and simulation analysis of the suspension behavior of large (5-30 mm) nonspherical particles in vertical pneumatic conveying,” Powder Technology, vol. 354, pp. 442–455, 2019.
View at: Publisher Site | Google Scholar
D. L. Yang, B. S. Xing, J. P. Li et al., “Experimental study on the injection performance of the gas-solid injector for large coal particles,” Powder Technology, vol. 354, pp. 442–455, 2019.
View at: Publisher Site | Google Scholar
J. Li, F. Zhou, D. Yang, B. Yu, and Y. Li, “Effect of swirling flow on large coal particle pneumatic conveying,” Powder Technology, vol. 362, pp. 745–758, 2020.
View at: Publisher Site | Google Scholar
T.-H. Shih, W. W. Liou, A. Shabbir, Z. Yang, and J. Zhu, “A new k-ϵ eddy viscosity model for high Reynolds number turbulent flows,” Computers & Fluids, vol. 24, no. 3, pp. 227–238, 1995.
View at: Publisher Site | Google Scholar
A. Li and G. Ahmadi, “Dispersion and deposition of spherical particles from point sources in a turbulent channel flow,” Aerosol Science and Technology, vol. 16, no. 4, pp. 209–226, 1992.
View at: Publisher Site | Google Scholar
P. G. Saffman, “The lift on a small sphere in a slow shear flow,” Journal of Fluid Mechanics, vol. 22, no. 2, pp. 385–400, 1965.
View at: Publisher Site | Google Scholar
A. M. Reynolds, “On the formulation of Lagrangian stochastic models for heavy-particle trajectories,” Journal of Colloid and Interface Science, vol. 232, no. 2, pp. 260–268, 2000.
View at: Publisher Site | Google Scholar
D. L. Yang, J. P. Li, C. L. Du, H. X. Jiang, and K. H. Zheng, “Injection performance of a gas-solid injector based on the particle trajectory model,” Advances in Materials Science and Engineering, vol. 23, 2015.
View at: Google Scholar
C. D. Dritselis and N. S. Vlachos, “Large eddy simulation of gas-particle turbulent channel flow with momentum exchange between the phases,” International Journal of Multiphase Flow, vol. 37, no. 7, pp. 706–721, 2011.
View at: Publisher Site | Google Scholar
J. K. Edwards, B. S. McLaury, and S. A. Shirazi, “Evaluation of alternative pipe bend fittings in erosive service,” Proceedings of ASME Fluids Engineering Summer Meeting, vol. 22, 2000.
View at: Google Scholar
D. Aitaliyahia, G. Baruzzi, W. G. Habashi et al., “Anisotropic mesh adaptation: towards user-independent, mesh-independent and solver-independent CFD. Part II. Structured grids,” International Journal for Numerical Methods in Fluids, vol. 39, no. 8, pp. 675–702, 2002.
View at: Publisher Site | Google Scholar
O. Rybdylova, M. A. Qubeissi, M. Braun et al., “A model for droplet heating and its implementation into ANSYS Fluent,” International Communications in Heat and Mass Transfer, vol. 76, pp. 265–270, 2016.
View at: Publisher Site | Google Scholar
M. Thiruvengadam, Y. Zheng, and J. C. Tien, “DPM simulation in an underground entry: Comparison between particle and species models,” International Journal of Mining Science and Technology, vol. 26, no. 3, pp. 487–494, 2016.
View at: Publisher Site | Google Scholar
M. Azimian, M. Lichti, and H.-J. Bart, “Investigation of particulate flow in a channel by application of CFD, DEM and LDA/PDA,” The Open Chemical Engineering Journal, vol. 8, no. 1, pp. 1–11, 2014.
View at: Publisher Site | Google Scholar
A. T. Adema, Y. Yang, and R. Boom, “Discrete element method-computational fluid dynamic simulation of the materials flow in an iron-making blast furnace,” ISIJ International, vol. 50, no. 7, pp. 954–961, 2010.
View at: Publisher Site | Google Scholar
K. Zheng, C. Du, J. Li, B. Qiu, L. Fu, and J. Dong, “Numerical simulation of the impact-breakage behavior of non-spherical agglomerates,” Powder Technology, vol. 286, pp. 582–591, 2015.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2020 Daolong Yang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

1179

Downloads

1072

Citations