Research Article  Open Access
Extension of Weighted Sum of Gray Gas Data to Mathematical Simulation of Radiative Heat Transfer in a Boiler with GasSoot Media
Abstract
In this study an expression for soot absorption coefficient is introduced to extend the weightedsumofgray gases data to the furnace medium containing gassoot mixture in a utility boiler 150 MWe. Heat transfer and temperature distribution of walls and within the furnace space are predicted by zone method technique. Analyses have been done considering both cases of presence and absence of soot particles at 100% load. To validate the proposed soot absorption coefficient, the expression is coupled with the Taylor and Foster's data as well as Truelove's data for CO_{2}H_{2}O mixture and the total emissivities are calculated and compared with the Truelove's parameters for 3term and 4term gray gases plus two soot absorption coefficients. In addition, some experiments were conducted at 100% and 75% loads to measure furnace exit gas temperature as well as the rate of steam production. The predicted results show good agreement with the measured data at the power plant site.
1. Introduction
One of the most important modes of heat transfer in a boiler furnace of a large power plant is radiation. Thus determination of the radiative properties of combustion products is vital to predict the temperature and heat flux distribution. The main products of combustion in an enclosure contain carbon dioxide and water vapor [1]. Some studies have been done on modeling of gas mixture in furnaces. Different models for calculation of the radiative properties of real gases have been proposed by researchers such as statistical narrow band (SNB), the fullspectrum correlatedk (FSCK) distribution, and the weightedsumofgray gases (WSGG) [2–6]. Among these models, the last one is more reasonable in engineering calculations in view of the accuracy and computing time. The WSGG model was developed by Hottel and Sarofim [1]. Modest [7] stated that this model can be applied for any solution methods for the transport equation. Based on WSGG concept, Taylor and Foster [8] utilized a “three gray plus one clear gas” model. A threeterm mixed gray gas model with third order polynomial for weighting factors was employed by Smith et al. [9]. Soufiani and Djavdan [10] proposed a fiveorder polynomial for gas combustion, where . After that a lot of efforts were given to model the radiative properties from semitransparent media containing a mixture of nongray gases and soot [11, 12] which was applied for analyzing the radiative heat transfer inside the furnace [13]. Among the different methods of modeling the combusting environments such as inverse methodologies, hybrid method, discreteordinates method, and finite volume method [14–17], the zone method is more practical and many attempts on radiative heat transfer analysis were conducted based on this method [18–20]. Zone method was originally developed by Hottel and Cohen [21] for an absorbing, emitting, nonscattering gray gas with constant absorption coefficient. Later Hottel and Sarofim [1] extended it to deal with threedimensional problems. Also Larsen and Howell [22] presented a method for calculations related to the direct exchange areas in zonal analysis based on lastsquares smoothing. Tucker [23] conducted a numerical integration and suggested an exponential expression for exchange areas which covers a range of optical thickness from 0 to 18. Lawson [24] proposed an improved method for smoothing approximate exchange areas. To achieve the total exchange areas, Noble [25] presented the explicit matrix relations. Batu and Selçuk [26] analyzed the radiative heat transfer in the freeboard of a fluidized bed combustor by using the zone method. Bordbar and Hyppanen [27] employed the zone method for predicting temperature and heat flux on the water walls of a steam boiler furnace. Recently, Méchi and coworkers [19] proposed a radiative model to extend the zonal method to semitransparent inhomogeneous composed of nongray gas and soot. Also Moghari et al. [18] used the zone method to predict thermal radiation behavior in the Dtype watercooled steam boiler furnace. Crnomarkovic et al. [28] used the simple gray gas (SGG) and WSGG to model the radiative properties of the twophase mixture composed of gas and particles inside the lignite fired furnace.
In this study a new expression for soot absorption coefficient has been presented depending on temperature which could be coupled with nonluminous flame data containing several gray gases and one clear gas. The results are based on the suggested soot absorption coefficient coupled with the data generated by Taylor and Foster. The validity of the calculated soot absorption coefficient is confirmed by comparison with the obtained total emissivities and the calculated values from Truelove’s models. In addition, for reconfirmation of the results, the soot expression is utilized in the zone method to model the furnace of a utility boiler 150 MWe. The temperature and heat flux distributions are discussed for 2 cases (with and without soot particles) at 100% load. Furthermore the furnace exit gas temperature and amount of steam production by considering the effect of soot for the loads of 100%; 75% are presented and compared with the captured data from the site.
2. Mathematical Model
2.1. The Weighted Sum of Gray Gas
The weightedsumofgray gases (WSGG) is one of the accurate techniques for modeling the radiative behavior of combustion gases. The total emissivity of real gas can be represented mathematically by a mixture of gray gases [1]: where, , and represent the absorption coefficient for the th gray gas, sum of the partial pressure of all radiating gases in the mixture, and effective path length, respectively, andis weighting factors [5, 9, 10, 12] of various commonly used correlations for the mixture of combustion products which have been reported by Taylor and Foster [8], Smith et al. [9], and Soufiani and Djavdan [10].
In fact WSGG is an appropriate tool which could be applied in the modeling of media containing CO_{2}, H_{2}O, and soot and, in this subject, some approaches have been developed to consider the effect of soot particles [12, 29].
Based on the suggestion of Truelove [12] for the gassoot mixture, the two absorption coefficients (gas mixture and soot) are contributed in the calculations. The expressions for emissivity of the combustion productsoot mixturecan be presented by: where is the soot concentration.
In order to determine the soot absorption coefficient, a relationship from the wavelength dependence of derived from experimental investigations is [30] as follows: where and .
By integrating over wavelength we have Introducinginto above equation, the soot absorption coefficient becomes where = 3.742 × 10 W·μm^{4}/m^{2} and = 1.4388 × 10^{4} μm·K are the first and second Planck function constants, respectively, = 5.669 × 10^{8} W m^{−2 }K^{−4} is StephaneBoltzman constant, is Gamma function, and is Rieman zeta function. The above temperature dependence as determined from (5) is expressed by following the simple polynomial equation: where is the temperature of the radiation source in Kelvin.
2.2. Zonal Method
In zone method the enclosure is subdivided into surfaces and volumes zones which could be assumed isothermal [31]. Then by using the gas flow and combustion pattern, the mass flow rate from/to each volume zone, generated heat by combustion and convection coefficients are obtained. A steady state energy balance is considered for each zone and then a set of simultaneous equations based on the temperatures and heat fluxes are produced. By solving these equations, the temperature and heat flux distributions are obtained.
Calculation of the Direct Exchange Areas and Total Exchange Areas. For finding the radiative heat transfer between two zones, the first step is to calculate the direct exchange areas (DEA) and then the total exchange areas (TEA). There are three types of DEAs: surfacesurface, volumesurface, and volumevolume. For instant, the volumesurface direct exchange area as illustrated in Figure 1 can be determined as follows:
The DEAs obey the reciprocity definitions where and . Direct numerical integration can be applied to calculate the respective areas.
For the gray gas the total flux between two zones and must be proportional to and the proportionality constant, called the total exchange area, is indicated by [1]. All of these terms are calculated by using the methods that have been reported by Hottel and Sarofim [1] and Modest [7].
Direct Flux Areas. The radiant energy between any two zones is proportional to the aweighted summation of the total exchange areas for each gas. For example, the net flux between zones and is given by [1]: whereandare replacing the terms in the brackets. These are called directedflux areas [1, 29]. Similarly expression for surfacesurface transfer is And for gasgas transfer it is expressed by
Total Energy Balance. For a volume zone , the total energy balance can be stated by whereand are the number of volume and surface zones, respectively. is the number of gases in the model, is the convection heat transfer to all surfaces in contact with the volume zone, andis the total sensible heat presented by whereis the mass flow rate of gas entering the zone from a neighboring zone , andis representing the total mass flow rate of gas leaving the zone . Also is heat released due to combustion plus the heat content in the combustion air, so this term can be expressed [27]:
On the other hand for a surface zone , the total energy balance could be represented by where is heat transfer rate to water walls.
Finally, the energy balance for the total number of volume and surface zones generates a series of nonlinear algebraic equations. These equations should be solved by the iterative techniques in order to achieve the temperature distribution in zones. In this study, the surfaces have been assumed gray and the combustion is complete in the zones in front of burners.
3. Experimental Facility
The experimental data were obtained from the furnace of a 150 MWe utility boiler. Schematic of boiler is shown in Figure 2. The dimensions of the boiler furnace arewhich is equipped with 9 natural gas fired burners in three rows of three. They are located at the left side wall of the furnace chamber. The operating conditions of the boiler and fuel characteristics are mentioned in Table 1. The experiments were conducted to measure the furnace exit gas temperature by a thermocouple with reasonable accuracy (0.05% of reading) located adequately far from the last raw of burners at the furnace outlet at 100% and 75% loads.

