Course Outline. Boundary Layer Flashback Core Flow Flashback and Combustion Induced Vortex Breakdown

Similar documents
Flame Stretch, Edge Flames, and Flame Stabilization

Lewis number and curvature effects on sound generation by premixed flame annihilation

Burning Rate, Kinetic Coupling, and Mechanism Reduction

THE ACOUSTIC RESPONSE OF BURNER-STABILIZED PREMIXED FLAT FLAMES

The prediction of turbulence intensities in unsteady flow

Applying Laminar and Turbulent Flow and measuring Velocity Profile Using MATLAB

Control-relevant Modelling and Linear Analysis of Instabilities in Oxy-fuel Combustion

ρ u = u. (1) w z will become certain time, and at a certain point in space, the value of

An Experimental Facility for Studying Flame Propagation at Elevated Pressures

STATIC, STAGNATION, AND DYNAMIC PRESSURES

CFD-Simulation thermoakustischer Resonanzeffekte zur Bestimmung der Flammentransferfunktion

IJAET International Journal of Application of Engineering and Technology ISSN: Vol.1 No.1

Measurements of Hydrogen Syngas Flame Speeds. at Elevated Pressures

Experimental Study of an Impinging Round Jet

MEASUREMENT OF TURBULENCE STATISTICS USING HOT WIRE ANEMOMETRY

Turbomachinery Lecture Notes

Nonlinear analytical flame models with amplitude-dependent time-lag distributions

FEA Solution Procedure

Pressure limit of hydrogen spontaneous ignition in a T-shaped channel

Fluid-induced rotordynamic forces produced by the uid in an annular seal or in the leakage

Two identical, flat, square plates are immersed in the flow with velocity U. Compare the drag forces experienced by the SHADED areas.

Large Eddy Simulation Of Flow Past A Two-dimensional Hill

Axial Compressor Design Parameters

DILUTE GAS-LIQUID FLOWS WITH LIQUID FILMS ON WALLS

Interrogative Simulation and Uncertainty Quantification of Multi-Disciplinary Systems

EDEXCEL NATIONAL CERTIFICATE/DIPLOMA. PRINCIPLES AND APPLICATIONS of FLUID MECHANICS UNIT 13 NQF LEVEL 3 OUTCOME 3 - HYDRODYNAMICS

Fluid Dynamics. Type of Flows Continuity Equation Bernoulli Equation Steady Flow Energy Equation Applications of Bernoulli Equation

High speed analysis of high pressure combustion in a constant volume cell

The Influence of Aerodynamic Characteristics of the Elements of the Flow Range of the Radial Two-row Microturbine on its Dynamic Characteristics

4 Exact laminar boundary layer solutions

Wall treatment in Large Eddy Simulation

Discontinuous Fluctuation Distribution for Time-Dependent Problems

5.1 Heat removal by coolant flow

1 JAXA Special Pblication JAXA-SP-1-E Small-scale trblence affects flow fields arond a blff body and therefore it governs characteristics of cross-sec

Evaluation of the Fiberglass-Reinforced Plastics Interfacial Behavior by using Ultrasonic Wave Propagation Method

Consistent Numerical Model for Wind Buffeting Analysis of Long-Span Bridges

Homotopy Perturbation Method for Solving Linear Boundary Value Problems

Modeling the transition to turbulence in shear flows

UNIT V BOUNDARY LAYER INTRODUCTION

Flood flow at the confluence of compound river channels

School of Aerospace Engineering. Course Outline

Lecture 15. The Turbulent Burning Velocity

Applicability Limits of Operational Modal Analysis to Operational Wind Turbines

ON THE SIMILARITIES OF THE ENGINEERING AND ATMOSPHERIC BOUNDARY LAYERS

Prof. Byoung-Kwon Ahn. College of Engineering, Chungnam National University. flow in pipes and the analysis of fully developed flow.

Spring Semester 2011 April 5, 2011

The Bow Shock and the Magnetosheath

