Influence of Reactant Mixing in a Laminar Flow Reactor: The Case of

A theoretical study of the influence of mixing on reburning under laminar flow non-premixed conditions has been performed. Two models with a different...
0 downloads 0 Views 248KB Size
3528

Ind. Eng. Chem. Res. 2007, 46, 3528-3537

Influence of Reactant Mixing in a Laminar Flow Reactor: The Case of Gas Reburning. 2. Modelling Study F. Liesa, M. U. Alzueta,* A. Millera, and R. Bilbao Arago´ n Institute of Engineering Research, Department of Chemical and EnVironmental Engineering, UniVersity of Zaragoza, Campus Rı´o Ebro, 50018 Zaragoza, Spain

A theoretical study of the influence of mixing on reburning under laminar flow non-premixed conditions has been performed. Two models with a different level of detail of the fluid dynamics description have been developed and adapted. Models have been validated against experimental results. The models are able to reproduce the experimental trends of the influence of the studied variables: temperature, stoichiometry, velocity ratio between the streams, and reactant feeding configuration. Modeling and experimental results have been compared with plug-flow reactor (PFR) results, showing significant differences. Reactant injection mode involves a change in local stoichiometry and residence time distribution in the reactor, which seems to be the main influencing factor. For low temperatures, a higher NO reduction is obtained when local fuel-lean conditions have a longer residence time. For high temperatures, a higher NO reduction is obtained when local fuel-rich conditions have a longer residence time. Introduction Reburning is a NOx control technology that involves the injection of a secondary fuel as a reducing agent (i.e., reburning fuel) into the so-called “reburning zone” downstream of the main combustion chamber of the furnace. The hydrocarbon radicals generated in the reburning zone provide the chemical pathways to reduce NOx to N2. The application of reburning process as a NOx control technology has been successfully demonstrated at full scale; therefore, it is a valuable technology to be further studied and optimized. In this process, reactants are not premixed and, in most cases, mixing and reaction times are of the same order of magnitude.1 Therefore, reburning is a process influenced by mixing effects. In fact, several authors2-8 have suggested that mixing between primary combustion products and reburning fuel is one of the main factors that influence NOx reduction performance through reburning. Thus, despite the valuable results obtained by means of reactors in which mixing effects are not considered, the performance of experiments to study the influence of mixing and the development of models to simulate mixing conditions are required. Most of the developed models include the influence of mixing by means of one of these two different approximations. Reburning fuel is entrained into the bulk flow coming from primary combustion, direct mixing,3,9,10 or the opposite situation, inverse mixing1,7-9,11,12 (i.e., the main flow is entrained into the reburning fuel jet). Both of these considerations improve the description of pilot- and real-scale reburn performance, compared to models that do not take into account mixing effects. Among these considerations, inverse mixing, which is a morerealistic description of physical processes, gives modeling results more approximated to experimental data.9,12 Despite the simplicity of this approximation, it takes into account the segregation of the reactants in the mixing area of the reactor, which may be the reason for its success.7 * To whom correspondence should be addressed. Tel.: +34 976 761876. Fax: +34 976 761879. E-mail: [email protected].

