Flame Stretch, Edge Flames, and Flame Stabilization

Size: px
Start display at page:

Download "Flame Stretch, Edge Flames, and Flame Stabilization"

Transcription

1 Flame Stretch, Edge Flames, and Flame Stabilization Natarajan et al. Combstion and Flame 2007 Flame Anchoring -1

2 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 Introdctory Concepts Flame Stretch Edge Flames Flame Stabilization in Shear Layers Flame Stabilization by Stagnation Points Flame Anchoring -2

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 Introdctory Concepts Flame Stretch Edge Flames Flame Stabilization in Shear Layers Flame Stabilization by Stagnation Points Flame Anchoring -3

4 Flame Stabilization and Blowoff Flame stabilization reqires: A point where the local flame speed and flow velocity match: s d = Typically fond in regions with: Low flow velocity» Aerodynamically decelerated regions (VBB) High shear» Locations of flow separation (ISL & OSL) Figres: Natarajan et al. Combstion and Flame 2007 Petersson et al. Applied Optics 2007 Flame Anchoring -4

5 Premixed Flame Stabilization: Basic Effects Flame stabilization: Balance combstion wave propagation with flow velocity Brning velocity, edge speed, atoignition front? Sggests that stable flames are rare However, flames have self-stabilizing mechanisms Shear layer stabilized: Upstream flame propagation increases wall qenching Aerodynamic stabilization velocity profiles ( cm s ) Flow δ q c s d b a x Lewis & Von Elbe 1987 Flame Anchoring -5

6 Flame Anchoring Locations Complex flows can provide mltiple anchor locations II: VBB/ISR I: VBB IV: OSR III: ISR Flame Anchoring -6

7 Flame Stabilization and Blowoff Stabilization locations determine location/spatial distribtion of flame Flame Shape Flame Length Combstor operability, drability, and emissions directly tied to these fndamental characteristics affecting: Heat loadings to combstor hardware Combstion instability bondaries Blowoff limits Flame Anchoring -7

8 Review of the Idealized Premixed Flame Simplest possible premixed flame configration 1-dimensional, planar Adiabatic T0 b = T ad Adiabatic Flame Temperatre ρ = ρ S0 L = ρ b b Mass Brning Flx What are the controlling parameters? Thermal and mass diffsivities Reaction rates Temperatre of reactants Pressre Exothermicity of fel/oxidizer S0 L Fndamental property of fel/oxidizer mixtre Flame Anchoring -8

9 What Happens if a Flame isn t Flat? Flame Anchoring -9

10 What Happens if Flow Field isn t 1-D? Common theme to 2 problems: misaligned convective and diffsive flxes Flame Anchoring -10

11 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 Introdctory Concepts Flame Stretch Edge Flames Flame Stabilization in Shear Layers Flame Stabilization by Stagnation Points Flame Anchoring -11

12 Overview of Flame Stretch Consider a stationary flame; focs on a C.V. intersection of a streamtbe, and the flame. Steady state energy balance (no viscos effects, no body forces) ρ h T = q Enthalpy Convection EnthalpyDiffsion Constittive relation of enthalpy flx (Fickian diffsion, no radiative heat transfer, no Soret or DFor effects) N q = kt T ρ hid i Yi i= 1 Heat Flx Mass Flx Flame Anchoring -12

13 Overview of Flame Stretch Flame stretch effects: is there a net enthalpy loss/gain or change in composition inside the C.V. becase of diffsion flxes throgh its lateral srface? Two mechanisms: Lewis Nmber effects Le =α /D Diffsion of mass and of heat are nbalanced. Differential diffsion effects D Fel /D ox Lighter species diffse faster than heavier species: eqivalence ratio or dilent/reactant ratio inside the C.V. can change. Negative stretch Positive stretch Flame Anchoring -13

14 Lewis nmber effects If Le = 1 then α =D : no net enthalpy loss throgh the lateral srface Energy Eq.: ρ ht = ht ( ρd ) If Le > 1 then α >D : Heat flx > Mass flx o Positive stretch: net enthalpy flx ot of the C.V. o Negative stretch: net enthalpy flx into the C.V. Negative stretch If Le < 1 then α <D : Heat flx < Mass flx o Positive stretch: net enthalpy flx into the C.V. o Negative stretch: net enthalpy flx ot of the C.V. Positive stretch Flame Anchoring -14

15 Differential diffsion effects Consider the temperatre at the tip of a Bnsen Flame T tip : CH 4 /Air: CH 4 is lighter than O 2, ths diffses faster (D Fel >D Ox ) and its concentration inside the C.V. decreases. o If overall eqivalence ratio is lean, then locally in the C.V. φ is made leaner: T tip is lower then T ad calclated at the overall φ. o If overall eqivalence ratio is rich, then locally in the C.V. φ plled toward stoichiometric: T tip is then higher then T ad calclated at the overall φ. C 3 H 8 /Air: C 3 H 8 is heavier than O 2, ths diffses more slowly (D Fel <D Ox ) and its concentration inside the C.V. increases. o The dependence of T tip on the overall φ is opposite to that of CH 4 /Air. Flame Anchoring -15

16 Bnsen Tip Flame Temperatre Data From C.K. Law, Combstion Physics Flame Anchoring -16

17 Example: Tips of Bnsen Flames Propane(C 3 H 8 ) Methane(CH 4 ) d=10mm φ = Ref: Mizomoto, Asaka, Ikai and Law, Proc. Combst. Inst. 20, 1933 (1984) Slide cortesy of J. Seitzman Flame Anchoring -17

18 Stretched Flames: Non-nity Lewis Nmbers Φ=0.27 Φ=0.27 Φ=0.37 T =298K P=1atm Le i >1 Le i 1 Le i <<1 Bell et al. Proceedings of the Combstion Institte (2007) Flame Anchoring -18

19 Mathematical School of Aerospace Engineering expressions of stretch κ κ 1 da A dt Williams (1975) Flame stretch rate is defined as the normalized differential change with respect to time of an infinitesimal flame srface area element Flame stretch rate qantifies the degree of stretch imposed on a differential flame srface element Lagrangian qantity Units of flame stretch are 1/s i.e. an inverse time scale Flame Anchoring -19

20 Mathematical expressions of stretch κ Expression for κ in terms of flow velocity,, and flame sheet velocity, v F : t1 t2 κ = + + ( vf n)( n) = t t + ( vf n)( n) t t 1 2 Hydrodynamic stretch κ a : variation of tangential flow velocity in the tangential direction (t 1, t 2 ) or, eqivalently (by continity), variation of normal flow velocity in the direction normal (n) to the flame. κ a κ b κ = Hydrodynamic stretch 1 da A dt Unsteady Crvatre stretch κ b : non-stationary flames o Positive κ: divergent tangential velocities or expanding flame flame area increases o Negative κ: convergent tangential velocity or contracting flame flame area decreases Unsteady crvatre stretch Flame Anchoring -20

21 Mathematical expressions of stretch κ Stationary flames vf = 0 κ = t t = n n since = n n ( ) ( ) Hydrodynamic stretch can be interpreted as a variation of the angle between flow velocity and flame normal or, eqivalently, as a variation in tangential flow magnitde along the flame srface. Alternative flame stretch expression t κ = nn : + s ( n) = n S n+ s ( n) 1 κ S κ crv T S = ( + ) Flow Strain, s n = ( vf ) 2 o κ s : stretch de to flow non niformities o κ crv : stretch of a crved flame in a niform approach flow κ can be non zero also when the flow strain is zero Flame Anchoring -21

22 Unsteady Effects Motion of Crved Flames Crvatre is present bt flow velocities align with flame srface normal Stationary spherical flame wold be stretchless t1 t2 κ = + + ( vf n)( n) = t t + ( vf n)( n) t t 1 2 κ a κ b Combstion Physics by C.K. Law (Cambridge University Press, 2006) Flame Anchoring -22

