Lecture 5 Emission and Low-NOx Combustors

Similar documents
Lecture 17: Kinetics of Phase Growth in a Two-component System:

Lecture 18: Kinetics of Phase Growth in a Two-component System: general kinetics analysis based on the dilute-solution approximation

Lecture 26: Leapers and Creepers

, on the power of the transmitter P t fed to it, and on the distance R between the antenna and the observation point as. r r t

Two-dimensional Effects on the CSR Interaction Forces for an Energy-Chirped Bunch. Rui Li, J. Bisognano, R. Legg, and R. Bosch

Chapter 7. Interference

Orthotropic Materials

The Global Trade and Environment Model: GTEM

Today - Lecture 13. Today s lecture continue with rotations, torque, Note that chapters 11, 12, 13 all involve rotations

Lecture 22 Electromagnetic Waves

Exponential and Logarithmic Equations and Properties of Logarithms. Properties. Properties. log. Exponential. Logarithmic.

Chapter 13 Homework Answers

A Numerical Hydration Model of Portland Cement

Low-complexity Algorithms for MIMO Multiplexing Systems

An Automatic Door Sensor Using Image Processing

Lecture-V Stochastic Processes and the Basic Term-Structure Equation 1 Stochastic Processes Any variable whose value changes over time in an uncertain

L. Yaroslavsky. Image data fusion. Processing system. Input scene. Output images

MEEN 617 Handout #11 MODAL ANALYSIS OF MDOF Systems with VISCOUS DAMPING

Chapter Finite Difference Method for Ordinary Differential Equations

An Open cycle and Closed cycle Gas Turbine Engines. Methods to improve the performance of simple gas turbine plants

ÖRNEK 1: THE LINEAR IMPULSE-MOMENTUM RELATION Calculate the linear momentum of a particle of mass m=10 kg which has a. kg m s

Applied Optoelectronics in Medicine

The Production of Polarization

Pressure Vessels Thin and Thick-Walled Stress Analysis

Lecture 5. Chapter 3. Electromagnetic Theory, Photons, and Light

Molecular Evolution and Phylogeny. Based on: Durbin et al Chapter 8

CS 188: Artificial Intelligence Fall Probabilistic Models

Computer Propagation Analysis Tools

Sections 3.1 and 3.4 Exponential Functions (Growth and Decay)

MATHEMATICAL FOUNDATIONS FOR APPROXIMATING PARTICLE BEHAVIOUR AT RADIUS OF THE PLANCK LENGTH

General Non-Arbitrage Model. I. Partial Differential Equation for Pricing A. Traded Underlying Security

Sharif University of Technology - CEDRA By: Professor Ali Meghdari

FINITE DIFFERENCE APPROACH TO WAVE GUIDE MODES COMPUTATION

Servomechanism Design

7 Wave Equation in Higher Dimensions

ENGI 4430 Advanced Calculus for Engineering Faculty of Engineering and Applied Science Problem Set 9 Solutions [Theorems of Gauss and Stokes]

The k-filtering Applied to Wave Electric and Magnetic Field Measurements from Cluster

Q & Particle-Gas Multiphase Flow. Particle-Gas Interaction. Particle-Particle Interaction. Two-way coupling fluid particle. Mass. Momentum.

PHYS GENERAL RELATIVITY AND COSMOLOGY PROBLEM SET 7 - SOLUTIONS

ME 304 FLUID MECHANICS II

Support Vector Machines

DEVELOPMENT OF A DESIGN AND PERFORMANCE PREDICTION TOOL FOR THE GROUND SOURCE HEAT PUMP SYSTEM

Measures the linear dependence or the correlation between r t and r t-p. (summarizes serial dependence)

Phys-272 Lecture 17. Motional Electromotive Force (emf) Induced Electric Fields Displacement Currents Maxwell s Equations

TP B.2 Rolling resistance, spin resistance, and "ball turn"

Module 4: Time Response of discrete time systems Lecture Note 2

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

KINEMATICS OF RIGID BODIES


Probabilistic Models. CS 188: Artificial Intelligence Fall Independence. Example: Independence. Example: Independence? Conditional Independence

Feedback Couplings in Chemical Reactions

156 There are 9 books stacked on a shelf. The thickness of each book is either 1 inch or 2

u(t) Figure 1. Open loop control system

Variance and Covariance Processes

GEOGRAPHY PAPER

UT Austin, ECE Department VLSI Design 5. CMOS Gate Characteristics

Gravity. David Barwacz 7778 Thornapple Bayou SE, Grand Rapids, MI David Barwacz 12/03/2003

Overview. Overview Page 1 of 8

Combinatorial Approach to M/M/1 Queues. Using Hypergeometric Functions

Numerical Dispersion

Chapter 14 Homework Answers

Control Volume Derivation

one primary direction in which heat transfers (generally the smallest dimension) simple model good representation for solving engineering problems

8.5 Circles and Lengths of Segments

where the coordinate X (t) describes the system motion. X has its origin at the system static equilibrium position (SEP).

Lecture Angular Momentum

( ) c(d p ) = 0 c(d p ) < c(d p ) 0. H y(d p )

