Design of experiments and empirical models for up to date burners design for process industries

Similar documents
DEVELOPMENT OF CFD MODEL FOR A SWIRL STABILIZED SPRAY COMBUSTOR

Heat Exchangers for Condensation and Evaporation Applications Operating in a Low Pressure Atmosphere

ADVANCED DES SIMULATIONS OF OXY-GAS BURNER LOCATED INTO MODEL OF REAL MELTING CHAMBER

Advanced Turbulence Models for Emission Modeling in Gas Combustion

ASSET INTEGRITY INTELLIGENCE. Featured Article. ACHIEVING A COMPREHENSIVE FIRED HEATER HEALTH MONITORING PROGRAM By Tim Hill, Quest Integrity Group

Characterisation of Deposits on Membrane Walls of Steam Generators by Heat Flux Density Measurement

MODELLING OF SINGLE-PHASE FLOW IN THE STATOR CHANNELS OF SUBMERSIBLE AERATOR

Experimental Study on the Non-reacting Flowfield of a Low Swirl Burner

The influence of the spectrum of jet turbulence on the. stability, NOx emissions and heat release profile of. pulverised coal flames. Ph.D.

Jet Aircraft Propulsion Prof. Bhaskar Roy Prof A M Pradeep Department of Aerospace Engineering Indian Institute of Technology, Bombay

NUMERICAL INVESTIGATION OF THE EFFECT OF THE INSULATION THICKNESS ON THE DEGREE OF NON-UNIFORMITY OF THE BILLET TEMPERATURE

COMBUSTION OF FUEL 12:57:42

XXXVIII Meeting of the Italian Section of the Combustion Institute

CHARACTERISTICS OF ELLIPTIC CO-AXIAL JETS

Lecture 9 Laminar Diffusion Flame Configurations

Overview of Turbulent Reacting Flows

Well Stirred Reactor Stabilization of flames

Numerical Simulation of a Complete Francis Turbine including unsteady rotor/stator interactions

On Annular Impinging Jets - Experimental Data Analysis

3. Write a detailed note on the following thrust vector control methods:

Rocket Propulsion Prof. K. Ramamurthi Department of Mechanical Engineering Indian Institute of Technology, Madras

Lecture 28. Key words: Heat transfer, conduction, convection, radiation, furnace, heat transfer coefficient

Direct fired strip preheating

EVALUATION OF FOUR TURBULENCE MODELS IN THE INTERACTION OF MULTI BURNERS SWIRLING FLOWS

Model of Mass and Energy Transfer in a Clinker Rotary Kiln

The Effect of Magnetic Roller Field Uniformity on Print Quality in Electrophotography

Advanced Lab Course. X-Ray Photoelectron Spectroscopy 1 INTRODUCTION 1 2 BASICS 1 3 EXPERIMENT Qualitative analysis Chemical Shifts 7

TOPICAL PROBLEMS OF FLUID MECHANICS 97

CFD study of gas mixing efficiency and comparisons with experimental data

The Fire Induced Ceiling Jet Revisited. Ronald L. Alpert, ScD, FSFPE Editor, Journal of Fire Protection Engineering Rockport, Cape Ann, MA USA

Reacting Flow Modeling in STAR-CCM+ Rajesh Rawat

Analysis of Trace (mg/kg) Thiophene in Benzene Using Two-Dimensional Gas Chromatography and Flame Ionization Detection Application

ANSYS Advanced Solutions for Gas Turbine Combustion. Gilles Eggenspieler 2011 ANSYS, Inc.

INTERNAL FLOW IN A Y-JET ATOMISER ---NUMERICAL MODELLING---

The Critical Velocity and the Fire Development

THERMAL ANALYSIS OF SECOND STAGE GAS TURBINE ROTOR BLADE

57:020 Mechanics of Fluids and Transfer Processes OVERVIEW OF UNCERTAINTY ANALYSIS

Best Practice Guidelines for Combustion Modeling. Raphael David A. Bacchi, ESSS

of the heat is dissipated as a result of the flow of electrical current in various conductors. In order to predict temperatures

Flame / wall interaction and maximum wall heat fluxes in diffusion burners

Laboratory Studies of Fire Whirls (preliminary)

A New Approach to CFD Research: Combining AVL s Fire Code with User Combustion Model

Flow Structure Investigations in a "Tornado" Combustor

Introduction. Mathematical model and experimental

EFFECT OF DISTRIBUTION OF VOLUMETRIC HEAT GENERATION ON MODERATOR TEMPERATURE DISTRIBUTION

