Numerical Study on Combustion Behavior of Semi-Coke in Blast Furnace Blowpipe-Tuyere-Combustion Zone

: Injecting low-cost semi-coke is critical to reducing blast furnace production costs. The combustion behavior of the co-injection of semi-coke and undersized coke powder is still rarely studied. In this work, a three-dimensional CFD mathematical model was established to simulate the gas velocity, temperature, composition distribution and burnout of semi-coke during the combustion of pulverized semi-coke. The inﬂuence of mass fraction of semi-coke on the composition and burnout in the combustion zone of blast furnace was also studied. The results show that the maximum concentrations of CO and H 2 in the combustion zone are 36% and 8%, respectively. With the decrease of the semi-coke ratio in the blended coal, the ﬁxed carbon content and the caloriﬁc value of the blended coal increase, but the burnout of the blended coal reduces. When bituminous coal is single injected, the burnout reaches 70%, which is higher than that of semi-coke. In actual production, for the semi-coke and coke powder injecting into the blast furnace, a proportion of bituminous coal can be appropriately added to improve the burnout rate of the coal blends and increase the H 2 content in reducing gas.


Introduction
Blast furnace pulverized coal injection is an important technology to reduce production cost and stabilize operations [1,2].After the pulverized coal is injected into the blast furnace combustion zone from the tuyere, a violent and complex multiphase chemical reaction occurs in the limited space with high temperature and high pressure.At present, blast furnaces in China are mainly based on injected bituminous coal and anthracite.However, production is faced with the problem that high-quality coal resources are increasingly depleted and some local governments are restricting the total coal consumption of enterprises [3].On the other hand, semi-coke is mainly produced from low-temperature pyrolysis of weakly cohesive or non-cohesive coal with high volatile content [4][5][6].Due to the abundance of low-rank coal resources in China and the fact that local governments have not included semi-coke in the restricted sequence of corporate energy products [7], the efficient utilization of semi-coke in blast furnaces has become the focus of attention.
In the past, the comprehensive properties of semi-coke as a blast furnace injection fuel have been extensively studied.Zou et al. [8] investigated the inherent characteristics and injection effect of semi-coke, and determined the basic conditions which satisfy the blast furnace production.Li et al. [9] carried out the experimental analysis of semi-coke used as blast furnace injection fuel.Bi et al. [10] pointed out the optimal ratio of blast furnace injection of semi-coke to replace part of the bituminous coal and anthracite from the perspective of combustibility.In addition, some iron and steel enterprises have carried out industrial tests and production practices of blast furnace injection of semi-coke, proving its feasibility and cost superiority [11,12].Recently, considering the advantages of numerical simulation in terms of cost and efficiency, especially the fact that it can provide detailed furnace phenomena, it has become an important tool to study semi-coke injection in blast furnace [13][14][15][16][17]. Goto et al. [18] established a three-dimensional model to simulate the combustion process of pulverized coal-plastic and coal-gas in the blowpipe.Du [19,20] and Gu et al. [21] constructed a three-dimensional mathematical model of the blowpipe-tuyere area to describe the combustion behavior of pulverized coal under different operating conditions.However, they did not pay attention to the shape of the raceway and the development process of the pulverized coal after leaving the tuyere.Shen et al. [22][23][24][25] developed a series of mathematical models and conducted a detailed study on the combustion behavior of blast furnace pulverized coal injection, which provided a theoretical basis for guiding and optimizing the blast furnace pulverized coal injection process.Recently, Hu et al. [26] established a blast furnace raceway model and compared the differences between the individual injection of semi-coke and the mixed injection of semi-coke and pulverized coal.The previous studies provided effective information for guiding the blast furnace injection of semi-coke.However, the combustion behavior of the injection of semi-coke and undersized coke powder is still rarely studied, and especially the effect of different semi-coke ratios on the combustion behavior needs to be clarified.
In this work, a three-dimensional model of blast furnace blowpipe-tuyere-combustion zone based on computational fluid dynamics (CFD) is established to explore the velocity, temperature, component concentration distribution and burnout rate of semi-coke injection in blast furnace.The effect of the proportion of semi-coke on the composition of the combustion zone and the burnout rate is investigated, and the differences between injection semi-coke and bituminous coal are analyzed.The findings of this work will be useful for the control and optimization of blast furnace injection of semi-coke.