Effects of MHD Laminar Flow Between a Fixed Impermeable Disk and a Porous Rotating Disk

Lecture Notes: Finite Element Analysis, J.E. Akin, Rice University

MODELLING OF TURBULENT ENERGY FLUX IN CANONICAL SHOCK-TURBULENCE INTERACTION

Appendix A: The Fully Developed Velocity Profile for Turbulent Duct Flows

Chapter 9 Flow over Immersed Bodies

called the potential flow, and function φ is called the velocity potential.

FLUCTUATING WIND VELOCITY CHARACTERISTICS OF THE WAKE OF A CONICAL HILL THAT CAUSE LARGE HORIZONTAL RESPONSE OF A CANTILEVER MODEL

ON THE LIMITS OF THE BETZ S EFFICIENCY COEFFICIENT OF WIND TURBINES WITH THE HORIZONTALLY SHAFT

Inertial Instability of Arbitrarily Meandering Currents Governed by the Eccentrically Cyclogeostrophic Equation

WEAR PREDICTION OF A TOTAL KNEE PROSTHESIS TIBIAL TRAY

4 Primitive Equations

Kragujevac J. Sci. 34 (2012) UDC 532.5: :537.63

Entropy ISSN

Primary dependent variable is fluid velocity vector V = V ( r ); where r is the position vector

Large Eddy Simulation of Flame Flashback by Combustion Induced Vortex Breakdown

Lateral Load Capacity of Piles

Quantitative Characterization of Pressure-related Turbulence Transport Terms using Simultaneous Nonintrusive Pressure and Velocity Measurement*

NATURAL CONVECTION No mechanical force to push the fluid pump, fan etc. No predefined fluid flowrate and velocity can t prescribe Reynolds

arxiv: v1 [physics.flu-dyn] 4 Sep 2013

5. The Bernoulli Equation

THE REDUCTION IN FINESTRUCTURE CONTAMINATION OF INTERNAL WAVE ESTIMATES FROM A TOWED THERMISTOR CHAIN

A Decomposition Method for Volume Flux. and Average Velocity of Thin Film Flow. of a Third Grade Fluid Down an Inclined Plane

Latent Differential Equation Modeling with Multivariate Multi-Occasion Indicators

SOIL NON-LINEAR BEHAVIOR AND HYSTERETIC DAMPING IN THE SPRING-DASHPOT ANALOG

Transient Approach to Radiative Heat Transfer Free Convection Flow with Ramped Wall Temperature

Computational Fluid Dynamics Simulation and Wind Tunnel Testing on Microlight Model

FUZZY BOUNDARY ELEMENT METHODS: A NEW MULTI-SCALE PERTURBATION APPROACH FOR SYSTEMS WITH FUZZY PARAMETERS

STUDY OF AC ELECTROOSMOTIC FLOW DEVELOPED BY CO-PLANAR MICROELECTRODE ARRAY IN A SLIT MICROCHANNEL

STUDY OF THE NON-DIMENSIONAL SOLUTION OF DYNAMIC EQUATION OF MOVEMENT ON THE PLANE PLAQUE WITH CONSIDERATION OF TWO-ORDER SLIDING PHENOMENON

INTERNATIONAL JOURNAL OF MATHEMATICS AND COMPUTERS IN SIMULATION Volume 10, 2016

IMECE UNSTEADY VISCOUS FLOWS AND STOKES S FIRST PROBLEM

Similarity Solution for MHD Flow of Non-Newtonian Fluids

A Single Species in One Spatial Dimension

SIMULATION OF PRECESSION IN AXISYMMETRIC SUDDEN EXPANSION FLOWS

EXERCISES WAVE EQUATION. In Problems 1 and 2 solve the heat equation (1) subject to the given conditions. Assume a rod of length L.

Analysis of lift-off height and structure of n-heptane tribrachial flames in laminar jet configuration

Modeling Effort on Chamber Clearing for IFE Liquid Chambers at UCLA

