chemkin 使用C3H8/空气在对冲里燃烧

复现文章的计算部分。
文章链接：

https://www.sciencedirect.com/science/article/pii/S0010218019300380

计算方面：

首先用chemkin pro里的opposed -flow，气相的化学反应机理是KAUST-Aramco_mech1.0,这个机理目前没有。AramcoMech1.3用于基本的化学反应以及添加物的形成苯环的反应。
PAH的形成路径是来自KAUST-PAH-Mechanisms

碳烟模型

碳烟模型包括：

1. 碳烟粒子的成核
2. 通过化学反应以及PAHs之间缩合的表面生长
3. 碳烟粒子与粒子之间的合并
4. 碳烟粒子的氧化

下面首先是：

一、粒子的成核
碳烟粒子的成核是8个大的PAHs,(A4-A8)碰撞形成的。成核反应有不同的反应速率，成核反应一共有36个，他们的反应速率用PAH（与PAH的大小有关）碰撞效率来形容。
成核模型是向着更大的分子还有中等大小的PAH的共价键迈进

https://pubs.rsc.org/is/content/articlehtml/2018/cp/c7cp07803j?page=search

这个过程中共振稳定自由基还有辅助侧链反应目前备受关注；

https://science.sciencemag.org/content/361/6406/997

碳烟粒子一旦成核，那么就会在质量以及大小上增长，质量。这个表面增长是通过反应以及粒子的合并。而表面反应一般就是脱氢加乙炔（HACA）机理，还有氢的提取。

重点来了

模型参数参考下面的文献：

http://xueshu.baidu.com/usercenter/paper/show?paperid=0d6215843c947b5fff0dde736bd41376&site=xueshu_se

但是呢乙炔添加反应的反应速率是不一样的C2H2:C_soot*+C2H2=C_soot-H+H_*（这个反应应该也是要参考论文的）
要和实验里的数据匹配（what?不太懂)

二、物理的PAH聚合碰撞，

包括八个PAH分子与碳烟分子在碳烟颗粒的表面发生碰撞，具体的细节的PAH聚合反应和聚合速率：主要就是PAH的聚合形成碳烟核的过程，参考文献为:

https://www.sciencedirect.com/science/article/pii/S0010218014002636

三、碳烟粒子之间的聚合

碳烟粒子与粒子之间的碰撞效率为2.2

http://xueshu.baidu.com/usercenter/paper/show?paperid=c41eed276c5258771460375e86ac25bd&site=xueshu_se

四、碳烟粒子的氧化

碳烟粒子的氧化考虑了OH和O₂，其中OH的碰撞效率是0.13.
碳烟粒子的和O₂的氧化主要是考虑和O2和A4。引用：

https://www.sciencedirect.com/science/article/pii/S1540748908001235

碳烟粒子的聚合还有碎片化没有考虑

碳烟粒子的动力学

先把文章贴出来

https://www.tandfonline.com/doi/abs/10.1080/02786829708965459

碳烟粒子首先被分成25份，每份是假设是球形并且，每份的密度是1.8g/cm3
模型中最小的粒子是0.878nm，对应的是芘的二聚物（C₁₆H₁₀），
空间系数是2，什么是空间系数？

空间系数就是两个连续的部分的质量体积比

粒子代表性的大小是0.878-224.8nm。碳烟的属性不随着粒子包数的增加而变化
所有的计算用的是mixture-average formulation，计算的快易收敛。
同时由于温度的梯度导致的质量的扩散也包含在内。
一些物质的热辐射也是包含在内的，CO₂，CH₄，CO，H₂O，还有碳烟。
粒子的热塑性也包括在内，为什么包括热泳效应呢？
因为粒子的扩散速度受到热泳的影响。

软件操作

打开chemkin创建新工程

然后选择对冲火焰模块，双击。然后就在右边出现了一个对冲火焰并且还有一个进气，但是呢，还需要一个进气

然后再选择出气就可以了

接下来更新工程

然后选择工作路径和反应机理

如果新拿到一个机理会有四个文件已KM2为例子

这四个例子是直接导不进chemkin的需要进行新建chemkin格式的机理，也就是.cks后缀的文件，很简单，点击new chemistry set 然后在下面的四个文件导入处导入KM2里对应的文件就可以了
然后导入.cks文件
然后点击run-pre processor
然后进行设置

reactor Physical property 页面的设置如下

然后选择求解能量方程，利用的是sectional method。关于这个方法，我不太懂，需要看帮助文件。然后利用了plateau profile for the initial guess，这个也不太懂为什么选这个，不过就是默认设置，，然后论文里明确说了要包含thermal diffusion也就是Soret effect；使用mixture average transport，听师兄说是因为算的快。然后气体的热辐射就没有考虑

grid property 设置

选圆柱坐标系，然后选择specific grid profile，为什么这么选，我也不太清楚，师兄经验。其中那个ending axial position是喷口之间的距离，8cm

气相的碳烟设置（disperse phase:carbon)

basic如下，第一个表面反应的系数是默认为1
使用sectional方法将粒子分为25份，所以number of sections是25
section的空间系数是2

空间系数是指连续的两个section的质量\体积比

section minimum size是指最小部分里的原子数这里取32，我理解的是两个A4的碳原子数，最后碳烟粒子之间的碰撞效率为2.2。

碳烟的Aggregation忽略。
粒子的辐射模型打开，辐射率按照默认的为700

气相的热辐射考虑CO₂，CO，H₂O和CH₄，实际上上边还考虑了碳烟热辐射，如上图。
最后计算