Research Article  Open Access
Daolong Yang, Jianping Li, Changlong Du, Hongxiang Jiang, Kehong Zheng, "Injection Performance of a GasSolid Injector Based on the Particle Trajectory Model", Advances in Materials Science and Engineering, vol. 2015, Article ID 871067, 8 pages, 2015. https://doi.org/10.1155/2015/871067
Injection Performance of a GasSolid Injector Based on the Particle Trajectory Model
Abstract
Gassolid injectors are widely used feeding equipment in pneumatic conveying systems. The performance of a gassolid injector has a significant influence on the type of application it can be employed for. To determine the key factors influencing the injection performance and address clogging problems in a gassolid injector during a pneumatic conveying process, the particle trajectory model has been utilised as a means to perform simulations. In the particle trajectory model, the gas phase is treated as a continuous medium and the particle phase is treated as a dispersed phase. In this work, numerical and experimental studies were conducted for different nozzle positions in a gassolid injector. A gassolid injector testbed was constructed based on the results of the simulations. The results show that the nozzle position is the key factor that affects the injection performance. The number of extrusive particles first increases and then decreases with the change in the nozzle position from left to right. Additionally, there is an optimum nozzle position that maximises the injection mass and minimises the number of particles remaining in the hopper. Based on the results of this work, the injection performance can be significantly increased and the clogging issues are effectively eliminated.
1. Introduction
A gassolid injector is an important piece of feeding equipment in pneumatic conveyers, which are widely used in the petroleum and chemical industries, material conveying, power stations, and other departments. They possess a simple structure, no moving parts, and concatenate conveniently with other pieces of equipment [1, 2]. Many methods have been used to simulate gassolid fluidisation. Some examples include the twofluid smoothed particle hydrodynamics method (TFSPH) [3], the multilattice deterministic trajectory model (MLDT) [4], the multifluid model (MFM), the twofluid model (TFM) [5], and the discrete phase model (DPM) [6]. In addition, particle trajectory models (PTMs) are also widely used because they require fewer partial differential equations when addressing the dispersed phase and because they provide exact results when forecasting the particle distribution [7–9]. The particle trajectory model is mainly used to forecast the movement of dilute particle phases in a turbulent flow. The Lagrangian solution model is adopted to address the particle phase. The particleparticle and particlewall collisions are calculated using a statistical method that is not restricted by the number of particles [10, 11].
Practically, gassolid injectors may jam depending on the particle diameter and as the volume being conveyed increases. These issues restrict their application and development when dealing with largesized particles and largemass flow pneumatic conveying. Many scholars have studied the conveying properties and static pressure distributions in gassolid injectors and agree that the location of the driving nozzle and the angle of the converging section have an obvious influence on the maximum achievable mass flow rate [12–15]. The goal of this paper is to improve the injection performance of a gassolid injection (to increase the conveying properties) and determine the key factors that affect the injection performance. The particle trajectories were simulated in a uniform flow field using the PTM and the injection performance at different driving nozzle locations was obtained by analysing the number of injected particles. Finally, a gassolid injector testbed was built to investigate the injection performance and verify the simulation models and results.
2. Mathematical Model
The movement of the particles is determined by the interactions between the solid and gas phases, which are determined in two ways. The first is by the exchange of mass, momentum, and energy. The second is through particleparticle and the particlewall collisions. The particle trajectories are obtained by the PTM, which uses different handling methods for the gas and particle phases. The gas phase is treated as a continuous medium and the variables of the gas phase are obtained by solving the gas control equation. The particle phase is treated as a dispersed phase and the Lagrangian method is used to track the particle trajectories. The mass, momentum, and energy of the gas and particle phases are exchanged during the iterative process [16].
2.1. Gas Control Equation
The gas phase is treated as a continuous medium in the PTM. Therefore, the gas phase continuity and momentum equations are based on the law of conservation of mass and Newton’s second law [17, 18]. In addition, the interactions between the gas and particles are considered. For larger particle diameters of up to 5 mm, additional interactions are added to the gas control continuity equation:where is the density of gas phase, is the coordinate of direction, is the velocity component of the gas phase in the direction, is the particle phase volume fraction in the gassolid twophase mixture, is the number density of the particle phase, and is the single particle mass. The momentum equation is defined aswhere is the dynamic viscosity of the gas phase, is the coordinate of direction, is the velocity component of the gas phase in the direction, is the density of the particle phase in the gassolid twophase mixture, is the velocity component of gas phase in the direction, and is the diffuse relaxation time of the particle phase. Equations (1) and (2) are obtained from [17, 18].
The above control equations are all for the gas phase. There are a total of four control equations when we consider the momentum equations and the three directions in the coordinate system. Therefore, the unknown variables for the gas phase are , , , , and . These are addressed by solving the closed form of the turbulence model [19]. The turbulent kinetic energy and turbulent momentum dissipation equations are shown respectively:Equations (3) and (4) are obtained from Yakhot and Orszag’s work [20].
2.2. Trajectory Equation of Particle Phase
The force analysis on the particle phase was performed using Newton’s second law in the PTM and the particle trajectory equations were acquired by integrating twice.
2.2.1. Motion Equation
The particle motion equation (Figure 1) is shown in the following equation, which is obtained from Newton’s second law:where is the particle velocity, is the average particle mass, is the fluid drag force on the particle, is the force of gravity on the particle, and represents the other forces on the particle.
The other particle forces include the Saffman force, the Brown force, and the Magnus force. In this study, the particle diameter is 5 mm, meaning that the Brown force and the Magnus force have little effect on the particle movement, as described in [21–24]. Therefore, these forces are ignored when considering the particle state of motion.
The fluid drag force is defined as [25]where is the particle density, is the particle diameter, is the fluid drag coefficient, which is defined as when the particle diameter is 5 mm, is the gas velocity, and is the particle velocity.
The Saffman force is defined aswhere is the Saffman force constant, which is defined as and is the deformation tenser.
The Saffman force is generated by the different velocity gradients in the fluid. Equation (7) is obtained from Li and Ahmadi [25], which determines the succession of the Saffman expression [26].
2.2.2. Particle Trajectory
When the particles are in a uniform flow field and do not collide with the walls or other particles, as described in (5), the particle trajectory equations may be defined as follows: direction: direction: direction: assuming that the initial particle velocities in the  and directions are zero. By substituting (6) and (7) into (8)~(10), the particle phase trajectory equations in the , , and directions may be obtained by integrating twice: direction: direction: direction: where , , and are the initial particle coordinates in the , , and directions, is the computational coefficient of fluid drag, is the computational coefficient of the Saffman force, is the coefficient of the simplified equation, , , and are the gas velocities in the , , and directions, is the difference in the initial gas and particle velocities, and and are the coefficients generated when simplifying the integral computation. These coefficients are shown in
The particle velocity recovery factor for the collisions with the wall is determined from Alister’s experiment [27]. There is a close relationship between the recovery factor and the particle impact angle, which is defined as where is the normal recovery factor, is the tangential recovery factor, and is the particle impact angle.
3. Simulations
The flow field of the gas phase in the gassolid injector is calculated using a colocated grid with the SIMPLE method. The particle coordinate positions are solved using the equations of motion. Additionally, the solid volume percentage of the colocated grids is also calculated. Then, the calculated value is returned for the calculation of the flow field. In this method, the coupled solution for the determination of the continuous fluid and particle phases may be alternated.
3.1. Simulation Model and Initial Parameters
The gassolid injector consists of a driving nozzle, feed opening, mixing chamber, contraction section, and delivery pipe (Figure 2). The inner diameter of the driving nozzle is 20 mm, the diameter of the delivery pipe is 14 mm, and the inclination angle is 8°. The nozzle locations addressed by the four models are , , , and , respectively. The initial parameters for the simulation model are shown in Table 1. The feeding time was set to be 0.2 s and 0.4 s, respectively, to prevent the material from becoming plugged in the gassolid injector.

