MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi 10.1155/2019/7609748 7609748 Research Article Numerical Simulation of Air Entrainment Performance of a Foam Generator Used for Dust Control in Coal Mines https://orcid.org/0000-0002-0713-5835 Wang Hetang 1 2 https://orcid.org/0000-0001-6257-3769 Tang Yan 3 Chmielewski Adrian 1 State Key Laboratory of Coal Resources and Safe Mining China University of Mining and Technology Xuzhou 221116 China cumt.edu.cn 2 School of Safety Engineering China University of Mining and Technology Xuzhou 221116 China cumt.edu.cn 3 School of Environmental and Safety Engineering Changzhou University Changzhou 213164 China cczu.edu.cn 2019 482019 2019 27 02 2019 10 06 2019 01 07 2019 482019 2019 Copyright © 2019 Hetang Wang and Yan Tang. 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.

Mine dust is one of the most serious environmental hazards in the coal mining process. This paper introduces a numerical simulation of a novel foam generator used for dust control in coal mines. The amount of foam generated by this device significantly depends on the amount of air entrainment. Therefore, a computational fluid dynamics (CFD) method was used to study three influencing factors, namely, throat-nozzle distance, mixing throat length, and the contraction angle of the suction chamber. The predicted values by the CFD simulation proved to be in good agreement with the experimental data. The results revealed that the air entrainment reached its maximum when the ratio of throat-nozzle distance to mixing throat length was 2/3. The optimum values of the throat ratio (its length to diameter) and the contraction angle of the suction chamber were obtained at 20 and 5°, respectively. This research provides essential guidance in the geometric parameter design of the self-suction type foam generator, which has the advantage of negating the need for compressed-air pipelines and having high reliability, compared to traditional foam generators.

China Postdoctoral Science Foundation 2018T110574 Independent Research Project of Sate Key Laboratory of Coal Resources and Safe Mining, CUMT SKLCRSM19X004 Jiangsu Planned Projects for Postdoctoral Research Funds 1601234C
1. Introduction

Mine dust is one of the most serious hazards in coal mines; it severely pollutes the air in underground working environments [1, 2], leads to dust explosions and coal workers’ pneumoconiosis (CWP), and thus can result in casualties and economic loss in coal mines . Among the 18 major coal-mine accidents recorded worldwide, each resulting in over 300 deaths, 16 were caused by dust or dust-methane explosions . CWP is a global occupational disease existing in nearly all major coal-mining countries. In the USA, 8% of the workers who served in underground coal mines for over 25 years were diagnosed with CWP; 69,337 workers died because of CWP during the period 1970-2004 . Thus, there is an urgent need for the development and application of an efficient dust suppression technology.

There are various technologies for mine dust control, such as water sprays , air-curtains [13, 14], dust collectors , wetting agents (surfactants) [16, 17], and foams [18, 19]. Among them, foam technology plays an extremely significant role. Compared with water mist, foam has a larger surface area and stronger ability in dust adhesion and wetting [20, 21]. In particular, foam has a higher efficiency in suppressing respirable dust with small diameters [6, 18]. Hence, foam technology has been widely investigated and adopted in mine dust control. The production of foam is to the core of foam technologies. Almost all foam generators produce foam with compressed air; however this presents some common shortcomings: (a) extra compressed air pipelines need to be paved, taking up a lot of space and increasing the complexity of the system; (b) higher costs due to the large consumption of water and foaming agent; (c) backflow of pressure water which cannot be adjusted easily due to the pressure energy of compressed air being so strong. In order to overcome such drawbacks, a novel foam generator with air self-suction has been developed . In this generator, the pressurized water is the power and used as the driven fluid, and the air and foaming agent are automatically sucked into the induction chamber by the negative pressure formed by the water jet. The water, air, and foaming agent are mixed in the mixing chamber, producing the primary foam, which will be enhanced by the swirler in the foaming chamber.

This foam generator with air-suction was shown to achieve positive dust suppression effect in field applications . The suction performance of foaming agent was also investigated at length . However, it is evident that there is still a lack of in-depth studies on the relationship between the air entrainment performance and the geometric parameters of this foam generator. This has hindered the high volume and low cost production of dust-suppression foam, since the foam production quantities significantly depend on the air entrainment performance of the foam generator.

Therefore, in order to get the optimum geometric parameters with which the air entrainment of the foam generator would be maximum, three influencing factors of air entrainment were investigated, namely, throat-nozzle distance, mixing throat length, and contraction angle of the suction chamber (or induction chamber). Two dimensionless parameters were used in the study, the ratio R td (throat-nozzle distance to mixing throat length) and the ratio R ld (the throat’s length to its diameter). A computational fluid dynamics (CFD) method was adopted due to the simulation’s lower-cost, reliability, and ease of use .