To describe the entrainment of the bulk flow into the reburning fuel jet, several approximations have been used, including instantaneous mixing, but with a distributed addition of the reactants in the reaction zone during a fixed period of time1,3,7,9-11 (mixing time), and the application of the two-stage Lagrangian (TSL) model.8,12 Using this model, the mixing time can be determined within the model itself. Moreover, it includes the reaction process, which occurs in an additional reaction zone (the flame sheet). The flame sheet temperature and composition can be different from the main reactor temperature and composition and, therefore, the reaction rate can also be different. When using a distributed addition of the reactants, mixing time is determined either empirically, or by means of the use of detailed fluid flow field calculation (e.g., by means of CFD calculations, or entrainment relations). Experimental work under pilot- and full-scale conditions have contributed to obtain interesting conclusions about reburning process. However, in these experimental facilities, the control and the study of how each variable (i.e., temperature, stoichiometry, reburning fuel, residence time, residence time distribution) influence the entire process, independent of the other variables, is quite complex. Moreover, the manner in which mixing occurs affects the local conditions in the reaction zones mainly, temperature, local stoichiometry, and residence time distributionsand affect the entire process. Because of the complexity of these reactive systems, a detailed description of both fluid dynamics and reaction processes is too computationally costly and most of the time, it is not feasible. Thus, the realization of a study of reburning process at laboratory scale, under well-controlled conditions to get a better understanding on how mixing affects the process, was considered interesting. Therefore, an experimental and modeling study of the influence of mixing under turbulent conditions was performed.13 In that study, it was stated that mixing has an effect on the local stoichiometry, which influences NOx reduction. Moreover, it was shown that simplified mixing models, as an inverse mixing model with a distributed addition of the reactants, did not reproduce the experimental behavior.

10.1021/ie060944i CCC: $37.00 © 2007 American Chemical Society Published on Web 04/10/2007

Ind. Eng. Chem. Res., Vol. 46, No. 11, 2007 3529

Mixing times (35-55 ms) in the experiments under turbulent conditions were shorter than in several pilot-scale7,14,15 and fullscale facilities. That is the reason why it was thought that the influence of mixing could be visualized much better by working under laminar flow non-premixed conditions, because the mixing times under these conditions are longer. Therefore, an experimental study of the influence of the different variables under non-premixed laminar conditions was performed.16 In this work, additionally, the result tendencies obtained from the experimental installation were compared with others obtained under plug-flow reactor (PFR) conditions. To analyze the results of those experiments,16 several models have been used and adapted. Similarly, these models allow us to study the influence of the mixing process between the reburning fuel and the primary combustion gases in the NO reduction through reburning from a theoretical viewpoint. The models used here do not require adjustable parameters. Modelling of Reburning Process under Laminar Mixing Conditions The theoretical description of the reburning process requires modeling of the chemical processes and mixing conditions that are occurring in the reactive system. Chemical processes are modeled by means of reaction mechanisms, and mixing conditions are determined by the fluid dynamics in the reactor as well as the conditions of the reactants at the entrance (premixed or non-premixed). Regarding the chemical kinetics of the reburning process, mechanisms with a different level of complexity can be found in the bibliography. The level of detail in the kinetic description of the reburning process determines the existence of reaction mechanisms that are classified as detailed mechanisms,17-21 simplified mechanisms,22-24 and reduced mechanisms.22,25,26 The different level of complexity of the mechanisms determines not only the accuracy of the results, but also their applicability to the modeling of systems, which are, more or less, complex, from the fluid-dynamics viewpoint. Thus, detailed mechanisms have a very wide application range, with regard to operational conditions (mainly, temperature and stoichiometry). Nevertheless, the computational cost requirements in complex systems are high and this significantly limits the use of the detailed schemes. Reduced mechanisms have a greater applicability, because of their simplicity and the low computational cost, compared with detailed mechanisms. However, in most cases, their application is restricted to a reduced interval of operational conditions. Because of computational advances, the reduced mechanisms currently developed27,28 incorporate a greater number of reactions and chemical species, because their application to fluid-dynamics complex systems without an excessive computational time is possible. These mechanisms, the so-called “augmented reduced mechanisms”, have an interval of validity of the variables that is higher, compared to the reduced mechanisms, which represents a significant advantage. In this work, the Augmented Reduced Mechanism (ARM) of Giral and Alzueta28 has been used to describe the kinetics of the process. The results obtained with this model for a wide interval of operation conditions were compared with those obtained from the detailed mechanism of Glarborg et al.20 and similar tendencies and results were obtained. The complexity of the fluid dynamics of the reactive system allows the use of this model with a reasonable computational cost, which makes possible to get trustworthy results, from the kinetic viewpoint. The ARM28 consists of 15 lumped reactions and 19 chemical species, whereas the initial mechanism20 consists of 438