23 Weak stretch effects Asymptotic analysis shows that in the linear limit of weak stretch the effect of varios type of stretch (κ, κ a, κ b, κ S, κ crv ) on flame characteristics (s, δ F, ) is the same. For the flame speed measred in a reference frame attached to the nbrned gases, s, we can write: M Definition of Markstein length is not niqe bt depends on the isosrface sed to define it; e.g. for the flame speed measred in a reference frame attached to the brned gasses s b we have: Dimensionless qantities Markstein nmber Ma Karlovitz nmber Ka s,0 s = s ( κ) s κ= 0 + κ= 0κ = s δmκ κ δ Markstein length s s δ κ δ δ b b,0 b b M M M 0 M δfκ s 0,0,0 F s s δ Ma = Ka = = 1 Ma Ka δ Flame Anchoring -23

24 Weak stretch effects The Markstein nmber contains all the stretch effects described previosly Ma= Ma(Le, D Fel /D Ox, φ, ) For lean mixtres of fels lighter than air (e.g. H 2, CH 4 ) and rich mixtres of fels heavier than air (e.g. C 3 H 8 ), Ma < 0 Conversely, for lean mixtres of heavier than air fels and rich mixtres of lighter than air fels, Ma > 0 Ma vales are sensitive to the position at which flow and flame speeds are measred For a specific mixtre, measred Ma nmbers in the literatre can vary widely among different investigators. Flame Anchoring -24

25 Application Stoichiometry Effects C 3 H 8 /Air p = 1 atm, T = 300 K Tseng et al., Comb. Flame, 95(2), 1993 Flame Anchoring -25

26 Application Fel Effects n-alkanes/air p =1atm, T =300K n-c 8 H 18 /air Ref: Tseng et al., Comb. Flame, 95(2), 1993 Ref: Halter et al. Comb and Flame, 157 (2010) Flame Anchoring -26

27 Application Pressre Effects atm 10 atm 15 atm 5 atm κ (1/s) atm 1 atm 10 atm 15 atm Ka=δ κ s 0,0 F H 2 /CO 30/70 (by vol.) in air, T = 300 K Conterflow twin flame φ adjsted at different p to maintain s,0 = 34 cm/s = const Flame Anchoring -27

28 Strong Stretch Effects PREMIXED FLAME CONCEPTS FLAME STRETCH AND FLAME EXTINCTION Flame Anchoring -28

29 Displacement Speed s d and Consmption Speed s c Displacement speed s d : speed at which the flame is moving along its normal relative to the flow The vale of s d depends on the reference srface: Low κ : the approach flow varies weakly pstream of the flame and the iso-srface choice is not to problematic High κ : velocity gradients occrs on a scale comparable to the flame thickness; s d definition becomes ambigos Displacement speed can also become negative when diffsive flxes are strong enogh to contrast the blk convection in the opposite direction Flame Anchoring -29

30 Displacement Speed s d and Consmption Speed s c Consmption speed s c : spatial integral of chemical rates Can be obtained integrating along a streamline the heat prodction rate and normalizing by the total change in sensible enthalpy across the flame Alternatively, the integrated qantity can be a reactant species consmption rate normalized by the total change in reactant species mass density Example: 1D steady flame sensible enthalpy balance (no viscos and body forces) N d( ρh x sens ) x d 0 b b dx = qdx ( ρ hf, iyi x, D, i ) dx ρ sc ρ sc qdx hsens, hsens, dx + dx dx = = i= 1 b b ρ hsens, ρ hsens, = ρ s b b c h sens h sens ρ = ρ (,, ) dq Similarly from species eq. sc = w idx ρ Yi Yi = 0 (,, ) ( ) Flame Anchoring -30

31 Displacement School of Aerospace Engineering Speed s d and Consmption Speed s c s c0 =s d 0 bt their vales differ at non-zero κ vales H 2 /CO 30/70 (by vol.) in air φ = 0.75, T = 300 K, p = 5 atm s (cm/s) s d s c ( ) H 2 s c (enthalpy ) κ (1/s) Flame Anchoring -31

32 Extinction stretch rate κ ext κ ext : maximm stretch that a flame can sstain before extingishing H 2 /CO 30/70 (by vol.) in air φ = 0.75, T = 300K atm 1 atm 10 atm 15 atm Ka=δ κ s 0,0 F Flame Anchoring -32

33 Example: Pressre Effects Most of the available data are for steady symmetric opposed flow flames : o κ ext depends on flame chemical time τ chem =δ F /s,0 (eg. p effects at fixed s,0 ) H 2 /CO 30/70 (by vol.) in air T = 300 K φ adjsted at different p to maintain s,0 = 34 cm/s = const Flame Anchoring -33

34 Example: Fel and Stoichiometry Effects Ref: Jackson et al., Comb. Fl., 1994, 25(1) Flame Anchoring -34

35 Example: Preheat Effects At high diltion/preheating levels, the flame does not "extingish increases in reactant temperatre are eqivalent to a redction in dimensionless activation energy Example: calclation of CH 4 /air flame stagnating against hot prodcts, whose temperatre is indicated on the plot (cm/s) [cm/s] s c K 1400 K 1450 K (cm/s) [cm/s] s d K 1400 K 1450 K Strain Strain Rate rate [1/s] (1/s) Strain rate Rate (1/s) [1/s] Flame Anchoring -35

36 Caveats on Stretch Sensitivity of Highly Stretched flames Stretch sensitivities and κ ext vales are not intrinsic to mixtre bt also depend on manner in which stretch is applied Example: it depends also on flame geometry and configration: o Velocity profile across the flame thickness (eg. κ ext for opposed flow flames depends on jets distance) o Type of stretch (κ, κ a, κ b, κ S, κ crv ) o Length and time scale of flame/stretch interaction H 2 /air, φ = 0.37, T = 298 K, p = 1 atm Flame Anchoring -36

37 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 Introdctory Concepts Flame Stretch Edge Flames Flame Stabilization in Shear Layers Flame Stabilization by Stagnation Points Flame Anchoring -37

38 Overview Real flames have edges strctre is different than continos flames previosly considered: Non-premixed flames do not propagate, bt their edges do; Premixed flame edge velocity is different from the laminar brning rate (e.g., can be negative) Applications Stabilization of non-premixed (a) and premixed flames (b) Propagation of an ignition front (c) Flame propagation after local extinction (d) Edge flame can be advancing, retreating or stationary Attention has to be paid to the observer reference frame. Flame Anchoring -38

39 Edge flame examples Piloted Bnsen flame A retreating flame edge that is stationary in lab coordinate High φ Low φ v flow v F Rajaram et al., Comb. Sci.Tech,(175) 2003 Flame Anchoring -39

40 Edge flame examples Premixed blff body stabilized flame near blow-off Chadhry et al., Comb. Flame (158), 2011 v flow v F Flame Anchoring -40

41 Cabra brner Edge flame examples Dnn et. al, Comb. Flame,(151) 2007 Vitiated coflow 1500K, 0.8m/s 200 Central jet velocity, (m/s) Downstream extinction observed Continos flame Natral gas/air Eqivalence Ratio,φ 0.8 Flame Anchoring -41

42 Edge Flame Concepts Illstrative Model Problem Flame Anchoring -42

43 Bckmaster s Edge flame model problem Generalize the one dimensional nonpremixed chambered flame (z + ). z x Oxidizer T a Y Ox, a, Flame edge Nonpremixed flame T b Y Fel b Fel 2 2 T T T T ρcp kt = k 2 T ρc 2 p x Q w Fel t z x x ρy Transverseflxes = B MW Fel Ea RT YOxe Fel z = 0 v F L Flame Anchoring -43

44 Bckmaster s Edge flame model problem Generalize the one dimensional nonpremixed chambered flame (z + ). z x Oxidizer T a Y Ox, a, Flame edge Nonpremixed flame T b Y Fel b Fel 2 2 T T T T ρcp kt = k 2 T ρc 2 p x Q w Fel t z x x ρy approximate as k T ( T T ) L b 2 = B MW Fel Ea RT YOxe Fel -Approximate transverse flx terms by convective loss-like term -L characterizes the scale of gradients normal to the flame, sch as de to strain. z = 0 v F L Flame Anchoring -44