2. Simulation Methods 2.1. Mathematical Model

The progress of jet suction in the foam generator can be classified into a gas-liquid two-phase flow (in this paper, the suction of the foaming agent was not studied as the amount of air entrainment was the main influencing factor contributing to foam generation). There are two methods for the numerical calculation of multiphase flows: the Euler-Lagrange method and the Euler-Euler method . In the Euler-Lagrange method, a particle is taken as the research object and the particle trajectory equations are solved in the Lagrangian framework, whilst in the Euler-Euler method, different phases are treated as interpenetrating continua which is more suitable for the jet suction process used in this paper. As the control equations for air and water are similar, water phase is taken as the example and its control equations are detailed as follows [29, 30].

The continuity equation: (1) t α 0 ρ 0 + · α 0 ρ 0 u 0 = m ˙ g 0 - m ˙ 0 g

where α 0 represents the fraction of water volume and u 0 is the velocity of the phase water, ρ 0 is the density of water, m ˙ g 0 characterizes the mass transfer from gas to water, and m ˙ 0 g characterizes the mass transfer from water to gas.

The conservation of momentum equation: (2) t α 0 ρ 0 u 0 + · α 0 ρ 0 u 0 u 0 = - α 0 p + · τ ̿ 0 + α 0 ρ 0 g + R g 0 + m ˙ g 0 u g 0 - m ˙ 0 g u 0 g + F 0 + F l i f t , 0 + F ν n , 0

where τ ̿ 0 is the stress-strain tensor of water, F 0 is an external body force, F l i f t , 0 is a lift force, F ν n , 0 is a virtual mass force, R g 0 is an interaction force between phases, and p is the pressure shared by two phases.

The energy conservation equation: (3) t α 0 ρ 0 h 0 + · α 0 ρ 0 u 0 h 0 = α 0 p 0 t + u 0 · τ ̿ 0 - · q 0 + S 0 + Q 0 g + m ˙ g 0 h g 0 - m ˙ 0 g h 0 g

where h 0 is the specific enthalpy of water phase, q 0 is the heat flux, S 0 is the source item that contains the enthalpy, Q 0 g is the strength of heat conduction between water and air, and h g 0 is the enthalpy value of the interphase.

The Realizable k-epsilon model is employed to model the turbulent behavior of flow in the generator.

2.2. Geometric Model

The geometric model was established using CAD and Gambit software. Subsequently, meshes were generated and the boundary conditions were set.

In the meshing process, the suction chamber was separated from the whole model and meshed with denser, unstructured “Tgrid” elements which have strong adaptability, whilst the other chambers were meshed with the “Hex/Wedge” elements. The total number of the grids varied as the structure parameters changed in different models. For instance, when the throat-nozzle distance was 12 mm, the mixing throat length was 360 mm and the contraction angle of the suction chamber was 20° and the total number of mesh elements was 219617, among which the number of the suction chamber was 74277.

Figure 1 shows the boundary conditions of the model. The pressurized water jet can be considered as incompressible fluid flow. Thus the velocity inlet was taken as the boundary condition for the flow of water when the volumetric flow rate of water at the nozzle was fixed. The air inlet was set as the pressure inlet whilst the ejector outlet was considered the pressure outlet. Their gauge pressures were considered to be 0 Pa since they were both open to atmosphere. The solution was initialized by taking the ejector geometry as being completely filled with air.

Schematic diagram of boundary conditions setting:   α contraction angle of suction chamber;   L m length of throat;   d m diameter of throat;   L tn throat-nozzle distance.

2.3. Solution Method

The calculations were conducted on an AMD6272 server (128 CPU core, 2.1 GHz dominant frequency and 256G internal storage) with FLUENT 14.5 software. The simulation process was divided into three steps: (a) simulation of air entrainment under different throat-nozzle distances; (b) simulation of air entrainment under different lengths of mixing throat; (c) simulation of air entrainment under different contraction angles of the suction chamber.

3. Results and Discussion

The foam generator studied in this paper is essentially a special jet pump which can safely and reliably mix several kinds of material. Based on the design theory of a jet pump, throat-nozzle distance, mixing throat length, and contraction angle are important structural parameters affecting the suction performance . Thus, this paper chose these parameters as a measure of the performance of the foam generator.

3.1. Validation of Model