Figure 1. Detail of the laminar flow reactor. Dimensions given in millimeters.

reactions and 65 chemical species. It is noteworthy to note that the used ARM only describes NO reduction through reburn reactions, because a parallel experimental work16 showed that this was the only significant mechanism for NO decrease, under the conditions studied. The fluid dynamics of the process is determined by the characteristics of the reactor and the operational conditions of the system. The reactor used in the experiments16 is presented in Figure 1, together with its dimensions. It includes two separate co-flowing streams (labeled as “1” and “2” in Figure 1) that are fed into the reaction zone coaxially. Both streams flow without mixing through the heating zone (labeled as “4” in Figure 1). Therefore, the reactor operates in non-premixed conditions. After the streams penetrate in the feature labeled as “5” in Figure 1 they start to mix to each other, mainly by diffusion, and react. In this reactor, to be different from ideal

3530

Ind. Eng. Chem. Res., Vol. 46, No. 11, 2007

Figure 2. Physical scheme and boundaries representation of the modeled system. Table 1. Experiments of Liesa et al.16 set

inner tube

O2 (ppm)

λ

total flow rate (NmL/min)

tra (s)

R

1 2 3 4 5 6 7 8

CH4,C2H6 NO, O2 CH4,C2H6 CH4,C2H6 CH4,C2H6 CH4,C2H6 NO, O2 NO, O2

1174 1136 1080 1925 3789 1081 1140 1151

0.60 0.60 0.57 1.00 1.97 0.57 0.59 0.59

6986 7073 4622 4616 4618 4567 7080 6966

814/T[K] 804/T[K] 1230/T[K] 1231/T[K] 1231/T[K] 1245/T[K] 803/T[K] 816/T[K]

8 2.8 1 1 1 8.3 1 7.9

a The average residence time (t ) is dependent on the temperature, as r listed. The temperature range of study is 900-1500 K.