45 Edge flame model problem Dimensionless eqation (Le=1) E z T Y Q τ = = = = = = R T L T ct τ a Fel, b E, z, T, Q, Da b b p b chem v F 2 dt d T = F( T, Da) 2 dz dz ρcl 2 flow p T 1B k z x Oxidizer T a Y Ox, a, Flame edge z = 0 v F L Nonpremixed flame T b Y Fel b Fel 2 (1 T + Q ) E where F( T, Da) = 1 T + YFel, bda exp and v F = Q T v F ρcl k T p Flame Anchoring -45

46 Soltion Limit for Edge-less Flame Steady state soltion with no z- direction variation: F( T, Da ) = 0 the problem becomes the same as the steady well stirred reactor; recover the same S-crve behavior. 2 (1 T + Q ) E F( T, Da) = 1 T + YFel, bda exp = 0 Q T E Da = 534 Da = 12, 070 I = 14, Q = 3 II We will focs on Da I <Da<Da II range, where propagating flame edges can occr Three possible behaviors: propagates (v F >0) retreats (v F <0) stays stationary (v F =0) Flame Anchoring -46

47 Edge flame model problem Edge flame velocity Thigh 2 dt v F = F( T, Da) dt dz dz T low E = 14, Q = 3 Flame Anchoring -47

48 Edge flame model problem Edge flame velocity T high 2 dt v F = F( T, Da) dt dz dz T low v F sign depends on the nmerator: v F >0 for Da>Da V and v F <0 for Da<Da V E = 14, Q = 3 v = F 0 Da V Physical meaning: Bi-stable steady state region frther divided into 2 regions if an edge exists for Da>Da V the flame edge acts as an ignition sorce; if the flame locally develops an hole this will close; for Da<Da V the nbrned gasses qench the flame; if the flame develops an hole this will spread ntil the whole flame is extingished. Flame Anchoring -48

49 Edge Flame Concepts Edge Strctre and Velocity Flame Anchoring -49

50 Edge flames strctre Non-premixed flames: Advancing (ignition front, v F >0) have often a triple flame strctre; Retreating (extinction front, v F <0) generally consist of a single edge. v F >0 v F ~ 0 v F < 0 Compted Reaction rate contors Dao et al., Proc. Comb. Inst.,(29) 2002 Premixed flames: generally have a single edge bt can have significant hook-like strctres Image conterflow brner Li et al., Comb. Sci. Tech.,(144) 1999 Compted Reaction rate contors Verdarajan et al., Comb. Flame,(114) 1998 Also stretch and heat losses inflence the flame edge strctre Low Stretch High Stretch False color images Cha et al., Comb. Flame,(146) 2006 Flame Anchoring -50

51 Flame edge velocity Gas Expansion Effects Velocity of edge flames v F depends on density ratio σ ρ (=ρ b /ρ ) across the flame, Damkohler nmber Da, heat losses, Le, D Fel /D Ox and Z st. Dependence on σ ρ : for a convex flame the deceleration of the approach flow in front of the flame monotonically increases with σ ρ (to be discssed later) Retsch et al. s nonpremixed flame scaling for Da v F Da s σ ρ Flame Anchoring -51

52 Flame edge velocity Heat Loss Effects Heat losses may be important process in edge flames stabilized near metal srfaces Heat losses can case v F <0 at high Da Example: Non-Premixed conterflow flame Cha, Ronney, Comb. Fl., 2006, 146 (1-2) Flame Anchoring -52

53 Conditions at the flame edge Consider a premixed flame sbject to a spatially varying stretch rate κ at t=0. Hole will form at points where κ >κ ext edges will retreat to points where κ=κ(v F =0)=κ edge κ κ ext κ edge time Spatial coordinate Flame Anchoring -53

54 Conditions at the flame edge From edge flame model problem κ κ e R T (1 T / T ) E 2 edge ad ext 0 ad a For T 0 /T ad =0.2 and E a /R T ad = <κ edge /κ ext <1 κ κ ext κ edge time Spatial coordinate Conclsion: once flame hole opens, it grows larger than wold be expected based on continos flame concepts Flame Anchoring -54

55 Conditions at the flame edge Experiments show additional physics: Tangential flows of hot gases (e.g., in a lab stationary, retreating edge flame) can increase v F and case κ edge >κ ext CH 4 Air, φ = eqilibrim prodcts ( T = 450K p = 1 atm) When a hole forms in a premixed flames, reactants and prodcts can mix: Mass brning rate increases becase of presence of radicals and increased initial temperatre. Peak heat release rate changes little across different diltion levels; The flame looses its S-crve character. Flame Anchoring -55

56 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 Introdctory Concepts Flame Stretch Edge Flames Flame Stabilization in Shear Layers Flame Stabilization by Stagnation Points Flame Anchoring -56

57 Stretch Effects on Shear Layer Flames Flame will extingish when flame stretch rate exceeds κ ext As expected, higher flow velocities reslt in flame extinction occrring at higher vales of κ ext Flame Anchoring -57

58 Shear Layer Stabilized Flames In high speed flows, althogh locally low velocities exist within the shear layer, flame extinction typically leads to liftoff/blowoff Limited by the amont of flame stretch which they can withstand before extinction e.g, 50 m/s jet with 1 mm shear layer thickness shear~d z /dx ~ s -1 Mch greater than typical extinction strain rates, κ ext How is flame stretch related to flow strain in a shear layer? Flame Anchoring -58

59 Sorces of Flame Stretch 1. Flame crvatre 2. Unsteadiness in flame and flow 3. Hydrodynamic strain: κ ( δ nn ) = s, strain ij i j x i j Q. Zhang et al. J. Eng. for Gas Trbines & Power 2010 For reference, flid strain rate given by tensor: S ij 1 i = + 2 xj x i j Flame Anchoring -59

60 Flame Stretch De to Flid Strain General expression can be redced by assming 2-D flow and incompressibility (pstream of flame): 2 2 z x κ s = ( nx nz ) nn x z + z z x z Normal Strain Contribtion Shear Strain Contribtion Flame Anchoring -60

61 Flame Stretch De to Flid Strain Shear Strain Contribtion κ x = nn + z s, shear x z x z Flame strain occrs de to variations in tangential velocity Leads to positive stretch Flame Anchoring -61

62 Flame Stretch De to Flid Strain Normal Strain Contribtion κ ( 2 2 n n ) = s, normal x z z Jet flows typically decelerate prodcing normal strain Leads to negative stretch z Flame Anchoring -62

63 Flame Stretch De to Flid Strain For high speed flows: κ z θ 1 nx Also assme: x z << z x 2 3 = 1 + O( θ ) nz = θ + O( θ ) x = nn + z s, shear x z κ s d s, shear & θ x x z z n z Expressions for shear and normal stretch can be simplified as follows: κ n x n θ 2 2 ( n n ) = s, normal x z κ s, normal z z z z Flame Anchoring -63

64 Stretch from Shear Strain κ s, shear θ x z If θ s d z then κ s, shear s d z x z κ ~ s, shear sd δ sh Stretch rate de to shear ~indep. of? increasing shear bt θ bt can inflence shear layer: δ sh -1/2 Shearing flow velocity profile Flame Anchoring -64

65 Stretch from Normal Strain κ s, normal ~ z z ~ L z char Obtain traditional flow time scaling approach Decelerating flow velocity profile Flame Anchoring -65

66 Total Stretch from Strain Total stretch de to strain κ s z z + θ x z Shearing flow velocity profile Opposite signs if decelerating flow which term dominates - θ small? x z even if << x z vs. θ?? z x z x Decelerating flow velocity profile Flame Anchoring -66

67 Flame Stretch Distribtion in Swirl Flame Dimensionless vale of 1= 4,125 1/s Shows dominance of deceleration term in first 10 mm Shear term dominates farther downstream Normalized Strain Rate (1/s) Axial Location (mm) 22 n-n * / r r x Zhang, Q., Shanboge, S., Shreekrishna,O Connor, J., Liewen, T., Strain Characteristics Near the Flame Attachment Point in a Swirling Flow, Combstion Science and Technology, Vol. 83, 2011, x 2 z 2 n-n z * v/ z z z n-n x r *n nz z * / z x z n-n r *n z * v/ r x nz z x κ κ s s x Flame Anchoring -67