Two-media boundary layer on a flat plate

Two Phase Flow Analysis in Electro-Chemical Machining using CFD

UNCERTAINTY FOCUSED STRENGTH ANALYSIS MODEL

Numerical Simulation of Three Dimensional Flow in Water Tank of Marine Fish Larvae

Prandl established a universal velocity profile for flow parallel to the bed given by

Variability sustained pattern formation in subexcitable media

Incompressible Viscoelastic Flow of a Generalised Oldroyed-B Fluid through Porous Medium between Two Infinite Parallel Plates in a Rotating System

Theoretical study on impingement heat transfer with single-phase free-surface slot jets

PLANETARY ORBITS. According to MATTER (Re-examined) Nainan K. Varghese,

Momentum Equation. Necessary because body is not made up of a fixed assembly of particles Its volume is the same however Imaginary

Evolution and transition mechanisms of internal swirling flows with tangential entry

Fundamentals of magnetohydrodynamics

This is a publisher-deposited version published in : Eprints ID : 18469

MAE 320 Thermodynamics HW 4 Assignment

Chem 4501 Introduction to Thermodynamics, 3 Credits Kinetics, and Statistical Mechanics. Fall Semester Homework Problem Set Number 10 Solutions

Process Modeling of Wellbore Leakage for GCS Risk Assessment

Transcription:

Corse Otline A) Introdction and Otlook B) Flame Aerodynamics and Flashback C) Flame Stretch, Edge Flames, and Flame Stabilization Concepts D) Distrbance Propagation and Generation in Reacting Flows E) Flame Response to Harmonic Excitation Bondary Layer Flashback Core Flow Flashback and Combstion Indced Vortex Breakdown Copyright T. Liewen 2013 Unathorized Reprodction Prohibited 1

Flashback and Flameholding Flashback: Upstream propagation of a premixed flame into a region not designed for the flame to exist Occrs when the laminar and/or trblent flame speed exceeds the local flow velocity Reference flow speed and brning velocity? Flameholding: Flame stabilizes in an ndesired region of the combstor after a flashback/atoignition event Problem has hysteretic elements Wall temperatre effects Bondary layer and swirl flow stability effects Copyright T. Liewen 2013 Unathorized Reprodction Prohibited 2

Flashback and Flameholding Mechanisms Flashback in the bondary layer Flame propagation into core flow We ll focs on swirl flows Combstion instabilities Strong acostic plsations lead to nearly reverse flow Note: p /p~ /c=m / i.e. /=(1/M)p /p Show video Heeger et al. Exp. In Flids 2010 Significance of above mechanisms is a strong fnction of: Fel composition Operating conditions Flid mechanics Kröner et al. CST 2007 Copyright T. Liewen 2013 Unathorized Reprodction Prohibited 3

Corse Otline A) Introdction and Otlook B) Flame Aerodynamics and Flashback C) Flame Stretch, Edge Flames, and Flame Stabilization Concepts D) Distrbance Propagation and Generation in Reacting Flows E) Flame Response to Harmonic Excitation Bondary Layer flashback Core Flow Flashback and Combstion Indced Vortex Breakdown Copyright T. Liewen 2013 Unathorized Reprodction Prohibited 4

Bondary Layer Flashback-Classical Treatment Neglects effects of Heat release (changes approach flow) Stretch (changes brning velocity) Flashback occrs if flame speed exceeds flow velocity at distance,, from the wall δ q ( y = δ ) = s ( y = δ ) x q d q Expanding velocity in a Taylor series, establish flashback condition: gδ x q ( y = δ ) δ = 1 y sd x q q g Assming, δq δf, define flashback Karlovitz nmber gδ F Ka = s d x δ q x y y = 0 flame Copyright T. Liewen 2013 Unathorized Reprodction Prohibited 5

