(整理)Chemkin模型学习读书笔记.

Chemkin模型学习读书笔记

一、模型总体介绍

大型气相动力学计算软件包Chemkin(chemical kinetics）可以用来解决带有化学反应的流动问题，是燃烧领域中普遍使用的一个模拟计算工具。该软件是1980 年美国Sandia 国家实验室Kee R. J. 等人开发并推出的，经几次完善发展，至今已开发出了第6个版本CHEMKIN 4.0.2。chemkin有多种针对不同模型的应用程序，在4.0版本中共有23种计算模型，分6大类：

○1封闭的0维反应器：包括封闭的内燃发动机模型(closed internal combustion engine simulator)，封闭的同质反应器(closed homogeneous batch reactor),封闭的部分搅拌反应器(closed partially stirred reactor)和封闭的等离子反应器(closed plasma reactor)。顾名思义，此类模型没有出入反应流，只根据反应器的初状态计算其末状态的参数。

○2开放的0维反应器：包括良搅拌反应器PSR(perfectly stirred reactor)，等离子良搅拌反应器(plasma PSR)和部分搅拌反应器(partially stirred reactor)。此类模型需要定义入流的流量、种类和温度等信息，计算后会给出出口的状态参数。

○3流动反应器：包括栓塞流反应器(plug-flow reactor)、等离子栓塞流反应器(plasma plug-flow reactor)、平面层流反应器(planar shear flow reactor)、圆柱形通道内的层流反应器(cylindrical shear flow reactor)和蜂窝整料反应器(honeycomb monolith reactor)。此类模型考虑流动中的化学反应，主要是表面反应。

○4火焰模拟反应器：包括预混层流燃烧器－稳定的火焰(premixed laminar burner-stabilized flame)、预混层流火焰－火焰速度计算(premixed laminar

flame-speed calculation)、和扩散/预混对撞火焰(diffuseion or premixed

opposed-flow flame)。

○5多晶片沉积滞留CVD反应器：停滞流CVD反应器(stagnation flow reactor)和旋转盘CVD反应器(rotating disk CVD reactor。

○6爆管反应器：通用瞬间爆轰反应器(normal incident shock)和通用反射爆轰反应器(normal reflected shock)。

此外还可独立进行化学平衡和相平衡的计算及对反应机理模型进行分析。

下面就其中一部分模型进行详细介绍。

二、部分模型的详细介绍

1、PSR模型

良搅拌反应器（Perfectly Stirred Reactor，PSR）计算模型常应用于化学反应基础研究中。PSR是连续理想混合流动模型的理想混合反应器，物料以稳定的流量进入反应器后，瞬间就在整个反应器内分散均匀并与器内原存留的物料完全混合，因此反应物转化为生成物的速率由化学反应速率控制而不是混合过程。这样的假设减小了计算强度，容器内的燃烧过程能够用详细化学反应机理来描述。

理想混合流动的特点是：（1）器内以及出口物料的组成和温度等参数均匀一致，且不随时间、空间而变化；（2）各物料微元在器内的停留时间不尽相同，存在停留时间分布。如图1所示。

图1 理想混合式流动

质量控制方程为：

()*0

(1,2,,)k k k k m Y Y VW k K ω--== （1-1）

能量控制方程为： ()**10K

k k k k k m Y h Y h Q =-+=∑ （1-2）

在上述方程中，m

表示反应气体的质量流率(g/s)；k ω 表示第k 种组分的摩尔生成速率(mol/(cm 3s))；Y k 表示第k 种组分的质量分数；W k 表示第k 种组分的摩尔质量(g/mol)；V 是搅拌器的容积(cm 3)；h k 是第k 种组分的比焓(J/g)；Q 是反应器的散热速率(J/s)。上标*表示进口参数。

滞留时间（Residence time ）τ由搅拌器的体积和气体的质量流速确定：

V m

ρτ= (s) （1-3） 其中质量密度ρ通过理想气体状态方程计算得到：

PW RT

ρ= （1-4） 上述公式中P 表示压力；R 是通用气体常数(J/mol ⋅K)；T 是温度(K)；W 是混合物的平均摩尔质量(g/mol)。

2、Plug-flow 模型

Plug-flow 是连续理想排挤流动模型，理想排挤反应器，物料以稳定的流量由反应器的一端流入反应器后，各物料微元沿流动方向齐头并进，完全没有轴向混合与扩散，就好像活塞在气缸里向前平推一样，如图2所示。

图2 理想排挤式流动

活塞流模型的特点是：（1）在与流动方向垂直的任意截面上各点，物料的流速、浓度、温度及停留时间等完全一样；（2）物料的浓度、温度等各参数沿流动方向递变；（3）在每一截面上物料各参数都不随时间而变化。

在管段流中，已被接受的理想模型是PFR ，其假设没有混合在轴向方向但在横向混合充分。它也被展示成没有轴向可实现的反应物转化最大。同样地，缺少横向梯度意味着没有质传递限度，再一次加强反应器性能。除了这些实用的优点外，PFR 反应器计算非常快捷，因为它是使用初命令ODE ’s 模型，不需要传递其他特性。

调节PFR 反应器性能方程式就是简化质传递、能量传递、动量传递的一半模式。他们可以衍生出很多在流向中微元段的平衡，规定在横向上没有变量，任何质量的轴向扩散相比于相关对流项都被忽略。这样，所有的气体质平衡（连续性方程）可表达为：

质量连续方程：

∑∑===++g K

k k m k M

m m i W s a dx d uA dx du A dx dA u 1,1, ρρρ （2-1） 这里ρ是质浓度，u 是轴向气流速，包括K g ；W k 是k 的摩尔重量，S k 表示所有的表观方应器中这种物质的摩尔产率，质量A 和a i,m 是截面面积和材料单位长度的有效内表面积。A 和a i,m 可作为x 的随机函数。式2-1简单描述了气体质流速在反应器中可以作为生成或者消耗结果反应所有物质。相似方程也可以独立列出。

气相转化方程：

⎪⎪⎭

⎫ ⎝⎛+=+∑∑∑===g g K k k m i m k k K k k m k M m m i k k A a s W W s a Y dx dY uA 1,,1,1,ωρ （2-2） 这里Y k 是物质k 质量分数，ωk 是同类气体反应的摩尔产率。这些反应不能改变总的气体质量，但可改变组分。

能量平衡：

k K K k k m k M m m i e e K k k m k M m m i K k k k P K k k k h W s a q a W s a u Y h dx du u dx dT C dx dY h uA g b g g g ∑∑∑∑∑∑======-=⎪⎪⎭⎫ ⎝⎛++⎪⎪⎭⎫ ⎝

⎛++',1,1,1,12121 ρ（2-3）

h K 是物质k 的比热焓，是单位气体质量的平均热容，T 是相对气体温度。在正确的求和中s k,m 是固体物质k 在物质临界面的体积摩尔产率。大量的与临界的物质区别在1.1章和第4章被讨论。式2-3描述了流动气体总能量变化引发从外墙环境的热通量变化Qe ，和散落固体的焓聚集。如果式2-3没包含临界面物质的焓，那么这个公式就无任何意义。

气体的动量方程表示压力、惯性的、粘滞力和临界反应动量之间的平衡。因而，动量方程：