One-Dimensional Isentropic Flow

ROLLER BEARING FAILURES IN REDUCTION GEAR CAUSED BY INADEQUATE DAMPING BY ELASTIC COUPLINGS FOR LOW ORDER EXCITATIONS

Regenerator design study combining numerical simulations and statistical tools

Low Inductance Ceramic Capacitor (LICC)

IMPROVED EVALUATION OF RECOVERY BOILER WATER CIRCULATION DESIGN WITH THE HELP OF STATE-OF-THE- ART CFD-BASED HEAT FLUX DATA

Contents. Preface... xvii

CFD Analysis of High Temperature and High Velocity System: Plasma Torch

Aerodynamics of Centrifugal Turbine Cascades

PROCESS CONTROL BASIS FOR A COST-EFFECTIVE SELECTIVE SOLDERING PROCESS

AAE COMBUSTION AND THERMOCHEMISTRY

EXPERIMENTAL DESIGN TECHNIQUES APPLIED TO STUDY OF OXYGEN CONSUMPTION IN A FERMENTER 1

Experimental and numerical investigations of heat transfer and thermal efficiency of an infrared gas stove

Efficient runner safety assessment during early design phase and root cause analysis

Environmental chemistry

JPM. 1. Introduction. Excerpt from the Proceedings of the 2016 COMSOL Conference in Munich. requirement is the. through. efficient design of

mywbut.com Hydraulic Turbines

Large-eddy simulation of an industrial furnace with a cross-flow-jet combustion system

DARS overview, IISc Bangalore 18/03/2014

The integration of response surface method in microsoft excel with visual basic application

White Paper FINAL REPORT AN EVALUATION OF THE HYDRODYNAMICS MECHANISMS WHICH DRIVE THE PERFORMANCE OF THE WESTFALL STATIC MIXER.

MECHANICAL ENGINEERING (ME)

Liquid Feed Injection in a High Density Riser

Numerical Simulation and Experimental Analysis of an Oxygen-Enriched Combustion Fiberglass Furnace

Washington University in St. Louis Chemistry Tournament Sample Problems for Individual Round #3: Kinetics, Electrochemistry, and Thermodynamics

Onset of Flow Instability in a Rectangular Channel Under Transversely Uniform and Non-uniform Heating

OXY-COAL COMBUSTION: EFFECTS OF P O2 ON COAL JET STABILITY IN O 2 /CO 2 ENVIRONMENTS

Prediction of Transient Deflector Plate Temperature During Rocket Plume Impingment and its Validation through Experiments

A NUMERICAL STUDY ON THE FLOW AND HEAT TRANSFER FOR THE INSIDE OF A NEW DIVERSION-TYPE LNG HEATING DEVICE

NUMERICAL SIMULATION OF LDI COMBUSTOR WITH DISCRETE-JET SWIRLERS USING RE-STRESS MODEL IN THE KIVA CODE

Analysis of NO-Formation for Rich / Lean - Staged Combustion

Hydrogen addition to the Andrussow process for HCN synthesis

Atomization. In Flame Emission

CFD Modeling Ensures Safe and Environmentally Friendly Performance in Shell Claus Off-Gas Treating (SCOT) Unit

CHAPTER 4 OPTIMIZATION OF COEFFICIENT OF LIFT, DRAG AND POWER - AN ITERATIVE APPROACH

HETEROGENEOUS CONDENSATION FOR SUBMICRONIC PARTICLE ABATEMENT: EXPERIMENTS AND MODELLING

4.1 Derivation and Boundary Conditions for Non-Nipped Interfaces

Numerical Modeling of Anode Baking Process with COMSOL Multiphysics

THE APPLICATION OF PROCESS MASS SPECTROMETRY TO FUMED SILICA PRODUCTION

CHAPTER 3 MODELLING AND ANALYSIS OF THE PACKED COLUMN

STATOR/ROTOR INTERACTION

EffectofBaffleDesignontheOff-bottomSuspensionCharacteristics ofaxial-flowimpellersinapilot-scalemixingvessel

Evaluation of the Air-Demand, Flame Height, and Radiation from low-profile flare tips using ISIS-3D

Heat and Mass Transfer Prof. S.P. Sukhatme Department of Mechanical Engineering Indian Institute of Technology, Bombay

COMSOL Multiphysics Simulations of the Electric Field and Gas Flow in a Microwave Axial Injection Torch