Bondary Layer Flashback Flashback Karlovitz nmber approach is well validated for open flames, sch as Bnsen brners Performed detailed kinetics calclations to determine flame speed and thickness for several data sets Shows how prior brning velocity, flame thickness tendencies can be sed to nderstand tendencies Pressre Preheat temperatre Stoichiometry g (1/s) CH 4 -air CH 4 -air (solid) and C 3 H 8 -air (dashed) Data for figres obtained from: Grmer Ind. & Eng. Chem. 1954 Dgger Ind. & Eng. Chem 1955 Copyright T. Liewen 2013 Unathorized Reprodction Prohibited 6

Bondary Layer Flashback Trblent Bondary Layers Mlti-zoned Near wall laminar sblayer, Basic scaling developed for laminar flows holds if: Most literatre data shows g, trblent 3 g, laminar Significant space-time variation dring flashback Images sggest flame interactions with bondary layer instabilities C. Eichler Exp. In Flids 2012 Show video Copyright T. Liewen 2013 Unathorized Reprodction Prohibited 7

Copled Effects of Flame Crvatre and Gas Expansion Flame blging into reactants Approach flow decelerates Streamlines diverge Adverse pressre gradient Implications: Bondary layers adverse pressre gradients lead to separation Swirl flows adverse pressre gradients can lead to vortex breakdown Triple flames flame can propagate into region with velocity that is higher than flame speed Flame stability flame spontaneosly develops wrinkles Copyright T. Liewen 2013 Unathorized Reprodction Prohibited b x,0 x,0 A ξ 2π k Solid thick contors: positive pressre flctations 0.8 0.6 0.4 0.2 0-0.2-0.4 Pressre Velocity -0.6-1 -0.8-0.6-0.4-0.2 0 kx 8

Heat Condction Inflences on Bondary Layers 4 3 2 Reacting 1 Nonreacting Important implications for Scaling velocity gradients in shear layers Flame stretch rates 0 0 0.2 0.4 0.6 0.8 1 max Shear layer instability freqencies acostic sensitivities Copyright T. Liewen 2013 Unathorized Reprodction Prohibited 9

Heat Release and Stretch Effects Heat release modifies approach flow Stretch modifies brning velocity C. Eichler Trbo Expo 2011 Inner lip Unconfined Confined Copyright T. Liewen 2013 Unathorized Reprodction Prohibited 10

Heat Release and Stretch Effects Particlarly important in explaining flameholding phenomenon Once a flashback event has occrred, difficlt to expel flame from combstor Leading point of advancing flashback event sbject to positive crvatre ( ) 5 10 1/ s g 1.5 0 0 01 confined 0 0.5 0 0 nconfined 0 0.25 0.5 0.75 1 φ Effect of gas expansion de to heat release on local flow velocity Copyright T. Liewen 2013 Unathorized Reprodction Prohibited C. Eichler Trbo Expo 2011 11

Stretch Effects Leading point of advancing flashback event sbject to positive crvatre For Ma < 0, this can case: s ( y = δ ) >> s,0 d q d x x y y = 0 flame,0 s d can be a significant nderestimate of flame speed Positively crved flame δ q Copyright T. Liewen 2013 Unathorized Reprodction Prohibited 12

Heat Release Effects Gas expansion across a crved flame alters the approach flow -0.9-0.7-0.5-0.4 Reslting adverse pressre gradient ahead of flame decelerates flow In extreme cases, can case bondary layer separation Approach flow scks flame back into nozzle Flame Separated flow region 10 2 mm Figres: C. Eichler Trbo Expo 2011 Heeger et al. Exp. In Flids 2010 Copyright T. Liewen 2013 Unathorized Reprodction Prohibited Reactants Reactants 13

Corse Otline A) Introdction and Otlook B) Flame Aerodynamics and Flashback C) Flame Stretch, Edge Flames, and Flame Stabilization Concepts D) Distrbance Propagation and Generation in Reacting Flows E) Flame Response to Harmonic Excitation Bondary Layer flashback Core Flow Flashback and Combstion Indced Vortex Breakdown Copyright T. Liewen 2013 Unathorized Reprodction Prohibited 14