In order to validate the model and obtain high accuracy of the simulation results, related experiments were undertaken to compare the simulation results with measured results by using the experimental setup in . Figure 2 shows the comparison between the simulation values and the experimental data for air entrainment under different R td values. The throat length and the contraction angle were set at 360 mm and 20°, respectively (with a jet velocity of 89.75 m/s). It can be seen from Figure 2 that the simulation and experimental values showed similar variation trends and good agreement overall. Attributing reasons for errors between the experimental and simulation data include the following: (a) the control equations used in FLUENT software are applicable to the general situation but may not be totally accurate to the specific problem presented in this paper; (b) the sensors used in the experiments may have inaccuracies/sensitivities and their locations may be not have been the most appropriate; (c) the initial conditions (such as the turbulence intensity) are time-averaged in the simulation but are distributed unevenly in reality. Taking all of these aspects into consideration, the simulation results can be considered credible.

Simulation and experimental data of air entrainment under different R td.

3.2. The Effect of Throat-Nozzle Distance on Air Entrainment

The throat length and the contraction angle of the suction chamber were set at 360 mm and 20°, respectively. The throat-nozzle distance was viewed as a variable and the simulations were conducted under the values of 1.0 mm, 3.0 mm, 6.0 mm, 9.0 mm, 12.0 mm, and 15.0 mm. As the throat diameter d m was 18 mm, the ratio ( R t d ) of throat-nozzle distance to mixing throat length was thus 1/18, 1/6, 1/3, 1/2, 2/3, and 5/6, respectively. There were 8 different values of jet velocity v n under each ratio, 52.97 m/s, 55.96 m/s, 65.55 m/s, 75.13 m/s, 83.94 m/s, 89.75 m/s, 96.51 m/s, and 104.22 m/s.

Figure 3 shows the velocity nephogram of the water and air phases on section Z=0mm when v n is 104.22 m/s and R td is 2/3. As water ejects from the nozzle at high velocity, the outside air is induced into the suction chamber and moves together with the water.

Velocity nephogram of (a) water phase and (b) air phase on section Z=0 mm when v n is 104.22 m/s and R td is 2/3.

Figure 4 shows the air entrainment Q g under different R td values. For each v n , the air entrainment was highest when the R td was 2/3 and decreased to its lowest point when the R td was 1/18. As the R td transitioned from 1/18 to 1/3, the air entrainment showed an increase and reached a high point when the R td was 1/3. With a jet velocity of 83.94m/s and R td at 1/3, Q g was 1.49 times greater than that when R td was at 1/18. As R td increased from 1/3 to 2/3, air entrainment first decreased, then raised, and reached the peak when R td was 2/3. Furthermore, the air entrainment showed a decreasing trend as R td increased from 2/3 to 5/6. Thus for this device, the optimal value for R td was found to be 2/3; hence, the optimal throat-nozzle distance was 12 mm when the throat length was 18 mm.

Change of Q g with R td when the L m and α are 360 mm and 20°, respectively.

3.3. The Effect of the Length of Throat on Air Entrainment

The throat-nozzle distance and contraction angle of the suction chamber were set to a value of 12mm and 20°, respectively, and the simulations were conducted under different mixing throat lengths of the 540mm, 450 mm, 360 mm, 270 mm, and 180 mm. The ratios for R l d (the length to diameter of throat) were 35, 30, 25, 20, 15, and 10 respectively. Furthermore, there were 8 values of jet velocity v n tested under each ratio: 47.79 m/s, 50.77 m/s, 56.75 m/s, 65.71 m/s, 77.18 m/s, 85.83 m/s, 94.00 m/s, and 102.18 m/s.

Figure 5 shows the changes in air entrainment Q g under 6 different values of R ld. As illustrated in the figure, the air entrainment reaches its highest point when R ld is 20 and falls to its lowest point when the R ld is 35. As R ld increases from 10 to 20, the air entrainment shows a corresponding increase and reaches a peak when the R td is 20. As R ld increases from 20 to 30, air entrainment experiences a slow decrease, and then it falls significantly when R ld is 35. With a jet velocity of 65.71m/s, the minimum value of air entrainment decreases by 41% as compared to peak values. Thus for this device, the optimal value of R ld was found to be 20, with the optimal throat length of 360 mm and throat diameter of 18 mm.

Change of Q g with R ld when L tn and α are 12mm and 20°, respectively.

3.4. The Effect of the Contraction Angle on Air Entrainment

The throat-nozzle distance and throat length were set with values of 12mm and 360mm, respectively, and simulations were conducted under different contraction angles α of 3°, 5°, 10°, 15°, 20°, 25°, 45°, and 60°. There were 8 values for the jet velocity v n under each α: 45.00 m/s, 55.96 m/s, 65.55 m/s, 75.13 m/s, 83.94 m/s, 89.75 m/s, 96.51 m/s, and 104.22 m/s.

