Next Article in Journal
Comparison of Different Approaches to Predict the Performance of Pumps As Turbines (PATs)
Next Article in Special Issue
Value-Added Performance and Thermal Decomposition Characteristics of Dumped Food Waste Compost by Pyrolysis
Previous Article in Journal
Numerical Research of the Effect of Surface Biomimetic Features on the Efficiency of Tidal Turbine Blades
Previous Article in Special Issue
Pyrolysis of Grape Marc from Tunisian Wine Industry: Feedstock Characterization, Thermal Degradation and Kinetic Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Advanced One-Dimensional Entrained-Flow Gasifier Model Considering Melting Phenomenon of Ash

1
Department of Chemical Engineering, Pohang University of Science and Technology (POSTECH), Pohang 37673, Korea
2
Graduate Institute of Ferrous Technology, Pohang University of Science and Technology (POSTECH), Pohang 37673, Korea
*
Author to whom correspondence should be addressed.
Energies 2018, 11(4), 1015; https://doi.org/10.3390/en11041015
Submission received: 3 March 2018 / Revised: 17 April 2018 / Accepted: 17 April 2018 / Published: 21 April 2018
(This article belongs to the Special Issue Biomass Chars: Elaboration, Characterization and Applications Ⅱ)

Abstract

:
A one-dimensional model is developed to represent the ash-melting phenomenon, which was not considered in the previous one-dimensional (1-D) entrained-flow gasifier model. We include sensible heat of slag and the fusion heat of ash in the heat balance equation. To consider the melting of ash, we propose an algorithm that calculates the energy balance for three scenarios based on temperature. We also use the composition and the thermal properties of anorthite mineral to express ash. gPROMS for differential equations is used to solve this algorithm in a simulation; the results include coal conversion, gas composition, and temperature profile. Based on the Texaco pilot plant gasifier, we validate our model. Our results show good agreement with previous experimental data. We conclude that the sensible heat of slag and the fusion heat of ash must be included in the entrained flow gasifier model.

1. Introduction

Generally, an entrained flow gasifier (EFG) uses finely pulverized coal with steam and oxygen co-current to make syngas. This design forms a uniform internal temperature, and has a residence time of only a few seconds [1]. The coal conversion reaches approximately 100%, because the gasifiers use pulverized coal at high temperature. An EFG is not affected by the rank of the coal [2]. Currently many commercial EFGs are operated by enterprises such as (General Electric) GE, Shell, Siemens, CB&I, MHI and ThyssenKrupp [3]. These types of gasifiers are operated at a temperature higher than the AFT. The ash, which is a coal residue, is discharged in the form of molten slag. The slag that is discharged to the bottom has a considerable amount of sensible heat. The design of the gasifier should consider this heat, because it affects the internal maximum temperature.
Existing EFG models [4,5,6] focus on calculating the composition of the gas. To improve the model, representations of the reaction mechanism of coal have been improved. Previous models have proposed various reaction kinetic models, such as random pore model [7], shrinking core model [8,9], and shrinking sphere model [10]. In addition, equilibrium models have been suggested calculating the reaction between gases [11,12,13,14]. Most of the studies [7,8,9,10,11,12,13,15] have focused on developing reaction models and adjusting parameters. The energy balance is not considered sufficiently to find the optimal gasifier design. They have developed heat balance with two variables: (i) input and output heat flow; and, (ii) reaction heat.
The Texaco pilot plant is a typical EFG. Several studies (Table 1) [4,5,6,16,17] have attempted to model it based on experimental data [18] acquired from the Electric Power Research Institute (EPRI). Wen et al. (1979) proposed a model that uses three reaction zones: (i) pyrolysis and volatile combustion zone; (ii) combustion and gasification zone; and, (iii) only gasification zone; the model applies a Stokes’ law approximation instead of momentum balance [16]. Govind and Shah (1984) used the same kinetics as those of Wen et al., but neglected the momentum balance [17]. Vamvuka et al. (1995) used thermogravimetric analysis data to develop the kinetics based on bituminous coal, but does not consider the momentum balance [4,5]. Hwang et al. (2015) expressed two reaction zones without considering the ‘pyrolysis and volatile combustion’ zone; this model applies the Stokes’ law approximation, and adjust parameters, such as outer wall temperature and the reaction rate constant [6].
All of these models of the EFG have limitations. For energy balance, they all consider only input and output heat flow, and reaction enthalpy; they neglect energy that is absorbed by the melting of ash, and therefore do not accurately represent the inside of the gasifier in the real world. As a result, the calculated temperature is too high. For this purpose, previous papers introduce an additional term that represents heat loss to the outer wall. This calculation of heat loss requires assumptions about variables such as wall temperature, overall heat transfer coefficient, and thermal conductivity. These assumptions decrease the accuracy of the models.
The objective of this study is to improve the existing one-dimensional EFG model by including the ash-melting phenomenon instead of approximating it as heat loss at the outer wall. We propose a model to increase the accuracy of the temperature profile. The resulting model can predict the composition change of the product gas. We apply a shrinking sphere model to consider the combustion reaction, and then suggest reaction kinetics to calculate the amount of ash. We also design a new algorithm to consider the melting phenomenon of ash to improve the accuracy of the predicted temperature profile. We discuss the energy balance equation in three cases, according to the temperature: (i) temperature is lower than AFT; (ii) temperature in the first cell is higher than AFT; and, (iii) temperature in the second cell is higher than AFT. Finally, we compare the simulation results with the experimental results.