Governing Equations
In the simulation of the combustion process of the injection semi-coke, the fluid phase was treated as an incompressible ideal gas.The Euler method was adopted to solve the three-dimensional Reynolds-averaged Navier-Stokes equation.The governing equations include mass conservation, momentum conservation, energy conservation and component conservation equations, as shown in Table 1.Particularly, the gas velocity at the inlet of the blowpipe reaches more than 160 m/s, and the turbulence is fully developed.To accurately and effectively predict the turbulent mixing and diffusion behavior of dust and pulverized coal particles in high-speed airflow, realizable k-ε equations were employed to describe the turbulent diffusion [27].The equations of turbulent kinetic energy k and turbulent dissipation ε can be written as: where G k and G b are the turbulent kinetic energy terms generated by laminar velocity gradient and buoyancy, respectively.σ k and σ ε are the turbulent Prandtl number of the k equation and ε equation, respectively.Particularly, The pulverized coal was treated as a discrete phase and solved by the Langrangian method in the simulation.Considering that the pulverized coal was rapidly diluted by the gas phase after entering the blowpipe from the coal lance, the interaction between pulverized coal particles was therefore ignored [28].The trajectory of the pulverized coal particles was solved by the random trajectory model, and its motion followed Newton's second law, given by, where F D (u − u p ) is the drag force per unit particle mass.u p and u are the velocities of the particles and the gas phase, respectively.ρ p and ρ are the densities of the particles and the gas phase, respectively.F x is the sum of other forces on the particle.F D can be written as, Table 1.Governing equations of gas phase in the model.
The mass, momentum and energy equations of the particle phase are listed in Table 2.The temperature variation of the particle phase occurs in three ways, namely convective heat transfer, conduction heat transfer and radiation heat transfer.Among them, radiation heat transfer was solved by the P-1 radiation model, where the emissivity was estimated by the Gray Gas Weighted Average Model (WSGGM) [4].

Chemical Reaction Model
The combustion of pulverized coal in the tuyere and raceway is a complex reaction process.Generally, the combustion of pulverized coal undergoes four stages, namely preheating, releasing of volatile matter (VM), combustion of VM, oxidation and gasification of residual carbon in the turbulent gas phase.The chemical reaction models for these stages are as follows.

Releasing of Volatile Matter
A two-step competitive reaction model was adopted to describe the releasing process of VM.In other words, two chemical reactions with different reaction rates and volatilization rates compete with each other, and the expression is shown in Equation (5).
where VM 1 and VM 2 are volatile matters, Char 1 and Char 2 are residual carbon, and α 1 and α 2 are the stoichiometric coefficients of the volatile pyrolysis reaction of the pulverized coal, respectively.α 1 is the content of VM in the industrial analysis of pulverized coal, and the value of α 2 is related to α 1 .The chemical reaction rate constants k 1 and k 2 are expressed by the Arrhenius formula, given by, where A is the pre-exponential factor constant and E is the activation energy.

Combustion Reaction of Gas Phase
The VM generated by the pyrolysis of raw coal react with O 2 to produce CO 2 and H 2 O.This reaction was regarded as a one-step reaction and solved by the Finite Rate/Eddy Dissipation model.The effects of turbulent diffusion and chemical reaction rates were both considered in this model, and as a result, the reaction rate was expressed by the Arrhenius formula and the turbulent diffusion equation.The final net reaction rate depends on the smaller of the two, given by,

.3. Oxidation and Gasification Models of Residual Carbon
Another product after the pyrolysis of pulverized coal is considered to be pure carbon, that is, residual carbon, which will participate in the relevant oxidation and gasification reactions after being completely separated from the raw coal.Heterogeneous reactions involving carbon residues mainly include C (s) -O 2 , C (s) -CO 2 and C (s) -H 2 O, as shown in Equations ( 10)-( 12).
C (s) + 0.5O 2 → CO (10) In the simulations, a multi-interface chemical reaction model, which assumed that particle-phase surface species could be consumed or produced stoichiometrically by particle surface reactions, was employed to express the chemical reaction rates of these processes.The surface material of pulverized coal particle was regarded as pure carbon and participated in oxidation and gasification reactions.In addition, the kinetic rate k r of the chemical reaction was calculated by the Arrhenius formula, given by, The apparent order of chemical reaction r was 1; as a result, the carbon consumption rate R j,r on the surface of pulverized coal particles was calculated by Equation ( 14).
where A p and η r are particle surface area and effective factor, respectively.Y j denotes the mass fraction of material j on the particle surface.P n is the partial pressure of the gas phase.D 0,r is the diffusion rate coefficient of the chemical reaction r, given by where C 1,r is the molar concentration of substance j in chemical reaction r.The kinetic parameters associated with the above heterogeneous chemical reactions are summarized and listed in Table 3.
Table 3. Kinetic parameters of heterogeneous chemical reactions.