3.2. Simulation Results
One group of steadystate simulations and two groups of transient simulations were conducted for different feeding times and particle trajectories in the gassolid injector to analyse the injection mass.
The coordinates at each time step are solved using the PTM. The particle trajectories are shown in Figure 3 after the particle coordinates are projected into the gassolid injection model. The nozzle locations at , , , and are shown in Figure 3, , , and , respectively.
The injection block not only increases the use cost but also reduces the reliability of the pneumatic conveying system. Reducing the number of particles in the injector may alleviate material blockage issues and achieve a balance between the feeding mass and the extrusive mass. A statistical analysis of the extrusive number and the residual number of particles is conducted separately for each of the four models. The statistical results are shown in Figures 4–7.
The number of extrusive particles and particles remaining in the hopper for a 0.2 s feeding time is shown in Figures 4 and 5, respectively, where the horizontal axis represents the simulation time and the vertical axis represents the number of particles remaining in the hopper.
The number of extrusive particles and the number of particles remaining in the hopper for a 0.4 s feeding time are shown in Figures 6 and 7, respectively, where the horizontal axis represents the simulation time and the vertical axis represents the remaining number of particles in the hopper.
4. Experiment
4.1. Experimental Equipment
The gassolid injector testbed consists of a gas inlet, a bolt, flange, support frame, injector, and hopper (Figure 8). One side of the gas inlet connects the gasholder and the other side inserts into the injector through the flange and is fixed by the bolt, which can move up and down in the axial direction.
4.2. Test Objective and Method
A mixture of sediment and stone was used with the gassolid injector testbed to explore the influence of the nozzle position on the injection performance and verify the simulation results. The material supply pipe had an inner diameter of 80 mm and is divided into eight equal parts along the axial direction of mixture chamber. These parts are denoted by A to I from left to right. The nozzle was placed in a position and the mixture was continually dropped into the hopper. The injected material was gathered every five minutes after a relatively stable injection mass is achieved. The results from the tests are shown in Table 2.