2. EFG Model

2.1. Basic Assumptions

To build the model, we assumed that:
  • The inside of the gasifier is cylindrical; this assumption is suitable for modeling the Texaco EFG [17].
  • Coal and gas mass flow rates are constant.
  • Temperature and gas concentration are uniform in the radial direction.
  • Each cell is perfectly mixed.
  • The reactor consists of equally divided cells (Figure 1). The reaction rate depends on each cell’s conditions.
  • Spherical coal particles react from the outer surface while moving through the cell from front to end. As the reaction of the coal progresses, the size of the particles decreases.
  • All of the coal grains in the same cell are uniform.
  • Ash changes to slag after temperature exceed an AFT.
  • Ashes are inert.

2.2. Reaction Kinetics

The reaction type can be divided into a heterogeneous reaction and a homogeneous reaction. The heterogeneous reaction indicates that the coal particles react with the gas. Coal reacts with oxygen, carbon dioxide, steam, and hydrogen (Table 2). The water gas shift reaction (WGS) and CO oxidation were considered as major reactions. The WGS reaction proceeds rapidly at a high temperature, and was therefore considered to be at equilibrium. CO oxidation was regarded as irreversible; as a consequence, gasification and combustion reaction could be expressed without dividing the reaction zone. Therefore, we solved reaction kinetics based on the single reaction zone; i.e., the six reaction schemes that are presented below were considered equally from the first cell to the last cells. The EFG does not produce much methane [19], so the methane-steam reforming reaction was not considered in this model.
To incorporate the coal components in reaction kinetics, the ultimate analysis was applied. The reaction kinetics (Table 3) of coal used in this study was based on Hwang et al.’s work. One difference is that we have attempted to quantify ash. We introduce subscript f, which means that the ash component contained in 1 g of coal, and is the result of the proximate analysis. The subscripts a through e represent elements contained in 1g of coal. For heterogeneous reactions, the reaction rate of coal was proportional to both the surface area of the coal particles and the partial pressure of each gas. In homogeneous reactions, all of the gases were assumed to be ideal.

3. Solving Procedure

3.1. Mass Balance

For mass balance, only production and consumption due to the reaction were calculated. The kinetics of the coal reaction was based on one coal grain. The number of coal particles per unit volume came from the following Equation (1):
N v = W c , 0 m c v c A
The four coal reactions caused the coal conversion when the coal particles passed through the unit cell. ∆z is unit cell length. The change of coal mass was calculated as
Δ W c = N v · A i = 1 4 R i · Δ z
Mass balance of the ash could be expressed in the same way as that of coal. The only difference is that the stoichiometric coefficients f are added. The stoichiometric coefficients that correspond to ashes were all equal to f (Equation (3)). For slag, the mass balance was not considered separately. When the temperature was higher than AFT, ash was regarded as slag.
Δ W a = N v · A i = 1 4 R i · f · Δ z
For gas mass balance, R5 and R6 were also included. R6 is irreversible, so it was applied in the same manner as the heterogeneous reactions. α means the degree of deviation from equilibrium. To calculate α, we used Equation (4)
K   =   ( F C O 2 α ) · ( F H 2 α ) ( F C O + α ) · ( F H 2 O + α )
After converting this to an explicit form, we took a small root of the quadratic equation. The reason was that large roots make the mole flow negative.
Equation (5) is the mass balance of gases where ξ i is a stoichiometric coefficient for each component. Heterogeneous and homogeneous reactions were considered in different forms. The terms on the right side of the equation mean in order: (i) initial value; (ii) heterogeneous reactions; (iii) CO oxidation; and, (iv) each cell reaches equilibrium.
F i = F i , 0 ( N v · A · Δ z · j = 1 4 ξ i , j · R j M W i ) + ( R 5 · ξ i , 5 ) ± α

3.2. Energy Balance