4. Results and Discussion
The soot absorption coefficient suggested in (6) can be coupled with available models for nonluminous flam with gray gases and one clear gas to obtain the total emissivities of gassoot mixture and to be applied in the zone method.
Validation of Presented Soot Absorption Coefficient with Truelove’s Model. The total emissivities of gassoot mixture are obtained by coupling the calculated by using (6) and Taylor’s data. The evaluated total emissivities at the different temperatures with soot concentration of 0.0001 Kg/m^{3} for gas combustion are well compared against benchmark data which are “threegray gas plus twosoot” and “fourgray gas plus twosoot” models suggested by Truelove (Figure 3). Also to show the suitability of using the present expression in other models, the suggestedis coupled with Truelove’s models without soot which are threeterm (twogray plus one clear) and fourterm (threegray plus one clear) gas models and the results are presented in Figure 3 as well.
Table 2 presents the discrepancies between the computed total emissivities by using the coupled models and benchmarks.

Figure 4 demonstrates the comparison between the calculated gassoot mixture total emissivity and benchmark data. The total emissivities are plotted versus path length on a logarithmic axis with two different soot concentrations (0.0001 Kg/m^{3} and 0.005 Kg/m^{3}) for gas, at 800, 1600, and 2400°K.
(a)
(b)
(c)
It is seen that the obtained results coincide with the benchmarks; deviations are acceptable specially in comparison of two cases of models: (i) present model coupled with Taylor’s data and (ii) present model coupled with 3term Truelove’s data with the Truelove’s model (threegray gas plus 2soot). The errors are not greater than 4.8 percent at > 0.1 m for the first case and 2 percent at all path lengths for second case for small soot concentration (0.0001 Kg/m^{3}) at 1600°K.
Application in the Zone Method. The proposed soot absorption coefficient coupled with Taylor’s data is used for radiative heat transfer analysis inside the boiler furnace. The emissivity of tube wall is considered 0.8. The zones of furnace are obtained by dividing the height ( direction) into five equal sections, the length ( direction) into 2 equal sections, and the width ( direction) into 2 equal sections; thus the furnace has been divided into 46 surface and 18 volume zones, as shown in Figure 5.
The results are based on the effects of existence and absence of soot particles. Tables 3 and 4 show the temperature distribution on height direction of front wall at 100% load for case 1 (excluding soot particles) and case 2 (including soot particles). It is noticeable that in the results of case 2, the value is fixed at an average of zones temperature which is obtained from results by excluding soot particles. The amount of excess air is 5% and soot concentration is 0.00005 Kg/m^{3}.


