Stanford University Research Program Shock Tube/Laser Absorption Studies of Chemical Kinetics. Ronald K. Hanson

Transcription

1 Stanford University Research Program Shock Tube/Laser Absorption Studies of Chemical Kinetics Ronald K. Hanson Dept. of Mechanical Engineering, Stanford University Experimental Methods Butanol Kinetics Methyl Ester Kinetics Future Work CEFRC Second Annual Conference August 17 19, 2011 Some of the work presented here is unpublished. Please check with RKH/DFD before regarding the data as final or importing it into publications. We would also value feedback from team members regarding our data and how it might be modeled. 1

2 Please Note Some of the work presented here is unpublished. Please check with RKH/DFD before regarding the data as final or importing it into publications. We would also value feedback from team members regarding our data and how it might be modeled. 2

3 Experiment Types Goal: High quality databases to validate detailed mechanisms Ignition delay times provide global targets Shock tubes provide constant volume data over wide pressure range Species time histories provide strong constraints on mechanisms Laser absorption can provide species time histories for: OH, CO, CO 2, CH 2 O, H 2 O, CH 3, CH 4, C 2 H 4, fuel, Direct determination of elementary reaction rate constants For reactions where estimates/theory are not sufficient 3

4 Kinetics Shock Tube 1 Kinetics Shock Tube 2 Experimental Approach Aerosol Shock Tube High Pressure Shock Tube 4

5 Advances in Shock Tube Methodology Tailored driver gas/new driver geometry provide extended test time for access to low T Driver inserts provide highly uniform reflected shock conditions approaching constant U/V Aerosol shock tube provides access to low vapor pressure fuels Gasdynamic modeling of shock tube flows to account for facility effects and energy release 5

6 Access to Low Temperatures Longer driver length and tailored gas mixtures can provide longer test times (> 40 ms) 2x Driver Extension 7 Shock tubes now can overlap with RCMs 6

7 Improvement in Temperature Uniformity Problem: Ignition delay times are artificially shortened by non ideal facility effects!! dp/dt 0 Solution: Driver Inserts Results: Near ideal constant volume performance!! dp/dt 0 V RS P 5 T 5 7

8 Aerosol Shock Tube for Low Vapor Pressure Fuels Dump Tank Diagnostics: Pressure Droplet scattering Fuel time history OH* emission Aerosol Tank Driver Section Driven Section Evaporated Aerosol Ultrasonic Nebulizers Does not require heated shock tube Eliminates fuel cracking and partial distillation Provides access to low vapor pressure fuels: large methyl esters, bio diesel surrogates 8

9 Current Laser Capabilities for Species Detection for real time, in situ sensing Ultraviolet CH nm NO 225 nm O nm HO nm CH 2 O 305 nm OH 306 nm NH 336 nm Visible CN CH NCO NO 2 NH 2 HCO 388 nm 431 nm 440 nm 472 nm 597 nm 614 nm Infrared CO 2.3 m H 2 O 2.5 m CO m Fuel 3.4 m NO 5.2 m MeOH 9.2 m MF 9.2 m C 2 H m Coherent MIRA Ti Sapphire Spectra Physics 380 Ring NovaWave Mid IR DFG

10 Kinetics Shock Tube 1 Kinetics Shock Tube 2 Butanol Kinetics 10

11 Overview of Butanol Studies Ignition delay times: atm, K Butanol isomers, high and low pressure Species time histories: N Butanol pyrolysis: OH, H 2 O, CH 2 O, C 2 H 4, CH 4, CO N Butanol oxidation: OH, H 2 O, C 2 H 4 Direct determinations of elementary rxn. rate constants: Butanol+OH=products, all isomers Butene+OH= Products, all isomers 11

12 Survey of Ignition Delay Times: Butanol Isomers at Low Pressure 1538 K 1429 K 1333 K 1250 K Variation in ignition delay times - tert-butanol slowest - 1-butanol fastest 1000 Lines - MIT (2011) 1-butanol 2-butanol i-butanol t-butanol t-but 2-but i-but 1-but MIT (2011) mechanism t ign [us] - Fair agreement with t- & i-butanol data 4% O 2 /Argon - Poorer agreement with 1- & 2-butanol data atm, =1.0 Data of sufficient quality to refine reaction mechanisms /T 5 [1/K] What happens at high pressure? 12

13 Survey of Ignition Delay Times: 2 Butanol Variation with Pressure 1538 K 1333 K 1177 K 1053 K Very low scatter ( 5-10 %) consistent with uncertainty in T Lines - MIT (2011) 1.5 atm 3 atm 19atm 43atm Ignition delay times scale approximately as P -0.7 MIT (2011) model P-dependence consistent with 2-butanol data t ign [us] Butanol 4% O 2 /Argon ~1.0 Data of sufficient quality to refine reaction mechanisms /T 5 [1/K] Next step: need for species time histories! 13

14 N Butanol Pyrolysis First shock tube/laser absorption speciation study of n butanol pyrolysis OH (306 nm) H 2 O (2.5 microns) CH 2 O (305 nm) CO (4.6 nm) C 2 H 4 (10.5 microns) CH 4 (3.4 microns) n Butanol 14

15 N Butanol Pyrolysis: Species Time Histories OH and H 2 O OH & H 2 O time histories reveal large variation in model performance Clear opportunity for model refinement 15