Heat flow of input and output and reaction enthalpy were considered equally in all cells. In some existing models [5,6,16,17], oxygen, steam, and coal temperature were set differently. When the temperatures of coal and gas are separately calculated, the temperature difference between them is generally ~10 K [22]. In this study, we simplified the problem by assuming that the gas temperature was the same as the coal temperature.
Coal residue is ash or slag depending on temperature condition (Figure 1). To consider ash melting, energy balance was calculated by dividing it into three cases. At temperature <AFT, Equation (6) was applied. It did not incorporate the fusion heat of ash.
Δ { ( W c · c p , c + W a · c p , a + W s · c p , s + i F i · c p , i ) T } =   N v · A · j = 1 6 R j · Δ H j H l o s s
Cells with temperature >AFT were considered using Equations (7) and (8). Equation (7) was applied only to the first cell, which has temperature >AFT. This equation calculates all of the latent heat of the accumulated ashes.
Δ { ( W c · c p , c + W a · c p , a · W s · c p , s + i F i · c p , i ) T } =   N v · A · j = 1 6 R j · Δ H j   +   W a · Δ H a s h H l o s s
Equation (8) was applied to the remaining cells. It considers the heat of fusion of ash that occurs as it passes through each cell.
Δ { ( W c · c p , c + W a · c p , a · W s · c p , s + i F i · c p , i ) T } =   N v · A · j = 1 6 R j · Δ H j   +   W a · Δ H a s h H l o s s
Calculation of the heat loss must consider radiation, convection, and conduction. In addition, the thermal conductivity varies depending on the material of the refractory. If these cannot be reasonably estimated, the overall heat transfer coefficient will produce large errors. In a previous study, the heat transfer coefficient was indirectly estimated instead of calculating the heat loss [6,16,17]; the authors claimed that 30% of the reaction heat was lost to the outer wall, and that the corresponding calorific value is about 7–10% of the heating value of coal. However, heat losses typically considered in an EFG range from 1 to 4% [23,24]. This value is dependent on the gasifier scale. When compared with industrial gasifiers, pilot-scale gasifiers have higher heat losses. Therefore, we assumed a heat loss of 4%.

3.3. Solving Algorithm

The length of the gasifier was divided into 1650 cells (500 parts/m). Each cell was a system. Mass and energy balance were calculated sequentially. The computational algorithm was terminated when it reached the end of the reactor, and we obtain a sufficiently smooth curve in the simulation result. The information calculated in each cell was as follows:
  • Information about coal: coal conversion, coal mass flow rate.
  • Molar flow rate of product gas.
  • Temperature profile.
The proposed algorithm (Figure 2) was designed to apply the energy balance of three cases according to temperature. gPROMS simulation software (Process Systems Enterprise Ltd., London, UK) was used to perform the algorithm. The DASolver was applied; it could be used to solve differential equations in the steady state. Both the relative tolerance and absolute tolerance were fixed as 1.0 × 10−5.

4. Required Information for Simulation

4.1. Operating Variables and Reactor Size

Required inputs, including operating variables and particle sizes (Table 4), were those of the existing pilot gasifier. These results were compared with those of previous studies. The experimental conditions were presented in the work of Govind and Shah, which is the same as for their simulation conditions. Experimental information on the outer wall temperature of the gasifier was not provided. The existing model assumed that the temperature of the outer wall started at 2100 K, and then decreased linearly [16,17]. We indirectly estimated the amount of heat loss from the coal heating value, instead of assigning an initial outer wall temperature.

4.2. Coal Properties

This study considered Illinois no. 6 coal, which is the same type using previous modeling studies. The results (Table 5) of coal elemental analysis were acquired from the EPRI [18]. Density and specific heat capacity of the coal were the same as in previous work [16,17]. We considered a smaller particle size (41 μm) of coal than did previous work, for two reasons. First, 41 μm is closer to the size of the particles that were used in the EFG in real industry; particularly, when entering a dry feed, the particle size is < 100 μm [25]. Second, the reaction of coal is directly related to particle size and kinetic parameters, so particle size is an important factor in kinetics. Therefore, the kinetic parameter and particle size must be taken from the same reference. We used the same kinetics parameters and the coal particle size as Hwang’s model [6].

4.3. Ash and Slag Properties

The EPRI provided only the coal analysis data; it does not provide information about ash. One of way to remedy this lack of information is to assume the ash component. However, the physical properties differ depending on the ash component. Calculate of thermal properties that consider its components is a difficult task. In some documents [26,27], a formula for predicting thermal properties according to the temperature have been proposed. However, they were out of the range of gasifier operating temperature.
In this study, different approaches were used. Coal ashes are composed of minerals (Figure 3). About 70–90% of the constituents are CaO, Al2O3, and SiO2. They also include such as Fe2O3 and MgO [28,29]. In this work, we adopt the thermal properties of gehlenite (Figure 3) to represent ash.
The properties (Table 6) of ash and slag used in this study were obtained from the literature because the ash was considered to be a well-known mineral, and its thermal properties were set accordingly.

5. Results and Discussion

5.1. Model Valdation