laboratory reactors, the mixing time is not short enough to be considered negligible, compared with the reaction time; thus, both mixing and reaction occur simultaneously and must be taken into account to describe the process theoretically. The experiments were performed under highly diluted conditions, using nitrogen as a carrier gas. The streams fed through the two different inlets consist of H2O, O2, and NO, to be representative of the products exiting the primary combustion zone of interest in a reburn configuration for one of the streams, and the reburn fuel, which is a synthetic mixture prepared with ∼90% CH4 and ∼10% C2H6, to be representative of a natural gas composition for the second stream. The global concentration of hydrocarbons diluted in N2 has been 800 ( 25 ppm of CH4 and 95 ( 5 ppm of C2H6 and the stoichiometry values (λ) used were in the range of 0.57-1.98. The temperatures used in the experiments were in the range of 1000-1500 K. The maximum total flow rate used in the experiments was ∼7000 NmL/min and the minimum was 2600 NmL/min that implies average residence times of 0.5-2 s. The specific experimental conditions considered in the present work are summarized in Table 1. Because of the experimental conditions used, the following approximations can be adopted: (1) The system can be considered axisymmetric, and, therefore, bidimensional. Taking into account this fact, the modeled system and its boundaries, where the boundary conditions are applied, is presented in Figure 2. (2) The reactor works under stationary conditions. (3) The reactor works under laminar flow conditions, because the Reynolds number is 1325 K, a greater NO reduction is obtained with the BLM and LMM models than with the PFR Model. For fuellean conditions (λ ) 2), this behavior is shifted to lower temperatures (∼1250 K); above this temperature, a greater reduction of NO is obtained with the LMM and BLM models, compared to the PFR model. This behavior may be attributed to the effect of non-premixed conditions that are modeled by LMM and BLM and not in PFR. Non-premixed conditions lead to the existence of fuel-rich zones within the reactor, despite the globally fuel-lean conditions. These fuel-rich zones, which do not exist under PFR conditions, favor NO reduction. Under fuel-lean and PFR conditions, the hydrocarbons oxidation process is too fast and does not allow hydrocarbon radicals to interact with NO effectively. Besides the differences previously explained, there is another difference between LMM, BLM, and PFR predictions for all the values of λ studied. There is a maximum of NO reduction with temperature that has been predicted by the PFR model that is not observed in any of the following cases: BLM model, LMM model, and experiments. Figures 5 and 7 show that the LMM and BLM models describe the experimental behavior better than the PFR model for fuel-rich, stoichiometric, and fuel-lean conditions in a wide range of temperatures. Moreover, results obtained with the LMM and BLM models are very similar to each other. As an example of the influence of stoichiometry under different mixing conditions, in Figure 8, there is a comparison of NO reduction modeling results obtained at 1425 K with PFR and LMM models for different λ values. Under PFR conditions, the evolution in the NO reduction, as the value of λ increases, is sharper than under non-premixed laminar conditions (LMM).

It can be observed that under non-premixed conditions and laminar mixing (LMM model), the maximum of NO reduction is shifted to leaner conditions, compared to PFR. Moreover, for λ > 0.7, a higher reduction of NO is obtained under laminar mixing conditions. This is a probe of evidence that a delayed mixing may result in an improvement of reburning performance. Miller et al.43 mentioned similar observations in a pilot plant, in which they obtained a higher reburn performance in a lean atmosphere while optimum reburn conditions are usually considered to be fuel-rich. They already noted that the explanation for those results could rely on the mixing effects. Moreover, Glickert and Houy44 quantified significant NO reductions under fuel-lean reburning conditions in a real installation. Despite the globally fuel-lean conditions in all cases, because of nonpremixed reactant feeding, there are fuel-rich zones within the reactor that make it possible to obtain a significant NO reduction. To obtain different local stoichiometries under laminar flow while maintaining constant global conditions, the amount of nitrogen added through the inner and outer streams has been varied, but the total flow rate introduced in the reactor is maintained constant. Thus, the local concentration of the reactants added change and, therefore, the local stoichiometry also changes. Moreover, the variation of the nitrogen added through each stream modifies the initial velocity profile at the entrance of the reactor and, therefore, the velocity flow field in the reactor also changes. To take into account the nitrogen introduced through each stream, the velocity ratio (R) has been defined as follows:

R)

velocity flow of the inner stream velocity flow of the outer stream

Figure 9 shows the theoretical and experimental influence of the velocity ratio (R) on NO reduction, as a function of temperature, when reburning fuel is added through the inner stream and NO and O2 through the outer stream. An example of the concentration field in the reactor for this configuration is shown in Figure 3. It is worth mentioning that, applying the PFR model, the simulation of both velocity ratios is indeed the same case of study and, consequently, the same results are obtained. The BLM and LMM models show slight differences between the results obtained for both velocity ratios, which may indicate that the influence of the initial velocity profile, which is only taken into account by the BLM model to calculate the velocity flow field in the reactor (velocity flow field is fixed in LMM), has a low impact on the results.

Ind. Eng. Chem. Res., Vol. 46, No. 11, 2007 3535

Figure 9. Experimental and theoretical influence of the velocity ratio (R). Reburning fuel is introduced through the inner tube (experiments 3 and 6, in Table 1).

Despite the differences being small, a greater NO reduction is obtained for R ) 1 and high temperatures. These differences are even lower at temperatures of 0.7, a greater reduction of NO is obtained under laminar mixing conditions. This is a probe of evidence that a delayed mixing may result in an improvement of reburning performance. Moreover, it is worth mentioning that, under PFR conditions, the evolution in the NO reduction, as a function of the global value of λ, is sharper than that under non-premixed laminar conditions. Stoichiometry and residence time distribution in the reactor have an important role in reburning performance, even with the same global stoichiometry. For low temperatures, a higher NO reduction is obtained when local fuel-lean conditions have a longer residence time (reburning fuel introduced through the inner stream); however, for high temperatures, a higher NO reduction is obtained when local fuel-rich conditions have a longer residence time (reburning fuel introduced through the outer stream). LMM and BLM follow these trends, which are also shown by the experimental results. Acknowledgment The authors express their gratitude to the EU (Contract No. ENK5-CT-2000-00324) and the MCYT (Project No. PPQ20014913-E) for financial support, and to Mr. Jose´ Manuel Martos and Tinus Cording for their cooperation in the early parts of the present investigations. Nomenclature Cl ) molar concentration of the lth species (kmol/m3) Dkm ) diffusion coefficient of species k in the mixture (kg m-1 s-1) ki ) rate constant of the ith reaction MWk ) molecular weight of species k (kg/kmol) p ) pressure (N/m2) R ) radius (m) Vx ) velocity in the x-direction (m/s) Yk ) mass fraction of species k in the mixture Greek Symbols µ ) viscosity (kg m-1 s-1)

