Chemkin 软件模拟低压乙烯预溷火焰.docxVIP

Chemkin模拟低压乙烯预混火焰——吴武华 S、Chemkin简介CHEMKIN软件是美国Sandia国家实验室开发的大型气相化学反应动力学软件，可以用来解决带有化学反应的流动问题，是燃烧领域中普遍使用的一个模拟计算工具。CHEMKIN是一种非常强大的求解复杂化学反应问题的软件，常用于对燃烧过程、催化过程、等离子体及其他化学反应的模拟。CHEMKIN以气相动力学、表面动力学、传递过程这三个核心软件包为基础，提供了对20多种常见化学反应模型及后处理程序。三个核心程序模块为：1) 气相动力学（Gas-Phase Kinetics）：是所有程序计算的基础，包括气相成分组成、气相化学反应与相关的Arrhenius数据等信息。2) 表面动力学（Surface Kinetics）：多相反应中，如催化反应、化学气相沉积等，Surface Kinetics提供两相反应所需的各种信息，如表面结构、表面成分组成及热力学数据、表面化学反应等。3) 传递（Transport）：提供气相多组分粘度、热传导系数、扩散系数和热扩散系数等。其中Surface Kinetics和Transport必须以Gas-Phase Kinetics为基础，因为它们中出现的成分都必须在Gas-Phase Kinetics中已定义。Gas-Phase Kinetics、Surface Kinetics和Transport提供了化学反应的基本信息，生成动态链接库供后续程序调用。如下图给出了 CHEMKIN软件模拟层流预混火焰时的计算结构。反应解释器首先读取化学反应机理和热力学数据，输出的连接文件中包含了机理中元素、物种和反应的信息，并进入输运解释器与输运数据一起计算出各物种和体系平均的粘度系数、扩散系数、热传导率和热扩散系数等传热传质参数。最后由反应程序库与输运程序库利用输入参数求解能量方程、动量方程、浓度矩阵等燃烧方程。二、乙烯火焰的模拟2.1 实验条件：表1乙烯火焰的实验条件实验压力/Torr喷嘴孔径/μm当量Ar流量/SLMO2流量/SLMC2H4流量/SLM304001.70.7501.1170.6331502401.70.7501.1170.6332.2 模拟部分在此，我所使用的版本为CHEMKIN-PRO版本，下面，以模拟30 Torr压力下乙烯预混火焰为例，简述乙烯火焰的模拟。模拟步骤：新建工程新建一个工程，取名为C2H4 wwh01。反应器模型选择在models面板中选择Premixed Laminar Burner-Stabilized Flame 模型，添加出口气体反应器后，点击Update project按钮。设定反应机理在Open Projects面板，双击Pre-Processing，进入机理设置面板。在Gas Phase Kinetic File中，选择“USTC(11.11.21) .inp”文件。在Thermodynamics Data File中，选择“ThermUSTC.DAT”文件。在Gas Transport Data File中，选择“TrandatUSTC.DAT”文件。注：以上文件是我们实验室自己编写的文件。气相反应动力学文件“USTC(11.11.21) .inp”是由丁烯火焰机理修改而来，文件中，包含176种物种，642个反应。设定完之后，点击“Run Pre Processor”。运行成功后可进行初始参数的设定。设定初始参数面板中，点开Pre-mixed_Burner（C1），进入初始参数设定界面。双击C1_Pre-mixed Burner，设置炉子的参数。在Reactor Physical Properties栏中，Problem Type选择“Fix Gas Temperature”，Temperature Profile中，点击Browse图标按钮，选择温度曲线文件“C2H4 Phi_2.2 - 2 -93-temperature.ckprf”（温度曲线为实验中测出的）。Pressure中输入30 Torr。在Grid Properties栏中，Ending Axial Position中，输入3 cm。其他参数默认。双击面板中C1_Iniet1，设置进气参数。在Stream Properties Data中，入口流量输入0.0041 g/cm2-sec。Species-specific Properties中，添加进气的反应物种的摩尔分数：Ar气，0.3；C2H4，0.2532；O2，0.4468。其他参数默认。双击面板中的Solver，对求解器设置。其中，Basic栏设置收敛条件，绝对误差限和相对误差限，大小可以参看chemkin帮助说明。Advanced栏设置当收敛遇到问题时，时间作出让步