r r r r r EE334 Electromagnetic Theory I Todd Kaiser

Numerical Integration

Monochromatic Wave over One and Two Bars

Reinforcement learning

ME 391 Mechanical Engineering Analysis

The shortest path between two truths in the real domain passes through the complex domain. J. Hadamard

Theoretical background and the flow fields in downhole liquid-liquid hydrocyclone (LLHC)

The sudden release of a large amount of energy E into a background fluid of density

Sterilization D Values

ECE 2100 Circuit Analysis

Ferent equation of the Universe

University of Toledo REU Program Summer 2002

Chapter 7: Inverse-Response Systems

Mobility of atoms and diffusion. Einstein relation.

AB for hydrogen in steel is What is the molar flux of the hydrogen through the steel? Δx Wall. s kmole

18.03SC Unit 3 Practice Exam and Solutions

Lecture 4 Kinetics of a particle Part 3: Impulse and Momentum

Physics 2001/2051 Moments of Inertia Experiment 1

CHEMICAL KINETICS: 1. Rate Order Rate law Rate constant Half-life Temperature Dependence

Faraday s Law (continued)

Representing Knowledge. CS 188: Artificial Intelligence Fall Properties of BNs. Independence? Reachability (the Bayes Ball) Example

Introduction to SLE Lecture Notes

r P + '% 2 r v(r) End pressures P 1 (high) and P 2 (low) P 1 , which must be independent of z, so # dz dz = P 2 " P 1 = " #P L L,

Electromagnetism Physics 15b

Growth of structure in an expanding universe The Jeans length Dark matter Large scale structure simulations. Large scale structure data

Physics 240: Worksheet 16 Name

Do not turn over until you are told to do so by the Invigilator.

Physics 2212 GH Quiz #2 Solutions Spring 2016

To become more mathematically correct, Circuit equations are Algebraic Differential equations. from KVL, KCL from the constitutive relationship

The fundamental mass balance equation is ( 1 ) where: I = inputs P = production O = outputs L = losses A = accumulation

Program: RFEM 5, RSTAB 8, RF-DYNAM Pro, DYNAM Pro. Category: Isotropic Linear Elasticity, Dynamics, Member

ECE 2100 Circuit Analysis

MATHEMATICS PAPER 121/2 K.C.S.E QUESTIONS SECTION 1 ( 52 MARKS) 3. Simplify as far as possible, leaving your answer in the form of surd

Transcription:

Lecue 5 Emiion and Low-NOx Combuo Emiion: CO, Nox, UHC, Soo Modeling equiemen vay due o diffeence in ime and lengh cale, a well a pocee In geneal, finie-ae ineic i needed o pedic emiion Flamele appoach ill ue ineic! Reduced ineic ucceful fo hea eleae and global dynamic many no wo fo emiion Accuacy in PPM i needed fo eliable pedicion Compuaional co fo finie-ae! Soo phyic i elaively unnown

Emiion nea LBO in DOE-HAT Combuo 40 NOx a 5% exce O2, ppmv 30 20 0 NO 200 pi 280 pi 00 pi Emiion Pobe 0 0.4 0.45 0.5 0.55 0.6 0.65 Equivalence aio NO T inle 700K P inle 3.8 am S 0.6 0.4 < φ < 0.53

DOE-HAT Seup and Condiion Tinle 700K,Pinle 3.8 am S 0.6, 0.4 < φ < 0.53 * 85 x 75 x 8 cylindical gid * 85 x 24 x 24 inne Caeian gid * O(2-4) inepolaion LES-LEM only in he flame zone eolve he flame * Load balancing o achieve peedup

Simplified appoach o pedic emiion Polluan (CO, NO, UHC) aced a he LES level ~ Y x ~ Y ~ Y m ( D D ) w! m ~ u ~ m j m T j x j Slow chemiy, eacion ae obained fom CHEMKIN G-equaion appoach ued o ac flame in LES Hea eleae in enegy equaion a a hin zone Tubulen flame peed model in he LES G-equaion S S ( uʹ, S ), uʹ obainedfom LDKM T T L LES-LEM appoach Global finie-ae ineic ued in he ubgid o obain lamina flame peed and flame ucue Tubulen flame peed acually pediced

CO Pedicion Thee mechanim modeled CO poducion a he flame fon Teaed a a jump diconinuiy Rae obained uing CHEMKIN Equilibium beween CO oxidaion and CO2 diociaion CO O CO 2 2 Fowad/bacwad ae obained fom CHEMKIN CO poducion via UHC oxidaion UHC fomed due o local flame quenching UHC oxidaion o CO modeled a an Aheniu ae O

UHC Pedicion Local quenching of flame due o ech effec (Meneveau and Poino, 9) Unbun fuel eleaed on he poduc ide Efficiency coefficien deemine poion of quenched flame Efficiency Eddy Size Size oo lage o peneae ino he eacion zone Size oo mall o uvive long enough o quench he flame

NO Pedicion Two mechanim included Fomaion a he flame fon Obained fom CHEMKIN Fomaion via he Zeldovich mechanim N N 2 O O 2 N O2 2 NO NO O and N aumed o be in equilibium N O 2NO

