Premix CH 4/Air flame with scheme GRI12
Contents
Objectives
The Stochastic_GRI12 test case describes a reduction of the GRI1.2 scheme for a 1D premixed flame. Starting with 32 species and 177 reactions, we reduce to 14 species and 26 reactions.
Key parameters
The target species considered for this test case are O2, CO and CO2. The characteristics of the premixed flame are displayed below :
listFlames.push_back(new PremixedFlames(/*T_fuel*/ 300,
/*T_ox*/ 300,
/*Pressure*/ 1E+05,
/*Ratio*/ 0.75,
/*Y_fuel*/ "CH4:1.0",
/*X_fuel*/ "",
/*Y_ox*/ "O2:0.21, N2:0.79",
/*X_ox*/ "",
/*reference_flame_path*/ "flames/flame__Phi_0_75__P_100000__T_300.cantera",
/*flame_path*/ "flames/flame"));
Results
DRGEP step
While running the DRGEP species step, your terminal should display the following information :
scheme : GRI12
number of species : 32
number of reaction : 177
reference flame : flames/flame__Phi_0_75__P_100000__T_300.cantera
Pressure : 1 bar
equivalence ration : 0.75
Composition : <O2:0.202024>
<CH4:0.0379821>
<N2:0.759994>
followed by the the species associated with their rank :
4.17578965277519e-189 AR 1.81320667696424e-06 HCCOH 3.33272520534669e-05 C2H 0.000136594690602079 CH2CO 0.000228724937160335 C 0.00039194258246205 C2H2 0.000909271897890726 HCCO 0.00324304055657054 CH 0.00391278892730799 C2H3 0.00395833642581853 CH3OH 0.00517918655958378 CH2OH 0.0158390004661296 C2H4 0.0158390004661296 C2H5 0.0160151369051444 CH2(S) 0.0209164664682532 C2H6 0.0317122146887158 CH2 0.0395596174129288 H2O2 0.0796917798682698 H2 0.0949180222420726 CH3 0.108929211050051 HCO 0.123697044495215 CH2O 0.154156160901222 H2O 0.165970760424094 CH3O 0.222439341655666 O 0.292243743837601 HO2 0.501413374426018 OH 0.501449580770414 H 1 O2 1 CH4 1 CO 1 CO2 1 N2
From a detailed 32 species scheme, we obtain reduced schemes (31 to 14 species). and we choose the one with 14 species for the next step.
Fig1 : Comparison between the reference trajectories of the target species, the temperature and the flame speed (in black), and the trajectories computed with the reduced mechanism with 14 transported species and 42 reactions, after a DRGEP species reduction (in red).
From the reduced scheme with 14 species, the DRGEP reaction step displays the associated 42 reactions (forward, reverse and global) with their rank :
Forward reaction 1 2.00309523982856e-05 Forward reaction 11 3.18147185453261e-05 Forward reaction 2 5.32566156092659e-05 Forward reaction 22 0.000208990760221458 Forward reaction 12 0.000233718026775223 Forward reaction 23 0.00168454263548541 Forward reaction 6 0.00205742859762978 Forward reaction 18 0.00214115102965165 Forward reaction 31 0.00256497802215946 Forward reaction 36 0.00423851339874799 Forward reaction 38 0.00628487407561224 Forward reaction 10 0.00709110759390656 Forward reaction 34 0.00727675620154213 Forward reaction 7 0.00816841804675756 Forward reaction 32 0.0105674795412119 Forward reaction 8 0.0113926379537138 Forward reaction 37 0.0121416499150591 Forward reaction 24 0.0125158192768498 Forward reaction 29 0.0283685658437263 Forward reaction 13 0.0346772420107361 Forward reaction 9 0.0421005353808225 Forward reaction 5 0.0438925336438239 Forward reaction 19 0.052214327287295 Forward reaction 14 0.0559978780179485 Forward reaction 40 0.117011315958426 Forward reaction 3 0.119098963953184 Forward reaction 35 0.134916815372161 Forward reaction 26 0.139922784638101 Forward reaction 21 0.156522233152828 Forward reaction 30 0.159817663946469 Forward reaction 20 0.170817390119866 Forward reaction 39 0.188379388541834 Forward reaction 4 0.216840759493952 Forward reaction 15 0.220017919211627 Forward reaction 42 0.233730982227189 Forward reaction 27 0.289635965581729 Forward reaction 25 0.297660981150591 Forward reaction 16 0.313186383157947 Forward reaction 17 0.408165884440031 Forward reaction 33 0.436800577700256 Forward reaction 41 0.466945029846515 Forward reaction 28 0.46798620169842 Reverse reaction 36 4.31457435977878e-12 Reverse reaction 10 4.07393245251796e-10 Reverse reaction 31 5.46767994862199e-10 Reverse reaction 37 6.94717166789094e-10 Reverse reaction 1 1.62679455171629e-09 Reverse reaction 2 3.36952444783841e-09 Reverse reaction 4 3.9394494587627e-09 Reverse reaction 6 1.07567882185921e-07 Reverse reaction 8 1.67295386056467e-07 Reverse reaction 18 1.80412812846361e-07 Reverse reaction 7 2.29730484442612e-07 Reverse reaction 29 1.02774798716307e-06 Reverse reaction 42 1.59774748772963e-06 Reverse reaction 38 1.60488656493203e-06 Reverse reaction 34 2.78747206613877e-06 Reverse reaction 32 4.79372054090737e-06 Reverse reaction 22 5.15258170490348e-06 Reverse reaction 21 9.09619375620443e-06 Reverse reaction 19 1.6264917798335e-05 Reverse reaction 33 2.2111561567654e-05 Reverse reaction 11 2.30018419933948e-05 Reverse reaction 14 3.91565530202607e-05 Reverse reaction 12 6.04970176547156e-05 Reverse reaction 9 0.000139162041703192 Reverse reaction 3 0.00016265208315078 Reverse reaction 41 0.000350182801894106 Reverse reaction 13 0.000499139457311589 Reverse reaction 20 0.000549518915338207 Reverse reaction 30 0.00069364782602912 Reverse reaction 26 0.000732044791986802 Reverse reaction 5 0.000811514966114664 Reverse reaction 16 0.000822421784788142 Reverse reaction 27 0.00399298992388133 Reverse reaction 15 0.00452951884981948 Reverse reaction 24 0.031077090852003 Reverse reaction 40 0.0404069054905756 Reverse reaction 35 0.0422888795532929 Reverse reaction 39 0.067009846281901 Reverse reaction 23 0.196527603907256 Reverse reaction 25 0.300502412047135 Reverse reaction 28 0.339570121629816 Reverse reaction 17 0.4037115227148 Reaction 11 1.76836114590199e-05 Reaction 1 2.00306712724843e-05 Reaction 2 5.32560412611922e-05 Reaction 12 0.000199370652543748 Reaction 22 0.000203838178516555 Reaction 6 0.00205742513026953 Reaction 18 0.00214112144276709 Reaction 31 0.00256497802215946 Reaction 36 0.00423851339874798 Reaction 38 0.00628440315477067 Reaction 10 0.00709110759390656 Reaction 34 0.00727675619860293 Reaction 7 0.00816841794021552 Reaction 32 0.0105674795412113 Reaction 8 0.011392637948932 Reaction 37 0.0121416498847893 Reaction 25 0.0223800727572329 Reaction 24 0.0282453396822729 Reaction 29 0.0283685654460358 Reaction 13 0.0346772420107361 Reaction 9 0.0420990341989184 Reaction 5 0.0436190317079101 Reaction 19 0.052214327287295 Reaction 14 0.0559978780179485 Reaction 40 0.116797747919766 Reaction 3 0.119098963953184 Reaction 35 0.124498803353466 Reaction 26 0.139190739846114 Reaction 21 0.156522233136118 Reaction 30 0.159814723012322 Reaction 20 0.170817390107279 Reaction 39 0.1880142273317 Reaction 23 0.196350171114342 Reaction 4 0.216840758956933 Reaction 15 0.216904850055577 Reaction 42 0.233730982227189 Reaction 17 0.288086975041949 Reaction 27 0.28850598819427 Reaction 16 0.313186383157947 Reaction 33 0.436800577699977 Reaction 41 0.466945029846514 Reaction 28 0.467986201698403
After this step, by comparing the reference trajectories with the new ones, we choose to delete 16 reactions, so the next step is performed with 14 species and 26 reactions.
Fig2 : Comparison between the trajectories of the target species, the temperature and the flame speed computed with the reduced 14-species 42-reaction mechanism (in black) , and the trajectories computed with the reduced mechanism with 14 transported species and 26 reactions after a DRGEP reaction reduction (in red).
QSS step
The QSS step provides the QSS criteria of each species and the different links with the others species :
Species H 0.0581710482364302 Species O 0.0160759870785582 Species O2 0.277560360139677 Species OH 0.0144373005142915 Species H2O 0.314791474636433 Species HO2 0.010252726793236 Species CH3 0.100336627880574 Species CH4 0.79524736833312 Species CO 0.694064507028327 Species CO2 0.692897407440143 Species HCO 0.00243878401227475 Species CH2O 0.0772557856780877 Species CH3O 0.0207969879289909 Species N2 0
Interactions with species H with QSS Criteria 0.0581710482364302 HO2:2 CH3:1 CH2O:2 CH3O:1 N2:1
Interactions with species O with QSS Criteria 0.0160759870785582 OH:1 HO2:1 CH3:1 CH2O:1 CH3O:1
Interactions with species OH with QSS Criteria 0.0144373005142915 O:1 OH:2 HO2:1 CH3:2 HCO:2 CH2O:1 CH3O:1
Interactions with species HO2 with QSS Criteria 0.010252726793236 H:2 O:1 OH:1 CH3:1 CH2O:1 N2:1
Interactions with species CH3 with QSS Criteria 0.100336627880574 H:1 O:1 OH:2 HO2:1
Interactions with species HCO with QSS Criteria 0.00243878401227475 OH:2
Interactions with species CH2O with QSS Criteria 0.0772557856780877 H:2 O:1 OH:1 HO2:1
Interactions with species CH3O with QSS Criteria 0.0207969879289909 H:1 O:1 OH:1
Interactions with species N2 with QSS Criteria 0 H:1 HO2:1
We choose to put the species CH3O and HCO in QSS hypothesis due to their low QSS coefficient. The species HO2 could have been a good candidate too but it is linked to himself (non linearity).
In order to obtain the trajectories of the 14-species 26-reaction reduced scheme with these 2 species in QSS hypothesis, we run the getQSSfile step and we obtain the following graphs.We observe that the final state and the shape of the trajectories are conserved. It is under these conditions that the optimisation step will be efficient.
Fig3 : Comparison between the reference trajectories of the target species, the temperature and the flame speed (in black), and the trajectories computed with the reduced mechanism with 12 transported species and 26 reactions after a DRGEp reaction reduction (in red).