68 Piloting or Flow Recirclation Effects Flame stabilization can be enhanced throgh: Pilot flames Recirclation zones Transport hot prodcts to the attachment point of a flame Flame Anchoring -68

69 Diltion/Liftoff Effects At high diltion/preheating levels, the flame does not "extingish" Increases in reactant temperatre are eqivalent to a redction in dimensionless activation energy Example: calclation of CH 4 /air flame stagnating against hot prodcts with indicated temperatre (cm/s) [cm/s] s c K 1400 K 1450 K (cm/s) [cm/s] s d K 1400 K 1450 K Strain Strain Rate rate [1/s] (1/s) Strain rate Rate (1/s) [1/s] Flame Anchoring -69

70 Blowoff of Blff Body Stabilized Flames Stages of blowoff Stage 1: Flame is continos bt marked by local extinction events Stage 2: Changes in wake dynamics as large scale strctres become visible Stage 3: Blowoff of flame Da Approach: Scaling captres onset of flame extinction events Ability to captre blowoff dependent on link between extinction events and blowoff physics Nair J. Prop. Power 2007 Images cortesy of D. R. Noble Flame Anchoring -70

71 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 Introdctory Concepts Flame Stretch Edge Flames Flame Stabilization in Shear Layers Flame Stabilization by Stagnation Points Flame Anchoring -71

72 Flame Anchoring Locations and Flame Shapes in Swirling Flows (c) IRZ & Oter Nozzle (d) Centerbody & Oter Nozzle Flame Anchoring -72

73 Flame Anchoring Locations and Flame Shapes in Annlar Nozzle Geometries Flame stabilizes in front of stagnation point of vortex breakdown bbble Stagnation point apparently precesses, probably also moves p and down i.e., flame anchoring position highly nsteady, in contrast to stabilization at edges/corners Under what circmstances can sch flames exist? Not always observed; flames may blowoff directly withot reverting to a free floating configration Flow mst have interior stagnation point From Kmar and Liewen Flame Anchoring -73

74 Flame Anchoring -74

75 Vortex Breakdown in Annlar Geometries Natre of centerbody wake/ VBB changes with geometry, swirl #, and Reynolds # Recirclation zone with vortex tbe Recirclation zone with bbble-like breakdown above Merged recirclation zones Sheen et al., Phys. Flids, 1996 Swirl nmber/ Centerbody Diameter Flame Anchoring -75

76 Vortex Breakdown in Annlar Geometries Natre of centerbody wake/ VBB changes with geometry, swirl #, heat release parameter, and Reynolds # Flame Anchoring -76

77 Vortex Breakdown in Annlar Geometries Natre of centerbody wake/ VBB changes with geometry, swirl #, and Reynolds # Flame Anchoring -77

78 Flame Stabilization Inflenced by Downstream Bondary Conditions Flid rotation introdces inertial wave propagation mechanism sb- and spercritical flow distinction Exit bondary condition has significant inflence on vortex breakdown bbble topology Emara et. al 2009 Flame Anchoring -78

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

Course Outline. Boundary Layer Flashback Core Flow Flashback and Combustion Induced Vortex Breakdown 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)

More information

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

Lewis number and curvature effects on sound generation by premixed flame annihilation Center for Trblence Research Proceedings of the Smmer Program 2 28 Lewis nmber and crvatre effects on sond generation by premixed flame annihilation By M. Talei, M. J. Brear AND E. R. Hawkes A nmerical

More information

Measurements of Hydrogen Syngas Flame Speeds. at Elevated Pressures

Measurements of Hydrogen Syngas Flame Speeds. at Elevated Pressures Paper # A16 Topic: Laminar Flames 5 th US Combstion Meeting Organized by the Western States Section of the Combstion Institte and Hosted by the University of California at San Diego March 25-28, 27. Measrements

More information

Burning Rate, Kinetic Coupling, and Mechanism Reduction

Burning Rate, Kinetic Coupling, and Mechanism Reduction Brning Rate, Kinetic Copling, and Mechanism Redction Yigang J Mechanical and Aerospace Engineering Princeton University 27 AFOSR MURI Kick-Off Meeting Generation of Comprehensive Srrogate Kinetic Models

More information

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

ρ u = u. (1) w z will become certain time, and at a certain point in space, the value of THE CONDITIONS NECESSARY FOR DISCONTINUOUS MOTION IN GASES G I Taylor Proceedings of the Royal Society A vol LXXXIV (90) pp 37-377 The possibility of the propagation of a srface of discontinity in a gas

More information

Second-Order Wave Equation

Second-Order Wave Equation Second-Order Wave Eqation A. Salih Department of Aerospace Engineering Indian Institte of Space Science and Technology, Thirvananthapram 3 December 016 1 Introdction The classical wave eqation is a second-order

More information

UNIT V BOUNDARY LAYER INTRODUCTION

UNIT V BOUNDARY LAYER INTRODUCTION UNIT V BOUNDARY LAYER INTRODUCTION The variation of velocity from zero to free-stream velocity in the direction normal to the bondary takes place in a narrow region in the vicinity of solid bondary. This

More information

Radiation Effects on Heat and Mass Transfer over an Exponentially Accelerated Infinite Vertical Plate with Chemical Reaction

Radiation Effects on Heat and Mass Transfer over an Exponentially Accelerated Infinite Vertical Plate with Chemical Reaction Radiation Effects on Heat and Mass Transfer over an Exponentially Accelerated Infinite Vertical Plate with Chemical Reaction A. Ahmed, M. N.Sarki, M. Ahmad Abstract In this paper the stdy of nsteady flow

More information

4 Exact laminar boundary layer solutions

