^{1}

^{2}

^{1}

One of the main concerns in particle pneumatic conveying process is the possibility of hazards for operation safety due to the electrostatic charge generation as a result of collisions between particles and the equipment wall. Indeed, the electrostatic discharge can occur in the equipment leading to fire or explosion. Simulation of these kinds of processes plays an important role in understanding the various aspects of the system in order to production loss prevention. This paper deals with the simulation of particle pneumatic conveying process inside an inclined tube using a particular method. In this method, the electrification of particles inside the tube is influenced by the vertical collision velocity against the tube wall. Simulation of the particle movements inside the tube, generation of electrostatic charges over the particle surfaces as well as the possibility of fire as a result of discharging the electrostatic energy are investigated. The possibility of fire is investigated by comparing the amount of electrostatic energy with minimum ignition energy (MIE) of the particles. The effect of particle properties including the size and mechanical ones in the simulation is studied. Finally, several solutions are proposed to manage the risk of fire and explosion. As results, the electrostatic energy (E) is beyond the MIE, and the electrostatic discharge can occur leading to explosion for the diameters more than 2 mm and also for elasticity constants lower than 140 MPa. Eventually, there is no hazard of fire and explosion, since all calculated electrostatic energy for the change of Particle Poisson’s ratio varying from 0.1 to 0.9 is less than the MIE value for the air flow rate of 10 m
^{3}/h.

One of the common processes in the chemical industry is the pneumatic conveying process. There is always a possible ignition source in such a process for fire and explosion, and thus, a hazard for operation safety in the process industry. Indeed, as a result of collisions between the solid particles and the wall inside the tanks, tubes and other equipment, the electrostatic charge generates over the solid particles surfaces. Afterward, the charged particles are transported, and the electrostatic discharge can occur which sometimes leads to fire and explosion.

One of the least costly and best approaches for understanding these types of risky processes is using process modeling and simulation. Having a complete model of the process in hand, various aspects of the process can be investigated leading to proper and safe design to minimize the risk of hazards.

