耦合ch基础介绍

1、CHENKIN基础介绍2009.9.主要内容o CHEMKIN介绍o CHEMKIN用户界面o CHEMKIN化学设置o CHEMKIN模型o CHEMKIN后处理o CHEMKIN算例化学模拟提品和服务o 市场范围:nnnnnn汽车能源航空化工材料微电子CHEMKIN应用于工程问题o设计方向nnnn改进反应器的影响工业过程的参数运行范围确定过程的可替代性反应过程时间尺度上的可性o反应器优化和改进nnn工况运行条件变化的影响过程参数扰动的敏感性 反应器的生产能力评估o厂商nn排放过程变化对下游的影响o模型假设和验证nn化学机理的发展与简化实验结果的CHEMKIN有效的求解问题o快速模型的发展n

2、nn分级近似能够把不同的理想反应器模型连接在一起不需要建立网格o有效准确的求解技术nnn良好的刚性方程组求解方法高效的体系结构继续计算和重新开始计算的功能o计算结果容易分析和解释nn内置的数据可视化数据管理的敏感性分析CHEMKIN 4.xo 基于java的操作界面。o 操作界面不再是基于计算模块，而是反应器模型o 全新的后处理功能o 项目图形化o 全新的参数研究功能o 颗粒追踪模块主要内容o CHEMKIN介绍o CHEMKIN用户界面o CHEMKIN化学设置o CHEMKIN模型o CHEMKIN后处理o CHEMKIN算例探究用户界面o 更改用户参数选择nnnn检查 “培训” 和“例子