Figure 6 shows the air entrainment Q g under 8 different contraction angles of the suction chamber. The air entrainment reached its highest point when α was at 5° and fell to its lowest point when α was 60°. When α increased from 3° to 5°, the air entrainment showed a slight increase. As α increases from 5° to 20°, the air entrainment decreased. The air entrainment remained relatively steady when α varied from 20° to 45°, after which it dropped dramatically when α was 60°. With respect to the jet velocity of 89.75m/s, the air entrainment decreased by 39.4% and was only 60% of maximum values.

Change of Q g with α when L tn and L m are 12mm and 360mm, respectively.

Since the secondary fluid air is induced by the negative pressure generated by the jet, the effect of α can be explained by the driving force of negative pressure in the chamber. Figure 7 shows the pressure cloud chart of the suction chamber on section Y=0 mm and Z=0 mm. The jet velocity was 77.18 m/s, the contraction angle of the suction chamber was 20°, throat length was 450 mm, and the throat-nozzle distance was 12mm. The high negative pressure mainly existed in the area around the nozzle exit and the entrance of mixing chamber.

Pressure distribution nephogram on section (a)Y=0mm and (b) Z=0mm.

Figure 8 describes the static pressure values along the central axis which begin from the nozzle exit with a jet velocity of 96.51 m/s and α value of 20°. The static pressure first decreased sharply and then increased rapidly until it is about 5 mm from the nozzle exit and reached the minimum point at a value of -6200 Pa. Subsequently, the static pressure went through a fluctuation phase which was not as intense. This result was consistent with the phenomenon of shock train in ejectors which is driven by the combined effect of interphase drag and turbulent eddies . In this case, the vacuum degree, which is the absolute value of negative pressure, was used to describe the strength of the negative pressure; correspondingly, the maximum vacuum degree is 6200 Pa.

Static pressure on central X axis.

Figure 9 indicates the maximum vacuum degrees along the central axis under different contraction angles. The maximum vacuum degree was highest when α was 5° and fell to its lowest point when α was 60°. As α increased from 3° to 5°, the maximum vacuum degree increased and reached a peak when α was 5°. The maximum vacuum degree exhibited a tendency to decrease as α increased from 5° to 60°. Therefore, the optimal value of α was 5° which illustrated a corresponding larger vacuum degree.

Change of maximum vacuum degree with α.

In conjunction with Figure 6, a synchronous variation can be found between the vacuum degree and the air entrainment. This phenomenon can be explained by the following equation which indicates that for the determined parameters, the air entrainment Q g increases with the vacuum degree in the suction chamber : (4) Q g = π d a 2 4 × 2 p a t m - p m / ρ g μ 1 + μ 2 + k d a 4 / d m 4 where d m is the throat diameter; d a is the diameter of air suction hole; ρ g is air density, kg/m3; p a t m and p m denote the atmospheric pressure and air pressure, respectively, at the throat inlet, MPa; μ 1 and μ 2 denote the local resistance coefficient at the inlet and outlet of the suction hole, respectively; k is the kinetic energy correction coefficient of the throat section.

3.5. Field Application

With adoption of the parameters mentioned above, this foam generator with air self-suction can have good foam-production ability and high dust-suppression efficiency in field applications and have been utilized in several coal mines in the Huaibei mining area such as Zhuxianzhuang Coal Mine [22, 23].

4. Conclusions

The air entrainment performance of the foam generator presented in this paper was shown to be significantly affected by the ratio of the throat-nozzle distance to the mixing throat length (R td), the length-diameter ratio of the mixing throat (R ld), and the contraction angle of the suction chamber (α). The following conclusions are presented.

The foam generator’s air entrainment reaches a peak when the R td is 1/3 and 2/3. The minimum and maximum values for air entrainment are obtained when R td is 1/18 and 2/3, respectively. With the optimal value of R td at 2/3 in this paper, the optimal throat-nozzle distance is 12 mm when the throat length is 18 mm.

Raising the length-diameter ratio, R ld, increases the amount of air suction initially but this begins to decrease after surpassing a certain value. The maximum and minimum values of air entrainment are obtained when R ld is 20 and 35, respectively. Hence, the optimal value of R ld is 20, with the optimal throat length at 360 mm, when the diameter of throat is 18 mm.

As the contraction angle of the suction chamber (α) increases from 3° to 60°, the air entrainment shows a trend of increasing initially and subsequently decreasing after surpassing a specific value. A peak value is obtained when α is 5° and this falls significantly when α is 60°. Thus, the optimal value of α for the device is 5°.

The optimum parameters have a corresponding larger vacuum degree. Thus, there is a positive relationship between the vacuum degree and the air entrainment.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the China Postdoctoral Science Foundation funded project (2018T110574), the Independent Research Project of Sate Key Laboratory of Coal Resources and Safe Mining, CUMT (SKLCRSM19X004), and the Jiangsu Planned Projects for Postdoctoral Research Funds (1601234C).