2.查找某物质的热力学数据并转换为chemkin数据格式2023-05-06

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

1|0查找热力学数据

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

2|0在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可以识别的热力学文件。

__EOF__

本文作者：宋培荣
本文链接：https://www.cnblogs.com/bbhnqxx/p/17376613.html
关于博主：评论和私信会在第一时间回复。或者直接私信我。
版权声明：本博客所有文章除特别声明外，均采用 BY-NC-SA 许可协议。转载请注明出处！
声援博主：如果您觉得文章对您有帮助，可以点击文章右下角【推荐】一下。您的鼓励是博主的最大动力！