Several researchers have studied the pneumatic conveying process both numerically and experimentally. Korevaar et al. [

In this study, a particular method is used to predict the electrostatic charge generation during the particle movement in the pneumatic conveying process inside an inclined tube. The hydrodynamic model describing the flow behavior of the gas-solid phase is coupled with the equations governing the particle charging phenomenon. In the developed model, it is shown that the calculation of the electrostatic charge over the particles is dominantly affected by the contact area as a result of colliding the particles and inner pipe wall. The equations used in contact mechanic especially those of based on Hertzian theory [

It is worth underlying that in order to better understand the flow dynamics and electrostatic charging of the particles during the process, the commercially available finite element software COMSOL Multiphysics is used to simulate the motion of charged particles under different operating parameters.

The rest of the paper is organized as follows. In the next section, the developed model focusing on the main governing equations and assumptions are described. The results and discussion about the effect of particle properties are presented in Section 3, while Section 4 concludes the paper.

At the first step, the outline of the simulation is briefly described in this section. As shown in

The model is developed and then solved to study the variation of two important variables, including the generated electrostatic charges and the electrostatic energy of the particles. The variation of these variables is influenced by several physical phenomena, and thus, the model consists of several differential and algebraic equations, the most important of which are conservation equations and equations related to Hertzian Theory.

At the first step, the equations related to the variation of the fluid and particle velocity are solved. Then, with the particle collision velocities toward the tube wall in hand, those parts of the equations responsible for calculating the electrostatic charges and electrostatic energy come to the solution sequence of the model, and finally the simulation is performed.

The developed Multiphysics model would be too complex without considering some assumptions. The key assumptions and the main mathematical equations governing the model are described in the following.

1) The momentum equations

The turbulent flow regime (Compressible flow with Mach number < 0.3) is assumed. Moreover, it is assumed that particle transferring has no effect on the velocity profiles of the particles in the tube.

The momentum equation of the air flow is as follows:

ρ ( u ⋅ ∇ ) u = ∇ ⋅ [ − p l + ( μ + μ T ) ( ∇ u + ( ∇ u ) T − 2 3 ( μ + μ T ) ( ∇ u ) l − 2 3 ρ k l ] + F + ρ g (1)

where p is the pressure, l is the length, µ is the viscosity, ρ is the density, and F denotes to the forces acting on the fluid. Using the “k-ε” model [

∂ ( ρ k ) ∂ t + ∂ ( ρ k u i ) ∂ x i = ∂ ∂ x j ( μ t σ k ∂ k ∂ x j ) + 2 μ t E i j E i j − ρ ε (2)

∂ ( ρ ε ) ∂ t + ∂ ( ρ ε u i ) ∂ x i = ∂ ∂ x j ( μ t σ ε ∂ ε ∂ x j ) + 2 C ε 1 ε k μ t E i j E i j − C ε 2 ρ ε 2 k (3)

where u_{i} the velocity in the corresponding direction is, E_{ij} denotes the rate of deformation and µ_{T} is the eddy viscosity. µ_{T} is calculated as follows:

μ T = ρ C μ k 2 ε (4)

Equations (2) to (4) also consist of some adjustable constants, the values of which have been estimated by numerous iterations of data fitting as follows:

C μ = 0.09 , σ ε = 1.3 , σ k = 1 , C ε 1 = 1.44 , C ε 2 = 1.92

In addition, the particle momentum equation is:

m p d 2 x d t 2 = F C + F D + F G (5)

where m_{p} is the mass of the particle. Moreover, F_{C}, F_{D}, and F_{G} are the contact force (calculated by the particle collision), the drag force and the gravity force, respectively.

The drag coefficient can be calculated by Equation (6).

C D = max { 24 R e ( 1 + 0.15 R e 0.687 ) , 0.44 } (6)

Note that the grinding process as a result of the particle-particle collision is neglected.

2) The equations of electrostatic charge (Q) and electrostatic energy (E) of the particles

The magnitude of the particle charge (as a result of particle and tube wall collision) depends upon the nature of the wall and particles, actual contact area, and the speed of particle towards the wall [

S = π r t 2 = { S max + ( 8 ρ p 3 P m ) π 2 R 4 ( V cos ( θ ) − V e l 2 ) } 0.5 (7)

S_{max} is the maximum contact area (

S max = 5.41 ( 1 − σ p 2 E p − 1 − σ w 2 E w ) 0.4 ρ p 0.4 R 2 ( V cos θ ) 0.8 (8)

In addition, V is the particle velocity, and V_{el} is the velocity at which the elastic limit is reached. R is the radius of the sphere and r_{t} denotes the radius of the contact area. P_{m} denotes the mean pressure over the contact area, and ρ_{p} is the density of the particle. Furthermore, E_{p} and E_{w} are the elasticity constants of the particle and wall. σ_{p} and σ_{w} are also the Poisson’s ratio of the particle and wall, respectively.

Note that, since the particles are much smaller than the tube diameter, it is assumed that the particles collide with a flat plate instead of the curved plane of the tube wall.

The main goal is to estimate the particle charge generation (Q) and electrostatic energy (E) of the particles. Typically, the charging amount of the particle can be described by its charge density. Having the velocity of each particle in hand and with regard to the contact area of each particle calculated by Equation (7), the Q is estimated using the following equation:

Q = S × q (9)

In which, q is the charge density of the particle. Note that the charge density of the particle is the function of its material [

For calculating the electrostatic energy, first, the special capacity (C) is calculated by the following equation:

C = κ ε 0 γ R ′ (10)

In which R ′ (resistance) can be obtained using the following equation:

R ′ = γ L A (11)

where L is the characteristic length and A denotes to the area of the particles.

Note that the nominator in Equation (10) is indeed the charge relaxation time of the particles defined as follow:

τ = κ ε 0 γ (12)

In the final step, the electrostatic energy can be calculated using the following equation:

E = 1 2 Q 2 C (13)

It should be noted that by comparing the electrostatic energy with minimum ignition energy (MIE) of the particles, the possibility of fire can be evaluated.

For conditions presented in ^{3}/s, the air pressure at the outlet is 1.013 × 10^{5} Pa, the initial values of the velocity of the particles are all zero and the tube wall is stationary. Moreover, the particle number transporting into the tube form the feeder is considered 500. The tube wall material is Stainless Steel 360 and the properties of which are utilized in the equations. The analysis of the effect of particle diameter, particle Poisson’s ratio and the particle elasticity constant on the simulation is performed, and the results are presented in next sub-sections.

Variable name | Value | Variable name | Value |
---|---|---|---|

Tube diameter (m) | 0.031 | Particle Young module (MPa) | 7.3 × 10^{2} |

Tube length (m) | 1.91 | Particle density (kg/m^{3}) | 800 |

Number of 90˚ elbow | 1 | Particle Poisson’s ratio | 0.3 |

Simulation time step | 6.84 × 10^{−5} | Particle volumetric resistivity (Ω.m.) | 10^{10} |

Dielectric constant | 6.8 | MIE (mJ) | 25.2 |

Finite element method is used as the numerical solution in this paper. Three-dimensional meshes are created for finite element analysis. Non-uniform meshing is applied to improve the computational efficiency, where the mesh is fine for great changes of variables, elsewhere the mesh is larger (

Note that the MIE is not a fixed value and is influenced by the gas phase condition like moisture content and the presence of flammable vapor. These factors are kept constant in the simulation.

It is worth underlying that the simulation is in dynamic mode, and the variables are the function of both time and the location. In fact, all the variables are changing dynamically and have different values at each time and each position along the tube. Accordingly, only the maximum values of the electrostatic charge and electrostatic energy in each case (which indeed are the worst condition) are shown, that can be considered as the basis of the safe design.

As the first scenario, the effect of particle diameter is studied.

As depicted, by increasing the particle diameter, the amount of particle charges enhances due to the possibility of increasing the contact area of the particles colliding with the tube wall. For the diameters more than 2 mm, the electrostatic energy (E) is beyond the MIE, and the electrostatic discharge can occur leading to fire or explosion. One solution to manage the risk of fire and explosion can be a reduction in the volumetric flow of air leading to reduction in particle velocity, which in turn, decreases the particle electrostatic charge.

As the second scenario, the elasticity constant of the particles is increased linearly, and its effect is studied (

The higher values of the elasticity constants mean the higher stiffness of the particle. The more the stiffness of the particles, the less the contact area after the collision of particles and pipe wall is made. Accordingly, as shown in

Quantitatively speaking, the results demonstrate that an increase of elasticity constant by near 500%, causes an about 100% reduction to particle electrostatic charge (for the range of values considered herein).

As the final scenario, we focused on the influence of particle Poisson’s ratio. Accordingly, the particle Poisson’s ratio is left to vary while other parameters of the simulation are kept constant. The simulation is performed for an air flow rate of 10 m^{3}/h.

As shown in the results, it can be concluded that the particle Poisson’s ratio plays an important role in charge generation over the surface of particles inside the tube. The increase of Poisson’s ratio leads to the reduction of particle electrostatic charge generation. For the air flow rate of 10 m^{3}/h, there is no hazard of fire and explosion, since all calculated electrostatic energy for the change of Particle Poisson’s ratio varying from 0.1 to 0.9 is less than the MIE value. However, this is mainly due to the low value of the air flow rate. As shown in ^{3}/h, there is a hazard of electrostatic discharge for Poisson’s ratios less than 0.7. When working with high air flow rate and with the particles with low Poisson’s ratios, reduction in the size or the number of the particles may reduce the risk of fire and explosion.

This work dealt with the simulation of particle pneumatic conveying process inside an inclined tube, focusing on prediction of electrostatic charge generation over the particles in various conditions.

In the developed model, the velocity of the particles inside the tube was obtained in the first step. Then, using the Hertzian contact theory, the contact area as a result of the collision between the particles and the pipe wall was calculated. Having calculated the contact area, the amount of charge of the particles at a single collision was calculated regarding the charge density of the particles. Afterward, the electrostatic energies of the particles as a result of the accumulation of the electrostatic charge were estimated. By comparing the amount of electrostatic energy with MIE of the particles, the possibility of the fire was studied. Finally, the effect of important parameters on the generation of particle charge, the

amount of electrostatic energy as well as the possibility of fire and explosion were studied. It was shown that by increasing the particle diameter, the electrostatic energy was linearly increased. On the other hand, the linear increase in the stiffness of the particles brings about the reduction of electrostatic charge generation. Indeed, the particles with higher stiffness had less contact area when colliding with the pipe wall, and consequently less generated values of electrostatic charge.

As results, for the diameters more than 2 mm, the electrostatic energy (E) is beyond the MIE, and the electrostatic discharge can occur leading to fire or explosion. Also, for elasticity constants lower than 140 MPa the electrostatic energy is more than MIE, and the possibility of fire is high. Also, for the air flow rate of 10 m^{3}/h, there is no hazard of fire and explosion, since all calculated electrostatic energy for the change of Particle Poisson’s ratio varying from 0.1 to 0.9 is less than the MIE value. However, by increasing the air flow rate from 10 to 30 m^{3}/h, there is a hazard of electrostatic discharge for Poisson’s ratios less than 0.7. Finally, when working with high air flow rate and with the particles with low Poisson’s ratios, reduction in the size or the number of the particles may reduce the risk of fire and explosion.

The authors declare no conflicts of interest regarding the publication of this paper.

Fathollahi, M., Ahmadi, M.E. and Hosseini, S.M.J. (2019) Electrostatic Charge Generation in Pneumatic Conveying Process: Effect of Particle Properties. Journal of Power and Energy Engineering, 7, 1-11. https://doi.org/10.4236/jpee.2019.76001