QIRT th International Conference on Quantitative InfraRed Thermography

The Effect of Endplates on Rectangular Jets of Different Aspect Ratios

The Dynamical Loading of the WWER440/V213 Reactor Pressure Vessel Internals during LOCA Accident in Hot and Cold Leg of the Primary Circuit

EXPERIMENT 6 Empirical Formula of a Compound

THERMO-FLOW CHARACTERISTICS OF A PIN-FIN RADIAL HEAT SINKS ACCORDING TO THEIR FIN HEIGHT PROFILE

NGSS. Science Items Grade 5 Middle School High School. Table of Contents. Grade Middle School... 5 High School... 10

Thermoacoustic Instabilities Research

Introduction to the course ``Theory and Development of Reactive Systems'' (Chemical Reaction Engineering - I)

THINK Fluid Dynamix THINK FLUID DYNAMIX. CFD for Aeration Applications

AN OPTIMAL RELATIONSHIP BETWEEN CASTING SPEED AND HEAT TRANSFER COEFFICIENTS FOR CONTINUOUS CASTING PROCESS

UQ in Reacting Flows

Transcription:

Design of experiments and empirical models for up to date burners design for process industries Proceedings of European Congress of Chemical Engineering (ECCE-6) Copenhagen, 16-20 September 2007 Design of experiments and empirical models for up to date burners design for process industries Petr Belohradsky, Vit Kermes, Petr Stehlik Institute of Process and Environmental Engineering, Faculty of Mechanical Engineering, Brno University of Technology (UPEI VUT), Technická 2, 616 69 Brno, Czech Republic Tel.: +420 541 142 368, Fax: +420 541 142 372, E-mail: belohradsky@upei.fme.vutbr.cz Abstract At present period the application of complex models based on computational fluid dynamics (CFD) increases. At the same time, the cost of full experiments with high number of factors grows rapidly. In some cases, however, experimental approach cannot be fully substituted by computational modeling, which applies also to tuning of design of new burners. Due to high cost of experiments therefore grows the importance of properly designed experiments and statistical data analysis. Statistical experimental design and analysis is an indispensable tool for experimenters and one of the core topics in a statistics curriculum. The present work deals with the design of experiment from its preparation, through experiment planning, that is followed by preliminary considerations about the model, implementation of the experiment and subsequent evaluation. Furthermore, we point out the problems of the performing of full experiment with high number of factors and outline the possibilities of creating the basis for fractional experiment. The experiment presented in this work is carried out at an up to date experimental facility for testing diffusion gas burners. The investigation focuses on the influence of design parameters on the formation of nitrogen oxides. Priority is however always given to the functionality of the burner (namely flame stability). This paper shows a convenient way how to set up the plan of expensive and time consuming experiments in testing industrial burners. Importance for industrial practice, equipment and consequently processes is demonstrated. Keywords: design of experiments, statistical data analysis, empirical model, formation of nitrogen oxides 1. Introduction Experiments are performed by investigators in virtually all fields of inquiry, usually to discover something about a particular process or system. A common goal in many - 1 -

Petr Belohradsky et al. types of experiments is to characterize the relationship between a response and a set of quantitative factors of interest of the researcher. This is accomplished by constructing a model that describes the response over the applicable ranges of the factors of interest. However, a little attention is paid to efficient carrying out of experiments during construction of the empirical models. Badly or insufficiently planned experiment can cause incorrect evaluation of acquired experimental data by employing any of statistical data analysis. Strictly speaking, it is necessary to set up such conditions in order to achieve the range of the experiment as small as possible and simultaneously high-quality information. This requirement is most important in experiments with high number of factors (independent variables). The design of experiment depends on particular solved problem, namely on the number of factors, the goal of experiment, the homogeneity or variability of experimental conditions. It is required to take into account the economical criteria, too. All the goals are fulfilled by properly designed experiments as well as by appropriate technical carrying out of experiments. Experimental design is an important tool in the engineering world for improving the performance of process and it also has extensive application in the development of new processes resulting in e.g. improved process yields or reduced overall costs. Experimental design methods also play a major role in engineering design activities including for example evaluation and comparison of basic design configurations or determination of key product design parameters that impact product performance. The development of efficient, low polluting combustion systems is a major goal of combustion researchers. Burner is always one of the key components of any combustion system. There are many factors related to the design of a burner that have significant impact on its emissions from its flame. There are three predominant mechanisms for nitrogen oxides formation from combustion processes (Baukal, 2004): thermal NO x, prompt NO x, fuel NO x, where the thermal NO x formation is the dominating source in the combustion of natural gas. Nitric oxide is the primary NO x compound emitted from industrial combustion sources. It is even emitted if natural gas or hydrogen is used because it can be formed from air nitrogen and oxygen at high temperatures. The dependence of NO x emissions concentration on various operating parameters is generally known. It was shown that the main factors which to a great extent influence NO x formation are the excess air and the temperature of preheated air (Carvalho et al., 1990, 1992). NO x increases with increasing combustion temperature and with increasing average furnace chamber temperature. Most NO x reducing techniques like staged combustion, gas re-circulation or water injection try to cut off peak - 2 -

