Fundamentals of Thermal Radiation For Energy Utilization in Fuel Combustion
Fundamentals of Thermal Radiation For Energy Utilization in Fuel Combustion
Fundamentals of Thermal Radiation For Energy Utilization in Fuel Combustion
Fundamentals of
Thermal Radiation for
Energy Utilization in
Fuel Combustion
Advanced Topics in Science and Technology
in China
Volume 67
Zhejiang University is one of the leading universities in China. In Advanced Topics
in Science and Technology in China, Zhejiang University Press and Springer jointly
publish monographs by Chinese scholars and professors, as well as invited authors
and editors from abroad who are outstanding experts and scholars in their fields.
This series will be of interest to researchers, lecturers, and graduate students alike.
Advanced Topics in Science and Technology in China aims to present the latest
and most cutting-edge theories, techniques, and methodologies in various research
areas in China. It covers all disciplines in the fields of natural science and
technology, including but not limited to, computer science, materials science, the
life sciences, engineering, environmental sciences, mathematics, and physics.
This book series is indexed by both the Scopus and Compendex databases.
If you are interested in publishing your book in the series, please contact Violetta Xu
(Email: [email protected]) and Mengchu Huang (Email: mengchu.huang
@cn.springer.com).
Shiquan Shan · Zhijun Zhou · Yanwei Zhang
Fundamentals of Thermal
Radiation for Energy
Utilization in Fuel
Combustion
Shiquan Shan Zhijun Zhou
College of Energy Engineering College of Energy Engineering
Zhejiang University Zhejiang University
Hangzhou, Zhejiang, China Hangzhou, Zhejiang, China
Yanwei Zhang
College of Energy Engineering
Zhejiang University
Hangzhou, Zhejiang, China
This Springer imprint is published by the registered company Springer Nature Singapore Pte Ltd.
The registered company address is: 152 Beach Road, #21-01/04 Gateway East, Singapore 189721,
Singapore
Preface
At present, fuel combustion is still the significant way for energy utilization. A large
amount of thermal energy released in combustion is converted to electrical energy
through various thermal power cycles. Electrical energy drives related devices to
obtain work, heat, light, sound, and other human terminal needs. The process is
a mode with direct combustion plus a simple thermal power cycle. However, the
single thermal energy conversion neglects the matching of various energy qualities,
resulting in low-efficiency power generation. The energy conversion process involves
the energy quantity and energy-carrying form conversion. The energy-carrying forms
can be photons, electrons, and phonons. For example, the radiation energy, with
photons as carriers, is ignored in the traditional energy conversion. The combus-
tion process releases photon and thermal energy simultaneously, but the energy
quality differs. In high-temperature combustion, the cascade conversion of photon
and thermal energy aligns with the principle of matching energy quality and poten-
tial in the energy conversion process. This matching of energy quality and potential
can theoretically break through individual thermal energy conversion limitations,
effectively improving energy conversion efficiency.
Based on the academic idea of the photo-thermal energy cascade conversion in
fuel combustion, it is crucial to discuss the release characteristics of photo-thermal
energy in combustion. Among them, photo energy refers to the radiative energy which
transfers in the form of electromagnetic waves during combustion process. In heat
transfer field, there is a relatively complete system on the research of thermal radia-
tion, and many radiative characteristic models and calculation methods for radiation
transfer have been proposed. However, few studies on thermodynamics and thermal
radiation exist from the perspective of direct utilization of thermal radiative energy
in fuel combustion. This involves the theory of energy quality for spectral radiation,
as well as the transfer and prediction theory of radiative energy with different qual-
ities under complex combustion reaction conditions. Therefore, it is an extremely
valuable and interesting topic that studying thermal radiation from the perspective
of energy utilization.
v
vi Preface
When I studied for my doctorate, my supervisors, Prof. Zhijun Zhou and Prof.
Yanwei Zhang, provided vital inspiration on combustion thermal radiation and high-
lighted the significance of radiation energy use in combustion, which helped me start
my scientific research career and form the academic idea of photo-thermal energy
cascade conversion. In recent years, several innovative results obtained in this field
have made me fortunate to receive the Outstanding Doctoral Dissertation Award of
Zhejiang University and be selected for the National Postdoctoral Innovative Talent
Support Program of China. Therefore, I plan to summarize our scientific research
results and academic thoughts and write a monograph for academic dissemination.
Thus, this book, Fundamentals of Thermal Radiation for Energy Utilization in Fuel
Combustion, was published. This monograph focuses on the thermal radiation in
fuel combustion based on the theoretical ideas of photo-thermal energy cascade
conversion. The monograph introduces thermal radiation theories developed by the
author’s group based on radiative energy use, and a theoretical system is initially
formed, providing a basis for subsequent research and inspiration for scholars in
related fields.
While drafting this book, Prof. Zhang and Prof. Zhou reviewed the full text and
provided valuable suggestions. Moreover, during my Ph.D. study, I was fortunate
to be sponsored by the China Scholarship Council to study at the Georgia Insti-
tute of Technology with Prof. Zhuomin Zhang, an expert in thermal radiation and
heat transfer, and benefited a lot. Moreover, Prof. Kefa Cen, the academician of
Chinese Academy of Engineering, gives key guidance and strongly supports to this
project. I would like to express my gratitude to them. The National Postdoctoral
Innovative Talent Support Program of China (BX2021254) and the China Postdoc-
toral Science Foundation (2021M702793) primarily funded this monograph. Besides,
the Ningxia Provincial Key Research and Development Program (2018BCE01004)
partially funded it. Moreover, three graduate students in our group, Ms. Siqi Jia, Ms.
Qiaoli Li, and Ms. Jialu Tian, translated and proofread some manuscripts. I express
my gratitude to them all.
1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.1 Photo-Thermal Energy Cascade Conversion for Fuel
Combustion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.2 Key Issues of Combustion Thermal Radiation for Energy
Utilization . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
1.2.1 Theory of Radiation Energy Quality Theory for Fuel
Combustion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
1.2.2 Thermal Radiation Model of Non-gray Gas Medium
in Combustion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.2.3 Radiative Energy Characteristics of High-Temperature
Oxy-Fuel Combustion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
1.3 Objective and Organization of This Book . . . . . . . . . . . . . . . . . . . . . . 9
References . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
2 Spectral Radiation Thermodynamic Theory for Combustion . . . . . . . 15
2.1 Exergy of Blackbody Radiation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
2.1.1 Several Representative Views . . . . . . . . . . . . . . . . . . . . . . . . . . 15
2.1.2 Discussion on the Reversibility . . . . . . . . . . . . . . . . . . . . . . . . . 16
2.1.3 Differences Among the Three Formulas . . . . . . . . . . . . . . . . . 18
2.1.4 Radiation Engine Model for Discussing the Exergy
of Blackbody Radiation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
2.2 Exergy of Monochromatic Radiation . . . . . . . . . . . . . . . . . . . . . . . . . . 21
2.2.1 Discussion of the Exergy of Monochromatic Photons . . . . . . 21
2.2.2 Discussion About the Formula of Equivalent
Temperature . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
2.2.3 Discussion About Radiation from the Combustion
Flame . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
2.3 Entropy of Spectral Radiation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
2.4 Summary . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
References . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
vii
viii Contents
Appendix . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 107
Chapter 1
Introduction
The primary way humans use energy is through the generation of thermal energy
through fuel combustion. Through high-temperature fuel combustion, a large amount
of thermal energy is released, which is then used to generate electric energy
through various thermal power cycles. However, the energy production process based
on the thermal power cycle does not consider energy quality matching, and the
energy released by fuel combustion is only utilized in the form of thermal energy.
However, the highest parameters of the thermal power cycle limit the system’s energy
conversion efficiency.
During the combustion process, the chemical energy in the fuel is converted into
photo energy and thermal energy. The term photo energy refers to both visible and
invisible light (e.g., ultraviolet and infrared) [1], i.e., radiative energy. The heat energy
mentioned here primarily refers to non-radiative heat energy transferred during
the combustion process. From a thermodynamic standpoint, the energy released
by combustion is not limited to having both a quantitative and qualitative char-
acter. Combustion produces both photo and thermal energy, but their qualities differ.
Additionally, photo energy has spectral characteristics, and spectral radiation with
different wavelengths has different energy quality. Short-wavelength radiation has
high energy quality and can directly excite the photovoltaic conversion device. Thus,
carrying out cascade and orderly conversion of photo-thermal energy of various qual-
ities is necessary: the high-grade radiative energy from high-temperature combustion
should be converted step-by-step into electrical energy following the frequency, and
the low-grade thermal energy should be further utilized through thermal power gener-
ation to establish a new photo-thermal energy cascade conversion system (Fig. 1.1),
which will greatly improve energy conversion efficiency compared with traditional
thermal power cycle [2, 3].
Flame radiation results from high-temperature combustion, and high-temperature
energy has higher quality. Thus, the technical route of cascade conversion of photo-
thermal energy conforms to thermodynamic principles and is an innovative idea to
overcome the limitation of the efficiency of single thermal power conversion. Melino
et al. [4, 5] used the radiative energy generated by the high combustion temperature
to power photovoltaic conversion, and the remaining thermal energy in flue gas is
used to power the Organic Rankine Cycle. Filters can also control and shape the
Fig. 1.2 The photo-thermal cascade conversion system that directly utilizes the spectral radiative
energy of combustion flame [3]. Reprinted from Shan et al. [3], Copyright (2019), with permission
from John Wiley and Sons
4 1 Introduction
than conventional air combustion under the conditions of increased oxygen concen-
tration or even directly using pure oxygen [8]. If the thermal power cycle will be
continued to be directly used for electricity generation in high-temperature oxy-fuel
combustion, the high energy quality will be wasted due to the thermal power cycle’s
maximum parameters. Furthermore, the high-temperature of oxy-fuel combustion
generates a significant amount of radiative energy. Thus, a photo-thermal energy
cascade conversion system based on oxy-fuel combustion can be built by prioritizing
radiative energy cascade conversion, followed by thermal power cycle conversion.
It not only controls CO2 emissions in fuel combustion, in line with the strategic
goal of carbon neutrality, but it also effectively matches energy quality and potential,
theoretically improving energy conversion efficiency.
Given the preceding concept and the prevalence of energy conversion through
fuel combustion, the development of combustion thermal radiation theory based on
the background of photo-thermal energy cascade utilization is not only important
for guiding the combustion-based photo-thermal energy cascade conversion system
but can also enrich the theoretical system of thermal radiation. These primarily
include the radiative energy quality theory, an efficient radiation calculation model
for complex combustion conditions, and the mechanism of radiative energy release
and distribution in the combustion process. Section 1.2 will go over these in greater
detail.
Radiative energy in fuel combustion has spectral characteristics, and spectral radi-
ation with different wavelengths have different energy qualities. To use radiative
energy more efficiently, analyzing its quality from the standpoint of thermodynamics
is necessary, thus quantitatively evaluating the value of radiative energy and guiding
the cascade conversion of high-temperature radiative energy. Exergy is used in ther-
modynamics to describe the quality of energy, which refers to the usable part, also
known as available energy. When compared with energy analysis, exergy analysis
can describe not only the quantity but also the quality of energy. Consequently, the
analysis of thermal radiation processes using the second law of thermodynamics
is of unique significance. On the one hand, radiation is the primary heat transfer
mode in high-temperature equipment, and second law analysis of radiation aids in
investigating the entropy production and exergy loss of high-temperature combustion
equipment, allowing the high-temperature energy conversion system to be evaluated
and optimized more rationally. The second law, on the other hand, can be used to
accurately analyze exergy in radiation transfer process within different wavebands,
1.2 Key Issues of Combustion Thermal Radiation for Energy Utilization 5
which further guides energy management with different energy grades to achieve
efficient exergy conversion.
Radiation, as an electromagnetic wave, has wavelength and frequency. Moreover,
as a way of heat transfer, it has entropy and exergy properties similar to thermal
energy. Radiative energy has its distinct characteristics, which distinguishes it from
thermal energy. Wright [9] identified the difference between the entropy of radiation
heat transfer and the entropy of heat conduction processes. Liu and Chu [10] inves-
tigated the expression of entropy production of the heat transfer process and discov-
ered that the conventional formula of entropy production rate for heat conduction is
ineffective for radiative heat transfer.
The blackbody radiation exergy theory proposed by Petela at the University of
Silesia in Poland in his early years [11], Spanner’s theory based on solar radiation
research [12], and the Carnot efficiency theory proposed by Jeter of Georgia Institute
of Technology [13] are representative views for the exergy of blackbody radiation.
The theoretical core of the three is the coefficient of radiation exergy-to-exergy ratio,
which is summarized in Table 1.1. Bejan [14] proposed some thought-provoking
discussions on all three, which will be introduced and discussed in depth in Chap. 2
of this book. Currently, scholars generally accept Petela’s blackbody radiation exergy
formula. Candau [15] also proposed a derivation of the intensity of spectral radiation
exergy based on Planck’s formula for entropy intensity of monochromatic radiation
and the concept of radiation temperature for spectral radiation exergy. Based on
the concept of equivalent temperature, Chen et al. [16] proposed a corresponding
calculation formula for monochromatic photon exergy. These studies ignored the
distinction between radiative and thermal energy.
Lior [17] previously reviewed the exergy calculation in the radiation transfer
process, but the radiative exergy problem was rarely considered in early studies.
Caldas and Semiao [18] were the first to consider entropy production in radiative
transfer calculations and they separately analyzed entropy production by emission,
absorption, and scattering processes in the media radiative transfer. Liu and Chu [19]
investigated the radiative exergy transfer in the media on this basis. These studies’
models are compatible with a standard method for solving radiative transfer equa-
tions. Based on these hypotheses, Zhang et al. [20, 21] investigated various entropy
generations, including thermal radiation in conditions (e.g., mixed flue gas and hydro-
carbon fuel combustion) and discovered that soot has a greater effect on radiative
entropy generation. These studies used the radiation temperature concept [22] and
expressed the spectral radiative entropy using a formula similar to the Carnot thermal
efficiency, which is intended for calculating the entropy generation of the radiation
heat transfer. In fact, the exergy of thermal energy in the form of Carnot’s law is not
suitable for characterizing radiative entropy and exergy, and the radiation tempera-
ture is a relatively macroscopic concept that is not used to directly characterize the
quality of monochromatic radiation energy. Thus, the monochromatic spectral radi-
ation itself must be considered to better characterize the monochromatic radiation
exergy.
model is typically 5–50 cm−1 ), with the typical narrow band models being a statis-
tical narrowband model (SNB) [25] and statistical narrowband correlated-k model
(SNBCK) [26]. Although the SNB model is computationally intensive, results from
the SNB model based on the EM2C Lab database are very similar to those of the
accurate LBL model based on the HITEMP2010 database, and thus it can be used
as a benchmark [27]. In general, the wideband model interval is 100–1000 cm−1 .
The most common are the exponential wideband (EWB) model [28], the wideband
k-distribution (WBK) model [29], and so on.
The global model is primarily concerned with radiative heat flux, with less
emphasis placed on spectrum characteristics. The global model can significantly
reduce calculation time while also meeting the accuracy requirements for engineering
calculations, making it suitable for CFD calculations. Typical global models include
the full spectrum k-distribution (FSK) model [30], the weighted-sum-of-gray-gases
(WSGG) model [31], and the spectral line-based WSGG model (SLW) [32]. The
basic idea behind the WSGG model is to substitute several gray gases for the actual
non-gray gas. This method’s principle is simple and its calculation efficiency is very
high; it can be combined with various radiation transfer equation solution methods,
giving it a wide range of adaptability. After considering calculation time and accuracy,
the WSGG model is the best gas radiation model for current combustion calculation
and research [33].
Smith et al. [31] developed a set of WSGG model parameters based on the
SNB model for air–fuel combustion conditions. However, the flue gases in oxy-
fuel combustion contain a high concentration of triatomic gases (e.g., CO2 ) which
have different radiation properties than air–fuel combustion. Thus, the traditional
WSGG model applied to air–fuel combustion is inapplicable to oxy-fuel combus-
tion [34]. In recent years, researchers have developed various WSGG models for
oxy-fuel combustion conditions. Based on the EWB model, Khare [35] established
8 1 Introduction
WSGG model parameters for oxy-fuel conditions and discovered that the gas emis-
sivity in oxy-fuel conditions is underestimated when using the traditional WSGG
model parameters in [31]. Yin [36] developed a new WSGG model using the EWB
model as a benchmark. Additionally, Johansson et al. [37] developed a new WSGG
model using the SNB model of the EM2C database as a benchmark. Later, Kang-
wanpongpan et al. [38], Dorigon et al. [39], and Bordbar et al. [40], presented their
respective models based on the more accurate HITEMP2010 database. The models
presented above were all designed for oxy-fuel combustion at atmospheric pres-
sures. Future studies should focus on developing WSGG models for gas radiation
that covers a wider range of pressures, temperatures, and atmospheres for complex
reaction conditions (e.g., pressurized oxy-fuel combustion and combustion in aero
engines).
is believed to contribute more to the O2 /CO2 atmosphere than the known contribu-
tion of elevated partial pressure of CO2 . The presence of CO2 not only increases
gas radiation but also has a significant impact on the formation of soot and its flame
radiation. Following that, Andersson et al. [44] continued to investigate the radiation
characteristics of oxy-fuel combustion flame of lignite on this furnace, discovering
that the flame temperature and total radiation intensity will increase as the oxygen
proportion in the atmosphere increases from 25 to 29%. When oxy-fuel combustion
is used instead of air–fuel combustion, the ratio of gas radiation to total radiation
increases slightly. Thus, the authors concluded that the majority of the radiation in
lignite flames is emitted by flame particles. In recent years, researchers have devel-
oped numerical cases to simulate the physical field in a furnace using the radiation
transfer equation. Gronarz et al. [45] investigated the coupling relationship between
gas radiation and particle radiation in an oxy-fuel combustion boiler for pulverized
coal. The findings show that particle emissivity had a significant impact on radiation
intensity. Adams et al. [46] established a two-dimensional furnace case for pulver-
ized coal combustion and calculated the influence of pressure on the radiative heat
flux of oxy-fuel combustion. The findings indicate that increasing pressure within
10 atm can increase radiative heat flux.
The preceding studies contribute to a better understanding of the radiation char-
acteristics of fuel combustion. The results show that oxy-fuel combustion can effec-
tively improve the proportion of radiation energy in combustion, and the radiative
heat flux can be effectively managed by adjusting the combustion parameters, which
is conducive to efficient radiation energy utilization. The majority of current studies
focus on the effects of fuel type and oxygen concentration on the radiative heat
transfer properties of oxy-fuel combustion. The spectral radiation characteristics
in combustion processes and their influence mechanism are important issues to
investigate in the future in terms of radiation energy utilization.
The goal of this book is to present, for the first time, academic thoughts on photo-
thermal energy cascade conversion of fuel combustion. Based on this, with a focus
on thermal radiation in fuel combustion, this book will primarily introduce a series
of thermal radiation theories for radiative energy utilization on recent studies of the
authors’ group, forming an initial theory system. This serves as the foundation for
future studies while also inspiring scholars in related fields with academic thought
and scientific theory.
The following key scientific problems on combustion thermal radiation must be
solved to achieve the goal of efficient and scientific utilization of high-temperature
combustion energy. (1) No reasonable method exists to evaluate radiation energy
quality, particularly the energy quality theory of monochromatic radiation for fuel
combustion, so expanding the radiation thermodynamics theory is critical. (2) A lack
of an efficient model suitable for a wide parameter range was noted in the calculation
10 1 Introduction
Fig. 1.4 The relationship diagram of the main chapters of this book
References 11
models), and the WSGG model structure is improved to make it compatible with
parameters in a wider pressure range. Furthermore, based on the benchmark of
the SNB model in the EM2C laboratory, a new WSGG model and coefficients
for complex combustion conditions are developed. The new WSGG model and
coefficients can produce accurate and efficient calculation results.
Chapter 4 discusses the calculation method and characteristics of radiation exergy
in participating medium of fuel combustion. Firstly, the theory of monochromatic
radiation exergy proposed in Chap. 2 is combined with the radiation transfer equa-
tion, establishing the monochromatic radiation entropy and radiation exergy transfer
equations and verifying the calculation method based on the second law of thermody-
namics. Following that, a one-dimensional furnace case is simulated to calculate the
spectral radiation energy characteristics of the participating medium in fuel combus-
tion, providing a foundation for numerical calculation of combustion radiation for
energy utilization.
Chapter 5 is the radiation energy characteristics of solid fuel combustion. Based
on the experimental platforms of the tubular furnace and the small flat-flame burner,
it introduces the radiation energy characteristics of oxy-fuel combustion with various
solid fuels, including pulverized coal and semi-coke. The effect of reaction condi-
tions on the radiation energy and spectral properties is investigated. Furthermore, the
radiation exergy and its spectral properties are investigated using the radiation ther-
modynamics theory proposed in Chap. 2. The findings are useful in understanding
and controlling radiative energy in fuel combustion.
References
37. Johansson R, Leckner B, Andersson K et al (2011) Account for variations in the H2 O to CO2
molar ratio when modeling gaseous radiative heat transfer with the weighted-sum-of-gray-gases
model. Combust Flame 158(5):893–901
38. Kangwanpongpan T, França FHR, da Silva RC et al (2012) New correlations for the weighted
sum of gray-gases model in oxy-fuel conditions based on HITEMP 2010 database. Int J Heat
Mass Transf 55(25–26):7419–7433
39. Dorigon LJ, Duciak G, Brittes R et al (2013) WSGG correlations based on HITEMP2010 for
computation of thermal radiation in non-isothermal, non-homogeneous H2 O/CO2 mixtures.
Int J Heat Mass Transf 64:863–873
40. Bordbar MH, W˛ecel G, Hyppänen T (2014) A line by line based weighted sum of gray
gases model for inhomogeneous CO2 –H2 O mixture in oxy-fired combustion. Combust Flame
161(9):2435–2445
41. Nozaki T, Takano S, Kiga T et al (1997) Analysis of the flame formed during oxidation of
pulverized coal by an O2 -CO2 mixture. Energy 22(2–3):199–205
42. Andersson K, Johnsson F (2007) Flame and radiation characteristics of gas-fired O2 /CO2
combustion. Fuel 86(5–6):656–668
43. Andersson K, Johansson R, Johnsson F et al (2008) Radiation intensity of propane-fired oxy-
fuel flames: implications for soot formation. Energy Fuels 22(3):1535–1541
44. Andersson K, Johansson R, Hjärtstam S et al (2008) Radiation intensity of lignite-fired oxy-fuel
flames. Exp Thermal Fluid Sci 33(1):67–76
45. Gronarz T, Schulze J, Laemmerhold M et al (2017) Quantification of the influence of parameters
determining radiative heat transfer in an oxy-fuel operated boiler. Fuel Process Technol 157:76–
89
46. Adams BR, Hosler TR (2019) Pressure and particle property impacts on radiation in oxy-coal
combustion. Fuel 239:667–676
Chapter 2
Spectral Radiation Thermodynamic
Theory for Combustion
Revised from Zhou et al. Exergy of Blackbody Radiation and Monochromatic Photon, International
Journal of Thermophysics, Copyright (2017), with permission from Springer Nature.
Petela [1] proposed the following formula for blackbody radiation exergy at
temperature T in 1964
( )[ ]
4 1 W
bBR = σ T − T0 T 3 + T04
4
(2.1)
3 3 m2
where T 0 is the ambient temperature. Moreover, radiative exergy can be calculated via
multiplying energy by the radiative exergy coefficient, which is defined as the ratio of
exergy to energy. The blackbody radiation intensity of is σ T 4 and the corresponding
radiative exergy coefficient can be expressed as
( )
4 T0 1 T0 4
ηP = 1 − + (2.2)
3T 3 T
This is called Petela’s formula. According to Jeter [2], the radiative exergy coef-
ficient is equal to the Carnot efficiency as shown in Eq. (2.3). Some researchers have
also recognized this [3, 4].
T0
ηJ = 1 − (2.3)
T
Spanner proposes a radiation exergy coefficient based on the study of solar
radiation [5]
4 T0
ηS = 1 − (2.4)
3T
Petela’s study on radiation exergy is earlier and more persistent. Although each
formula is supported by some scholars, Eq. (2.1) is the most widely used formula,
and many researchers have amended or expanded on it [6–11].
In this section, the differences and relationships among the three radiation exergy
formulas are analyzed. At the same radiation temperature T and ambient temperature
T 0 (T > T 0 ), ηS < ηP < ηJ can be observed. The Carnot efficiency is the highest
reversible conversion efficiency for thermal energy. Thus, Jeter’s formula appears to
be reasonable if radiation is regarded as thermal energy. Furthermore, some scholars
believe that the reason why ηP < ηJ is because Petela’s model for calculating ηP is
irreversible. Bejan [12], for example, proposed a special model in 1987 that includes
three reversible processes: first, photon gas at temperature T 1 fulfills the space volume
V 1 ; second, it is connected with a reversible Carnot engine and outputs work to the
outside at temperature T 2 ; and third, the photon gas is exhausted and the volume
2.1 Exergy of Blackbody Radiation 17
T0 ( )
T0
W = 1− dQ (2.5)
T
T
The differences between the formulas of Petela and Spanner are investigated in [15],
pointing out that Spanner the absolute work, instead of useful work, is regarded as
exergy in Spanner’s formula. According to energy conservation law, after absorbing
the energy Q from the outside, the initial internal energy of the medium increases
from U 1 to U 2 and performs the absolute work W 1–2 . Furthermore, the absolute work
W 1–2 includes two parts, which are useful work W u , as well as the work W e spent
for the compression of the environment, where W e = p0 (V 1 − V 0 ). Exergy means
the useful energy, thus, it makes sense that and using useful work as radiative exergy.
Thus, when the exergy is represented by useful work W u , Petela’s theory is deduced.
In solar energy engineering, since the radiation temperature is as high as 5800 K, the
calculation results from Spanner’s theory and Petela’s theory are very similar.
Figure 2.2 is plotted to discuss Jeter’s formula, which is expressed as Carnot effi-
ciency. Moreover, Fig. 2.2a showed that the radiation from a blackbody at a high
temperature of T is directly absorbed by an ideal Carnot engine. This radiation could
do the most work with Carnot efficiency. Nevertheless, thermal energy should be
the input of Carnot engine rather than radiative energy. Otherwise, the process of
radiation will be overlooked and will not be able to obtain the exact radiative exergy.
Thus, it is necessary to set a receiver to convert radiative energy into thermal energy,
and then, the Carnot engine could convert thermal energy into work as shown in
Fig. 2.2b. Here the temperature of thermal energy T ' is always less that the tempera-
ture of blackbody T; otherwise, there is no radiation heat transfer. The Carnot engine
absorbs thermal energy to output work, and its efficiency ηb is 1 − T 0 /T ' , with 1 −
T 0 /T ' < 1 − T 0 /T. This implies that the Carnot efficiency should be greater than the
radiative exergy coefficient. Furthermore, the radiative exergy should be calculated
based on the quasi-equilibrium process as shown in Fig. 2.1, and the result conforms
to Petela’s theory, which is less than the efficiency deduced by Carnot engine. Conse-
quently, Jeter’s theory implies that radiation enters the Carnot engine directly, this is
the difference between Petela’s theory and Jeter’s theory. Thus, it’s more rational to
represent the exergy of blackbody radiation through Petela’s formula.
In this section, a radiation engine model, similar to the Carnot heat engine model,
is presented to discuss the exergy formula of blackbody radiation. The nature of the
radiation could thus be better explained. The radiation engine, as shown in Fig. 2.3,
is an ideal black box converter; the input energy is the blackbody radiative energy
flow Rin , with temperature T and entropy S in . The output energy is the blackbody
radiative energy flow Rout , with temperature T 0 and entropy S out ; it is balanced with
the ambient conditions. The radiation engine’s entire process produces work W,
2.1 Exergy of Blackbody Radiation 19
which is the maximum amount of work. The heat conduction with the environment
is Q, and its entropy is S.
The entropy balance equation can be expressed using the second law of
thermodynamics as follows
where ΔS is the entropy production; for a reversible process with the highest conver-
sion efficiency, ΔS = 0. The input energy flow is expressed as Rin = σ T 4 , the
/ flow3is expressed as Rout = σ T0 , the inlet entropy flow
4
outlet energy / is expressed
as Sin = 4 3σ T , the outlet entropy is expressed as Sout = 4 3σ T03 , and/ the
entropy flow of heat conduction with the environment is expressed as S = Q T0 .
The formula below is obtained when the expression for entropy is substituted into
Eq. (2.7)
4 ( 3 )
Q = (Sin − Sout )T0 = σ T − T03 T0 (2.8)
3
The formula below is obtained by substituting Eq. (2.8) and the expressions for
energy flow into Eq. (2.6)
This is the coefficient from Petela’s exergy formula. When compared to thermal
energy, the entropy of radiation can be expressed as Q/T. Furthermore, the inlet and
outlet entropy flow is Sin = σ T 3 and Sout = σ T03 , respectively. Then, it is substituted
into (2.9) and obtain
W σ T 4 − σ T 3 T0 T0
η= = =1− (2.12)
Rin σT 4 T
Equation (2.12) represents exactly the Carnot efficiency proposed by Jeter; thus,
the reason why the coefficient of Petela’s exergy formula is less than the Carnot
efficiency is not the irreversibility of the radiation system but the difference between
radiation and heat. The distinction implies that the entropy-to-energy ratios of radia-
tion and heat conduction differ. In the case of reversible conversion of radiation energy
conversion, the radiation engine should be left in the form of blackbody radiation at
the ambient temperature. This is consistent with Wright’s [17] viewpoint.
In addition to using classical thermodynamic theory to analyze the radiation
exergy, Badescu [18] demonstrated Petela’s formula using the quantum method of
statistical thermodynamic theory. Additionally, Markvart and Bauer [19] discussed
the exergy of photovoltaics and obtained Petela’s formula when they calculated the
exergy over the entire wavelength range. All of the preceding discussions and anal-
yses indicate that Petela’s formula essentially clarifies the differences between radia-
tive energy and thermal energy. Moreover, it can accurately represent radiative exergy
and is widely accepted by scholars.
T0
ηu = 1 − (2.14)
Tr
v ( )
T0
Wu = h 1 − dv (2.15)
T(v)
v0
f Tm
k
T(λ) = (2.16a)
λn
24 2 Spectral Radiation Thermodynamic Theory for Combustion
In this section, the accuracy of the equivalent temperature formula Eq. (2.16) is
observed and the values of constants m, n, k, and f are discussed. It can be deduced
from Eq. (2.16) that
f 1/ k T m / k n / k
T(v) = v (2.17)
cn / k
T0 T(v0 ) ( v ) nk
0
= = (2.18)
T(v) T(v) v
k k−n n n
Wu = hv − hv k v0k + hv0 . (2.19)
k−n k−n
c1 λ−5
i b,λ = (2.21)
e(c2 / λT ) − 1
/
where c1 = 2π hc02 = 3.74 × 10−16 W m2 , c2 = hc0 k = 1.4388 × 10−2 m K, c1
denotes the first Planck constant, and c2 denotes the second Planck constant. Thus,
the following relationship exists
∞ [ ( ) ] [ ( )4 ]
k T0 n T0 k 4 T0 1 T0
i b,λ 1− + dλ = σ T 4 1 − +
k − n T(λ) k − n T(λ) 3T 3 T
0
(2.22)
Equation (2.22) can be used to calculate equivalent temperature, and the left side
of Eq. (2.22) can be written as
∞ [ ( ) ]
k T0 n T0 k
i b,λ 1− + dλ
k − n T(λ) k − n T(λ)
0
∞ [ ( ) ]
c1 λ−5 k T0 n T0 k
= 1− + dλ
ec2 / λT − 1 k − n T(λ) k − n T(λ)
0
∞ ∞
k c1 T0 λ−5+n / k n c1 T0k λ−5+n
= σT − 4
dλ + dλ
k − n f n/ k T m / k ec2 / λT − 1 k − n f Tm ec2 / λT − 1
0 0
(2.23)
/
Set x = c2 λT ; then, d x = − λc22T dλ. Thus,
∞
k c1 T0 λ−5+n / k
dλ
k − n f n/ k T m / k ec2 / λT −1
0
0 ∞
k c1 T0 T 4−n / k−m / k x 3−n / k k c1 T0 T 4−n / k−m / k x 3−n / k
=− d x = dx
k−n f n/ k c /
4−n k ex − 1 k−n f n/ k c /
4−n k ex − 1
2 ∞ 2 0
(2.24)
and
∞
n c1 T0k λ−5+n
dλ
k − n f Tm ec2 / λT − 1
0
26 2 Spectral Radiation Thermodynamic Theory for Combustion
0 ∞
n c1 T0k T 4−n−m x 3−n n c1 T0k T 4−n−m x 3−n
=− d x = d x (2.25)
k−n f c24−n ex − 1 k−n f c24−n ex − 1
∞ 0
∞ 3−n/ k ∞ 3−n
Make a = 0 xex −1 d x; b = 0 exx −1 d x.
Thus, Eq. (2.22) can be expressed as
k c1 a n c1 b k 4−n−m
σT4 − T0 T 4−n / k−m / k + T T
k − n f n/ k c 4−n / k k − n f c24−n 0
2
4 1
= σ T 4 − σ T0 T 3 + σ T04 (2.26)
3 3
[ ] [ ]
[ 4According to the]theory of dimensional analysis, T0 T 3 = T0 T 4−n / k−m / k and
] [ k 4−n−m
T0 = T0 T exist. Thus, k = 4, m = 3, and n = 1 can be deduced; Eq. (2.16)
can be written as:
f T3
4
T(λ) = (2.27)
λ
From Eq. (2.20), the corresponding coefficient of exergy to energy ηu can be
obtained
( )
4 T0 1 T0 4
ηu = 1 − + (2.28)
3 T(λ) 3 T(λ)
4 c1 a 1 c1 b 4 4 1
σT4 − T0 T 3 + T0 = σ T 4 − σ T0 T 3 + σ T04 (2.29)
3 f 1/ 4 c /
15 4 3
3 f c2 3 3
2
∞ 11/ 4 ∞ 2
where σ is the Stefan–Boltzmann constant and a = 0 exx −1 d x; b = 0 exx−1 d x.
Using MATLAB software to calculate the integral, a = 4.8732 and b = 2.4041 can
be obtained.
Then, Eq. (2.29) is discussed under various conditions. When T > T 0 , the third
term in Eq. (2.29) accounts for a small proportion, and thus the value of exergy is
determined by the second term, i.e., the second term is primarily considered and
c1 a = f 1/ 4 c2 / σ ; the calculated result is f = 4.5638 × 10−3 m K.
15 4
When T < T 0 , the second term accounts for a small proportion of the total, and thus,
the value of exergy is determined by the third term, i.e., the third term is primarily
considered and c1 b = f c23 σ ; and the calculated result is f = 5.33 × 10−3 m K.
In this case, when T < T 0 , the formula can be expressed as 5.33×10−3 T 3 = λT(λ) 4
,
and an error ε2 will be noted after the introduction of Eq. (2.16) for the exergy
calculation of blackbody radiation.
( )
4 c1 a
− σ T0 T 3 0.24
f 1/ 4 c2 /
3 15 4
ε2 = = ( / ) ( / )3 (2.31)
σ T 4 − 43 σ T0 T 3 + 13 σ T04 3 T T0 + T0 T − 4
The ambient temperature T 0 is set at 300 K. The errors ε1 and ε2 are calculated
at various temperatures. The results show that the error introduced by the formula
is very small in the case of high or low temperatures. Furthermore, under common
industrial radiation conditions, with a temperature of ~1000 K, the error is <0.1%;
for a solar radiation temperature of ~5800 K, the error is <10−8 and thus negligible.
The equivalent temperature formula applies to both industrial and solar radiation.
The expression form of the equivalent temperature in Eq. (2.16) in this book is
approximate. However, the integral form of Eq. (2.15) for radiation photon exergy
is still accurate. Furthermore, the preceding analysis shows that the monochromatic
photon exergy of blackbody radiation calculated using the equivalent temperature
expression is applicable in most situations, particularly in high-temperature engi-
neering. Furthermore, the exergy-to-energy coefficient of monochromatic photons
calculated using the equivalent temperature expression is similar to Petela’s coef-
ficient. Thus, it is simple, easy to remember, and can be used to calculate
non-blackbody radiation in most engineering fields.
Equivalent temperature T (λ) can represent the energy quality of monochromatic
radiation photons. Additionally, the higher the value of T (λ) , the better the energy
quality. Thus, the exergy of monochromatic photons can be calculated using equiva-
lent temperature. The exergy produced by converting monochromatic radiation into
electricity, mechanical energy, or other forms of energy is no greater than the exergy
of monochromatic radiation.
Fuel combustion can generate a high temperature of about 2000 K, and for oxy-fuel
combustion conditions, it can generate a high temperature of more than 2500 K.
Based on the formula of equivalent temperature Eq. (2.27), Fig. 2.5 shows the
28 2 Spectral Radiation Thermodynamic Theory for Combustion
4500 0.95
(a) Equivalent temperature (b) Exergy-to-energy coefficient
4000
0.90
3500
2000 K 2000 K
2500 K 0.85 2500 K
3000 K 3000 K
3000
0.80
2500
0.75
2000
1500 0.70
1000 0.65
0 2 4 6 8 10 0 2 4 6 8 10
Fig. 2.5 The monochromatic radiation equivalent temperature T (λ) and exergy-to-energy coeffi-
cient ηu versus the wavelength λ at combustion radiation temperatures
The study of photon entropy in monochromatic radiation has been going on for
a long time. Planck [29] gave statistical average photon entropy in an equilibrium
state with fixed frequency ν and specific blackbody radiation temperature T. Further-
more, Ito [30] believed that the entropy of a monochromatic photon is hν/T for the
photon emitted from a radiation source with temperature T and frequency ν. All
of this indicates that the entropy of monochromatic radiation is related to radia-
tion temperature T and the frequency of the photon ν. On the other hand, Kirwan
[31] investigated the inherent entropy of a single photon using the second law of
thermodynamics, Planck’s equation, and the Gibbs–Duhem equation and discovered
that photon entropy is constant regardless of wavelength. Furthermore, Chen [26]
believed that the entropy of a monochromatic photon is hν/T (ν ) based on the concept
of equivalent temperature T (v) . Additionally, a conclusion that photon entropy is a
2.3 Entropy of Spectral Radiation 29
constant was obtained based on the relationship between equivalent temperature and
photon frequency in this book. Most scholars, however, find it difficult to accept the
notion that photon entropy is a constant. It is perhaps doubtful that an independent
photon is directly given the macroscopic parameter temperature T, as Kirwan’s paper
suggests. Due to the difference between radiative energy and thermal energy, it is also
questionable whether the entropy expression analogous to the Q/T (ν ) form is given
directly using equivalent temperature T (ν ) . Thus, the infinite-staged Carnot engine
model proposed in this chapter is used to analyze the entropy of a photon.
The infinite-staged Carnot heat engine is used to analyze the entropy of a photon
as shown in Fig. 2.6. A monochromatic photon has the energy hν and the equivalent
temperature T (ν ) . The absorbing model absorbs a photon, converting a small portion
of the energy converted into thermal energy dq and emitting the remainder as radiation
with a slightly lower frequency ν ' ; thus, its energy is hν ' . The entropy of the thermal
energy transformed from radiation can be expressed in the form Q/T (v) ; thus, dq/T (v)
= hdv/T (v) . All of the photon’s radiation energy is converted into thermal energy
through infinite processes like this until the frequency reaches zero. Because each
step of the process is reversible and generates no entropy, the entropy of a photon
with frequency ν should be expressed as the sum of thermal energy entropy in each
process.
The expression of monochromatic photon entropy should be
v
h
sv = dv (2.32)
T(v)
0
hc1/ 4 3
sv = v4 (2.33)
f 1/ 4 T 3/ 4
30 2 Spectral Radiation Thermodynamic Theory for Combustion
where c is the speed of light and f is the coefficient in the equivalent temperature equa-
tion. Equation (2.33) should also be approximated because the relationship between
radiation equivalent temperature and frequency in Eq. (2.27) is approximated in this
chapter. However, nothing is wrong with qualitative analysis. As a result of this anal-
ysis, the entropy of a monochromatic photon can be concluded to be proportional to
its frequency ν and radiation temperature T.
Moreover, after substituting the Eq. (2.32) into the formula for the relationship
between exergy and entropy, the expression for monochromatic photon exergy can
be obtained
( )
Wu = H − H0 − T0 (S − S0 ) = hv − hv0 − T0 sv − sv0
⎛ v ⎞
v v0
h h
= h dv − T0 ⎝ dv − dv ⎠
T(v) T(v)
v0 0 0
v v v ( )
h T0
=h dv − T0 dv = h 1 − dv (2.34)
T(v) T(v)
v0 v0 v0
This is just the exergy expression in Eq. (2.15), and it also indicates that the
exergy and entropy of monochromatic radiation photons proposed in this chapter are
intrinsically harmonious.
2.4 Summary
Distinctions exist between radiative energy and thermal energy, and the mechanism
of radiation heat transfer differs from that of heat conduction or convection. Thus,
radiative energy cannot be considered synonymous with thermal energy. In several
views of blackbody radiation exergy, Petela’s formula is reversible in derivation
and distinguishes the differences between radiative energy and thermal energy when
compared to the Carnot efficiency formula proposed by Jeter. Therefore, Petela’s
formula accurately represents radiation exergy.
Based on the concept of monochromatic photon equivalent
v ( temperature
/ ) T (v) , it
is believed that the formula with integral form Wu = v0 h 1 − T0 T(v) dv should
be used to represent the exergy of monochromatic photons, which can reflect the
difference between radiative energy and thermal energy.
The equivalent temperature of monochromatic photon T (λ) is proportional to radi-
ation temperature T and inversely proportional to wavelength λ. The higher the
equivalent temperature T (λ) , the greater the exergy of monochromatic photons.
An approximate formula exists to calculate monochromatic photon equivalent
temperature T (λ) from the radiation temperature T and wavelength λ is given. When T
> T 0 and T < T 0 , 4.5638×10−3 T 3 = λT(λ) 4
and 5.33×10−3 T 3 = λT(λ)
4
, respectively.
Based on this equivalent temperature expression, the formula of exergy-to-energy
References 31
( )4
ratio ηu for a monochromatic photon can be obtained as ηu = 1 − 43 TT(λ)0 + 13 TT(λ)0 ,
which is similar to Petela’s coefficient formula and is simple and easy to apply. The
calculation of the exergy of blackbody radiation shows that the error of this formula
is negligible at high temperatures, being no more than 10−8 . The exergy-to-energy
ratio of monochromatic radiation in blackbody surface radiation agrees with previous
studies. Using the emission spectrum of engineering radiation, this formula can be
used to conveniently calculate radiative exergy in most engineering fields.
The entropy
v expression
/ of monochromatic radiation photons should be the integral
form sv = 0 h T(λ) dv based on the concept of monochromatic photon equivalent
temperature. It can reflect the distinctions between radiative energy and thermal
energy. The entropy and exergy of monochromatic photons presented in this book
satisfy the thermodynamic relation, indicating that they are essentially harmonious.
References
20. Duysens LNM (1958) The path of light in photosynthesis. Brookhaven Symp Biol 11:18–25
21. Knox RS (1977) Photosynthetic efficiency and exciton transfer and trapping. Primary Process
Photosynthesis 2:55–97
22. Lems S, Van Der Kooi HJ, De Swaan AJ (2010) Exergy analyses of the biochemical processes
of photosynthesis. Int J Exergy 7(3):333–351
23. Juretić D, Županović P (2003) Photosynthetic models with maximum entropy production in
irreversible charge transfer steps. Comput Biol Chem 27(6):541–553
24. Lavergne J, Joliot P (2000) Thermodynamics of the excited states of photosynthesis. Energy
transduction in membranes. WA Cramer, Rockville Biophysical Society
25. Meszéna G, Westerhoff HV, Somsen O (2000) Reply to comment on non-equilibrium
thermodynamics of light absorption. J Phys A: Math Gen 33(6):1301
26. Chen Z, Mo S (2007) Effective temperature and exergy of monochromic blackbody radiation.
Prog Nat Sci 17(10):1250–1254
27. Petela R (2009) Radiation spectra of surface. Int J Exergy 7(1):89–109
28. Modest MF (2003) Radiative heat transfer. Academic Press
29. Planck M (1914) The theory of heat radiation. Philadelphia Blakistons
30. Ito E, Komatsu T, Suzuki H (1998) The entropy generation in visual-pigment system by the
absorption of light. Biophys Chem 74(1):59–70
31. Kirwan AD (2004) Intrinsic photon entropy? The darkside of light. Int J Eng Sci 42(7):725–734
Chapter 3
Gas Radiation Model Under Complex
Combustion Conditions
The actual combustion chambers (e.g., pressurized boilers, gas turbines, internal
combustion engines, and aircraft engines) are all pressurized. Additionally, many
triatomic gases exist in the flue gas during oxy-fuel combustion (CO2 ). Thus, inves-
tigating gas radiation models under complex combustion conditions (e.g., high-
pressure environments and oxy-fuel atmospheres) is critical. The first section of this
chapter summarizes the steps and methods for developing a new efficient and accurate
WSGG model based on the benchmark model. Second, the EM2C laboratory’s SNB
model is used as a benchmark to improve the WSGG model and establish a new set
of correlations for a broader range of parameters for H2 O/CO2 mixtures, adapting
to radiation heat transfer calculations under complex combustion conditions. The
new WSGG model and its parameters are suitable for pressure, temperature, path-
length, and molar ratio ranges of 1–30 bar, 500–2500 K, 0.001–60 m, and 0.125–4
for H2 O/CO2 , respectively. Thus, the new WSGG model is suitable for a wide range
of fuels under pressurized conditions. When compared to the benchmark model, the
one-dimensional case results show that the new WSGG model can accurately predict
the radiation heat transfer of the gas mixture under pressurized conditions. Further-
more, the new model is used to investigate the effect of pressure on the radiation heat
transfer of the gas mixture. The findings show that increasing the pressure within a
certain range can improve the radiation heat transfer of a gas mixture, and this effect
is also related to temperature, molar ratio, and path-length. In conclusion, the new
WSGG model can be used to calculate and study actual combustion conditions as
well as provide a theoretical foundation for the study of radiative energy utilization.
Revised from Shan et al., New pressurized WSGG model and the effect of pressure on the radiation
heat transfer of H2 O/CO2 gas mixtures, International Journal of Heat and Mass Transfer, Copyright
(2018), with permission from Elsevier.
Revised from Shan et al., New weighted-sum-of-gray-gases model for typical pressurized oxy-
fuel conditions, International Journal of Energy Research, Copyright (2017), with permission from
John Wiley and Sons.
The SNB model is a band model that was proposed by Malkmus in 1967 [1]. It is
the model closest to the LBL model among all gas radiation characteristics models
and is frequently used as a reference model for comparison. The SNB model divides
the spectrum into many small wavenumber intervals Δν; it is assumed that each
wavenumber interval Δν contains N spectral lines with a random distribution of the
intensity or location. Thus, the probability of a spectral line appearing in any position
of Δν is the same. For the wavenumber interval Δν, it is assumed that the intensity
of blackbody radiation is constant in the wavenumber interval Δν; thus, Δν should
be sufficiently small. However, Δν cannot be too small to ensure the validity of the
statistical spectral distribution. The wavenumber interval Δν is generally selected to
be from 5 to 50 cm−1 .
The SNB model parameters developed by EM2C Laboratory (France) in 1997
[2] have been demonstrated in many SNB model databases to be closest to those of
the LBL model based on the HITEMP2010 database, and this SNB model has been
recommended as a benchmark [3, 4]. Furthermore, the EM2C laboratory developed
new SNB parameters in 2012 [5]; the CO2 and H2 O parameters of this model are
based on the CSCD-4000 and HITEMP databases, respectively. However, recent
studies showed that when the two EM2C SNB models are used under pressurized
oxy-fuel combustion conditions, the early SNB models are much closer to the LBL
model based on the HITEMP2010 database [6].
The SNB model of the EM2C divides the spectrum into narrow bands of 25 cm−1 ,
and the transmittance of an isothermal homogeneous gas mixture in a band is
expressed as
⎡ ⎛/ ⎞⎤
2γ ⎝ X P Lkδ
τ v = exp⎣− 1+ − 1⎠⎦ (3.1)
δ γ
where the EM2C database provides the parameters k and δ, and the average spectral
half-width γ is calculated by solving the following equations:
⎧ ( )0.5 [ ( )]⎫
P Tref Tref 0.0792 1 − X CO2 − X O2
γ H2 O = 0.462 X H2 O + ×
Pref T T +0.106X CO2 + 0.036X O2
(3.2)
( )
P Tref 0.7 [ ( ) ]
γ CO2 = 0.07X CO2 + 0.058 1 − X CO2 − X H2 O + 0.1X H2 O (3.3)
Pref T
3.1 Non-gray Gas Radiation Models 35
where Pref = 1 bar and T ref = 296 K. An increase in pressure primarily influences
the radiation characteristics by affecting the average spectral half-width γ and the
pressure term P in Eqs. (3.2) and (3.3). Thus, the gas transmittance decreases as the
pressure increases.
The WSGG model was first proposed by Hotel [7], and the relevant radiation heat
transfer calculation was performed using the WSGG model in conjunction with
the regional method. Modest [8] pointed out that the WSGG model considers the
absorption coefficient and is not limited to calculating total emissivity. This model
can also be used to solve the radiation transfer equation and is suitable for the arbitrary
solution method of the radiation transfer equation. The core idea of the WSGG model
is to use the n gray gases to represent the actual non-gray gas; the total emissivity
or the radiative heat flux of each gray gas is multiplied by the weight, and the total
value equals their sum [9]. The emissivity on the path-length L through the WSGG
model can be expressed as
Σ
n
[ ]
ε= ai 1 − exp(−κi X P L (3.4)
i=1
where X is the mole fraction of the gas mixture„ expressed as the sum of the molar
fractions of H2 O and CO2
X = X H2 O + X CO2 (3.5)
Σ
n
ai = 1 (3.6)
i=0
The WSGG model is concerned with calculating the absorption coefficient κ i and
the weighting factor ai of each gray gas.
36 3 Gas Radiation Model Under Complex Combustion Conditions
The radiation transfer equation (RTE) describes the process of radiation transfer in
a medium; for a medium with absorption, emission, and non-scattering properties,
the radiation transfer equation is [9]
d Iν
= κν Ibν − κν Iν (3.7)
ds
where I ν and I bν are the spectral and blackbody spectral emission intensity, respec-
tively. The left side of the RTE represents the change of radiation intensity at a
point along the path-length in the medium. The first term on the right represents the
emission intensity at a point along the path-length in the medium, and the second
represents the radiation intensity absorbed at a point along the path-length in the
medium.. The RTEs for the entire spectrum as well as all spatial directions must be
calculated to obtain the radiation transfer of the entire space.
Numerous approaches were noted to solve the RTE. The discrete coordinate
method (DOM) was used in this book because it had broad applicability [10].
Following that, the DOM for the RTE coupled with the WSGG and SNB models
is represented in one- and two-dimensional settings.
The RTE of each gray gas in the WSGG model was solved for the radiation
transfer problem under the one-dimensional condition and the results of four gray
gases were then summed. For each gray gas, the RTE can be written as:
d Ii (x, ŝ)
= κi ai Ib (x) − κi Ii (x, ŝ) (3.8)
ds
/
where the blackbody emission intensity is expressed as Ib (x) = σ T (x)4 π ; the
items in Eq. (3.8) are determined based on the state at position x. Figure 3.1a depicts
the solution diagram of the one-dimensional RTE for the DOM. The one-dimensional
transfer directions are divided into positive and negative series during the solution
process, the cosine of the angle between each direction and the x-axis is ul , and the
direction weight is wl . Thus, the RTE in each direction can be expressed as follows:
Fig. 3.1 Diagram of DOM calculation process: a one-dimensional case and b two-dimensional
case
The one-dimensional axis was divided into some computing units; therefore, the
partial differential equation was written as a difference equation for each unit, and the
latter format was selected as the differential way. The discrete form of the difference
equation is as follows:
After calculating the intensity of the emission along the path-length, the radiative
heat flux q and the radiation source q̇ can be written as:
ΣΣ
qm = 2π μl wl [Ii+ (xm , l) − Ii− (xm , l)] (3.12a)
i l
ΣΣ
q̇m = {2π κi wl [Ii+ (xm , l) + Ii− (xm , l)] − 4π κi wl Ib,i (xm )} (3.12b)
i l
The radiation heat flux q and the radiation source term q̇ in the calculation of
radiation transfer were mainly investigated. Additionally, the radiation source term
q̇ can be added to the energy equation of the CFD as a part of the source term for the
numerical simulation of a practical problem.
As demonstrated in Sect. 3.1.1, the SNB model can be used to calculate the trans-
mittance of a given path-length. Additionally, the emissivity can be determined using
the relationship between transmittance and emissivity. Thus, the average emissivity
is expressed as
ε L = 1 − exp(−κ X P L) (3.13)
d Iν (x, ŝ)
= κν aν Ib (x) − κν Iν (x, ŝ) (3.14)
ds
The DOM was used to solve the RTE for each band of the SNB model, and the
results were then summed. The radiation heat flux q and the radiation source term q̇
were then calculated, and the detailed processes are similar to those of the WSGG
model.
For the radiation transfer problem under the two-dimensional condition, the
WSGG model was used to describe the radiation characteristics. The radiation
transfer equation for DOM can be written as shown in Eq. (3.15), where μl or ηl is
the cosine of the angle between the discrete direction and the x- or y-axis. Moreover,
Fig. 3.1b depicts the calculation diagram:
The calculation region was divided into a square lattice; thus, Δx = Δy. The
difference format is as follows:
Similarly, the positive and negative radiation transfer equations of ŝ2 are
analogous.
When the SNB model is used under the two-dimensional condition, the DOM
with the SNB model is similar to that with the WSGG model. However, unlike the
WSGG model, which needs only the bands of four gray gases for calculations, the
SNB model needs the radiation transfer equations of all the narrow bands to be solved
and the results to be added together. Thus, the RTE can be written as follows:
Σ
n
ε= ai [1 − exp(−κ p,i X L)] (3.19)
i=1
There are two types of WSGG models for oxy-fuel combustion: those with three
gray gases and those with four gray gases. Johansson’s study showed that the WSGG
model with four gray gases and one white gas is better than that with three gray gases
and one white gas [11]. Currently, the majority of the WSGG models proposed by
scholars for oxy-fuel combustion have four gray gases. Thus, the model containing
four gray gases and one white gas is chosen to investigate the new WSGG model
correlations for pressurized oxy-fuel combustion; thus, n = 4. In the WSGG model,
the weight value ai is mainly dependent on temperature; thus, ai can be expressed as
the polynomial of the dimensionless temperature term with m-th power
Σ
m ( )m− j
T
ai = ci, j (3.20)
j=0
Tref
In Eq. (3.20), T ref is the reference temperature, which can make the polynomial
coefficient ci,j dimensionless, so as to improve the accuracy of the polynomial [11,
12]. The recommended value of T ref is 2000 K in the literature [13] for higher model
accuracy, and this value is also used here. As shown in Eq. (3.21), the polynomial
coefficient ci,j and the absorption coefficient κ i with same molar ratio M and different
pressures are expressed as the polynomial of pressure:
κi = K 2i P 2 + K 1i P + K 0i (3.21b)
3.2 Development of WSGG Model Under Complex Combustion Conditions 41
M = X H2 O / X CO2 (3.22)
To improve the accuracy of the polynomial in Eq. (3.22), piecewise fitting within
different pressure ranges is used.
In this section, the detailed fitting process of WSGG model correlations under
pressurized conditions is elaborated.
(1) First, the SNB model of the EM2C database was used to calculate the emissivity
of a gas mixture over a pressure range of 1–30 bar; then, the emissivity table was
established as a database. Under each pressurized condition, the molar ratio range
was 0.125–4, the temperature range was 500–2500 K, and the path-length range was
0.001–60 m.
When calculating emissivity under various conditions with the SNB model, the
transmittance τ ν in a narrow band was first calculated using Eq. (3.1), and the emis-
sivity in the narrow band of the gas mixture was then calculated using Kirchhoff’s
law and the relationship between the gas transmittance and the absorptivity:
εν = α ν = 1 − τ ν (3.23)
The emission intensity in the narrow band I v was then calculated, followed by
the sum of all emission intensities to obtain the emission intensity I for the entire
spectral range. The emission intensity in the narrow band I v is
Iν = Ibν × ε ν (3.24)
ε = I /Ib (3.25)
Ib = σ T 4 /π (3.26)
certain path-length, all emissivity values with different temperatures were added to
determine the sum of emissivity εL , which is only associated with the path-length L:
Σ
εL = εt,L (3.27)
t
Σ
n
εL = Ai [1 − exp(−κi X L)] (3.28)
i=1
where Ai is the sum of the single gray gas weights ai at different temperatures:
3.2 Development of WSGG Model Under Complex Combustion Conditions 43
Σ
Ai = ai,t (3.29)
t
In order to make the new model suitable to the pressure range of 1–30 bar, piecewise fit
was used to improve polynomial accuracy in different pressure ranges,. Considering
most combustion chambers are in the temperature range of 500–2500 K, the new
model uses this temperature range [14–17] Additionally, the temperature range of
many published WSGG models is included in it [13, 18, 19]. The path-length range
is 0.001–60 m, which covers nearly the entire size range of the combustion chambers
[20–22]. The model is also suitable for a broad molar ratio range of 0.125–4, which
can represent various organic fuels, where M is 0.5 for acetylene and benzene, M is
1 for olefins and many oil fuels [23, 24], and M is 2 for the combustion of methane
or methanol. In addition, the molar ratio M less than 1 can represent the dry cycle
conditions of oxy-fuel combustion, while the M more than or equal to 1 can represent
the wet cycle conditions oxy-fuel combustion [11–13].
44 3 Gas Radiation Model Under Complex Combustion Conditions
1.0 1.0
(a) (b)
0.8 0.8
T=1500 K T= 1500 K
0.6 P=5 bar 0.6 P=5 bar
Emissivity
Emissivity
M=0.15 M=1.5
0.4 0.4
SNB
SNB
WSGG (Specific table)
WSGG (Specific table)
WSGG (Linear interpolation)
WSGG (Linear interpolation)
0.2 0.2
0.0 0.0
0 10 20 30 40 50 60 0 10 20 30 40 50 60
Path-length L (m) Path-length L (m)
Fig. 3.3 Comparison between two methods for applying new WSGG model for a given MR
(1500 K, 5 bar): a MR = 0.15 and b MR = 1.5
3.3 Model Calculation Results 45
research. Although a minor error exists, it is not visible. Thus, the specific table
method in Table 3.1 is recommended first for a given M in actual CFD calculations.
In order to investigate the applicability of the new model under high pressure condi-
tions and the effect of pressure on the radiation heat transfer of the gas mixture,
the one-dimensional infinite parallel plates cases are used in this study because
one-dimensional calculation is not only simple but also enough to investigate the
accuracy of the model and the effect of pressure on the radiation heat transfer. Thus,
many studies used one-dimensional cases [11–13, 18, 19]. In the course of inves-
tigating model accuracy, non-isothermal homogeneous cases and non-isothermal
non-homogeneous cases were chosen for study because non-isothermal conditions
are more consistent with real conditions. Furthermore, the two boundary plates are
blackbody walls with temperatures of 1000 K, the temperature distribution between
the plates can be found in [13, 25, 26], and this temperature distribution is similar to
that of a real furnace
( )
2π x
T = 1400 − 400 cos (K) (3.30)
L
where x is the position of the one-dimensional length coordinate, and L is the distance
between the two plates. For non-isothermal, homogeneous conditions (case 1 series),
five cases with different M of 0.25, 0.5, 1, 2, and 4 were selected. The total molar
fraction was set to 0.9, which covers both air–fuel and oxy-fuel conditions [27]. Table
3.2 displays the specific information.
For the non-homogeneous conditions, similar to literature [18], the molar fractions
of H2 O and CO2 are varying, while the molar ratio remains constant. In this study,
three different cases with the M of 0.5, 1, and 2 were selected; furthermore, the molar
fraction distributions of H2 O in these three cases are expressed as Eq. (3.31), and the
specific information of the cases is shown in Table 3.3. Furthermore, when solving
RTE with DOM, the S 8 discrete format is used.
( )
2π x
X H2 O = 0.25 + 0.05 cos (3.31a)
L
( )
2π x
X H2 O = 0.4 + 0.08 cos (3.31b)
L
( )
2π x
X H2 O = 0.5 + 0.1 cos (3.31c)
L
In this section, the accuracy of new model in the emissivity calculation is investigated
under pressurized oxy-fuel conditions, and the influence of the pressure, temperature,
and molar ratio on the emissivity of the gas mixture is also investigated.
First, the results of the benchmark in this study were compared with the results
published in [27] using the same model (SNB model of EM2C) to verify the accuracy
of the benchmark model used in this study. Figure 3.4 shows the comparison of the
emissivity from the SNB model with the published emissivity results. The results
of the SNB model used in this paper were consistent with previously published
benchmark results [27], indicating that the SNB model of EM2C in this paper is
accurate and can be used as benchmark.
First, the new WSGG model was used to predict the emissivity of the gas mixture,
and the results were compared to the benchmark solution. In addition, Fig. 3.5 depicts
the emissivity of a gas mixture calculated using the WSGG model and the benchmark
SNB model under various conditions.
Figure 3.5a depicts the emissivity of the gas mixture at different pressures of 1
(atmospheric), 2, 5, 15, and 30 bar with a path-length range of 0.001–60 m with a
temperature of 1500 K and molar ratio of 1. It can be seen that the results of the
WSGG model are in good agreement with that of the SNB model, indicating that
the new WSGG model has suitable accuracy when calculating the emissivity at a
3.3 Model Calculation Results 47
0.8
1000 K
0.6
Emissivity
1800 K
0.4
SNB in calculation
Benchmark
0.2
0.0
0 10 20 30 40 50
Path length L (m)
Fig. 3.4 Comparison of the emissivity calculated using the SNB model with the published
emissivity results in literature
1.0 1.0
(a) (b)
0.8 0.8
0.6 0.6
P=5 bar
Emissivity
Emissivity
M=1
SNB 1bar WSGG 1bar
0.4 0.4 SNB 600K WSGG 600K
T=1500 K SNB 2bar WSGG 2bar
SNB 900K WSGG 900K
SNB 5bar WSGG 5bar
M=1 SNB 15bar WSGG 15bar
SNB 1200K WSGG 1200K
SNB 1500K WSGG 1500K
0.2 SNB 30bar WSGG 30bar 0.2 SNB 1800K WSGG 1800K
0.0 0.0
0 10 20 30 40 50 60 0 10 20 30 40 50 60
Path-length L (m) Path-length L (m)
1.0 1.0
(c) (d)
0.8 0.8
Emissivity
P=5 bar
SNB
0.4 0.4 T=1500 K WSGG
SNB M=0.125 WSGG M=0.125
SNB M=0.5 WSGG M=0.5
P=5 bar WSGG-J
WSGG-Y
0.2
SNB M=1 WSGG M=1
0.2
M=1 WSGG-S
SNB M=4 WSGG M=4
0.0 0.0
0 10 20 30 40 50 60 0 10 20 30 40 50 60
Path-length L (m) Path-length L (m)
Fig. 3.5 Emissivity of the gas mixture a at different pressures (1500 K, M = 1), b at different
temperatures (5 bar, M = 1), c with different molar ratios (5 bar, 1500 K), and d calculated by
different WSGG models (5 bar, 1500 K, M = 1)
48 3 Gas Radiation Model Under Complex Combustion Conditions
high pressure. In addition, the emissivity of the gas mixture was greater when the
pressure was raised; furthermore, the emissivity increased quickly when the pressure
was below 5 bar while slowly when the pressure was higher than 5 bar. This result
indicates that elevating the pressure can enhance the emissivity of the gas mixture
and that the effect of pressure on the emissivity is reduced as the pressure increases.
Furthermore, there is a significant effect of pressure on the emissivity when the
pressure is below 5 bar.
The emissivity at different temperatures with a molar ratio of 1 at a pressure of
5 bar is shown in Fig. 3.5b, and it can be found that the results of the WSGG model
were consistent with that of the SNB model. As shown in Fig. 3.5b, the emissivity
of the gas mixture was reduced when the temperature increased from 600 to 1800 K;
furthermore, the decrease of emissivity from 600 to 1200 K was more than that
from 1200 to 1800 K. This indicates that the emissivity of the gas mixture decreases
with the temperature and that the effect of temperature is reduced as the temperature
increases under high-pressure conditions; this is consistent with the results in the
literature [13].
Figure 3.5c shows the emissivity of the gas mixture with different molar ratios at
a temperature of 1500 K and pressure of 5 bar. Through comparing the cases with
different molar ratios of 0.125, 0.5, 1, and 4, it can be seen that the increase of molar
ratio could enhance the emissivity of the gas mixture, moreover, the increment of
emissivity from the M of 0.125 to 0.5 was more than that from the M of 0.5 to 4. This
indicates that there is greater emissivity for higher molar ratio under high-pressure
conditions and that the effect of M is reduced as the M increases; thus, the higher
the ratio of H/C is in the fuel, the greater the emissivity of flue gas. Furthermore, for
pressurized oxy-fuel combustion, the flue gas emissivity of the wet cycle condition
is higher than that of the dry cycle condition, and the molar ratio has a greater effect
on the emissivity of the flue gas in the dry cycle condition than in the wet cycle
condition.
Figure 3.5d further shows the emissivity calculated through different WSGG
models with a molar ratio of 1 at a temperature of 1500 K and a pressure of 5 bar.
The emissivity was calculated through the new pressurized WSGG model (WSGG)
proposed in this paper and three WSGG models proposed by Smith (WSGG-S)
[28], Yin (WSGG-Y) [12] and Johansson (WSGG-J) [11]. Subsequently, the emis-
sivity results were compared with the benchmark from the SNB model. As shown
in Fig. 3.5d, the new WSGG model can accurately predict the emissivity of the
gas mixture in the path-length range of 0.001–60 m under high-pressure conditions.
Furthermore, the results of the WSGG models proposed by other researchers were
quite different from the benchmark, which further indicates that the results are not
ideal when the atmospheric WSGG model is applied in high-pressure conditions;
thus, it is necessary to develop new pressurized WSGG model [6]. In addition,
Johansson’s WSGG model is closer to the benchmark than Yin’s WSGG model;
this is because Yin’s model is based on the EWB model, while Johansson’s model
is based on the SNB model of EM2C, which is the benchmark used in this study.
Smith’s model was developed in an early year and based on an old database of the
wide band model; furthermore, its structure is not optimal, and its path-length range
3.3 Model Calculation Results 49
is limited. Therefore, Smith’s model has been completely unsuitable for the radiation
characteristic calculation of the gas mixture under high-pressure conditions.
In this section, the results of one-dimensional cases were mainly discussed to verify
the adaptability of the new WSGG model under high-pressure conditions. The
non-isothermal, homogeneous cases and non-isothermal, non-homogeneous cases
described in Sect. 3.3.1 were investigated. Furthermore, each case was evaluated
at four different high pressures: 2, 5, and 30 bar; furthermore, three different path-
lengths of 1, 5, and 10 m were used in each case for comparison. In the course
of numerical calculation, one-dimensional space was divided into 100 grids. The
radiation heat transfer results of one-dimensional cases mainly include the radiation
heat flux and the radiation source term, which were calculated by the new WSGG
model and the SNB model, respectively. Subsequently, the errors were analyzed to
determine the accuracy of the new model. In this study, the errors of the radiation
heat flux and radiation source term can be defined as [18]:
|qWSGG − qSNB |
δq = × 100% (3.32)
max |qSNB |
|q̇WSGG − q̇SNB |
δ q̇ = × 100% (3.33)
max |q̇SNB |
where, symbol ‘max’ represents the maximum absolute value of the radiation heat
flux or the radiation source term calculated by the SNB model. The radiation heat
flux error at the wall δqwall and the radiation source error at the midpoint δ q̇half were
primarily considered in this study because the heat flux error δqwall and source term
error δ q̇half are typically the largest and are frequently of concern to researchers [13,
18]. Furthermore, in each case, the average heat flux error δqavg and the average
source term error δ q̇avg were investigated. The errors of the new model in these cases
with varying high pressures and path-lengths are listed in Table 3.4.
Table 3.4 was used to investigate the new model errors at various high pressures.
In the 5-bar pressurized condition, the maximum wall heat flux error was 14.85%
(case 1.5, 10 m), and the average heat flux error was not > 5.9% (case 1.5, 10 m).
Additionally, the maximum midpoint source term error was 3.93% (case 1.5, 1 m),
and the average source term error was no more than 3.26% (case 1.5, 10 m). For non-
homogeneous conditions, the maximum wall heat flux error was 11.28% (case 2.3,
10 m), and the average heat flux error was not > 4.28% (case 2.3, 10 m). Additionally,
the maximum midpoint source term error was 3.18% (case 2.1, 5 m), and the average
source term error was not > 3.33% (case 2.3, 1 m). Thus, the maximum average
error at 5 bar pressure is no more than 6%. Furthermore, the average error was no
more than 5% at the pressure of 2 bar, this indicates that the errors change little
50 3 Gas Radiation Model Under Complex Combustion Conditions
from 2 to 5 bar. However, on the whole, there are more errors for higher pressure.
At a pressure of 30 bar, for homogeneous condition, the maximum of average heat
flux error increment was 7.79% (case 1.5, 10 m), and the maximum of average
source term error increment was 3.66% (case 1.5, 5 m) when compared with the
5 bar pressure conditions; while for non-homogeneous conditions, the maximum of
average heat flux error increment was 6.41% (case 2.3, 10 m), and the maximum of
average source term error increment was 3.23% (case 2.1, 10 m). It can be seen that
there is not significant increment of the average error, and this indicates that the error
3.3 Model Calculation Results 51
of the new model will increase with the pressure but remains within a certain range;
therefore, the new model can ensure the accuracy of calculation.
In addition to the pressure, it can be found in general that the errors would be
slightly more for the cases with a higher M; it implies to some extent that the model
accuracy is slightly worse when the weighted sum of gray gases method is used
to predict the radiation characteristics of the flue gas with a high moisture ratio. In
addition, the errors will also increase with the path-length, especially for the heat
flux. Table 3.4 also shows that the average errors in the cases with a path-length of
5 m were very low (below 7%), this indicates that the new model has a practical
significance since the size of most actual combustion chambers is no more than 5 m
[17, 22]. Moreover, it also can be found that the errors of the non-homogeneous
conditions were less than that of the homogeneous conditions; therefore, the new
model is more suitable for the actual problems with non-homogeneous condition.
Figure 3.6a, b intuitively reflect the new WSGG model results of the homogeneous
cases (case 1.2) with M = 0.5 under different pressures (L = 1 m). Figure 3.6c, d
show the new WSGG model results of the homogeneous cases (case 1.2) with M
= 0.5 for different path-lengths (P = 5 bar). Furthermore, the new WSGG model
results of the homogeneous cases (P = 5 bar, L = 1 m) with different molar ratios
are shown in Fig. 3.6e, f.
The radiation heat flux curves at different pressures were similar as shown in
Fig. 3.6a, b, and increasing pressure can increase the absolute value of the radiation
heat flux. When calculating the heat flux and source term in high-pressure conditions,
e.g., 2, 5, and 30 bar, the new WSGG model agrees well with the benchmark. Further-
more, the error of the midpoint source term is slightly higher in the 30-bar pressure
condition than in the low-pressure conditions. The error results in Table 3.4 show that
the new WSGG model has good accuracy under different pressures, and in general,
the higher the pressure, the greater the source term errors. Figure 3.6c, d show that
the new WSGG model results for different path-lengths of 1, 5, and 10 m are consis-
tent with the benchmark. Furthermore, the wall heat flux error with a path-length of
10 m is slightly greater than that with a path-length of 1 m, indicating that the heat
flux error will be greater in general for longer path-lengths. Furthermore, as shown
in Table 3.4, the new WSGG model errors for different path-lengths are not large,
indicating that the new model is accurate. Figure 3.6e, f show that the new model
results agree well with the benchmark for various molar ratios. Furthermore, Table
3.4 shows that the errors will be slightly higher for higher molar ratios. The results
shown in Fig. 3.6 and Table 3.4 concluded that the new WSGG model is in good
agreement with the reference model for calculating the radiation heat transfer of the
flue gas mixture, and that, while different conditions (e.g., pressure, path-length, and
molar ratio) will affect model accuracy, the errors are still not very high and within a
certain range in general. Thus, the new WSGG model has reasonable accuracy and
can be used to calculate the radiation heat transfer of a gas mixture under pressurized
conditions.
52 3 Gas Radiation Model Under Complex Combustion Conditions
150 600
(a) Case 1.2 (b) Case 1.2
100 M=0.5 M=0.5
L=1 m L=1 m
50 0
0
SNB 2 bar
WSGG 2 bar
SNB 5 bar SNB 2 bar
-50 -600 WSGG 2 bar
WSGG 5 bar
SNB 30 bar SNB 5 bar
WSGG 30 bar WSGG 5 bar
-100 SNB 30 bar
WSGG 30 bar
-150 -1200
0.0 0.5 1.0 0.0 0.5 1.0
x (m) x (m)
150 500
(c) (d) Case 1.2
Case 1.2
100
M=0.5 Radiation source term (kW/m3) M=0.5
Radiation heat flux (kW/m2)
0 SNB 1m
WSGG 1m
SNB 5m SNB 1m
-50 WSGG 5m -500 WSGG 1m
SNB 10m SNB 5m
WSGG 10m WSGG 5m
-100 SNB 10m
WSGG 10m
-150 -1000
0.0 0.5 1.0 0.0 0.5 1.0
Dimensional coordinate x* Dimensionless coordinate x*
150 500
(e) (f) P=5 bar
100 P=5 bar L=1 m
Radiation source term (kW/m3)
Radiation heat flux (kW/m2)
L=1 m
50 0
0
SNB M=0.25
WSGG M=0.25
SNB M=0.25
-50 SNB M=0.5 -500
WSGG M=0.25
WSGG M=0.5
SNB M=0.5
SNB M=1
WSGG M=0.5
-100 WSGG M=1 SNB M=1
WSGG M=1
-150 -1000
0.0 0.5 1.0 0.0 0.5 1.0
x (m) x (m)
Fig. 3.6 New WSGG model results of homogeneous cases: a radiation heat flux and b radiation
source term under different pressures (case 1.2, M = 0.5, L = 1 m); c radiation heat flux and d
radiation source term with different path-lengths (case 1.2, MR = 0.5, P = 5 bar, x* = x/L); e
radiation heat flux and f radiation source term with different molar ratios (M = 0.25 for case 1.1,
M = 0.5 for case 1.2, M = 1 for case 1.3, P = 5 bar, L = 1 m)
3.3 Model Calculation Results 53
This section investigates the effect of pressure on the radiation heat transfer of the
gas mixture. The isothermal homogeneous cases are discussed in order to investigate
the effect of pressure on the radiation heat flux of the gas mixture with various
molar ratios, temperatures T, and path-lengths L. The isothermal, homogeneous case
is convenient for comparative study since it has explicit characteristic parameters;
thus, the isothermal, homogeneous results are of general significance.
The effect of pressure on the average radiation heat flux of a gas mixture with
different molar ratios is shown in Fig. 3.7. In this case, the average radiation heat flux
is defined as the average value of all the radiation heat fluxes at every grid point along
the path-length in a one-dimensional case. In the isothermal, homogeneous cases,
five kinds of gas mixtures have different M of 0.25, 0.5, 1, 2 and 4, respectively; the
total molar fraction is 0.9; the temperature is 1500 K; and the path-length is 1 m. In
the pressure range of 1–30 bar, the effect of pressure on the average radiation heat
flux was investigated.
40 90
(a) T=1500 K (b)
80
Average radiation heat flux (kW/m2)
Average radiation heat flux (kW/m2)
L=1 m M=1
1200 K
1350 K
35 70
L=1 m 1500 K
1650 K
60 1800 K
50
30
M=0.125 40
M=0.25
M=1 30
25 M=2
M=4 20
10
20
0 5 10 15 20 25 30 0 5 10 15 20 25 30
Pressure P (bar) Pressure P (bar)
40
(c)
T=1500 K
Average radiation heat flux (kW/m2)
35 M=1
30
25
1m
3m
5m
20 6m
10 m
15
0 5 10 15 20 25 30
Pressure P (bar)
Fig. 3.7 Effect of pressure on the average radiation heat flux of the gas mixtures a with different
molar ratios, b at different temperatures, and c with different path-lengths
54 3 Gas Radiation Model Under Complex Combustion Conditions
The average radiation heat flux of the gas mixture increased rapidly at first, then
slowly, and finally began to decrease in each case from 1 to 30 bar pressure (Fig. 3.7a).
It implies that increasing the pressure in a specific range can improve the radiation
heat transfer of the gas mixture. Additionally, the enhancing effect will gradually
diminish as the pressure increases. In contrast, the radiation heat transfer is reduced
if the pressure exceeds a certain threshold. Thus, an optimal pressure was noted at
which the average heat flux of a gas mixture is greatest.
The optimum pressure to improve heat transfer was ~13, ~5, and ~3 bar for the
M of 0.25, 1, and 4, respectively (Fig. 3.7a). This means that as the molar ratio of
the gas mixture increases, the pressure at which the radiation heat flux begins to
decrease is lower. Thus, the optimal pressure for enhancing heat transfer decreases.
Increasing the H2 O concentration in the gas mixture results in the best radiation
heat transfer enhancement at lower pressure. In practical industrial applications, the
optimal pressure will be lower for the fuel with a higher H/C ratio. Thus, the optimal
pressure for enhancing the radiation heat transfer of methane combustion will be
lower than that of fuel oil combustion. Additionally, the optimal pressure for the
radiation heat transfer of wet flue gas cycle oxy-fuel combustion is lower than that
of dry flue gas cycle oxy-fuel combustion. Furthermore, Fig. 3.7a shows that if the
molar ratio was lower, the variation range of radiation heat flux over the pressure
range was larger. This implies that pressure has a greater influence on radiation heat
transfer of gas mixtures with a lower M. Thus, pressure has a greater influence on
radiation heat transfer of the CO2 molecule than the H2 O molecule.
All of the curves in Fig. 3.7a roughly intersected at a pressure of about 9 bar,
implying that there is a pressure boundary dividing the entire pressure range into
two sections. The radiation heat flux of the gas mixture generally increases with
the molar ratio when the pressure is below the boundary pressure, and it decreases
with the molar ratio when the pressure is greater than the boundary pressure. Thus,
increasing the molar ratio of the gas mixture can improve its radiation heat transfer
under lower pressure, but decreasing the molar ratio of the gas mixture is required
to improve its radiation heat transfer under higher pressure.
The effect of pressure on the average radiation heat flux of the gas mixture at
different temperatures is shown in Fig. 3.7b. The five gas mixtures in the isothermal,
homogeneous cases have temperatures of 1200, 1350, 1500, 1650, and 1800 K.
Additionally, the cases have a molar ratio of 1 and a path-length of 1 m. In the
pressure range of 1–30 bar, the effect of pressure on the average radiation heat flux
was investigated. The average radiation heat flux of the gas mixture increased first
and then decreased with pressure at different temperatures as shown in Fig. 3.7b.
Optimal pressure for heat transfer enhancement was noted, as discussed previously
in Fig. 3.7a. Moreover, Fig. 3.7b also shows that the optimal pressure was ~2, ~4, and
~8 bar when the temperature was 1200, 1500, and 1800 K, respectively, indicating
that the optimal pressure increases with temperature. Furthermore, Fig. 3.7b shows
that the higher the temperature is, the more obvious change of the radiation heat
flux caused by pressure, which indicates that pressure has a greater effect on the
radiation heat transfer of the gas mixture at higher temperatures. Thus, for high-
temperature combustion chambers, the pressurization method improves radiation
3.4 Summary 55
heat transfer more clearly, and the pressure range to improve radiation heat transfer
is wider. Furthermore, Fig. 3.7b shows that the average radiation heat flux is greater
at higher temperatures, implying that temperature has a significant influence on the
radiation heat transfer of the gas mixture. The radiation intensity is quite high in
combustion chambers with high temperatures and pressures, such as those found in
rocket engines.
The effect of pressure on the average radiation heat flux of gas mixtures with
different path-lengths is shown in Fig. 3.7c. In the isothermal, homogeneous case,
the path-lengths of the five gas mixtures are 1, 3, 5, 6, and 10 m. Additionally, the
molar ratio of cases is 1 and the temperature is 1500 K. In the pressure range of
1–30 bar, the effect of pressure on the average radiation heat flux was investigated.
The variations in average radiation heat flux with pressure for gas mixtures with
different path-lengths were different as shown in Fig. 3.7c. For gas mixtures with
path-lengths no more than 5 m, the average radiation heat flux increased first and
then decreased with pressure. Additionally, when the radiation heat flux began to
decrease, the pressure was about 5, 2, and 1.5 bar for the path-lengths of 1, 3, and
5 m, respectively. When the path-length was more than 5 m, the average radiation
heat flux of the gas mixture decreased without any initial increase, e.g., the path-
lengths of 6 and 10 m, which are shown in Fig. 3.7c. This suggests that the variation
of radiation heat transfer of the gas mixture with pressure can be affected by the
path-length. In general, for a longer path-length, the optimal pressure for enhancing
radiation heat transfer is lower. Nonetheless, the radiation heat flux decreases with
pressure and does not undergo any process of first increasing when the path-length
exceeds a certain value. Thus, the curves in Fig. 3.7c have no peak, and no optimum
pressure was noted for the cases with a longer path-length. When evaluating the
effect of pressure on radiation heat transfer in practical applications, considering the
size of the combustion chambers and the actual path-length in addition to the specific
flue gas characteristics is necessary. Furthermore, Fig. 3.7c shows that the variation
of radiation heat flux in the pressure range of 1–30 bar is not very different for the
gas mixtures with different path-lengths.
The curves in Fig. 3.7c intersected at approximately 3 bar. Thus, a boundary
pressure exists that divides the pressure range into two sections, with different trends
in the radiation heat flux curves in these two sections. The average radiation heat
flux decreased with path-lengths when the pressure was greater than the boundary
pressure as shown in Fig. 3.7c. The size of the high-pressure combustion chamber
should not be too large to ensure adequate radiation heat transfer intensity.
3.4 Summary
In this chapter, a new WSGG model and its parameters suitable for H2 O/CO2 mixtures
at pressures ranging from 1 to 30 bar were established based on the SNB model. In
addition, the new model was used to calculate the emissivity and one-dimensional
radiation heat transfer of the gas mixture. The results were also compared with the
56 3 Gas Radiation Model Under Complex Combustion Conditions
benchmark. Finally, the new WSGG model was used to investigate the effect of
pressure on the radiation heat transfer of a gas mixture. The following conclusions
can be drawn:
(1) The gas mixture emissivity calculated by the new WSGG model is consis-
tent with the benchmark under high-pressure conditions. Additionally, no ideal
results were noted for other atmospheric WSGG models. Thus, the development
of a new pressurized WSGG model is of great significance.
(2) The results of the new WSGG model in one-dimensional cases are in good
agreement with the benchmark. The maximum average heat flux and maximum
average source term errors are around 10% and 6%, respectively. According to
the most actual pressurized combustion chambers with small size, it can be seen
that the average error is no more than 4% for the cases with a path-length of
1 m. Thus, the new WSGG model can provide a reasonably accurate prediction
of radiation heat transfer for the gas mixtures under high-pressure conditions
and can be used for numerical calculation or design of pressurized combustion
chambers.
(3) In a certain pressure range, increasing the pressure can improve H2 O/CO2
mixture radiation heat transfer. Thus, pressurized oxy-fuel combustion can
improve flue gas radiation heat transfer when compared to atmospheric oxy-fuel
combustion. Both the emissivity and radiation heat transfer intensity of the oxy-
fuel gas mixture increased with pressure. Furthermore, as pressure increases,
this effect gradually diminishes; when pressure exceeds this range, the radiation
heat flux decreases. Thus, an optimum pressure exists to improve the radiation
heat transfer of the gas mixture.
(4) The optimal pressure for radiation heat transfer enhancement can be influenced
by the molar ratio of H2 O/CO2 , temperature, and path-length of the gas mixture.
The higher the molar ratio is, the lower the optimal pressure will be; nevertheless,
the optimum pressure will be higher for a higher temperature; moreover, there
is a lower optimum pressure for the longer path-length.
(5) Pressure has a greater effect on the radiation heat transfer of a gas mixture with
a lower molar ratio as well as a greater influence on the radiation heat transfer
characteristics of CO2 than H2 O. Thus, an increase in pressure can effectively
improve the radiation heat transfer of the dry flue gas recycle condition in
oxy-fuel combustion to wet flue gas recycling. Furthermore, pressure has a
greater influence on the radiation heat transfer of a gas mixture with a higher
temperature. However, the effect of pressure on the radiation heat transfer of
gas mixtures with different path-lengths is similar.
References
1. Malkmus W (1967) Random Lorentz band model with exponential-tailed S−1 line-intensity
distribution function. JOSA 57(3):323–329
References 57
2. Soufiani A, Taine J (1997) High temperature gas radiative property parameters of statistical
narrow-band model for H2 O, CO2 and CO, and correlated-K model for H2 O and CO2 . Int J
Heat Mass Transf 40(4):987–991
3. Goutiere V, Liu F, Charette A (2000) An assessment of real-gas modelling in 2D enclosures. J
Quant Spectrosc Radiat Transfer 64(3):299–326
4. Chu H, Liu F, Zhou H (2011) Calculations of gas thermal radiation transfer in one-dimensional
planar enclosure using LBL and SNB models. Int J Heat Mass Transf 54(21–22):4736–4745
5. Rivière P, Soufiani A (2012) Updated band model parameters for H2 O, CO2 , CH4 and CO
radiation at high temperature. Int J Heat Mass Transf 55(13–14):3349–3358
6. Chu H, Gu M, Consalvi J et al (2016) Effects of total pressure on non-grey gas radiation transfer
in oxy-fuel combustion using the LBL, SNB, SNBCK, WSGG, and FSCK methods. J Quant
Spectrosc Radiat Transfer 172:24–35
7. Hottel H, Sarofim A (1967) Radiative transfer. McGraw-Hill, New York
8. Modest M (1991) The weighted-sum-of-gray-gases model for arbitrary solution methods in
radiative transfer
9. Modest M (2013) Radiative heat transfer, 3rd edn. Academic press, Amsterdam
10. Chandrasekhar S (1960) Radiative transfer. Dover Publications, New York
11. Johansson R, Andersson K, Leckner B et al (2010) Models for gaseous radiative heat transfer
applied to oxy-fuel conditions in boilers. Int J Heat Mass Transf 53(1–3):220–230
12. Yin C, Johansen L, Rosendahl L et al (2010) New weighted sum of gray gases model appli-
cable to computational fluid dynamics (CFD) modeling of oxy-fuel combustion: derivation,
validation, and implementation. Energy Fuels 24(12):6275–6282
13. Kangwanpongpan T, França F, Silva R et al (2012) New correlations for the weighted-sum-of-
gray-gases model in oxy-fuel conditions based on HITEMP 2010 database. Int J Heat Mass
Transf 55(25–26):7419–7433
14. Rashedul H, Kalam M, Masjuki H et al (2014) Numerical study on convective heat transfer of
a spark ignition engine fueled with bioethanol. Int Commun Heat Mass Transfer 58:33–39
15. Arghode V, Gupta A (2010) Effect of flow field for colorless distributed combustion (CDC)
for gas turbine combustion. Appl Energy 87(5):1631–1640
16. Zhang J, Dai B, Meng Y et al (2015) Pilot-scale experimental and CFD modeling investigations
of oxy-fuel combustion of Victorian brown coal. Fuel 144:111–120
17. Bahramian A, Maleki M, Medi B (2017) CFD modeling of flame structures in a gas turbine
combustion reactor: velocity, temperature, and species distribution. Int J Chem Reactor Eng
15(4)
18. Dorigon L, Duciak G, Brittes R et al (2013) WSGG correlations based on HITEMP2010 for
computation of thermal radiation in non-isothermal, non-homogeneous H2 O/CO2 mixtures.
Int J Heat Mass Transf 64:863–873
19. Bordbar M, W˛ecel G, Hyppänen T (2014) A line by line based weighted sum of gray gases model
for inhomogeneous CO2 –H2 O mixture in oxy-fired combustion. Combust Flame 161(9):2435–
2445
20. Yin C (2016) Effects of moisture release and radiation properties in pulverized fuel combustion:
a CFD modelling study. Fuel 165:252–259
21. Gicquel L, Staffelbach G, Poinsot T (2012) Large eddy simulations of gaseous flames in gas
turbine combustion chambers. Prog Energy Combust Sci 38(6):782–817
22. Li Y, Han J, Chen X et al (2017) Numerical study of the internal flow field of a dual pulse
solid rocket motor including conjugate heat transfer. Proc Inst Mech Eng Part G: J Aerosp Eng
231(8):1535–1549
23. Westbrook C, Pitz W, Herbinet O et al (2009) A comprehensive detailed chemical kinetic
reaction mechanism for combustion of n-alkane hydrocarbons from n-octane to n-hexadecane.
Combust Flame 156(1):181–199
24. Seshadri K, Frassoldati A, Cuoci A et al (2011) Experimental and kinetic modeling study of
combustion of JP-8, its surrogates and components in laminar premixed flows. Combust Theor
Model 15(4):569–583
58 3 Gas Radiation Model Under Complex Combustion Conditions
This chapter combines the radiation transfer process with the thermodynamic second
law to obtain more accurate results for energy quality and its variability in the spec-
tral radiation transfer process. First, the fundamental concepts of the monochromatic
radiative exergy theory are reviewed and discussed. This theory is then combined
with the radiation transfer equation, and the spectral radiative entropy and the radia-
tive exergy transfer equations are established and verified using the second law of
thermodynamics. Finally, one-dimensional furnace case calculations are performed
to determine the spectral distribution and variation of radiative energy and radia-
tive exergy in the combustion medium. Furthermore, the effects of the temperature
field, gas molar ratio, and particle number density on spectral radiation characteris-
tics are investigated. This is useful for studies into flame energy spectrum-splitting
conversion systems.
the concept of equivalent temperature and the differences between the radiative and
thermal energy to derive the monochromatic photon exergy as shown in Chap. 2.
The monochromatic photon exergy derived from the infinite-staged Carnot heat
engine can be expressed as an integral form:
(ν ( )
T0
Eν = h 1 − dν (4.1)
Tν
ν0
f T3
Tλ4 = (4.2)
λ
In addition, the corresponding spectral radiative exergy–energy coefficient can be
calculated by comparing it with the photon energy:
( )
4 T0 1 T0 4
ηλ = 1 − + (4.3)
3 Tλ 3 Tλ
when T > T 0 : 4.5638 × 10−3 T 3 = λTλ4 ; when T < T 0 : 5.33 × 10−3 T 3 = λTλ4 .
Furthermore, the entropy of a monochromatic photon was studied using the
infinite-staged Carnot heat engine model, and the entropy of the photon with
frequency ν is:
(ν
hdν
Sν = (4.4)
Tν
0
After substituting Eq. (4.2) into Eq. (4.4), the coefficient of the spectral radiative
entropy to the energy can be obtained using a comparison with the photon energy:
Sν 4
ξν = = (4.5)
hν 3Tλ
Following that, some remarks on the radiative exergy theory are made. First,
the Gouy–Stodola theorem was used to demonstrate that the relationship between
monochromatic photon exergy and entropy represented by Eqs. (4.1) and (4.4) was
proved to be in agreement with:
E = H − H0 − T0 (S − S0 ) (4.6)
4.1 Thermodynamic Theory of Spectral Radiation Transfer in Combustion … 61
Thus, the proposed entropy and exergy theories are intrinsically harmonious.
Furthermore, the new spectral radiative exergy theory focuses on monochromatic
radiation’s ability to do work, which is greater for shorter wavelength radiation.
The radiative exergy in the entire blackbody spectrum corresponds to the blackbody
radiative exergy expression proposed by Petela [1].
Second, it is evident that the equivalent temperature [2, 3] is explicitly defined to
characterize the spectral radiative energy quality; this is different from the concept
of radiation temperature, which is only related to the spectral radiation intensity.
Furthermore, the spectral radiative exergy coefficient derived from the infinite-staged
Carnot heat engine model differs from the Carnot efficiency, reflecting the differences
between the radiation and thermal energy. Thus, the method proposed in Eq. (4.3)
appears to be more appropriate for the theoretical characterization of the spectral
radiative exergy.
Finally, a special relationship between the monochromatic photons and spectral
radiation is investigated from the perspective of radiative exergy. According to the
infinite-staged Carnot engine model, for a photon of hν, hν 0 is the reference state for
calculating its radiative exergy. Furthermore, for blackbody spectral radiation I b,λ ,
I b,λ (T 0 ) is the reference state for calculating its radiation exergy [2, 3]. The macro-
and microscales should have a corresponding relationship. Consequently, Eqs. (4.1)
and (4.2) indicates the following relationship between photons and spectral radiation:
The following is a spectral radiation transfer equation with emission, absorption, and
scattering [4]:
(
d Iλ (s) ( ) κs,λ ( ) ( )
= − κa,λ + κs,λ Iλ (s) + κa,λ Ib,λ + Iλ s' Φ s' , s dΩ' (4.8)
ds 4π
4π
where κ a,λ is the spectral absorption coefficient, κ s,λ is the scattering coefficient, and
Φ(s' , s) is the scattering phase function.
The spectral radiative exergy transfer equation can be directly obtained using the
spectral radiative exergy theory:
62 4 Thermodynamic Calculation of Radiative Energy in Combustion Medium
( )
d E λ (s) ( ) 4 T0 T04
= − κa,λ + κs,λ E λ (s) + κa,λ Ib,λ 1 − +
ds 3 Tλ 3Tλ4
(
κs,λ ( ) ( )
+ E λ s' Φ s' , s dΩ' (4.9)
4π
4π
The first term on the right side of the equation represents the exergy intensity
reduction caused by absorption and scattering, the second term is the exergy of the
emission, and the third term is the exergy enhancement caused by spatial scattering
in the direction of radiation transfer. The spatial propagation of spectral radiative
exergy is described in Eq. (4.9). The transfer process for the radiative exergy is very
similar to the transfer process for radiative energy shown in Fig. 4.1. The exergy is
generated for radiation transfer in a participating medium, and it varies primarily in
the radiation field, as well as in the media [3]. Thus, the exergy loss in the participating
media radiation transfer consists primarily of the exergy loss in the radiation field
and the exergy loss in the media as a result of radiation.
To determine the variation of the radiative exergy in the radiation field, Eq. (4.9)
is integrated over the entire solid angle of the space to obtain the local net exergy
loss of the radiation with a spectrum interval dλ after passing through the differential
unit dV.
( [ ( )]
4 T0 T04
d E R,λ = −κa,λ d V dλ
V
E λ (s) − Ib,λ 1 − + dΩ
3 Tλ 3Tλ4
4π
⎡ ⎤
( (
1 ( ') ( ' ) '
− κs,λ d V dλ ⎣ E λ (s) − E λ s Φ s , s dΩ ⎦dΩ (4.10)
4π
4π 4π
Due to the radiation transfer in the participating media, the media also obtain or
lose energy through radiation absorption or emission; thus, the exergy of the media
varies [3]. As shown in Fig. 5.1b, the local net exergy increment of the media due to
the radiation heat flux can be calculated according to the definition of thermal exergy
[3]:
4.1 Thermodynamic Theory of Spectral Radiation Transfer in Combustion … 63
( ) ( )(
T0 T0 [ ]
V
d E M,λ = d Qλ 1 − = κa,λ d V dλ 1 − Iλ (s) − Ib,λ dΩ (4.11)
TM TM
4π
Thus, Eqs. (4.10) and (4.11) are combined. After the radiation of wavelength λ
passes through a differential unit dV within the spectral interval dλ, the radiation
exergy loss in the entire system is:
d E RV M,λ = −d E R,λ
V
− d E M,λ
V
⎡ ( ) ⎤
T0
( ⎢ E λ (s) − Iλ (s) 1 − ⎥
⎢ TM ⎥
= κa,λ d V dλ ⎢ ( 4 )⎥dΩ
⎣ T0 4 T0 T0 ⎦
4π −Ib,λ − +
TM 3 Tλ 3Tλ4
⎡ ⎤
( (
1 ( ') ( ' ) '
+ κs,λ d V dλ ⎣ E λ (s) − E λ s Φ s , s dΩ ⎦dΩ (4.12)
4π
4π 4π
T4
Because Ib,λ 3T0 4 = 13 Ib,λ (T0 ) according to Eq. (4.7), which is constant for a
λ
radiation with wavelength λ in a spectrum of dλ. That is, its value does not change
after passing through a differential unit dV; therefore:
( (
T4 1
d V dλ Ib,λ 0 4 dΩ = d V dλ Ib,λ (T0 )dΩ = 0 (4.13)
3Tλ 3
4π 4π
d E RV M,λ = −d E R,λ
V
− d E M,λ
V
( [ ( ) ( )]
T0 T0 4 T0
= κa,λ d V dλ E λ (s) − Iλ (s) 1 − − Ib,λ − dΩ
TM TM 3 Tλ
4π
⎡ ⎤
( (
1 ( ) ( )
+ κs,λ d V dλ ⎣ E λ (s) − E λ s' Φ s' , s dΩ' ⎦dΩ (4.14)
4π
4π 4π
The first term on the right side of the equation is the radiative exergy loss caused
by the absorption and emission process and the second term is the exergy loss caused
by the scattering process.
64 4 Thermodynamic Calculation of Radiative Energy in Combustion Medium
The first term on the right side of the equation represents the entropy intensity
reduction due to absorption and scattering, the second term is the entropy of the
emission, and the third term is the entropy enhancement in the direction of radia-
tion transfer from the spatial scattering. Equation (4.15) reflects the spatial transfer
process of radiative entropy. It can be seen in Fig. 4.1 that the transfer process for the
radiative entropy is very similar to the radiative energy and exergy transfer processes.
Therefore, the entropy generation in the participating media radiation transfer also
includes the entropy generation in the radiation field and the entropy generation in
the medium caused by the radiation [5, 6].
For the radiative entropy generation in the radiation field, the local net entropy
increment of the radiation with the spectrum interval dλ after passing through the
differential unit dV can be obtained by integrating Eq. (4.15) over the entire solid
angle of the space:
( [ ]
4Ib,λ
V
d S R,λ = −κa,λ d V dλ Sλ (s) − dΩ
3Tλ
4π
⎡ ⎤
( (
( ') ( ' ) '
− κs,λ d V dλ ⎣ Sλ (s) − 1 Sλ s Φ s , s dΩ ⎦dΩ (4.16)
4π
4π 4π
Due to the radiation transfer in the participating media, the media also obtain or
lose energy through radiation absorption or emission; thus, the entropy of the media
varies [5, 7]. As shown in Fig. 4.1c, the local net entropy increment of the media
due to the radiation heat flux can be calculated according to the definition of thermal
entropy [5, 7]:
(
d Qλ Iλ (s) − Ib,λ
V
d S M,λ = = κa,λ d V dλ dΩ (4.17)
TM TM
4π
As a result, Eqs. (4.16) and (4.17) are combined. After the radiation of wavelength
λ passes through a differential unit dV within the spectral interval dλ, the radiative
entropy generation in the entire system is:
d S RV M,λ = d S R,λ
V
+ d S M,λ
V
4.1 Thermodynamic Theory of Spectral Radiation Transfer in Combustion … 65
( [ ]
4Ib,λ Iλ (s) Ib,λ
= −κa,λ d V dλ Sλ (s) − − + dΩ
3T λ TM TM
4π
⎡ ⎤
( (
1 ( ) ( )
− κs,λ d V dλ ⎣ Sλ (s) − Sλ s' Φ s' , s dΩ' ⎦dΩ (4.18)
4π
4π 4π
The first term on the right side of the equation is the radiative entropy generation
caused by the absorption and emission processes and the second term is the entropy
generation caused by the scattering process.
The spectral radiative exergy transfer equation (Eq. 4.9), the entropy transfer
equation (Eq. 4.15), and the exergy loss or entropy generation in the radiation
transfer process can be combined with various radiation transfer numerical calcula-
tion methods, such as the discrete ordinate method (DOM), finite volume method
(FVM), discrete transfer method, and spherical harmonics method [4]. In this study,
some examples are analyzed using the widely applicable DOM.
According to the Gouy–Stodola theorem in Eq. (4.6), the relationship between the
radiative entropy and exergy could be verified.
First, for the spectral radiative entropy transfer equation and exergy transfer
equation, according to Gouy–Stodola theorem, there is:
[ ]
E λ = Iλ − Ib,λ (T0 ) − T0 Sλ − Sb,λ (T0 ) (4.19)
and:
[ ] [ ]
d Eλ d Iλ − Ib,λ (T0 ) d Sλ − Sb,λ (T0 )
= − T0 (4.20)
ds ds ds
After substituting the radiation transfer equation (Eq. 4.8) and the spectral radia-
tive entropy transfer equation (Eq. 4.15) into the right side of the Eq. (4.20), the
following is obtained:
[ ] [ ]
d Iλ − Ib,λ (T0 ) d Sλ − Sb,λ (T0 )
− T0
( ds ){ ds [ ]}
= − κa,λ + κs,λ Iλ (s) − Ib,λ (T0 ) − T0 Sλ (s) − Sb,λ (T0 )
( [ ])
4Ib,λ 4Ib,λ (T0 )
+ κa,λ Ib,λ − Ib,λ (T0 ) − T0 −
3Tλ 3T0
(
κs,λ { ( ' ) [ ( ') ]} ( )
+ Iλ s − Ib,λ (T0 ) − T0 Sλ s − Sb,λ (T0 ) Φ s' , s dΩ' (4.21)
4π
4π
66 4 Thermodynamic Calculation of Radiative Energy in Combustion Medium
d E RV M,λ = d I R,λ
V
+ d Q VM,λ − T0 d S RV M,λ (4.25)
Using Eq. (4.8) for the radiative intensity and Eq. (4.18) for entropy and
considering Eq. (4.7), the following is obtained:
V
d I R,λ + d Q VM,λ − T0 d S RV M,λ
(
[ ]
= −κa,λ d V dλ Iλ (s) − Ib,λ − Iλ (s) + Ib,λ dΩ
4π
⎡ ⎤
( (
1 ( ') ( ' ) '
− κs,λ d V dλ ⎣ Iλ (s) − Iλ s Φ s , s dΩ ⎦dΩ
4π
4π 4π
( [ ]
4T0 T0 T0
− κa,λ d V dλ T0 Sλ (s) − Ib,λ − Iλ (s) + Ib,λ dΩ
3Tλ TM TM
4π
⎡ ⎤
( (
1 ( ') ( ' ) '
− κs,λ d V dλ ⎣T0 Sλ (s) − T0 Sλ s Φ s , s dΩ ⎦dΩ
4π
4π 4π
4.1 Thermodynamic Theory of Spectral Radiation Transfer in Combustion … 67
( [ ( ) ( )]
T0 T0 4 T0
= κa,λ d V dλ E λ (s) − Iλ (s) 1 − − Ib,λ − dΩ
TM TM 3 Tλ
4π
⎡ ⎤
( (
1 ( ) ( )
+ κs,λ d V dλ ⎣ E λ (s) − E λ s' Φ s' , s dΩ' ⎦dΩ = d E RV M,λ (4.26)
4π
4π 4π
It can be seen that the entropy generation and exergy loss in the radiation process
agree with the laws of thermodynamics, demonstrating that the entropy generation
and the exergy loss are uniform and in internal harmony in the radiation transfer
process.
The numerical verification method proposed by Liu and Chu [7] was used to validate
the new radiative entropy transfer equation. A one-dimensional case is used with
equally high temperatures of the two boundary walls; the wall temperatures are set
to 3500 K, 2500 K, and 1750 K, respectively. The temperature of the absorbing and
emission media is constant at 350 K, the absorption coefficient is set to 0.5, and the
length is 1 m. The DOM is used for the calculation and the discrete format is S 8 .
This method has been applied in the calculation of the one-dimensional radiative
heat flux.
By integrating Eq. (4.18) on a one-dimensional space and the whole spectrum,
we can obtain the entropy generation of the radiation transfer in a one-dimensional
space [7]:
( (
ΓG = d S RV M,λ dλd V (4.27)
V λ
Γ G TL
ΠG = (4.28)
Q
0.8
3500 K
0.7 2500 K
1750 K
0.6
0.0 0.2 0.4 0.6 0.8 1.0
Coordinate x (m)
The results of the theoretical calculation (Eq. 4.29) indicate that the dimensionless
entropy generation in the one-dimensional space is 0.9, 0.86, and 0.8 for the wall
temperatures of 3500 K, 2500 K, and 1750 K, respectively. A comparison between
theoretical results and numerical results in Fig. 4.2 shows that the numerical entropy
production in the whole space is very similar to the theoretical results. According to
Ref. [7], this approach verifies the new spectral radiative exergy theory and numerical
calculation method of the second law calculation for the spectral radiation transfer.
Here, the one-dimensional infinite parallel plate cases are used for calculation, which
represents the pulverized-coal furnace. Therefore, the one-dimensional radiation
media, consists of gas and particle mixture, has the absorption, emission, and scat-
tering properties. The characteristics of radiative energy and radiative exergy are
investigated under the conditions with different flue gas temperature, molar ratio,
and particle characteristics. Three furnace temperature distribution characteristics
are shown in Eq. (4.30), and the corresponding curves are shown in Fig. 4.3.
) (
2π x
T = 1400 − 400 cos (K) (4.30a)
L
( )
2π x
T = 1200 − 300 cos (K) (4.30b)
L
4.2 Radiation Energy Characteristics in Combustion Medium of 1-D … 69
2000
Fig. 4.3 Temperature
1400 K
profiles of the 1800 1200 K
one-dimensional cases 1000 K
1600
Temperature (K)
1400
1200
1000
800
600
0.0 0.2 0.4 0.6 0.8 1.0
x/L
( )
2π x
T = 1000 − 200 cos (K) (4.30c)
L
Four oxy-fuel atmospheres with different molar ratio are investigated; they are all
non-isothermal and homogeneous medium and the molar ratio is shown in Table 4.1.
The absorption and scattering factors for particles need to be calculated from the
Mie theory according to the fuel properties. Generally, the average absorption and
scattering factors are used, that is, the spectral distribution curve is drawn according
to the Mie theory, and then the Planck average is taken. In this book, we do not
discuss and perform the specific calculation of Mie theory. Based on the research
published results, appropriate radiation characteristic parameters are selected for
sensitivity analysis. The particle absorption factor Qabs,coal of coal is taken as 0.9, the
scattering factor Qsca.coal is taken as 1.1, the absorption factor Qabs,ash of ash is taken
as 1, and the scattering factor Qsca,ash is taken as 1.5 [8]. The projected area per unit
volume of coal, Ap,coal , is 0.75 m−1 , and the projected area per unit volume of ash,
Ap,ash , is 0.125 m−1 [9]. The calculation methods of the absorption coefficient and
the projected area per unit volume are:
κ = Q abs A p (4.31a)
70 4 Thermodynamic Calculation of Radiative Energy in Combustion Medium
σ = Q sca A p (4.31b)
A p = (π/4)nd 2 (4.32)
where n denotes the particle number density, which remains constant throughout the
reaction process:
m par t m gas
n= / (4.33)
ρ par t V par t ρgas
( (∞
( )
e= E λ+ (s) − E λ− (s) dλdΩ (4.35)
4π 0
The positive (+) and negative (−) symbols indicate the two directions of transmis-
sion, i.e., right and left, respectively; this is similar to the calculation of the radiative
heat flux.
According to the blackbody radiation distribution, the majority of the radiation
in the high-temperature furnace occurs in the near-infrared region, and the radiation
intensity is weak for the wavelength above 50 μm. Thus, the radiation intensity in
the 0–50 μm wavelength range is used to represent the total radiation intensity in the
calculation. The wavelength interval is 25 nm. The spectral radiative exergy distri-
bution in the 0–50 μm wavelength range is investigated according to the distribution
of the spectral radiation intensity of fuel combustion [10, 11].
4.2 Radiation Energy Characteristics in Combustion Medium of 1-D … 71
Absorption, emission, and particle scattering states exist simultaneously in gas and
particle mixed media. The WSGG model, developed in Chap. 3, and the SNB model
are used to calculate the gas radiation properties to obtain the properties of one-
dimensional spectral radiation heat flux. The H2 O/CO2 gas molar ratio M is 0.25 in
the calculation, and the total proportion of radiation gas is 0.9; the temperature curve is
Eq. (4.30a) with an average temperature of 1400 K, the pressure is 1 bar; additionally,
the two sets of absorption coefficient and scattering coefficient of the particles are 0.63
and 0.77 and 1.35 and 1.65, respectively. The SNB model is a spectral band model,
which is more accurate in calculating spectral radiative energy, whereas the WSGG
model is a global model, which was developed in principle to calculate total heat flux.
However, Fig. 4.4 shows that the spectral characteristics of heat flux calculated by the
WSGG model are very similar to those calculated by the SNB model. This is because
particles with absorption and scattering properties play a significant role in radiation
[12]. The radiation properties of the particles are similar to the gray-body radiation,
and the calculation uses spectrally averaged absorption and emission coefficients.
Thus, under such conditions, the non-gray radiation characteristics of the gas are
weakened, and the calculation results of the WSGG model are very similar to those
of the SNB model. It also shows that the WSGG model can be used for gas radiation
calculation when calculating the spectral radiative properties of combustion media
containing particles. Given the computational efficiency of the WSGG model, the
results in Sect. 4.2 are all computed using the WSGG model introduced in Chap. 3.
45 45
(a) P=1 bar; M=0.25 SNB
(b)
40 40 P=1 bar; M=0.25 SNB
κs=0.63; σ=0.77 WSGG
κs=1.35; σ=1.65 WSGG
35 35
Energy proportion (%)
30 30
25 25
20 20
15 15
10 10
5 5
0 0
0-2 2-3 3-4 4-6 6-10 10-50 0-2 2-3 3-4 4-6 6-10 10-50
Waveband (μm) Wavelength (μm)
Fig. 4.4 Comparison of spectral radiative energy calculated by WSGG model and SNB model
72 4 Thermodynamic Calculation of Radiative Energy in Combustion Medium
45 45
40
(a) Energy 1400 K 40
(b) Exergy 1400 K
P=1 bar; M=0.25 1200 K P=1 bar; M=0.25 1200 K
35 1000 K 35 1000 K
Energy Proportion (%)
25 25
20 20
15 15
10 10
5 5
0 0
0-1 1-2 2-3 3-4 4-6 6-10 10-50 0-1 1-2 2-3 3-4 4-6 6-10 10-50
Waveband (μm) Waveband (μm)
Fig. 4.5 Characteristics of spectral proportion of participating media under different temperature
conditions: a average radiative energy flux and b average radiative exergy flux
consistent with Petela’s black-body radiation exergy efficiency Eq. (2.2), so the radi-
ation energy is more valuable. This is mainly related to the increase in the proportion
of high-frequency radiation energy. Thus, increasing the temperature of the combus-
tion media is an important technical measure in terms of photovoltaic conversion and
radiation energy utilization.
Figure 4.7 further depicts the one-dimensional total radiative energy flux and
radiative exergy flux along the path-length. The radiative energy flux and the radia-
tive exergy have a similar trend, which is similar to the previous analysis; i.e., the
calculation method of the radiative exergy in the media is similar to that of radiative
energy, and the radiative exergy at each point is calculated by the radiative energy.
Furthermore, Fig. 4.7 still directly indicates that the radiative energy and radiation
exergy increase significantly as the average temperature rises, with the maximum
energy flux increasing from 25 kW/m2 at 1000 K to 150 kW/m2 at 1400 K, and the
maximum exergy flux increasing from about 20 kW/m2 to about 120 kW/m2 .
different temperature 75
conditions
70
65
60
55
50
45
0.0 0.2 0.4 0.6 0.8 1.0
x coordinate (m)
74 4 Thermodynamic Calculation of Radiative Energy in Combustion Medium
200 150
(a) Energy (b)
150 Exergy
M=0.25 100 M=0.25
100
50
50
0 0
-50
1400 K -50 1400 K
1200 K
-100 1200 K
1000 K
-100 1000 K
-150
-200 -150
0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0
x coordinate (m) x coordinate (m)
Fig. 4.7 a Radiative energy flux and b radiation exergy flux along the path length under different
temperature conditions
The effect of the gas molar ratio on the radiation properties of participating media
is discussed. The average temperature of the one-dimensional radiation media is
1400 K, and the temperature distribution follows Eq. (4.30a), the gas pressure is
1 bar, and the length of the one-dimensional case is 1 m. The particle radiation char-
acteristics are the same as those discussed above, and the H2 O/CO2 molar ratio M
of the media is chosen from the four cases in Table 4.1. Since the spectral character-
istics of radiation energy and radiation exergy are almost the same, Fig. 4.8a shows
the spectral energy proportion of the radiative media under different molar ratios.
The spectral energy proportion is very similar in all four cases. As the molar ratio
increases from 0.25 to 2, the proportion of radiation in the waveband below 3 μm
decreases slightly, while that in the waveband above 3 μm increases slightly, but
the variation does not exceed 1% point, that is, the molar ratio has little effect on
the radiation spectral characteristics. In fact, this is because, the proportion of H2 O
increases as the molar ratio, and the radiation of H2 O has its special spectrum in
the waveband of 5.6–7.6 μm, so the spectral proportion increases in the waveband
above 3 μm. However, in the presence of particle media with considerable absorption
and scattering coefficients, the particle radiation plays a decisive role in the overall
radiation characteristics, and the effect of individual gas radiation is less. Thus, gas
molar ratio variation has little effect on the radiation spectrum characteristics.
Figure 4.8b depicts the effect of the molar ratio on the average radiative energy
flux and radiative exergy flux. When it increases from 0.25 to 2, the average radia-
tive energy flux increases from 112 to 118.7 kW/m2 . This is because the increase
in the proportion of H2 O helps to improve the radiative heat transfer. Thus, the
corresponding average radiative exergy also increases from 89.2 to 94.4 kW/m2 . In
addition, the total exergy-to-energy ratio shows a downward trend in the range of
79.7–79.5%, which can be regarded as almost no variation. The effect of molar ratio
on exergy-to-energy ratio can be ignored in the presence of particle radiation.
4.2 Radiation Energy Characteristics in Combustion Medium of 1-D … 75
45 125 79.75
(a) Energy (b)
40 M=0.25
P=1 bar; T=1400 K
M=0.5 115
35 79.70
M=1
Energy flux
M=2
10 79.55
75
5
0 65 79.50
0-1 1-2 2-3 3-4 4-6 6-10 10-50 0.0 0.5 1.0 1.5 2.0 2.5
Wavelength (μm) Molar ratio
Fig. 4.8 a Spectral proportion characteristics of media radiation under different molar ratios; b the
influence of molar ratio on average radiative energy flux and radiative exergy flux
In order to consider the influence of particle characteristics, the coal particle diam-
eter is selected as 100 μm and the ash particle diameter is 30 μm. In the previous
calculation, Ap,coal = 0.75 and Ap,ash = 0.125 were selected following Eq. (4.33); the
ratio of particle number density of coal and ash is ncoal /nash = 0.54, and the propor-
tion of coal particles number is 35%. Assuming that the total number of coal and
ash particles remains constant, the radiation characteristics of participating media
are investigated when the coal particle number ratio r varies from 20 to 60%. The
burnout area in the combustion chamber is represented by a smaller proportion of
coal particles, while the combustion area is represented by a larger proportion, and
the larger the proportion, the closer it is to the combustion center.
Figure 4.9a depicts the spectral radiation characteristics of participating media
as the coal particle concentrations vary. The spectral radiation characteristics are
still consistent. When the proportion of coal particle increases from 20 to 60%,
the proportion of spectral radiation in each waveband changes less. This is because
although there are variations in the absorption and scattering coefficients of coal
and ash particles after the number density of coal particles changes, the average
absorption and scattering coefficients are used in the calculation, while the spectral
coefficients are not used. The radiation characteristics are less affected by the average
radiation parameters. However, it can still conclude that when the spectral radiation
parameters of the particles are relatively stable, the particle number density has little
effect on the spectral radiation characteristics of the participating media.
The effect of coal particle concentration on the average radiative energy flux and
radiative exergy flux is further reflected in Fig. 4.9b. The average radiative energy
flux decreases from 99.6 to 121.2 kW/m2 as the coal particle concentration increases
from 20 to 60%, indicating that the radiative heat transfer is greater when the coal
particle concentration is high, i.e., the closer to the combustion flame center, the
greater the radiative heat transfer. The corresponding average radiative exergy also
increases from 79.5 to 96.6 kW/m2 . Furthermore, the radiative energy-to-energy ratio
76 4 Thermodynamic Calculation of Radiative Energy in Combustion Medium
45 130 79.80
(a) Energy (b)
40 T=1400 K 20%
120
M=0.25 30%
35 40%
50% 110
20 90
15 79.70
80
10 Energy flux
Exergy flux
70
5 Exergy-to-energy ratio
0 60 79.65
0-1 1-2 2-3 3-4 4-6 6-10 10-50 20 30 40 50 60
Waveband (μm) Ratio of coal partical concentration (%)
Fig. 4.9 a Spectral proportion characteristics of media radiation under different coal particle
concentrations; b influence of coal particle concentration on average radiative energy flux and
radiative exergy flux
flow is almost stable at about 79.7%. These results indicate that under the combustion
conditions of solid fuels such as pulverized coal, the radiation heat transfer is greater
if the participating media is closer to the center of the combustion area, while the
proportion of ash particles is greater if the media is farther from the combustion center,
and the radiation intensity is smaller. If the combustion temperature is similar, the
change of various parameters has little effect on the radiative exergy-to-energy ratio
of the participating media. Elevating the temperature is the preferred solution to the
radiative exergy in combustion media.
4.3 Summary
This chapter combines the radiation transfer equation with the theory of the
monochromatic photon exergy developed in Chap. 2 by considering the second
law of thermodynamics. Based on the second law, the spectral radiative entropy
transfer equation and radiative exergy transfer equation are established; this process
can be used to numerically analyze the variability of the energy quality during the
spectral radiation transfer process more appropriately and accurately. The radiative
entropy and the radiative exergy are found to be consistent with the thermodynamic
Gouy–Stodola law.
Numerical calculations of the spectral radiative exergy transfer equation can be
used to determine the distribution and variation of the radiation spectral exergy flux
during the radiation transfer process. Calculating the spectral radiative exergy distri-
bution, which is very important for the analysis of the flame energy spectrum-splitting
conversion system in high-temperature engineering, is also possible. The results of
the one-dimensional cases simulating the furnace temperature field show that the
proportion characteristics of spectral radiative energy and spectral radiative exergy
are consistent, implying that the proportion characteristics of spectral energy can be
References 77
used to predict those of spectral radiative exergy. Temperature is the primary factor
influencing spectral radiative energy. A higher temperature indicates that the propor-
tion of radiative energy in the waveband of < 2 μm is greater, and more available
radiative energy exists. At the same temperature, the combustion atmosphere and
coal particle concentrations have little effect on the proportional characteristics of
spectral radiative energy. These parameters will influence the radiative energy charac-
teristics in the actual combustion process by influencing the combustion temperature.
Thus, additional experiments are required to investigate the characteristics of spectral
radiative energy.
References
Revised from Shan et al., Spectral energy characteristics of radiation in oxy-coal combustion for
energy utilization, Fuel, Copyright (2021), with permission from Elsevier.
Revised from Shan et al., Spectral radiation characteristics in semi-coke jet flame for energy
utilization, Fuel, Copyright (2021), with permission from Elsevier.
A tube furnace platform is constructed for the solid fuel combustion experiment.
Figure 5.1 shows the experimental system, which is similar to that used in a previous
study [1]. The tube furnace was 60 cm long. A high-temperature-resistant quartz tube
(140 cm long and 6.5 cm in diameter) was used, with one side of the tube opened
for measurement. Furthermore, the sample was placed at the center of the quartz
tube, and a radiation sensor was placed 97 cm from the sample. In the experiment, a
CT-20 high-temperature radiation sensor (Captec Co., France) was used to measure
the radiative heat flux at the outlet of the quartz tube. The sensor can measure the
total radiation power in the 0.3–50 μm spectrum band and has a response time of
0.25 s. The experimental data were recorded using a Daq-5300 eight-channel data
logger (Fourier System Co., Israel), with an acquisition interval of 1 s. An infrared
filter (Spectrogon Co., Sweden), a standard shortwave pass filter that allows for the
passage of short-wavelength radiation, was placed in front of the radiation sensor
to separate the radiation spectrum. The specific cutoff band λc and transmittance
τ of the different filters are λc = 1125 ± 25 nm, τ = 94.98% (filter no. 1); λc =
4100 ± 100 nm, τ = 89.29% (filter no. 2); and λc = 10500 ± 200 nm, τ = 76.6%
(filter no. 3). Although the separation method using a filter is relatively unrefined, it
could still provide a preliminary reference based on the experimental conditions [2].
The experiment used two Chinese coals with different coal ranks: Ximeng lignite
(XM) and Shenhua bituminous coal (SH). The particle size of the coal samples ranged
from 75 to 150 μm and were dried in a 105 °C oven for 2 h. Table 5.1 presents the
proximate and ultimate analyses.
In the experiment, the furnace was heated to the desired temperature, and the
reaction atmosphere was then passed through the tube at 1.5 ± 0.01 L/min for 30 min
to ensure the experimental atmosphere in the tube. An empty crucible was then
quickly pushed into the center of the reaction zone, and the data were subsequently
recorded as a calibration set. After calibration, a similar process was performed for
the crucible with a coal sample weighing 0.8 ± 0.005 g and evenly placed on the
bottom of the crucible. The crucible was specially designed to conduct better radiation
measurements (Fig. 5.1). Similar experiments were conducted for each filter. The
radiative heat flux was represented by the difference between the experimental and
calibration results. The experiment was conducted in a dark laboratory to minimize
the effect of the environment as much as possible.
5.1 Multi-scale experiments 81
In the experiment, oxygen ratios of 21–60% for the O2 /CO2 atmosphere at 800 °C
and a temperature range of 600–1000 °C in a 30% O2 /CO2 atmosphere were investi-
gated. These conditions are chosen because a 30% O2 /CO2 atmosphere is often used
as an alternative air-combustion condition [3]. Moreover, the selection of temper-
ature range mainly considers the temperature of the furnace flue gas for oxy-fuel
boiler is mainly around 1000 °C [4]. Therefore, this is still meaningful to the study
of oxy-fuel combustion characterized by flue gas radiation. In this study, the control
variable method is used to achieve the general significance of the comparison between
different conditions, and it has been discovered that the radiation flux changes in a
tube furnace [1]. To provide a reference for an actual combustion furnace, the average
value of the estimated radiation flux during the period from 80 to 140 s was used, as
the radiation flux variation was relatively stable in this stage.
The instrument uncertainty (3%) of the radiation sensor is considered when eval-
uating the uncertainty factors; thus, the measured value is multiplied by 3% as the
instrument error δ 1 . In addition, the standard deviation δ 2 is used to represent the
degree of fluctuation in the experimental data because the data is investigated over a
1-min period. Furthermore, the instrument uncertainty /(δ 1 and statistical
)
uncertainty
δ 2 are used to calculate the synthetic uncertainty δ = δ1 + δ2 [2], which is shown
2 2
Chap. 2:
( )
4 T0 1 T0 4
ηλ = 1 − + (5.1)
3 Tλ 3 Tλ
where T λ denotes the equivalent temperature and can be used to characterize radiation
exergy. Spectral radiation energies with different frequencies have different exergies,
so T λ is related to the wavelength and the radiation temperature. It is expressed as
f T3
Tλ4 = (5.2)
λ
Here, when T > T 0 , 4.5638 × 10−3 T 3 = λTλ4 ; when T < T 0 , 5.33 × 10−3 T 3 = λTλ4 .
In the experiment, radiative energy is divided into four different wavebands
through the three filters. When Eq. (5.1) is used to calculate the radiative exergy
for each waveband, a black-body equivalent wavelength λm is designed, and its
equivalent temperature is T λm at a certain temperature (Eq. 5.2). The exergy effi-
ciency expressed by T λm (Eq. 5.1) can be used to calculate the blackbody radiative
exergy E λ1−λ2 (T ) in the wavelength range of λ1 − λ2 :
[ ( ) ]
4 T0 1 T0 4
E λ1 −λ2 (T ) = L λ1 −λ2 (T ) × 1 − + (5.3)
3 Tλm 3 Tλm
where L λ1−λ2 (T ) denotes the blackbody radiation flux in this spectral waveband and
λm ranges between λ1 and λ2 . After obtaining the values of λm for the four wavebands
of 0.3–1.125, 1.125–4.1, 4.1–10.5, and 10.5–50 μm, Eq. (5.3) is used to calculate
the radiative exergy flux in each waveband based on the experimental radiation flux,
and the total radiative exergy flux is the sum of all fluxes.
In the laboratory of the current study, a Hencken flat-flame burner for diffusion
combustion was built based on the design in [5, 6]. The flat-flame burner uses gas
fuel combustion to produce flue gas with a specific composition and temperature to
support solid fuel combustion. The particle jet is ignited using the high-temperature
flue gas produced by gas fuel combustion and the particle heating rate under this
experimental condition are similar to that of actual furnace [6]. The detailed descrip-
tion of the Hencken flat-flame burner is contained in Sect. 5.1.2.2. Moreover, Fig. 5.2
shows the experiment system, which mainly includes a Hencken burner, a powder-
feeding device, a gas-supplying subsystem, and an optical measurement subsystem.
84 5 Radiative Energy Characteristics of Solid Fuel Combustion
Fig. 5.2 Diagram of the experimental system based on Hencken flat-flame burner
Under the effect of the ultrasonic vibrator and the stirring motor, the pulverized fuel is
carried by air from the powder bin through the stainless-steel tube to the central burner
outlet, which is surrounded by gas outlets of carbon monoxide/propane evenly. The
high-temperature flue gas produced by gas combustion provides required temperature
and atmosphere for the combustion of solid fuel.
All experimental gases are controlled by mass-flow meters (Alicat Co.). The fuel
gas and combustion-supporting atmosphere are mixed by the gas mixer and then
enters the burner. Moreover, the AvaSpec-ULS2048XL spectrometer is used to obtain
the radiation spectrum of the target flame with a wavelength range of 200−1160 nm
(spectral resolution is 2.4 nm), and the diameter of receiving area is 1 mm for the
fiber probe. The spectrometer is calibrated by the standard blackbody furnace signal
before this experiment, so it can be used directly to obtain the radiation spectrum.
To facilitate spectrum measurement at different heights, the optical fiber probe is
connected to a collimator and fixed on the optical bracket, which is fixed on the slide
rail. The test platform is also calibrated to maintain its horizontal position, and the
experiment is conducted in a dark room to avoid interference from other visible light.
Before data collection, the probe position must be adjusted using laser alignment to
determine the flame center position. Continuous monitoring is conducted when the
flame is stable, and 50 consecutive sets of spectral data are chosen for averaging to
obtain the results. To achieve real-time spectral data processing, the probe is linked
to the computer via optical fiber. The signals at 7, 10, 14, 17, and 20 cm from the
burner outlet are measured in the experiment. Because the measuring points are at
5.1 Multi-scale experiments 85
actual height, the measured results can represent the radiation intensity in the same
radiation area, which is better for comparing different cases.
The Hencken burner can generate a thin-plane flame (2–3 mm) for gas combustion
and is also flexible for adjusting the flue gas composition and temperature. Further-
more, it is more suitable for optical measurement because of the absence of the need
for electric or other heating devices.
The cross-sectional size of the self-built Hencken burner is about 55 × 55 mm
(Fig. 5.2), and it can form a gas flame plane with a corresponding area. The burner
is primarily made of upper and lower chambers. The lower chamber has a height of
about 100 mm for fuel gas, while the upper chamber has a height of about 200 mm
for combustion-supporting gas. The upper and lower chambers are separated by heat-
resistant silica-gel plates. Two holes were noted at the bottom of the lower chamber
for the gas introduction. The gas fuel in the experiment is carbon monoxide. The upper
space of the gas inlets is filled with color-changing colloidal silica for adjusting and
drying the gas flow. In the upper chamber, 652 stainless-steel capillary tubes serve as
fuel channels, penetrating the silicone plate into the lower chamber to deliver fuel gas
to the burner outlet. Two symmetrically designed holes were noted in the lower part
of the upper chamber to input combustion-supporting gas (e.g., air), which enters
the upper chamber and flows outside the capillary tubes until merging with the fuel
gas at the burner outlet. The water-cooling groove is located on the periphery of
the upper chamber, and a thicker powder-feeding tube (outer diameter, 3 mm; inner
diameter, 2.5 mm; and total length, about 330 mm) is installed through the entire
central section of the burner. A quartz glass tube is placed as a protective cover in the
upper space of the burner outlet to prevent the flue gas temperature from dropping
too quickly and to suppress the influence of external factors during the experiment.
Furthermore, the Hencken burner powder-feeding device consists primarily of a
powder bin, a stepping motor, and an ultrasonic vibrator with a controller. The solid
fuel is carried to the inlet of the powder-feeding tube at the bottom of the burner
by adjusting the carrier air. A small amount of propane is added with a flow rate of
20 mL/min (SCCM) because the gas fuel CO is difficult to burn consistently at high
gas flow rates.
The hot flue gas produced by gas combustion in the flat-flame burner can provide the
temperature and atmosphere for solid fuel combustion. The atmosphere produced
by gas combustion is referred to as secondary airflow. The flow rate of each gas
entering the burner is reversely deduced based on secondary air conditions to obtain
the designed experimental conditions of the secondary airflow (temperature, oxygen
ratio, and gas flow rate). The appropriate method is as follows:
86 5 Radiative Energy Characteristics of Solid Fuel Combustion
Since the main fuel gas in this study is CO and the atmosphere is assumed to
be O2 /N2 , the following reaction equation was used to calculate the adiabatic flame
temperature
where a, b, and c are the flow rates of CO, O2 , and N2 , respectively, in L/min (SLM).
According to the reaction in Eq. (5.4), if the flow rates of each reactant are known,
the adiabatic flame temperature T ad can be calculated by referring to the specific
heat capacity, enthalpy, and temperature of gases. Furthermore, the conservation of
material mass before and after the reaction can be demonstrated as follows:
where ρ is the density of the mixed gas, in kg/m3 ; V is the gas velocity on the burner
section, in m/s; S is the effective area of the burner outlet and can be estimated by S
= 0.8D2 , where D = 55 mm and is the side length of the outlet section. Using the
ideal state equation Pm = NρRT and the conservation equation, the gas velocity can
be calculated as follows:
( / )
a 2 + b + c Tad
V = (5.6)
ST0
b − a/2
YO2 = × 100% (5.7)
a/2 + b + c
Table 5.3 Proximate and ultimate analysis of the semi-coke from bituminous coal
Proximate analysis % Calorific value Qnet (MJ/kg) Ultimate analysis %
Aad V ad FCad C ad H ad Oad N ad S t,ad
9.47 5.76 84.77 28.71 85.04 0.19 4.51 0.16 0.63
Note ad air dry basis; Qnet net calorific value
The radiation spectrum was collected by spectrometer after the solid fuel flame was
stable.
The best method for measuring the temperature of the solid fuel flame is spec-
troscopy. When compared with contact measurement, this spectroscopy method does
not disturb the flame shape, which prevents the experimental results from being influ-
enced. According to the spectral signal directly measured by the spectrometer, the
two-color method coupled with the emissivity ratio model developed by the labo-
ratory of the current study is used to calculate the temperature of the semi-coke
flame. Temperature measurement is not affected by the measuring distance with this
method, and it includes the variation of spectral radiation with wavelength, yielding
a more accurate result of the temperature of burning particles. The specific principle
can be observed in [7].
This study investigates the flame radiation exergy-based experimental data and
the radiation thermodynamics theory, which provides a reference for the energy
cascading conversion system of fuel combustion, to characterize the available radi-
ation energy. Spectral radiation exergy can be obtained by multiplying the radiative
energy by a radiative exergy-to-energy ratio ηλ as shown in Eqs. (5.1) and (5.2).
The spectrum for particle radiation in solid fuel combustion is continuous and can
be approximately regarded as a gray body [8, 9] in the waveband of < 1.1 μm. In
the waveband of < 2.25 μm, no significant variation was also noted in the spectral
emissivity with temperature [10]. Spectrum continuity and emissivity stability in the
long waveband for solid vegetation combustion were also demonstrated in [11]. In a
preliminary assessment, the linear extrapolation method is used in [12] to predict the
spectral emissivity of the infrared waveband for solid fuel combustion. In this study,
88 5 Radiative Energy Characteristics of Solid Fuel Combustion
The radiative energy and exergy flux for the two coals in the O2 /CO2 atmospheres
with different oxygen ratios are shown in Fig. 5.3a, b. For the two coals, the radiative
energy flux increases from about 2 to 9 W/m2 with an increase in the oxygen ratio
from 21 to 60%; this occurs because the combustion rate increases with the oxygen
ratio, releasing more radiative energy. Furthermore, higher local temperatures caused
by higher oxygen ratios increase radiation flux. The radiative energy flux of the
bituminous coal is greater than that of lignite. The radiative energy flux may be
affected by both volatilization and calorific value. Bituminous coal and lignite burn
more easily and have higher radiation fluxes because they have higher volatilization.
Furthermore, its radiative energy flux is slightly higher than that of lignite because
bituminous coal has the highest calorific value. In addition to the radiative energy,
Fig. 5.3a, b show that the radiative exergy flux increases from approximately 1.5
to 6 W/m2 as the oxygen ratio increases from 21 to 60%. The radiative exergy flux
of the bituminous coal is greater than the lignite; this trend is similar to that of the
radiative energy flux. Since the furnace remains at 800 °C, the results indicate that
the trend of the radiative exergy flux as a function of the oxygen ratio is consistent
with that of the radiative energy flux when the combustion temperature is constant.
A comparison of the exergy-to-energy ratio illustrates that it is stable between 0.61
and 0.63. This indicates that the radiative exergy-to-energy ratio is generally similar
at a specific combustion temperature. Theoretically, as the oxygen ratio increases,
the proportion of radiative energy in the short-wavelength range increases, and thus
the total radiative exergy increases.
Figure 5.3c, d show the radiative energy and exergy flux of two different coals
at temperatures ranging from 600 to 1000 °C in a 30% O2 /CO2 atmosphere. The
temperature has a significant impact on the radiation flux. The radiative energy flux
of the bituminous coal and lignite increases from approximately 2–14 W/m2 as the
temperature increases from 600 to 800 °C. According to Stefan-Boltzmann’s law, the
radiative power is generally proportional to the fourth power of the temperature; thus
the radiative energy flux increases rapidly with the temperature [13]. The radiation
flux of the bituminous coal is similar to that of the lignite. This is due to the use
5.2 Experimental Results and Discussion 89
12 0.630 12 0.630
(a) Lignite (XM) (b) Bituminous (SH)
10 Energy flux 10 Energy flux
Exergy flux 0.625 Exergy flux 0.625
Exergy-to-energy ratio
Exergy-to-energy ratio
Radiation flux (W/m2)
6 0.620 0.620
6
4
0.615 4 0.615
2
2
0 0.610 0.610
20 30 40 50 60 20 30 40 50 60
Oxygen ratio (%) Oxygen ratio (%)
18 0.70 18 0.70
16
(c) Lignite (XM) 0.68 16
(d) Bituminous (SH) 0.68
Energy flux 0.66 Energy flux 0.66
14 14
Exergy flux
Exergy-to-energy ratio
Exergy-to-energy ratio
Radiation flux (W/m2)
Bituminous-Exergy Lignite-Energy
0.6 Bituminous-Energy
0.4
Lignite-Exergy
0.3 Bituminous-Exergy
0.4
0.2
0.2
0.1
0.0 0.0
20 30 40 50 60 600 700 800 900 1000
Oxygen ratio (%) Temperature (°C)
Fig. 5.3 Total radiative energy and exergy flux during oxy-coal combustion for a lignite and b
bituminous coal with different oxygen ratios and c lignite and d bituminous coal with different
temperatures. Radiation fraction and radiative exergy fraction during oxy-coal combustion with e
different oxygen ratios and f different temperatures
of highly volatile bituminous coal, which, like lignite, has a high burning intensity.
Thus, their local temperatures are higher, and the radiative energy flux is greater;
additionally, the combustion temperature influences the radiation flux significantly.
Figure 5.3c, d also show that the radiative exergy flux of the lignite and bituminous
coal increases from about 1 to 10 W/m2 with the increase in temperature from 600 to
1000 °C. The radiative exergy flux follows the same trend as the radiative energy as
temperature rises. A comparison of the exergy-to-energy ratio shows that while the
magnitudes of the radiative energy and exergy flux differ, the exergy-to-energy ratios
of the two coals increase almost linearly from 0.54 to 0.68 with the temperature, and
little difference exists between the two coals. This result indicates that the radiative
exergy-to-energy ratio is essentially similar at a specific combustion temperature.
90 5 Radiative Energy Characteristics of Solid Fuel Combustion
Furthermore, the coal ranks have little influence on the exergy-to-energy ratio with
the temperature being the primary influencing factor.
The radiation and radiative exergy fractions are also investigated in this study. The
radiation fraction is defined as the ratio of the radiation power released during the
combustion process to the input power of the fuel, and the radiative exergy fraction
is defined as the ratio of the radiative exergy power to the input exergy power of the
fuel. Thus, the radiation and radiative exergy fractions can be calculated following
the definition, and their variations are then investigated under different conditions.
The chemical exergy of coal is calculated using the approximate formula shown
below:
[ ]
n(H) n(O) n(N)
E fuel = L H V 1.0064 + 0.1519 + 0.0616 + 0.0429 (5.8)
n(C) n(C) n(C)
where LHV is the low heat value, and w(C), w(H), w(O), and w(N) are the mass
fractions of the elements in the coal. The equivalent reaction time is set to 160 s
because the tube furnace is fed with coal at one time, and the input powers of different
coals are calculated according to the LHV. Assuming that the radiation basically
conforms to Lambert’s law, it can be supposed that the received radiation flux is
consistent in the entire hemisphere with the radius of the distance from the sample to
the sensor. Therefore, the output radiation power of fuel combustion can be calculated
by multiplying the hemisphere area by the detected radiation heat flux.
Despite the relevant assumptions, it is believed that the estimated results could still
reflect the relative variation of the radiation fraction under different conditions, which
differs from the radiative power variation shown in Fig. 5.3a–d. Figure 5.3e, f show
how the radiation fraction and radiative exergy fraction change with oxygen ratio
and temperature during the combustion process. As the oxygen ratio and temperature
increased, so did the radiation fraction and radiative exergy fraction. These findings
suggest that increasing the oxygen ratio and combustion temperature in the coal
combustion process can increase the proportions of radiative energy and exergy to
the total energy, which is beneficial to the graded conversion of radiative energy
during combustion. Furthermore, the variation characteristics of the radiation and
radiative exergy fractions are similar in different atmospheres. The radiative exergy
fraction is less than the radiation fraction, indicating that the quality of fuel chemical
energy is higher than that of radiative energy, and combustion is an energy-losing
process. A more efficient energy cascade conversion system must take chemical
energy conversion into account. Furthermore, lignite has been discovered to have
higher radiation and radiative exergy fractions than bituminous coal. This indicates
that lignite, which burns more easily, emits more radiation during the combustion
process. Bituminous coal is commonly used in power plant boilers because lignite
has a lower calorific value and anthracite is difficult to burn. Highly volatile solid
fuels (e.g., lignite) could be considered for combustion energy cascade conversion
systems.
5.2 Experimental Results and Discussion 91
The spectral radiation results are discussed in this section. The spectral radiation
proportion is defined as the ratio of spectral radiation flux in a waveband to total
radiation flux. Similarly, the spectral radiative exergy proportion is defined as the
ratio of spectral radiative exergy flux in a waveband to the total radiative exergy flux.
Thus, they can be computed using the definition. Figure 5.4a, b show how the spectral
radiative energy and exergy change with the oxygen ratio during combustion in an
O2 /CO2 atmosphere at 800 °C temperature. The radiative energy ratio is highest in
the waveband of 1.125–4.1 μm and increases from about 40 to 50% with the oxygen
ratio. Additionally, the radiative energy ratio is approximately 40% in the waveband
of 4.1–10 μm, whereas it is about 10% in the spectrum of > 10 μm. The experimental
data for radiation of < 1.125 μm are close to zero due to the low energy proportion
in this spectrum. Furthermore, fluctuations affect the results. Thus, these data should
only be used as a reference.
The participating media influence the spectral radiation distribution. The energy
proportion in the far-infrared spectrum of the actual radiation is greater than that of
gray-body radiation due to the presence of CO2 and H2 O in the combustion process
[13]. The radiation proportion in the waveband of 1.125–4.1 μm is found to increase
60 60
(a) Lignite (XM) 1.125-4.1 µm (b) Bituminous coal (SH) 1.125-4.1 µm
Spectral radiation proportion (%)
50 50
40 40
energy 0.3-1.125 µm exergy 0.3-1.125 µm energy 0.3-1.125 µm exergy 0.3-1.125 µm
energy 1.125-4.1 µm exergy 1.125-4.1 µm energy 1.125-4.1 µm exergy 1.125-4.1 µm
30 30
energy 4.1-10.5 µm exergy 4.1-10.5 µm energy 4.1-10.5 µm exergy 4.1-10.5 µm
energy 10.5-50 µm exergy 10.5-50 µm energy 10.5-50 µm exergy 10.5-50 µm
10 10
0 0
20 30 40 50 60 20 30 40 50 60
Oxygen ratio (%) Oxygen ratio (%)
70 70
(c) Lignite (XM) (d) Bituminous coal (SH)
Spectral radiation proportion (%)
60 1.125-4.1 µm 60 1.125-4.1 µm
50 50
40 40
0 0
600 700 800 900 1000 600 700 800 900 1000
Temperature (°C) Temperature (°C)
Fig. 5.4 Spectral radiative energy and exergy during oxy-coal combustion versus the oxygen ratio
for a lignite and b bituminous coal versus the temperature for a lignite and b bituminous coal
92 5 Radiative Energy Characteristics of Solid Fuel Combustion
with the increase in the oxygen ratio, whereas it decreases in the 4.1–10.5 and 10.5–
50 μm wavebands. This occurs because a higher oxygen ratio increases the burning
rate and the local combustion temperature, causing the spectral distribution to shift
to shorter wavelengths following Wien’s displacement law [13]. Furthermore, the
CO2 decrease in the atmosphere caused by an increase in the oxygen ratio reduces
the shift to longer wavelengths. Bright flames are observed during the combustion
of lignite and bituminous coal. This is attributed to the fact that bituminous coal and
lignite burn intensely and have higher burning temperatures [14]; thus, Fig. 5.4a,
b reveal higher radiation proportions in the waveband of < 4.1 μm for lignite and
bituminous coal. Furthermore, radiation accounted for a greater proportion of the
waveband of > 4.1 μm during lignite combustion because lignite contains more
moisture than bituminous coal; however, this effect is weak. In addition, Fig. 5.4a, b
show that the trend of the radiative exergy flux with increasing oxygen ratio is similar
to the trend of the radiative energy flux in various bands. This is primarily because
radiative exergy is derived from radiative energy. Moreover, the radiative exergy ratio
was also discovered to be greater than the radiative energy ratio in the waveband of
< 4.1 μm, but less than the radiative energy ratio in the waveband of > 4.1 μm. This
indicates that the exergy-to-energy ratio is greater for shorter-wavelength radiation
of < 4.1 μm, implying that shorter wavelengths are more useful.
Figure 5.4c, d depict the spectral radiative energy and exergy flux versus temper-
ature in a 30% O2 /CO2 atmosphere. In terms of the spectral radiative energy flux
characteristics, the radiative energy ratio in the waveband of 1.125–4.1 μm increases
from about 30 to 60% with the increase in the temperature from 600 to 1000 °C,
whereas that in the waveband of 4.1–10.5 μm decreases from approximately ~50 to
35%. Furthermore, the radiative energy ratio in the spectrum of > 10.5 μm falls from
approximately 15 to 5%. The radiation proportion in the waveband of < 1.125 μm
should theoretically increase with temperature. This trend, however, is not observed
in Fig. 5.4c, d because the energy ratio in this band is small and the experimental
fluctuation has a large influence on the results. These findings suggest that temper-
ature is the most important factor influencing the spectral distribution of radiative
energy in pulverized coal combustion; this is explained by Wien’s displacement law,
which states that as temperature rises, spectral radiation shifts to a shorter wavelength
[13]. Generally, radiation transfer occurs with participating media in pulverized coal
combustion. Because of the presence of soot, CO2 , and H2 O, it cannot be simply
regarded as gray-body radiation [13]. Nonetheless, the spectral radiation distribu-
tion characteristics still follow the basic law of radiation; thus, in this experiment,
the proportion of radiative energy increases in the spectrum range of < 4.1 μm and
decreases in the spectrum range of > 4.1 μm. Furthermore, as shown in Fig. 5.4c,
d, the different coals do not exhibit different influences on the spectral radiation as
temperature increases. Coal combustion is complete at temperatures of > 800 °C
[15]; so the spectral radiation characteristics during combustion for the two coals at
high temperatures are essentially similar. The trend of the radiative exergy flux with
increasing temperature is similar to that of the radiative energy in each waveband.
Furthermore, the radiative exergy ratio is greater than the radiative energy ratio in
the waveband of < 4.1 μm, while it is lower than the radiative energy ratio in the
5.2 Experimental Results and Discussion 93
waveband of > 4.1 μm; this is consistent with the results in Fig. 5.4a, b. These
findings suggest that radiation with short wavelengths of < 4.1 μm has a higher
exergy-to-energy ratio and thus has greater practical significance. Thus, increasing
the short-wavelength radiative energy ratio in combustion radiation can improve
overall radiative exergy, which is important for energy conversion and management.
Table 5.4 displays the temperature measurements made using the improved two-color
method based on spectral data. The spectral radiation that peaks near 589 and 769 nm
in the flame radiation spectra are primarily caused by the characteristic radiation of
the Na and K elements [15]. The waveband affected by characteristic radiation is thus
removed when calculating the temperature. The measured temperature is observed
to deviate from the design temperature near the burner outlet. At the measuring point
of 7 cm, for example, it is 1519.37 and 1602.71 °C in cases 2 and 3, which differ
from the design temperature of 1400 °C. This is primarily due to the deviation of
secondary airflow from ideal conditions and the temperature increase in the particles
caused by the semi-coke combustion. In addition, it may also be related to the signal-
to-noise ratio of the data collection and the model-matching degree of the two-color
method.
The measured temperature at 20 cm above the burner outlet is similar to the
design temperature, indicating that the combustion reaction has a significant impact
on the flame temperature. Despite being influenced by the experimental conditions,
the temperature results are still representative of the 10- to 14-cm middle section.
The temperature of case 2 is similar to that of case 3, close to the setting temperature
of 1400 °C, whereas the temperature of case 1 is lower due to the lower setting
temperature. Moreover, the temperature in case 4 is higher than that in case 3 by
approximately 100 °C, indicating a higher actual temperature for a higher oxygen
ratio when the setting temperature is similar.
Figure 5.5 depicts the total radiative energy and exergy along the flame under
various experimental conditions. The measurement’s spectral range of total radiative
energy is 500–1100 nm, covering visible and near-infrared wavebands. Although it
Table 5.4 Measured temperature by the two-color method based on radiation spectrum
Temperature °C 7 cm 10 cm 14 cm 17 cm 20 cm
Case 1 1517.36 1501.52 1425.56 1360.51 1218.18
Case 2 1519.37 1482.20 1463.70 1450.99 1428.38
Case 3 1602.71 1464.27 1450.05 1420.93 1420.25
Case 4 1877.65 1690.17 1597.56 1527.61 1451.74
94 5 Radiative Energy Characteristics of Solid Fuel Combustion
accounts for a small proportion of the total thermal radiation spectrum, it is extremely
useful due to its short-wavelength. This waveband is important for studies because
it is suitable for photovoltaic conversion of commercial Si cells. Moreover, Fig. 5.5
shows that the radiative energy intensity under the four cases increases first and
then decreases from the burner outlet upwards. A slight fluctuation exists in case 3
as shown in Fig. 5.5c, which is due to measurement error, but it is roughly in line
with the overall trend. These findings suggest that the radiation intensity distribution
along the flame is related to the combustion reaction intensity, showing a trend
of first increasing and then decreasing and presenting a cone-shaped flame shape.
Furthermore, the highest radiation intensity is found to be at 17, 10, 17, and 15 cm
in cases 1, 2, 3, and 4, respectively. This indicates that ignition occurs first in case
2, second in case 4, and finally in case 1 or 3. The radiation intensity peak appears
first in case 2 because it has the highest setting temperature (1400 °C). Furthermore,
the diffusion rate of O2 molecules is faster in the O2 /N2 atmosphere than in the
O2 /CO2 atmosphere, resulting in a faster oxidation rate of semi-coke and an earlier
ignition. Because the ignition occurs late in the O2 /CO2 atmosphere, the radiation
intensity peak in cases 3 and 4 occurs later than in case 2. Furthermore, because
the O2 ratio is higher in case 4, the radiation intensity peak arises earlier than that
in case 3; additionally, the radiation intensity peak in case 1 appears later because
its design temperature is lower. This analysis shows that the reaction atmosphere,
temperature, and oxygen ratio all affect the radiation intensity peak and ignition. The
O2 /CO2 atmosphere retards ignition, whereas high-temperature and high oxygen
ratio conditions promote early ignition and cause the radiation intensity to peak
earlier.
Furthermore, the highest radiation intensity in case 2 is higher than that in case
1 by about 700 W/m2 , indicating that high-temperature directly increases radiation
intensity; whereas the highest radiation intensity in case 4 is higher than that in
case 3 by about 350 W/m2 , indicating that increasing oxygen ratio will also increase
flame radiation intensity. Furthermore, the peak radiation intensity in case 2 is about
600–950 W/m2 higher than in cases 3 or 4, indicating that when the parameters
(e.g., oxygen ratio or flue gas temperature) are similar under different atmospheres,
the radiation intensity in the O2 /CO2 atmosphere is lower than that in the O2 /N2
atmosphere. This is primarily due to the lower combustion intensity caused by the
slower diffusion of O2 in the CO2 and the low local temperature caused by the higher
specific heat capacity of CO2 . Furthermore, the gasification of semi-coke and the
slow oxidation reaction rate caused by the addition of CO2 will affect the properties
and distribution of combustion potentially increasing medium absorption along the
path length and resulting in lower measured radiation intensity.
Finally, Fig. 5.5 shows that the total radiation exergy curve along the flame path
length is consistent with the radiation energy curve, and its absolute value is less than
the radiation energy. Furthermore, higher radiation intensity results in a greater differ-
ence between radiative energy and radiative exergy. Further investigation reveals that
the exergy-to-energy ratios of the four cases are all around 82%, where it is higher
and lower in cases 4 (82–85%) and case 1 (80–82.5%), respectively. The measured
temperature of each point in case 4 is higher while it is lower in case 1. The variation
5.2 Experimental Results and Discussion 95
450
(c) Case 3 800
(d) Case 4
85.0
50 81.5 0 82.0
6 8 10 12 14 16 18 20 22 6 8 10 12 14 16 18 20 22
Measuring point height (cm) Measuring point height (cm)
Fig. 5.5 Total radiative energy and radiative exergy along the flame under different cases
curve of the radiation exergy-to-energy ratio under the four cases is similar to the
measured temperature. It implies that the radiation temperature is the most important
factor influencing radiation exergy, whereas the influence of the spectral waveband
is not obvious in Fig. 5.5.
Particle radiation dominates in solid fuel combustion, and a continuous spectrum was
noted [16, 17]. Based on the measured emissivity trend in the near-infrared spectrum,
the spectral emissivity is linearly extrapolated to the wavelength of 15 μm. Figure 5.6
show the spectral radiative energy and exergy characteristics in a broad spectrum
(cases 2 and 4 are selected). The spectral radiative energy is primarily concentrated
in the 1.1–3 μm band at the current flame temperature (~ 1200–1800 °C). The spec-
tral radiation in the wavelength of < 3 μm is useful in actual photovoltaic conversion
because the Si cell match the spectral radiation of < 1.1 μm and the III–V semi-
conductor cells (e.g., InAsSbP) matches the spectral radiation of < 3 μm [18]. High
temperatures are advantageous for the utilization of radiative energy in combustion.
The spectral radiative energy in the waveband of < 3 μm along the path length is
similar to the variation trend in Fig. 5.5 (0.3–1.1 μm) as shown in Fig. 5.6. The
spectral radiative energy in the waveband of > 3 μm shows an increasing trend at
17 and 20 cm. The radiation in waveband of > 3 μm accounts for more because the
temperatures at 17 and 20 cm are lower. The spectral distributions of the radiation
96 5 Radiative Energy Characteristics of Solid Fuel Combustion
exergy under different conditions are similar to those of radiative energy because
the radiation exergy is calculated from the radiation energy. Figure 5.6e depicts the
exergy-to-energy ratio of spectral radiation in cases 2 and 4, demonstrating that
the spectral exergy-to-energy ratio is greater for shorter wavelengths in each case
which could be between 76 and 85% for the waveband of < 3 μm. The measured
temperature influences the difference in the exergy-to-energy ratio at different points.
Furthermore, the exergy-to-energy ratio is greater at higher temperatures. These find-
ings suggest that short-wavelength and high-temperature radiation energy are more
valuable.
Even when professional fixtures (e.g., optical brackets and moving rails) are used
to assist spectrometer measurement, a deviation in the measured radiation value
14 14
(a) Spectral energy in Case 2 (b) Spectral exergy in Case 2
Spectral radiation power (kW/m2)
10 10
7cm
7cm 10cm
8 10cm 14cm
8
14cm 17cm
6 17cm 6 20cm
20cm
4 4
2 2
0 0
0--1.1 1.1--3 3--6 6--10 >10 0--1.1 1.1--3 3--6 6--10 >10
Waveband (µm) Waveband (µm)
7 7
(c) Spectral energy in Case 4 (d) Spectral exergy in Case 4
Spectral radiation power (kW/m2)
6 6
7cm
5 7cm 5 10cm
10cm 14cm
4 4
14cm 17cm
3 17cm 3 20cm
20cm
2 2
1 1
0 0
0--1.1 1.1--3 3--6 6--10 >10 0--1.1 1.1--3 3--6 6--10 >10
Waveband (µm) Waveband (µm)
0.90 0.90
(e) Case 2 0--1.1 µm 1.1--3 µm (f) Case 4 0--1.1 µm
3--6 µm
1.1--3 µm
6--10 µm
3--6 µm 6--10 µm
Spectral exergy-to-energy ratio
0.80 0.80
0.75 0.75
0.70 0.70
0.65 0.65
0.60 0.60
6 8 10 12 14 16 18 20 6 8 10 12 14 16 18 20
Measuring point height(cm) Measuring point height (cm)
Fig. 5.6 a–d Spectral radiative energy and exergy distribution and e and f spectral exergy-to-energy
ratio for different cases
5.2 Experimental Results and Discussion 97
under different conditions still exists due to the influence of the measurement envi-
ronment. The absolute value of the spectral radiative energy in Fig. 5.6 provides
only a limited amount of information. The proportion of spectral radiative energy
and exergy in various cases is discussed further. The proportion characteristics have
more significance for reference than the absolute value (Fig. 5.7). In different cases,
the largest proportion can be observed to be in a waveband of 1.1–3 μm. The propor-
tion of spectral radiative energy of < 3 μm wavelength decreases with distance
from the burner outlet, whereas the proportion of > 3 μm wavelength increases. The
spectral radiation distribution shifts to the long-wavelength direction according to
Planck’s law because the flame temperature decreases with distance; thus, the propor-
tion of the long-wavelength radiation increases while that of the short-wavelength
radiation decreases. Emissivity decreases with temperature in the experiment; this
phenomenon is more visible in the measured short-wavelength and theoretically
weakens the increasing short-wavelength proportion. Indeed, the results show that
increasing the temperature does not affect the short-wavelength proportion of the
measuring point. However, the proportion of short-wavelength radiation (< 3 μm) is
still greater than that of long-wavelength radiation for a given measuring point spec-
trum. This means that the effect of temperature on spectral emissivity is insufficient
to change the overall trend of spectral distribution with temperature. Furthermore,
Fig. 5.7 shows that the proportion of the spectral radiative energy under different
cases is relatively stable, and the variation in each waveband is within 5%. Case
4 has a higher proportion of radiative energy in the waveband of 0.5–1.1 μm than
the other cases because its measured temperature is higher, and case 2 has a higher
proportion of energy in the wavelength of < 3 μm than case 1. These characteris-
tics are especially noticeable at 20 cm, indicating that increasing the combustion
temperature may increase the proportion of short-wavelength radiation. Case 4 has
a higher proportion of energy in the wavelength of < 3 μm than in case 3, indicating
that increasing the oxygen ratio may increase the proportion of short-wavelength
radiation. Cases 2 and 3 have very similar spectral energy distribution characteris-
tics, indicating that the oxy-fuel atmosphere has little effect on the spectral radiation
proportion under similar temperature conditions. Particle radiation dominates solid
fuel combustion in an oxy-fuel atmosphere, while radiation from gases and other
media accounts for a smaller proportion [16, 17]. CO2 gasification and its high
specific heat could affect particle combustion, affecting local particle temperature
and particle emissivity. Because temperature is the most influential factor influ-
encing combustion radiation, increasing the oxygen ratio to raise the temperature is
beneficial for increasing short-wavelength radiation, which is important for radiative
energy utilization. The results of radiation exergy proportion are consistent with that
of radiative exergy in various cases as shown in Fig. 5.7e, f. Thus, the radiative energy
proportion can be approximated by the radiation exergy proportion. The proportion
of radiative exergy in the waveband of < 3 μm is slightly higher than that of radiative
energy, whereas it is slightly lower in the waveband of > 3 μm. This is because
high-frequency radiative energy has a higher equivalent radiation temperature and
its radiative exergy is relatively higher.
98 5 Radiative Energy Characteristics of Solid Fuel Combustion
70 70
(a) Spectral energy in Case 1 (b) Spectral energy in Case 2
60 60
50 7cm
50
7cm 10cm
40 10cm 40 14cm
14cm 17cm
30 17cm 30 20cm
20cm
20 20
10 10
0 0
0--1.1 1.1--3 3--6 6--10 >10 0--1.1 1.1--3 3--6 6--10 >10
Waveband (µm) Waveband (µm)
70 80
(c) Spectral energy in Case 3 (d) Spectral energy in Case 4
60 70
7cm 60
50
10cm 50 7cm
40 14cm 10cm
17cm 40
14cm
30 20cm 30 17cm
20 20cm
20
10 10
0 0
0--1.1 1.1--3 3--6 6--10 >10 0--1.1 1.1--3 3--6 6--10 >10
Waveband (µm) Waveband (µm)
70 80
(e) Spectral exergy in Case 2 (f) Spectral exergy in Case 4
Spectral radiation proportion (%)
60 70
60
50 7cm
7cm
50 10cm
40 10cm
14cm 40 14cm
30 17cm 17cm
30 20cm
20cm
20
20
10 10
0 0
0--1.1 1.1--3 3--6 6--10 >10 0--1.1 1.1--3 3--6 6--10 >10
Waveband (µm) Waveband (µm)
Fig. 5.7 Proportion distribution of radiative energy and radiative exergy under different experi-
mental cases
One-dimensional cases are built for media radiation transfer calculation based on
the parameters of the tube furnace experiment (e.g., furnace temperature and reac-
tion conditions). The spectral radiative energy characteristics are predicted via one-
dimensional cases, and the calculated and experimental results are compared. The
length of the quartz tube determines the size of the one-dimensional case, which
is set at 60 cm. The thermocouple temperature curve in the furnace has a trape-
zoidal distribution. The temperature formulas in Table 5.5 are used for the cases to
match the temperature distribution in the furnace. The calculated temperature has a
5.3 Numerical Verification 99
similar curve bend to the measured temperature. This study is primarily concerned
with the spectral radiation proportion; thus, it should ensure that the highest and
lowest temperatures of the one-dimensional temperature field are consistent with
the measured values and that the designing temperature curve is also as close to the
measured value as possible along the entire length. This ensures the model’s ratio-
nality and calculation precision, allowing the model to be applied to the calculation
cases.
The calculation cases are based on bituminous coal (SH) combustion and a 30%
O2 /CO2 atmosphere. The effective reaction time is set to 160 s based on the experi-
mental results. The flue gas atmosphere of the cases is calculated using the parame-
ters of airflow and coal characteristics under the assumption of complete combustion.
The calculated H2 O/CO2 molar ratio is 0.022 and the proportion of triatomic gas is
0.75. The gas radiation properties are described using the WSGG model proposed in
Chap. 3. Coal combustion generates a large number of particles and is an important
radiation medium; thus, particle radiation is taken into account in the calculation
cases. The parameters of the particles during combustion are determined using the
results reported in the existing literature and summarized in Table 5.6.
In Table 5.6, absorption coefficient κ and scattering coefficient σ is calculated
according to the formulas [17, 19]:
κ = Q abs Ap (5.9)
σ = Q sca Ap (5.10)
For the total particles, κ part = κ coal + κ ash = 0.8 and σ part = σ coal + σ ash = 1.
The radiation parameters are also highest because the particle concentration is
highest at the combustion center. Furthermore, the radiation parameters are zero
because the minimum concentration at both ends is zero. Consequently, the radiation
parameter curve in the one-dimensional case is expressed by the following equations.
( )
2π x
κpart = 0.4 − 0.4 cos (5.11)
L
( )
2π x
σpart = 0.5 − 0.5 cos (5.12)
L
where g is set as 0.86 based on a previous study [19]. The radiation characteristics of
different temperature conditions in a 30% O2 /CO2 atmosphere are considered in the
case study, and the calculation is similar to that used in Chap. 3. The proportional
characteristics of spectral radiative heat flux are investigated here. The radiation spec-
tral proportion is unaffected by the absolute radiation intensity and can be compared
to experimental results.
Furthermore, the same analysis methods and parameters used in the flat-flame
burner experiment for semi-coke particle parameters are used here. The semi-coke
particles radiation parameters are κ part = 0.8, σ part = 1. The particle concentration
is assumed to be the greatest at the center of the flame, so the radiation parameter is
greatest, and the concentration at the two walls is the lowest to zero, so the radiation
parameter is also zero. Moreover, Eqs. (5.11) and (5.12) continue to represent the
radiation parameter distribution curves. The scattering phase function also chooses
the Henyey–Greenstein distribution as expressed in Eq. (5.13), with g = 0.86.
First, the comparison of the spectral radiation proportions calculation results in one-
dimensional cases with the tube furnace experimental results is discussed. The vari-
ations in the predicted spectral radiation proportion and the experimental results
versus the temperature under a 30% O2 CO2 atmosphere are similar, and the calcu-
lated results agree well with the experimental results for various cases at different
temperatures. The maximum difference between calculated and experimental results
is ~ 3% for the case at the temperature of 600 °C, whereas the differences are less
with a maximum of ~ 2%, for the three respective cases at temperatures of 700, 800,
and 900 °C. Furthermore, the difference is greater for the case with a temperature of
5.3 Numerical Verification 101
1000 °C, with a maximum difference of about 5%. This could be due to the difficulty
of the experimental control at high temperatures as well as the large fluctuations in the
reaction after feeding the coal. Figure 5.8a, b depict the typical cases at 700 °C and
1000 °C, respectively. These results indicate that the simulated results are reasonable
in the 600–1000 °C temperature range. Furthermore, the corresponding simulated
temperature curve has a larger gap with the measured temperature field along the
path length when the temperature difference is > 1000 °C. Thus, the theoretical
accuracy of the prediction results will be slightly reduced. Thus, the calculation and
the experimental results are relatively close to the gray-body radiation spectrum.
This also suggests that the radiation spectrum characteristics of combustion media
in the presence of particles are similar to gray-body radiation, indicating that particle
radiation has the greatest impact in solid fuel combustion [19].
The predictions of the oxyfuel conditions for cases 3 and 4 are better with small
errors up to ~ 2% for the numerical verification results of small flat-flame burner
experiment results. However, the calculation results for the air-combustion condi-
tions in cases 1 and 2 are poor, with errors of ~ 6%. Figure 5.8c, d depict typical
cases 2 and 4. The calculation results also show that the calculated radiative heat flux
and radiative source terms under oxy-fuel conditions are several times higher than
those under air-combustion conditions, indicating that the oxy-fuel combustion atmo-
sphere has stronger radiative heat transfer. The error generated when analyzing the
60 70
(a) 700°C Case Tube furnace experiment (b)
Spectral radiation proportion (%)
1000°C Case
Spectral radiation proportion (%)
0 0
0-1.125 1.125-4.1 4.1-10.5 10.5-50 0-1.125 1.125-4.1 4.1-10.5 10.5-50
Waveband (µm) Waveband (µm)
80 80
(c) Case 2 Flat-flame burner experiment (d) Case 4
Spectral radiation proportion (%)
70 70
Flat-flame burner experiment
60 60 Calculation
Calculation Experiment
50 Experiment 50
40 40
30 30
20 20
10 10
0 0
0--1.1 1.1--3 3--6 6--10 >10 0--1.1 1.1--3 3--6 6--10 >10
Waveband (µm) Waveband (µm)
Fig. 5.8 Comparison of the calculation results of the spectral radiation proportion in one-
dimensional cases a and b for tube furnace experiment and c and d for the small flat-flame burner
experiment
102 5 Radiative Energy Characteristics of Solid Fuel Combustion
radiative energy proportion may be smaller for a higher radiative heat flux. Further-
more, regardless of the oxy-fuel or air-combustion condition, both the calculation
and experimental results of the proportional distribution are similar, indicating that
particle radiation plays a significant role in the combustion medium. The results
show that a one-dimensional numerical calculation can accurately predict the spec-
tral proportion when compared to experimental results. Thus, further studies could
adopt numerical simulation to investigate spectral radiative energy characteristics in
flame as a reasonable supplementary.
5.4 Summary
This chapter introduced the spectral radiative energy and exergy characteristics in the
oxyfuel combustion of solid fuels under different atmospheres in a fixed-bed furnace
and small Hencken flat-flame burner from the perspective of energy utilization. The
effects of various fuels, temperature, atmosphere, and oxygen ratio are investigated.
The following summary can be drawn.
1. In the tube furnace experiment, the radiative energy of coal combustion increases
with oxygen ratio and temperature, with temperature having a greater influence.
The radiative exergy follows a similar pattern as the radiative energy. Temperature
influences the radiative exergy-to-energy ratio more than coal rank and oxygen
ratio. The radiative exergy-to-energy ratios of various coals are comparable at a
given combustion temperature.
The spectral radiation ratio increases in the waveband of < 4.1 μm as the oxygen
ratio and combustion temperature increase, with the temperature being the most
influential factor. At temperatures of > 800 °C, the coals have no discernible effect
on the spectral radiation distribution. The trends of radiative exergy as temperature
and oxygen ratio increase in each waveband are similar to those of the radiative
energy. Increased short-wavelength radiation ratio during combustion can improve
total radiative exergy, which is advantageous for energy conversion and management.
The radiative exergy fraction is less than the radiation fraction during combustion,
indicating that the chemical energy in the fuel is of a higher quality than that of
radiative energy. Thus, a more efficient energy cascade conversion system must
consider chemical energy conversion. Lignite has a higher radiation fraction and
radiative exergy fraction than bituminous coal. Thus, lignite or other highly volatile
solid fuels could be used for combustion in photo-thermal energy cascade conversion
systems.
other products may affect the emission and absorption of the medium along the
path length. The total radiation exergy curve along the flame height is similar to
that of radiative energy, and its absolute value is less than that of the radiative
energy. The radiation exergy-to-energy ratio curve along the path length is similar
to that of the measured temperature in different cases. Thus, the temperature is
assumed to be the primary factor affecting radiation exergy.
The spectral radiative energy has the highest proportion in the 1.1–3 μm waveband,
reaching > 60%. In different cases, the variation of spectral radiation proportion in
each waveband is within 5%. The proportion of the short-wavelength energy can
be increased by increasing the combustion temperature and oxygen ratio, whereas
the O2 /CO2 atmosphere has little effect on the proportion of spectral radiation. The
spectral radiative energy proportion can be used to approximate the spectral radiation
exergy. Short-wavelength and high-temperature radiation energy is more important
for utilization because the spectral exergy-to-energy ratio could reach between 76
and 85% in the waveband of < 3 μm.
3. The one-dimensional calculation radiation transfer cases built on combustion
characteristics can accurately simulate the proportional characteristics of spec-
tral radiation. As a reasonable supplement, further studies could use numerical
simulation to investigate spectral radiative energy characteristics in flame.
References
1. Shan S, Zhou Z, Wang Z et al (2019) Radiative energy flux characteristics and model analysis
for one-dimensional fixed-bed oxy-coal combustion. J Zhejiang Univ Sci A 20(6):431–446
2. Kirdyashkin A, Orlovskii V, Sosnin E et al (2010) Energy and spectral characteristics of
radiation during filtration combustion of natural gas. Combust Explosion Shock Waves
46(5):523–527
3. Bejarano P, Levendis Y (2008) Single-coal-particle combustion in O2 /N2 and O2 /CO2
environments. Combust Flame 153(1–2):270–287
4. Zhang Z, Li X, Luo C et al (2018) Investigation on the thermodynamic calculation of a 35
MWth oxy-fuel combustion coal-fired boiler. Int J Greenhouse Gas Control 71:36–45
5. Yuan Y, Li S, Xu Y et al (2017) Experimental and theoretical analyses on ignition and surface
temperature of dispersed coal particles in O2 /N2 and O2 /CO2 ambients. Fuel 201:93–98
6. Yuan Y, Li S, Li G et al (2014) The transition of heterogeneous–homogeneous ignitions of
dispersed coal particle streams. Combust Flame 161(9):2458–2468
7. Shao L, Zhou Z, Chen L et al (2018) Study of an improved two-colour method integrated with
the emissivity ratio model and its application to air-and oxy-fuel flames in industrial furnaces.
Measurement 123:54–61
8. Sun Y, Lou C, Zhou H (2011) A simple judgment method of gray property of flames based on
spectral analysis and the two-color method for measurements of temperatures and emissivity.
Proc Combust Inst 33(1):735–741
9. Cai X, Cheng Z, Wang S (2007) Flame measurement and combustion diagnoses with spectrum
analysis. AIP Conf Proc Am Inst Phys 914(1):60–66
10. Graeser P, Schiemann M (2020) Investigations on the emissivity of burning coal char particles:
influence of particle temperature and composition of reaction atmosphere. Fuel 263:116714
104 5 Radiative Energy Characteristics of Solid Fuel Combustion
11. Johnston J, Wooster M, Lynham T (2014) Experimental confirmation of the MWIR and LWIR
grey body assumption for vegetation fire flame emissivity. Int J Wildland Fire 23(4):463–479
12. Schiemann M, Gronarz T, Graeser P et al (2019) A correlation between char emissivity and
temperature. Fuel 256:115889
13. Modest M (2013) Radiative heat transfer, 3rd edn. Elsevier, Amsterdam
14. Moon C, Sung Y, Ahn S et al (2013) Thermochemical and combustion behaviors of coals of
different ranks and their blends for pulverized-coal combustion. Appl Thermal Eng 54(1):111–
119
15. He Y, Zhu J, Li B et al (2013) In-situ measurement of sodium and potassium release during
oxy-fuel combustion of lignite using laser-induced breakdown spectroscopy: effects of O2 and
CO2 concentration. Energy Fuels 27(2):1123–1130
16. Johansson R, Leckner B, Andersson K et al (2013) Influence of particle and gas radiation in
oxy-fuel combustion. Int J Heat Mass Transf 65:143–152
17. Gronarz T, Schulze J, Laemmerhold M et al (2017) Quantification of the influence of parameters
determining radiative heat transfer in an oxy-fuel operated boiler. Fuel Process Technol 157:76–
89
18. Ferrari C, Melino F, Pinelli M et al (2014) Thermophotovoltaic energy conversion: analytical
aspects, prototypes and experiences. Appl Energy 113:1717–1730
19. Kez V, Consalvi J, Liu F et al (2019) Investigation of gas and particle radiation modelling in
wet oxy-coal combustion atmospheres. Int J Heat Mass Transf 133:1026–1040
Chapter 6
Outlook
From the perspective of energy cascade utilization, this book points out that the
traditional power generation method is limited by the maximum parameters of the
thermal power cycle, so it cannot effectively utilize high-temperature combustion
energy. Against this, it is proposed the academic idea of photo-thermal energy cascade
conversion for fuel combustion. Several thermal radiation theories for fuel combus-
tion are introduced based on radiative energy utilization with the main innovations
of these theories are summarized as follows:
1. A radiation engine model is proposed to explain the exergy expression of black-
body radiation using the classical thermodynamic analysis method. The infinite-
staged Carnot heat engine model is innovatively created based on the concept
of equivalent temperature, and the expression of the monochromatic radiative
exergy is deduced. Consequently, a relatively complete theory for evaluating the
energy quality of spectral radiation is established. Combining the monochro-
matic radiation exergy theory with the radiation transfer equation yields the
monochromatic radiation transfer equations for radiative entropy and radiative
exergy.
2. The WSGG model for gas medium radiation characteristics is improved to
accommodate the pressurized combustion conditions. Several parameters suit-
able for the radiation characteristics of complex combustion conditions are devel-
oped based on the improved WSGG model, particularly the WSGG model
parameters in the high-pressure range and oxy-fuel combustion conditions.
Furthermore, their accuracy and necessity are verified.
3. From the perspective of energy utilization, the radiative energy characteristics
and spectral energy distribution characteristics of oxy-fuel combustion are inves-
tigated using various furnaces, with special emphasis placed on the spectral radia-
tive exergy characteristics. The findings of the study serve as the foundation for
the design and development of a photo-thermal energy cascade conversion system
for fuel combustion.
An initial academic research topic was noted on the theory and application of
thermal radiation in fuel combustion based on energy utilization, which inspires
academic ideas and scientific theories in related fields, as well as a foundation for
subsequent studies. Future studies should focus on the following areas:
1. The combustion flue gas contains CO2 , H2 O, CO, OH, and other complex gas-
phase components due to the complex combustion environment and various fuel
components. Furthermore, the products contain nitrogen oxide when ammonia
fuel combustion is considered. Consequently, relevant experiments and models
must be performed to accurately predict the radiation characteristics of flue gas
under complex combustion conditions.
2. Particle radiation is dominant in fuel combustion. Thus, the radiation character-
istics of oxyfuel combustion under high concentration particle conditions should
be investigated. The characteristics of high-temperature particle radiation should
be experimentally measured, and a particle radiation calculation model should
be developed.
3. The efficient calculation model of spectral radiation should be considered for
radiative energy utilization. A WSGG model suitable for the radiation character-
istics of multiple wavebands can be developed for gas media, and a database of
radiation characteristics of different wavebands can be established for particles.
4. The spectral range tested by the spectrometer in the research on the radiative
energy characteristics of combustion is limited due to the limitations of the
current test equipment. In the future, instruments with a wider spectrum, partic-
ularly those with radiative characteristics in the infrared waveband, should be
considered for experiments. The radiative energy characteristics of larger-scale
combustion equipment should be considered in the future because the current
experiments are on a small scale. For high-power combustion equipment (e.g.,
boilers), CFD simulation can be used to study the characteristics of combustion
radiative energy.
5. Using micro/nanostructured materials, the spectral radiative energy can be
adjusted based on the radiation characteristics of the combustion media.
Consequently, future studies should focus on the coupling relationship
between combustion reaction conditions and thermal radiation regulation
of micro/nanostructured materials. A test bench should be built especially
for the photovoltaic conversion of combustion flame radiation based on
the current measurement experiment for combustion radiation energy. The
micro/nanostructured emitter is used to reshape the flame spectrum to improve
the photovoltaic conversion efficiency based on fuel combustion.
Appendix
Table A.5 M = 1
P 1–2 (bar) 2–4 (bar)
κi K2 K1 K0 K2 K1 K0
κ1 − 0.081448825 − 0.000864595 − 0.074800399
0.023790214 0.002443296 0.00648264
κ2 − 0.366457306 0.21847772 0.000731319 0.035172062 0.482394599
0.098932084
κ3 − 2.057737171 2.944584676 − 0.973859027 3.943178351
0.410099882 0.117809229
κ4 − 42.34449532 59.57355714 − 6.92563965 94.1426706
9.89783392 0.83068445
cij C2 C1 C0 C2 C1 C0
c14 0.146388406 0.370016121 0.217756312 − 0.91945684 0.239554131
0.133781409
c13 − − − 0.473798549 − −
0.083120074 2.421659879 0.443538656 3.083571614 1.347389678
c12 − 4.521505523 0.085494128 − 3.567226694 2.406334374
0.485981366 0.589052013
c11 0.55684035 − 0.281589319 0.291446386 − −
3.031216099 1.589363583 1.540539857
c10 − 0.519086447 0.218241077 − 0.142323742 0.593255647
0.132370278 0.037742568
c24 0.21711866 − − 0.079427312 − −
0.580025964 0.454280171 0.189851541 0.683863625
c23 − 0.9647052 2.383235586 1.516783791 − 0.277513938 2.859544568
0.24753457
c22 1.529733352 − − 0.254529042 0.287916284 −
3.54582877 1.920137248 4.486810116
c21 − 2.293561238 0.881870488 − − 2.78127746
1.00260602 0.104030225 0.453293838
c20 0.226845976 − 0.192609708 0.010886816 0.128406199 −
0.539168934 0.278703918
cij C2 C1 C0 C2 C1 C0
c34 − 0.556659211 − 0.002172357 − −
0.320809314 0.147289776 0.104710652 0.116476734
c33 0.958765864 − 0.194782685 − 0.541639001 0.428506742
1.300978382 0.020973842
c32 − 0.796693448 0.246518783 0.071137741 − −
1.036536552 1.143489416 0.303812661
c31 0.507634118 − − − 0.971597734 −
0.076728817 0.576484461 0.072739954 0.351641276
c30 − 0.1158029 0.096535014 0.348535756 0.014396399 − 0.432090554
0.205640983
(continued)
Appendix 125
Table A.6 M = 2
P 1–2 (bar) 2–4 (bar)
κi K2 K1 K0 K2 K1 K0
κ1 0.004183992 − 0.068147347 0.001049754 − 0.072554919
0.010373304 0.006308614
κ2 0.101532064 − 0.650404666 − 0.06265385 0.461490718
0.24119527 0.003164009
κ3 0.28644159 − 4.391868256 − 1.04745419 3.943379962
0.021764799 0.136045831
κ4 8.0948341 − 75.08422843 − 5.649844395 67.07187336
15.57011057 0.512054615
cij C2 C1 C0 C2 C1 C0
c14 − 4.952590649 − − 1.125891931 0.112474264
1.243515114 3.230723556 0.16596521
c13 4.024870484 − 9.80373301 0.547003288 − −
16.03976646 3.494118478 1.376094169
c12 − 18.38733328 − − 3.61377363 3.128227747
4.646282902 10.31072653 0.619241648
c11 2.213411654 − 4.398299147 0.27025991 − −
8.583966607 1.344376897 2.308273297
c10 − 1.141064195 − − 0.064988686 0.767119204
0.312189098 0.251363117 0.028771924
c24 1.022651866 − 2.072909162 0.041428824 0.103572986 −
3.552667663 1.314679965
c23 − 10.26320718 − − − 5.156706239
3.034288688 5.451248065 0.082650091 0.944047166
c22 3.092021234 − 4.144747873 0.00323339 2.096348852 −
9.902423135 7.497644722
c21 − 3.74565114 − 0.055065157 − 4.523253504
1.283581888 0.781372714 1.583956059
c20 0.170701212 − 0.172432498 − 0.357434561 −
0.428002848 0.021821394 0.628351897
cij C2 C1 C0 C2 C1 C0
c34 − 0.491749416 − 0.141233332 − 1.693283319
0.169652676 0.359970633 1.156649576
c33 0.442677898 − 0.796605744 − 3.994202978 −
0.933503277 0.479159192 5.371458409
c32 − − − 0.610477647 − 6.221918426
0.280246772 0.272482378 0.053479685 5.191630272
c31 − 1.103029518 − − 3.002659825 −
0.028079824 0.868840988 0.342844492 3.409042933
c30 0.003524402 − 0.513453948 0.062865637 − 0.962470141
0.235812241 0.579002806
(continued)
Appendix 129
Table A.7 M = 4
P 1–2 2–4
κi K2 K1 K0 K2 K1 K0
κ1 − 0.0118807 0.046532356 0.018045543 0.000754349 − 0.06752071
0.003475325
κ2 0.016125356 0.065470566 0.394837816 − 0.086330646 0.447485527
0.007466612
κ3 − 5.792356922 − − 1.048119719 4.003185597
1.416839776 0.383808258 0.141469638
κ4 − 100.6233201 − − 9.80541848 52.42165954
26.39996476 27.1921835 0.89447472
cij C2 C1 C0 C2 C1 C0
c14 − 1.899480863 − − 1.497825229 −
0.361394306 0.579049571 0.220416154 0.339650914
c13 1.292774778 − 1.54464316 0.701108077 − −
6.574744127 4.520110356 0.197957579
c12 − 7.991516215 − − 4.53672888 2.177164989
1.650468142 1.18528327 0.76368654
c11 0.879092358 − 0.293455506 0.321411648 − −
3.947963937 1.647388273 2.076972979
c10 − 0.494250895 0.321309732 − 0.103367952 0.747901899
0.124607614 0.035814185
c24 − 0.281809353 − 0.103422307 − −
0.106869686 1.101283986 0.327830563 0.723172126
c23 0.275276188 − 3.744270311 − 0.316857359 3.354365269
0.954798952 0.263075707
c22 − 1.271623902 − 0.17996892 0.856379971 −
0.208410976 4.888539825 5.611571547
c21 0.01913329 − 2.749736151 − − 3.775703659
0.671570019 0.011633992 1.123019209
c20 0.006765944 0.125801296 − − 0.296651647 −
0.269812754 0.012641271 0.533884597
cij C2 C1 C0 C2 C1 C0
c34 1.212447928 − 3.59044676 0.114052604 − 1.001655337
4.276503486 0.785317125
c33 − 13.94028924 − − 2.885346145 −
3.896163398 11.31184997 0.408054214 3.15440051
c32 4.371944412 − 12.47918026 0.554057703 − 3.660562599
16.10193081 4.05684856
c31 − 7.533934773 − − 2.555029124 −
1.930912478 5.708099259 0.331516476 2.147871972
c30 0.2068702 − 0.985055123 0.065144754 − 0.765549326
0.928642764 0.535438972
(continued)
Appendix 133