Geometric Model
The geometric model of this work was the same as the actual size of the blast furnace air supply pipe in a steel plant, including the blowpipe, the tuyere, the combustion zone and the coal injection lance, as shown in Figure 1a.The coal lance with a diameter of 35 mm was obliquely inserted into the blowpipe and the tuyere, and the outlet tip of the coal lance intersected the centerline of the blowpipe with an included angle of 9.8 • .The blowpipe was a tapered pipe structure with a length of 1050 mm, and its inlet and outlet diameters were 150 mm and 130 mm, respectively.The right side of the blowpipe was connected to a tuyere with a length of 570 mm and the included angle was 5 • .The combustion zone was set as a gradually expanding area with a length of 1740 mm and an expansion angle of 3 • .In the simulations, the pipe wall was treated with adiabatic no-slip boundary conditions, and the pressure at the inlet of blowpipe was set to 437 kPa.The outlet of the combustion zone was set as pressure outlet.The geometric model was divided in the form of mixed mesh with a mesh number of 120,000.Note that mesh refinement was performed in the coal injection lance section, as shown in Figure 1b.Particularly, the mesh-independence tests were conducted before the simulation.Four different mesh resolutions with cell numbers of about 60,000, 120,000, 180,000 and 240,000 were tested.The simulation results with a cell number of 60,000 are significantly different from the other conditions.To balance calculation accuracy and speed, the cell number of 120000 was thus selected in the simulations.

Simulation Conditions
The pulverized coal in the simulations is mainly composed of 90% semi-coke and 10% coke powder.Table 4 lists the data of industrial analysis and elemental analysis of semi-coke and coke powder.Particularly, the industrial analysis process was performed according to the national standards GB/T 2001-2013.The element analysis of the fuel was conducted in an YX-CHN5000 analyzer.The blended coal had a uniform particle size of 60 µm, a density of 1400 kg/m 3 , and a calorific value of 30,472.301kJ/kg.All simulation cases were performed on FLUENT 18.0.The SIMPLE algorithm was used to calculate the pressure-velocity coupled equations.The gradient and pressure in spatial discretization were calculated by least-square cell-based and PRESTO algorithms, respectively.The convergence standard required the residual curves of each index to be less than 1 × 10 −5 , and especially the residual curves of energy were less than 1 × 10 −6 .In the simulations, the operational parameters were set according to the actual production conditions of the blast furnace, which are listed in Table 5.
Table 5. Operational parameters adopted in the simulations.