Petr Belohradsky et al temperatures, keep the residence time in high temperature areas low and avoid high oxygen concentration in these areas. 2. Description of the testing facility Experimental data are presently obtained from testing a low-no x burner in an up to date experimental facility intended for testing and improving industrial burners. The main apparatus of the burner experimental facility is horizontal water cooled combustion chamber with a maximum heat load approximately 1.8 MW (internal diameter of 1 meter, length of 4 meters). The cooling of the combustion chamber allows simulating conditions similar to conditions in the heating furnaces of the process industry. The cooling shell is split up to seven sections that gives the possibility to evaluate intensity of the heat flux along the flame. The next positive effect of the splitting in the sections is an improved water circulation in the cooling shell, thus avoiding local boiling caused by poor water circulation. There is one inspection hole (window) per each section that enables detailed observation of the flame and installation of measuring instrumentation. The experimental facility is equipped with automatic data collection (temperatures, pressures, flows). The picture of the experimental facility is shown in Figure 1. There was a possibility to measure concentrations of NO, NO 2, CO and O 2 by flue gas analyzer. Other measured data are as follows: flue gas temperature, temperature of combustion air, flow rate of combustion air and flow rate of fuel. Experimental data is measured during testing of a gas burner with forced inflow of combustion air. Nominal duty of the gas burner was 0.93 MW. Gas injector stabilizing burner is used for flame stabilization. Figure 1: Experimental facility. - 3 -

Petr Belohradsky et al. 3. Design of experiment By statistical design of experiment, we refer to the process of planning the experiment so that appropriate data that can be analyzed by statistical methods would be collected, resulting in valid and objective conclusions. The presented design of experiment was employed with respect to the basic principles of experimental design, namely to replication, randomization, and blocking (Montgomery, 1991; Wu et al., 2000; Mason et al., 2003). The following chapters summarize the important steps that we had to address during the designing of experiment. 3.1. Objective of experiment The objective of the experiment is to investigate the influence of burner design parameters on the formation of nitrogen oxides and to find the relationship between the design parameters and the response. The tested gas burner is of a diffusion-mixed type. Fuel staging was used, which means that secondary injectors in the burner are used to inject a portion of the fuel downstream of the root of the flame. Non-preheated atmospheric air is used as the oxidizer. Natural gas was used as a fuel. 3.2. Choice of factors and levels. Response variable. A review of the input factors (design parameters) which are varied in the experiment and the ranges of these factors may be found in Table 1. Limits of the factors were chosen to meet practical constraints. Response variable is the concentration of nitrogen oxides. Factors Lower level Upper level Burner duty [kw] 765 1120 Diameter of swirl generator [mm] 240 280 Pitch angle of swirl generator blades [ ] 35 55 Diameter of primary gas throttle [mm] 5,5 6,5 Air equivalence ratio [-] 1,1 1,2 Pitch angle of secondary nozzles [ ] 20 40 Tangential orientation of secondary nozzles [ ] 0 45 Radial position of secondary nozzles [mm] 0 50 Axial position of secondary nozzles [mm] 0 80 Table 1: Burner design parameters. Besides the response variable (nitrogen oxides concentration) are monitored also characteristic features of the flame like its diameter, length and shape, for each of the measurements. The character of these observations is however rather qualitative. - 4 -