4 Exact laminar boundary layer solutions 4 Eact laminar bondary layer soltions 4.1 Bondary layer on a flat plate (Blasis 1908 In Sec. 3, we derived the bondary layer eqations for 2D incompressible flow of constant viscosity past a weakly crved

More information

Study on Combustion Characteristics of the Blast Furnace Gas in the Constant Volume Combustion Bomb

Study on Combustion Characteristics of the Blast Furnace Gas in the Constant Volume Combustion Bomb Proceedings of the 7th WEA International Conference on Power ystems, Beijing, China, eptemer 15-17, 2007 115 tdy on Comstion Characteristics of the Blast Frnace Gas in the Constant Volme Comstion Bom IU

More information

An Experimental Facility for Studying Flame Propagation at Elevated Pressures

An Experimental Facility for Studying Flame Propagation at Elevated Pressures Paper # 7LT-21 Topic: Laminar and Trblent Flames 8 th U. S. National Combstion Meeting Organized by the Western States Section of the Combstion Institte and hosted by the University of Utah May 19-22,

More information

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

Fluid Dynamics. Type of Flows Continuity Equation Bernoulli Equation Steady Flow Energy Equation Applications of Bernoulli Equation Tye of Flows Continity Eqation Bernolli Eqation Steady Flow Energy Eqation Alications of Bernolli Eqation Flid Dynamics Streamlines Lines having the direction of the flid velocity Flids cannot cross a

More information

Turbomachinery Lecture Notes

Turbomachinery Lecture Notes Trbomachinery Lectre Notes KTH Corse MJ49/MJ41 Trbomachinery for Compressible Flids DRAFT version Compressors Damian M. Vogt KTH Heat and Power Technology Trbomachinery Lectre Notes 1 007-09-19 Axial Compressor

More information

Experimental Study of an Impinging Round Jet

Experimental Study of an Impinging Round Jet Marie Crie ay Final Report : Experimental dy of an Impinging Rond Jet BOURDETTE Vincent Ph.D stdent at the Rovira i Virgili University (URV), Mechanical Engineering Department. Work carried ot dring a

More information

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

IJAET International Journal of Application of Engineering and Technology ISSN: Vol.1 No.1 IJAET International Jornal of Application of Engineering and Technology ISSN: 395-3594 Vol1 No1 ANALYSIS OF SUPERSONIC FLOWS IN THE E -LAVAL NOZZLE AT 1 INTO A SUENLY EXPANE UCT AT L/=WITH CAVITY ASPECT

More information

THE ACOUSTIC RESPONSE OF BURNER-STABILIZED PREMIXED FLAT FLAMES

THE ACOUSTIC RESPONSE OF BURNER-STABILIZED PREMIXED FLAT FLAMES Proceedings of the Combstion Institte, Volme 29, 2002/pp. 115 122 THE ACOUSTIC RESPONSE OF BURNER-STABILIZED PREMIXED FLAT FLAMES K. R. A. M. SCHREEL, R. ROOK and L. P. H. DE GOEY Eindhoven University

More information

Discontinuous Fluctuation Distribution for Time-Dependent Problems

Discontinuous Fluctuation Distribution for Time-Dependent Problems Discontinos Flctation Distribtion for Time-Dependent Problems Matthew Hbbard School of Compting, University of Leeds, Leeds, LS2 9JT, UK meh@comp.leeds.ac.k Introdction For some years now, the flctation

More information

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

Lecture Notes: Finite Element Analysis, J.E. Akin, Rice University 9. TRUSS ANALYSIS... 1 9.1 PLANAR TRUSS... 1 9. SPACE TRUSS... 11 9.3 SUMMARY... 1 9.4 EXERCISES... 15 9. Trss analysis 9.1 Planar trss: The differential eqation for the eqilibrim of an elastic bar (above)

More information

Axial Compressor Design Parameters

Axial Compressor Design Parameters Trbomachinery Lectre Notes 007-09-9 Axial Compressor Design Parameters Damian Vogt Corse MJ49 Nomenclatre Sbscripts Symbol Denotation Unit c Absolte velocity m/s h Enthalpy J/kg m& Mass flow rate kg/s

More information

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

Momentum Equation. Necessary because body is not made up of a fixed assembly of particles Its volume is the same however Imaginary Momentm Eqation Interest in the momentm eqation: Qantification of proplsion rates esign strctres for power generation esign of pipeline systems to withstand forces at bends and other places where the flow

More information

STATIC, STAGNATION, AND DYNAMIC PRESSURES

STATIC, STAGNATION, AND DYNAMIC PRESSURES STATIC, STAGNATION, AND DYNAMIC PRESSURES Bernolli eqation is g constant In this eqation is called static ressre, becase it is the ressre that wold be measred by an instrment that is static with resect

More information

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

Kragujevac J. Sci. 34 (2012) UDC 532.5: :537.63 5 Kragjevac J. Sci. 34 () 5-. UDC 53.5: 536.4:537.63 UNSTEADY MHD FLOW AND HEAT TRANSFER BETWEEN PARALLEL POROUS PLATES WITH EXPONENTIAL DECAYING PRESSURE GRADIENT Hazem A. Attia and Mostafa A. M. Abdeen

More information

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

High speed analysis of high pressure combustion in a constant volume cell High speed analysis of high pressre combstion in a constant volme cell.j.m. Frijters *, R.J.H. Klein-Dowel, S.S. Manski, L.M.T. Somers and R.S.G. Baert Section Combstion Technology Eindhoven University

More information

MEASUREMENT OF TURBULENCE STATISTICS USING HOT WIRE ANEMOMETRY

MEASUREMENT OF TURBULENCE STATISTICS USING HOT WIRE ANEMOMETRY MEASUREMENT OF TURBULENCE STATISTICS USING HOT WIRE ANEMOMETRY Mrgan Thangadrai +, Atl Kmar Son *, Mritynjay Singh +, Sbhendra *, Vinoth Kmar ++, Ram Pyare Singh +, Pradip K Chatterjee + + Thermal Engineering,

More information

Microscale physics of fluid flows

Microscale physics of fluid flows Microscale physics of flid flows By Nishanth Dongari Senior Undergradate Department of Mechanical Engineering Indian Institte of Technology, Bombay Spervised by Dr. Sman Chakraborty Ot line What is microflidics

More information

MODELLING OF TURBULENT ENERGY FLUX IN CANONICAL SHOCK-TURBULENCE INTERACTION

MODELLING OF TURBULENT ENERGY FLUX IN CANONICAL SHOCK-TURBULENCE INTERACTION MODELLING OF TURBULENT ENERGY FLUX IN CANONICAL SHOCK-TURBULENCE INTERACTION Rssell Qadros, Krishnend Sinha Department of Aerospace Engineering Indian Institte of Technology Bombay Mmbai, India 476 Johan

More information

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

Transient Approach to Radiative Heat Transfer Free Convection Flow with Ramped Wall Temperature Jornal of Applied Flid Mechanics, Vol. 5, No., pp. 9-1, 1. Available online at www.jafmonline.net, ISSN 175-57, EISSN 175-645. Transient Approach to Radiative Heat Transfer Free Convection Flow with Ramped

More information

5.1 Heat removal by coolant flow

5.1 Heat removal by coolant flow 5. Convective Heat Transfer 5.1 Heat removal by coolant flow Fel pellet Bond layer Cladding tbe Heat is transferred from the srfaces of the fel rods to the coolant. T Temperatre at center of fc fel pellet

More information

Burgers Equation. A. Salih. Department of Aerospace Engineering Indian Institute of Space Science and Technology, Thiruvananthapuram 18 February 2016

Burgers Equation. A. Salih. Department of Aerospace Engineering Indian Institute of Space Science and Technology, Thiruvananthapuram 18 February 2016 Brgers Eqation A. Salih Department of Aerospace Engineering Indian Institte of Space Science and Technology, Thirvananthapram 18 Febrary 216 1 The Brgers Eqation Brgers eqation is obtained as a reslt of

More information

5. The Bernoulli Equation

5. The Bernoulli Equation 5. The Bernolli Eqation [This material relates predominantly to modles ELP034, ELP035] 5. Work and Energy 5. Bernolli s Eqation 5.3 An example of the se of Bernolli s eqation 5.4 Pressre head, velocity

More information

Well Stirred Reactor Stabilization of flames

Well Stirred Reactor Stabilization of flames Well Stirred Reactor Stabilization of flames Well Stirred Reactor (see books on Combustion ) Stabilization of flames in high speed flows (see books on Combustion ) Stabilization of flames Although the

More information

The Bow Shock and the Magnetosheath

The Bow Shock and the Magnetosheath Chapter 6 The Bow Shock and the Magnetosheath The solar wind plasma travels sally at speeds which are faster than any flid plasma wave relative to the magnetosphere. Therefore a standing shock wave forms

More information

DILUTE GAS-LIQUID FLOWS WITH LIQUID FILMS ON WALLS

DILUTE GAS-LIQUID FLOWS WITH LIQUID FILMS ON WALLS Forth International Conference on CFD in the Oil and Gas, Metallrgical & Process Indstries SINTEF / NTNU Trondheim, Noray 6-8 Jne 005 DILUTE GAS-LIQUID FLOWS WITH LIQUID FILMS ON WALLS John MORUD 1 1 SINTEF

More information

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

EDEXCEL NATIONAL CERTIFICATE/DIPLOMA. PRINCIPLES AND APPLICATIONS of FLUID MECHANICS UNIT 13 NQF LEVEL 3 OUTCOME 3 - HYDRODYNAMICS EDEXCEL NATIONAL CERTIFICATE/DIPLOMA PRINCIPLES AND APPLICATIONS of FLUID MECHANICS UNIT 3 NQF LEVEL 3 OUTCOME 3 - HYDRODYNAMICS TUTORIAL - PIPE FLOW CONTENT Be able to determine the parameters of pipeline

More information

CFD-Simulation thermoakustischer Resonanzeffekte zur Bestimmung der Flammentransferfunktion

CFD-Simulation thermoakustischer Resonanzeffekte zur Bestimmung der Flammentransferfunktion CFD-Simlation thermoakstischer Resonanzeffekte zr Bestimmng der Flammentransferfnktion Ator: Dennis Paschke Technische Universität Berlin Institt für Strömngsmechanik nd Technische Akstik FG Experimentelle

More information

Study of the diffusion operator by the SPH method

Study of the diffusion operator by the SPH method IOSR Jornal of Mechanical and Civil Engineering (IOSR-JMCE) e-issn: 2278-684,p-ISSN: 2320-334X, Volme, Isse 5 Ver. I (Sep- Oct. 204), PP 96-0 Stdy of the diffsion operator by the SPH method Abdelabbar.Nait

More information

The prediction of turbulence intensities in unsteady flow

The prediction of turbulence intensities in unsteady flow University of Wollongong Research Online Faclty of Engineering and Information Sciences - Papers: Part A Faclty of Engineering and Information Sciences 24 The prediction of trblence intensities in nsteady

More information

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

Control-relevant Modelling and Linear Analysis of Instabilities in Oxy-fuel Combustion Control-relevant Modelling and Linear Analysis of Instabilities in Oxy-fel Combstion Dagfinn Snarheim, Lars Imsland and Bjarne A Foss Abstract Semi-closed gas trbine cycles based on oxyfel combstion have

More information

Applying Laminar and Turbulent Flow and measuring Velocity Profile Using MATLAB

Applying Laminar and Turbulent Flow and measuring Velocity Profile Using MATLAB IOS Jornal of Mathematics (IOS-JM) e-issn: 78-578, p-issn: 319-765X. Volme 13, Isse 6 Ver. II (Nov. - Dec. 17), PP 5-59 www.iosrjornals.org Applying Laminar and Trblent Flow and measring Velocity Profile

More information

Chapter 1: Differential Form of Basic Equations

Chapter 1: Differential Form of Basic Equations MEG 74 Energ and Variational Methods in Mechanics I Brendan J. O Toole, Ph.D. Associate Professor of Mechanical Engineering Howard R. Hghes College of Engineering Universit of Nevada Las Vegas TBE B- (7)

More information

Rich Spray-Flame Propagating through a 2D-Lattice of Alkane Droplets in Air

Rich Spray-Flame Propagating through a 2D-Lattice of Alkane Droplets in Air Rich Spray-Flame Propagating throgh a D-attice of Alkane Droplets in Air Colette Nicoli, Brno Denet, Pierre Haldenwang o cite this version: Colette Nicoli, Brno Denet, Pierre Haldenwang. Rich Spray-Flame

More information

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

arxiv: v1 [physics.flu-dyn] 4 Sep 2013 THE THREE-DIMENSIONAL JUMP CONDITIONS FOR THE STOKES EQUATIONS WITH DISCONTINUOUS VISCOSITY, SINGULAR FORCES, AND AN INCOMPRESSIBLE INTERFACE PRERNA GERA AND DAVID SALAC arxiv:1309.1728v1 physics.fl-dyn]

