Chemical Reactor Network modeling B. Rosati Prediction offrom emissions combustiontems using sys- 0Dreacting flow and models 1D ehiceUiesti Delft Universiteit Technische
Prediction of emissions from combustion systems using 0D and 1D reacting flow models Prediction of emissions from combustion systems using 0D and 1D reacting flow models Chemical Reactor Network modeling
by B. Rosati
to obtain the degree of Master of Science at the Delft University of Technology, to be defended publicly on Monday December 14, 2015 at 10:00 AM.
Student number: 4259343 Project duration: November 15, 2014 – December 14, 2015 Thesis committee: Prof. dr. D.J.E.M. Roekaerts, TU Delft, supervisor Dr. A. Gangoli Rao, TU Delft, co-supervisor Dr. A. Bhat, Honeywell, co-supervisor Dr. H. B. Eral, TU Delft
Faculty of Mechanical, Maritime and Materials Engineering (3mE). Process & Energy Department. Delft University of Technology
An electronic version of this thesis is available at http://repository.tudelft.nl/.
Abstract
Emission prediction is a complex problem involving the coupling between the flow field and chemistry. Most of the time CFD is the preferred modeling approach, yielding predictions with varying degrees of accuracy. But because of a high computational cost, CFD investigations are often limited to the use of reduced chemical mechanisms. In this work the specific features of chemical reactor networks are exploited to build a fast and reliable emission estimator. The main advantage of this modeling approach is a much lower computational cost than CFD, hence offering the potential for relatively fast predictions while allowing the use of detailed chemistry.
This methodology has been applied to three different combustion systems, with mixed results. It may not be the most suitable modeling technique to obtain emissions from a lifted jet flame, but a successful estimator has been designed for flameless furnaces. It is based solely on analytical sub– models, giving it the potential to predict the emissions from any type of flameless furnace installation. For three different experimental setups, the correct trends were reproduced as well as the right order of magnitude for 푁푂 and 퐶푂 emissions, if not within experimental measurements uncertainty. Finally the emissions from a lean-premixed gas turbine combustor burning cryogenic fuel have been successfully modeled and this investigation has brought out the major sensitivities of this system.
Lastly, despite some promising results, several developments have been suggested to improve the accuracy and stability of the flameless furnace estimator. The combustor estimator, for one, can be used as basis to investigate the behavior of the more comprehensive hybrid combustion system it has originally been designed for: the dual combustion chamber of the AHEAD hybrid engine (Advanced Hybrid Engines for Aircraft Development).
Keywords: Emissions, Reactor networks, Flameless combustion, Strong–Jet / Weak–Jet, Lean- premixed combustion.
iii
Contents
List of Figures ix List of Tables xiii 1 Introduction 1 1.1 Energy & Environment ...... 1 1.2 Design & Modeling ...... 2 1.3 Thesis objectives ...... 2 1.4 Thesis outline ...... 3 2 Literature : scientific background 5 2.1 Emissions ...... 5 2.1.1 Carbon monoxide (퐶푂) ...... 5 2.1.2 Nitrogen oxides (푁푂 ) ...... 6 2.1.3 Quantification of emissions ...... 8 2.2 Chemical mechanisms ...... 9 2.2.1 GRI-Mech 3.0...... 10 2.2.2 퐶2_푁푂 Mechanism ...... 10 2.3 Flameless combustion ...... 10 2.3.1 One technology, different names ...... 10 2.3.2 Theory and characteristics ...... 11 2.3.3 Desirable consequences of flameless combustion ...... 14 2.3.4 Operation and stability domain...... 14 2.3.5 Burner and furnace technology ...... 15 2.4 Jet modeling ...... 17 2.4.1 Free jet structure ...... 17 2.4.2 Entrainment models ...... 17 3 Emission modeling of a lifted jet flame 21 3.1 Introduction ...... 21 3.2 Additional literature ...... 23 3.2.1 Turbulent jet flame emission modeling : a 2–reactor model ...... 23 3.2.2 Diffusion flame stability...... 24 3.3 Reactor network ...... 26 3.3.1 Architecture ...... 26 3.3.2 Model inputs ...... 28 3.3.3 Simulation cases ...... 30 3.3.4 Calibration of the reactor network ...... 30 3.3.5 Conclusion ...... 30 3.4 Results ...... 32 3.5 Conclusion ...... 34 4 Emission modeling of flameless furnaces 35 4.1 Introduction ...... 35 4.2 Additional literature ...... 36 4.3 The ”Bending Model” ...... 37 4.3.1 Setup and equations ...... 37 4.3.2 The 3 jets problem ...... 38 4.3.3 Implementation of the dual–hole configuration ...... 39 4.3.4 Validation against the reported performance ...... 40 4.3.5 Validation against the SJ/WJ model predictions ...... 41
v vi Contents
4.3.6 Final validation: the IFRF furnace ...... 45 4.3.7 Final validation: the University of Mons furnace ...... 45 4.3.8 Final validation: the University of Adelaïde furnace ...... 45 4.3.9 Conclusion ...... 48 4.4 Jet entrainment modeling in a SJ/WJ configuration...... 49 4.4.1 Mutual ”shielding” of the jets ...... 49 4.4.2 Shielding factor ...... 49 4.4.3 Conclusion ...... 51 4.5 CRN modeling of a flameless furnace ...... 52 4.5.1 General architecture ...... 52 4.5.2 Role of the analytical submodels ...... 53 4.5.3 Parameters of particular interest...... 54 4.5.4 Conclusion: model parameters classification ...... 55 4.6 Prediction of emissions from the IFRF furnace ...... 56 4.6.1 Inputs ...... 56 4.6.2 Simulation cases ...... 57 4.6.3 Base case results...... 57 4.6.4 Parametric study ...... 62 4.6.5 Conclusion ...... 63 4.7 Prediction of emissions from the Mons furnace ...... 64 4.7.1 Inputs ...... 64 4.7.2 Simulation cases ...... 66 4.7.3 Results ...... 66 4.7.4 Conclusion ...... 68 4.8 Prediction of emissions from the TU Delft furnace ...... 68 4.8.1 Furnace characteristics ...... 68 4.8.2 Necessary approximations ...... 69 4.8.3 Inputs ...... 70 4.8.4 Simulation cases ...... 71 4.8.5 Results ...... 71 4.9 General conclusion ...... 73 5 Emission modeling of a jet engine lean–premixed combustor 75 5.1 Introduction ...... 75 5.2 The cryo–combustor of the AHEAD engine ...... 76 5.3 Network architecture ...... 80 5.4 Model inputs ...... 81 5.5 Calibration ...... 83 5.5.1 Heat loss calibration ...... 83 5.5.2 Recirculation intensity calibration ...... 83 5.6 Results: comparison with experimental data ...... 84 5.7 Parametric study ...... 89 5.7.1 Influence of the overall residence time: ̇푚 ...... 89 5.7.2 Influence of the heat loss intensity...... 89 5.7.3 Influence of the heat loss distribution ...... 91 5.7.4 Influence of 푁 ...... 91 5.7.5 Conclusion ...... 93 5.8 Predictions from the real combustor in operating conditions ...... 93 5.8.1 Differences with the experimental installation ...... 93 5.8.2 Simulation cases and parameters ...... 94 5.8.3 Emission prediction of the operating combustor ...... 95 5.9 Conclusion ...... 97 6 Conclusion 99 6.1 Emission modeling of a lifted jet flame ...... 99 6.2 Emission modeling of a flameless furnace ...... 100 6.3 Emission modeling of an innovative gas turbine combustor ...... 101 Contents vii
A Guide of good practices using Cantera 103 A.1 Miscellaneous ...... 103 A.2 The Plug–Flow Reactor in Cantera ...... 104 B Calibration of the lifted jet flame reactor network 105 B.1 Calibration of 퐿 ...... 105 B.2 Calibration of 푛 ...... 105 B.3 Calibration of the radiant fraction ...... 107 C Bending Model development 109 C.1 The dual–hole configuration: 5 possible integration modes ...... 109 C.1.1 The slope approximation ...... 109 C.1.2 Discussion on the initial conditions ...... 110 C.1.3 Conclusion ...... 111 C.2 Comparison with 10° diverging injection angle trajectories...... 111 C.3 Comparison with 20° diverging injection angle trajectories...... 113 C.4 Trajectory prediction comparison with the SJ/WJ model ...... 113 C.5 Convergence point prediction comparison with the SJ/WJ model ...... 113 C.6 Conclusion: selection of the preferred integration mode ...... 113 D Shielded Entrainment Model development 117 D.1 Shielding factor: method A...... 117 D.2 Shielding factor: method B...... 118 D.3 Adoption of a combined criteria ...... 121 D.4 One weakness of the Shielded Entrainment Model...... 121 E Details on the development of Mancini’s reactor network and additional simulations 123 E.1 Differences with the designed analytical estimator ...... 123 E.2 Additional simulations: impact of the mixer entrainment definition...... 124 E.3 Sensitivity analysis results ...... 125 E.4 Conclusion ...... 125 F Alternative heat loss implementation in the cryogenic combustor network 127 G User guide: flameless furnace emissions estimator 131 G.1 Program structure ...... 131 G.2 Modeling a new flameless furnace installation ...... 132 H User guide: gas turbine combustor emissions estimator 133 H.1 Program structure ...... 133 H.2 Guidelines for a proper use of the reactor networks ...... 134 Bibliography 135
List of Figures
2.1 Schematic representation of the flameless combustion technique principle ...... 11 2.2 Conventional and flameless combustion of heavy fuel oil ...... 12 2.3 Temperature field comparison between conventional and flameless combustion ..... 13 2.4 Temperature fluctuation comparison between conventional and flameless combustion . 13 2.5 Radiative heat flux and wall temperature profiles in a flameless combustion furnace .. 13 2.6 Stability diagram for different combustion regimes ...... 15 2.7 Different views of the REGEMAT® CD 200 burner ...... 16 2.8 Different views of the CGRI burner ...... 16 2.9 Han and Mungal’s free jet entrainment coefficient ...... 18 2.10 Effect of buoyancy on jet entrainment ...... 19
3.1 Schematic representation of the combustion duct of Fujimori et al. experimental setup . 22 3.2 Change of liftoff height and evolution of EINOx of pure methane jets with coflowing air temperature ...... 22 3.3 Sketch of the TSL model ...... 23 3.4 Schematic of the 2–reactor model ...... 24 3.5 Features of a lifted jet diffusion flame ...... 25 3.6 Images of lifted jet flames from the DJHC setup ...... 25 3.7 Sketch of the combustion duct discretization in 4 regions ...... 27 3.8 Sketch of the designed reactor network for the lifted jet flame in hot coflow system ... 28 3.9 Typical time evolution of the igniter mass flow ...... 28 3.10 Entrainment characteristics for 푇 = 1310 퐾 / ℎ = 65 푚푚 ...... 29 3.11 Reactor network predictions vs experimental data for the lifted jet flame ...... 33 3.12 Predicted temperature profiles along the combustion duct ...... 33 3.13 Detailed view of the peak temperature region ...... 33
4.1 The reactor network built by Mancini ...... 35 4.2 Parametrization of the Bending Model dual–hole configuration ...... 37 4.3 Parametrization of the Bending Model 3 jets configuration ...... 38 4.4 Validation of the isolated SJ hypothesis ...... 40 4.5 Trajectories predictions for 10° diverging injection ...... 42 4.6 Trajectories predictions for 20° diverging injection ...... 42 4.7 Trajectories predictions for 10° diverging injection with modified entrainment ...... 43 4.8 Trajectories predictions for 20° diverging injection with modified entrainment ...... 43 4.9 Comparison between predictions from the implemented Bending Model and the SJ/WJ model ...... 44 4.10 NFK regenerative burner of the IFRF furnace ...... 46 4.11 Convergence point prediction for the IFRF furnace ...... 46 4.12 Vertical cut of the Mons furnace base in the injectors plane ...... 46 4.13 Convergence point prediction for the Mons furnace (10% excess air) ...... 47 4.14 Convergence point prediction for the Mons furnace (20% excess air) ...... 47 4.15 Adelaïde furnace dimensions (a) and two burner configurations (b) ...... 48 4.16 Convergence point prediction comparison for the Adelaïde furnace ...... 48 4.17 Strong jet entrainment data comparison with the Ricou & Spalding model ...... 50 4.18 Weak jet entrainment data comparison with the Ricou & Spalding model ...... 50 4.19 Strong jet entrainment data comparison with the modified Ricou & Spalding model ... 50 4.20 Weak jet entrainment data comparison with the modified Ricou & Spalding model .... 51 4.21 Weak jet entrainment data comparison with the Ricou & Spalding model modified with methods A and B ...... 51
ix x List of Figures
4.22 General architecture of the designed flameless furnace reactor network ...... 52 4.23 Volume determination of a unit PSR ...... 54 4.24 Jets temperature evolution comparison ...... 60 4.25 Jets 푂 concentration evolution ...... 61 4.26 Jets 푁푂 concentration evolution ...... 61 4.27 Influence of 푛 on CPU time ...... 63 4.28 푁푂 emissions from the Mons furnace at 10% excess air ...... 67 4.29 푁푂 emissions from the Mons furnace at 20% excess air ...... 67 4.30 Sketch of the TU Delft multiple burner flameless furnace ...... 69 4.31 Prediction of 퐶푂 emission from the TU Delft flameless furnace ...... 72 4.32 Prediction of 푁푂 emission from the TU Delft flameless furnace ...... 72
5.1 ACARE vision for Europe ...... 77 5.2 Schematic of the AHEAD hybrid engine ...... 77 5.3 General layout of the hydrogen combustor ...... 77 5.4 Schematic of the burner ...... 78 5.5 Evolution of 푥 with 퐽 for various 푇 and 휙 ...... 79 5.6 Flame shape and cone angle issuing from the designed lean–premixed burner ..... 79 5.7 Basic architecture modeling the cryo combustor of the AHEAD engine ...... 80 5.8 Determination of the reactors volumes ...... 81 5.9 Comparison of 푁푂 predictions with experimental data for various 휙 and 푇 , at constant ̇푚 = 130 푘푔/ℎ ...... 85 5.10 Temperature and 푁푂 distribution along the combustor for 푇 , = 453 퐾, 휙 = 0.8 and ̇푚 = 130 푘푔/ℎ (Mixer: 0, Flame: 1, Recirculation: 2) ...... 85 5.11 푁푂 pathway analysis for 푇 , = 453 퐾 and ̇푚 = 130 푘푔/ℎ on the flame reactor 88 5.12 Influence of ̇푚 on 푁푂 emissions ...... 89 5.13 Heat loss intensity influence on the overall 푁푂 emissions ...... 90 5.14 Influence of the heat loss distribution on 푁푂 emissions ...... 91 5.15 Influence of the heat loss distribution on 푁푂 emissions for all cases ...... 92 5.16 Influence of the PFR discretization on 푁푂 emissions ...... 92 5.17 Architecture of the reactor network modeling the integrated combustor in a gas turbine . 94 5.18 푁푂 formation pathways for take–off conditions ...... 96