Flow Stability and Vortex Breakdown The degree of swirl in the flow, S, has profond inflences on the flow strctre Most prominent featre of high swirl nmber flows is the occrrence of vortex breakdown, which is manifested as a stagnation point followed by reverse flow Stagnation points Billant et al., JFM, 1998 Sarpkaya, JFM, 1971 Copyright T. Liewen 2013 Unathorized Reprodction Prohibited 15

Prominent Featres of Swirling Flows with Vortex Breakdown: Precessing Vortex Core The flow does not instantaneosly rotate abot the geometric centerline The location of zero azimthal velocity is referred to as the precessing vortex core (PVC) The freqency of rotation of the precessing vortex core scales with a Strohal nmber based on axial flow velocity and diameter Leads to a helical pattern in instantaneos axial flow velocity Important to differentiate the PVC from the other helical shear flow strctres which may also be present Copyright T. Liewen 2013 Unathorized Reprodction Prohibited Sorce: Syred, Prog. Energy and Comb. Sci., 2006 16

Prominent Featres of Swirling Flows: Shear Layer Instability Shear layers exist in both span- and streamwise directions Can be axisymmetric or helical Hang and Yang, Proc. Comb. Inst., 2005 Copyright T. Liewen 2013 Unathorized Reprodction Prohibited 17

Flow Stability and Vortex Breakdown Vortex breakdown can be described as a fold catastrophe Bifrcation of the possible steady state soltions to the Navier-Stokes eqations In high Re flows, there is an intermediate swirl nmber range where flow is bi-stable and hysteretic i.e., either vortex breakdown or no vortex breakdown flow state possible Sorce: Lopez, Physics of Flids, 1994 Swirl Nmber Vortex Breakdown Either No Breakdown Vortex Core Size/Nozzle Radis Copyright T. Liewen 2013 Unathorized Reprodction Prohibited 18

Flow Stability and Vortex Breakdown: Example calclation Vortex breakdown can be predicted for given velocity profile x,0 b,0 Q-vortex" velocity profile: 3 2.5 2 1.5 1 0.5 1 0.8 0.6 0.4 0.2 rθ,0 b,0 x,0 2 b,0 b,0 2χ 5 r = 1+ exp( ( ) ) 1 χ 4 r 5 r rθ,0 Sv 4 rc = ( r/ r ) (1 exp( 5 / 4)) c 2 (1 exp( ( ) )) c 0 0 0.5 1 1.5 2 0 rr c Axial and azimthal velocity profiles sed for vortex breakdown calclation, sing S v =0.71 for θ,0 plot. χ = a,0 b,0 + a,0 b,0 Copyright T. Liewen 2013 Unathorized Reprodction Prohibited 19

Flow Stability and Vortex Breakdown: Example Calclation Wake Jet 1.2 1.0 0.8 Breakdown 2.5 2 1.5 Breakdown S v 0.6 S v 1 0.4 0.2 No breakdown 0.5 No breakdown 0 0 0.2 0.4 0.6 0.8 r c /a 1.0 0-0.5-0.25 0 0.25 0.5 χ Following Z. Rsak Copyright T. Liewen 2013 Unathorized Reprodction Prohibited 20

Core Flow Flame Propagation A B Stable Flame U f CIVB r z Swirl Generator Mixing Tbe Combstion Chamber Image reprodced from T. Sattelmayer Vortex breakdown flame interaction Can occr even if flame speed everywhere less than flow speed Gas expansion across a crved flame: 1. Adverse pressre gradient & radial divergence imposed on reactants 2. Low/negative velocity region generated pstream of flame 3. Flame advances frther into reactants 4. Location of vortex breakdown region advances pstream De to bi-stable natre of vortex breakdown bondaries CIVB itself not necessarily bi-stable Copyright T. Liewen 2013 Unathorized Reprodction Prohibited 21