More information

Two-media boundary layer on a flat plate

Two-media boundary layer on a flat plate Two-media bondary layer on a flat plate Nikolay Ilyich Klyev, Asgat Gatyatovich Gimadiev, Yriy Alekseevich Krykov Samara State University, Samara,, Rssia Samara State Aerospace University named after academician

More information

4 Primitive Equations

4 Primitive Equations 4 Primitive Eqations 4.1 Spherical coordinates 4.1.1 Usefl identities We now introdce the special case of spherical coordinates: (,, r) (longitde, latitde, radial distance from Earth s center), with 0

More information

Modeling Effort on Chamber Clearing for IFE Liquid Chambers at UCLA

Modeling Effort on Chamber Clearing for IFE Liquid Chambers at UCLA Modeling Effort on Chamber Clearing for IFE Liqid Chambers at UCLA Presented by: P. Calderoni own Meeting on IFE Liqid Wall Chamber Dynamics Livermore CA May 5-6 3 Otline his presentation will address

More information

Computational Fluid Dynamics Simulation and Wind Tunnel Testing on Microlight Model

Computational Fluid Dynamics Simulation and Wind Tunnel Testing on Microlight Model Comptational Flid Dynamics Simlation and Wind Tnnel Testing on Microlight Model Iskandar Shah Bin Ishak Department of Aeronatics and Atomotive, Universiti Teknologi Malaysia T.M. Kit Universiti Teknologi

More information

School of Aerospace Engineering. Course Outline

School of Aerospace Engineering. Course Outline Course Outline A) Introduction and Outlook B) Flame Aerodynamics and Flashback C) Flame Stretch, Edge Flames, and Flame Stabilization Concepts D) Disturbance Propagation and Generation in Reacting Flows

More information

FEA Solution Procedure

FEA Solution Procedure EA Soltion Procedre (demonstrated with a -D bar element problem) EA Procedre for Static Analysis. Prepare the E model a. discretize (mesh) the strctre b. prescribe loads c. prescribe spports. Perform calclations

More information

PROBLEMS

PROBLEMS PROBLEMS------------------------------------------------ - 7- Thermodynamic Variables and the Eqation of State 1. Compter (a) the nmber of moles and (b) the nmber of molecles in 1.00 cm of an ideal gas

More information

Setting The K Value And Polarization Mode Of The Delta Undulator

Setting The K Value And Polarization Mode Of The Delta Undulator LCLS-TN-4- Setting The Vale And Polarization Mode Of The Delta Undlator Zachary Wolf, Heinz-Dieter Nhn SLAC September 4, 04 Abstract This note provides the details for setting the longitdinal positions

More information

Thermal balance of a wall with PCM-enhanced thermal insulation

Thermal balance of a wall with PCM-enhanced thermal insulation Thermal balance of a wall with PCM-enhanced thermal inslation E. Kossecka Institte of Fndamental Technological esearch of the Polish Academy of Sciences, Warsaw, Poland J. Kośny Oak idge National aboratory;

More information

Sources of Non Stationarity in the Semivariogram

Sources of Non Stationarity in the Semivariogram Sorces of Non Stationarity in the Semivariogram Migel A. Cba and Oy Leangthong Traditional ncertainty characterization techniqes sch as Simple Kriging or Seqential Gassian Simlation rely on stationary

More information

Material Transport with Air Jet

Material Transport with Air Jet Material Transport with Air Jet Dr. István Patkó Bdapest Tech Doberdó út 6, H-1034 Bdapest, Hngary patko@bmf.h Abstract: In the field of indstry, there are only a very few examples of material transport

More information

Turbulence and boundary layers

Turbulence and boundary layers Trblence and bondary layers Weather and trblence Big whorls hae little whorls which feed on the elocity; and little whorls hae lesser whorls and so on to iscosity Lewis Fry Richardson Momentm eqations

More information

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

Prandl established a universal velocity profile for flow parallel to the bed given by EM 0--00 (Part VI) (g) The nderlayers shold be at least three thicknesses of the W 50 stone, bt never less than 0.3 m (Ahrens 98b). The thickness can be calclated sing Eqation VI-5-9 with a coefficient

More information

Reflections on a mismatched transmission line Reflections.doc (4/1/00) Introduction The transmission line equations are given by

Reflections on a mismatched transmission line Reflections.doc (4/1/00) Introduction The transmission line equations are given by Reflections on a mismatched transmission line Reflections.doc (4/1/00) Introdction The transmission line eqations are given by, I z, t V z t l z t I z, t V z, t c z t (1) (2) Where, c is the per-nit-length

More information

Variability sustained pattern formation in subexcitable media