To validate the model, several main variables were selected: final coal conversion rate; major product gas composition at exit; and, the hydrogen-to-carbon monoxide ratio (H2/CO). The results (Table 7) of this simulation are presented together with the results of the previous researchers and experimental data. We did not arbitrarily adjust the kinetic parameters that were used in this study. We applied the same reactor size and operating conditions of the existing pilot gasifier. For this reason, results of this work were similar to results of previous studies. In addition, we considered the melting phenomenon of ash that was not considered in the existing one-dimensional (1-D) model and reduced the estimated heat loss in the previous model. As a result, the modeling results are improved.

5.2. Simulation Results

5.2.1. Coal Conversion

Most of the reaction proceeded rapidly at the front of the gasifier (Figure 4). Only 10% of coal was reacted at around 0.15 m of the gasifier, but 80% had reacted at 0.24 m. This increase occurred because the combustion reactions were accelerated. The reaction rate of coal decreased as the size of coal particles decreased. Especially after oxygen was completely consumed, the reaction rates were very slow. This trend is characteristic of typical EFGs; it is consistent with the results in past research [5,6,16]. The final coal conversion was 98.8%.

5.2.2. Gas Composition

The mole fractions of the major gases changed over the reactor length (Figure 5). Oxygen was abruptly consumed near 0.25 m; this change is the result of rapid combustion. Syngas was generated in an oxygen-free environment. Steam was slightly generated in the front of the gasifier; afterwards, the proportion of steam was controlled by WGS equilibrium. After all of the oxygen was consumed, the gas composition did not change significantly. CH4, H2S, and N2 were < 1% of the product; they are not represented in Figure 5. Trends in the graph agreed with trends that were reported in previous research [5,6,16]. The exit gas consists mainly (mole fraction 98.9%) of carbon dioxide, hydrogen, steam, and carbon monoxide. All of these compositions are determined by WGS equilibrium, which is a function of temperature and is closely related to the outlet temperature of the gasifier. The attempt to calculate the temperature is the first step in predicting the composition of the gas, and this study is significant in that respect.

5.2.3. Temperature and Heat Flow

Simulations were used to calculate two temperature profiles (Figure 6). The heat balance of the first case (Figure 6, blue line) considered the sensible heat of the slag and the latent heat of the ash. The second case (Figure 6, red line) neglected these phenomena; both were set to 0. Below AFT, the temperatures of the two cases were the same, but the peak and outlet temperatures differed between the two cases. The same results were obtained in the temperature range below AFT; they are acceptable because the melting of ash had not yet been considered. The maximum temperature was calculated as 2112 K when sensible heat of the slag and the latent heat of the ash were considered, but 2155 K when they were neglected. The exit temperatures were 1464 K when the sensible heat of the slag and the latent heat of the ash were considered, and 1521 K when they were neglected; the difference between the two outlet temperatures was 57 K. The difference can be explained, as follows. The heat capacity of ash and slag is not taken into account in the system, and the corresponding energy was transferred to the gas, so the calculated temperature increased. Previous models have released heat to outside. EFGs generally use thick refractories. The assumption that a significant amount of heat escapes to the outside must be modified: these two temperature profiles show that the ash melting phenomenon must be considered when the internal temperature of the gasifier is calculated. This is a reasonable conclusion, given the fact that the internal peak temperature of the gasifier is higher than the AFT.
In this study, ash accounted for 15% of the mass of the coal and 7% of the total mass. The effect of sensible heat of slag is evident when the heat flows are divided into the solid and gas (Figure 7). Solids include coal, ash, and slag. After most of the reaction has proceeded, the solids are mainly in slag form. Especially at the exit condition, >99% of the solid was slag, which accounted for 5.0% of the total heat. In the coal that was used in this study, the ash was not negligible. As a result, our new attempt was meaningful. This phenomenon must be considered, especially for coal that has high ash content.

5.3. Consideration of Ash Melting Effect

The melting of ash has a similar effect on energy balance, as does heat loss from the outer wall. Because the ash melts and absorbs energy, the temperature of the system is lowered. In contrast, the loss to the outer wall lowers the temperature because heat escapes from of the system. These are different phenomena that occurred inside the gasifier, and were considered independently.
We quantitatively calculated the melting effect. This approach used the same method that was used to consider the heat loss. Based on the coal heating value, the latent heat of the ash and the sensible heat of the slag were calculated. The melting effect at the outlet was ~1%. This is not negligible given that heat loss is assumed to be 4% in this study. We calculated the melting effect at the exit as
S l a g   S e n s i b l e   H e a t + A s h   F u s i o n   H e a t C o a l   H e a t i n g   V a l u e ( H H V )   =   W s , L · c P , s , L · T L + W a , L · X L · Δ H a s h W c , 0   · H H V   o f   I l l i n o i s   N o . 6
In this way, we corrected the overestimate of heat loss that occurred in previous studies. In addition, we could explain some of the uncertainty of heat loss.