The heat flux distribution along walls (front wall, right side wall, and left side wall) for both cases with soot and without soot is shown in Figure 6. It is observed that existence of soot raised considerably the heat flux on the wall. Similar result was reported in the literature [19, 32–34]. In fact the existence of soot enhances the radiative intensity because of continuum radiation in the visible and infrared regions of the wavelength spectrum [35] and in this situation the radiative transfer is conducted towards the wall and as a result of this, temperature of the medium has reduced.
(a)
(b)
(c)
(d)
To confirm the obtained results and validate the applied mathematical model, Tables 5 and 6 present a comparison between the calculated data and measured practical data. Available results of measurements are furnace exit gas temperature and rate of steam production at 100% and 75% loads.


The table shows that there is a good agreement between the present results and the experimental data.
5. Conclusions
A new soot absorption coefficient proposed in this paper has been assessed though coupling with WSGG parameters suggested by Taylor. It has been utilized for modeling the radiation heat transfer in a utility boiler 150 MWe. The total emissivities are calculated and compared with the Truelove’s parameters for 3term and 4term gray gases plus twosoot absorption coefficients. In addition, some experiments were conducted at 100% and 75% loads to measure furnace exit gas temperature as well as the rate of steam production and the following results are obtained.(1)The soot absorption coefficient model is compatible with WSGG models containing gray gases and one clear gas.(2)The existence of soot particles leads to a decrease in gas temperature and an increase in wall heat flux.(3)The exhaust gas temperature and steam production could be estimated with reasonable accuracy at different loads.
Nomenclature
:  Surface area (m^{2}) 
:  Weighting factor 
:  Mass concentration (Kg/m^{3}) 
:  First Planck function constant (W μm^{4}/m^{2}) 
:  Second Planck function constant (μm K) 
:  Specific heat capacity at constant pressure (kJ/kg K) 
:  Black body emissive power (W/m^{2}) 
Direct exchange area, total exchange area, and flux exchange area for volume to volume , (m^{2})  
:  Direct exchange area, total exchange area, and flux exchange area for volume to surface , (m^{2}) 
:  Extinction coefficient of the medium (m^{−1}) 
:  Soot absorption coefficient (m^{2}/Kg) 
:  Effective path length (m) 
:  Number of gases 
:  Partial pressure of gas (atm) 
:  Heat transfer rate (W) 
:  Distance between two zones (m) 
:  Stoichiometric air/fuel volume ratio 
Direct exchange area, total exchange area, and flux exchange area for surface to surface , (m^{2})  
:  Temperature (K) 
:  Volume (m^{3}) 
:  Percentage excess air level. 
:  Gas absorptivity 
Gas emissivity  
:  Emissivity of gassoot mixture 
:  Angle between the beam joining the center points of two zones and the normal to one of the two zones (Rad) 
:  Gamma function 
:  Rieman zeta function 
:  Stephane Boltzman constant (W/m^{2} K^{4}) 
:  Density (Kg/m^{3}). 
Air  
:  Gas 
:  Surface or volume zone 
:  Number of volume 
:  Number of surface 
:  th gray gas 
:  Surface or soot. 
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The authors gratefully acknowledge High Impact Research Grant UM.C/HIR/MOHE/ENG/45 and UMRG Grant RP012D13AET, University of Malaya, Malaysia, for support to conduct this research work. The authors are also grateful to the persons in charge in power plant site for their technical support during the experiments.
References
 H. C. Hottel and A. F. Sarofim, Radiative Transfer, McGrawHill, 1967.
 M. N. Borjini, K. Guedri, and R. Saïd, “Modeling of radiative heat transfer in 3D complex boiler with nongray sooting media,” Journal of Quantitative Spectroscopy and Radiative Transfer, vol. 105, no. 2, pp. 167–179, 2007. View at: Publisher Site  Google Scholar
 V. Goutiere, F. Liu, and A. Charette, “An assessment of realgas modelling in 2D enclosures,” Journal of Quantitative Spectroscopy and Radiative Transfer, vol. 64, no. 3, pp. 299–326, 1999. View at: Google Scholar
 P. J. Coelho, “Numerical simulation of radiative heat transfer from nongray gases in threedimensional enclosures,” Journal of Quantitative Spectroscopy and Radiative Transfer, vol. 74, no. 3, pp. 307–328, 2002. View at: Publisher Site  Google Scholar
 I. H. Farag, “Nonluminous gas radiation: approximate emissivity models,” in Proceedings of the 7th International Heat Transfer Conference, vol. 2, pp. 487–492, Miinchen, Germany, 1982. View at: Google Scholar
 H. Q. Chu, F. S. Liu, and H. C. Zhou, “Calculations of gas radiation heat transfer in a twodimensional rectangular enclosure using the linebyline approach and the statistical narrowband correlatedk model,” International Journal of Thermal Sciences, vol. 59, pp. 66–74, 2012. View at: Google Scholar
 M. F. Modest, Radiative Heat Transfer, Academic Press, 2003.
 P. B. Taylor and P. J. Foster, “The total emissivities of luminous and nonluminous flames,” International Journal of Heat and Mass Transfer, vol. 17, no. 12, pp. 1591–1605, 1974. View at: Google Scholar
 T. F. Smith, Z. F. Shen, and J. N. Friedman, “Evaluation of coefficients for the weighted sum of gray gases model,” Journal of Heat TransferTransactions of the ASME, vol. V 104, no. 4, pp. 602–608, 1982. View at: Google Scholar
 A. Soufiani and E. Djavdan, “A comparison between weighted sum of gray gases and statistical narrow band radiation models for combustion applications,” Combustion and Flame, vol. 97, no. 2, pp. 240–250, 1994. View at: Publisher Site  Google Scholar
 A. T. Modak, “Radiation from products of combustion,” Fire Safety Journal, vol. 1, no. 6, pp. 339–361, 1979. View at: Google Scholar
 J. S. Truelove, A Mixed Grey Gas Model for Flame Radiation, Thermodynamics Division, AERE, 1976.
 R. Yadav, A. Kushari, A. K. Verma, and V. Eswaran, “Weighted sum of gray gas modeling for nongray radiation in combusting environment using the hybrid solution methodology,” Numerical Heat Transfer, Part B, vol. 64, no. 2, pp. 174–197, 2013. View at: Google Scholar
 S. Payan, S. M. H. Sarvari, and A. Behzadmehr, “Inverse estimation of temperature profile in a nongray medium with soot particles between two parallel plates,” Numerical Heat Transfer, Part A, vol. 63, no. 1, pp. 31–54, 2013. View at: Google Scholar
 F. C. Lockwood and N. G. Shah, “A new radiation solution method for incorporation in general combustion prediction procedures,” Symposium (International) on Combustion, vol. 18, no. 1, pp. 1405–1414, 1981. View at: Publisher Site  Google Scholar
 A. S. Jamaluddin and P. J. Smith, “Predicting radiativetransfer in rectangular enclosures using the discrete ordinates method,” Combustion Science and Technology, vol. 59, no. 4–6, pp. 321–340, 1988. View at: Google Scholar
 P. J. Coelho, “Fundamentals of a new method for the solution of the radiative transfer equation,” International Journal of Thermal Sciences, vol. 44, no. 9, pp. 809–821, 2005. View at: Publisher Site  Google Scholar
 M. Moghari, S. Hosseini, H. Shokouhmand, H. Sharifi, and S. Izadpanah, “A numerical study on thermal behavior of a Dtype watercooled steam boiler,” Applied Thermal Engineering, vol. 37, pp. 360–372, 2012. View at: Publisher Site  Google Scholar
 R. Méchi, H. Farhat, K. Guedri, K. Halouani, and R. Said, “Extension of the zonal method to inhomogeneous nongrey semitransparent medium,” Energy, vol. 35, no. 1, pp. 1–15, 2010. View at: Publisher Site  Google Scholar
 H. Ebrahimi, A. Zamaniyan, J. S. S. Mohammadzadeh, and A. A. Khalili, “Zonal modeling of radiative heat transfer in industrial furnaces using simplified model for exchange area calculation,” Applied Mathematical Modelling, vol. 37, no. 1617, pp. 8004–8015, 2013. View at: Google Scholar
 H. C. Hottel and E. S. Cohen, “Radiant heat exchange in a gasfilled enclosure: allowance for nonuniformity of gas temperature,” AIChE Journal, vol. 4, no. 1, pp. 3–14, 1958. View at: Google Scholar
 M. E. Larsen and J. R. Howell, “Leastsquares smoothing of directexchange areas in zonal analysis,” Journal of Heat Transfer, vol. 18, pp. 239–242, 1986. View at: Google Scholar
 R. J. Tucker, “Direct exchange areas for calculating radiation transfer in rectangular furnaces,” Journal of Heat Transfer, vol. 108, p. 707, 1986. View at: Google Scholar
 D. A. Lawson, “An improved method for smoothing approximate exchange areas,” International Journal of Heat and Mass Transfer, vol. 38, no. 16, pp. 3109–3110, 1995. View at: Google Scholar
 J. J. Noble, “The zone method: explicit matrix relations for total exchange areas,” International Journal of Heat and Mass Transfer, vol. 18, no. 2, pp. 261–269, 1975. View at: Google Scholar
 A. Batu and N. Selçuk, “Modeling of radiative heat transfer in the freeboard of a fluidized bed combustor using the zone method of analysis,” Turkish Journal of Engineering and Environmental Sciences, vol. 26, no. 1, pp. 49–58, 2002. View at: Google Scholar
 M. H. Bordbar and T. Hyppanen, “Modeling of radiation heat transfer in a boiler furnace,” Advanced Studies in Theoretical Physics, vol. 1, no. 12, pp. 571–584, 2007. View at: Google Scholar
 N. Crnomarkovic, M. Sijercic, S. Belosevic, D. Tucakovic, and T. Zivanovic, “Numerical investigation of processes in the lignitefired furnace when simple gray gas and weighted sum of gray gases models are used,” International Journal of Heat and Mass Transfer, vol. 56, no. 1, pp. 197–205, 2013. View at: Google Scholar
 T. R. Johnson and J. M. Beer, “Radiative heat transfer in furnaces: further development of the zone method of analysis,” Symposium (International) on Combustion, vol. 14, no. 1, pp. 639–649, 1973. View at: Publisher Site  Google Scholar
 J. M. Rhine and R. J. Tucker, Modelling of GasFired Furnaces and Boilers and Other Industrial Heating Processes,, British Gas, 1991.
 R. Siegel and J. R. Howell, Thermal Radiation Heat Transfer, Taylor & Francis, 2002.
 F. Liu, H. A. Becker, and Y. Bindar, “A comparative study of radiative heat transfer modelling in gasfired furnaces using the simple grey gas and the weightedsumofgreygases models,” International Journal of Heat and Mass Transfer, vol. 41, no. 22, pp. 3357–3371, 1998. View at: Google Scholar
 N. W. Bressloff, “The influence of soot loading on weighted sum of grey gases solutions to the radiative transfer equation across mixtures of gases and soot,” International Journal of Heat and Mass Transfer, vol. 42, no. 18, pp. 3469–3480, 1999. View at: Publisher Site  Google Scholar
 Y.L. Hwang and J. R. Howell, “Local furnace data and modeling comparison for a 600MWe coalfired utility boiler,” Journal of Energy Resources Technology, vol. 124, no. 1, pp. 56–66, 2002. View at: Publisher Site  Google Scholar
 T. L. Farias, M. G. Carvalho, and Ü. Ö. Köylü, “Radiative heat transfer in sootcontaining combustion systems with aggregation,” International Journal of Heat and Mass Transfer, vol. 41, no. 17, pp. 2581–2587, 1998. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2014 Samira Gharehkhani et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.