LES-G veu LES-LEM Reoluion Iue Eddie lage han flame hicne eolved in LES-G and LES-LEM LES-G baely eolve flame hicne while LES-LEM ha aound 2 cell wihin flame hea eleae implemened in enegy equaion a a hin-zone Eddie of ize of flame pehea zone ae eolved in LES-LEM Flame boadening effec included in LES-G via a model Eddie of ize of flame eacion zone ae paially eolved in LEM No eolved in LES-G δ RR L Δ LES Δ LEM η δ F

Emiion pedicion: UHC (φ 0.4 ) Flame uface UHC io-uface Conou line ae CO ma facion UHC poducion localized in egion of high hea Oue bounday laye Flame lif-off Combuion in he diibued egime.

Emiion pedicion: CO φ 0.4 φ 0.53 CO ma facion (ppm) 8000 0 Flame locaion (oange) and CO ma facion (conou line) Low equivalence aio: long flame and low CO oxidaion High equivalence aio: ho flame and fa CO oxidaion

CO Pedicion fo φ 0.4 NO UHC UHC wih fa oxidaion in CO UHC wih low oxidaion in CO Dump plane Poducion of CO a he flame fon CO oxidaion Emiion pobe Poducion of CO via UHC oxidaion CO a equilibium UHC oxidaion ae i eenial o pedic CO emiion accuaely

CO emiion in he DOEHAT Combuo Indicae level Of inabiliy GLES LEMLES A GLES: Model can be uned o mach daa bu wih no phyic LEMLES: No paamee o adju o conol Noe: Boh imulaion employed he ame CO emiion model

How Soo i Fomed? Sep in Soo Fomaion Fomaion of pecuo Paicle Incepion Suface gowh paicle agglomeaion Paicle oxidaion Range of cale 0. 0 nm Spaially and empoally vaying in he domain

Modeled Soo Relaed Pocee Inenal pocee Nucleaion: Soo nuclei incepion by aceylene Coagulaion: Paicle coalece Suface gowh: Ma depoiion on paicle Agglomeaion: Fomaion of lage chain-lie ucue Oxidaion: Deucion by O2 and OH Exenal pocee Radiaion (opically hin model fo abopion by oo, CO 2 H 2 O gae (Kaplan 996) Themophoei Tanpo by Bownian diffuion Ohe unnown pocee

Radiaion Model Cuen Implemenaion Opically Thin model fo abopion by oo and CO 2, H 2 O gae (Kaplan 996) Moe deailed, bu elaively efficien FAST Coelaed appoach unde udy (Dembele and Wen, 2003) Ue 43 pecal band of vaiable widh fo H2O, CO2 and CO inead of many naow band 5 poin G-L quadaue (inead of 7 o 0 poin) Need moe wo o chec i applicabiliy wihin LEMMOM

Soo Kineic - Linded (994) Soo nucleaion Soo uface gowh Soo Oxidaion CH 2C H n 2 2 2 CH C H 2 2 2 2 C OH CO H C O2 CO 2 Baed on aceylene a a oo pecuo Suiable fo ubulen flame, wih low cabon conen fuel (CH 4 - C 2 H 4 )

Mehod of Momen Appoach The paicle ize diibuion (PSD) i unnown in advance Fo polydipee paicle i i vey had o pecify one ype of PSD (Fiedlande, 2000) Howeve, nowing he momen i equivalen o nowing he PSD (Hudon, 963) MOM wih Inepolaive Cloue (MOMIC) developed by Fenlach and Wang ued in LEMLES El-Aag e al. (Comb. Flame 2006, 2007) Ohe mehod being developed (Pich)

LEM-MOM Subgid Model LEM (Diffuion Reacion Hea Tanfe Soo anpo) T Ga pecie C 2 H 2, OH, O 2 Y -- f v Aea oo Soo Opical Pop. (Ligh caeing- Exincion- Radiaion) MOM

Subgid Combuion Model fo Sooing Flame Ti N p p N p p F q W h C x T x C x T Y V C C T.., ω κ ( ) Si c T Ki K F W x V V Y Y N F W x V Y Y ω ω.., i i i N m M M i S G R d dm Whee M i he h Momen of Paicle Size Diibuion (PSD) Funcion R Nucleaion ae G Coagulaion ae S Suface Gowh ae

LEM-MOM Subgid Model Rad F W h C x T x C x T Y V C C T Ti N p p N p p., ω κ Ki c Ki K F W Y N F W x V Y Y ω ω.., i i i N m M Mi F Ox S G R M Whee M i he h Momen of Paicle Size Diibuion Funcion R Nucleaion ae G Coagulaion ae S Suface Gowh ae Ox Oxidaion ae

Soo Popeie Fom MOM Δ M M d M A M M d Y f M Y M N p o o p v o 2 /3 6 ω π π N Soo Numbe Deniy Y Soo Ma Facion f v Soo Volume Facion d p Soo Paicle Diamee A Soo Suface Aea ω Souce Tem fo LEM M 0 Zeo Momen of PSD M Fi Momen of PSD