5.4. Applicability of High Ash Content Coal and Limitations

Coal containing a large amount of ash can be applied in the same manner. If information about the ash component is available, then it can be used to calculate the slag properties and incorporate thermal properties in the heat balance. When the ash content is high, modeling errors can be reduced by considering the fusion heat of ash and the sensible heat of the slag.
However, the reaction kinetic parameters that are cited in this study represent only coal combustion and gasification, so the model is only applicable to the entrained coal gasifier. In addition, lab-scale or pilot-scale data of other coal that has high ash content should be used to further assess applicability of this model.

6. Conclusions

This study suggests that a heat-balance model of an EFG must consider the effects of melting of ash. We attempted modeling based on kinetic models similar to those of previous researchers. Our proposed 1-D model is the only one that includes ash in the heat-balance and temperature-calculation algorithm. Gas production and coal conversion trends were similar to those of existing ones, and the results at the exit were mostly consistent with the experiments.
This result is meaningful in that it reflects actual phenomena occurring inside the EFG. Ash melts in any slagging-type gasifier. We can expect to calculate the internal temperature more accurately based on this study. These results can be used to guide the choice of design elements of EFG, such as material and thickness of the refractory wall.
One limitation of this model is that information on the ash component was not available. We used minerals to express ash components. This demerit must be eliminated. For further research, advanced modelling should use thermal properties based on ash analysis data. This model can be extended to account for radial direction temperature distribution.

Acknowledgments

This work was financially supported by the Pohang Iron & Steel Company Limited in South Korea (No. 2017Y005).

Author Contributions

Jinsu Kim programmed the gasifier model and wrote draft paper. Hyunmin Oh checked the coding results, analyzed the simulation results, and corrected minor errors. Seokyoung Lee did the validation work and wrote the figure. Young-Seek Yoon proposed idea, supervised research based on literature review and finished the paper.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

