复现文章的计算部分。
文章链接:
https://www.sciencedirect.com/science/article/pii/S0010218019300380
计算方面:
首先用chemkin pro里的opposed -flow,气相的化学反应机理是KAUST-Aramco_mech1.0,这个机理目前没有。AramcoMech1.3用于基本的化学反应以及添加物的形成苯环的反应。
PAH的形成路径是来自KAUST-PAH-Mechanisms
碳烟模型
碳烟模型包括:
- 碳烟粒子的成核
- 通过化学反应以及PAHs之间缩合的表面生长
- 碳烟粒子与粒子之间的合并
- 碳烟粒子的氧化
下面首先是:
一、粒子的成核
碳烟粒子的成核是8个大的PAHs,(A4-A8)碰撞形成的。成核反应有不同的反应速率,成核反应一共有36个,他们的反应速率用PAH(与PAH的大小有关)碰撞效率来形容。
成核模型是向着更大的分子还有中等大小的PAH的共价键迈进
https://pubs.rsc.org/is/content/articlehtml/2018/cp/c7cp07803j?page=search
这个过程中共振稳定自由基还有辅助侧链反应目前备受关注;
碳烟粒子一旦成核,那么就会在质量以及大小上增长,质量。这个表面增长是通过反应以及粒子的合并。而表面反应一般就是脱氢加乙炔(HACA)机理,还有氢的提取。
重点来了
模型参数参考下面的文献:
但是呢乙炔添加反应的反应速率是不一样的C2H2:Csoot*+C2H2=Csoot-H+H*(这个反应应该也是要参考论文的)
要和实验里的数据匹配(what?不太懂)
二、物理的PAH聚合碰撞,
包括八个PAH分子与碳烟分子在碳烟颗粒的表面发生碰撞,具体的细节的PAH聚合反应和聚合速率:主要就是PAH的聚合形成碳烟核的过程,参考文献为:
https://www.sciencedirect.com/science/article/pii/S0010218014002636
三、碳烟粒子之间的聚合
碳烟粒子与粒子之间的碰撞效率为2.2
四、碳烟粒子的氧化
碳烟粒子的氧化考虑了OH和O2,其中OH的碰撞效率是0.13.
碳烟粒子的和O2的氧化主要是考虑和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,对应的是芘的二聚物(C16H10),
空间系数是2,什么是空间系数?
空间系数就是两个连续的部分的质量体积比
粒子代表性的大小是0.878-224.8nm。碳烟的属性不随着粒子包数的增加而变化
所有的计算用的是mixture-average formulation,计算的快易收敛。
同时由于温度的梯度导致的质量的扩散也包含在内。
一些物质的热辐射也是包含在内的,CO2,CH4,CO,H2O,还有碳烟。
粒子的热塑性也包括在内,为什么包括热泳效应呢?
因为粒子的扩散速度受到热泳的影响。
软件操作
打开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
气相的热辐射考虑CO2,CO,H2O和CH4,实际上上边还考虑了碳烟热辐射,如上图。
最后计算