Numerical Simulation of the Tar Mist and Dust Movement Process in a Low-Temperature Dry Distillation Furnace

Zhang, Qiuli; Hui, Xiangrong; Yan, Long; Luo, Min; Feng, Wenru; Zhou, Jun; Lan, Xinzhe

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

Journal of Chemistry

On this page

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

Research Article | Open Access

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

Numerical Simulation of the Tar Mist and Dust Movement Process in a Low-Temperature Dry Distillation Furnace

Qiuli Zhang,¹Xiangrong Hui,¹Long Yan,²Min Luo,¹Wenru Feng,¹Jun Zhou,¹and Xinzhe Lan³

Academic Editor: Marco Anni

Received28 Nov 2019

Revised30 Dec 2019

Accepted03 Jan 2020

Published02 Mar 2020

Abstract

In the low-temperature dry distillation of low-rank coal, the important liquid product of coal tar is produced, but its quality and utilization rate are degraded by entrained dust. The movement of coal tar and dust in the furnace is a key factor in causing particles such as dust to mix with coal tar. Therefore, the Euler–Lagrangian method is used to simulate the two-phase motion process of gas, tar, and dust in a furnace. By considering the effects of tar particle size, dust particle size, gas velocity, tar density, and dust density, the motion process mechanism is revealed, enabling the dust content in coal tar to be reduced and the quality improved. The results indicate that tar particles with sizes less than 0.20 mm can be removed from the furnace by gas, and the smaller the particle size is, the shorter the time required for removal. Dust particles greater than 0.18 mm in size cannot be completely removed from the furnace. As the gas velocity increases, the time required for complete removal of the tar mist and dust gradually decreases. When the speed is 0.70 m/s, all tar mist is removed, although some particles remain. Tar mist with a density of more than 900 kg/m³ can be extensively removed, but dust with a density of more than 1400 kg/m³ is difficult to remove and remains in the furnace. Finally, particle size distribution experiments in the product were conducted to verify the accuracy of the numerical simulation.

1. Introduction

China’s coal resources are relatively abundant, and its reserves account for approximately 90% of the country’s mineral resource reserves and 95% of its fossil energy, with significant resource potential [1]. Moreover, the country’s coal resources include coal types from lignite to anthracite in various coalification stages, especially low-rank reserves that account for more than 40% of coal reserves and occupy a vital position in China’s coal resources [2, 3]. Coal accounts for approximately 72% of the country’s total energy consumption, and China’s coal-based energy structure is expected to remain stable for a long time [4].

China’s northern Shaanxi region is rich in coal resources, with proven reserves of 146 billion tons. These resources are dominated by low-rank long-flame coal, which has a high oil content (average 6% to 7%) and has the characteristics of “three lows and one high” (low ash, low sulfur, low phosphorus, and high calorific value). This coal is an ideal raw material for low-temperature retorting [5]. Coal dry distillation is a necessary process for coal conversion processes such as coal liquefaction [6], gasification [7], combustion, and coking [8]. In addition, coal dry distillation is an efficient, low-pollution coal classification and conversion technology [9–11]; the gas produced during the dry distillation process can be used as a fuel and industrial raw material [12], and semicoke can be used in the metallurgical industry, syngas production, and other fields [13]. The phenols and hydrocarbons rich in coal tar become important raw chemical materials after separation and purification [14]; moreover, this material can partially replace oil and optimize the energy allocation of the petrochemical industry. China’s crude oil storage is small, and China has long relied on imports. Therefore, improving the production and quality of coal tar during low-temperature dry distillation can partially alleviate China’s petroleum energy crisis and substantially augment China’s reserve oil.

In the low-temperature dry distillation process, the furnace type adopted in the northern part of Shaanxi is mainly an SJ-type internal heat low-temperature dry distillation furnace. This furnace type is a new process developed by a company based on the Lurgi-Spilgas three-stage furnace. It utilizes the advantages of related furnace types at home and abroad, combines local actual production experience, and designs according to the characteristics of low-rank coal. A moving bed reactor is used, and the burned gas is used as a gas heat carrier to heat coal in the bed to facilitate the pyrolysis reactions of the coal. The process is simple, and a high yield of tar can be obtained, making this process one of the preferred techniques for low-rank coal cryogenic dry distillation in China. However, the furnace gas produced by the process is cooled during the temperature rise, and part of the tar forms liquid tar mist. Moreover, dust such as pulverized coal and semicoke entrained in the furnace gas flow are also incorporated, which makes the separation and recovery of coal tar difficult, resulting in poor quality and greatly reducing the potential value of the tar. Therefore, the separation of coal tar and impurities has become a key problem. Based on actual observations of phase flow and the complexity of measuring industrial data, it is very difficult to solve this problem by experimental methods. However, with the rapid development of computers, the computational fluid dynamics (CFD) method has gradually become effective for studying complex multiphase flow processes in reactors. Many researchers have implemented the coupling of complex reaction kinetic models into CFD models for flow-phase simulation [15–18]. This work provides a technical basis for the simulation calculation of tar mist and dust movement processes by the CFD method.