Across section area of gasifier (m2)
Cmolar concentration (mole/m3)
cpspecific heat capacity (J/g·K for solid, J/mole·K for gas)
dccoal particle size (m)
Fmolar flow rate of gas (mole/s)
Henthalpy (J)
Hlossheat loss to outer wall (J/s)
HHVhigh heating value of coal (J/g)
Kequilibrium constant of water gas shift reaction (-)
kirate constant of ith reaction
Lreactor length (m)
mcmass of one coal particle (g/#)
MWmolecular weight (g/gmole)
Nvnumber of particles contained in coal (#/m3)
Pipartial pressure of gas i (atm)
Rjreaction rate of j (g/s for heterogeneous reaction, mole/s for homogeneous reaction)
Ttemperature (K)
vccoal velocity (m/s)
Waash mass flow rate (g/s)
Wccoal mass flow rate (g/s)
Wsslag mass flow rate (g/s)
wphi control function according to temperature (-)
XCoal conversion (-)
Zthe coordinates of the axis of the gasifier in the previous model (m)
zthe coordinates of the axis of the gasifier in this model (m)
Greek Characters
αdegree of deviation from equilibrium (mole/s)
φthe stoichiometric coefficient to adjust the complete, incomplete combustion
ξstoichiometric coefficient
Subscripts
athe weight of C contained in 1 g of coal
bthe weight of H contained in 1 g of coal
cthe weight of O contained in 1 g of coal
dthe weight of N contained in 1 g of coal
ethe weight of S contained in 1 g of coal
fthe weight of ash contained in 1 g of coal
igas species
jreaction
Lvalue at the reactor exit
0value at the reactor inlet

References

  1. Phillips, J. Different types of gasifiers and their integration with gas turbines. In The Gas Turbine Handbook; 2006; Volume 1. Available online: https://www.netl.doe.gov/File%20Library/Research/Coal/energy%20systems/turbines/handbook/1-2-1.pdf (accessed on 2 January 2018).
  2. Guo, X.; Dai, Z.; Gong, X.; Chen, X.; Liu, H.; Wang, F.; Yu, Z. Performance of an entrained-flow gasification technology of pulverized coal in pilot-scale plant. Fuel Process. Technol. 2007, 88, 451–459. [Google Scholar] [CrossRef]
  3. Wang, T.; Stiegel, G.J. Integrated Gasification Combined Cycle (IGCC) Technologies; Woodhead Publishing: Cambridge, UK, 2016; ISBN 978-0-08-100185-1. [Google Scholar]
  4. Vamvuka, D.; Woodburn, E.T.; Senior, P.R. Modelling of an entrained flow coal gasifier. 2. Effect of operating conditions on reactor performance. Fuel 1995, 74, 1461–1465. [Google Scholar] [CrossRef]
  5. Vamvuka, D.; Woodburn, E.T.; Senior, P.R. Modelling of an entrained flow coal gasifier. 1. Development of the model and general predictions. Fuel 1995, 74, 1452–1460. [Google Scholar] [CrossRef]
  6. Hwang, M.; Song, E.; Song, J. One-Dimensional Modeling of an Entrained Coal Gasification Process Using Kinetic Parameters. Energies 2016, 9, 99. [Google Scholar] [CrossRef]
  7. Bhatia, S.K.; Perlmutter, D.D. A random pore model for fluid-solid reactions: II. Diffusion and transport effects. AIChE J. 1981, 27, 247–254. [Google Scholar] [CrossRef]
  8. Ishida, M.; Wen, C.Y. Comparison of kinetic and diffusional models for solid-gas reactions. AIChE J. 1968, 14, 311–317. [Google Scholar] [CrossRef]
  9. Sadhukhan, A.K.; Gupta, P.; Saha, R.K. Modelling of combustion characteristics of high ash coal char particles at high pressure: Shrinking reactive core model. Fuel 2010, 89, 162–169. [Google Scholar] [CrossRef]
  10. Gupta, P.; Sadhukhan, A.K.; Saha, R.K. Analysis of the combustion reaction of carbon and lignite char with ignition and extinction phenomena: Shrinking sphere model. Int. J. Chem. Kinet. 2007, 39, 307–319. [Google Scholar] [CrossRef]
  11. Kong, X.; Zhong, W.; Du, W.; Qian, F. Three Stage Equilibrium Model for Coal Gasification in Entrained Flow Gasifiers Based on Aspen Plus. Chin. J. Chem. Eng. 2013, 21, 79–84. [Google Scholar] [CrossRef]
  12. Biagini, E.; Bardi, A.; Pannocchia, G.; Tognotti, L. Development of an Entrained Flow Gasifier Model for Process Optimization Study. Ind. Eng. Chem. Res. 2009, 48, 9028–9033. [Google Scholar] [CrossRef]
  13. Chen, C.; Horio, M.; Kojima, T. Numerical simulation of entrained flow coal gasifiers. Part I: Modeling of coal gasification in an entrained flow gasifier. Chem. Eng. Sci. 2000, 55, 3861–3874. [Google Scholar] [CrossRef]
  14. Eri, Q.; Wu, W.; Zhao, X. Numerical Investigation of the Air-Steam Biomass Gasification Process Based on Thermodynamic Equilibrium Model. Energies 2017, 10, 2163. [Google Scholar] [CrossRef]
  15. Liu, X.; Wei, J.; Huo, W.; Yu, G. Gasification under CO2–Steam Mixture: Kinetic Model Study Based on Shared Active Sites. Energies 2017, 10, 1890. [Google Scholar] [CrossRef]
  16. Wen, C.Y.; Chaung, T.Z. Entrainment Coal Gasification Modeling. Ind. Eng. Chem. Process Des. Dev. 1979, 18, 684–695. [Google Scholar] [CrossRef]
  17. Govind, R.; Shah, J. Modeling and simulation of an entrained flow coal gasifier. AIChE J. 1984, 30, 79–92. [Google Scholar] [CrossRef]
  18. Robin, A.M. Hydrogen Production from Coal Liquefaction Residues; Final Report; Montebello Research Lab, Texaco, Inc.: Montebello, CA, USA, 1976. [Google Scholar]
  19. Basu, P. Biomass Gasification and Pyrolysis: Practical Design and Theory; Academic Press: New York, NY, USA, 2010; ISBN 978-0-08-096162-0. [Google Scholar]
  20. Callaghan, C.A. Kinetics and Catalysis of the Water-Gas-Shift Reaction: A Microkinetic and Graph Theoretic Approach. Ph.D Dissertation, Worcester Polytechnic Institute, Worcester, MA, USA, 2006. [Google Scholar]
  21. Deng, Z.; Xiao, R.; Jin, B.; Huang, H.; Shen, L.; Song, Q.; Li, Q. Computational Fluid Dynamics Modeling of Coal Gasification in a Pressurized Spout-Fluid Bed. Energy Fuels 2008, 22, 1560–1569. [Google Scholar] [CrossRef]
  22. Breault, R.W. Gasification Processes Old and New: A Basic Review of the Major Technologies. Energies 2010, 3, 216–240. [Google Scholar] [CrossRef]
  23. Gazzani, M.; Manzolini, G.; Macchi, E.; Ghoniem, A.F. Reduced order modeling of the Shell–Prenflo entrained flow gasifier. Fuel 2013, 104, 822–837. [Google Scholar] [CrossRef]
  24. Tremel, A.; Becherer, D.; Fendt, S.; Gaderer, M.; Spliethoff, H. Performance of entrained flow and fluidised bed biomass gasifiers on different scales. Energy Convers. Manag. 2013, 69, 95–106. [Google Scholar] [CrossRef]
  25. Krishnamoorthy, V.; Pisupati, S.V. A Critical Review of Mineral Matter Related Issues during Gasification of Coal in Fixed, Fluidized, and Entrained Flow Gasifiers. Energies 2015, 8, 10430–10463. [Google Scholar] [CrossRef]
  26. Richardson, M.J. The specific heats of coals, cokes and their ashes. Fuel 1993, 72, 1047–1053. [Google Scholar] [CrossRef]
  27. Rezaei, H.R.; Gupta, R.P.; Bryant, G.W.; Hart, J.T.; Liu, G.S.; Bailey, C.W.; Wall, T.F.; Miyamae, S.; Makino, K.; Endo, Y. Thermal conductivity of coal ash and slags and models used. Fuel 2000, 79, 1697–1710. [Google Scholar] [CrossRef]
  28. Mills, K.C.; Rhine, J.M. The measurement and estimation of the physical properties of slags formed during coal gasification: 1. Properties relevant to fluid flow. Fuel 1989, 68, 193–200. [Google Scholar] [CrossRef]
  29. Blissett, R.S.; Rowson, N.A. A review of the multi-component utilisation of coal fly ash. Fuel 2012, 97, 1–23. [Google Scholar] [CrossRef]
  30. Žigo, O.; Adamkovičová, K.; Kosa, L.; Nerád, I.; Proks, I. Determination of the heat of fusion of 2CaOAl2O3 SiO2 (gehlenite). Chem. Pap. 1987, 41, 171–181. [Google Scholar]
  31. Richet, P.; Fiquet, G. High-Temperature Heat Capacity and Premelting of Minerals in the System MgO-CaO-Al2O3-SiO2. J. Geophys. Res. 1991, 96, 445–456. [Google Scholar] [CrossRef]
  32. Berman, R.G.; Brown, T.H. Heat capacity of minerals in the system Na2O-K2O-CaO-MgO-FeO-Fe2O3-Al2O3SiO2-TiO2-H2O-CO2: Representation, estimation, and high temperature extrapolation. Contrib. Mineral. Petrol. 1985, 89, 168–183. [Google Scholar] [CrossRef]
Figure 1. Gasifier internal scheme that considers melting of ash. Red dot: ash; brown dot: slag.
Figure 1. Gasifier internal scheme that considers melting of ash. Red dot: ash; brown dot: slag.
Energies 11 01015 g001
Figure 2. Algorithm that includes the melting of ash.
Figure 2. Algorithm that includes the melting of ash.
Energies 11 01015 g002
Figure 3. Ternary phase diagram of Al2O3-CaO-SiO2. Brown: gehlenite.
Figure 3. Ternary phase diagram of Al2O3-CaO-SiO2. Brown: gehlenite.
Energies 11 01015 g003
Figure 4. The profile of the coal conversion and coal mass flow rate.
Figure 4. The profile of the coal conversion and coal mass flow rate.
Energies 11 01015 g004
Figure 5. The profile of the major product gas mole fraction.
Figure 5. The profile of the major product gas mole fraction.
Energies 11 01015 g005
Figure 6. Differences in temperature profile depending on whether melting of ash is considered or not.
Figure 6. Differences in temperature profile depending on whether melting of ash is considered or not.
Energies 11 01015 g006
Figure 7. Heat flow of solid and gas.
Figure 7. Heat flow of solid and gas.
Energies 11 01015 g007
Table 1. Texaco pilot plant entrained flow gasifier (EFG) model based on the Electric Power Research Institute (EPRI) data.
Table 1. Texaco pilot plant entrained flow gasifier (EFG) model based on the Electric Power Research Institute (EPRI) data.
ResearcherKineticsMomentumEnergy Balance
Wen et al.3 reaction zonesStokes’ law approximationThermal
Govind and Shah3 reaction zonesNot consideredThermal
Vamvuka et al.Parameter based on thermogravimetric analysisNot consideredThermal
Hwang et al.2 reaction zonesStokes’ law approximationThermal
Table 2. Homogeneous and heterogeneous reaction list used in this study.
Table 2. Homogeneous and heterogeneous reaction list used in this study.
Reaction TypeReactionDescription
Heterogeneous R 1 :   C a H b O c N d S e A f + ( a ϕ * + b 4 c 2 e 2 ) O 2 ( 2 1 ) a C O 2 + 2 ( 1 1 ) a C O + ( b 2 e ) H 2 O + e H 2 S + d 2 N 2 + f A Coal↔O2
R 2 :   C a H b O c N d S e A f + a C O 2 2 a C O + c H 2 O + ( b 2 c e ) H 2 + e H 2 S + d 2 N 2 + f A Coal↔CO2
R 3 :   C a H b O c N d S e A f + ( a c ) H 2 O a C O + ( a + b 2 c e ) H 2 + e H 2 S + d 2 N 2 + f A Coal↔H2O
R 4 :   C a H b O c N d S e A f + ( 2 a b 2 + c e ) H 2 a C H 4 + c H 2 O + e H 2 S + d 2 N 2 + f A Coal↔H2
Homogeneous R 5 :   C O +   H 2 O   C O 2 + H 2 WGS
R 6 :   C O +   1 2 O 2 C O 2 CO oxidation
* φ is a function of absolute temperature T, and adjusts the ratio of complete and incomplete combustion. =   2 w + 2 w + 2 ,   w = 2500   · exp ( 6249 T ) . [4].
Table 3. Reaction kinetics and parameters that were used in this study.
Table 3. Reaction kinetics and parameters that were used in this study.
action TypeReaction Rate (g/s)k (g/(m2·atm·s))Reference
Heterogeneous k 1 π d c 2 P O 2 6180 · e ( 10 , 233.9 T ) [5,6]
k 2 π d c 2 P C O 2 198 , 100 · e ( 10 , 233.9 T )
k 3 π d c 2 P H 2 O 198 , 100 · e ( 10 , 233.9 T )
k 4 π d c 2 P H 2 385 · e ( 17 , 451.7 T )
HomogeneousEquilibriumEquilibrium ConstantReference
K = P C O 2 P H 2 P C O P H 2 O 2.4198 + 0.003855 · T + 2180.6 T [20]
Reaction Rate (mole/s)k (m3/(mole · s))Reference
k 6 C C O C C O 2 3.09 · e ( 11 , 199 T ) [21]
Table 4. Comparison of operating variables and gasifier size between this and previous studies.
Table 4. Comparison of operating variables and gasifier size between this and previous studies.
Operating VariablesThis WorkWenGovind
Coal feed rate (g/s)507577
Steam to coal ratio (-)0.240.240.241
Oxygen to coal ratio (-)0.860.860.86
Feed coal temperature (K)900900505
Feed gas temperature (K)900900697(H2O), 298(O2)
Gasifier pressure (MPa)2.02.02.4
Gasifier wall temperature (K)Not considered2100–600(Z/L)2100–600(Z/L)
Gasifier SizeThis WorkWenGovind
Gasifier length (cm)330330330
Inner diameter (cm)152152152
Table 5. Element analysis (wt %), coal properties, and feed conditions that were used in this study.
Table 5. Element analysis (wt %), coal properties, and feed conditions that were used in this study.
ElementMeasurement
C74.05
H6.25
O1.32
N0.71
S1.77
Ash15.33
Coal PropertiesAssumed
Density (g/cm3)1.80
Specific heat capacity (cal/(g·K))0.45
Feed ConditionsSet
Feed particle size ( μ m . )41
Velocity (cm/s)50
Table 6. Thermal properties of ash and slag used in this study.
Table 6. Thermal properties of ash and slag used in this study.
Thermal PropertiesValueReference
Fusion heat of ash (J/g)627.6[30]
Ash fusion temperature (K)1863[30]
Specific heat capacity of ash (J/(g·K)) 0.976 + lnT +   279 T   1.094 · 10 5 T 2 +   1.169 · 10 7 T 3 [31]
Specific heat capacity of slag (J/(g·K))1.15[32]
Table 7. Simulation results and comparison of previous models with experimental data (Dry based).
Table 7. Simulation results and comparison of previous models with experimental data (Dry based).
Major VariablesThis WorkWen [16]Govind [17]Experiment [17]
Conversion (%)98.898.998.198.6
H2 (mole %)38.639.840.039.1
CO (mole %)57.756.655.557.6
H2/CO ratio0.6680.7030.7210.679
CO2 (mole %)2.562.923.952.95

Share and Cite

MDPI and ACS Style

Kim, J.; Oh, H.; Lee, S.; Yoon, Y.-S. Advanced One-Dimensional Entrained-Flow Gasifier Model Considering Melting Phenomenon of Ash. Energies 2018, 11, 1015. https://doi.org/10.3390/en11041015

AMA Style

Kim J, Oh H, Lee S, Yoon Y-S. Advanced One-Dimensional Entrained-Flow Gasifier Model Considering Melting Phenomenon of Ash. Energies. 2018; 11(4):1015. https://doi.org/10.3390/en11041015

Chicago/Turabian Style

Kim, Jinsu, Hyunmin Oh, Seokyoung Lee, and Young-Seek Yoon. 2018. "Advanced One-Dimensional Entrained-Flow Gasifier Model Considering Melting Phenomenon of Ash" Energies 11, no. 4: 1015. https://doi.org/10.3390/en11041015

APA Style

Kim, J., Oh, H., Lee, S., & Yoon, Y. -S. (2018). Advanced One-Dimensional Entrained-Flow Gasifier Model Considering Melting Phenomenon of Ash. Energies, 11(4), 1015. https://doi.org/10.3390/en11041015

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop