Chemkin 输入文件的规则.docxVIP

CHEMKIN的文件规则输入文件的一些通用规则：符号“！”无论出现任何位置，此行后面的文本将作为注释文本而被忽略。②输入文件每行不应超过80个字符（这个和Fortran读取txt文件每行的长度是一样的）③除了个别有极其严格规则的地方外（如热力学数据的定义等，均会特殊声明），空格作为分隔符，而且多个空格将被视为一个。④数字格式：可以为整数（如99）、浮点数如（99.99）、或E格式（如9.999E2，E大小写均可）。编写Gas-Phase Kinetics的输入文件。该文件包括四部分的内容：元素、组分、热力学数据、化学反应，如下例所示：! 例：Gas-Phase Kinetics输入文件ELEMENTS H O END! 元素定义SPECIES H2 H O2 O OH H2O END! 组分定义THERMO! 热力学数据(本例中只重新定义了“OH”的热力学数据)OH 121286O1H 1 G 0300.00 5000.00 1000.0010+02 0-02-0-050-09-0-14 20+05 0+02 0+02 0-02-0-05 30-07-0-11 0+05 0+01 4ENDREACTIONS! 反应方程及Arrhenius系数H2+O2=2OH 0.170E+14 0.00 47780OH+H2=H2O+H0.117E+10 1.30 3626O+OH=O2+H 0.400E+15 -0.50 0O+H2=OH+H 0.506E+05 2.67 62902OH=O+H2O 0.600E+09 1.30 0H+H+M=H2+M 0.100E+19 -1.00 0H2O/0.0/ H2/0.0/! 辅助数据H+H+H2=H2+H2 0.920E+17 -0.60 0H+H+H2O=H2+H2O 0.600E+20 -1.25 0END①元素(Elements)［规则］此部分以ELEMENTS（或者ELEM，两者等价）关键字开头；其后以空格为间隔符列出将在反应中出现的所有元素；最后以END关键字结束。例：ELEMENTSH O END! 元素定义用户要按周期表的元素名（两个字母均须大写）来定义元素，CHEMKIN可以辨认的元素如表2.1。如果用户定义同位素或新元素的话，可以以1~2个字母命名（不与元素表重名），并将其原子量以“/”括住跟随其后。如定义氢（H）的同位素氚（命名为HH）如下：ELEM HH / 3.0 /END表2 CHEMKIN的元素表H, HE, LI, BE, B, C, N, O, F, NE, NA, MG, AL, SI, P, S, CL, AR, K, CA, SC, TI, V, CR, MN, FE, CO, NI, CU, ZN, GA, GE, AS, SE, BR, KR, RB, SR, Y, ZR, NB, MO, TC, RU, RH, PD, AG, CD, IN, SN, SB, TE, I, XE, CS, BA, LA, CE, PR, ND, PM, SM, EU, GD, TB, DY, HO, ER, TM, YB, LU, HF, TA, W, RE, OS, IR, PT, AU, HG, TL, PB, BI, PO, AT, RN, FR, RA, AC, TH, PA, U, NP, PU, AM, CM, BK, CF, ES, FM, D, E其中，D为氢（H）的同位素氘（D），E为电子，如果化学反应中有离子参加，电子必须作为一个元素进行定义。②组分（Species）［规则］此部分以SPECIES（或SPEC）开头；其后以空格为分隔符列出将在反应中出现的所有组份；最后以关键字END结束。例：SPECIES H2 H O2 O OH H2O END! 组分定义CHEMKIN的热力学数据库therm.dat自带的组分有（778种）。我们要使用这些组分的热力学数据时，组分的命名就必须按照CHEMKIN数据库的命名，这样CHEMKIN才可以从数据库中自动获取该组分的元素组成及其热力学性质。当然，我们可以用自己定义组分（不超过16个字符且以字母开头），但组分的热力学数据同样需要自行定义。③热力学数据（Thermodynamic Data）［规则］此部分以THERMO开头；其后以空格为分隔符列出定义组份的名称和热力学数据；最后以关键字END结束。例：THERMO! 热力学数据