[求教]敏感性分析

sigou

在output control里选择了all rate of production,极限设为1.0e-4，输出的标准化的敏感性分析结果有正数有负数，请问应怎样选择敏感性大的反应？我的目的是想去掉一些不重要的反应得到骨架机理。

sigou

中科大的同仁，那位给小妹指点一下下。。。

zyw

Chemkin里好像没有这中自动功能吧？除非你自编程，不过既然都得到敏感性系数的大小了，选择其骨架反应应该不是很难。另外不知道是你写错了还是真的弄错了，敏感性系数应该是sensitivity，而rate of production是产率！

eric

各反应机理的敏感系数随反应过程是变化的，也就是说，“骨架机理”也是变化的。华中理工大学煤燃烧国家重点实验室好像也在做这方面的工作。

sigou

谢谢zyw。
chemkin4.0在output control里可以选输出敏感性分析结果。有三种选择，我现在试着选了rate of production，输出结果里多了对每一种物质生成率影响最大的前N个(N由设置的threshold决定)基元反应，影响有多大用正交标准化后的一个值来表征。
把对每种物质影响较大的基元反应收集在一起，去掉重复的就可以做到简化了吗？chemkin的基元反应都有A,n,E的，简化后的基元反应的A,n,E要不要变呢？听说浙大也有人做，按不变的来，我们导师觉得可能不太对，让我找一下资料倒底对不对，我还没找着，不知该找哪儿。

sigou

就是这个样子的结果，每个组分有一个。threshold设得大一点则列出的反应就少一点。
NORMALIZED AND ABSOLUTE RATE-OF-PRODUCTION COEFFICIENTS
               (surface reactions are normalized by the material surface area to volume ratio)
               Threshold normalized value for printing of ROP coefficients =  0.100E-01

  Gas Phase Species   1. H2                                                                  NORMALIZED  (MOLES/CC-SEC)
          Gas Phase Reaction   3.                                        O+H2<=>H+OH        -0.303     (-7.6106E-08)
          Gas Phase Reaction  41.                                    2H+H2O<=>H2+H2O         0.014     ( 7.1678E-10)
          Gas Phase Reaction  45.                                      H+HO2<=>O2+H2         0.082     ( 4.2158E-09)
          Gas Phase Reaction  49.                                        H+CH<=>C+H2         0.015     ( 7.5416E-10)
          Gas Phase Reaction  53.                                     H+CH4<=>CH3+H2         0.205     ( 1.0511E-08)
          Gas Phase Reaction  55.                                      H+HCO<=>H2+CO         0.019     ( 9.7214E-10)
          Gas Phase Reaction  58.                                    H+CH2O<=>HCO+H2         0.189     ( 9.7093E-09)
          Gas Phase Reaction  75.                                   H+C2H4<=>C2H3+H2         0.023     ( 1.2000E-09)
          Gas Phase Reaction  84.                                      OH+H2<=>H+H2O        -0.697     (-1.7507E-07)
          Gas Phase Reaction 126.                                      CH+H2<=>H+CH2         0.086     ( 4.4208E-09)
          Gas Phase Reaction 146.                                  CH2(S)+H2<=>CH3+H         0.034     ( 1.7310E-09)
          Gas Phase Reaction 174.                             C2H4(+M)<=>H2+C2H2(+M)         0.011     ( 5.6258E-10)
          Gas Phase Reaction 214.                                      HNO+H<=>H2+NO         0.049     ( 2.5046E-09)
          Gas Phase Reaction 284.                                     O+CH3=>H+H2+CO         0.206     ( 1.0579E-08)
          Gas Phase Reaction 288.                                    OH+CH3=>H2+CH2O         0.018     ( 9.4075E-10)
                     Total Rate-of-production  (moles/cc-sec) =  5.1278E-08
                     Total Rate-of-consumption (moles/cc-sec) =  2.5117E-07
                     Net   Rate-of-production  (moles/cc-sec) = -1.9989E-07