Ind. Eng. Chem. Res., Vol. 46, No. 11, 2007 3537

νi,k ) stoichiometric coefficient of the kth species in the ith reaction F ) density (kg/m3) ωk ) molar rate of formation of species k (kmol m-3 s-1) Superscripts f ) forward reaction r ) reverse reaction Literature Cited (1) Cha, C. M.; Kramlich, J. C. Modeling Finite-Rate Mixing Effects in Reburning Using a Simple Mixing Model. Combust. Flame 2000, 122, 151. (2) Chen, S. L.; McCarthy, J. M.; Clark, W. D.; Heap, M. P.; Seeker, W. R.; Pershing, D. W. Bench and Pilot Scale Process Evaluation of Reburning for In-Furnace NOx Reduction. Proc. Combust. Inst. 1986, 21, 1159. (3) Rota, R.; Bonini, F.; Servida, A.; Morbidelli, M.; Carra´, S. Modelling of the Reburning Process. Combust. Sci. Technol. 1997, 123, 83. (4) Takeya, R.; Nakamura, T. Numerical and Experimental Investigation on the Application of Natural Gas Reburning to Municipal Solid Waste Incinerators. In EnVironmental Control/Fuels and Combustion Technologies/ Nuclear Engineering; Proceedings of the 1997 International Joint Power Generation Conference, EC Series, Vol. 5; American Society of Mechanical Engineers (ASME): New York, 1997; p 389. (5) Adams, B. R.; Harding, N. S. Reburning Using Biomass for NOx Control. Fuel Process. Technol. 1998, 54, 249. (6) Weber, R.; Wecel, G.; Verlaan, A.; Breussin, F.; Dugue, J. Experimental and Numerical Studies on Reburn Jet Penetration and Mixing with Application to Boilers and Municipal Waste Incinerators. J. Inst. Energy 1998, 71, 94. (7) Lissianski, V. V.; Zamansky, V. M.; Maly, P. M.; Sheldon, M. S. Reburning Chemistry-Mixing Model. Combust. Flame 2001, 125, 1310. (8) Zarnescu, V.; Pisupati, S. V. The Effect of Mixing Model and Mixing Characteristics on NOx Reduction during Reburning. Energy Fuels 2001, 15, 363. (9) Alzueta, M. U.; Bilbao, R.; Millera, A.; Glarborg, P.; Østberg, M.; Dam-Johansen, K. Modelling Low-Temperature Gas Reburning. NOx Reduction Potential and Effects of Mixing. Energy Fuels 1998, 12, 329. (10) Zamansky, V. M.; Sheldon, M. S.; Maly, P. M. Enhanced NOx Reduction by Interaction of Nitrogen and Sodium Compounds in the Reburning Zone. Proc. Combust. Inst. 1998, 27, 3001. (11) Cha, C. M.; Kramlich, J. C.; Kosa´ly, G. Finite-Rate Mixing Effects in Reburning. Proc. Combust. Inst. 1998, 27, 1427. (12) Han, D.; Mungal, M. G.; Zamansky, V. M.; Tyson, T. J. Prediction of NOx Control by Basic and Advanced Gas Reburning Using the TwoStage Lagrangian Model. Combust. Flame 1999, 119, 483. (13) Liesa, F.; Alzueta, M. U.; Millera, A.; Bilbao. R. An Experimental and CFD Simulation Study of Reburning under Turbulent Mixing Conditions. Energy Fuels 2005, 19, 833. (14) Kolb, T.; Janshon, W.; Leuckel, W. Reduction of NOx Emissions in Turbulent Combustion by Fuel Staging. Effects of Mixing and Stoichiometry in the Reduction Zone. Proc. Combust. Inst. 1988, 22, 1193. (15) Ballester, J.; Fueyo, N.; Dopazo, C.; Herna´ndez; M.; Vidal, P. J. Influence of Operational Parameters on the Results of Reburning in Coal Combustion. In Proceedings of the 3rd International Conference on Technologies and Combustion for a Clean EnVironment, II, Lisbon, Portugal, 1995; pp 27.3/17. (16) Liesa, F.; Alzueta, M. U.; Millera, A.; Bilbao. R. Influence of Reactant Mixing in a Laminar Flow Reactor: The Case of Gas Reburning. 1. Experimental Study. Ind. Eng. Chem. Res. 2007, 46, 3520-3527. (17) Glarborg, P.; Miller, J. A.; Kee, R. J. Kinetic Modeling and Sensitivity Analysis of Nitrogen Oxide Formation in Well-Stirred Reactors. Combust. Flame 1986, 6, 177. (18) Miller, J. A.; Bowman, C. T. Mechanism and Modeling of Nitrogen Chemistry in Combustion. Prog. Energy Combust. Sci. 1989, 15, 287. (19) Kilpinen, P.; Glarborg, P.; Hupa, M. Reburning Chemistry: A Kinetic, Modeling Study. Ind. Eng. Chem. Res. 1992, 31, 1477. (20) Glarborg, P.; Alzueta, M. U.; Dam-Johansen, K.; Miller, J. A. Kinetic Modelling of Hydrocarbon/Nitric Oxide Interactions in a Flow Reactor. Combust. Flame 1998, 115, 1. (21) Smith, G. P.; Golden, D. M.; Frenklach, M.; Moriarty, N. W.; Eiteneer, B.; Goldenberg, M.; Bowman, C. T.; Hanson, R. K.; Song, S.; Gardiner, W. C.; Lissianski, V. V.; Win, Z. The GRI Mechanism, v. 3.0., 1999, http://www.me.berkeley.edi/gri_mech/.