Currently, the calculation methods of fluid-solid coupling are divided into three categories: the Euler-Eulerian method, the Euler–Lagrangian method, and the Lagrangian-Lagrangian method. Because the Euler–Lagrangian method can accurately reflect the airflow-dust coupling characteristics [19–22], this method is increasingly widely used in the simulation of fluid-solid two-phase flow. Torañod et al. [23] studied the migration behavior of airflow and air-dust flow through a 14.7-m² cross section on a fully mechanized excavation surface under strong exhaust ventilation conditions by the Euler–Lagrangian method. Kurnia et al. [24] used the Euler–Lagrangian method to study the diffusion and management of dust on the excavation surface. Patankar and Joseph [25] used the Lagrangian method to describe the movement of dust particles.

In previous studies, numerical simulations have been widely used to study the dust diffusion pollution behavior in coal, and some useful results have been achieved. For example, Lu et al. [26] used CFD to evaluate various possible scenarios for underground mine configurations. Zhou et al. [27] combined CFD simulation and onsite measurement to study the diffusion behavior of inhalable dust. Cai et al. [28] studied the effect of airflow on the diffusion mode of coal dust particles in a fully mechanized mining face and found the optimal airflow.

In the past 5 years, the model of pulverized coal-fired boilers has improved significantly. Modeling furnaces for coal-fired boilers is a complex task. Any model developed to predict critical performance parameters, such as furnace exit gas temperature and pollutant emission levels, will need to capture the multifarious nature of combustion in the boiler furnace. Despite the large amount of computationally expensive CFD codes, researchers have also developed simpler but effective reactor network models. Sankar et al. [29] reviewed the research status of calculation models for pulverized coal boilers from three perspectives: the development of coal combustion models, radiation models, and overall furnace models.

Nevertheless, there are very few simulation studies on the movement process of focusing tar oil mist and dust in a low-temperature dry distillation furnace. Through numerical simulation of the gas-solid two-phase flow state, the two-phase motion mechanism of the coal retorting process is revealed. The physical parameters (e.g., particle size of tar and dust, gas velocity, and density of tar and dust) can be optimized so that the tar can be quickly removed by the gas, while the dust, such as pulverized coal and semicoke, remains in the furnace to continue pyrolysis. This condition can reduce the waste slag in the coal tar, improve the quality of the coal tar, partially alleviate China’s petroleum energy crisis, and substantially augment China’s reserve oil. In addition, the raw materials can be brought closer to the theoretically superior value before production, thereby improving the efficiency of pyrolysis production and reducing energy consumption and costs.

2. Physically Simplified Model of Airflow-Particle Motion

In an SJ-type internal heat type low-temperature dry distillation furnace, the tar that is dry distillated by coal exists in a gaseous form in the high-temperature region. The tar and gas move upward in the space between the coal particles, and when it reaches the low-temperature region (below the boiling temperature of the tar component), the tar condenses into a tar mist. The main components and boiling points of tar are shown in Table 1.

The tar mist is brought out of the furnace by the rapidly rising gas. At the same time, the coal generates dust when it enters the furnace, and some of this dust is removed from the furnace by the gas. By simplifying the model, the process of tar mist and dust removal by gas is considered a pneumatic conveying process. Showing a characteristic length of L (unit is m) and an equivalent diameter of d (unit is m) of the microelement tube as the research object, a schematic diagram of the tar mist and dust movement is shown in Figure 1, where u is the apparent velocity in m/s and n is the void fraction. Due to the presence of scattered solid particles in the tube, which increases the difficulty of establishing the geometric model, we arrange the scattered solid particles in Figure 1(a) in an orderly manner to obtain Figure 1(b). When the arrangement is tight enough, we can ignore its existence in the tube and then express it by a void fraction. The movement process of tar mist and dust in the tube is shown in Figure 1(c).

3. Mathematical Model of Airflow-Particle Motion

In this study, gas is considered as the fluid phase, and tar mist and dust are considered as the particle phase. It is also assumed that the flow of gas in the model conforms to the continuous medium hypothesis, conservation of energy law, and Newton’s second law of motion [30]. On the basis of the standard k-ε model, we constructed the closed control equations to calculate the turbulence. This model only considered momentum transfer and ignored thermal conduction. Then, dust particles were treated as discrete media whose movements were calculated based on the Lagrangian method [27, 31, 32]. The associated mathematical models are as follows.

3.1. Mathematical Model of Air Movement

The continuity equation [33, 34]:

The Navier–Stokes equation:where ρ is the fluid density; t is time; x is the coordinate; u is velocity; p is the turbulent pressure; and μ_t is the coefficient of turbulent viscosity.