Variability sustained pattern formation in subexcitable media Variability sstained pattern formation in sbexcitable media Erik Glatt, Martin Gassel, and Friedemann Kaiser Institte of Applied Physics, Darmstadt University of Technology, 64289 Darmstadt, Germany (Dated:

More information

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

Chem 4501 Introduction to Thermodynamics, 3 Credits Kinetics, and Statistical Mechanics. Fall Semester Homework Problem Set Number 10 Solutions Chem 4501 Introdction to Thermodynamics, 3 Credits Kinetics, and Statistical Mechanics Fall Semester 2017 Homework Problem Set Nmber 10 Soltions 1. McQarrie and Simon, 10-4. Paraphrase: Apply Eler s theorem

More information

Chapter 9 Flow over Immersed Bodies

Chapter 9 Flow over Immersed Bodies 57:00 Mechanics o Flids and Transport Processes Chapter 9 Proessor Fred Stern Fall 01 1 Chapter 9 Flow over Immersed Bodies Flid lows are broadly categorized: 1. Internal lows sch as dcts/pipes, trbomachinery,

More information

Curves - Foundation of Free-form Surfaces

Curves - Foundation of Free-form Surfaces Crves - Fondation of Free-form Srfaces Why Not Simply Use a Point Matrix to Represent a Crve? Storage isse and limited resoltion Comptation and transformation Difficlties in calclating the intersections

More information

Homotopy Perturbation Method for Solving Linear Boundary Value Problems

Homotopy Perturbation Method for Solving Linear Boundary Value Problems International Jornal of Crrent Engineering and Technolog E-ISSN 2277 4106, P-ISSN 2347 5161 2016 INPRESSCO, All Rights Reserved Available at http://inpressco.com/categor/ijcet Research Article Homotop

More information

Pulses on a Struck String

Pulses on a Struck String 8.03 at ESG Spplemental Notes Plses on a Strck String These notes investigate specific eamples of transverse motion on a stretched string in cases where the string is at some time ndisplaced, bt with a

More information

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

Appendix A: The Fully Developed Velocity Profile for Turbulent Duct Flows Appendix A: The lly Developed Velocity Profile for Trblent Dct lows This appendix discsses the hydrodynamically flly developed velocity profile for pipe and channel flows. The geometry nder consideration

More information

UNIT IV BOUNDARY LAYER AND FLOW THROUGH PIPES Definition of bonary layer Thickness an classification Displacement an momentm Thickness Development of laminar an trblent flows in circlar pipes Major an

More information

CHEMICAL REACTION EFFECTS ON FLOW PAST AN EXPONENTIALLY ACCELERATED VERTICAL PLATE WITH VARIABLE TEMPERATURE. R. Muthucumaraswamy and V.

CHEMICAL REACTION EFFECTS ON FLOW PAST AN EXPONENTIALLY ACCELERATED VERTICAL PLATE WITH VARIABLE TEMPERATURE. R. Muthucumaraswamy and V. International Jornal of Atomotive and Mechanical Engineering (IJAME) ISSN: 9-8649 (int); ISSN: 18-166 (Online); Volme pp. 31-38 Jly-December 1 niversiti Malaysia Pahang DOI: http://dx.doi.org/1.158/ijame..11.11.19

More information

3. Flows in a Steam Path. Steam Turbine 3. Flow in Steam Path 1 / 112

3. Flows in a Steam Path. Steam Turbine 3. Flow in Steam Path 1 / 112 3. Flows in a Steam Path Steam Trbine 3. Flow in Steam Path / 3 4 5 6 7 Flid Dynamics Steam Trbine Flow Model 4 Thermodynamics and Flid Dynamics for Steam Trbines 35 Dimensionless Nmbers 49 Implse Trbine

More information

A Survey of the Implementation of Numerical Schemes for Linear Advection Equation

A Survey of the Implementation of Numerical Schemes for Linear Advection Equation Advances in Pre Mathematics, 4, 4, 467-479 Pblished Online Agst 4 in SciRes. http://www.scirp.org/jornal/apm http://dx.doi.org/.436/apm.4.485 A Srvey of the Implementation of Nmerical Schemes for Linear

More information

IMECE UNSTEADY VISCOUS FLOWS AND STOKES S FIRST PROBLEM

IMECE UNSTEADY VISCOUS FLOWS AND STOKES S FIRST PROBLEM Proceedings of IMECE 006 006 ASME International Mechanical Engineering Congress Chicago, Illinois, USA, November 5-0, 006 IMECE 006-430 UNSTEADY VISCOUS FLOWS AND STOKES S FIRST PROBLEM Y.S. Mzychka Faclty

More information

Steady State and Transient Thermal Analysis of Switched Reluctance Machine

Steady State and Transient Thermal Analysis of Switched Reluctance Machine Steady State and Transient Thermal Analysis of Switched Relctance Machine E. Annie Elisabeth Jebaseeli and S. Paramasivam Abstract This paper presents the two dimensional (-D) steady state and transient

More information

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

Two identical, flat, square plates are immersed in the flow with velocity U. Compare the drag forces experienced by the SHADED areas. Two identical flat sqare plates are immersed in the flow with velocity U. Compare the drag forces experienced by the SHAE areas. F > F A. A B F > F B. B A C. FA = FB. It depends on whether the bondary

More information

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

Inertial Instability of Arbitrarily Meandering Currents Governed by the Eccentrically Cyclogeostrophic Equation Jornal of Oceanography, Vol. 59, pp. 163 to 17, 3 Inertial Instability of Arbitrarily Meandering Crrents Governed by the Eccentrically Cyclogeostrophic Eqation HIDEO KAWAI* 131-81 Shibagahara, Kse, Joyo,

More information

International Journal of Scientific & Engineering Research, Volume 5, Issue 3, March ISSN

International Journal of Scientific & Engineering Research, Volume 5, Issue 3, March ISSN International Jornal of Scientific & Engineering Research, Volme 5, Isse 3, March-4 83 ISSN 9-558 Doble Dispersion effects on free convection along a vertical Wavy Srface in Poros Media with Variable Properties

More information

TRANSONIC EVAPORATION WAVES IN A SPHERICALLY SYMMETRIC NOZZLE

TRANSONIC EVAPORATION WAVES IN A SPHERICALLY SYMMETRIC NOZZLE TRANSONIC EVAPORATION WAVES IN A SPHERICALLY SYMMETRIC NOZZLE XIAOBIAO LIN AND MARTIN WECHSELBERGER Abstract. This paper stdies the liqid to vapor phase transition in a cone shaped nozzle. Using the geometric

More information

Spring Semester 2011 April 5, 2011

Spring Semester 2011 April 5, 2011 METR 130: Lectre 4 - Reynolds Averaged Conservation Eqations - Trblent Flxes (Definition and typical ABL profiles, CBL and SBL) - Trblence Closre Problem & Parameterization Spring Semester 011 April 5,

More information

Large Eddy Simulation Of Flow Past A Two-dimensional Hill

Large Eddy Simulation Of Flow Past A Two-dimensional Hill Large Eddy Simlation Of Flow Past A Two-dimensional Hill Sankara N.Vengadesan ) and Akihiko Nakayama ) ) Research Associate, Email: vengades@kobe-.ac.jp, ) Professor, Email: nakayama@kobe-.ac.jp Gradate

More information

Reduction of over-determined systems of differential equations

Reduction of over-determined systems of differential equations Redction of oer-determined systems of differential eqations Maim Zaytse 1) 1, ) and Vyachesla Akkerman 1) Nclear Safety Institte, Rssian Academy of Sciences, Moscow, 115191 Rssia ) Department of Mechanical

More information

Laminar Premixed Flames: Flame Structure

Laminar Premixed Flames: Flame Structure Laminar Premixed Flames: Flame Structure Combustion Summer School 2018 Prof. Dr.-Ing. Heinz Pitsch Course Overview Part I: Fundamentals and Laminar Flames Introduction Fundamentals and mass balances of

More information

Lecture 6 Asymptotic Structure for Four-Step Premixed Stoichiometric Methane Flames

Lecture 6 Asymptotic Structure for Four-Step Premixed Stoichiometric Methane Flames Lecture 6 Asymptotic Structure for Four-Step Premixed Stoichiometric Methane Flames 6.-1 Previous lecture: Asymptotic description of premixed flames based on an assumed one-step reaction. basic understanding

More information

Technical Note. ODiSI-B Sensor Strain Gage Factor Uncertainty

Technical Note. ODiSI-B Sensor Strain Gage Factor Uncertainty Technical Note EN-FY160 Revision November 30, 016 ODiSI-B Sensor Strain Gage Factor Uncertainty Abstract Lna has pdated or strain sensor calibration tool to spport NIST-traceable measrements, to compte

More information

AN ISOGEOMETRIC SOLID-SHELL FORMULATION OF THE KOITER METHOD FOR BUCKLING AND INITIAL POST-BUCKLING ANALYSIS OF COMPOSITE SHELLS

