甲烷燃烧的数值模拟
微燃烧室内甲烷燃烧的数值模拟
Ｖ Ｏ ． Ｎ Ｏ．５ ２０１ １ ４ｌ ２
童
加
熟
Ｉ ＮＤＵＳ ＴＲＩ ＡＬ ＨＥＡＴＩ ＮＧ
１ ９
阜｝乍｝乍｝乍｝ｅ｝年毫譬｝譬 － 亡碌
誊 热能研 究 ｌ
｝ ｝ ｝ ｛ 毒 毒 亨 兔
ＤＯ ： １ ． ６ ￣ ｓｎ ０ ２１ ３ ． １．５０ ８ Ｉ ０３ ９ ．ｉ ．１ ０ —６ ９２ ２０ ． ９ ｓ ０ ０
１１ 燃 烧器 结构 ．
微燃烧器采用的燃烧 室长为 １ １ｍ，宽为 ６ ０Ｔ ｉ ｍｍ，高 为 ０ ．ｍｍ，计算 中采用均匀网格，如 图 １ ８ 所示 。
和二氧化碳的质量百分 比浓度分布 ，图 ３给 出了不同流
量下微燃烧 室内的温度分布 。对应不 同的流量 ，初始 条 件相 同的 甲烷和 氧气的混合气在 微燃烧室 内的燃烧情 况
Ｋｅ ｙ ｗｏｒ ｄｓ： ｍ ｉｒ ｃ ｃ ｏ ｏｍ ｂ ｔ ｒ ｃ ｍ ｂｕ ｔｏ ｕｓｏ ； ｏ ｓｉ ｎ； ＦＬＵ ＥＮ Ｔ ｏｆｗ ａｅ；ｎｕ ｅ ｉａ ｉ ｕ ａ ｉｎ ｓ ｔ ｒ ｍ ｒｃ ｌｓｍ ｔ ｔｏ
自２ ０世纪 ８ ０年代以来， 微机 电系统 （ ＭＳ ｍｉ ｏ ＭＥ ： ｃ — ｒ ｅ ｃ ｏｍｅｈｎｃｌｙｔ ）迅速发展，在航 空、航天 、汽 ｌ ｔ — ｃａ ｉ ｓ ｍｓ ｅｒ ａｓ ｅ
＿ ） ＋） ｝Ｅ 丧［ Ｅｐ＝ （＋ （ 】 ｐ ｐ
ａｆ ａ Ｔ ，
一
ｌ （ ３ ）
＋Ｒ ｉ
组分守恒方程 ：
（１ ＝ 一 ｐ， ）
气体状 态方程 ：
＝ ｐ （） ５
式 中：Ｐ 为密度；ｆ 为时间；Ｕ 方 向的速度；Ｘ 方 向坐 ｉ 为ｉ ｉ 为ｉ
基于fluent的甲烷燃烧二维模拟分析
基于fluent的甲烷燃烧二维模拟分析陈飞 1434422(同济大学汽车学院，上海)摘要: 目前，由于环境污染和排放法规的日趋严格，降低排放已经成为了汽车工业的重点，而寻求清洁的替代燃料是一种较为有效的解决办法。
甲烷属于可再生气体燃料，可以实现与空气的良好预混，利用fluent进行甲烷燃烧的数值模拟进行仿真分析。
Fluent提供了涡耗散模型用于求解燃料快速燃烧，整体反应速率由湍流控制的单步或双步总包燃烧反应。
根据模拟结果分析甲烷作为车用替代燃料的可行性。
关键词: 替代燃料；燃烧的数值模拟；甲烷燃烧；fluent仿真1. 引言燃烧是燃料与氧化剂之间的发光发热的化学反应，根据反应前各组分的分布，可以分为预混燃烧，扩散燃烧和部分预混燃烧。
其中预混燃烧较多的应用于汽车工业的车用汽油发动机。
目前，由于环境污染和排放法规的日趋严格，降低排放已经成为了汽车工业的重点，而寻求清洁的替代燃料是一种较为有效的解决办法。
1.1. 燃烧的数值模拟燃烧的数值模拟是通过CFD软件实现对实际燃烧过程的仿真模拟，求解流畅流动特性及其混合特性，温度场、组分浓度场以及颗粒和污染物排放等，从而提供实际燃烧过程的参考，对于产品研发，科学研究都有很大的意义。
燃烧的数值模型主要运用模拟软件根据燃烧模型进行仿真，目前可用于燃烧数值模拟的软件有FLUENT,STAR-CD,CHEMKIN,KIVA等。
燃烧模型主要根据不同燃烧的特点设置求解参数，包括如下内容：稀疏相模型、输运控制方程、燃烧模型、辐射换热模型、污染物模型。
Fluent提供了涡耗散模型用于求解燃料快速燃烧，整体反应速率由湍流控制的单步或双步总包燃烧反应。
其中对于反应r中的物质i的产生速率由下面两个式子给出：（1.1）（1.2）式中，——任何一种产物的质量组分；——某种产物的质量组分；——经验常数4.0；——经验常数0.5。
1.2. 甲烷性质介绍甲烷在自然界的分布很广，甲烷是最简单的有机物，是天然气，沼气，坑气等的主要成分，俗称瓦斯。
甲烷二次燃烧火焰传播速度的数值模拟研究
甲烷二次燃烧火焰传播速度的数值模拟研究史增凯;席文雄;金星;张玉坤【摘要】二次燃烧是常见热力设备和高速推进燃烧室内的重要现象,但目前对二次燃烧的研究较为匮乏.为了研究当量比、甲烷添加量对二次燃烧自点火火焰传播速度的影响情况,论文利用CHEMKIN软件中火焰速度反应器模拟研究了向不同当量比(0.3～0.9)甲烷/空气一次燃烧产物中添加不同量(摩尔分数0.02～0.24)甲烷时的自点火火焰传播速度.研究表明当量比主要是通过影响尾气温度和一次燃烧产物中富余氧气来影响火焰传播速度,而甲烷添加量会影响二次燃烧时局部当量比,局部当量比在1附近时火焰传播速度最大.%Secondary combustion is a common phenomenon in thermal equipment and high-speed propulsion combustion chamber, but the current studies on the secondary combustion are relatively deficient. In order to study the effects of equivalence ratio and added methane on the flame propagation velocity of secondary flame, this paper simulates the combustion process of secondary combustion at different equivalence ratios (0. 3～0. 9) and different molar fractions of added methane (0. 02～0. 24) by using the flame velocity reactor of CHEMKIN, and achieves the flame propagation velocity. The results show that the equivalence ratio affects the flame propagation velocity mainly through influencing the temperature and remaining oxygen of the primary combustion products. Moreover, the amount of added methane can affect the local equivalence ratio of secondary combustion, and the maximum of flame propagation velocity occurs when the local equivalence ratio is around 1.【期刊名称】《机电产品开发与创新》【年(卷),期】2018(031)002【总页数】3页(P52-54)【关键词】二次燃烧;自点火;当量比;火焰传播速度【作者】史增凯;席文雄;金星;张玉坤【作者单位】航天工程大学激光推进及其应用国家重点实验室,北京 101416;航天工程大学激光推进及其应用国家重点实验室,北京 101416;航天工程大学激光推进及其应用国家重点实验室,北京 101416;航天工程大学激光推进及其应用国家重点实验室,北京 101416【正文语种】中文【中图分类】V231.20 引言二次燃烧自点火是常见热力设备和高速推进燃烧室内的重要现象，在能源、材料、推进领域中有着广泛的应用[1～3]。
甲烷燃烧的数值模拟
甲烷燃烧的数值模拟及分析主要分为三个部分，第一部分讲解甲烷在炉膛内燃烧的模型建立的方法；第二部分对甲烷燃烧模型的数值模拟结果进行分析和比较；第三部分为结论。
一、模型建立1、在Gambit中建立计算区域在本例中建立圆柱形炉膛，并研究甲烷和空气在炉膛内的燃烧反应。
物理模型如下：甲烷入口直径为10mm；空气入口直径为50mm；炉膛为直径为500mm；长度为1200mm的圆柱形。
如图1。
图1圆柱形炉膛模型图2、绘制网格图2进口网格分布甲烷和空气进口的网格元素选择四边形网格，如图2。
炉膛表面的网格也是四边形网格，如图3。
图3炉膛表面网格分布图4炉膛表面网格分布图5炉膛出口网格分布图6炉膛内部网格分布3、指定边界条件图7炉膛边界条件Inlet1为甲烷入口，边界条件为速度入口；Inlet2为空去入口，边界条件为速度入口；Outlet为炉膛出口，边界条件为自由流；其他炉膛壁面为墙体，边界条件为墙体。
4、导入fluent具体信息如下：54440mixed cells,zone2,binary.326quadrilateral wall faces,zone3,binary.1900quadrilateral wall faces,zone4,binary.350quadrilateral wall faces,zone5,binary.218quadrilateral outflow faces,zone6,binary.204quadrilateral velocity-inlet faces,zone7,binary.18quadrilateral velocity-inlet faces,zone8,binary. 108880triangular interior faces,zone10,binary.11144nodes,binary.11144node flags,binary.缩放信息如下图：图8缩放信息图5、选择计算模型图9定义求解器图10考虑能量方程图11考虑粘性模型图12考虑辐射模型图12考虑燃料模型图13燃烧物质和炉膛材料6、操作环境的设置图14操作环境（压力场和重力场）7、设置边界条件图15空气入口边界条件空气入口的速度为8m/s，温度为300K，入口空气中氧气的含量为21%。
甲烷燃烧温度随过剩空气系数变化数值模拟
甲烷燃烧温度随过剩空气系数变化数值模拟李振;李佳璇【摘要】采用Fluent软件的ED燃烧模型对不同过剩空气系数α下圆柱空腔内甲烷燃烧进行二维稳态模拟,得出其温度场的变化规律.圆柱空腔长2m,直径为0.6m,保持燃烧器额定功率230 kW不变,即保持燃气的喷射速度不变,保持燃气与空气的温度不变,使α在1.05～1.80范围内变化.结果表明:α对于甲烷燃烧的最高燃烧温度具有明显的影响,最高燃烧温度随α增大呈线性逐渐降低;改变α后最高燃烧温度均出现在x=1.2m左右,在此之前燃烧温度呈二次函数增长,在此之后燃烧温度呈线性缓慢减小并保持平稳,这表明α只影响炉膛内最高燃烧温度的值及高温区的范围,不能影响炉膛内温度的变化趋势.【期刊名称】《煤气与热力》【年(卷),期】2019(039)001【总页数】5页(P38-42)【关键词】过剩空气系数;圆柱空腔燃烧;温度场;二维稳态数值模拟【作者】李振;李佳璇【作者单位】山东省冶金设计院股份有限公司,山东济南250101;山东建筑大学热能工程学院,山东济南250101;山东济华燃气有限公司,山东济南250101【正文语种】中文【中图分类】TK224.11 概述目前，除了实验方法研究燃烧，还有一种比较好的方法，就是利用计算流体力学进行数值模拟。
计算流体力学(CFD)是在计算机上求解描述流体运动、传热和传质的偏微分方程组，并且对上述现象进行过程模拟。
CFD可用来进行流体力学的基础研究、复杂流动结构的工程设计，了解在燃烧过程中的化学反应，分析实验结果等[1-2]。
其基本特征是数值模拟和计算机实验。
Fluent是用于模拟具有复杂外形的流体流动以及热传导的计算机程序，是一种比较常用的CFD软件。
Fluent软件对预混及非预混等各种复杂的燃烧问题有比较好的研究效果，尤其是内置的化学反应模型自诞生以来一直占据着非常重要的地位[3-6]。
2 圆柱空腔燃烧的物理模型本数值模拟原本是三维几何模型，但是，考虑到以下原因：该三维模型为轴对称图形且具有旋转轴，Fluent软件提供了Axisymmetric Swirl的二维代替三维模拟的计算方式，二维模拟可以生成质量更高、数量更少的计算网格，我们采用二维对称轴旋转成三维的几何模型代替三维几何模型，获得了同等精度范围的模拟结果。
Fluent 模拟燃烧
3
混合状态
反应机制
甲烷在空CH4 +3O2 =2CO+4H2O 2CO +O2 = CO2
甲烷完全燃烧 甲烷不完全燃烧
模拟过程中，假设燃料完全燃烧成CO2和H2O
流动条件
甲烷在空气的燃烧
层流
· 各项参数（速度等）稳定 · 低雷诺数
湍流
· 局部参数脉动 · 高雷诺数
甲烷在空气的燃烧
混合状态
非预混火焰：
有限速率化学反应
求解过程中采用的方程为涡耗散模型
访谈结果与析
☞ 模拟结果
燃烧器内，甲烷从开始点燃到趋于稳定过程中温度的变化
☞ 模拟结果
空气：0.5m/s，300k
■
甲烷 :80m/s,300k
甲烷含量监测点
☞ 模拟结果
监测点处，甲烷浓度的变化值
总结与分析
实例概述
图中所示为甲烷火焰燃烧器，
主要用于处理污水厌氧过程中 产生的沼气. 甲烷燃烧器多为圆柱型，甲烷 从中间喷口进入.
模型建立
模拟 计算 区域
空气：0.5m/s，300k
甲烷 :80m/s,300k
网格模型
中间区域及左侧喷嘴附近的区域 在计算过程中需要较密的网格
☞ 模拟机理
1
反应机制
2
流动条件
甲烷燃烧器的优化
• 燃烧器尺寸的优化
不足
• 模拟过程中将燃烧器的桶壁考虑成 绝热，计算过程中，器内部温度要 大于实际过程中的温度.
•
进气速度的确定
谢谢！
微燃烧器内甲烷催化燃烧的数值模拟
吸收受颗粒大小 、 颗粒密实程度等多种因素的影响 , 要实现煤质的在线分析 ,还有大量工作要做 。 参考文献 :
量的能量储备 ,还是在能量转换与利用效率上都明 显地优于传统的化学电池 。例如 , 碳氢燃料提供的 储能密度约为 40～50 MJ / kg ,而目前普遍使用的 、 性 能最好的锂离子电池只有 0. 4 ～ 0. 5 MJ / kg 。因此 , 即使燃料的化学能只有 10 %转变为电能 , 碳氢燃料 提供的储能密度也超过电池 10 倍以上 。正是由于 碳氢燃料的廉价和高储能密度的优点 , 微型发电系 统和微型推进系统的研究已经成为最近几年美国国 防和能源高技术研究的重点投资之一 。麻省理工学 院、 加州大学伯克利分校 、 普林斯顿大学以及南加州 大学等分别在美国 NASA 、 国防部 ( DOD ) 和能源部 (DOE) 的支持下进行了微型火箭 , 微型透平机和以 燃烧方式直接发电等动力微机电系统 ( Power MEMS) 的研 究 , 并 取 得 了 初 步 的 成 果 。日 本 也 把 动 力 MEMS 的研究列入新能源开发机构 ( NEDO) 的研究 计划之中 。 基于燃烧的动力 MEMS 中 , 如何实现微尺度条 件下燃料的持续稳定燃烧是制约系统成败的关键因 素之一 。在微尺度燃烧中 ,由于燃烧器尺度的减小 ,
k =1 j =1 K s N +N
g
2 物理模型和数学模型
计算中采用的物理模型是一个二维的有逆流换 热的 U 形通道反应器 , 计算中采用均匀网络 ( 见图 1) 。反应器长度为 10 mm ,通道宽度为 2 mm ,催化表 面为外径 2 mm 的外圆柱面 , 其上镀上一层 Pt 催化 剂 。催化表面置于燃气和烟气通道的交汇处 , 燃气 在催化表面着火燃烧 。本文联合使用 FLUENT[1 ] 和 DETCHEM[2 ] 数值模拟了甲烷/ 空气混合物在通道宽 度为 2 mm 、 长为 10 mm 、 布置不同催化剂表面微燃烧 器的燃烧过程 。
微通道壁面传热对甲烷催化燃烧影响的数值模拟
微通道壁面传热对甲烷催化燃烧影响的数值模拟近年来，由于石油资源的减少、环境污染的加剧以及全球能源安全的日益受到重视，研究甲烷催化燃烧技术的发展受到越来越多的关注。
由于甲烷催化燃烧的反应机理十分复杂，研究者们对甲烷催化燃烧反应的实验研究仍然处于探索阶段。
将甲烷催化燃烧及其影响机制以数值模拟方式进行探究，可以进一步深入地了解其反应过程及详细的物理机制，从而更有效地设计出更先进的催化剂和工艺条件，结合试验数据，系统地实现对燃烧机理的研究。
其中，微通道壁面传热对甲烷催化燃烧的影响是相当重要的，它影响着反应的动力学，可以改变反应的温度及燃烧比例，控制燃烧方式及产物组成，因此就有必要以数值模拟的方式来进行研究。
本文将以《微通道壁面传热对甲烷催化燃烧影响的数值模拟》为标题，详细地阐述甲烷催化燃烧微通道壁面传热的影响，并且采用数值模拟方法进行研究，为今后甲烷催化燃烧技术的发展提供理论指导。
首先，本文介绍了微通道燃烧的一般原理及基本数学方程，其中包含了热传导方程、质量守恒方程以及动量守恒方程。
然后，详细地介绍了微通道壁面传热对甲烷催化燃烧影响的数学模型，包括壁面传热系数的求解，燃烧速率与温度的关系，及壁面传热对反应活化能等参数的影响。
此外，文章还将数值模拟方法应用于甲烷催化燃烧微通道壁面传热影响的研究，并以一个简化的催化燃烧实验数据为示例，模拟了反应过程的温度分布，模拟结果表明，壁面传热对甲烷催化燃烧反应的影响极其重要，可以有效地提高燃烧速率及改善温度分布。
最后，文章以数值模拟的方法深入研究了甲烷催化燃烧微通道壁面传热对燃烧机制的影响，并给出了详细的模拟结论。
本文研究提示，微通道壁面传热对甲烷催化燃烧有重要影响，应当在设计及优化催化剂和反应条件时得到充分的考虑，以期获得良好的燃烧性能。
因此，本文的数值模拟研究结果可以为甲烷催化燃烧技术的发展提供理论指导，以及为今后研究者研究甲烷催化燃烧机理提供有价值的参考。
本文最后总结了围绕微通道壁面传热对甲烷催化燃烧机理影响的数值模拟研究，指出了微通道壁面传热对甲烷催化燃烧性能的影响，以及此类模拟研究将会为甲烷催化燃烧技术的发展提供有效的科学依据。
甲烷气体的一维冲击反应动力学数值模拟
甲烷气体的一维冲击反应动力学数值模拟甲烷是一种常见的天然气体，具有广泛的应用。
然而，甲烷的燃烧过程中会产生大量的热量和压力波，因此对于其反应动力学的研究具有重要的意义。
本文将介绍甲烷气体的一维冲击反应动力学数值模拟的研究成果。
首先，需要说明的是，一维冲击反应动力学数值模拟是一种基于计算机数值方法的研究手段。
该方法可以模拟气体在一维空间内的流动和化学反应过程，从而研究气体的反应动力学特性。
在甲烷气体的一维冲击反应动力学数值模拟中，研究者首先需要确定气体的初态。
在本研究中，气体的初态为压力为1 atm，温度为300 K的甲烷气体。
然后，研究者需要在气体中引入一个冲击波，以模拟气体在外界作用下的反应过程。
在模拟中，研究者使用了化学动力学模型来描述甲烷气体的反应过程。
该模型包含了甲烷气体的反应机理和速率常数等参数。
通过对模型的求解，可以得到气体中各种反应物和产物的浓度随时间的变化情况。
同时，还可以计算出气体的温度、压力和速度等物理量随时间的变化情况。
为了验证模拟结果的准确性，研究者还进行了大量的实验对比。
实验结果表明，模拟结果与实验结果具有较高的一致性，证明了数值模拟方法的有效性。
通过对甲烷气体的一维冲击反应动力学数值模拟，研究者发现，甲烷气体的反应过程具有明显的非平衡特性。
在冲击波作用下，气体中的甲烷、氧气和氮气等组分会发生复杂的反应过程，产生大量的热量和压力波。
同时，气体的速度和密度等物理量也会随着反应过程的进行而发生变化。
总之，甲烷气体的一维冲击反应动力学数值模拟是一种有效的研究手段，可以帮助我们更好地理解甲烷气体的反应动力学特性。
未来，我们还可以将该方法应用于其他气体的研究中，以深入探究气体反应动力学的本质。
甲烷燃烧器流场数值模拟
zቤተ መጻሕፍቲ ባይዱ
xx 2 yy 2 zz 2
u div(u ) ， x v div(v) ， y w div( w) ， z u v )， y x w u ) x z w v ) y z
1 流体数学模型
1.1 连续方程
质量守恒方程也称为连续性方程。根据 连续性假设， 粘性可压流体其表达式可以表 达为： ( u ) ( v ) ( w ) 0 （1.1） t x y z 对于不可压缩流体， 密度为常数时可以将上 式化简为： u v w 0 x y z （1.2） 速度矢量 V= { u ,v, w } ，式 1.2 为粘性 不可压流体运动的连续方程。
前[1]。李兵臣，宋景慧等对二次风旋流场进 行了数值模拟， 发现不同的调节范围对流场 [2] 的影响较大 。 苏亚欣,汪文辉,邓文义对燃烧 的流场进行了数值模拟， 研究了流场温度分 布情况与燃烧产物的燃烧情况[3]。 张力,王炯 等对低热值煤层气燃烧器结构进行了数值 模拟与优化， 研究了甲烷体积分数为 30%时 的不同比例的直流风和旋流风对燃烧过程 的影响[4]。 冯明杰,李德立等对火焰长度可调 的燃烧器的燃烧过程建立了数学模型， 采用 fluent 软件对流场的热负荷及中心燃料总流 量的比值对燃烧室温度场进行了数值研究 , 结果表明改变燃气流量比， 能够实现火焰长 [5] 度的调整 。钟北京 , 洪泽恺联合 CFD 和
甲烷燃烧器流场数值模拟
李丽，白锦川
（1.酒泉职业技术学院 化学工程系 ， 甘肃 酒泉 735000
摘要：本文通过采用 Workbench Design Molder 对燃烧器的流体进行物理建模，用 Fluent 软 件对空气流速在 10m/s、 20m/s、 30m/s 和 40m/s 情况下的甲烷燃烧情况进行了数值模拟研究。 分析结果表明当甲烷流速一定时随着空气流速的增大甲烷的燃烧程度越来越充分， 供给的热 量越来越多。另外，软件计算结果与实际的产热效果很接近，表明了 Fluent 对燃烧场的模 拟研究可以应用实际工程， 根据本文的计算思路和方法对燃烧的结构优化优化和降低能源消 耗提供一定参考作用。 关键词：甲烷；燃烧器、Fluent、数值模拟
管道内甲烷_空气预混爆炸燃烧的数值模拟fluent
应的输运方程为 :
99ρt +
9 (ρkvx) 9x
=
9 9x
[
(μ
+
μi σk
)
9k 9r
]
+
Gk
+
Gb
-
ρε-
YM
+ Sk
(4)
3 项目资助 :山西省留学人员科研资助项目 。
9 (ρε) 9t
+
9 (ρεvr) 9r
=
9 9r
[
(μ+
μt σε
)
99εx ] +
ε C1ε k ( Gk +
C3εGb)
(6)
初始时刻 t ( t0) = 0 ;初始速度 v ( t0) = 0 。
已燃区 : T = 2 000 K, P0 = 101 325 Pa , YCH4 = 0 , YO2 = 0 , YCO2 = 0. 151 4 , YH2O = 0. 123 9 。
未燃区 : T = 300 K, P0 = 101 325 Pa , YCH4 = 0. 055 , YO2 = 0. 22 , YCO2 = 0 , YH2O = 0 。
损失 ,因 此 假 设 壁 面 为 绝 热 。v ( R , t )
= 0,
9ρ( R , t) 9x
= 0,
9T( R , t) 9x= 0,Biblioteka 9Yi ( R , t) 9x
= 0;绝热壁面热流密度为零 ,即
qw
= 0 ,壁面无滑移 。
管道两端边界条件 :
x =0和
x = 800
mm ,0 ≤t ≤∞处 , v (0 , t)
甲烷超声速燃烧过程的数值模拟
（ 科学院力学研究所．北寐 ｌ （ｏ 中 ０｝） ｏ８
摘
要 ：提 出 ｒ 种 ｃ ／ 证 超 声 速 流 － 燃 烧 的 ６ 程 化 学 反 血 模 、并 且 用 在 超 声 速 高 温 空 － 横 向 喷 Ｏ Ｃ ｐ 方 Ｃ Ｊ ｛
ｌ 引 言
研究 吸 热 碳 氧 燃 料 在 超 声 速 气 流 中 燃 烧 的 机
理 是 ｊ｛ 究 的 热 ｆ 课 题 Ｉ ｆ 订 ｊ 液 态 碳 氯 燃 料 燃 烧 时 的 化学 反 应 和 热 物 理 过 很 复 杂 ， 如 ， Ｈ 例 Ｃ 和 Ｏ 的 完 化 学 反 应 包 括 ２ ０多 个 组 元 、０ ２０多 个 化 学 反 应
ｇｈ 一 ａ ｌ ｎ ｆ ｔ ｅ ｎｅ］ａｎ ｏ ｇ ｓ ｓ ｒ （ｔｏ ｏ ｈ ｔ １ｌ ｅ— ｘｙ ｅｎ ｙ ｋｍ ｔ
ｍ ｄ ｉ ］ ｈ ｄ ｌｏ ｓｎ ａ ｅｗ ｆ ｔｅ ｍｏ ｅ ｆ ｉ． １
Ｋｅ ｒ ｓ Ｓ ｐ ｒｏ ｉ ｏ ｕ ｔ ｎ ｒ， ｔｅ ｇｔ ｙ ｗｏ ｄ ： ｕ ｅｓ ｎｃｃ ｍｂ ｓ ｏ ａ ｎ ｉ ｉ ￣ｅ ｍ；Ｓ ｐ ｒｏ ｉ ｃ ｎ ｂ ｓ ｏ ｕ ｅｓ ｎｃ ｏ ｖ ｕ ｔ ｎ；Ｍｅ ／ｎ ｉ ｌ ｆａ ｅ；Ｎ ｍｅｉａ ｉ ｌ ｉｎ ｕ ｒｃ ｌ ｍｕａ ｏ ｓ ｔ
通过 直接 求 解 雷 诺 方 程 及 补 充 完 全化 学 反 应 的模 型 ，
已 超 ｍ ｉ 的 计 算 能 力 ， ｒｒ前 因此 需 采 用 简 化 的 动 力 学 搂 型 找 到 种 既 能 基 率 反 映 燃 烧 过 程 的 化 学 动 学 特 扯 叉 能 为 现 有 的计 算 能力 所 接 受 的化 学 反 应 模 足 解决 问题 的关 键 。 文 献 ： ］ 完 全 反 应模 型 的 基 ２在 础上， 建 了甲烷 超声 速燃 烧 ４组 元 二 个 基元 反应 的 － 简 化化 学 动 力 学 模 喇 ； 献 ： ］ 用 ｒ乙 烯 燃 料 ５组 文 ３采
基于fluent的甲烷燃烧二维模拟分析
基于fluent 的甲烷燃烧二维模拟分析陈飞 1434422(同济大学汽车学院，上海)摘要: 目前，由于环境污染和排放法规的日趋严格，降低排放已经成为了汽车工业的重点，而寻求清洁的替代燃料是一种较为有效的解决办法。
甲烷属于可再生气体燃料，可以实现与空气的良好预混，利用fluent 进行甲烷燃烧的数值模拟进行仿真分析。
Fluent 提供了涡耗散模型用于求解燃料快速燃烧，整体反应速率由湍流控制的单步或双步总包燃烧反应。
根据模拟结果分析甲烷作为车用替代燃料的可行性。
关键词: 替代燃料；燃烧的数值模拟；甲烷燃烧；fluent 仿真1. 引言燃烧是燃料与氧化剂之间的发光发热的化学反应，根据反应前各组分的分布，可以分为预混燃烧，扩散燃烧和部分预混燃烧。
其中预混燃烧较多的应用于汽车工业的车用汽油发动机。
目前，由于环境污染和排放法规的日趋严格，降低排放已经成为了汽车工业的重点，而寻求清洁的替代燃料是一种较为有效的解决办法。
1.1. 燃烧的数值模拟燃烧的数值模拟是通过CFD 软件实现对实际燃烧过程的仿真模拟，求解流畅流动特性及其混合特性，温度场、组分浓度场以及颗粒和污染物排放等，从而提供实际燃烧过程的参考，对于产品研发，科学研究都有很大的意义。
燃烧的数值模型主要运用模拟软件根据燃烧模型进行仿真，目前可用于燃烧数值模拟的软件有FLUENT,STAR-CD,CHEMKIN,KIVA 等。
燃烧模型主要根据不同燃烧的特点设置求解参数，包括如下内容：稀疏相模型、输运控制方程、燃烧模型、辐射换热模型、污染物模型。
Fluent 提供了涡耗散模型用于求解燃料快速燃烧，整体反应速率由湍流控制的单步或双步总包燃烧反应。
其中对于反应r 中的物质i 的产生速率r i R ,由下面两个式子给出：⎪⎪⎭⎫ ⎝⎛=R w r R R R i w r i r i M v Y k A M v R ,',,',,min ερ （1.1） ∑∑=N j jw nr j p p i w r i r i M v Y k AB M v R ,,,',,ερ （1.2） 式中，p Y ——任何一种产物的质量组分；R Y ——某种产物的质量组分；A——经验常数4.0；B——经验常数0.5。
Fluent大作业
Fluent大作业——圆筒燃烧器内甲烷燃烧的数值模拟引言：根据公安部消防局的统计数据，2010年因火灾死亡的人数为1205人，其中多数人是因为火灾产生的有毒有害高温气体而死，因此研究火灾中有毒有害气体的分布有着重要意义。
下面以一个简单的模型，对一个圆筒燃烧器内的甲烷和空气的混合物的流动与燃烧过程进行研究，模拟其中的温度场、有害气体的分布情况。
问题描述：长为2m、直径为0.45m的圆筒燃烧器结构如下图所示，燃烧器壁上嵌有三块厚为0.005m，高0.05m的薄板，以利于甲烷与空气的混合。
燃烧火焰为湍流扩散火焰。
在燃烧器中心有一个直径为0.01m、长0.01m、壁厚为0.002m的小喷嘴，甲烷以60m/s的速度从小喷嘴注入燃烧器。
空气从喷嘴周围以0.5m/s的速度进入燃烧器。
总当量比约为0.76（甲烷含量超过空气约28%），甲烷气体在燃烧器中高速流动，并与低速流动的空气混合，基于甲烷喷口直径的雷诺数约为5.7X103。
图1燃烧器结构使用通用的finite-rate化学模型分析甲烷-空气混合与燃烧过程。
同时假定燃料完全燃烧并转换为CO2和H2O。
反应方程为CH4+2O2→CO2+2H2O反应过程是通过化学计量系数、形成焓和控制化学反应率的相应参数来定义的。
计算结果：图2采用恒定的Cp值（1000J/kg·K）计算的温度分布图3采用mixing-law计算的温度分布从上面两图可以看出，当Cp值恒定为1000J/kg·K时，最高温度超过2900K。
火焰温度的计算结果偏高，可以通过一个更真实的依赖于温度和组分热容模型来修正。
比热对温度和组分的依赖性将对火焰温度的计算结果有着明显的影响。
Mixing-law会得到基于全部组分质量分数加权平均的混合比热。
在Fluent中，还有一个Fluent物性数据库随温度变化的Cp（T）多项式，可以启动组分比热随温度的变化特性。
设置后的计算结果如图2，可以看出最高温度已经降低到大约2200K。
甲烷旋流预混燃烧的数值模拟研究
中国工程热物理学会 燃烧学 学术会议论文 编号：084177甲烷旋流预混燃烧的数值模拟研究杨朝乐，林伟荣，原 鲲*，于溯源清华大学核能与新能源技术研究院，北京 100084(*摘要：本文采用数值模拟手段，对甲烷旋流射流预混火焰进行了研究，重点考察了射流出口旋流数、反应当量比和射流出口速度对火焰形貌，以及NOx 和CO 排放的影响。
研究结果包括：NOx 和CO 的排放随着当量比的降低而下降；旋流可使CO 排放大幅下降，但主要发生在弱旋工况；对于弱旋流火焰，预混射流速度的提高增加了CO 排放，而在中强旋工况下，情况正好相反。
本研究结果为以后的实验设计和开展提供了理论指导。
关键词：预混，旋流，燃烧0 前言随着我国社会经济的快速发展，对包括天然气在内的清洁能源的需求正逐年上涨。
目前我国天然气占一次能源消费比重不到3.5%，尚不足国际水平（25%左右）的1/7。
但年度消费量增长迅速，近年来每年增幅在10%以上。
天然气已经成为一些城市客车的重要交通燃料，也成为北方许多城市冬季采暖的首选燃料。
与燃煤相比，天然气没有二氧化硫和粉尘排放，被称为清洁燃料，但是在氮氧化物（NOx ）减排方面效果有限。
天然气主要成分为甲烷，其燃烧产生的NOx 污染正逐渐引起人们关注。
天然气燃烧产生氮氧化物的机理主要是热力型氮和快速型氮，其中快速型氮的生成过程相对复杂，但热力型氮是最主要的生成产物，而且排放水平和火焰温度高度相关。
采用贫燃料预混火焰，可以大幅度降低火焰温度，从而将NOx 排放控制在很低的水平[1]。
预混火焰越接近贫燃极限，NOx 排放越低。
但是在接近贫燃极限的工况下，如何使预混火焰得到稳定成为问题的瓶颈。
作为一种稳定火焰的方式，旋流燃烧的火焰稳定能力已得到实践检验[2]，传统上还主要用来稳定扩散火焰。
旋流燃烧的主要功能包括：1）在火焰中心形成回流区以实现火焰稳定；2）强化燃气混合以提高燃烧效率；3）缩短火焰长度以节省燃烧室尺寸。
微细直管内甲烷燃烧的数值模拟研究
ｅ ｔｍ ａ ｅ ｔｌｓｅ ｆｅ ｔｄｔｅａ ｉｌ ｅ ｔ ｒ ｓｅ ｓｄ ｅｗａｌ ａ ｉｌ ａ ａ ｆｒｉ ｅｆ ｉ，ｆ ｅ ｇ ｉｏ ｄｓａ ｉ ｔ ａ ｅ ｕｔｅ， ｘｅ ｌ ａ ｓ ｓａ ｃｅ ｘａ ａ ａ ｆｒｉ ｉｅｔ ｌ ｈ ｏ ｈ ｈ ｔｎ ｎ ｈ ，ｒｄａ ｔｒ ｅ ｔ ｕｄ ｕ ｌ ｎｔ ｎａ ｔｂｌ ｆ ｅ ｂ ｔｎ ｎｈ ｌ ｉ ｉ ｎ ｉ ｌ ｙ ｍ ，ｆｒｈ ｒ
ｂ ｓｉ ｎ ｃ ａ ａ ｔ ｒｓｉｓｗｅ ｅｓｕ ｉ ｄ ｈ ｓ ｌ ｆｓ ｕ ｔ ｈ ｒ ｃ ｅ ｉｔ ｒ ｔ ｄ ｅ ．Ｔ ｅｒ ｕ ｔ ｏ ｍｕ ａｉ ｎ ｉ ｄ ｃ ｔ ｄ ｔ ａ ｏ ｕ ｔ ｒ ｉ ｎ ｉ ｎ ，ｈ ａ ｏ ｄ ｃ ｉ ｉ ｎ ｉｋ ｅ ｓｏ ｗａ ｌ ｏ ｃ ｅ ｓ ｉ ｌ ｔ ｉａ ｅ ｔ ｍｂ ｓ ｍｅ ｓｏ ｓ ｅ ｔ ｎ ｕ ｔ ｔ ａ ｄ ｔ ｃ ｎ ｓ ｆ ｌ ｏ ｎ ｈ ｃ ｏ ｄ ｃ ｖｙ ｈ ，
Ａｂ ｔ ａ ｔ ｓｒ ｃ ：Ａ ｏ ｕ ａ ｉ ｎ ｌ ｕ ｄｄ ｎ ｍｉ ｓ ｓ ｍｐ ｏ ｅ ａ ｃ ｌｔ ｅｃ ｍｂ ｓｉ ｎ ｏ ｐ ｅ ｘ ｄｍｅ ａ ｅ ａ ｒ ｎｔ ｅ ｗｏ ｄ ｍｅ ｓｏ ａ ｃ ｏ ｃ ｍｐ ｔ ｔ ａ ｉ ｙ ａ ｃ ｏ ｌ ｆ ｗａ ｅ ｌ ｙ ｄｔ ｃ ｌ ｕ ａｅｔ ｏ ｏ ｈ ｕ ｔ ｆ ｒ ｍｉ ｅ ｔ ｎ － ｉ ｉ ． ｉ ｎ ｉｎ ｌ ｏ ｈ ｈ ｔ ｍｉ ｒ ．
ｉｄ ｃ ａ ｘｉｃｉｎ ｎ ｕ ｅｆ ｍｅｅ ｔ ｔ ．Ｗ ａ ｕ ｈ ｅ ｓｅ ｈ ｎ ｅ ｕ ｄ ｄ ｓ ｒ ｅ ｎ ｅ ｔ ｘ ｈ ｇ ｅ ｅ ｎ ｆ ｉ ｄｗａ１ ｎ ｕ ｎ ｅ ｏ ｕ ｔ ｎ ｓ ｂ ｌ ｙ ｌ ｎ ｏ Ｕｒ ｇ ｎ ｓ ｎ ａ ｃ ｄ ｆ ｉ ｉｏ ｄ ｒａ ｄ ｈ ａ ｃ ａ ｅｂ ｔ ｅ ｕ ｄ ａ ｌ ｏ ｌ ｅ ｎ ｗ ｌ ｎ ．ｉ ｆ ｅ ｃ ｄ ｃ ｍｂ ｓｉ ｔ ｉ ｔ ． ｌ ｏ ａ ｉ
各向同性多孔介质燃烧器的低溶度甲烷燃烧二维数值模拟
各向同性多孔介质燃烧器的低溶度甲烷燃烧二维数值模拟一、模拟对象描述由于圆柱状燃烧器的对称性，可简化为二维模拟。
如图1.1所示，在进口处有一个小喷嘴，甲烷体积溶度为20%的甲烷-氮气混合气体以80m/s 的速度从喷嘴成圆锥状喷出，入口管道内径为0.01m 。
周围空气以0.5m/s 的速度流入燃烧器，通风管道内径为0.076m 。
过量空气系数为 1.28，甲烷混合气流动的雷若数为5700，入口燃气、入口空气和绝热管壁温度皆为300k ，甲烷与空气的反应采用最常见的单步总包反应，如式1.1所示。
422222CH O CO H O +→+ （1.10）图1.1 二维多孔介质燃烧器二、数学模型。
1、多孔介质的基本参数1）代表性单元体在多孔介质的一定范围内选取一个控制体，用这个控制体的局部容积求得的任何参数的平均值具有代表性，同时，从形态上看，这个控制体包含了足够多的多孔介质微结构体。
这个控制体称为代表性单元体（REV），也叫表征体元，图2.10就是一个REV。
图2.1就是一个REV。
REV是人为选定的控制体，它强调的是“代表性”，即该单元体必须在固体骨架的形态及结构特性、孔隙尺度特性以及流体在孔隙中的流动和传热特性等方面能基本代表整个计算区域，其相关基本参数的平均值能逼近真值。
REV选定的要点是：（1）REV的尺度远小于整个计算区域的宏观尺度。
（2）REV应包含足够多的固体骨架和孔隙，以致单元体内的孔隙体积远大于单个孔隙空间。
（3）在计算区域内，REV所表征的基本参数随空间坐标的变化幅度很小。
2）孔隙率孔隙率是指在多孔介质的REV中固体骨架间孔隙的总体积V与该单元体多f孔介质总体积V之比，即图2.1 多孔介质代表性单元体结构示意图100%fV V φ=⨯ （2.10）当多孔介质中被流体充满（单相或多相）时，f V 也就等于流体所占据的空间体积。
孔隙率φ是多孔介质研究和应用中非常重要的一个参数。
3）比面比面是指多孔介质固体骨架的总面积s A 与多孔介质总容积V 之比，即 s f A a V= （2.11） 在相同容积的REV 中，f a 正比于s A ，而对于细小骨架的多孔介质，其骨架表面积要大于粗骨架的多孔介质，因此，骨架或颗粒越细，所对应的比面f a 也就越大。
微通道内甲烷催化燃烧的数值模拟研究的开题报告
微通道内甲烷催化燃烧的数值模拟研究的开题报告
一、研究背景和意义：
微通道内的催化燃烧技术作为一种高效低污染的能源转化方式，在化工、航天、汽车等领域得到了广泛的应用。
与传统的反应器相比，微通道反应器内的催化燃烧反应速率更快，反应效率更高，且所需催化剂量更少，经济性更好。
本文将研究微通道内甲烷催化燃烧的数值模拟研究。
该研究的意义在于，通过数值模拟研究，可以更全面地了解微通道内甲烷催化燃烧的反应过程和应对方法，为微通道催化燃烧技术的发展提供理论依据。
二、研究内容：
1. 对微通道内甲烷催化燃烧反应机理进行建模，分析微通道内流场和组分分布的特点，探究反应过程中的化学反应和传质过程。
2. 基于CFD方法，对微通道内甲烷催化燃烧反应的数值模拟进行研究。
建立微通道内甲烷催化燃烧的反应动力学模型，采用Fluent等数值模拟软件进行模拟分析。
3. 通过模拟分析，探究微通道内甲烷催化燃烧的影响因素，如反应温度、流速、催化剂种类和浓度等因素对反应效率的影响，并根据模拟结果提出相应的优化方案和建议。
三、研究方法和流程：
1. 文献调研和资料收集，建立微通道内甲烷催化燃烧反应机理的理论模型。
2. 建立CFD数值模拟模型，对微通道内甲烷催化燃烧反应进行数值模拟分析，并进行模拟结果的验证。
3. 在模拟分析的基础上，设计相应的实验并进行实验验证。
4. 分析模拟结果和实验数据，提出优化方案和建议。
四、预期结果和意义：
通过数值模拟的研究，可以更全面地了解微通道内甲烷催化燃烧的反应机理和流场特点，提高微通道催化燃烧技术的效率和可行性，对于推进微通道燃烧技术的应用和发展具有重要的理论和应用价值。
基于化学动力学的甲烷燃烧性能数值模拟
法等 ； 文献 ［ ８ ］ 研 究 了 甲烷 在 不 同 水 油 混 合 比下 的燃烧 特性 ； 文献 ［ ９ ］ 致 力 于 掺 氢 来 改 善 天 然 气
［ ５ ］ 利用 三维 Ｃ Ｆ Ｄ软件 对 甲烷 的燃烧 和排 放进 行
了研 究 ， 得 出了各排 放物 在缸 内的分布 隋况 ； 文献 ［ ６ ］ 针 对 甲烷燃 烧 过 程 中 的碳 氢 排 放 ， 提 出 了一
种快 速有 效 的 测 试 碳 氢 排 放 量 的 方 法 ； 文献 ［ ７ ］
量密度 和抗爆 性 ， 还 具有 节 能减 排 的特 点 ， 因此 在
发 动机 的运行 特 性 ， 研 究 了不 同掺 氢 比例 下 甲烷
燃 烧 的变化 状况 。 笔 者在 甲烷 燃 烧 的化 学 动力 学 过 程 的 基 础
替代 能源 中受 到 重视 。但天 然气 在 均质 压 燃 的过 程 中会 出现相应 的缺点和不 足 ， 在此背 景下 笔者对 天 然气 中的主要 成分 甲烷 的燃烧性 能进行 了研究 。 天然气 的 主要 成 分 为 甲烷 。其 主 要 优 点 有 ：
从 表 ｌ可 以看 出 ， 甲烷和 正庚 烷在 引燃 温度 、
分子量 和 热值方 面完 全 不 同 ， 属 于两 种 完 全 不 同 的燃料 。 正庚 烷 具 有 与 柴 油 相 当 的 十 六 烷 值 （ 均 为
献［ ４ ］ 分别 以敏感度分析 和最优简化 法为基础 ，
对Ｇ Ｒ Ｉ —ＭＥ Ｃ Ｈ ３ ． ０机理进 行 简化 ， 并 与 原 有机 理 进 行 了对 比 ， 得 到 了简化后 的甲烷反 应机 理 ； 文献
燃 烧 产 生 的污 染 物少 、 使 用 特 性好 、 抗 爆性 好 、 与 空 气混 合气 着火 界 限宽 、 经济性 好 。
甲烷加湿燃烧的数值模拟
第３８卷 第５期２０２０年９月 石化技术与应用ＰｅｔｒｏｃｈｅｍｉｃａｌＴｅｃｈｎｏｌｏｇｙ＆Ａｐｐｌｉｃａｔｉｏｎ Ｖｏｌ．３８ Ｎｏ．５ Ｓｅｐ．２０２０研究与开发（２９５～２９９）甲烷加湿燃烧的数值模拟王珂（长江大学石油工程学院，湖北武汉４３００１０）摘要：对空气氛围下甲烷的加湿燃烧进行了数值模拟，在保持甲烷与空气流量恒定的条件下，考察了加湿率对燃料燃烧特性、燃烧速率以及主要污染物（碳烟、氮氧化物、一氧化碳）排放的影响。
模拟结果表明：随着加湿率的增加，燃料的燃烧速率提高，各污染物的排放均减少，一氧化碳和碳烟排放分别在加湿率为０～２０％，０～１０％时下降明显，氮氧化物排放则呈线性降低；综合考虑节能环保与燃烧器的安全运行，燃料的最佳加湿率不宜超过２０％。
关键词：甲烷；加湿燃烧；空气；加湿率；碳烟；氮氧化物；一氧化碳；数值模拟中图分类号：ＴＱ２２１．１＋１ 文献标志码：Ｂ 文章编号：１００９－００４５（２０２０）０５－０２９５－０５ 加湿燃烧是在燃料燃烧时向其中加入水或水蒸气来提高燃烧效率、减少污染物排放的技术。
根据加湿方式可分为燃料加湿燃烧和氧化剂加湿燃烧，根据氧化剂不同可分为空气加湿燃烧和富氧加湿燃烧。
目前，关于加湿燃烧的研究主要集中在氧化剂加湿燃烧过程中对污染物排放的分析［１－７］，加湿燃烧对抑制碳烟和氮氧化物（ＮＯｘ）生成的效果已成为行业共识，但关于加湿燃烧对ＣＯ生成影响的研究还不多见。
且研究的氧化剂环境多为富氧条件或如Ｈ２Ｏ／Ｏ２氛围的非常规氧化剂［７－８］，工业应用成本较高。
因此，研究常规空气氧化剂氛围下的加湿燃烧对节约燃烧成本、降低技术难度方面具有工业应用价值。
本工作以小型圆筒燃烧器为研究对象，对常规空气氛围下甲烷的加湿燃烧进行了数值模拟，考察了常规空气氛围下甲烷的加湿燃烧性质及加湿率对主要污染物排放的影响，确定了最佳加湿率，可为工业锅炉的环保、高效、安全的运行提供借鉴。
１ 模型建立与数值模拟方法选择１．１ 模型建立与网格划分模型采用圆柱燃烧器，建立简化后的几何模型，其外形与尺寸如图１所示。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
甲烷燃烧的数值模拟及分析
主要分为三个部分，第一部分讲解甲烷在炉膛内燃烧的模型建立的方法；第二部分对甲烷燃烧模型的数值模拟结果进行分析和比较；第三部分为结论。
一、模型建立
1、在Gambit中建立计算区域
在本例中建立圆柱形炉膛，并研究甲烷和空气在炉膛内的燃烧反应。
物理模型如下：
甲烷入口直径为10mm；空气入口直径为50mm；炉膛为直径为500mm；长度为1200mm的圆柱形。
如图1。
图1圆柱形炉膛模型图
2、绘制网格
图2进口网格分布
甲烷和空气进口的网格元素选择四边形网格，如图2。
炉膛表面的网格也是四边形网格，如图3。
图3炉膛表面网格分布
图4炉膛表面网格分布
图5炉膛出口网格分布
图6炉膛内部网格分布
3、指定边界条件
图7炉膛边界条件
Inlet1为甲烷入口，边界条件为速度入口；
Inlet2为空去入口，边界条件为速度入口；
Outlet为炉膛出口，边界条件为自由流；
其他炉膛壁面为墙体，边界条件为墙体。
4、导入fluent
具体信息如下：
54440mixed cells,zone2,binary.
326quadrilateral wall faces,zone3,binary.
1900quadrilateral wall faces,zone4,binary.
350quadrilateral wall faces,zone5,binary.
218quadrilateral outflow faces,zone6,binary.
204quadrilateral velocity-inlet faces,zone7,binary.
18quadrilateral velocity-inlet faces,zone8,binary. 108880triangular interior faces,zone10,binary.
11144nodes,binary.
11144node flags,binary.
缩放信息如下图：
图8缩放信息图
5、选择计算模型
图9定义求解器
图10考虑能量方程
图11考虑粘性模型
图12考虑辐射模型
图12考虑燃料模型
图13燃烧物质和炉膛材料
6、操作环境的设置
图14操作环境（压力场和重力场）
7、设置边界条件
图15空气入口边界条件
空气入口的速度为8m/s，温度为300K，入口空气中氧气的含量为21%。
图16甲烷入口边界条件
甲烷入口的速度为19.1m/s，温度为300K。
图17炉膛出口边界条件
图18炉膛壁面边界条件
炉膛壁面的边界条件设置为steel钢。
一般fluent软件中避免材料默认设置为铝，考虑到铝的熔点温度为660℃，甲烷在炉膛中燃烧的温度可能很高，炉膛存在融化的危险，为了让模拟条件尽量接近真实的条件，所以选用钢材料，普通钢的熔点为1500℃。
8、求解方法设置及控制
求解参数设置压力速度耦合求解方式选择“SIMPLE”，Discretization对应的压力离散方式选择“Standard”，其他离散方式选择“QUICK”。
QUICK格式一般用于四边形/六面体时具有三阶精度，用于杂交网格或三角形/四面体时只具有二阶精度，比second-order upwind有更高的精度。
图19求解方法设置
将入口的甲烷温度和空气温度设为初始值为300K，开始迭代，迭代1000步。
9、模拟结果
可以看到模拟结果较符合甲烷在炉膛内燃烧时的温度分布，火焰成本生灯状。
并在下面列出速度分布图和各个成分的摩尔质量图。
图20温度分布图
图21速度分布图
图22甲烷摩尔质量
图23氧气摩尔质量
图24氮气摩尔质量
图25水的摩尔质量
图26二氧化碳的摩尔质量
二、分析与对比
1、改变计算模型
将原来的稳态模型，改为非稳态模型，非稳态方程选择默认的1st-Order Implicit。
可以从下图中看到，整个炉膛内的温度呈均匀分布，没有燃烧反应进
行。
由于没有燃烧反应进行，炉膛中的水和二氧化碳都没有生成。
图27计算模型选项图
图28温度分布图
图29速度分布图
图30甲烷摩尔质量
图31氧气摩尔质量
图32氮气摩尔质量
图33二氧化碳摩尔质量
2、改变求解方式
保持求解参数设置压力速度耦合求解方式选择“SIMPLE”，Discretization 对应的压力离散方式选择“Standard”。
而将其他离散方式改为“First Order Upwind”一阶迎风。
已知Quick格式一般用于四边形/六面体时具有三阶精度，用于杂交网格或三角形/四面体时只具有二阶精度，比second-order upwind有更高的精度。
First Order Upwind格式具有较好的收敛性但是精度低，在模型简单的结构网格计算上，first-order解算精度与second-order区别不大。
所以本模型尝试使用First Order Upwind来计算。
下图将对比使用First Order Upwind格式计算与使用Quick格式计算的结果，燃烧温度分布和速度分布没有太大的不同，所以主要比较氧气、氮气、水和二氧化碳生成的模拟结果。
可以看到使用First Order Upwind格式计算，燃烧生成物扩散的更大。
图34求解方式选项图
图35氧气摩尔质量对比图（右边为first-order格式计算）
图36氮气摩尔质量对比图（右边为first-order 格式计算）
图37水的摩尔质量对比图（右边为first-order 格式计算）
图38二氧化碳摩尔质量对比图（右边为first-order 格式计算）
First Order Upwind
Mass Flow Rate (kg/s)
----------------------------------------------------
air-inlet
0.018060772outlet -0.01897534
------------------------------------
Net -0.00091456872
三、结论
从上面的模拟结果可以看出，氧气和甲烷在燃烧的前半段基本反应完毕，氮气的生成主要在火焰的中间区域，二氧化碳和水的生成区域相似。
从比较和分析中可以看出，对大多数的简单模型，First Order Upwind格式的计算已经足够了，选用其他的格式计算并不一定能够提高精度和速度。
在Fluent的模拟计算中，选择合适的计算模型和求解方式对最后的模拟结果有很大的影响。
又由于这种影响的随机性喝不确定性，使模拟的可信度降低，只有和实验多做对比，更多的累积经验，才能提高模拟的准确性和可信度。