The energy equation [35]:where h is the specific enthalpy, is the ratio of effective viscosity and the Prandtl number, is the dissipation function, and is the source of enthalpy generated on a volumetric source basis. can further be calculated as

Standard k-ε equation:

The k-equation [36, 37]:

The ε-equation:

In these equations, k is the fluid’s turbulent kinetic energy; ε is the viscous dissipation induced by a vortex; G_k is the turbulent kinetic energy; and C_1ε, C_2ε, and C_μ are the empirical constants with values of 1.44, 1.92, and 0.09, respectively.

3.2. Force Analysis of the Fluid Phase

As shown in Figure 1, in the microelement length L, the intracavity void ratio is n, and it is assumed that the in-tube free cross section fraction is equal to n.

The fluid phase is subjected to a frictional resistance F_{f, g}, the force of the pressure drop gradient on the fluid phase F_{Δp, g}, the self-gravity, and the particle drag F_{D, g}. All other forces are negligible:where Δp_f is the pressure difference generated by the fluid motion, and the unit is Pa.

3.3. Mathematical Model of Tar Mist and Dust Movement

In the Lagrangian coordinate system, the motion trajectory of discrete-phase particles was calculated by integrating the differential equation of the interaction among particles [30, 38–40].

Time-averaged particle velocity is calculated as the sum of particle convective velocity and diffusive velocity [41]. The convective velocity of particles is

The diffusive velocity:

The additional forces in the particle force balance, mainly including the pressure gradient force, virtual mass force, Saffman lift force, and Magnus lift force, are important under special circumstances [42].

The pressure gradient force:

The virtual mass force:where is the virtual mass factor with a fixed value of 0.5.

The Saffman lift force [43, 44]:where K = 2.594; d_ij, d_lk, and d_kl refer to the deformation tensor; and is the kinematic viscosity.

The Magnus lift force:where is the rotational angular velocity.

4. Geometry Model Establishment and Boundary Condition Setting

Figure 2 is a schematic diagram of the structure of a 75,000-ton SJ-type internal heat low-temperature dry distillation furnace used by a coal chemical company in Shaanxi. The furnace adopts a cavity design. From top to bottom is a drying section, a dry distillation section, and a cooling section. Coal blocks gradually descend from the coal box through the cavity. The gas and air are mixed into the furnace from below, and the gas and oxygen burn to form a high-temperature heat carrier to heat the coal seam layer by layer. The tar produced by pyrolysis is removed from the riser by rising gas. According to this process, a geometrical model of the riser is constructed, which is divided into a square grid with a total of 15,320 cells, including 40 cells on the X-axis and 383 cells on the Y-axis. As shown in Figure 3, the mesh with quality higher than 0.3 accounts for more than 97% of all meshes, indicating that the mesh model is reasonable.

In practice, it is a continuous process for tar mist and dust to be removed from the furnace by the rising gas, but it is difficult to simulate this process. Therefore, in the numerical simulation, it is assumed that the lower part of the riser is filled with particles with a volume fraction of 30%, and the gas is then blown upward from the bottom of the tube at a certain speed to observe the rise of particles. Figure 4 shows the state of the particles in the tube filling the bottom (red portion) before iteration after setting the particle volume fraction to 30% in ANSYS Fluent.

Because the tar mist particles only exist in the riser, once they are removed from the furnace, they will agglomerate to form liquid tar. It is difficult to directly measure during the production process, so a large number of simulation experiments were performed in the early stage to determine the particle size scope. Finally, the research scope is determined as follows: the tar mist particle size is 0.1 mm–1.0 mm, and the dust particle size is 0.05 mm–0.20 mm. According to the density scope of tar and dust provided by the coal chemical company, numerous simulation experiments were carried out in the early stage. The final determination scope is as follows: tar mist density 900 kg/m³–1100 kg/m³ and dust density 1100 kg/m³–1500 kg/m³. The scope of gas velocity provided by the coal chemical company is 0.675–0.723 m/s. In this study, 0.68, 0.69, 0.70, 0.71, and 0.72 m/s were selected as the gas velocities. Table 2 lists the boundary conditions for the geometric model.

5. Results and Discussion

5.1. Effect of the Tar Mist Particle Size

When the gas velocity is 0.70 m/s and the density of the tar mist is 900 kg/m³, the particle size of the tar mist is changed to study the two-phase motion. Figure 5 shows the volume fraction of the tar mist at different times.

(a)

(b)

(c)

(d)

(e)

It can be seen from Figures 5(a) and 5(b) that the distribution area of the particles in the furnace increases rapidly with time, and the entire area is filled within 1.8 s. Then, the distribution of the particles in the working area is continuously reduced until the particles are completely removed by the gas. When the particle size is 0.10 mm and the time is 6.1 s, only particles at a volume fraction of 0.03 are present in the furnace, and the remaining positions are filled with the gas phase, indicating that 97% of the tar particles are carried out by the gas. When the particle size is 0.20 mm and the time is 0.9 s, the particles are brought to a height of approximately 4/6 in the tube; when the time reaches 22.0 s, the particles are all carried out of the working area by the gas.