4.3. Experiment Results
The jet gas flow caused a significant number of impacts and there was a great deal of noise when the nozzle was located at point A. The material erupted in reverse from the hopper, leading to the failure of the test. Only a handful of the material was injected into the mixture chamber because the material could not be mixed with the highspeed gas when the nozzle was located at point I. The average injection mass was adjusted by eliminating the data at nozzle locations A and I. The fitted curve is shown in Figure 9 and the relationship between the average injection mass and the nozzle location is expressed in which achieved an squared value of .
5. Discussion
The particle trajectories in the mixing chamber of the gassolid injectors, which are shown in Figure 3, are sparser compared to the other three models when using the PTM, which indicates less particle deposition and particlewall collisions. In contrast, the particle trajectories shown in Figures 3 and 3 are denser, indicating more particle deposition and energy consumption, an unstable flow field, and a greater degree of wear in mixing chamber. The particle trajectories in the delivery pipe, shown in Figure 3, exhibit a greater degree of linearity compared to the other three models when using the PTM, which indicates less collisions between the particles and wall. The other three models exhibit clear particleparticle or particlewall collisions. In general, a greater number of collisions lead to a greater degree of wear and higher energy consumption in the pipes. As a result, the material conveying process becomes more difficult.
The particles are injected over 0.15 s and the number of extrusive particles increases from 0.18 s to 0.22 s. From 0.25 s to 0.34 s, the particle count increases steadily before beginning a downward trend after 0.35 s and sharply declines during the 0.38 s~0.50 s period (Figure 4). The number of extrusive particles is largest when the nozzle is at . Additionally, a small increase occurs from 0.50 s to 0.75 s when the nozzle is at and . This illustrates that, during this time, the sedimentary particles at the bottom of injector reenter the flow field and are ejected out of the mixing chamber.
The number of remaining particles is less using the PTM than that for the other three injector models when the nozzle location is at after 0.22 s (Figure 5). Remarkably, the number of particles for all four models peaks at 0.30 s (after the feeding stops). The curves are similar for nozzle locations at and after 0.30 s. The curves are also similar when the nozzle location is at and . However, the number of particles remaining in the hopper in the latter is less than the former after 0.70 s to the end of the injection period because the flow field experienced acute fluctuations when the sedimentary particles quickly reentered the flow field.
As the feeding time increases, the number of extrusive particles reaches a maximum when the nozzle location is at (Figure 6). The number of extrusive particles increases steadily from 0.25 s to 0.50 s but increases slowly when the nozzle position is at , indicating that the feeding and injecting masses are balanced. The number of extrusive particles tends to decline after 0.55 s and then declines sharply from 0.58 s to 0.62 s. This is the same trend observed in Figure 3 for a nozzle position of . The number of extrusive particles increases slightly after 0.70 s. However, this does not appear when the nozzle position is . The number of remaining particles in the four models is clearly different from each other; the number of remaining particles in the PTM is less than that for the other three models for a nozzle position of (Figure 7). After 0.5 s, the material feed stops and the number of extrusive particles linearly decreases in all of the models.
The simulation results show that the number of extrusive particles is maximised when the nozzle is positioned at , indicating that the injection performance is very good and the number of remaining particles is then minimised. However, the particlewall collisions are a significant issue in the delivery pipe. The linearity of particle trajectories is better and fewer particlewall collisions are observed compared to that for the other three models for a nozzle location of . However, the particle trajectories are complex and the number of remaining particles is maximised. Consequently, a nozzle location of is chosen for the gassolid injector testbed to verify the simulation results regarding the injection performance.
The experimental results show that nozzle location is the key factor that influences the injection performance, especially the injection mass. The injection mass increases first and then decreases as the nozzle position moves from left to right (or from A to I), indicating that there is an optimum nozzle position and thus verifying the simulation.
When the nozzle was in position A, the jet gas flow caused serious impacts, leading to significant noise and resulted in little mass being ejected. The material erupted in reverse from the hopper and blocked the mixing chamber during continuous feeding.
When the nozzle was in position D, the injection mass was maximised and the average maximum injection mass was 3.497 m^{3}/h with a peak injection mass of 3.66 m^{3}/h. These injection masses were determined to have met the requirements for the engineering application.
When the nozzle was in position I, the injection air did not mix with the material. Only a small amount of material was ejected out the mixing chamber and most of the material did not leave the hopper.
Due to the limitations of the experimental conditions, this paper only addresses the injection performance from the view of the nozzle position. There are other parameters that influence the injection properties of a gassolid injector that are not examined here, such as the geometric construction of the mixing chamber, the wall conditions, and the material properties. Future studies are expected to examine these aspects.
6. Conclusions
Based on the theoretical analysis, the conclusions obtained from the simulation and experimental studies regarding key factors that affect the injector performance are as follows:(1)The nozzle location is the key factor affecting the gassolid injector performance. There is an optimal nozzle location at which the number of extrusive particles is maximised and the number of particles remaining in the hopper is minimised.(2)The simulation results indicate that injection performance peaks when the nozzle is positioned at 30 mm. However, the particlewall and particleparticle collisions are more intense in the delivery pipe at this nozzle position. The particle trajectories are more linear and the particlewall collisions decrease for a nozzle position of 90 mm. The particle trajectories at this nozzle position are more complicated, and the greatest amount of remaining particles in the hopper is observed for this nozzle position.(3)The results of the experimental study indicate that the nozzle location has a clear influence on the injection performance. The injection performance improves remarkably in the optimum nozzle position of the gassolid injector when the structural parameters do not change.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This project is supported by National HighTech Research and Development Program of China (863 Program) (no. 2012AA062102), the Graduate Student Innovation Training Project in Jiangsu Province (nos. KYLX_1379 and CXLX13_936), and the Priority Academic Program Development of Jiangsu Higher Education Institutions.
References
 W.Y. Yu, R. Tian, and S.Y. Yan, “Numerical simulation on the flow structure and performance of the ejector,” Journal of Engineering Thermophysics, vol. 33, no. 11, pp. 1881–1883, 2012. View at: Google Scholar
 H. J. Chen, W. Lu, C. Cao, and L. Yang, “Applicability analysis on ejectors with cylindrical and conical mixing chambers,” Journal of Chemical Industry and Engineering, vol. 64, no. 6, pp. 2043–2049, 2013. View at: Publisher Site  Google Scholar
 L. J. Deng, Y. N. Liu, W. Wang, W. Ge, and J. H. Li, “A twofluid smoothed particle hydrodynamics (TFSPH) method for gassolid fluidization,” Chemical Engineering Science, vol. 99, pp. 89–101, 2013. View at: Publisher Site  Google Scholar
 H. Zhang, M. Liu, B. Wang, and X. Wang, “Dense gasparticle flow in vertical channel by multilattice trajectory model,” Science in China Series E: Technological Sciences, vol. 55, no. 2, pp. 542–554, 2012. View at: Publisher Site  Google Scholar
 K. Hong, Z. Shi, W. Wang, and J. Li, “A structuredependent multifluid model (SFM) for heterogeneous gassolid flow,” Chemical Engineering Science, vol. 99, pp. 191–202, 2013. View at: Publisher Site  Google Scholar
 K. Liu, J.S. Bai, and P. Li, “Application of particle trajectory model in 1D planar ejection,” Journal of Central South University of Technology, vol. 15, no. 1, pp. 149–154, 2008. View at: Publisher Site  Google Scholar
 M. A. I. Rasel, M. A. Taher, and H. D. Kim, “A study on the gassolid particle flows in a needlefree drug delivery device,” Journal of Thermal Science, vol. 22, no. 4, pp. 340–344, 2013. View at: Publisher Site  Google Scholar
 Z. G. Feng, E. E. Michaelides, and S. L. Mao, “Simulation of particlewall collisions in a viscous fluid using a resolved discrete particle method,” in Proceedings of the Asme Fluids Engineering Division Summer Conference, vol. 1, pp. 57–62, 2010. View at: Google Scholar
 H. Sakaino, “A semitransparencybased opticalflow method with a point trajectory model for particlelike video,” IEEE Transactions on Image Processing, vol. 21, no. 2, pp. 441–450, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 C. F. You, H. L. Zhao, Y. Cai, H. Y. Qi, and X. C. Xu, “Experimental investigation of interparticle collision rate in particulate flow,” International Journal of Multiphase Flow, vol. 30, no. 9, pp. 1121–1138, 2004. View at: Publisher Site  Google Scholar
 Z. Z. An, H. Zhou, Y. Yang, and X. L. Shi, “A study on particle trajectory of particle swarm optimization,” Applied Mechanics and Materials, vol. 433435, pp. 662–666, 2013. View at: Publisher Site  Google Scholar
 Y.Q. Xiong and M.Y. Zhang, “Experimental study on static pressure distribution properties of gassolid injector under pressure,” Journal of Combustion Science and Technology, vol. 10, no. 2, pp. 130–134, 2004. View at: Google Scholar
 Y.Q. Xiong, M.Y. Zhang, Z.L. Yuan, and D.M. Xue, “Threedimensional numerical simulation on gassolid twophase flows in gassolid injector,” Proceedings of the Chinese Society of Electrical Engineering, vol. 25, no. 20, pp. 77–82, 2005. View at: Google Scholar
 Y. Q. Xiong, M. Y. Zhang, and Z. L. Yuan, “Threedimensional numerical simulation method for gassolid injector,” Powder Technology, vol. 160, no. 3, pp. 180–189, 2005. View at: Publisher Site  Google Scholar
 Y. Q. Xiong, Z. L. Yuan, and M. Y. Zhang, “Threedimensional numerical simulation on conveying properties of gassolid injector under pressurization,” Journal of Chemical Industry and Engineering, vol. 10, no. 55, pp. 1638–1643, 2004. View at: Google Scholar
 A. Levy, “Twofluid approach for plug flow simulations in horizontal pneumatic conveying,” Powder Technology, vol. 112, no. 3, pp. 263–272, 2000. View at: Publisher Site  Google Scholar
 D. Z. Zhang and W. B. VanderHeyden, “The effects of mesoscale structures on the macroscopic momentum equations for twophase flows,” International Journal of Multiphase Flow, vol. 28, no. 5, pp. 805–822, 2002. View at: Publisher Site  Google Scholar
 J. D. Franklin and J. S. Lee, “Momentum/continuity coupling with large nonisotropic momentum source terms,” International Journal for Numerical Methods in Fluids, vol. 61, no. 9, pp. 946–969, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 A. Mehdizadeh, B. Firoozabadi, and S. A. Sherif, “Particle trajectory study in submerged flows with baffles using $\stackrel{}{{v}^{2}}f$ and $k\epsilon $ turbulence models,” Journal of Fluids Engineering, vol. 132, no. 5, Article ID 051105, 10 pages, 2010. View at: Publisher Site  Google Scholar
 V. Yakhot and S. A. Orszag, “Renormalization group analysis of turbulence. I. Basic theory,” Journal of Scientific Computing, vol. 1, no. 1, pp. 3–51, 1986. View at: Publisher Site  Google Scholar  MathSciNet
 N. Huber and M. Sommerfeld, “Modelling and numerical calculation of dilutephase pneumatic conveying in pipe systems,” Powder Technology, vol. 99, no. 1, pp. 90–101, 1998. View at: Publisher Site  Google Scholar
 Z. C. Zheng, X. Y. Zou, X. F. Yang, and H. Cheng, “Simulation of effects of the Saffman force and the Magnus force on sand saltation in turbulent flow,” in Proceedings of the ASME International Mechanical Engineering Congress and Exposition (IMECE '11), vol. 6, pp. 501–506, Denver, Colo, USA, November 2011. View at: Google Scholar
 Z. Zhu and J. Lin, “Pressure gradient force, saffman lift, and magnus lift on the fiberlike particle in fluid,” Journal of Dong Hua University, vol. 17, no. 2, pp. 23–27, 2000. View at: Google Scholar
 Q. Li, V. Rudolph, F. Y. Wang, and M. Horio, “A study of particle packing compression under fluid drag force by DEM simulations,” Developments in Chemical Engineering and Mineral Processing, vol. 13, no. 56, pp. 693–708, 2005. View at: 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. Forder, M. Thew, and D. Harrison, “A numerical investigation of solid particle erosion experienced within oilfield control valves,” Wear, vol. 216, no. 2, pp. 184–193, 1998. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2015 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.