Petr Belohradsky et al 3.3. Experimental plan In many industrial applications the response surface methodology is an effective tool to find the adequate model. Since the curvature effects are expected, the goal is to fit the second-order (quadratic) model in k input factors 0 k k k 2 i i ij i j ii i i= 1 i< j i= 1 ε, y = β + β x + β xx + β x + where β i represents the linear effect of x i, β ij represents the bilinear interaction between x i and x j, and β ij represents the quadratic effect of x i, ε is an error term. An experimental design for fitting the second-order model must have at least three levels of each factor. The so-called central composite design is the most widely-used experimental design approach for fitting a second-order response surface. The design enables to construct the second-order model without needing to use the full three level factorial experiments. In order to carry out the full three level factorial experiment with high number of factors, extremely high number of test runs would be required, approximately 20000. This would be very time consuming and costly. Due to limitations given by nature of the present experiment, the central composite design was selected, which requires far less measurements than in a 3 k factorial experiment design. The actual design employed in the experiment is a special form of the central composite design in that the axial points are located on the centers of the faces. Thus, the design is called a face-centered central composite design and uses only three levels of each factor. The design consists of the following three parts: 2 1) 2 k fractional factorial layout, where k = 9 is the number of the input factors. 9 2 Thus the number of cube points was n f = 2 = 128. 2) 2k axial points along the coordinate axis. Each pair of axial points is denoted, using coded levels, as follows: ( ± α,0,0,,0), (0, ± α,0,,0), (0,0,0,, ± α), where α = 1for the face-centered central composite design. 3) n = 1 center points (0,0,,0). c k 2 The number of test runs in the plan is equal to N = 2 + 2k+ n c = 147without replication. In order to be able to estimate the experimental error and to increase the reliability of conclusion, the basic experiment was replicated three times. Thus, the total number of test runs in the experiment is 441. - 5 -

Petr Belohradsky et al. Since in general most new settings of factor levels for individual test runs would require shutdown and dismounting of the burner, the order of test runs in the experiment is not completely random. The order of test runs was arranged in such a way that the setting of each test run requires the smallest possible number of changes of factor levels (compared to the respective preceding run). The benefit of resulting time-savings in the factor-setting stages is that the temperature of the combustion chamber walls does not decrease too much between individual test runs. Thus further time and fuel savings result, as equilibrium conditions for new measurement are attained quickly. The technique of blocking is not applied to the plan of experiment even though the temperature of combustion air fluctuates during the whole experiment (depending on ambient conditions in the laboratory). The variation of temperature of combustion air is about ±10 C around the mean value. This however has a negligible influence on the measured response. 4. Preliminary results At the moment of completing this text, only a part of measurements is completed. Expected time interval of test runs, necessary to reach steady temperature of the combustion chamber walls in each of them has been confirmed and is approximately 20 minutes. This period is minimized by suitable ordering of test runs as described above. The preliminary results reveal that the concentration of NO x emissions has decreased with increasing radial position of secondary nozzles as well as with increasing clockwise tangential orientation of secondary nozzles (reference orientation is that the secondary nozzles are oriented to the burner axis). Next, the concentration of nitrogen oxides has increased with increasing pitch angle of secondary nozzles. The change of axial position of secondary nozzles has shown a small influence on the change of NO x concentration. 5. Conclusion The present work presents the proper way to create a smallest sufficient plan of experiment while respecting statistical requirements of the theory of experiment design. The principles of experiment design were implemented during the construction of a plan of measurements for testing a low-no x gas burner. Facecentered composite design of experiment was proposed that is sufficient for this kind of experiment in spite of small number of test runs. The presented design that was proposed is specific for the given type of burner. However, the methodology of proper experiment design can be applied to various types of burners. - 6 -

Petr Belohradsky et al 6. Future work Future work will consist of carrying out the measurements according to the designed plan of experiment. This will be followed by construction of the empirical model by employing statistical data analysis. The aim of the model is further improvement of the predictive capabilities of NO x in the field of burner design. References Montgomery, D. C., Design and Analysis of Experiments, John Wiley & Sons, Inc., USA (1991). Wu, C. F. J., Hamada, M., Experiments: Planning, Analysis, and Parameter Design Optimization, John Wiley & Sons, Inc., USA (2000). Mason, R. L., Gunst, R. F. and Hess, J. L., Statistical Design and Analysis of Experiments with Applications to Engineering and Science, John Wiley & Sons, Inc., USA (2003). Baukal, C. E., Industrial Combustion Pollution and Control, Marcel Dekker, Inc., USA (2004). Bebar L., Kermes V., Stehlik P., Canek J. and Oral J., (2002) Waste Management, 22, 443-451. Carvalho, M. G., Semiao, V. S., Coelho, J.P., (1192) ASME J. Engineering for Industry, 114, 514-523. Carvallho, M. G., Semiao, V. S., Lockwood, F. C., Papadopoulos, C., (1990) Journal of Inst. of Energy March, 39-47. - 7 -