AN ISOGEOMETRIC SOLID-SHELL FORMULATION OF THE KOITER METHOD FOR BUCKLING AND INITIAL POST-BUCKLING ANALYSIS OF COMPOSITE SHELLS th Eropean Conference on Comptational Mechanics (ECCM ) 7th Eropean Conference on Comptational Flid Dynamics (ECFD 7) 5 Jne 28, Glasgow, UK AN ISOGEOMETRIC SOLID-SHELL FORMULATION OF THE KOITER METHOD

More information

3 2D Elastostatic Problems in Cartesian Coordinates

3 2D Elastostatic Problems in Cartesian Coordinates D lastostatic Problems in Cartesian Coordinates Two dimensional elastostatic problems are discssed in this Chapter, that is, static problems of either plane stress or plane strain. Cartesian coordinates

More information

Lecture 8 Laminar Diffusion Flames: Diffusion Flamelet Theory

Lecture 8 Laminar Diffusion Flames: Diffusion Flamelet Theory Lecture 8 Laminar Diffusion Flames: Diffusion Flamelet Theory 8.-1 Systems, where fuel and oxidizer enter separately into the combustion chamber. Mixing takes place by convection and diffusion. Only where

More information

Fundamentals of magnetohydrodynamics

Fundamentals of magnetohydrodynamics Fndamentals of magnetohydrodynamics Part II Daniel Gómez 1, Email: dgomez@df.ba.ar Webpage: astro.df.ba.ar (1 Institto de Astronomía y Física del Espacio, CONICET, Argentina ( Departamento de Física, Universidad

More information

Classify by number of ports and examine the possible structures that result. Using only one-port elements, no more than two elements can be assembled.

Classify by number of ports and examine the possible structures that result. Using only one-port elements, no more than two elements can be assembled. Jnction elements in network models. Classify by nmber of ports and examine the possible strctres that reslt. Using only one-port elements, no more than two elements can be assembled. Combining two two-ports

More information

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

Evaluation of the Fiberglass-Reinforced Plastics Interfacial Behavior by using Ultrasonic Wave Propagation Method 17th World Conference on Nondestrctive Testing, 5-8 Oct 008, Shanghai, China Evalation of the Fiberglass-Reinforced Plastics Interfacial Behavior by sing Ultrasonic Wave Propagation Method Jnjie CHANG

More information

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

Pressure limit of hydrogen spontaneous ignition in a T-shaped channel 4 th International Conference on Hydrogen Safety, 12-14 September 2011, San Francisco, USA Pressre limit of hydrogen spontaneos ignition in a T-shaped channel Maim Bragin, Dmitriy Makarov, Vladimir Molkov

More information

Efficiency Increase and Input Power Decrease of Converted Prototype Pump Performance

Efficiency Increase and Input Power Decrease of Converted Prototype Pump Performance International Jornal of Flid Machinery and Systems DOI: http://dx.doi.org/10.593/ijfms.016.9.3.05 Vol. 9, No. 3, Jly-September 016 ISSN (Online): 188-9554 Original Paper Efficiency Increase and Inpt Power

More information

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

STUDY OF AC ELECTROOSMOTIC FLOW DEVELOPED BY CO-PLANAR MICROELECTRODE ARRAY IN A SLIT MICROCHANNEL Proceedings of the International Conference on Mechanical Engineering 211 (ICME211 18-2 ecember 211, haka, Bangladesh ICME 11 225 STUY O AC EECTROOSMOTIC OW EVEOPE BY CO-PANAR MICROEECTROE ARRAY IN A SIT

More information

1. INTRODUCTION. A solution for the dark matter mystery based on Euclidean relativity. Frédéric LASSIAILLE 2009 Page 1 14/05/2010. Frédéric LASSIAILLE

1. INTRODUCTION. A solution for the dark matter mystery based on Euclidean relativity. Frédéric LASSIAILLE 2009 Page 1 14/05/2010. Frédéric LASSIAILLE Frédéric LASSIAILLE 2009 Page 1 14/05/2010 Frédéric LASSIAILLE email: lmimi2003@hotmail.com http://lmi.chez-alice.fr/anglais A soltion for the dark matter mystery based on Eclidean relativity The stdy

More information

Advanced topics in Finite Element Method 3D truss structures. Jerzy Podgórski

Advanced topics in Finite Element Method 3D truss structures. Jerzy Podgórski Advanced topics in Finite Element Method 3D trss strctres Jerzy Podgórski Introdction Althogh 3D trss strctres have been arond for a long time, they have been sed very rarely ntil now. They are difficlt

More information

Design and Data Acquisition for Thermal Conductivity Matric Suction Sensors

Design and Data Acquisition for Thermal Conductivity Matric Suction Sensors 68 TRANSPORTATION RSARCH RCORD 1432 Design and Data Acqisition for Thermal Condctivity Matric Sction Sensors J. K.-M. GAN, D. G. FRDLUND, A. XING, AND W.-X. LI The principles behind sing the thermal condctivity

More information

Study on the impulsive pressure of tank oscillating by force towards multiple degrees of freedom

Study on the impulsive pressure of tank oscillating by force towards multiple degrees of freedom EPJ Web of Conferences 80, 0034 (08) EFM 07 Stdy on the implsive pressre of tank oscillating by force towards mltiple degrees of freedom Shigeyki Hibi,* The ational Defense Academy, Department of Mechanical

More information

An Interpretation of the Black Energy in Universe by Using a Hydro-Dynamical Analogy with Newton Gravity

An Interpretation of the Black Energy in Universe by Using a Hydro-Dynamical Analogy with Newton Gravity An Interpretation of the Black nergy in Universe by Using a Hydro-Dynamical Analogy with Newton Gravity Corneli BRBNT*,1, Sorin BRBNT *Corresponding athor *,1 POLITHNICA University of Bcharest, Faclty

More information

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

Quantitative Characterization of Pressure-related Turbulence Transport Terms using Simultaneous Nonintrusive Pressure and Velocity Measurement* Qantitative Characterization of Pressre-related Trblence Transport Terms sing Simltaneos Nonintrsive Pressre and Velocity Measrement* Xiaofeng Li and Joseph Katz Department of Aerospace Engineering San

More information

DELFT UNIVERSITY OF TECHNOLOGY

DELFT UNIVERSITY OF TECHNOLOGY DELFT UNIVERSITY OF TECHNOLOGY REPORT 02-01 The Krylov accelerated SIMPLER method for incompressible flow C. Vik and A. Saghir ISSN 1389-6520 Reports of the Department of Applied Mathematical Analysis

More information

L = 2 λ 2 = λ (1) In other words, the wavelength of the wave in question equals to the string length,

L = 2 λ 2 = λ (1) In other words, the wavelength of the wave in question equals to the string length, PHY 309 L. Soltions for Problem set # 6. Textbook problem Q.20 at the end of chapter 5: For any standing wave on a string, the distance between neighboring nodes is λ/2, one half of the wavelength. The

More information

Optimal Control of a Heterogeneous Two Server System with Consideration for Power and Performance

Optimal Control of a Heterogeneous Two Server System with Consideration for Power and Performance Optimal Control of a Heterogeneos Two Server System with Consideration for Power and Performance by Jiazheng Li A thesis presented to the University of Waterloo in flfilment of the thesis reqirement for

More information

Experiment and mathematical model for the heat transfer in water around 4 C

Experiment and mathematical model for the heat transfer in water around 4 C Eropean Jornal of Physics PAPER Experiment and mathematical model for the heat transfer in water arond 4 C To cite this article: Naohisa Ogawa and Fmitoshi Kaneko 2017 Er. J. Phys. 38 025102 View the article

More information

1.061 / 1.61 Transport Processes in the Environment

1.061 / 1.61 Transport Processes in the Environment MIT OpenCorseWare http://ocw.mit.ed 1.061 / 1.61 Transport Processes in the Environment Fall 008 For information abot citing these materials or or Terms of Use, visit: http://ocw.mit.ed/terms. Answer 8.1

More information