It can be seen from Figures 5(c) and 5(d) that within 1.8 s, the particle distribution area begins to decrease significantly; from 1.8 s to 120 s, 95% of the particles are below a height of 1/2, and the distribution area and state are similar. After 120 s, this state remains unchanged, indicating that the particles are difficult to remove. In Figure 5(e), the particle distribution area is always at a height of 1/6, and the area does not change very much with time, indicating that the particles are also hard to remove in this case.

A comparison of Figures 5(c)–5(e) indicates that the volume fraction of the particles at the bottom of the furnace is unevenly distributed, and the particles are fluidized. Moreover, as the particle size increases, the particle distribution area gradually decreases, as does the degree of fluidization of the particles. A comparison of the subfigures of Figure 5 reveals that as the particle size increases, the distribution area of the particles decreases at every time point, the volume at the bottom gradually increases, and the degree of fluidization gradually decreases. At the same gas velocity, the larger the particle size, the more difficult removal becomes.

5.2. Effect of the Dust Particle Size

When the gas velocity is 0.70 m/s and the dust density is 1100 kg/m³, the particle size of the dust is changed to study the two-phase motion. Figure 6 shows the volume fraction of dust at different times.

(a)

(b)

(c)

(d)

(e)

Figures 6(a)–6(d) show that the distribution area of the particles in the furnace increases rapidly with time, and the entire area is filled within 1.6 s. Then, the distribution of the particles in the working area is continuously reduced. Finally, only particles at a volume fraction of 0.015 remain, the rest having been removed by the gas. As the particle size of the dust increases, the time required for dust removal gradually increases. In Figure 6(e), when the particle size of the dust is 0.20 mm, the distribution area and the volume fraction at 23.0 and 120 s do not change much, and the particles are deemed difficult to remove.

In Figure 6, the particle distributions of 0.2 s and 0.4 s are the same, and only the proportions of the different volume fractions are different. Particles (with a volume fraction of 0.30) are gradually depleted, and the proportion of remaining particles gradually increases, producing a volume fraction between 0.27 and 0.03. At 0.8 s, the distribution of the upper part is unchanged, but the volume fraction is gradually reduced, the height of the particles is slightly lowered, and the volume fraction at the bottom is gradually increased. These effects are caused by the increase in particle size. At 1.6 s, the particles fill the entire working area; however, as the particle size increases, the upper distribution area decreases, and particles gradually gather at the bottom with a volume fraction of 0.30.

5.3. Experimental Verification

The above simulation results demonstrate that the influence of the particle size on the rise of oil mist and dust is important. Zhang et al. [42] used the AKFC-92A-type mine dust sampler to perform dust sampling 30 minutes after coal mining to verify the accuracy of the simulation results. To verify the accuracy of the simulation results, Cai et al. [28] set up multiple measurement points onsite and collected dust particle samples using a dust sampler to measure the dust concentration. However, due to the complicated structure of the coal low-temperature dry distillation furnace and the uninterrupted production process, it is difficult for the measuring instrument to be directly introduced into the rising pipe for measurement. The coal tar produced on-site by this type of a low-temperature dry distillation furnace in a factory is dissolved with an organic solvent such as toluene, acetone, or quinoline, followed by suction filtration and washing to obtain dust particles, and the particle size distribution of the dust is measured by the OMEC Topsizer Plus laser particle size analyzer. The obtained particle size distribution is shown in Table 3 and Figure 7.

As shown in Figure 7, the dust particles less than 0.18 mm accounted for 90%, and the largest particle size was not more than 0.339 mm, indicating that the dust particles less than 0.18 mm were easily taken out by the gas and affected the quality of the tar. It can be seen from Figure 6(e) that the 0.20 mm dust cannot be completely taken out, but it is not solidified like tar mist deposited on the bottom, or a small part will be taken out, which accounts for about 4.68% (Table 3). In summary, it was found that the particle size distribution could achieve the expected simulation effect in practice. By comparing the particle size change in experiments and simulations, the effects on the dust rising process were investigated, and the accuracy of the mathematical models used in the simulation process and the reliability of the simulation results were verified. The rationality of particle size parameter selection was further proven.

5.4. Effect of the Gas Velocity on the Tar Mist and Dust Movement

Two-phase motion is studied by changing gas velocity. The volume fractions of tar mist and dust are shown in Figures 8 and 9, respectively.

(a)

(b)

(c)

(d)

(e)

(a)

(b)

(c)

(d)

(e)