3、” 文件的位置设置“记事本” 作为 缺省编辑器设置IE作为 HTML 阅读器测试 Adobe Acrobat Reader 作为 PDF 阅读器o 打开 开始手册o 为“线性沉积率”设置缺省为 “microns/min”介绍图形界面o 注意:输出打开反应器 要求一个 输入 和一个GoodBad研究图像选项1.2.3.4.5.6.7.8.创建一个新方案创建一个PSRs与一个或多个输入的网在末端处添加一个管流反应更新方案试验在参数项间连线在 PSRs间增加热流并再次更新打开前处理面板为化学设置选择 grimech30.cks 文件 (据的化学设置系统数前处理程序的化学设置由工作路径可以把化学设置文

4、件放在不同的地方继续运行前先运行前处理程序工作路径是运行一个方案输入和输出文件的地方在保存的上方，化学设置的定义保存在一个选择“*.cks” f文件中反应输入面板了输入手册中“”将鼠标放在上面可以提供快速的输入信息只需要填入必须的信息; 剩下空白的为缺省可以局部地选择每个输入的输入面板包括参数和反应物成分输入流柱参数包括速率和温度/氧化剂平衡率可以输入反应物的成分用按氧化剂混合物定义氧化剂的分数,它们的和必须等于1.0按总的反应混合物定义附加的反应物分数, 它们的和必须 <1按混合物定义分数; 它们的和必须等于 1.0平衡率和氧化剂的相对量平衡率定义o 例如: CH4作为, O2 作为氧

5、化剂o CH4完全燃烧的平恒方程:CH4 +O2 ® CO2 +H2Oo 方程给出了CH4:O2的化学当量比 (1:2)Reactant Fractionso 平衡率与下面的有关:j = ( CH4 /O2 = CH4 /O2 CH4 /O2 )stoich0.5øO2CH40.91.00.451.01.00.51.11.00.551.21.00.61.31.00.65有平衡率选项,内在的执行这些计算22求解面板:通常可以设成缺省在输入手册中详细地描述了高级选项调整相对和绝对残差能够改进收敛性继续和新的运行选项o 继续选项n 稳态:某个计算结果可以用于另一个初始设定n 瞬态

6、: 从上一个求解结果开始新的计算o 包括 distance-marching 模块 (PFR, 剪切流反应器)o 新的运行选项nn开始运行每个符合要求的参数研究方法(pamameter study),即在同一个模型中设置不同的变量,并且同时处理继续 / 新的运行面板”为每一个继续(或新运行)通过指定新的参数值给运行增加 “案例选择继续或新的运行选项输出面板着火延迟定义选项敏感性分析选项运行面板当前正在运行程序; 同时显示一个监视器项面板在运行前为反应模型创建一个“”输入文件补充输入, 如元素, 可以增加到输入文件中指定文件名称当运行完全, 弹出图像后处理程序存档选项o Archive生成一个与

7、方案相关的压缩文件,含输入和输出文件nnn容易将方案总体转移创建与以前有关的工作存档文件将方案发给同事o unarchive生一个方案设置你在一个新的工作目录下再n不依赖原始文件路径 n 便于 RD 技术支持主要内容o CHEMKIN介绍o CHEMKIN用户界面o CHEMKIN化学设置o CHEMKIN模型o CHEMKIN后处理o CHEMKIN算例一个“化学设置”由描述体系的两个或的文件热力学数据文件 (必需)onn适合于热力学参数的多项式系数适用于许多反应物的预包装数据气相化学文件 (必需)onn元素、 气相反应物、气相反应对于用户感的问题必需或获得输运特性文件 (有时必需)onn适

8、合于气相反应物的基本参数适用于许多反应物的预包装数据表面化学文件 (可选)on表面状态、点反应物、固体反应物、表面反应n对于用户感的问题必需或取得热力学数据o数据是适应于温度的求解反应物特定热、焓和熵的多项式系数0cpT 2 +T 3 +T 4=+T +RH 0a2a3a4a5a6 234=a1 +T +TTTRTS 02345Ta3a4a5a7 234= a1 ln T + a2T +T+T+T+R234T输入是固定的格式(相关 NASA平衡代码)on灵活使用许多新扩充的o热力学数据包括反应物的元素成分o在反应机理文件中用户可以随意地包含反应物热力学数据a5a4a3a2a1热力学数据输入的例

9、子反应物名字符号反应物适用的温度界限!Species: Name: Source: Comment: H0(298K)AL2H6CASNumber: 12004-30-7DimericAluminum Trihydride,SNL fit to data generated from Pollard fit, 6/29/87R. Pollard, J. Crystal Grow., V.77, P.200 (1986)=21.3500 (Kcal/mole),S0(298K) =62.7500 (cal/mole-K)一个反应物的数据62987G12342.63488400e+008.8713

10、4600e+03-1.11958200e-072.13595200e-02 3.15415100e-07-7.68467400e-092.33583200e-121.03974700e-059.82751500e+00-6.80068100e+008.45915500e-11 1.06053700e+045.08074400e-025.55452600e+01a aAL2ME662987AL2C6H18G300.0001500.000600.00123417适用于 1st适用范围-2.19321200e-071.64414400e-10-3.51554600e+043.89076300e+01

11、(高温)注释提供了反应物的信息1.77314700e+01 4.93574700e-02 1.19685400e-06-1.63982600e-08 4.89086700e-12-3.85556000e+04-5.05329800e+01-7.15975000e-01 1.06710900e-01 2.11760500e-05! Species: AL2ME6CAS Number: 15632-54-9! Name:Trimethylaluminum, Dimeric! Source: SNL fit to data generated from Pollard fit, 6/29/87! C

12、omment: R. Pollard, J. Crystal Grow., V.77, P.200 (1986)! H0(298K) =-61.2000 (Kcal/mole), S0(298K) =131.050 (cal/mole-K)300.000 1500.000600.00AL 2H6AL2H6气相反应机理例子反应物名字必须包含在热力学数据中ELEMENTS HON OENDSPECIESH2 H O2OHHO2H2O2H2O N N2 NO ENDREACTIONSH2+O2=2OH OH+H2=H2O+H O+OH=O2+H O+H2=OH+H H+O2+M=HO2+MH2O/1

13、8.6/ OH+HO2=H2O+O2 H+HO2=2OH O+HO2=O2+OH 2OH=O+H2O0.170E+140.117E+100.400E+150.506E+050.361E+18 N2/1.26/0.750E+130.140E+150.140E+140.600E+090.001.30-0.502.67-0.72477803626062900!D-L&W JAM 1986 KLEMM,DIXON-LEWISH2/2.86/0.000.000.001.300107310730!D-LD-LD-LCOHEN-WEST.END固定系数: A, B, E反应 H2/空气火焰 3 种元

14、素 11 种反应物 23 个反应反应速率常数假定Arrhenius 表达式o 修正的Arrhenius 表达式有三个参数:expæ - Ei ök= ATç÷R TfiièøcoE的是 calories/mole (缺省)n 为全部或一个反应指定其他选项:o kcal/mole, Joule/mole, Kelvin, or eV (电子伏特)Bi反应速率常数假定Arrhenius 表达式oA的依赖于反应的类型必须是 mole/cm3-snn切记过程速率的u '依赖于KÕkXkk =1可以是 (moles,cm,s)

15、的组合 (缺省)n输入或 (molecules, cm, s)可逆反应逆反应速率系数o详细平衡的方法k fik=riKci平衡常数用浓度n 平衡常数源自于反应物的热力学参数:标准大气压Kuæ DS°DH ° öæ Påkik =ö1= expçi -i ÷ç atm ÷Kciç÷èRRTRTøèø反应的焓变反应的熵变o用户可以不考虑特征并明确的规定逆反应速率系数n 辅助反应数据:REV/ ArevBrev Erev/输运数据输

16、入的例子反应物名e/k136.500136.50071.40097.530232.400232.400209.000209.000224.700417.000s3.3303.3303.2983.6213.8283.8284.1004.1004.1623.690m0.0000.0000.0000.0000.0000.0000.0000.0000.0001.700a0.0000.0000.0001.7600.0000.0000.0000.0000.0000.000Zrot0.0000.0000.0004.0001.0001.0002.5002.5001.0002.000线性0001111122AR

17、AR* C C2 C2O CN2 C2HC2H2 C2H2OH CH2OHLennard-Jones直径极性转动松驰因子Lennard-Jones深度偶极力矩表面反应机理的例子SITE/SILICON/SDEN/1.66E-9/SI(S)ENDBULKSI(B)/2.33/ ENDTHERMO ALL300.600.1685.SI(S)J 3/67SI1000000000S300.0001685.000123412340.24753989E 01 0.88112187E-03-0.20939481E-060.42757187E-11 0.16006564E-13-0.81255620E 03-

18、0.12188747E 02 0.84197538E 000.83710416E-02-0.13077030E-040.97593603E-08-0.27279380E-11-0.52486288E 03-0.45272678E 01SI(B)J 3/67SI 100 000 000 0S300.000 1685.0000.24753989E 01 0.88112187E-03-0.20939481E-06 0.42757187E-11 0.16006564E-13-0.81255620E 03-0.12188747E 02 0.84197538E 00 0.83710416E-02-0.13

19、077030E-040.97593603E-08-0.27279380E-11-0.52486288E 03-0.45272678E 01 ENDREACTIONSSIH4 SI2H6 SIH2 SI2H2 2SI2H3 H2SISIH2 2SI2H52SIH32SIHSI H3SISIH SI2 SI3+SI(S)SI(S)SI(S)2SI(S)4SI(S)2SI(S)4SI(S)2SI(S)2SI(S)SI(S)2SI(S)2SI(S)3SI(S)=>=>=>=>=>=>=>=>=>=>=>=>=>SI(S)2S

20、I(S)SI(S)2SI(S)4SI(S)2SI(S)4SI(S)2SI(S)2SI(S)SI(S)2SI(S)2SI(S)3SI(S)+SI(B) + 2H21.05E174.55E263.9933E111.7299E206.2219E371.7007E206.1186E372.3659E202.4465E204.1341E111.7007E201.7607e208.6586E280.50.50.50.50.50.50.50.50.50.50.50.50.540000400000.00.00.00.00.00.00.00.00.00.00.02SI(B)SI(B)2SI(B)4SI(B)2S

21、I(B)4SI(B)2SI(B)2SI(B)SI(B)2SI(B)2SI(B)3SI(B)+ 3H2+H2 H2 3H22H25H23H2H2+2H2END表面反应物的热力学数据通常包含在表面动 力学输入文件中（不是“热力学数据中”）化学设置文件概要o 前处理程序创建连接文件，连接文件将指定的化学信息传递到反应模型ü始终要求(ü ) 有些问题要求反应类型气相动力学热力学数据表面动力学输运参数数据平衡üü(ü)0-D 封闭 PSRsüü(ü)管流üü(ü)PaSRüü

22、激波管üü剪切通道流üü(ü)ü1-D 火焰 模拟器üüüCVD 反应器üüüü机理分析üü(ü)(ü)LPCVD 反应炉üüüüLPCVD 热化学文件来源获得的 CHEMKIN 机理o从其他研究nhtml.CHEMKIN用户群:nn从文献， NIST 编辑系数oNIST 动力学数据库搜索 CHEMKINmlnn数据库:o由计算化学来计算系数Ab initio 方法动力学结合和聚合添加方

23、法nnno基于理论和经验的估计n类推到其它的化学体系RD公司: 机理集成服务o对于你的问题必须确实化学设置的正确性o研究数据源谁开发的?什么目的及开始条件 ?nn确定了解:onnn开发者介绍的应用范围用什么数据作的验证或校对确认研究中使用了什么热力学数据和输运数据o如果可能使用这个数据o敏感性分析来确定主要的反应关注主要反应的路径来源nn了解模型的不确定性主要内容o CHEMKIN介绍o CHEMKIN用户界面o CHEMKIN化学设置o CHEMKIN模型o CHEMKIN后处理o CHEMKIN算例真实的理想化近似H2O CO2YkH2 N2O2RPCONOCH4I2time化学和相平衡反

24、应器完全扰动反应器平推流反应器均质批反应器扩散火焰内燃机多数模块可包含表面化学YkRPItime柱塞和剪切流反应器完全扰动反应器xr转动圆盘 CVD 反应器LPCVD 热分析器和反应炉反应模块概括不发生化学反应的气体混合器化学反应和相平衡计算器机理分析器内燃机闭式均相反应器闭式 PaSR闭式等离子反应器完全扰动反应器 (PSR)等离子 PSR部分反应器 (PaSR)管流反应器蜂窝陶瓷反应器等离子体管流反应器平面剪切圆管剪切流预混燃烧器火焰速度计算器扩散或预混火焰相对滞流 CVD 反应器旋转盘 CVD 反应器入射激波反射激波混合标志外部源初始化出口 ( 4.1 新增)Gas flow split

25、ter ( 4.1 新增)平衡计算器o计算相和化学平衡n包括凝相¶G¶x= 0o完全基于反应物的热力学参数反应范围n不要求动力学数据ooo绝热火焰温度Chapman-Jouguet结果处理:CO2N2绘制结果 vs变参数图H2H2OnO2CO NCH4O2 n描述热力学循环Gibbs 能Gibbs 自由能的最小化闭式0-D 反应器o内燃机模型R容积 vs 扫过时间均匀混合物压缩点火PnnnYkIo一般 “间歇式反应器”模型nnn热传递选项限定容积或气体和表面化学timeReaction #4o闭式等离子模型Reaction #2nnn电子能量方程指定能量沉积气体和表面化学R

26、eaction #18o闭式部分扰动反应器time n限定容积或 n跟踪混合和动力学时间尺度Sensitivity开式0-D反应器开式等离子PSR反应器oo完全搅拌器nnnnnnn不涉及化学计算热力学状态稳态或瞬态模拟可成群气体和表面化学电能和等离子能沉积表面化学的离子碰撞能稳态或瞬态模拟可成群nno一般完全扰动反应器(PSR)o开式部分扰动反应器(PaSR)nnn气体和表面化学稳态或瞬态模拟可成群n湍力学相互作用开式0-D反应器群开式0-D反应器可形成的多反应器“群”o反应器间质量/热循环同时求解所有反应器, 提供封闭耦合部分扰动反应器不可用nnno其它反应器模型是单反应器“群”n PFRs

27、, 剪切反应器, CVD 反应器, 火焰o一个 “反应器网” 可包括一个或多个群nnn连续的群连接一个群外面不逆向循环使用连接或信息连接开式0-D反应器群o多样输入流n不同组分、温度、流率ooo质量和热循环用户指定循环分数和路径反应器间的热传递n 对流/ 传导和/或辐射Inlet BR31Inlet CR13123Inlet AOutletR32R111-D 塞式反应器o 一般平推流反应器nnn在横向忽略扩散方向假定反应器是均衡的RPYkI热损失和表面化学面积o 蜂窝陶瓷反应器n 几何上的活性表面面积& 水力直径距离o 等离子平推流反应器nn电子能量方程能量沉积轮廓每距离近似作为表面面

28、积全部包括敏感性分析、热传递选项、气相和表面化学平衡、0-D或塞式模型为复杂问题建立平台o 系统简单o 研究时间量程nnnn物理直觉较容易复杂因素很少停滞时间传输时间vs动力学时间o 研究化学问题o 利用结果比较好的初始nnnn忽略传输作用平衡vs动力学敏感性分析产物速率分析另外, 火焰计算收敛n1-D 剪切反应器o 圆柱剪切反应器nn圆柱坐标包含径向扩散o 平面剪切反应器nnn平面坐标包含横向扩散对称或非对称选项1-D 火焰反应器o 预混燃烧器nnnn和氧化剂预混合稳定燃烧火焰层流热的形成和反应速率敏感性分析1-D 火焰反应器o 预混火焰速度计算器nnnnn绝热火焰速度确定可燃限和火焰厚度

29、热损失和辐射传递的作用层流热的形成和反应速率敏感性分析1-D 火焰反应器o 对流扩散火焰类似减少2-D 平面或3-D 轴对称流场为1-D 模型的转化1-D边界条件假设扩散火焰可燃限为产生CFD平台进行火焰模拟热的形成和反应速率敏感性分析nnnnnShower-head CVD 反应器旋转盘x3-D的近似1-D转换:on 全部3-D流场n 轴向分量& 温度轮廓旋转盘或滞止 瞬态和稳态选项表面沉积 (在x=0处)rooo滞止xrShower-head CVD 反应器o 输运选项nn均匀或多组分混合热扩散选项o 假定nnn径向均匀边界条件确定径向范围没有到反应器壁面的径向损失标准激波模型o

30、入射激波模型nnnn激波波动前用户指定条件相关气体动力决定初始后激波条件瞬态动力学在激波通过后包括粘性作用o 反应激波模型nn反射激波后相关气体动力决定条件反射激波通过后开始瞬态动力学入射激波反射激波2125U2U2UrsUs机理分析器o 计算化学设置的初始信息化学机理热力学数据传输数据反应速率vs温度逆向反应速率平衡常数nnn反应速率.o 计算热力学参数vs温度热力和传输参数.nnn焓熵 比热机理分析器o 计算传输参数nnnn作为温度的函数粘性扩散系数热传导率指定“bath”气体on提取随而定的反应速率单独增加的两个专业微电子反应器可用o 分层式烘炉低压化学气相沉积(LPCVD) 模型n L

31、PCVD 热力分析模型o温度贯穿反应器n LPCVD 烘炉模型ooo利用热分析模型结果通过薄膜的沉积速率通过负载的沉积速率o 通过CHEMKIN中的可用菜单nn命令输入形式的用户菜单主要内容o Reaction Design 介绍o CHEMKIN介绍o CHEMKIN安装o CHEMKIN用户界面o CHEMKIN化学设置o CHEMKIN模型o CHEMKIN后处理o CHEMKIN算例CHEMKIN 4.1 新的可视化选项3-D 轮廓图交互式旋转镜像图对称坐标多重Y坐标首先运行一个例子问题来产生一个解算方法例子:在一个发中临界压缩比的确定 压缩比在着火发生处,其它所有参数固定 临界压缩比

32、 (CCR)可用于估算» 基于CCR和RON或MON间的关系的辛烷值(ON)» 用来探究混合物和替代的辛烷值 这个例子研究了丁烷-空气混合物(富氧)的CCRCHEMKIN IC 发模型模拟了一个发缸内的压燃汽封闭系统, 瞬态模拟进气门关闭到排气门打开前的一段时间容积 vs. 时间由以下决定 :LC/LA比压缩比 = 总容积 / 压缩容积总容积发速度其它输入:开始曲轴转角 (进气门关闭)初始热损失和温度现象:依靠动力学HCCI着火, 柴油发SI 燃烧敲缸速压机实验上已经观测了CCR和RON间的关系 RON = Research Octane Number 基于计算出的临界压缩

33、比可以调配出理想值的研究法辛烷数据来自Westbrook及其他人:* = 试验值其它符号 = 计算值例子:临界压缩比1.打开 training 案例:“trainingic_engine_ignition_compression_ratio.ckprj”观察前处理面板化学设置并运行前处理观察反应器面板设置2.3.模拟时间=0.04 sec (或 360 度曲轴转角)在参数研究中指定压缩比变化4.在 RunParameter 面板中, 选择 All cases 和 Run Selected更新运行状态在RumParameter Study 面板中输出文件可以被存取RumParameter Stu

34、dy 后处理结果创建温度 vs. 时间和 vs.压缩比的三维图5.选择后处理变量在图中改变显示的过滤器控制输入变量包含反应物的生产速率作为输出变量弹出后处理程序后处理程序面板在同一时间里可以加载多重解算文件镜像图用来显示对称结果可以应用以前保存的 “宏”不同图形的类型选项 (用户可以创建新的类型)创建图形HAND S ON用温度作为Z轴绘制一个 3-D 轮廓图放大/缩小旋转图: 文注释连续撤消(Cntrl-Z),对于所有操作!本框, 图画形状, 箭头记号临界压缩比 13.2用OH摩尔分数作为Z轴绘制一个 3-D轮廓图OH分数 vs时间、压缩比的轮廓图OH是燃烧的中心很容易输出绘制的图像以几种格

35、式输出图像 (bmp, gif, tiff, etc.)图像的清晰度和剪切/粘贴到其它的 Windows 应用程序为几个CCR值绘制一个温度 vs时间的线性图通过按下ctrl键进行多项选择数据区域工具快速地放大感区域的点击图激活数据区域点击数据区域图标缩放区域画距形框或者用箭头-/+ 来调整区域有很多选项可以图像格式和外观双击图像进入Visualization Browser窗口数据空间 (顶部): 轴缩放比例,三维轮廓线/线性: 颜色, 符号颜色, 线型, 宽度背景色,开/关 外形, 尺寸, 颜色, 填充,Z轴坐标轴: 线型, 符号开/关, 正文字体 , 记号注释: 图例文本字体 , 间距,

36、 轮廓改变字体，坐标轴和图例文字，线型并增加注释双击图像进入Visualization Browser在图的Window工具条中使用Annotation图标可以保存并在以后应用用户定义的图像类型方便使用设置首选字体轴外观线形 & 符号风格注释为将来绘图保存一个类型用 Edit ® Style ® Create Style 然后选择, 或 Edit®Style ® Styleeditor其他图像工作可由“Macros”自动完成为其它操记录并应用“macros” 可以任何鼠标点击或键盘敲击面板，Save/Reapply在后处理（保存/再运用）选项可以

37、重复案例Save / Reapply 绘图设置选项 选项可以一个接一个地保存工程可以后处理部分的数据，并可以输出到文本文件在一个单独的窗口所选择的数据设置为应用到其它程序如Microsoft Excel 或Tecplot,输出数据到一个文本文件用Excel画数据图 的例子温度vs时间图（输出到Excel）60Pressure (atm)807050403020CR=12.5 CR=12.9 CR=13.1 CR=13.2 CR=13.3 CR=13.6CR=14100150250350450550650Crank angle (degrees)数据也可以很容易地输入到CHEMKIN后处理中从例

38、如Microsoft Excel输入逗号/制表符/空格分隔符文本文件不核查柱标题在面板可利用新的绘图设置主要内容o CHEMKIN介绍o CHEMKIN用户界面o CHEMKIN化学设置o CHEMKIN模型o CHEMKIN后处理o CHEMKIN算例模拟的基本步骤o 确定建模的方法n能否用理想的反应器模拟实际过程o 建立反应器（群）o 寻找合适的机理，并进行预处理o 确定问题类型并给定初始条件、几何条件、物理参数o 创建输入文件，执行程序o 通过后处理程序分析结果问题求解技术o 对于任何问题的定义，可能要求几种方法来找到解决办法0-D模型适合于化学研究或者证实系统的效果氧化剂高级模型适合于

39、设计和处理优化问题rx对于复杂几何, 模拟通常要求分级进行o 总的化学描述o 从工程概念来设计完全3-D CFD模拟x1-D M(类似转换)化学复杂性几何复杂性CHEMKIN®CFD/KINeticschem.inptherm.datchem.asc绝热火焰温度o H2/airo 项目名：equilibrium_gas.ckprjo Reactor：equilibriumo Chem.inpn 3种元素：H、O、Nn 9种物质：H2、H、O2、O、OH、HO2、H2O、N2、H2O2o Therm.dat绝热火焰温度o 等压等焓o 初始温度：300K：1atmoo 预估温度解：2000Ko 反应物质：n N2：3.76，H2：2.0，O2：1.0o 连续计算：改变初温n 400K和500K稳态燃烧o 大气压下H2、N2和O2混合物的燃烧，燃烧温度随平衡比的变化o条件：n 温度：298Kn 平衡比：=1.0成分（摩