Operational Parameters Value
Blast flow rate (Nm

Model Validation
Blast furnace coal injection is a complex black-box process with high temperature and high pressure, and the combustion effect of pulverized coal in the actual blast furnace is often difficult to detect.To validate the model, the measurement results from a pilot scale test rig in terms of coal burnout along the tuyere centerline were used to compare with simulation results.Note that the coal properties in the reference are close to that in this work (PC-4 in Section 3.3).The burnout indicates the combustion efficiency of coal, which is defined as: where m 0 is the ash content in the original coal, and m 1 is the ash content in the burnt residual collected.Figure 2 shows the comparison of the predicted burnout with the experimental results [29].Obviously, the simulated burnout (69.96%) is very close to the experimental result (91.7%), and the relative error is about 2.7%.This means that the model used in this work can accurately simulate the combustion behavior of pulverized coal.

Transmission Phenomenon in the Blowpipe-Tuyere-Combustion Zone 3.2.1. Velocity Distribution
Figure 3 shows the velocity distribution in the symmetrical section of the blowpipetuyere-combustion zone.As the diameter of the blowpipe decreases from 150 mm to 130 mm, the gas velocity shows a gradual increasing trend, and the gas velocity at the outlet of the blowpipe is about 250 m/s.The velocity of the gas used for conveying pulverized coal is 22 m/s at the outlet of the injection lance.In addition, the two airflows from the blowpipe and injection lance converge at the tuyere with the maximum gas velocity of 260 m/s.Due to the existence of the pulverized coal gas stream, the gas velocity in the central region is significantly lower than that in the edge region.Since the combustion zone is an expansion area with an included angle of 3 • , the gas velocity gradually decreases when it approaches the end of the combustion zone.A coke bed is located around the raceway in the blast furnace, and the high-speed gas will circulate after encountering the packed coke bed area.This work focuses primarily on combustion and thus simplifies the process.Figure 4 plots the velocity distribution at cross sections A and B of the combustion zone (marked in Figure 1).Obviously, the velocity in the central region of the combustion zone is lower than that in the edge region, which is due to the impediment of the jet flow by the low-velocity pulverized coal, and the complex chemical reaction around the pulverized coal also affects the flow of the gas.This phenomenon gradually weakens with the backward development of the combustion zone, and the velocity distribution of the flow field tends to be uniform.

Temperature Distribution
The temperature distribution contour on the symmetrical section is shown in Figure 5.It can be seen that the gas with 1190 • C is evenly distributed at the inlet of the blowpipe.The low-temperature pulverized coal is continuously heated during the process of entering the tuyere zone from the injection lance.The temperature of the gas and powder phases increase continuously in the injection lance, but no combustion occurred in this process due to the fast flow rate of the pulverized coal.There is an enriched low temperature zone at the tip of the injection lance, which is due to the release of volatiles in the pulverized coal as an endothermic process.With the release of volatiles, they undergo combustion reaction at a higher O 2 concentration near the pipe wall and release a lot of heat.Particularly, the maximum temperature reaches more than 2400 K.In addition, the temperature area at the front and upper of the combustion zone is significantly larger.This is because the injection lance is inserted into the tuyere at a certain inclination angle, and the pulverized coal is concentrated in the lower part of the combustion zone, so there is sufficient reaction space and oxygen concentration on the upper side to satisfy the oxidation reaction.As shown in Figure 6, the high temperature zone on the cross section of the combustion zone is annularly distributed.The central pulverized coal stream is not exposed to enough oxygen and can only be continuously heated during the process, thereby releasing volatiles.On the other hand, a vigorous combustion reaction occurs where the pulverized coal stream is in sufficient contact with the oxygen, thus forming an annular high temperature zone.With the extension of the depth of the raceway, more volatiles and residual carbon in the pulverized coal participate in the combustion reaction; the distribution of high temperature zones in the radial direction is therefore continuously expanding.This verifies that the high temperature zone formed in the raceway can provide sufficient heat for the hearth.

Composition Distribution
After the pulverized coal leaves the lance, it is quickly mixed with the high-temperature blast injected by the blowpipe, so the pulverized coal particles will be heated to the temperature condition for the release of VM in a short period of time, as shown in Figure 7.The pulverized coal begins to release VM after it leaves the injection lance at about 200 mm.Due to the continuous supply of pulverized coal, an enriched region with the highest volatiles content of 28% appeared at 200-600 mm.On other hand, the concentration distribution of VM and O 2 are just opposite, which is caused by the oxidation reaction of volatiles and O 2 .Generally, O 2 is continuously consumed as the reaction proceeds, and there is essentially no O 2 at the end of the raceway.Figure 8 shows the gas composition distribution in the blowpipe-tuyere-combustion zone.CO 2 and H 2 O are the main products of VM combustion, so their concentrations are higher near the annular high temperature region.With the removal of VM from the pulverized coal, the remaining residual carbon will continue to undergo oxidation or gasification reactions.CO and H 2 are the main products of C(s)-O 2 , C(s)-CO 2 , C(s)-H 2 O reactions.Thus, the enriched regions of CO and H 2 exist in regions with lower oxygen concentrations, and the highest concentrations of CO and H 2 are 36% and 8%, respectively.It should be noted that the concentration of H 2 is slightly reduced at the end of the raceway; this is because of the dilution of the H 2 concentration by the surrounding CO or CO 2 gas.
To further clarify the reaction mechanism of each gas component in the blowpipetuyere-combustion zone, the gas distribution data along the direction of the coal lance tip were extracted, as shown in Figure 9. Obviously, there are only two gases, O 2 and N 2 , at the discontinuity of the coal lance, and the concentration of O 2 decreases quickly.This is because the VM released by the pulverized coal consume O 2 in the combustion reaction, and the concentration of the combustion products CO 2 and H 2 O continues to increase.According to the above discussion, it can be determined that the residual carbon after the removal of VM from the pulverized coal will undergo a gasification reaction with CO 2 and H 2 O, but the CO concentration has been rising at the beginning.The main reason is that the concentration of O 2 is still higher at the position closer to the coal lance and more generated residual carbon will burn with O 2 to produce CO.Therefore, the concentration of CO increases continuously throughout the process, and finally becomes conserved.Particularly, in the mixed gas at the outlet of the raceway, the reducing gases CO and H 2 account for about 45%, while CO 2 and H 2 O are less than 10%.

Motion Behavior of Pulverized Coal Particles
Pulverized coal injection technology has been widely used in blast furnace ironmaking due to its advantages of economical, operational and environmental benefits.However, the high pulverized coal injection rate often leads to low burnout rate, and the unburned pulverized accumulates at the edge of the raceway, reducing the permeability of the coke bed and the stability of the blast furnace.Improving the burnout rate of pulverized coal is crucial to increase the amount of pulverized coal injection and reduce the production cost.A comprehensive understanding of the pulverized coal combustion process in the blast furnace tuyere raceway is critical for optimal selection of injection rates and blast parameters.Hence, the particle trajectories are employed to explore the motion behavior of the pulverized coal particles in the tuyere raceway.
Figure 10 shows the velocity, temperature and residual time distribution of pulverized coal particles.The velocity distribution of pulverized coal particles is almost the same as the flow field in the model, which indicates that the high-speed hot wind flow and conveying gas flow determine the motion state of pulverized coal in the tuyere and the combustion zone.The pulverized coal particles are quickly diluted after entering the tuyere, so the interaction between the particles can be ignored.In addition, the temperature variation of pulverized coal particles also confirms the rationality of the distribution of gas components.The continuously heated pulverized coal undergoes the process of VM removal, residual carbon combustion and gasification, providing a large amount of reducing gas and heat for the raceway.The pulverized coal is covered by the high-speed hot air flow during the movement.Since there is no circulation flow in the model, the longest residence time of the pulverized coal particles is only 70 ms.Therefore, the efficient combustion of pulverized coal in a very short period of time is one of the keys to maintaining the stable operation of the blast furnace.To clarify the combustion of the blended coal in the raceway, the variation of the burnout of the pulverized coal particles during the movement process was calculated, as shown in Figure 11.It can be found that the low-temperature pulverized coal moves faster near the tip of the coal lance, and the combustion is not sufficient.The VM in the pulverized coal are gradually released with the increase of temperature, and more residual carbon participates in the oxidation and gasification reactions, resulting in a sharp increase in the burnout of the pulverized coal.Since the O 2 is consumed in the end of the raceway, the composition and temperature distribution of the gas tends to be uniform; the burnout of pulverized coal changes little and finally stabilizes at about 67%.

Effect of Semi-Coke Ratio
In this section, the effect of the ratio of semi-coke to coke powder on the tuyere raceway is explored, and the scheme is listed in Table 6.Among them, PC-1 to PC-3 are blended coals composed of 90% semi-coke + 10% coke powder, 80% semi-coke + 20% coke powder and 70% semi-coke + 30% coke powder, respectively.PC-4 is a bituminous coal commonly used in blast furnaces for comparison.It can be seen that with the increase of the semi-coke ratio, the content of fixed carbon in the blended coal increases gradually, while the content of volatiles decreases.As a comparison, bituminous coal has a VM content of 31.33%, which is about 20% higher than that of blended coal.Figure 12 shows the variation trend of VM volume fraction and temperature along the tip of the coal lance with different ratios of semi-coke.With the decrease of semi-coke ratio in the blended coal, the VM content in the blended coal decreases, resulting in a decrease in the VM content precipitated at the tip of the coal lance.Obviously, due to the high VM content of bituminous coal (PC-4), the highest concentration of VM precipitated from bituminous coal reaches more than 20%.The compensation heat and O 2 required for complete combustion of bituminous coal per unit mass are more than that of blended coal.Thus, the maximum allowable injection amount of bituminous coal is smaller than that of semi-coke under the same conditions, and the injection of the same flow of bituminous coal will lead to a decrease in the overall temperature of the raceway and a decrease in heat income.In fact, the injection of bituminous coal is more favorable for operation with high air temperature and high oxygen enrichment.Due to the variation of the pulverized coal composition, the volume fraction of gas in the combustion zone also changes, as shown in Figure 13.When the blended coal is injected, the fixed carbon content and the overall temperature in the combustion zone are higher, which is conducive to the gasification reaction of the residual carbon.Hence, compared with bituminous coal (PC-4), blended coal (PC-1 to PC-3) has a lower CO concentration and a higher CO 2 concentration.In addition, with the increase in the proportion of coke powder in the blended coal, the same rule is also presented.It should be noted that the CO concentration of bituminous coal (PC-4) is significantly lower.This is because the volatiles in bituminous coal contain a lot of H 2 , and the proportion of CO in the reducing gas is therefore reduced.On the other hand, the increase of the proportion of H 2 enhances the reducibility of the generated gas, which is more conducive to the indirect reduction in the blast furnace.Figure 14 shows the variation of pulverized coal burnout rate with different ratios of semi-coke.Bituminous coal has the characteristics of high VM content, low ignition point and flammability, so its burnout rate (about 70%) is significantly higher than that of blended coal.For the blended coal, with the increase of the proportion of coke powder, the burnout rate decreases.This indicates that although the calorific value of the blended coal increases, the combustion performance is poor, and it still cannot provide sufficient heat for the operation of the blast furnace.Therefore, it is necessary to control the blending ratio of coke powder.In practice, in addition to semi-coke and coke powder, bituminous coal can be appropriately added.On the one hand, it can improve the combustion performance of pulverized coal, and on the other hand, it can provide H 2 with higher reduction efficiency, which is also in line with the current needs of hydrogen metallurgy.

Conclusions
In this work, a three-dimensional mathematical model of blast furnace injection of semi-coke was established to investigate the velocity, temperature, component concentration distribution and the combustion behavior of pulverized coal in the blowpipetuyere-combustion zone.The effect of the semi-coke ratio on the combustion behavior was discussed.The main findings can be summarized as: (1) During the blast furnace injection of semi-coke, the velocity in the central region of the combustion zone is lower than that in the edge region, but with the development of the combustion zone, the velocity distribution of the flow field tends to be uniform.The high temperature zone of the combustion zone is annularly distributed, and its radial distribution continues to expand with the extension of the depth of the combustion zone.(2) The highest concentrations of CO and H 2 in the combustion zone reach 36% and 8%, respectively.The longest residence time of semi-coke particles in the combustion zone is 70 ms under the current simulation conditions.(3) With the decrease of the semi-coke ratio in the blended coal, the fixed carbon content and the calorific value of the blended coal increase, but the burnout of the blended coal reduces.The burnout of bituminous coal is about 70%, which is significantly higher than that of blended coal.Therefore, an appropriate amount of bituminous coal can be added to improve the combustion performance of pulverized coal in actual production.

Figure 1 .
Figure 1.Geometric model (a) and mesh (b) of the blowpipe-tuyere-combustion zone in the simulation.

Figure 2 .
Figure 2. Model validation of the simulation results with the experimental results data from [29].

Figure 4 .
Figure 4. Velocity distribution on the cross sections A (a) and B (b) of the combustion zone.

Figure 6 .
Figure 6.Temperature distribution on the cross sections A (a) and B (b) of the combustion zone.

Figure 7 .
Figure 7. Mass fraction distributions of the VM (a) and O 2 (b).

Figure 9 .
Figure 9. Gas distribution in the direction of the coal lance tip.

Figure 11 .
Figure 11.Evolution of pulverized coal burnout rate along the coal lance tip.

Figure 12 .
Figure 12.Variation trend of VM volume fraction (a) and temperature (b) along the tip of the coal lance with different ratios of semi-coke.

Figure 13 .
Figure 13.Variation trends of CO 2 (a) and CO (b) concentrations along the tip of the coal lance with different ratios of semi-coke.

Figure 14 .
Figure 14.Variation of pulverized coal burnout with different ratios of semi-coke.

Table 2 .
Governing equations of particle phase in the model.

Table 4 .
Industrial and elemental analysis of fuels.

Table 6 .
Composition of blended coal under different simulation conditions.