查找某物质的热力学数据并转换为chemkin数据格式

在使用CHEMKIN时，其库中的热力学数据有时无法满足需求，需要自行查询物质的热力学数据并编写热力学文件，本文将以B2O3为例子描述查找到转换热力学数据的整个过程。

查找热力学数据

• 进入 NIST的热力学数据库网站 http://webbook.nist.gov/chemistry/form-ser.html
• 输入化学式，使用Calorie-based单位，勾选气相数据，serch
  
• POLY NIST 将热力学参数已多项式拟合系数的方式给出，如下图所示，点击view table
  
• DATA NIST 将热力学参数罗列成表格如图，将表格数据复制
  

在CHEMKIN中转换热力学数据

这种数据是不被chemkin所接受的，需要进行转换。

• 创建一个文本文件，后缀.txt或者.inp都可以
• 内容如下，前面的内容自己根据自己情况进行更改，后面数据就是查到的数据复制的内容

SPEC B2O3(G)
ELEM O 3
ELEM B 2
H298 -199.8000
DATA NIST
PHAS G
TEMP 4000
3000. 25.41 118.5 97.27 63.65
3100. 25.43 119.3 97.97 66.19
3200. 25.46 120.1 98.65 68.74
3300. 25.48 120.9 99.31 71.28
3400. 25.50 121.7 99.96 73.83
3500. 25.52 122.4 100.6 76.39
3600. 25.54 123.1 101.2 78.94
3700. 25.55 123.8 101.8 81.49
3800. 25.57 124.5 102.4 84.05
3900. 25.58 125.2 103.0 86.61
4000. 25.59 125.8 103.5 89.16
4100. 25.60 126.5 104.1 91.72
4200. 25.61 127.1 104.6 94.29
4300. 25.62 127.7 105.2 96.85
4400. 25.63 128.3 105.7 99.41
4500. 25.64 128.8 106.2 102.0
4600. 25.65 129.4 106.7 104.5
4700. 25.65 130.0 107.2 107.1
4800. 25.66 130.5 107.7 109.7
4900. 25.67 131.0 108.1 112.2
5000. 25.67 131.5 108.6 114.8
5100. 25.68 132.1 109.0 117.4
5200. 25.68 132.6 109.5 119.9
5300. 25.69 133.0 109.9 122.5
5400. 25.70 133.5 110.4 125.1
5500. 25.70 134.0 110.8 127.6
5600. 25.71 134.5 111.2 130.2
5700. 25.71 134.9 111.6 132.8
5800. 25.72 135.4 112.0 135.4
5900. 25.72 135.8 112.4 137.9
6000. 25.73 136.2 112.8 140.5
END

• 保存
• 打开chemkin软件，我使用的是CHEMKIN PRO 2022 R1
• Utility中选择RUN Polynominal Equation Fitting
• 设置工作路径，也就是刚刚保存的txt文件的文件路径
• 在输入文件中选择保存的文件
• run fitdat
• 打开文件路径，找到.out文件，使用记事本打开

Species B2O3(G) Fit:
B2O3(G)                 O   3B   2          G  3000.000  6000.000 4000.00      1
 1.18494870E+01 6.07616309E-04-1.32582973E-07 1.31081583E-11-4.67876530E-16    2
-1.03891360E+05-1.85244692E+01 1.26110458E+01-4.20459661E-04 3.52889457E-07    3
-8.44612371E-11 6.67129341E-15-1.04287393E+05-2.29877776E+01                   4

以上数据就是CHEMKIN可以识别的热力学文件。