1 s2.0 S1876610212010946 Main
1 s2.0 S1876610212010946 Main
com
Abstract
During the oxy-fuel combustion process pulverized coal is burnt in an atmosphere consisting of pure O2 mixed with
recycled flue gas whereas during the conventional process air serves as the only oxidant. This entails specific conditions
regarding thermo-physical properties which impact both combustion characteristics and heat transfer. Accordingly,
adjustments within CFD codes are required in order to maintain accuracy and prediction quality criteria within simula-
tions of oxy-coal combustion. The CFD code AIOLOS was used to evaluate recent oxy-coal specific implementations
concerning the global chemistry mechanism and the heat transfer. For validation purposes extensive tests have been
carried out at IFK’s semi-industrial scale furnace (500 kWth ). Simulations have been performed for both, conventional
air-firing and oxy-coal combustion conditions with US bituminous coal, and a comparison of simulation results and
corresponding experimental data is given. In general, satisfactory agreement is observed.
©
c 2012
2011The Authors.byPublished
Published Elsevier by Elsevier Ltd. Selection and/or peer-review under responsibility of SINTEF Energi AS
Ltd.
Open access under CC BY-NC-ND license.
Keywords: oxy-fuel process, pulverized coal combustion, CFD modelling, validation
Nomenclature
t time [s]
uj j-th component of the velocity vector [m/s]
xj j-th component of the position vector [m]
SΦ source term of variable Φ [var.]
ΓΦ generalized diffusion coefficient kg/m/s
Φ generalized conservation variable [var.]
ρ density kg/m3
∗ Corresponding author
Email address: (Michael Müller)
1876-6102 © 2012 The Authors. Published by Elsevier Ltd. Selection and/or peer-review under responsibility of SINTEF Energi AS
Open access under CC BY-NC-ND license. doi:10.1016/j.egypro.2012.06.054
198 Michael Mu_ller et al. / Energy Procedia 23 (2012) 197 – 206
1. Introduction
Recent energy statistics reveal that fossil fuels in general and specifically coal are main contributors
to the worldwide energy consumption [1]. Furthermore, rapidly increasing energy demands are expected
due to emerging economies in developing countries. Hence, innovative solutions are required to reduce
greenhouse gas emissions resulting from fossil fuel utilization. The oxy-fuel combustion is seen as one of
the most promising options for CO2 capture to be realized in short-term. During the oxy-fuel combustion
process pulverized coal is burnt in a mixture of pure O2 and recycled flue gas. Flue gas recycling is required
to moderate flame temperatures. This leads to an exhaust gas with high CO2 concentration which is ready for
storage after further conditioning. A more detailed illustration of the oxy-coal process within the subsequent
steps in the chain of CO2 capture and storage (CCS) is shown in Fig. 1. Due to the additional auxiliary power
required for the O2 generation and the purification and liquefaction of the exhaust gas, the net efficiency will
be lowered by about 7 to 11 % [2]. Nonetheless, the oxy-coal combustion process is regarded as a relatively
cost effective CCS method [3, 4] and existing conventional power plants may also be retro-fitted to oxy-coal
operation without major complications.
Fig. 1: Simplified schematic of the oxy-coal process with its main stages: air separation, combustion with
flue gas recycling, exhaust gas conditioning, and CO2 separation (adapted from [2])
The use of simulation tools such as Computational Fluid Dynamics (CFD) leads to a detailed fundamen-
tal understanding of the complex processes in coal combustion. Thus, mathematical modelling provides an
important instrument for future development and design of the oxy-coal combustion technology. However,
the specific conditions of the oxy-coal process result in several major changes concerning thermo-physical
properties which influence both combustion characteristics and heat transfer as well as emission behaviour
[5–9]. Given that, adjustments in the AIOLOS code, which has been developed at the IFK, were required in
order to extend the simulation capabilities of the AIOLOS code for oxy-coal operation as well as conven-
tional air-firing operation. Details about the enhanced models are given in the following section focusing on
the modelling of the coal conversion process and the heat transfer.
The objective of this study was to evaluate the accuracy and prediction quality of the enhanced modelling
approach as well as the impact of the operating conditions. For validation purposes extensive tests have
been carried out at IFK’s semi-industrial scale furnace (500 kWth ) which has been retro-fitted for oxy-coal
combustion. Detailed experimental data were compiled by in-flame measurements as well as continuous
exhaust gas monitoring. Simulations have been performed for both, conventional air-firing and oxy-coal
combustion conditions. The results are presented by comparison of simulation results with corresponding
experimental data highlighting one test case for each operation mode.
Michael Mu_ller et al. / Energy Procedia 23 (2012) 197 – 206 199
Modelling of turbulent reacting flows, and in particular of pulverized coal combustion systems, has
to account for fluid flow, chemical reactions, and heat transfer phenomena. The mutual influence of the
involved physico-chemical processes requires a simultaneous numerical solution of a system of strongly
coupled differential equations, i.e. the governing equations of mass, momentum, energy, and the mean
mass fraction of each species which participates in the considered chemical reactions. The general transport
equation in direction j may be formulated as
∂(ρΦ) ∂(ρu j Φ) ∂ ∂Φ
+ = ΓΦ + SΦ (1)
∂t ∂x j ∂x j ∂x j
with ρ, t, u, x, ΓΦ , and S Φ denoting density, time, velocity, coordinate, diffusion coefficient, and source term,
respectively. Equation (1) describes the local change of the Favre-averaged variable Φ with the correspond-
ing transient, convective, diffusive, and source/sink terms. Furthermore, additional sub-models are required
to account for the broad range of physical and chemical processes within the combustion system.
Simulations have been performed using the CFD code AIOLOS which is based on the Finite Volume
method. Incompressible flow and stationary boundary conditions are assumed. The two-phase flow is
treated as a simplified Eulerian quasi-one-phase flow neglecting slip between gas and dispersed particle
phase. Pressure-velocity coupling is modelled by the SIMPLE method in combination with the interpolation
scheme from Date [10] for pressure correction. The standard k-ε model [11] is applied to describe turbulence
phenomena and the Eddy Dissipation Concept [12] accounts for turbulence-chemistry interactions. The
Discrete Ordinates method is utilized for modelling radiative heat transfer. In addition, the code has been
optimized for vectorization and parallelization and allows domain decomposition. Hence, efficient use of
high performance vector platforms has been enabled. More information concerning the AIOLOS code is
given elsewhere [13–15].
Further details of the developed chemical reaction scheme and the radiation model are illustrated in the
following sections. The coal combustion process can be split into a sequence of stages: devolatilization,
gasification and combustion of the remaining char, and combustion of the volatiles released during pyrol-
ysis. Due to the number of reactions involved, the modelling approach within a CFD framework has to be
simplified in order to maintain reasonable computational effort.
H + O2 O + OH . (R.1)
200 Michael Mu_ller et al. / Energy Procedia 23 (2012) 197 – 206
H + CO2 CO + OH (R.2)
is of similar relevance due to the specific CO2 -rich atmosphere [17]. Reactions (R.1) and (R.2) indicate
that in oxy-coal systems CO2 competes with O2 for the available H-radicals, which may cause a reduction
of flame speed and lower availability of O-radicals. Furthermore, high CO2 levels locally promote the
formation of CO in fuel-rich regions via reaction (R.2). Because detailed chemistry models which are
capable of modelling radical reactions are computationally prohibitive within engineering applications, the
standard global gas phase reaction model has to be extended to consider the specific oxy-coal conditions
[18].
In order to account for the chemical effects of high CO2 concentrations in the oxidizing atmosphere, the
homogeneous water-gas-shift reaction
H2 O + CO CO2 + H2 (R.3)
needs to be incorporated in the global reaction scheme as an equilibrium reaction [18]. The combustion of
light hydrocarbons Cn Hm , which are formed during the devolatilization stage, is modelled as
H2 + 0.5 O2 H2 O . (R.5)
This assumption has proven to be very important for high temperature flames with elevated O2 concentra-
tions which may occur under oxy-coal conditions, since the chemical equilibrium of reaction (R.5) is shifted
towards the educt side with temperatures exceeding 2000 K. Neglecting the reverse reaction would then lead
to a local over-prediction of the flame temperature. The calculation of the respective equilibrium constants
as a function of standard Gibbs free energy is achieved by polynomial fitting to the respective data from
JANAF tables [19]. The corresponding rate expressions and kinetic parameters of reactions (R.3)-(R.5) are
compiled in [18].
(1 + f ) C + O2 → (1 − f ) CO2 + 2 f CO (R.6)
C + CO2 → 2 CO (R.7)
Michael Mu_ller et al. / Energy Procedia 23 (2012) 197 – 206 201
C + H2 O → CO + H2 (R.8)
with f denoting a statistical mechanism factor within the range of 0 to 1. At combustion temperatures above
1000 ◦C f may be presumed to be constant with the main product of reaction (R.6) being CO, resulting in:
f ≈ 1 [20]. The Boudouard reaction (R.7) and the heterogeneous water-gas-shift reaction (R.8) are mostly
neglected in air combustion simulations. But in O2 -lean regions both reactions may have major impact
because partial pressures of CO2 and H2 O are generally higher in oxy-coal combustion systems, the latter
especially in case of wet recycle systems. In general, both reactions should be considered as equilibrium
reactions since product inhibition may occur [21]. Yet, regarding the equilibrium constants of each reaction
reveals that – at typical combustion temperatures and ambient pressure – the equilibrium is shifted towards
the product side. On this basis, reactions (R.7) and (R.8) may be considered irreversible. A comprehensive
analysis of char reaction models under oxy-coal combustion conditions including details about modelling
the char morphology is given in [18].
Fig. 2: Details of the swirl burner with inlets for air and oxy-coal operation (left: light grey denotes air;
right: light / dark grey denote recycled flue gas / pure oxygen).
generator is integrated in the outer annular cross-section. Due to the cylindrical shape of the furnace and the
burner geometry, the flame may be assumed to be axis-symmetric. The burner has been designed specifically
for oxy-coal combustion and because each of the inlet flows has a separate control system, this burner design
allows high flexibility concerning distribution and composition of the individual volume flows at the inlet.
The total inlet flow is divided into a concentric orifice where the pulverized coal is transported by a primary
stream, and four separate secondary streams with the compositions depending on the operation mode (see
Fig. 2). Bottled CO2 and air is used as carrier gas in the primary stream for oxy-coal and conventional air
operation, respectively. Similarly, pure CO2 or air serves as purge gas of the integrated flame detectors for
the different operating conditions.
Fuel properties are compiled in Table 3. The approximated particle size distribution of the pulverized
US bituminous coal is derived from sieve analysis resulting in ten discrete particle size classes. This exper-
imental setup has been operated successfully in a previous study with similar operating conditions utilizing
Lausitz lignite as fuel [29, 30].
Fig. 3: Axial profiles along the furnace centerline showing both benchmark test cases (oxy-fuel and air
operation).
(a) Level 2 @ 0.18 m from the burner (b) Level 2 @ 0.18 m from the burner
(c) Level 3 @ 0.33 m from the burner (d) Level 3 @ 0.33 m from the burner
Fig. 4: Radial profiles at furnace levels 2 and 3 showing the gas temperature and species concentrations of
the oxy-fuel case.
Michael Mu_ller et al. / Energy Procedia 23 (2012) 197 – 206 205
4. Conclusions
An efficient computational modelling framework for oxy-coal combustion was developed and has been
implemented into the CFD code AIOLOS. In order to account for the particular conditions of oxy-coal
combustion, adjustments within certain sub-models were carried out covering mainly reaction mechanisms
and gas radiation. Validation simulations were performed against benchmark tests at IFK’s 500 kWth facility
considering an oxy-fuel case as well as a conventional air-firing case.
The comparison of simulation results and experimental data reveals that the implemented models are
applicable for both operation modes, oxy-coal and air-firing. In general, the calculations agree quite well
with the measurements for both benchmark cases. Nevertheless, some deviations are identified – especially
in the near burner region where the simulations predict the ignition zone further downstream than detected
in the corresponding experiments. Those deviations appear to be caused by more general issues: on the
one hand regarding the modelling of the mixing behaviour and turbulence, and on the other hand regarding
partially uncertain boundary conditions. Moreover, additional validation simulations at various scales are
essential to ensure the reliability of the developed modelling framework.
Acknowledgement
The authors would like to thank ALSTOM Power Systems GmbH, EnBW Kraftwerke AG, and E.ON
Energie AG for funding the presented work. Computational resources have been provided by the High
Performance Computing Center Stuttgart (HLRS).
References
[1] IEA, World Energy Outlook 2010, Organization for Economic Co-operation and Development OECD, ISBN 978-9264086241,
2010.
[2] A. Kather, G. Scheffknecht, The oxycoal process with cryogenic oxygen supply, Naturwissenschaften 96 (2009) 993–1010.
[3] D. Singh, E. Croiset, P. L. Douglas, M. A. Douglas, Techno-economic study of CO2 capture from an existing coal-fired power
plant: MEA scrubbing vs. O2 /CO2 recycle combustion, Energy Conversion and Management 44 (19) (2003) 3073–3091.
[4] T. F. Wall, Combustion processes for carbon capture, Proceedings of the Combustion Institute 31 (2007) 31–47.
[5] B. J. P. Buhre, L. K. Elliott, C. D. Sheng, R. P. Gupta, T. F. Wall, Oxy-fuel combustion technology for coal-fired power generation,
Progress in Energy and Combustion Science 31 (2005) 283–307.
[6] S. P. Khare, T. F. Wall, A. Z. Farida, Y. Liu, B. Moghtaderi, R. P. Gupta, Factors influencing the ignition of flames from air-fired
swirl pf burners retrofitted to oxy-fuel, Fuel 87 (2008) 1042–1049.
[7] M. B. Toftegaard, J. Brix, P. A. Jensen, P. Glarborg, A. D. Jensen, Oxy-fuel combustion of solid fuels, Progress in Energy and
Combustion Science 36 (2010) 581–625.
[8] R. Stanger, T. Wall, Sulphur impacts during pulverised coal combustion in oxy-fuel technology for carbon capture and storage,
Progress in Energy and Combustion Science 37 (2011) 69–88.
[9] G. Scheffknecht, L. Al-Makhadmeh, U. Schnell, J. Maier, Oxy-fuel coal combustion—A review of the current state-of-the-art,
International Journal of Greenhouse Gas Control 5 (2011) S16–S35.
[10] A. W. Date, Complete pressure correction algorithm for the solution of compressible navier-stokes equations on a nonstaggered
grid, Numerical Heat Transfer, Part B: Fundamentals 29 (4) (1996) 441–458.
[11] B. E. Launder, D. B. Spalding, The Numerical Computation of Turbulent Flows, Computer Methods in Applied Mechanics and
Engineering 3 (2) (1974) 269–289.
[12] I. R. Gran, B. F. Magnussen, A Numerical Study of a Bluff-Body Stabilized Diffusion Flame. Part 2. Influence of Combustion
Modeling And Finite-Rate Chemistry, Combustion Science and Technology 119 (1) (1996) 191–217.
[13] U. Schnell, Numerical Modelling of Soild Fuel Combustion Processes using advanced CFD-based Simulation Tools, Progress in
Computational Fluid Dynamics 1 (4) (2001) 208–218.
[14] D. Förtsch, A Kinetic Model of Pulverised Coal Combustion for Computational Fluid Dynamics, Ph.D. thesis, Universität
Stuttgart, Stuttgart, 2003.
[15] J. Ströhle, Spectral Modelling of Radiative Heat Transfer in Industrial Furnaces, Energietechnik, Shaker Verlag, Aachen, 2003.
[16] L. Al-Makhadmeh, Coal Pyrolysis and Char Combustion Under Oxy-Fuel Conditions, Energietechnik, Shaker Verlag, Aachen,
2009.
[17] P. Glarborg, L. L. B. Bentzen, Chemical Effects of a High CO Concentration in Oxy-Fuel Combustion of Methane, Energy &
Fuels 22 (1) (2008) 291–296.
[18] S. Leiser, Numerical Simulation of Oxy-Fuel Combustion, Energietechnik, Shaker Verlag, Aachen, 2011.
[19] E. W. Lemmon, M. O. McLinden, D. G. Friend, NIST Chemistry WebBook, 2009.
[20] J. R. Arthur, Reactions between Carbon and Oxygen, Transactions of the Faraday Society 47 (1951) 164–178.
206 Michael Mu_ller et al. / Energy Procedia 23 (2012) 197 – 206
[21] D. G. Roberts, D. J. Harris, Char gasification in mixtures of CO2 and H2 O: Competition and inhibition, Fuel 86 (2007)
2672–2678.
[22] T. Wall, Y. Liu, C. Spero, L. Elliott, S. Khare, R. Rathnam, F. Zeenathal, B. Moghtaderi, B. Buhre, C. Sheng, R. Gupta, T. Yamada,
K. Makino, J. Yu, An overview on oxyfuel coal combustion – State of the art research and technology development, Chemical
Engineering Research and Design 87 (2009) 1003–1016.
[23] H. C. Hottel, A. F. Sarofim, Radiative Heat Transfer, McGraw-Hill, New York, 1967.
[24] R. Johansson, K. Andersson, B. Leckner, H. Thunman, Models for gaseous radiative heat transfer applied to oxy-fuel conditions
in boilers, Journal of Heat and Mass Transfer 53 (2010) 220–230.
[25] R. Johansson, B. Leckner, K. Andersson, F. Johnsson, Account for variations in the H2 O to CO2 molar ratio when modelling
gaseous radiative heat transfer with the weighted-sum-of-grey-gases model, Combustion and Flame 158 (2011) 893–901.
[26] B. Leckner, Spectral and total emissivity of water vapor and carbon dioxide, Combustion and Flame 19 (1) (1972) 33–48.
[27] N. Lallement, A. Sayre, R. Weber, Evaluation of emissivity correlations for H2 O-CO2 -N2 /air mixtures and coupling with solution
methods of the radiative transfer equation, Progress in Energy and Combustion Science 22 (6) (1996) 543–574.
[28] P. Mönckert, D. Reber, J. Maier, G. Scheffknecht, Operation of a Retrofitted 0.5 MWth PF Combustion Facility under Oxyfuel
Conditions, in: Clearwater Clean Coal Conference, 2007.
[29] S. Grathwohl, O. Lemp, U. Schnell, J. Maier, G. Scheffknecht, F. Kluger, B. Krohmer, P. Mönckert, G. N. Stamatelopoulos,
Highly Flexible Burner Concept for Oxyfuel Combustion, in: 1st Oxyfuel Combustion Conference, 2009.
[30] M. Müller, O. Lemp, S. Leiser, U. Schnell, S. Grathwohl, J. Maier, G. Scheffknecht, F. Kluger, P. Mönckert, Advanced Modeling
of Pulverized Coal Combustion Under Oxy-Fuel Conditions, in: Clearwater Clean Coal Conference, 2010.