From Figures 8 and 9, the distribution area of the particles in the furnace increases rapidly with time and fills almost the entire area within 2.0 s. At 4.0 s, the proportion of particles (with a volume fraction of 0.20 to 0.40) continuously decreases, the particles gradually move upward, and the bottom is gradually occupied by a volume fraction of less than 0.12. Finally, after approximately 20 s in the furnace, there is only a small amount of tar mist and dust particles with a volume fraction of 0.04. The simulation results reveal that as the gas velocity increases, the time over which the tar mist and dust are completely removed gradually decreases. The reason is that the gravity G of the mass per unit particle is constant in the case where the particle diameter and the density are the same. However, due to the difference in gas velocities, the resistance between the gas-solid phase and the additional force in the particle force balance may vary from equations (11)–(14). Following force balancing, an acceleration in the upward direction is finally obtained so that the time required for the belt to remove material decreases as the velocity increases.

A comparison of Figures 8 and 9 reveals that there is no change in the particle distribution shape, height, area, or volume fraction at 0.5 s and 1.0 s. The reason may be that the increase in the velocity is not very large; the particle phase is initially at the bottom of the furnace, and the acceleration in the upward direction needs to gradually increase with time. After 2.0 s, the acceleration starts to cause a visible change; the volume fraction at the bottom decreases, the particles spread upward, and increasingly more particles are removed, gradually reducing the total time required for complete removal. Taken together with the simulation results at the same velocity, these results reveal that the time required for the dust to be removed is more than that for the tar mist to be removed. Because the difference in particle size is not very large, the difference in density plays a large role. Density determines the gravity force, and the density of tar mist is less than that of dust, meaning that the gravity force on the tar mist is less than that of dust of the same mass. At the same velocity, the other additional forces are unchanged, so the time required for the tar mist to be carried out is less than that for the dust.

Since the tar mist and dust are mixed into the furnace, the gas velocity is the same for both phases. The main purpose of this study was to separate the two phases, letting the tar mist be carried out and allowing as much dust as possible to flow back into the furnace for further dry distillation, so it is important to determine the appropriate corresponding velocity. Figure 10 presents the time at which the particles are removed in Figures 8 and 9. The time difference between 0.68, 0.69, 0.71, and 0.72 m/s is small, and all values lie within a range of 0.7 s. At 0.70 m/s, the time difference is the largest, which is 1.7 s, which means that there is still more dust in the furnace when the tar mist is completely removed.

It can be seen from 5.1 and 5.2 that the above selected tar mist and dust particle size are the maximum values that can be removed by the gas. Select the two-phase particle with the smallest particle size when not removed and study the two-phase motion at the maximum velocity, as shown in Figures 11 and 12.

(a)

(b)

(a)

(b)

Figure 11 shows that the distribution area of particles in the furnace increases rapidly over time. At 1.8 s, the particle fills the area to 5/6, and the distribution area reaches its maximum. At 9.0 s, the particle distribution area is reduced to 4/6. At 120 s, it is reduced to 1/2, and the particles are always left in the furnace. A comparison of the simulation results at the same time reveals no change at 0.3 s. With increasing velocity at 0.6 s, the tar particles with a volume fraction of 0.35–0.45 aggregate in a higher location, and the particles with a bottom volume fraction of 0.20–0.30 are depleted. At 0.9 s, as the velocity increases, tar particles are concentrated with a volume fraction of 0.30 to 0.40 and 0.35 to 0.45 at the upper and bottom portions of the distribution area, respectively. At 9.0 s and 120 s, particles at a volume fraction of 0.35 to 0.50 disappear with increasing velocity, showing that although most of the particles remain in the furnace, some of them are still removed. It can be seen from Figure 12 that the distribution area of the particles in the furnace increases rapidly over time, and the entire area is filled in 1.6 s, with the distribution of particles at a volume fraction of 0.18–0.30 continuously decreasing in the working area. However, there are still some particles left in the furnace. A comparison of the simulation results at the same time reveals no change at 0.2 s and 0.4 s; at 0.8 s, particles (with a volume fraction of 0.21 to 0.27) are collected at a height of 4/6, and there is no significant difference in the remaining time except that the distribution area is expanded from 1/2 to 5/6.

Figures 8, 9, 11, and 12 verify that the gas velocity can change the length of time required for the two phases to be carried out, but the particle size of the particles is a key factor in determining whether removal is possible. The particle size of the tar mist should be less than 0.20 mm, and the particle size of the dust should be greater than 0.18 mm to achieve the aims of this study.

5.5. Effect of the Tar Mist and Dust Density

The density of tar mist and dust was changed, and the two-phase motion was studied. The motion simulation results are shown in Figures 13 and 14.

(a)

(b)

(c)

(d)

(e)

(a)

(b)

(c)

(d)

(e)

Figure 13 shows that the distribution area of the particles in the furnace increases rapidly over time. At 2.7 s, the particles fill almost the entire area, and the distribution of particles gradually decreases. Finally, at 33.5 s, there is only a small amount of tar particles in the furnace, with a volume fraction of 0.03. As the density increases, the time required to remove the particles increases because of the influence of gravity. There is little difference between the results at 0.3 s and 0.9 s. The distribution area at 2.7 s is at a height of 5/6; however, the particles (with a volume fraction of 0.06–0.12) are enriched with increasing density, and the particles are slowly transported to the bottom, with a volume fraction of 0.30. At 8.1 s, only particles at a volume fraction of 0.12 or less remain in the furnace. As the density increases, the distribution area first increases and then decreases. Specifically, only a few particles (with a volume fraction of 0.03 to 0.09) are left at ρ = 900 kg/m³; at ρ = 950 kg/m³ and 1000 kg/m³, the density increases, and it is difficult to remove many particles at the same time, resulting in an increase in the distribution area. When ρ = 1050 kg/m³ and 1100 kg/m³, the volume fraction of the particles increases significantly, and the particles collectively move downward, resulting in a decrease in the distribution area. As shown in Figure 14, the distribution of particles in the furnace increases rapidly with time. The area is almost filled at 2.7 s, and the distribution of particles in the area is then gradually reduced. The final density is between 1100 kg/m³ and 1400 kg/m³, and removal is difficult when the density is 1500 kg/m³. A comparison of the simulation results indicates that since the density is less than that of the dust, the time required for tar mist removal is less than that of dust removal. The tar mist can be completely removed in the density range of 900–1100 kg/m³. When the dust density is 1500 kg/m³, the dust starts to gather at the bottom, and the degree of fluidization of the particles gradually decreases. At 120 s, the particle fluidization degree is the lowest, and the dust is difficult to remove.

6. Conclusions

In this paper, a physical model of the two-phase motion of furnace gas, tar, and dust in the internal structure of the internal heat type low-temperature dry distillation furnace is established. The Euler–Lagrangian method is used to simulate the two-phase flow. The OMEC Topsizer Plus laser particle size analyzer is used to measure the particle size distribution of the dust to verify the accuracy of the mathematical model and the reliability of the simulation results. On this basis, the effects of gas velocity, tar particle size, and dust particle size on the movement of tar mist and dust are investigated. The following main conclusions are drawn:(1)The tar particles with a particle diameter of less than 0.20 mm are removed from the furnace by the gas. The smaller the particle size, the shorter the time required for removal, the lower the probability of secondary pyrolysis of the tar, and the higher the tar recovery rate. In actual production, efforts should be made to reduce the particle size of the tar.(2)Dust with a particle size greater than 0.18 mm cannot be removed from the furnace and remains in the furnace for pyrolysis, which not only increases the yield of the tar but also reduces the content of particulate impurities in the tar. In actual production, the coal powder generated by the breakage of coal in the furnace is reduced as much as possible.(3)As the gas velocity increases, the time for the tar mist and dust to be completely removed gradually decreases. When the velocity is 0.70 m/s, the tar mist can be completely removed, but the dust with a particle diameter larger than 0.18 mm cannot be completely removed.(4)Providing theoretical guidance for industrial production, the results indicate that the tar produced by low-temperature dry distillation can be quickly removed by the gas, while the dust such as pulverized coal and semicoke remains in the furnace. When the gas velocity is controlled to 0.70 m/s, the tar particle diameter at a density of 900 kg/m³ is 0.10 mm, and the particle size of the dust with a density of 1500 kg/m³ is 0.20 mm.

Data Availability

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

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (no. 21663034) and the Natural Science Foundation Program of Shaanxi Province for Joint fund project (no. 2019JLM-45).

References

S. Peng, Bo Zhang, and T. Wang, “China’s coal resources: octothorpe shaped distribution characteristics and sustainable development strategies,” Engineering Sciences, vol. 17, no. 9, pp. 29–35, 2015.
View at: Google Scholar
Y. Zhu and Z. Wang, Coal Chemistry, Chemical Industry Press, Beijing, China, 2013.
Z. Liu, H. Chen, and W. Huang, “The current situation of coal resources & prospecting straegy in China,” Coal Technology, vol. 24, no. 10, pp. 1-2, 2005.
View at: Google Scholar
Y. Wang, D. Guo, H. Zhang et al., “Spatial distribution of coal resource potential in China and its application,” Journal of Natural Resources, vol. 21, no. 2, pp. 225–230, 2006.
View at: Google Scholar
K. Xie, Research on the Development Strategy of Energy, Chemical Industry Press, Golden Triangle, Thailand, 2016.
T. Yao, Li Wang, H. Jingtuan et al., “Study on coal liquefied technology,” Guangzhou Chemical Industry, vol. 42, no. 15, pp. 34-35, 2014.
View at: Google Scholar
J. Chen and X. Hong, “Coal gasification technology,” Chemical Engineering & Equipment, vol. 4, no. 4, pp. 110-111, 2011.
View at: Google Scholar
Y. Sun, “Analysis of the development status and trend of coal chemical technology,” China Petroch2em, vol. 23, no. 23, pp. 115-116, 2016.
View at: Google Scholar
T. Bai, “Coal low-temperature pyrolysis multi-generation technology-the best way to achieve clean and efficient use of coal,” Coal Processing & Comprehensive Utilization, vol. 1, no. 12, pp. 6–10, 2014.
View at: Google Scholar
L. Hao, W. Song, X. Zhang et al., “Application and development of pyrolysis technology in low-carbon economy,” Journal of Engineering Studies, vol. 4, no. 3, pp. 231–236, 2012.
View at: Google Scholar
Y. Wang, Y. Wang, Y. Dong et al., “Development prospect of semi-coke coal tar and dry distillation gas produced by coal dry distillation,” Petrochemical Industry Application, vol. 27, no. 6, pp. 1–3, 2008.
View at: Google Scholar
K. Cen, M. Ni, Z. Luo et al., “The prospect of power generation technology based on coal staged conversion,” Engineering Sciences, vol. 17, no. 9, pp. 118–122, 2015.
View at: Google Scholar
Y. He, “Economic analysis of coal pyrolysis and tar processing technology,” Energy and Transportation, vol. 18, no. 18, pp. 24-25, 2010.
View at: Google Scholar
Y. Song and J. Tang, Coal Chemical Technology, Chemical Industry Press, Beijing, China, 2016.
G.-Q. Chen and Z.-H. Luo, “New insights into intraparticle transfer, particle kinetics, and gas-solid two-phase flow in polydisperse fluid catalytic cracking riser reactors under reaction conditions using multi-scale modeling,” Chemical Engineering Science, vol. 109, pp. 38–52, 2014.
View at: Publisher Site | Google Scholar
J. Gao, C. Xu, S. Lin, G. Yang, and Y. Guo, “Simulations of gas-liquid-solid 3-phase flow and reaction in FCC riser reactors,” AIChE Journal, vol. 47, no. 3, pp. 677–692, 2001.
View at: Publisher Site | Google Scholar
C. Wu, Y. Cheng, and Y. Jin, “Understanding riser and downer based fluid catalytic cracking processes by a comprehensive two-dimensional reactor model,” Industrial & Engineering Chemistry Research, vol. 48, no. 1, pp. 12–26, 2009.
View at: Publisher Site | Google Scholar
Y. Wu, Q. Zhang, W. Yang, and W. Blasiak, “Two-dimensional computational fluid dynamics simulation of biomass gasification in a downdraft fixed-bed gasifier with highly preheated air and steam,” Energy & Fuels, vol. 27, no. 6, pp. 3274–3282, 2013.
View at: Publisher Site | Google Scholar
H. R. Nazif and H. Basirat Tabrizi, “Development of boundary transfer method in simulation of gas-solid turbulent flow of a riser,” Applied Mathematical Modelling, vol. 37, no. 4, pp. 2445–2459, 2013.
View at: Publisher Site | Google Scholar
S. Akhshik, M. Behzad, and M. Rajabi, “CFD-DEM approach to investigate the effect of drill pipe rotation on cuttings transport behavior,” Journal of Petroleum Science and Engineering, vol. 127, pp. 229–244, 2015.
View at: Publisher Site | Google Scholar
F. Chaumeil and M. Crapper, “Using the DEM-CFD method to predict Brownian particle deposition in a constricted tube,” Particuology, vol. 15, pp. 94–106, 2014.
View at: Publisher Site | Google Scholar
I. Goldasteh, G. Ahmadi, and A. R. Ferro, “Monte Carlo simulation of micron size spherical particle removal and resuspension from substrate under fluid flows,” Journal of Aerosol Science, vol. 66, pp. 62–71, 2013.
View at: Publisher Site | Google Scholar
J. Toraño, S. Torno, M. Menéndez, and M. Gent, “Auxiliary ventilation in mining roadways driven with roadheaders: validated CFD modelling of dust behaviour,” Tunnelling and Underground Space Technology, vol. 26, no. 1, pp. 201–210, 2011.
View at: Publisher Site | Google Scholar
J. C. Kurnia, A. P. Sasmito, A. S. Mujumdar, and S. Arun, “Dust dispersion and management in underground mining faces,” International Journal of Mining Science and Technology, vol. 24, no. 1, pp. 39–44, 2014.
View at: Publisher Site | Google Scholar
N. A. Patankar and D. D. Joseph, “Modeling and numerical simulation of particulate flows by the Eulerian-Lagrangian approach,” International Journal of Multiphase Flow, vol. 27, no. 10, pp. 1659–1684, 2001.
View at: Publisher Site | Google Scholar
Y. Lu, A. Saad, A. P. Sasmito, and J. C. Kurnia, “Prediction of air flow, methane, and coal dust dispersion in a room and pillar mining face,” International Journal of Mining Science and Technology, vol. 27, pp. 657–662, 2017.
View at: Publisher Site | Google Scholar
G. Zhou, Q. Zhang, R. Bai, T. Fan, and G. Wang, “The diffusion behavior law of respirable dust at fully mechanized caving face in coal mine: CFD numerical simulation and engineering application,” Process Safety and Environmental Protection, vol. 106, pp. 117–128, 2017.
View at: Publisher Site | Google Scholar
P. Cai, W. Nie, D. Chen, S. Yang, and Z. Liu, “Effect of air flowrate on pollutant dispersion pattern of coal dust particles at fully mechanized mining face based on numerical simulation,” Fuel, vol. 239, pp. 623–635, 2019.
View at: Publisher Site | Google Scholar
G. Sankar, D. S. Kumar, and K. R. Balasubramanian, “Computational modeling of pulverized coal fired boilers-a review on the current position,” Fuel, vol. 236, pp. 643–665, 2019.
View at: Publisher Site | Google Scholar
G. Zhou, B. Feng, W. Yin, and J. Wang, “Numerical simulations on airflow-dust diffusion rules with the use of coal cutter dust removal fans and related engineering applications in a fully-mechanized coal mining face,” Powder Technology, vol. 339, pp. 354–367, 2018.
View at: Publisher Site | Google Scholar
B. Taraba, Z. Michalec, V. Michalcová, T. Blejchař, M. Bojko, and M. Kozubková, “CFD simulations of the effect of wind on the spontaneous heating of coal stockpiles,” Fuel, vol. 118, pp. 107–112, 2014.
View at: Publisher Site | Google Scholar
R. Zhou, N. Yang, and J. Li, “CFD simulation of gas-liquid-solid flow in slurry bubble columns with EMMS drag model,” Powder Technology, vol. 314, pp. 466–479, 2017.
View at: Publisher Site | Google Scholar
J. C. Kurnia, A. P. Sasmito, and A. S. Mujumdar, “CFD simulation of methane dispersion and innovative methane management in underground mining faces,” Applied Mathematical Modelling, vol. 38, no. 14, pp. 3467–3484, 2014.
View at: Publisher Site | Google Scholar
J. Cheng, S. Li, F. Zhang, C. Zhao, S. Yang, and A. Ghosh, “CFD modelling of ventilation optimization for improving mine safety in longwall working faces,” Journal of Loss Prevention in the Process Industries, vol. 40, pp. 285–297, 2016.
View at: Publisher Site | Google Scholar
A. M. Eaton, L. D. Smoot, S. C. Hill, and C. N. Eatough, “Components, formulations, solutions, evaluation, and application of comprehensive combustion models,” Progress in Energy and Combustion Science, vol. 25, no. 4, pp. 387–436, 1999.
View at: Publisher Site | Google Scholar
W. Li, S. Lu, Y. Liu, R. Wang, Q. Huang, and J. Yan, “CFD simulation of the unsteady flow of a single coal log in a pipe,” The Canadian Journal of Chemical Engineering, vol. 93, no. 11, pp. 2084–2093, 2015.
View at: Publisher Site | Google Scholar
M. K. Neuman, J. W. Boulton, and S. Sanderson, “Wind tunnel simulation of environmental controls on fugitive dust emissions from mine tailings,” Atmospheric Environment, vol. 43, pp. 520–529, 2009.
View at: Publisher Site | Google Scholar
Z. Jiang, L. Jiang, and J. Chen, “Numerical simulation of dust concentration distribution regularities of down-the-hole drilling in open-pit mine,” Journal of Shenzhen University Science and Engineering, vol. 30, no. 3, pp. 313–318, 2013.
View at: Publisher Site | Google Scholar
M. W. Korevaar, J. T. Padding, M. A. V. D. Hoef, and J. A. M. Kuipers, “Integrated DEM-CFD modeling of the contact charging of pneumatically conveyed powders,” Powder Technology, vol. 258, pp. 144–156, 2014.
View at: Publisher Site | Google Scholar
G. Zhou, W. Cheng, L. Chen, and W. Nie, “Numerical simulation and its application of dust concentration spatial distribution regularities in fully-mechanized caving face,” Journal of China Coal Society, vol. 35, pp. 2094–2099, 2010.
View at: Google Scholar
A. Milićević, S. Belošević, N. Crnomarković et al., “Mathematical modelling and optimisation of lignite and wheat straw cocombustion in 350 MW_e boiler furnace,” Applied Energy, vol. 260, p. 114206, 2020.
View at: Google Scholar
Q. Zhang, G. Zhou, X. Qian, M. Yuan, Y. Sun, and D. Wang, “Diffuse pollution characteristics of respirable dust in fully-mechanized mining face under various velocities based on CFD investigation,” Journal of Cleaner Production, vol. 184, pp. 239–250, 2018.
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

Copyright

Copyright © 2020 Qiuli Zhang 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