16 N Butanol Pyrolysis: Species Time Histories CH 2 O and CO Formaldehyde and CO uniformly underpredicted! Measured OH+H 2 O+CH 2 O+CO account for >90% of O atoms! Remaining O atoms likely in CH 3 CHO, CH 2 CO, 16

17 N Butanol Pyrolysis: Species Time Histories C 2 H 4 and CH 4 All models underpredict C 2 H 4 ; better agreement for CH 4 Measured C 2 H 4 +CH 4 +CO+CH 2 O account for >70% of C atoms! Remaining C atoms likely in C 3 H 6, C 2 H 6, C 2 H 2, 17

18 OH+Butanol Products First order removal of OH measured using laser absorption in 30 ppm TBHP/200ppm butanol/argon mixtures Rate Constant [cc/mole/s] 3E13 1E13 1E K 1000 K 909 K 1E /T [1/K] 1-but iso-but 2-but tert-but Sarathy (2011) Strong dependence on isomer MIT (2011) Model: Trends not consistent with data Sarathy (2011) Model: Consistent with data Overall rate dependent on product channel chemistry 18

19 Kinetics Shock Tube 2 Methyl Ester Kinetics Aerosol Shock Tube 19

20 Overview of Methyl Ester Studies Ignition delay times: Aerosol Shock Tube Studies Methyl Decanoate Methyl Oleate Species time histories during pyrolysis: Methyl Formate: CO, OH, C 2 H 4, CH 2 O, CH 3, CH 4, Me OH, MF Methyl Acetate/Propanoate: CO, CH 3, C 2 H 4 Methyl Butanoate: CO, CO 2, C 2 H 4, OH Reflected shock conditions: atm, K 20

21 Methyl Decanoate Ignition Delay Times: Aerosol Shock Tube Studies E A = 29.2 kcal/mol Westbrook model: 3500 species reactions E A = 42.5 kcal/mol Westbrook model correctly predicts ignition delay times, activation energies, and [O 2 ] dependence Can we apply AST method to larger esters? 21

22 Methyl Oleate Ignition Delay Times: Aerosol Shock Tube Studies Data reveal weak dependence on equivalence ratio Westbrook model: underestimates ignition delay times by about 50% captures reveal weak dependence on equivalence ratio Next step: Need for species time histories! 22

23 Methyl Formate Pyrolysis: Time Histories and Rate Data Species time histories: CO, OH, C 2 H 4, CH 2 O, CH 3, CH 4, MeOH, MF Rate constant determination for all three major decomposition channels (Princeton/NUI 2010) 1: MF CO+MeOH using CO (4.6 m) using MeOH (9.23 m) 2: MF CH 4 +CO 2 using CH 4 (3.4 m) using CO 2 (2.7 m) 3: MF 2CH2O using CH 2 O (306 nm) 23

24 Methyl Formate Decomposition: CH 3 OH+CO channel CO Mole Fraction [%] Measurement Dryer et al % MF/Ar 1.5 atm CO Time Histories 1488K 1376K 1285K 1202K Time [ s] 1607K k 1 [cm -3 mol -1 s -1 ] MF CH 3 OH+CO 1667K Current study Best fit Dryer et al (1.6 atm) Curran et al (1.6 atm) P = atm /T [K] Direct measurement of CH 3 OH+CO channel possible with CO laser Excellent agreement at with Dooley et al. (2010) particularly at lower T 24

25 Methyl Formate Decomposition: CH 4 +CO 2 channel CH 4 Time Histories 1667 K MF CH 4 +CO CH4 Mole Fraction [%] % MF/Ar 1.5 atm Current study Dryer et al Updated k K Time [ms] 1289K k 2 [cm 3 mol -1 s -1 ] Current study Best fit Dryer et al (1.45 atm) /T [K] P = atm Direct measurement of CH 4 +CO 2 channel possible with CH 4 or CO 2 laser Excellent agreement at with Dooley et al. (2010) 25

26 Methyl Formate Decomposition: CH 2 O+CH 2 O channel CH2O Mole Fraction [%] CH 2 O Time Histories Current study Dryer et al Updated k 1, k 2 3% MF/Ar 1.4 atm 1462K 1347K 1257K k 3 [cm 3 mol -1 s -1 ] 1667 K MF CH 2 O+CH 2 O P = atm Current study Best fit Dryer et al (1.4 atm) Time [ s] /T [K] Direct measurement of CH 2 O+CH 2 O channel possible with CH 2 O laser Excellent agreement at with Dooley et al. (2010) 26

27 Oxygen Balance in Methyl Formate Pyrolysis Oxygen Atom Balance CH 3 OCHO CH 4 (CO 2 ) T = 1420K, P = 1.5atm (0.2-3% MF/Ar) CO CH 3 OH CH 2 t = 300 s % Oxygen balance: 5.5% in MF 34.8% in MeOH 44.9% in CO 5.8% in CO2 7.2% in CH2O Time [ s] Total: 98.2% Laser data successfully tracks all major contributors to O atom balance Significant opportunity for mechanism validation 27

28 Future Work Butanol Kinetics Extend species time history database to all isomers Extend speciation database to oxidation Methyl Ester Kinetics Extend ignition delay time database to large bio diesel components Extend speciation studies to other methyl esters (small and large) Other Opportunities Collaborate with modelers Apply multi species methods to other fuels/surrogates 28