2011 us combustion meeting - kinetic modeling of methyl formate oxidation
DESCRIPTION
A kinetic model for methyl formate oxidation is generated using our open-source Reaction Mechanism Gen- erator (RMG) software, supplemented with high level quantum calculations and transition state theory (TST). New rate coefficients are calculated for the decomposition pathways of methyl formate, methoxy-formyl (CH3 OC · O), and formyloxy-methyl (C · H2 OCHO), and hydrogen abstractions from methyl formate by H and methyl radicals. We compare the predictions to experimental data including previously unpublished shock tube ignition delays over a wide range of T and P, as well as atmospheric-pressure laminar burning velocities and low-pressure flame species profiles from the literature. Using RMG we investigate the effect of changing the small molecule (C0−C1) “seed mechanism” and show that predictions of all the experiments are sensi- tive to these reactions. Until the small molecule chemistry is resolved it is impossible to have a conclusive mechanism for the fuel molecule oxidation.TRANSCRIPT
![Page 1: 2011 US Combustion Meeting - Kinetic Modeling of Methyl Formate Oxidation](https://reader034.vdocument.in/reader034/viewer/2022042715/55904fa51a28ab9f548b457f/html5/thumbnails/1.jpg)
Massachusetts Institute of Technology
Kinetic Modeling of Methyl Formate OxidationRichard H West, C Franklin Goldsmith, Michael R Harper, William H Green, Laurent Catoire, Nabiha Chaumeix
1a06
![Page 2: 2011 US Combustion Meeting - Kinetic Modeling of Methyl Formate Oxidation](https://reader034.vdocument.in/reader034/viewer/2022042715/55904fa51a28ab9f548b457f/html5/thumbnails/2.jpg)
We would like to model kinetics of biofuel combustion predictively.
•Automatically built models are more extensible than hand-curated models, and easier to update.
Hand curated models Automatically built models
Interpretation tool Design tool
Specific to validated fuels, T, P Extensible to new fuels, T, P
Chosen sub-models may be inconsistent with new data
Model can be rebuilt with new data, adding all new pathways
Accurate (where validated) Improving...
2
![Page 3: 2011 US Combustion Meeting - Kinetic Modeling of Methyl Formate Oxidation](https://reader034.vdocument.in/reader034/viewer/2022042715/55904fa51a28ab9f548b457f/html5/thumbnails/3.jpg)
Three steps towards predictive kinetic models of biofuels.
•Develop method to build kinetic models automatically.
•Collect or generate data and rules suitable for biofuels.
•Perform experiments to validate and check understanding.
3
![Page 4: 2011 US Combustion Meeting - Kinetic Modeling of Methyl Formate Oxidation](https://reader034.vdocument.in/reader034/viewer/2022042715/55904fa51a28ab9f548b457f/html5/thumbnails/4.jpg)
Three steps towards predictive kinetic models of biofuels.
•Develop method to build kinetic models automatically.
•Collect or generate data and rules suitable for biofuels.
•Perform experiments to validate and check understanding.
⇌���facebook.com/rmg.mitr m g . s o u rc e f o r g e . n e t
Reaction Mechanism Generator
•free and open source software
•version 3.3 released last month
3
![Page 5: 2011 US Combustion Meeting - Kinetic Modeling of Methyl Formate Oxidation](https://reader034.vdocument.in/reader034/viewer/2022042715/55904fa51a28ab9f548b457f/html5/thumbnails/5.jpg)
Methyl Formate is a useful moleculeto study for biofuel modeling.
•Biodiesel consists of methyl esters.
•Methyl formate is the smallest methyl ester.
•Rules that can predict methyl formatecan likely predict other esters.
��
Methyl Formate
Biodiesel
Diesel
4
![Page 6: 2011 US Combustion Meeting - Kinetic Modeling of Methyl Formate Oxidation](https://reader034.vdocument.in/reader034/viewer/2022042715/55904fa51a28ab9f548b457f/html5/thumbnails/6.jpg)
Overview of this study
•Use RMG to build a kinetic model
5
![Page 7: 2011 US Combustion Meeting - Kinetic Modeling of Methyl Formate Oxidation](https://reader034.vdocument.in/reader034/viewer/2022042715/55904fa51a28ab9f548b457f/html5/thumbnails/7.jpg)
Overview of this study
•Use RMG to build a kinetic model
•Improve some estimates
•thermochemistry data from quantum calculations
•decomposition kinetics of fuel molecule and radicals
6
![Page 8: 2011 US Combustion Meeting - Kinetic Modeling of Methyl Formate Oxidation](https://reader034.vdocument.in/reader034/viewer/2022042715/55904fa51a28ab9f548b457f/html5/thumbnails/8.jpg)
Overview of this study
•Use RMG to build a kinetic model
•Improve some estimates
•Use improved RMG to build new kinetic models
•use 3 different ‘seed’ mechanisms from the literature
7
![Page 9: 2011 US Combustion Meeting - Kinetic Modeling of Methyl Formate Oxidation](https://reader034.vdocument.in/reader034/viewer/2022042715/55904fa51a28ab9f548b457f/html5/thumbnails/9.jpg)
Overview of this study
•Use RMG to build a kinetic model
•Improve some estimates
•Use improved RMG to build new kinetic models
•Compare with experiments
•new dilute shock tube ignition delay times
•flame speed and low-P flame profile from Dooley et al. (Dryer group at Princeton).
8
![Page 10: 2011 US Combustion Meeting - Kinetic Modeling of Methyl Formate Oxidation](https://reader034.vdocument.in/reader034/viewer/2022042715/55904fa51a28ab9f548b457f/html5/thumbnails/10.jpg)
Overview of this study
•Use RMG to build a kinetic model
•Improve some estimates
•Use improved RMG to build new kinetic models
•Compare with experiments
•Conclude
•choice of small-molecule seed mechanism (C0-C1 sub-model) affects all other conclusions
9
![Page 11: 2011 US Combustion Meeting - Kinetic Modeling of Methyl Formate Oxidation](https://reader034.vdocument.in/reader034/viewer/2022042715/55904fa51a28ab9f548b457f/html5/thumbnails/11.jpg)
RMG generates secondary chemistry of byproducts, exploring paths with high flux.
•Reaction family templates generate all possible reactions.
•‘Edge’ species with highest flux are moved to ‘core’.
•Expansion continues until tolerance reached.
AB
CD
E
FG
H
AB
CD
E
F
10
![Page 12: 2011 US Combustion Meeting - Kinetic Modeling of Methyl Formate Oxidation](https://reader034.vdocument.in/reader034/viewer/2022042715/55904fa51a28ab9f548b457f/html5/thumbnails/12.jpg)
Thermochemistry for ~200 small molecules from high-level QCI//DFT calculations
!
!"
! !" #
!
!"
!
! !"
! !"#$
!" !
!"
!"#$ !"# !""""
!" !" !"# "#! !"#
! !" ! !"
!"#$ #$!
"
!"#$
! !"
!"" # !"" # !"#
!"#$"
!"#$ #$!
! !"
!"
!
!
!
!"
!
!"#
!
!
! !
!"""# !""#$
!""#$
!" !"#$
!"
#
! !"#$
%!
!"#$
%
!"#$
%
!"#$
%!
!"#$
% %!
!"
!
!"
!" #
!"#$
%
%!
% %! !"
#$ !"#$
%
%! ! !"#$ % !
"#$ %
%
!"#$ %!
!"#$
!"#$# !" # !" ! "#$%
!
!"# !"#
$#
!"# !"#
$
!
!
!"#
$
!"
!
!"#
$#
!"#$
%#
!
!"
! !#
!
! !"
!
!"
! !# ! !"
!
!
!"
!
!
!"
!
!
!"
!"
#$ !" !"#$ % ! !
"# !
!"
!
!"
!
!
!
!"
!"#
!"#
$
!"#
$!
!"#
$#
!"#
$ $!
!"" # !!
!" !"
!"!
"!
!
!"
!
!!" !
!
!" ! ! !"" !"" !"#
!"
!"# !"# $ !"# $ !" !#
!" #
$#
!" "
"
"#
!" !
! !#
!" "
#$!
"
"!
!" "!
!" #
#
$#
!"#$
%
% %! !"
"
!"#$
!"
!"
!"#! !" ! !"
!
!
!!" !
!
!
!!" !"
!"#$
!
!"
!
!"
!"" !" !" !" !" ! ! !" !# !
!"
!"
! !#
!"
!
!
! !"
!
"# !
!
!"! !"# !
! !"
! ! "#$ !
"#$!
!"#
$
!"
!
!"#$
%
!"
!
!
!
!
!"
!"! !" !!" !" !" ! !" !"
! !"
!"#
!"
!
!"#$
%
%!
!"#
$
$#
!"#$
% %!
!"
!"#$
%
!
!
!"
• conformer search:CBS-QB3
• geometry and frequencies:B3LYP/6-311++G(d,p)
• scan for hindered rotors:B3LYP/6-31+G(d,p)
• electronic energy:RQCISD(T)/CBS extrapolated from RQCISD(T)/cc-pVTZ and RQCISD(T)/cc-pVQZ
11
![Page 13: 2011 US Combustion Meeting - Kinetic Modeling of Methyl Formate Oxidation](https://reader034.vdocument.in/reader034/viewer/2022042715/55904fa51a28ab9f548b457f/html5/thumbnails/13.jpg)
Same QM methods used to calculate PES of methylformate & radical decomposition
•High-Pressure rate coefficients calculated with TST:
��� ��������������������������� ��������������������
"4�@3*�8.+624).+2/786=�,46�+').�4,�� �75+)/+7�8.'8�'55+'6+*�/3�'�56+:/497� ���-+3+6'8+*�0/3+8/)24*+1�4,��� 4</*'8/43� ;+�5+6,462+*�496�4;3� ���!��"����!�����$���������*�5��)'1)9�1'8/437�97/3-�8.+�2+8.4*414-=�4,��41*72/8. �����%�&��11���" '3*���!�����)'1)91'8/437�;+6+*43+�97/3-��'977/'3�%&� �11�����'3*����!��"� )'1)91'8/437�;+6+�*43+�97/3-����� �% &�
!/2/1'6�)'1)91'8/437�;+6+�5+6,462+*�43�8.+�86'37/8/43�78'8+7�4,�*+)42547/8/43�6+')8/437�4,�2+8.=1,462'8+� 2+8.4<=�,462=1�������
·��� '3*�,462=14<=�2+8.=1��� ·�������� '3*�4,�.=*64-+3�'(�786')8/437�,642�2+8.=1�,462'8+�(=�� '3*�2+8.=1�6'*/)'17� �46�8.+ ��� �� '3* ������ *+)42�547/8/43�6+')8/437� 8.+�7'2+�./-.�1+:+1� ���!��"� )'1)91'8/437�'7�,46�8.+�8.+624).+2/786=�;+6+5+6,462+*� �46���'(786')8/43�6+')8/437� 8.+���!�����+1+)8643/)�+3+6-=�;'7�97+*� �/3*+6+*�64846)'1)91'8/437�;+6+�86+'8+*�/3�8.+�7'2+�2'33+6�'7�%�&� ".+�./-.�56+7796+�1/2/8�,46�+').�6+')8/43�6'8+)4+,@)/+38��"'(1+ ��;'7�)42598+*�97/3-�86'37/8/43�78'8+�8.+46=�%�&�
��� �� � ������������������ ������������������������ k =A(T/1K)n exp(!Ea/RT )� ���������������� ������ ���
+')8/43 A �)23� 241� 7� n Ea �0)'1�241�#3/241+)91'6�*+)42547/8/43������!��"� ��������!!!"!! ��+����� 6.49"104 ���� � � �������!!!"!! ���+�� 6.88"106 �� � ����� ·�!!!"!! ���+� ·�� 2.34"109 �� ����+���� · !!!"!! ����� ·� 1.97"109 ��� ����� ·�+����!!!"!! � ·������ 5.03"103 �� � ����
��'(786')8/43�����!������� ·+�������!!!"!! ��+�����
·� 2.28"105 ��� ����� ·+�������!!!"!! ��+�
·������ 1.16"105 ���� ����� ·��+�������!!!"!! �� +����� ·� 6.34 ���� ���� ·��+�������!!!"!! �� +� ·������ 0.257 ���� ���
'*/)'1� '*/)'1�6+)42(/3'8/437��+78/2'8+*����� ·+�� ·�!!!"!! ������� 1.8"1013 � ·��+�
·���!!!"!! ������� 1.8"1013
��� ��� ����������
".6++�*/,,+6+38�7++*�2+).'3/727�'6+�97+*�/3�8./7�789*=� ;./).�'6+�6+,+66+*�84�'7�� ����1'6(46-� '3*�441+=� ".+�392(+6�4,�75+)/+7�'3*�6+')8/437�,46�+').�7++*�2+).'3/72�/7�7.4;3�/3�"'(1+ ��
��� ".+�� � 7++*�2+).'3/72�/7�)437869)8+*�,642�8.+�� ���+).����2+8.'3+�4</*'8/43�2+).�'3/72�%�&� �11�75+)/+7�'3*�6+')8/437�/3:41:/3-�3/864-+3�)438'/3/3-�)425493*7�'3*�8.+�1925+*���� 75+)/+7� ;./).� �� 24*+17�'7�8;4�*/78/3)8�75+)/+7� ;+6+�6+24:+*�
������ ".+��1'6(46-�7++*�2+).'3/72�)438'/37�'11 �!� 6+')8/437�,642�8.+��1'6(46-�-6495?7�'78+6��+).'3/72� 2478�6+)+381=�6+5468+*�(=��45+> ����� %�&� �11�75+)/+7�'3*�6+')8/437�/3:41:/3-3/864-+3�)438'/3/3-�)425493*7�;+6+�6+24:+*�
Methyl Formate PESFormyloxymethyl and Methyloxyformyl PES
12
kcal/mol
90
60
30
0
-30
O OCO + CH3OH
CH4 + CO2
kcal/mol
30
20
10
0
-10
40
-20
O OO O
OO
OO
OO +
O
OH +
CO + CH3O·
CH3 + CO2
80.8
-27.2
67.5
9.60.0
![Page 14: 2011 US Combustion Meeting - Kinetic Modeling of Methyl Formate Oxidation](https://reader034.vdocument.in/reader034/viewer/2022042715/55904fa51a28ab9f548b457f/html5/thumbnails/14.jpg)
Pressure-dependence calculated in RMG using reservoir state method
•RMG explores additional pathways and solves full Master Equation for chemically activated and fall-off reaction rates
�����
����
����
�����
�� ��
�����
�� ��������� ��
������
��� ������� ��
����
����� �������� ��� ���
���� ���������� ���
����
����� ������ ���
����
��� ������ ���
�����
�����
���
�������
�������
���
�����
���
��
���
����
���
���� ���
�������
����
���
����
�������������������
����
��������� ���������
�����
������ ������
����
������� �������
����
��������
����
��������
����
���������� �����
� �������� �������
���
���������� ������
�����
���������� ����
���
�����������
���
�� �������� ��������� ����
�� �������� ������
����
������������������
��
��������������������
���
��������������������
����
�� �������������
���
����������� �����
����������
���
������������� ����
������������� ���
�������������
��
�������� ��������
�����
����������
Methyl Formate PESFormyloxymethyl and Methyloxyformyl PES
13
O O
O OO O
![Page 15: 2011 US Combustion Meeting - Kinetic Modeling of Methyl Formate Oxidation](https://reader034.vdocument.in/reader034/viewer/2022042715/55904fa51a28ab9f548b457f/html5/thumbnails/15.jpg)
“Seed” mechanism placed in RMG’s core at start of mechanism generation
•Used for trusted rates of small-molecule chemistry.
•Rate expressions used without modification.
•Additional reactions generated by RMG.
AB
CD
E
F
AB
C
DE
HG
F
Seed Mechanism 1
Seed Mechanism 2
14
![Page 16: 2011 US Combustion Meeting - Kinetic Modeling of Methyl Formate Oxidation](https://reader034.vdocument.in/reader034/viewer/2022042715/55904fa51a28ab9f548b457f/html5/thumbnails/16.jpg)
“GRI”
•GRI-Mech 3.0
•(Methane
oxidation)
•Smith (1999)
Three Seed Mechanisms were taken from published literature.
������ &30��9960C�=00/�70.3,84=7�.98>,48=�,66 ��!� <0,.>498=�1<97�>30��<C0<�2<9?:F=�70>3C619<7,>0�9B4/,>498�70.3,84=7�-C��9960C ������*+� �>�,6=9�48.6?/0=�79=>�91�>30 �!�� =:0.40=��>304<<0,.>498=�,<0�0=>47,>0/�-C�$����
��� ��� ��� ��� � �����
&30<79.3074.,6�/,>,�A0<0�>,508�1<97��48�9</0<�91�:<010<08.0�� >30�$#��%��&������)" /,>,-,=0/0=.<4-0/�48�=0.>498 �� �$���0.3����� =074�07:4<4.,6�"���.,6.?6,>498=��19<�.C.64.�=:0.40=�986C�A4>3�3C/<9208�-98/�48.<0708>�.9<<0.>498=������ 19<�<,/4.,6=� 0=>47,>0=�-,=0/�98��08=98F=�2<9?:,//4>4@4>C�70>39/�A4>3���� 19<�<,/4.,6=� &<,8=:9<>�:,<,70>0<=�A0<0�>,508�1<97��$���0.3����A30<0�,@,46,-60� ,8/�0=>47,>0/�-C�$�� 9>30<A4=0� $0,.>498�<,>0�.901G.408>=�A0<0�>,508�1<97��489</0<�91�:<010<08.0�� >30�=:0.4G0/�=00/�70.3,84=7��/4110<08>�19<�0,.3�79/06� =00�=0.>498 ����A4>3>30�,==?7:>498�>3,>�:<0==?<0�/0:08/08.0�3,=�,6<0,/C�-008�,..9?8>0/�19<� >30 ��������� <,>0=�=39A848�&,-60 � >30��6,<-9<2��,=>0<�70.3,84=7�*�+� 0=>47,>0=�-,=0/�98�$��F=�<0,.>498�1,746C�5480>4.=0=>47,>498�<?60=� '4>3�>30�0B.0:>498�91�>39=0�1<97�>30�=00/�70.3,84=7=� ,66�<,>0�.901G.408>=�,<0,==?70/�>9�-0�<0:9<>0/�48�>30�3423�:<0==?<0�6474>�,8/��,=>0<��;?,>498�.,6.?6,>498=�,<0�:0<19<70/>9�.9<<0.>�19<�1,66�911�,8/�.3074.,6�,.>4@,>498� &30��,=>0<��;?,>498�.,6.?6,>498=�,<0�:0<19<70/98�,��,?==��30-C=30@�2<4/�91���>07:0<,>?<0=����E� �� �� ,8/� �:<0==?<0=�����E�� -,<��,8/�30-C=30@�:96C8974,6=�91�9</0<���48 T ,8/���48 P ,<0�G>>0/�>9�>30�<0=?6>=�19<�?=0�48������� ���66>30�=:0.40=�<0:9<>0/�48�>30�69A�:<0==?<0�H,70�0B:0<4708>=*�+�A0<0�,//0/�>9�>30�79/06�.9<0�A4>3D0<9�.98.08><,>498��>34=�08=?<0=�>30�=:0.40=�,::0,<�48�>30�G8,6�79/06� 0@08�41�>304<�<,>0�91�.<0,>498�4=:<0/4.>0/�>9�-0�69A�� %:0.40=�A4>3�79<0�>3,8���9BC208�,>97=�A0<0�19<-4//08� &34<>C�>A9�<0,.>498=C=>07=�A0<0�=47?6,>0/� 19<70/�-C�>30�:<9/?.>�91���>07:0<,>?<0=������ ���� ��� ������ " �:<0==?<0=������� � �� ��-,<� " ���!��< 796,<�<,>49=�������,8/���� -9>3 ! = 1�� &30��<0,.>498�=C=>07=�A0<0�0,.3�=47?6,>0/�19<��� �=0.98/=�?8/0<�4=9>30<7,6�4=9-,<4.�-,>.3�.98/4>498=A4>3�>30�0/20�H?B0=��>9�/0>0<7480�A3,>�80B>�>9�48.6?/0�48�>30�.9<0��0@,6?,>0/�,>�0@0<C�>470�=>0:91�>30�=>411�!�� =96@0<� &30�>0<748,>498�.<4>0<498�A,=���� �,6>39?23�89>�,66�<0,.>498�=C=>07=�<,8>9�.97:60>498� &30�G8,6�79/06=�A0<0�2080<,>0/�A4>3�,�<060,=0�.,8/4/,>0�91�$�� @0<=498������:<0.4=0�@0<=498�3,=3� �/, 1�� �9/06�2080<,>498�>995�E��39?<=�:0<�70.3,84=7�98�,�=48260��� ��D��8>06�(098�:<9.0==9<� ,6>39?23�79=>�91�>30 #�����"���.,6.?6,>498=�19<�.C.64.�=:0.40=3,/�,6<0,/C�-008�:0<19<70/��0B4=>482�<0=?6>=�,<0�<0�?=0/�� &30�=4D0=�91�>30�79/06=�,<0�=?77,<4D0/48�&,-60 �
�� �� ���������������"��� ����#����������$��������%���������� ����������������� ��� ��������� ����#����������$����� �������������������������������������� �������������������������!�� �� �!�������������
%00/�70.3,84=7 �9/06�=00/ �9/06�.9<0 �9/06�0/20=:0.40= <0,.>498= =:0.40= <0,.>498= =:0.40= <0,.>498=
�$� �� �� � ��� ��� ����6,<-9<2 � �� � ���� � �� �����9960C � �� ��� ��� ����
“Glarborg”
•Glarborg group’s “Master Mechanism”
•C0-C1 only
•Lopez (2009)
“Dooley”•Dryer group’s mechanism
•Used in methyl formate study
•C0-C1 only•Dooley (2010)
15
![Page 17: 2011 US Combustion Meeting - Kinetic Modeling of Methyl Formate Oxidation](https://reader034.vdocument.in/reader034/viewer/2022042715/55904fa51a28ab9f548b457f/html5/thumbnails/17.jpg)
The models are built using a hierarchy of data sources.
16
QCISD(T)//B3LYP database
Semi-empirical PM3
Seed MechanismGRI / Glarborg / Dooley
Methyl-formatecalculations
Glarborg Master Mechanism
Rule-based estimates
Benson group estimates
GRI-Mech
Group-based estimates
Kinetics Thermochemistry Transport
150
10
100
4000
![Page 18: 2011 US Combustion Meeting - Kinetic Modeling of Methyl Formate Oxidation](https://reader034.vdocument.in/reader034/viewer/2022042715/55904fa51a28ab9f548b457f/html5/thumbnails/18.jpg)
The models are built using a hierarchy of data sources.
16
QCISD(T)//B3LYP database
Semi-empirical PM3
Seed MechanismGRI / Glarborg / Dooley
Methyl-formatecalculations
Glarborg Master Mechanism
Rule-based estimates
Benson group estimates
GRI-Mech
Group-based estimates
Kinetics Thermochemistry Transport
k(T,P) calculated by Master Equation solution in RMG (+3000)ME
ME
ME
ME
150
10
100
4000
![Page 19: 2011 US Combustion Meeting - Kinetic Modeling of Methyl Formate Oxidation](https://reader034.vdocument.in/reader034/viewer/2022042715/55904fa51a28ab9f548b457f/html5/thumbnails/19.jpg)
The models are built using a hierarchy of data sources.
16
QCISD(T)//B3LYP database
Semi-empirical PM3
Seed MechanismGRI / Glarborg / Dooley
Methyl-formatecalculations
Glarborg Master Mechanism
Rule-based estimates
Benson group estimates
GRI-Mech
Group-based estimates
Kinetics Thermochemistry Transport
Radicals calculated by Hydrogen Bond Increment methodHBI
HBI
HBI
k(T,P) calculated by Master Equation solution in RMG (+3000)ME
ME
ME
ME
150
10
100
4000
![Page 20: 2011 US Combustion Meeting - Kinetic Modeling of Methyl Formate Oxidation](https://reader034.vdocument.in/reader034/viewer/2022042715/55904fa51a28ab9f548b457f/html5/thumbnails/20.jpg)
The models are built using a hierarchy of data sources.
16
QCISD(T)//B3LYP database
Semi-empirical PM3
Seed MechanismGRI / Glarborg / Dooley
Methyl-formatecalculations
Glarborg Master Mechanism
Rule-based estimates
Benson group estimates
GRI-Mech
Group-based estimates
Kinetics Thermochemistry Transport
Radicals calculated by Hydrogen Bond Increment methodHBI
HBI
HBI
k(T,P) calculated by Master Equation solution in RMG (+3000)ME
ME
ME
ME Constant across all models in this work
150
10
100
4000
![Page 21: 2011 US Combustion Meeting - Kinetic Modeling of Methyl Formate Oxidation](https://reader034.vdocument.in/reader034/viewer/2022042715/55904fa51a28ab9f548b457f/html5/thumbnails/21.jpg)
Shock tube experiments performed by CNRS in Orleans.
•99% Argon
•Equivalence ratios ɸ = 0.5, 1, 2
•Ignition defined as half of peak OH* emission
•P5 = 151 – 1446 kPa
•T5 = 1268 – 1812 K
High Pressure Low Pressure
P transducers
UV photomultiplier
! "!! #!! $!! %!! &!!!
'()*+, -.
!
!
Pressure
- OH* emission
delay
17
![Page 22: 2011 US Combustion Meeting - Kinetic Modeling of Methyl Formate Oxidation](https://reader034.vdocument.in/reader034/viewer/2022042715/55904fa51a28ab9f548b457f/html5/thumbnails/22.jpg)
Unscaled ignition delay times (IDT) are hard to interpret plotted vs 1/T
Phi = 10 p scaling
0.908841851 correl0.7425 optimum
0.688625 global optimum2.690715605
!""#
!"""#
!""""#
"$%"# "$&"# "$'"# "$("#
!"#$%&'()
*+++),)%)$-)
)*+,-./,01#
Phi = 20 p scaling
0.89880299 correl0.64875 optimum
0.688625 global optimum
!""#
!"""#
!""""#
"$%"# "$&"# "$'"# "$("#
!"#$%&'()
*+++),)%)$-)
)*+,-./,01#
Phi = 0.50 p scaling
0.883070764 correl0.689325 optimum0.688625 global optimum
!""#
!"""#
!""""#
"$%"# "$&"# "$'"# "$("#
!"#$%&'()
*+++),)%)$-)
)*+,-./,01#
•Each experiment performed at a different P5 pressure.
ɸ=0.5 ɸ=1 ɸ=2
18
![Page 23: 2011 US Combustion Meeting - Kinetic Modeling of Methyl Formate Oxidation](https://reader034.vdocument.in/reader034/viewer/2022042715/55904fa51a28ab9f548b457f/html5/thumbnails/23.jpg)
Phi = 20 p scaling
0.89880299 correl0.64875 optimum
0.688625 global optimum
!""#
!"""#
!""""#
"$%"# "$&"# "$'"# "$("#
!"#$%&'()
*+++),)%)$-)
)**+,-#./0#123,456,78#.+94:*4;#
Phi = 10 p scaling
0.908841851 correl0.7425 optimum
0.688625 global optimum2.690715605
!""#
!"""#
!""""#
"$%"# "$&"# "$'"# "$("#
!"#$%&'()
*+++),)%)$-)
)**+,-#./0#123,456,78#.+94:*4;#
Phi = 0.50 p scaling
0.883070764 correl0.689325 optimum0.688625 global optimum
!""#
!"""#
!""""#
"$%"# "$&"# "$'"# "$("#
!"#$%&'()
*+++),)%)$-)
)**+,-#./0#123,456,78#.+94:*4;#
ɸ=0.5 ɸ=1 ɸ=2
Unscaled ignition delay times (IDT) are hard to interpret plotted vs 1/T
•Each experiment performed at a different P5 pressure.
•Simulate each experiment with the corresponding T5, P5.
19
![Page 24: 2011 US Combustion Meeting - Kinetic Modeling of Methyl Formate Oxidation](https://reader034.vdocument.in/reader034/viewer/2022042715/55904fa51a28ab9f548b457f/html5/thumbnails/24.jpg)
Scale ignition delay times by P5α
to improve correlation with 1/T
Phi = 10 p scaling
0.908841851 correl0.7425 optimum
0.688625 global optimum2.690715605
!""#
!"""#
!""""#
"$%"# "$&"# "$'"# "$("#
!"#$%&'()
*+++),)%)$-)
)*+,-./,01#
Phi = 20 p scaling
0.89880299 correl0.64875 optimum
0.688625 global optimum
!""#
!"""#
!""""#
"$%"# "$&"# "$'"# "$("#
!"#$%&'()
*+++),)%)$-)
)*+,-./,01#
Phi = 0.50 p scaling
0.883070764 correl0.689325 optimum0.688625 global optimum
!""#
!"""#
!""""#
"$%"# "$&"# "$'"# "$("#
!"#$%&'()
*+++),)%)$-)
)*+,-./,01#
ɸ=0.5 ɸ=1 ɸ=2
•Each experiment performed at a different P5 pressure.
20
![Page 25: 2011 US Combustion Meeting - Kinetic Modeling of Methyl Formate Oxidation](https://reader034.vdocument.in/reader034/viewer/2022042715/55904fa51a28ab9f548b457f/html5/thumbnails/25.jpg)
Scale ignition delay times by P50.69
to maximize correlation with 1/T
Phi = 10.688625 p scaling
0.982161783 correl0.7425 optimum
0.688625 global optimum2.946534131
!""#
!"""#
!""""#
"$%"# "$&"# "$'"# "$("#
!"#$%&'()!*
+%,-.(/0123
4///353%3$+3
)*+,-./,01#
Phi = 20.688625 p scaling
0.982344879 correl0.64875 optimum
0.688625 global optimum
!""#
!"""#
!""""#
"$%"# "$&"# "$'"# "$("#
!"#$%&'()!*
+%,-.(/0123
4///353%3$+3
)*+,-./,01#
Phi = 0.50.688625 p scaling
0.982027468 correl0.689325 optimum0.688625 global optimum
!""#
!"""#
!""""#
"$%"# "$&"# "$'"# "$("#
!"#$%&'()!*
+%,-.(/0123
4///353%3$+3
)*+,-./,01#
ɸ=0.5 ɸ=1 ɸ=2
21
![Page 26: 2011 US Combustion Meeting - Kinetic Modeling of Methyl Formate Oxidation](https://reader034.vdocument.in/reader034/viewer/2022042715/55904fa51a28ab9f548b457f/html5/thumbnails/26.jpg)
Phi = 0.50.688625 p scaling
0.982027468 correl0.689325 optimum0.688625 global optimum
!""#
!"""#
!""""#
"$%"# "$&"# "$'"# "$("#
!"#$%&'()!*
+%,-.(/0123
4///353%3$+3
)**+,-#./0#.+123*24##567,289,:;#
Phi = 10.688625 p scaling
0.982161783 correl0.7425 optimum
0.688625 global optimum2.946534131
!""#
!"""#
!""""#
"$%"# "$&"# "$'"# "$("#
!"#$%&'()!*
+%,-.(/0123
4///353%3$+3
)**+,-#./0#.+123*24##567,289,:;#
Phi = 20.688625 p scaling
0.982344879 correl0.64875 optimum
0.688625 global optimum
!""#
!"""#
!""""#
"$%"# "$&"# "$'"# "$("#
!"#$%&'()!*
+%,-.(/0123
4///353%3$+3
)**+,-#./0#.+123*24##567,289,:;#
ɸ=0.5 ɸ=1 ɸ=2
•All too fast
•Glarborg seed fastest
•Dooley seed slowest
•Good agreement
•Glarborg seed slowest
•Dooley seed fastest
•All too slow
•Glarborg seed slowest
•GRI-Mech seed fastest
Shock tube ignition delays predicted OK, but vary with ɸ and seed mechanism.
22
![Page 27: 2011 US Combustion Meeting - Kinetic Modeling of Methyl Formate Oxidation](https://reader034.vdocument.in/reader034/viewer/2022042715/55904fa51a28ab9f548b457f/html5/thumbnails/27.jpg)
!"#!"$
%&'
(!'
(&'
)!'
)&'
*!'
*&'
!"+' !"#' %' %"(' %"*' %"+'
!"#$%"&'()&%$%*'+,-./$01
'2/#345'
67)$+"-,%/,'8"9.'
,-./01/2'3445',67'3445'811-49'3445':;<4/=>4?@'A811-49B'(!%!C'811-49'A(!%!C'>154-'
Laminar flame speeds all over-predicted, and also vary with seed mechanism.
•Experimental data and model from Dooley et al. (2010) Int. J. Chem. Kinet.
23
![Page 28: 2011 US Combustion Meeting - Kinetic Modeling of Methyl Formate Oxidation](https://reader034.vdocument.in/reader034/viewer/2022042715/55904fa51a28ab9f548b457f/html5/thumbnails/28.jpg)
Seed mechanism is mostly responsible for over-predicted flame speed.
•For model with GRI seed mechanism, the10 most sensitive reactions are all from the GRI seed.
24
Reac%on Sensi%vityH.0+0O20=0O0+0OH 0.22OH0+0CO0=0H.0+0CO2 0.15H.0+0HO20=0O20+0H2 90.09H.0+0HO20=0OH0+0OH 0.09H.0+0O20+0H2O0=0HO20+0H2O 90.08HC.O0+0H2O0=0H.0+0CO0+0H2O 0.08HC.O0+0m0=0H.0+0CO0+0m 0.07HC.O0+0O20=0HO20+0CO 90.07HO20+0CH30=0OH0+0CH3O. 0.06H.0+0HC.O0=0H20+0CO 90.05
Sensitivity of flame speed with respect to reaction A factor
!"#!"$
%&'
(!'
(&'
)!'
)&'
*!'
*&'
!"+' !"#' %' %"(' %"*' %"+'
!"#$%"&'()&%$%*'+,-./$01
'2/#345'
67)$+"-,%/,'8"9.'
,-.'/001'
2340567089':;<<=0>?'(!%!@'
![Page 29: 2011 US Combustion Meeting - Kinetic Modeling of Methyl Formate Oxidation](https://reader034.vdocument.in/reader034/viewer/2022042715/55904fa51a28ab9f548b457f/html5/thumbnails/29.jpg)
Minor species profiles poorly predicted in low-pressure flame. Seed has some effect.
•Low pressure (22 torr) stoichiometric (ɸ=1) flame data and temperature profile from Dooley et al. (2010) Comb. Flame.
•Overpredicted C2H4 in model with Glarborg seedis due to missing decomposition pathway.
!"#$
%&'()%($*
('()+((*
&'()%($*
,'()%(-*
,'&)%(-*
"'()%(-*
(* "* $* .* /* ,(*
!"#$%&'()*"
+%
,-./(+)$%&'"0%12'+$'%3%00%
01234536*078*9551:;*)<=:3>?:@A*
!"#
$%&'($')*'&'(+''*%&'($')*,&'($'#*,&%($'#*-&'($'#*-&%($'#*#&'($'#*#&%($'#*)&'($'#*
'* -* )* .* /* ,'*
!"#$%&'()*"
+%
,-./(+)$%&'"0%12'+$'%3%00%
01234536*078*9551:;*(<=:3>?:@A*
CH3 mole fraction C2H4 mole fraction
25
![Page 30: 2011 US Combustion Meeting - Kinetic Modeling of Methyl Formate Oxidation](https://reader034.vdocument.in/reader034/viewer/2022042715/55904fa51a28ab9f548b457f/html5/thumbnails/30.jpg)
Pressure-dependent ethyl formate reactions important to C2H4 profile
•Replacing high-pressure-limit ethyl formate reactions with chemically activated pressure-dependent network increases prediction of C2H4 and C2H2 concentrations.
C2H2 mole fractionC2H4 mole fraction
26
![Page 31: 2011 US Combustion Meeting - Kinetic Modeling of Methyl Formate Oxidation](https://reader034.vdocument.in/reader034/viewer/2022042715/55904fa51a28ab9f548b457f/html5/thumbnails/31.jpg)
Acknowledgments
•C Franklin Goldsmith
•Michael R Harper
•William H Green
•Laurent Catoire
•Nabiha Chaumeix
•Stephen Dooley
27
![Page 32: 2011 US Combustion Meeting - Kinetic Modeling of Methyl Formate Oxidation](https://reader034.vdocument.in/reader034/viewer/2022042715/55904fa51a28ab9f548b457f/html5/thumbnails/32.jpg)
Contributions
•Improved tools to automatically build kinetic models.
•Calculated new data for methyl esters.
•New shock tube ignition delay data for methyl formate.
•Information in hard-to-read fonts is better remembered than easier to read information (Diemand-Yauman et al., Cognition, 2011)
28
![Page 33: 2011 US Combustion Meeting - Kinetic Modeling of Methyl Formate Oxidation](https://reader034.vdocument.in/reader034/viewer/2022042715/55904fa51a28ab9f548b457f/html5/thumbnails/33.jpg)
Contributions
•Improved tools to automatically build kinetic models.
•Calculated new data for methyl esters.
•New shock tube ignition delay data for methyl formate.
•Choice of C0-C1 kinetic model affects all other conclusions. Hard to “validate” other parts of model when this remains unsure.
29