(22) Glarborg, P.; Lilleheie, N. I.; Byggstoyl, S.; Magnussen, F.; Kilpinen, P.; Hupa, M. A Reduced Mechanism for Nitrogen Chemistry in Methane Combustion. Proc. Combust. Inst. 1992, 24, 889. (23) Bilbao, R.; Alzueta, M. U.; Millera, A.; Duarte, M. Simplified Kinetic Model of the Chemistry in the Reburning Zone Using Natural Gas. Ind. Eng. Chem. Res. 1995, 34, 4540. (24) Marro, M. A.; Pivovarov, M. A.; Miller, J. H. Strategy for the Simplification of Nitric Oxide Chemistry in a Laminar Methane/Air Diffusion Flamelet. Combust. Flame 1997, 11, 208. (25) Hewson, J. C.; Bollig, M.; Byggstoyl, S.; Magnussen, F.; Kilpinen, P.; Hupa, M. Reduced Mechanisms for NOx Emissions from Hydrocarbon Diffusion Flames. Proc. Combust. Inst. 1996, 26, 2171. (26) Mallampali, H. P.; Fletcher, T. H.; Chen, J. Y. Evaluation of CH4/ NOx Global Mechanisms Used for Modelling Lean Premixed Turbulent Combustion of Gas Natural. Presented at the Fall Meeting of the Western States Section of the Combustion Institute, Pittsburgh, PA, 1996, Paper No. 96F-098. (27) Han, X.; Ru¨ckert, F.; Schnell, U.; Hein, K. R. G. Computational Modelling of the NOx Reduction Process by Hydrocarbon Reburning with Reduced Kinetics. In Proceedings of the 6th International Conference on Technologies and Combustion for a Clean EnVironment, I, Porto, 2001; pp 6.2/179-6.2/188. (28) Giral, I.; Alzueta, M. U. An Augmented Reduced Mechanism for the Reburning Process. Fuel 2002, 81, 2263. (29) Bird, R. B.; Stewart, W. E.; Lightfoot, E. N. Transport Phenomena; Wiley; New York, 1960. (30) Kee, R. J.; Dixon-Lewis, G.; Warnatz, J.; Coltrin, M. E.; Miller, J. A.; Moffat, H. K. A Fortran Computer Code Package for the Evaluation of Gas Phase Multicomponent Transport Properties; Report No. SAND868246. Sandia National Laboratories, Livermore, CA, 1986. (31) Liesa, F. Influencia de la mezcla en el proceso de reburning, Ph.D. Thesis, University of Zaragoza, Zaragoza, Spain, 2004. (32) Branch, M. C.; Kee, R. J.; Miller, J. A. A Theoretical Investigation of Mixing Effects in the Selective Reduction of Nitric Oxide by Ammonia. Combust. Sci. Technol. 1982, 29, 147. (33) Takagi, T.; Xu, Z. Numerical Analysis of Laminar Diffusion Flames-Effects of Preferential Diffusion of Heat and Species. Combust. Flame 1994, 96, 50. (34) Figueira da Silva, E.; Pe´rez Guerrero, J. S.; Cotta, R. M. Integral Transform Solution of Boundary Layer Equations in Stream Function-Only Formulation. Int. J. Non-Linear Mech. 1999, 34, 51. (35) Revuelta, A.; Sa´nchez, A. L.; Lin˜a´n, A. Laminar Mixing in Diluted and Undiluted Fuel Jets Upstream from Lifted Flames. Combust. Flame 2002, 128, 199. (36) Villadsen, J.; Michelsen, M. L. Solution of Differential Equation, Models by Polynomial Approximation; Prentice Hall: Engelwood Cliffs, NJ, 1978. (37) Villadsen, J. V.; Stewart, W. E. Solution of Boundary Value Problems by Orthogonal Collocation. Chem. Eng. Sci. 1995, 50, 3981. (38) Rice, G. R.; Do, D. D. Applied Mathematics and Modeling for Chemical Engineers; Wiley: New York, 1995. (39) Petzold, L. R. A Description of DASSL; Report No. SAND828637. Sandia National Laboratories, Livermore, CA, 1982. (40) Patankar, S. V. Numerical Heat Transfer and Fluid Flow; McGrawHill: New York, 1980. (41) Bilbao, R.; Millera, A.; Alzueta, M. U.; Prada, L. Dilution and Stoichiometry Effects on Gas Reburning: An Experimental Study. Ind. Eng. Chem. Res. 1997, 36, 2440. (42) Alzueta, M. U.; Glarborg, P.; Dam-Johansen, K. Low Temperature Interactions between Hydrocarbons and Nitric Oxide. An Experimental Study. Combust. Flame 1997, 109, 25. (43) Miller, C. A.; Touati, A. D.; Becker, J.; Wendt, J. O. L. NOx Abatement by Fuel-Lean Reburning: Laboratory Combustor and Pilot-Scale Package Boiler Results. Proc. Combust. Inst. 1998, 27, 3189. (44) Glickert, R. W.; Houy, S. P. Application of Fuel Lean Gas Reburn Technology at Commonwealth Edison’s Joliet Generating Station 9. In 1998 Proceedings, Spring Technical Conference, ASME Internal Combustion Engine DiVision: Engine Emissions and EnVironmental Issues, Vol. 1, Ft. Lauderdale, FL, April 26-29, 1998; p 75.

ReceiVed for reView July 20, 2006 ReVised manuscript receiVed February 12, 2007 Accepted March 6, 2007 IE060944I