CHEMKIN学习提高入门基础指南.doc
《CHEMKIN学习提高入门基础指南.doc》由会员分享,可在线阅读,更多相关《CHEMKIN学习提高入门基础指南.doc(18页珍藏版)》请在得力文库 - 分享文档赚钱的网站上搜索。
1、CHEMKIN 4.0.1 入门指南入门指南燃烧学燃烧学 1辅助教程辅助教程文中如有任何错误,敬请指出,以便不断改进;如有任何文中如有任何错误,敬请指出,以便不断改进;如有任何问题,欢迎提出,共同探讨问题,欢迎提出,共同探讨助教博士生:卢智恒助教博士生:卢智恒联系方式:热能系系馆办公室联系方式:热能系系馆办公室 201(O)62782108(H)2005.3一、一、CHEMKIN 的安装和简介的安装和简介11 安装 CHEMKINChemkin 最早的版本始于 1980,由美国 Sandia 实验室的 Kee RJ 等人编写,经过多年 的不断发展日趋完善。后来由 Reaction Design
2、 公司收购并继续开发,目前最新版为 4.0.1。 由于学习和科研需要,我们花费 12000$向 ReactionDesign 公司订购了一套最新版本的 CHEMKIN 4.0.1,其中包括可供 20 人同时在线计算的 license,用于燃烧学课程的学习。【安装安装】请登录 ftp:/combustion:combustion166.111.56.155 下载相关文件,其中 chemkin401_pc_setup.exe 为 CHEMKIN 的安装程序,chemkin.lic 为网络认证文件,详细的 安装信息可以参看 ftp 上的“安装说明.txt”文件。安装完后会自动在桌面及开始菜单建立 快
3、捷方式。 【注意】1、本套教学用的 CHEMKIN 软件采用网络认证的方式,请确保电脑已经联 网(校内) ,否则无法计算。 2、建议采用 1024768 的分辨率,否则某些界面将无法完全显示。12 CHEMKIN 简介CHEMKIN 是一种非常强大的求解复杂化学反应问题的软件包,常用于对燃烧过程、 催化过程、化学气相沉积、等离子体及其他化学反应的模拟。CHEMKIN 以气相动力学、 表面动力学、传递过程这三个核心软件包为基础,提供了对 21 种常见化学反应模型及后处 理程序。 三个核心程序模块为: 1) 气相动力学(气相动力学(Gas-Phase Kinetics):是所有程序计算的基础,包括
4、气相成分组成、 气相化学反应与相关的 Arrhenius 数据等信息。 2) 表面动力学(表面动力学(Surface Kinetics) 。很多反应过程包括多相反应,如催化反应、化学 气相沉积、固体腐蚀等。在这些反应里,Surface Kinetics 提供两相反应所需的各种信息, 如表面结构、表面和体内的成分组成及热力学数据、表面化学反应等。 3) 传递(传递(Transport) 。提供气相多组分粘度、热传导系数、扩散系数和热扩散系数等。其中 Surface Kinetics 和 Transport 必须以 Gas-Phase Kinetics 为基础,因为它们中出现 的成分都必须在 Ga
5、s-Phase Kinetics 中已定义。 Gas-Phase Kinetics、Surface Kinetics 和 Transport 提供了化学反应的基本信息,生成动 态链接库供后续程序调用。用户可以自己编写程序调用它们来实现指定的功能,但最方便 的是使用 CHEMKIN 自带的反应模型,共分 7 类,包括 21 个反应模型和 2 个应用程序, 将在下一节介绍。二、二、CHEMKIN 的简单使用入门的简单使用入门21 CHEMKIN 的用户界面CHEMKIN 的用户界面如图 1 所示。除了传统的 Windows 菜单和按钮以外, CHEMKIN 的工作界面分为任务栏、窗口栏、消息栏三部
6、分。图 1 CHEMKIN 4.0.1 的用户界面 其中任务栏包括 Open Projects 和 Models 两个选项卡,Models 选项卡中列出了 CHEMKIN 提供的 21 种反应模型和 2 种应用程序,分 7 类排列,如表 1 所示。有关 Open Projects 选 项卡的使用将在下文的例子中介绍。 表 1 CHEMKIN 提供的 7 类 21 种反应器与 2 种应用程序External Source of Inlet Gas,添加入口气体源,一般在有多路 气体输入时才使用Non-Reactive Gas Mixer,计算不反应气体组分的混合过程Chemical and Ph
7、ase Equilibrium Calculation,计算化学平衡和 相平衡Miscellaneous (杂项)(杂项)Mechanism Analyser,分析气相和表面化学反应系统中的热化 学、传递及动力学数据Closed Internal Combustion Engine Simulator,模拟封闭的内 燃机内的燃烧情况Closed 0-D Reactors (封(封闭闭 0 维反应器)维反应器)Closed Homogenous Batch Reactor,模拟封闭的全混同性反应器,包括定压、定容反应器Closed Partially Stirred Reactor(PaSR),
8、模拟封闭的部分混合 的反应器Closed Plasma Reactor,模拟封闭的等离子体反应器Perfectly Stirred Reactor(PSR),模拟稳态及瞬态的充分混合反 应器Plasma PSR,模拟充分混合等离子体反应器。Open 0-D Reactors (开(开口口 0 维反应器)维反应器)Partially Stirred Reactor(PaSR),模拟部分混合反应器Plug Flow Reactor,一维柱塞流反应器Plasma Plug Flow Reactor,等离子体一维柱塞流反应器Planar Shear Flow Reactor,平板层流反应器Cylind
9、rical Shear Flow Reactor,圆柱层流反应器Flow Reactors (流动反应器)(流动反应器)Honeycomb Monolith Reactor,蜂窝结构反应器Premixed Laminar Burner-stabilized Flame,模拟层流预混的 稳态火焰Premixed Laminar Flame-speed Calculation,层流预混火焰传 播速度的计算Flame Simulators (火焰模拟)(火焰模拟)Diffusion or Premixed Opposed-flow Flame,模拟扩散或预混 的对冲火焰Stagnation Flow
10、 CVD Reactor,模拟用于化学气相沉积的滞流 反应器CVD Reactors (化学气相沉(化学气相沉积反应器)积反应器)Rotating Disk CVD Reactor,模拟用于化学气相沉积的转盘流 反应器Normal Incident Shock,模拟入射激波的化学动力学Shock Tube Reactors (激(激波管道反应器)波管道反应器)Normal Reflected Shock,模拟反射激波的化学动力学22 CHEMKIN 的求解过程1、Gas Phase Kinetics(气相动力学)的处理(气相动力学)的处理 Gas Phase Kinetics 的前处理器(Pr
11、e-processor)读取用户编写的气相动力学输入文件和自 带的热力学数据库(therm.dat),生成包含元素、组分、热力学数据反应信息的 Gas-Phase Kinetics 连接文件。Gas-Phase Kinetics 提供子程序库处理该连接文件。 2、Surface Kinetics(表面动力学)和(表面动力学)和 Transport(传递过程)的处理(传递过程)的处理 如果化学反应包含表面反应或传递过程,则需要相应地执行这两个核心程序块。 Surface Kinetics 的前处理器读取用户编写的表面动力学输入文件,生成包含表面反应 信息的 Surface Kinetics 连接
12、文件,Surface Kinetics 提供子程序库处理该连接文件。Transport 的前处理器根据 Gas-phase Kinetics 连接文件中的信息,自动从 CHEMKIN 自 带的传递数据库(tran.dat)读取相应的数据,然后生成包含传递信息的 Transport 连接文件, Transport 提供子程序库处理该连接文件。 3、反应模型求解、反应模型求解 根据问题需要,CHEMKIN 读取模型输入文件确定求解方法。在 4.0 及其以上版本的 CHEMKIN 中,新加入了通过相应的模型设置窗口中设置模型的参数的功能。用户设置完 成后点击 Create Input File 即可
13、生成模型输入文件,然后用户可以通过 View Input File 按钮 查看输入文件的内容。而对于 4.0 以前的版本,用户需要通过手动编辑模型输入文件。但 对于文件的内容,都是相同的,都是采用关键字的形式声明模型的功能调用和参数设置。 有关常用的平衡计算模型和全混反应器模型的关键字列于附录 2、3。 完成所有设置后,即可 Run Model 进行计算,CHEMKIN 会自动调用上述 Gas-Phase Kinetics、Surface Kinetics、Transport 各自的子程序来读取反应信息,调用模型输入文件控 制模型求解过程。程序计算结束后,会生成一数据文件 xxxxxx.out
14、 供用户查阅数据,以及 一动态连接文件 XMLdata.zip 供后处理(绘图)使用。 4、后处理(、后处理(Post-Process) CHEMKIN 提供了统一的后处理器,用于对应用程序的结果进行分析和绘图。23 CHEMKIN 的用户操作步骤下面结合实际例子介绍用户操作的步骤。 【例例 1】计算化学当量的计算化学当量的 H2 与空气的定压绝热燃烧温度。与空气的定压绝热燃烧温度。 1运行 CHEMKIN,点击菜单 Project-New,输入项目名称,这里我们定为 H2-air。 2. 决定问题的性质,选择适当的反应模型。由于绝热燃烧温度的计算是一个相平衡过 程,不涉及具体的反应过程,于是
15、我们很容易地想到用 Chemical and Phase Equilibrium Calculation 模型。点击任务栏中的 model 选项卡,点击 Chemical and Phase Equilibrium Calculation 图标,此时在窗口栏的 Diagram View 窗口中将看到新加入一个相平衡计算的 模型,最后点击窗口右下角黄色的 Update Project 按钮,如下图所示。通常情况下,一般的化学反应问题通过适当的假设和简化,都可以对应到某一种 CHEMKIN 包含的反应模型,有时一个问题还可以有多种选择。选择恰当的应用程序是求 解问题的第一步。但是如果实在不幸没有一
16、种模型可以很好的解决你的问题,或者你想要 更完美地解决一些问题,就只有自己编写程序,调用 CHEMKIN 里的子程序库进行计算了。3此时左侧任务栏会自动切换至 Open Projects 选项卡,双击 Pre-Processing,窗口栏 出现 Pre-Processing 的参数窗口。在窗口中的 Working Dir 一项中填入你希望的保存路径, 或者通过右侧的 Browse 按钮点击选取。然后按 New Chemistry Set 按钮,点击 Gas-Phase Kinetics Files 项右端的编辑按钮,如下图所示。在弹出的窗口中选择刚才的工作路径,输入文件名 chem.inp,按
17、 Open/Create 按钮后即 可开始编辑气相动力学输入文件了。-【CHEMKIN 的文件规则】- CHEMKIN 的输入文件有它自己的规则,用户在编辑输入文件的时候应该遵守这些规 则。在介绍 Gas-Phase Kinetics 输入文件之前,先介绍输入文件的一些通用规则: 注释符号“!” 。符号“!”无论出现任何位置,此行后面的文本将作为注释文本 而被忽略。 输入文件每行不应超过 80 个字符 除了个别有极其严格规则的地方外(如热力学数据的定义等,均会特殊声明) ,空 格作为分隔符,而且多个空格将被视为一个。 数字格式:可以为整数(如 99) 、浮点数如(99.99) 、或 E 格式(
18、如 9.999E2,E 大小写均可) 。 下面介绍如何编写 Gas-Phase Kinetics 的输入文件。该文件包括四部分的内容:元素、 组分、热力学数据、化学反应,如下例所示: ! 例:Gas-Phase Kinetics输入文件 ELEMENTS H O END! 元素定义 SPECIES H2 H O2 O OH H2O END! 组分定义 THERMO! 热力学数据(本例中只重新定义了“OH”的热力学数据) OH 121286O 1H 1 G 0300.00 5000.00 1000.00 1 0.02882730E+02 0.10139743E-02-0.02276877E-05
19、 0.02174683E-09-0.05126305E-14 2 0.03886888E+05 0.05595712E+02 0.03637266E+02 0.01850910E-02-0.16761646E-05 30.02387202E-07-0.08431442E-11 0.03606781E+05 0.13588605E+01 4 END REACTIONS! 反应方程及Arrhenius系数 H2+O2=2OH 0.170E+14 0.00 47780 OH+H2=H20+H0.117E+10 1.30 3626 O+OH=O2+H 0.400E+15 -0.50 0 O+H2=OH
20、+H 0.506E+05 2.67 6290 2OH=O+H2O 0.600E+09 1.30 0 H+H+M=H2+M 0.100E+19 -1.00 0 H2O/0.0/ H2/0.0/! 辅助数据 H+H+H2=H2+H2 0.920E+17 -0.60 0 H+H+H2O=H2+H2O 0.600E+20 -1.25 0 END 元素元素(Elements)规则此部分以 ELEMENTS(或者 ELEM,两者等价)关键字开头;其后以空格为间隔符列出将在反应中出现的所有元素;最后以 END 关键字结束。例:ELEMENTS H O END! 元素定义 用户要按周期表的元素名(两个字母均须
21、大写)来定义元素,CHEMKIN 可以辨认的 元素如表 2.1。如果用户定义同位素或新元素的话,可以以 12 个字母命名(不与元素表重 名) ,并将其原子量以“/”括住跟随其后。如定义氢(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,
22、 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 为电子,如果化学反应中有离子参加,电子 必须作为一个元素进行定义。 组分(组分(
23、Species)规则此部分以 SPECIES(或 SPEC)开头;其后以空格为分隔符列出将在反应中出现的所有组份;最后以关键字 END 结束。例:SPECIES H2 H O2 O OH H2O END! 组分定义 对于我们所要用到的组分,都可以从 CHEMKIN 的热力学数据库 them.dat 找到。 CHEMKIN 数据库自带的组分名称(778 种)列入附录 1 中。我们要使用这些组分的数据, 组分的命名就必须按照 CHEMKIN 的规则来,这样 CHEMKIN 才可以从数据库中自动获取 该组分的元素组成及其热力学性质。 当然,我们也可以用不超过 16 个字符而且以字母开头的任意字符串来
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- CHEMKIN 学习 提高 入门 基础 指南
限制150内