第六章%20预混层流火焰
chemkin模拟稳态一维层流
稳态一维层流燃烧实验 一实验目的 一维层流火焰在预混燃料-氧化剂混合物中传播是最简单的燃烧现象之一。在此火焰中，化学动力学以及能量和组分扩散输运起重要作用。描述一维预混火焰的方程组是：压力为常数的条件下的质量守衡，能量守恒，组分守恒以及理想气体状态方程。虽然守恒方程和状态方程提供了缓燃的未燃气体和已燃状态之间的关系，但不能唯一确定层流火焰速度，其准确解只有通过数值积分才能获得。本实验考察层流火焰的传播速度以及与燃烧参数如燃料类型、化学配比、压力及未燃气体的温度的关系。 二实验方法 采用Chemkin自带的实例flame_speed_freely_pripagating.ckprj（甲烷-空气火焰传播速度）。模拟绝热、大气压力、自由传播、化学当量混合甲烷-空气预混合燃烧火焰的传播速度。模拟计算中不考虑NO x形成，仅采用甲烷-空气骨架燃烧机理。火焰用详细轴向温度分布做定温计算。设置火焰温度（在入口温度到峰值温度间），通过调节反应器内部的计算区域，来获得预热到反应完整过程，保证初始温度变化曲线足够平坦（温度梯度为0），计算报表反馈火焰传播速度。 三实验步骤 ?启动Ckemkin ?点击Open Project ?双击samples ?单击flame_speed_freely_propagating.ckprj ?单击Select按钮 ?双击左侧浏览器中的Pre-Processing选项 ?在弹出的新窗口中，点击Run Pre-Processon按钮，①View Results...按钮
可用；②左侧浏览器中出现Run Model选项 ?(可选)点击View Results...按钮，可查看甲烷的气相反应机理和气相传递 数据。 ?双击右侧浏览器中的Run Model选项，出现Run Model(flame_speed_freely_propagating)窗口。 ?点击Create Input File按钮，Run Model按钮可用。 ?点击Run Model按钮，计算甲烷-空气层流燃烧。 四查看和分析实验结果 1）查看实验结果 ?打开工作目录下的flame_speed_freely_propagating.out文件，编辑/查找，在 查找对话框内输入“cm/s”，查找该文件中最后一个速度栏。在该栏下的第一个数值就是层流火焰的传播速度，为41.01cm/s. ?点击Run Model窗口中的Run Post Processor按钮，弹出Select Results to Import to Plotting Package窗口。 ?点击Solution Sets选项卡，选择最后一个计算结果： Solution_variables_vs_distance_for_Soln_No_3 ?点击Species/V ariables选项卡，首先点击None按钮，清除所有被选中的选项， 然后选择所需查看的参数，点击OK按钮显示二维曲线图。 ①压力pressure-distance曲线 图1 压力-距离曲线
实验二-层流火焰传播速度的测定实验
实验二-层流火焰传播速度的测定实验
实验二层流火焰传播速度的测定实验 一、预备知识 1、火焰传播和化学反应 燃烧发生了一系列化学反应，在这些反应中，燃料在一些自由基例如O、OH、H碰撞下发生反应，产生更多的H或者是分解成更小的碎片。 例如，CH4被连续地转化成CH3，CH2，CH。最初形成的各种氧化的中间产物与燃料中的碳结合而首先变为CO，并且燃料中的氢基变为H2，所有的中间产物将接着进一步氧化，再一次通过自由基的作用，而变为 CO2和H2O。总热量的一大部分释放都是发生在第二阶段。这个次序使燃烧具有自持性，且只能够发生在高温下（如1500K以上）。因为只有在高温下，才能是自由基产生的速率比消耗的速率快，而这对燃料完全变形 以及中间产物的氧化是有必要的。 当点燃预混燃料时，局部温度将提高到一个非常高的值，提高了反应速率，从而也引起燃料的燃烧，并且释放出热量。通过热传导把热量 引导到了未燃的相邻区域，相邻区域的温度以及反应率都提高了，因此 燃烧就在那里发生了。我们知道，热量的扩散是火焰传播的原因，燃烧 波传播的速度取决于燃烧后的温度以及未燃混合物的热扩散性。为了把 高温区域的自由基传递到与之接触的低温的未燃混合物中，质量扩散也 是很重要的；通常质量和热扩散率是相同的。 在本实验中，未燃混合物的压力和温度与环境大气一致。火焰传播速度只依赖于混合物中的燃料/氧化剂的数量，它们反过来又控制着火焰的温度。贫油(Φ<1)和富油(Φ>1)的火焰温度比化学恰当比(Φ=1)时更低因 为偏离化学恰当比时多余的物质吸收了由可燃燃料燃烧所产生的热量。 实际上，温度最大值出现在当量比比1稍大一些的地方，因为产物的比热容比化学恰当比时稍低。 如果混合物过贫，燃气温度将太低，而不能产生大量的自由基，因此火焰传播变得不可能。如果混合物过富，大量的燃料将吸收自由基， 因此使燃烧第二阶段不能进行。因此，火焰传播只在某个当量比范围内 才有可能，这被称为可燃极限。对于甲烷—空气混合物，其贫燃极限是 Φ=0.53，其富燃极限是Φ=1.6。 2、火焰稳定性
本生灯法测定燃气法向火焰传播速度测试装置
本生灯法测定火焰法向传播速度实验指导书 一、实验目的 1．巩固火焰传播速度的概念，掌握本生灯法测量火焰传播速度的原理和方法。 2．测定液化石油气的层流火焰传播速度。 3．掌握不同的气/燃比对火焰传播速度的影响，测定出不同燃料百分数下火焰传播速度的变化曲线。 二、实验原理 层流火焰传播速度是燃料燃烧的基本参数。测量火焰传播速度的方法很多，本试验装置是用动力法即本生灯法进行测定。 正常法向火焰传播速度定义为在垂直于层流火焰前沿面方向上火焰前沿面相对于未燃混合气的运动速度。在稳定的Bensun 火焰中，内锥面是层流预混火焰前沿面。在此面上某一点处，混合气流的法向分速度与未燃混合气流的运动速度即法向火焰传播速度相平衡，这样才能保持燃烧前沿面在法线方向上的燃烧速度(图1)，即 0sin s u u α=? 式中：u s －混合气的流速（cm/s ）； α－火焰锥角之半。 或 0318 q u = 式中：q v －混合气的体积流量（L/s ）； h －火焰内锥高度（cm ）； r －喷口半径（cm ）。 上式是使用本生灯火焰高度法测定可燃混合气体的层流火焰传播速度0u 的计算式。在我们的实验中，可燃混合气体的体积流量v q 是用湿式流量计分别测定燃气与空气的体积流量而得到的，内锥焰面底部圆的半径r 可取本生灯喷口半径；内焰锥高度h 可由测高尺测量。 三、实验设备结构 实验台由本生灯、旋涡气泵、湿式气体流量计、U 型管压差计、测高尺等组成。旋涡气泵产生的空气通过泻流阀、稳压罐、湿式气体流量计、调压阀后进入本生灯，燃气经减压器、湿式气体流量计、防回火器、调压阀后进入本生灯与空气预混合，点燃后通过测量内焰锥高度计算火焰的传播速度。 四、实验步骤 1、启动旋涡气泵，调节风量使本生灯出口流速约为0.6m/s ，并由湿式流量计读出空气流量。 2、由以上空气流量，可粗略地估算出一次空气系数1α约为0.8、0.9、1.0、1.1、1.2时的燃气流量。 3、开启燃气阀，调整燃气流量分别为上述5个 计算值的近似值（流量值由流量计读出）。
第15章 预混燃烧模拟
第十五章预混燃烧模拟FLUENT有一个预混湍流燃烧模型，基于反应过程参数方法。有关这一模型的内容按以下节次给出： ●15.1 概述和限制 ●15.2 预混燃烧模型 ●15.3 使用预混燃烧模型 15.1 概述和限制 15.1.1 概述 在预混燃烧中，燃料和氧化剂在点火之前进行分子级别的混合。火焰前锋传入未燃烧的反应物产生燃烧。预混燃烧的例子有吸气式内燃机，稀薄燃气轮机的燃烧器，气体泄露爆炸。 预混燃烧比非预混燃烧更难以模拟。原因在于（亚音速）预混燃烧通常做为薄层火焰产生，并被湍流拉伸和扭曲。火焰传播的整体速率受层流火焰速度和湍流涡旋控制。层流火焰速度由物质和热量逆流扩散到反应物并燃烧的速率决定。为得到层流火焰速度，需要确定内部火焰结构以及详细的化学动力学和分子扩散过程。由于实际的层流火焰厚度只有微米量级或更小，求解所需要的开销是不可承受的。 湍流的影响是使传播中的层流火焰层皱折、拉伸，增加了薄层的面积，并因此提高了火焰速度。大的湍流涡使火焰层皱折，而小的湍流涡，如果它们比层流火焰的厚度还小，将会穿过火焰层并改变层流火焰结构。 与之相比，非预混燃烧可以极大地简化为一个混合问题（例如，14.1节中介绍的混合物组分方法）。预混燃烧模拟的要点在于捕获湍流火焰速度，它受层流火焰速度和湍流的影响。 在预混火焰中，燃料和氧化剂在进入燃烧设备之前已经紧密混合。反应在燃烧区发生，这一区域将未燃烧的反应物和燃烧产物隔开。部分预混火焰具有预混和扩散火焰两方面的性质。它们发生在有额外的氧化剂或燃料气流进入预混系统，或是当扩散火焰离开燃烧器以在燃烧前产生某些预混的情况。 预混和部分预混火焰FLUENT的有限速率公式(见13章)模拟。还可以参阅16章了解更多有关FLUENT部分预混燃烧模型方面的信息。如果火焰是完全预混合的，则只有一股具有单一混合比的气流进入燃烧器，可以使用预混燃烧模型。 15.1.2 限制 在使用预混燃烧模型时有以下限制： ●必须使用非耦合求解器。预混燃烧模型在两种耦合求解器中都不能得到。 ●预混燃烧模型只对湍流、亚音速模型有效。这一类型的火焰成为爆燃。在爆炸中， 可燃混合物被冲击波后面的热量点燃，这一类型的燃烧可以使用非耦合和耦合求解 器用有限速率模型模拟。有关限速率模型见13章。 ●预混燃烧模型不能和污染物（如碳烟和NOx）模型一起使用。但完全预混系统可以 用部分预混模型（见16章）模拟。 ●不能用预混燃烧模型模拟反应的离散相粒子。只有惰性粒子可以使用预混燃烧模 型。 15.2 预混燃烧理论 湍流预混燃烧模型基于Zimont等人的工作[275，276，278]，涉及求解一个关于反应过
球形传播火焰法测量层流火焰速度及Markstein长度的方法研究 吴蕴超
球形传播火焰法测量层流火焰速度及 Markstein 长度的方法研究 摘要：球形传播火焰法目前被广泛地用来测量不同燃料的层流火焰速度和Markstein 长度。本文研究了用于球形传播火焰的各种线性与非线性模型。首先，本文基于详细模型推导出了关于拉伸率与曲率的线性与非线性模型；接下来研究了各种模型的准确度及其在球形传播火焰法中的应用情况。研究结果表明：对于层流火焰速度和Markstein 长度，基于曲率的非线性模型而预测的结果最为准确。本文研究的结果对球形传播火焰法测量层流火焰速度及Markstein 长度有一定的指导意义。 关键词：球形传播火焰法；层流火焰速度；Markstein 长度；非线性模型 0 前言 层流火焰速度被定义为一维、绝热、平面火焰相对于未燃预混气体的速度[1]。层流火焰速度是影响燃料燃烧状况和效率的最重要参数之一。同时，层流火焰速度是验证燃料化学反应机理的重要参数[2, 3]。另一个反映燃料燃烧特性的重要参数是Markstein 长度。Markstein 长度表征了由于拉伸率的存在而引起的火焰速度变化，它决定了预混火焰传播的稳定性[2, 4]。同时层流火焰速度与Markstein 长度也都是湍流预混燃烧模型中的重要参数[5]。近五十年来，由于层流火焰速度与Markstein 长度的重要性，各种实验方法被用来测量这两个物理量。球形传播火焰法[6-9]由于有着火焰结构简单以及拉伸率定义准确等特点，目前被广泛地用来测量不同条件下不同燃料的层流火焰速度和Markstein 长度。在国内，西安交通大学黄佐华教授及其团队利用球形传播火焰法系统地研究了多种碳氢燃料的层流燃烧特性[9-15]。 目前关于球形传播火焰法的研究重点之一是如何提高测量精度。例如，点火[16]、非稳定性[17]、非球对称[18]、热辐射[19]、气体压缩[19, 20]等影响因素已经得到了系统的研究。然而，在文献中，不同实验者对同种燃料在同一条件下测得的层流火焰速度和Markstein 长度仍有较大的偏差，特别是Markstein 长度的测量相对偏差甚至可以达到300%[20-22]。造成这些偏差的原因有待被量化地解释。一个重要原因在于对实验数据的处理，尤其是那些在理论模型不成立的范围内的实验数据的使用会导致结果偏差很大。 对于球形传播火焰，在压力变化可被忽略时，可以认为已燃气体处于静止状态。因此，球形火焰面移动的速度就等于相对于已燃气体的层流火焰速度，即S b =dR f (t )/dt 。当拉伸率较小时，火焰速度与拉伸率成线性关系[6-9]： K L S S b b b -=0 (1) 其中S b 0和L b 分别是相对于已燃气体的层流火焰速度和Markstein 长度，K =(2/R f )(dR f /dt ) 是球形火焰的拉伸率。根据式(1)，S b 0和L b 可以通过对实验测得的S b 和K 进行线性拟合而得到。未燃烧气体的平面层流火焰速度S u 0则可以通过质量守恒得到：S u 0=σS b 0，其中 基金项目：国家自然科学基金(编号50976003)与内燃机国家重点实验室开放课题(编号K2010-02)
使用非预混燃烧模型
《数值计算与工程仿真》专刊—FLUENT HELP 算例精选中文版（二）
算例 13
引言
使用非预混燃烧模型
煤粉燃烧的模拟包括气相连续流场的建模和它与煤粒非连续相的作用的建 模。穿过气体的煤粒会挥发燃烧并成为与气相反应的燃料源。反应可以用组份 输运模型（the species transport）或模型（the non-premixed combustion）模拟， 在本指南中你将用非预混燃烧模型模拟简单煤粉燃烧炉中的化学反应。 在本指南中你将学会： 1.怎样用 prePDF 预处理程序为煤粉燃料准备 PDF 表格。 2.怎样为非预混燃烧化学模型定义输入条件。 3.怎样定义煤粒的非连续相。 4.怎样解决包含非连续相煤粒的反应的模拟。 非预混燃烧模型用这样的一种建模方法：用一个或二个守恒量，即混合分 数求解输运方程。多种化学组份，包括基团和中间产物组份可能被包含在对问 题的定义当中，而且它们的浓度将来至于混合分数分布的预测。组份的特性参 数是通过化学数据库获得。湍流化学反应是用 Beta 或者双 delta 概率密度函数 来模拟的。关于非预混燃烧模拟方法的更多细节请参看使用手册。
前提条件
本指南是建立在你已经熟悉 FLUENT 的菜单结构并且已经做完指南 1 的基 础上的。因此在建立过程中的一些步骤和解决过程将被省略。
问题描述
本指南中用的煤燃烧系统为一简单的 10m*1m 的二维管道， 如图 13.1 所示。 因为是对称的，所以只模拟宽度方向上的一半区域。2D 管道的进口分为两股流 动。管道中心附近的高速流速度为 50m/s，宽度为 0.125m。另一股流的速度为 15m/s， 宽度为 0.375m.两股流都为 1500K 的空气。 煤粒在高速流的附近以 0.1kg/s
—151 —
https://www.360docs.net/doc/fe1827030.html,
湍流燃烧模型-PDF
PDF 模型 概率密度函数PDF方法以随机的观点来对待湍流问题，对解决湍流化学反应流的问题具有很强的优势。在湍流燃烧中存在一些非输运量( 如反应速率, 密度, 温度及气相体积分数等) 的湍流封闭问题。尽管这些量没有输运方程, 但它们常常是输运变量的已知函数。平均或者过滤高度非线性的化学反应源项会引起方程的封闭问题。因此，用PDF的方法来解决这些非输运量的湍流封闭问题显然是一个既简单又直接的途径。 PDF方法是一种较为流行的湍流燃烧模型，能够较为精确的模拟任何详细的化学动力学过程, 适用于预混、非预混和部分预混的任何燃烧问题。目前, 确定输运变量脉动概率密度函数的方法有输运方程和简化假定两种, 分别称之为输运方程的PDF和简化的PDF。前者建立输运变量脉动的概率密度输运方程，通过求解该方程来获得输运变量脉动的概率分布。后者假定输运变量脉动的概率密度函数的具体形式, 通过确定其中的一些待定参数来获得输运变量脉动的概率分布。湍流燃烧中, 后者应用最为普遍和广泛。在简化的PDF 中, 输运变量脉动的概率密度函数常常采用双 D 分布、截尾高斯分布和B 函数分布等形式。 PDF在理论上可以精确考虑任意详细的化学反应机理，但是其具体求解时需借助其它的模型和算法，而且计算量相对较大。PDF的方程是由N-S方程推导而来，其中的化学反应源项是封闭的，但压力脉动梯度项以及分子粘性和分子扩散引起的PDF的分子输运项是不封闭的，需要引入模型加以封闭。例如，在速度- 标量-湍流频率PDF中，必须采用小尺度混合模型、随机速度模型和湍流频率模型加以封闭。 模化后的输运方程难以用有限容积、有限差分和有限元等方法来求解，比较可行的一种方法是蒙特卡洛(MonteCarlo)方法,在该方法中输运方程被转化为拉格朗日(Lagrangian)方程,流体由大量遵循Lagrang ian方程的随机粒子的系统来描述, 最后对粒子作统计平均得到流场物理量和各阶统计矩。另有与有限容积法相结合的蒙特卡洛法。 PDF 模型的发展 1969年Lungdren首先推导、计算了速度的联合PDF运输方程，避免了对梯度扩散模型进行模拟，对很简单的流动过程得到了简析解[1]。
层流预混火焰稳定浓度界限测定
层流预混火焰浓度稳定界限的测定 实验成员：徐俊卿 郑仁春 韩超 一、实验的理论基础 现代工业要求尽可能高效地利用燃料，这就要求在稳定燃烧条件下，使燃料与空气以适当的比例混合燃烧。因此对燃料燃烧特性及流体力学和热力学方面需要有详细了解。 预混可燃气燃烧时，如果预混气体的速度在火焰锋面上的法向分量大于火焰传播速度，火焰将向下游移动，最后完全熄灭，称为吹脱或吹熄。反之，如果预混气的法向速度小于火焰传播速度，火焰将逆流向上游移动，进入燃烧器内部，即出现回火现象。 在燃烧过程中，出现回火和熄灭都是不允许的。回火会引起爆炸。熄灭使动力机械停止工作，并向周围扩散有毒气体，有中毒和爆炸的危险。 回火现象只能出现在预混燃烧过程中。在扩散燃烧中，燃料和空气是分别送入燃烧室的，在燃烧器内两者并不接触，因此没有回火现象。熄灭或吹脱现象在预混燃烧和扩散燃 煤气 甲烷 （l/min ） 图1-1 甲烷/空气 煤气/空气 稳定曲线 要维持正常的稳定燃烧，就需要避免出现回火或熄灭现象，因而要求知道燃料的稳定燃烧范围。这一稳定界限与燃料/空气比和环境的温度和压力有关。图1—1是甲烷/空气煤气/空气在常温常压下预混燃烧的火焰稳定范围。图中下部是富油燃烧区，稳定范围宽，不易回火和吹脱。但火眼呈黄色，发光。下面一根曲线为火焰尖顶出现黄色的界限。图中上方为火焰吹脱熄灭界限。当预混气处于该曲线上方时，火焰被吹脱。 左边回火界限呈斜山立形，预混气的状态处于山丘内时出现回火现象。回火界限与预混气的空气气/燃气比有关。当空气/燃气比接近化学计量值时，火焰传播速度快，稳定范围小，容易出现回火。但这时燃烧效率最高。 二、实验原理 本实验用一套预混火焰装置（本生灯）进行实验，既可以定性观察预混火焰的各种现象，又可以定量地测定火焰浓度稳定的界限，火焰传播速度等各种燃烧特性。 ～220V 图1-2实验系统图示意图 图1—2为实验装置系统图，小型压气机供给压缩空气，通过一个稳压减压阀保持气压
第九章 层流扩散火焰 参考答案
第九章 层流扩散火焰 9-1 解: 由连续方程及7.48 ()0)(=??+??x r v x r rv r ρρ g r r v r x v x rv rv x v v x rv rv x v x x r r x x x x x )()()()()()(ρρμρρρρ-+????=??+??+??+??∞ g r r v r x rv x v x r v x v rv rv x v x r x x x r x x )()()()()(ρρμρρρρ-+????=????????+??+??+??∞ 由于)(x f r ≠，const v ==ρμ, 忽略浮力 )()(r v r x r x v v v x v x x r x x ????=??+??υ得证 9-7 证明：由于 ∑=1i Y , Y Pr =1-Y F , )1(F s F Y f Y f -+= 带入Y F s s F f f Y f +-=)1( s s F f f f Y --=1 得证 9-8 解： (1) 222226276.33276.3aN O H CO aN aO H C ++→++ a=x+y/4=2+6/4=3.5 A/F|S =4.76aMW air /MW fuel =4.76(3.5)28.85/30.069=15.98 059.098 .1511,=+=stoich F Y stoich F F f Y Q D L ,83π= Sc=1 D=ν=μ/ρ ρ=ρair =1.16kg/m 3 26/105.96 2m S N H C -?=-μ 27/106.184m S N air -?=-μ s m D /1005.12]2 105.9106.184[16.11/266 7---?=?+?===ρμν s m m s m D v R v Q e e f /10927.34 )01.0()/05.0(4362 22-?====πππ m L f 66.0059.010927.3) 1005.12(836 6=??=--π (2) 关联式 ) /11ln()/(1330exp ,S T T Q L F F t f +=∞
第六章 预混层流火焰
第六章 层流预混火焰传播 §6-1 火焰速度和火焰结构 一维层流火焰在预混燃料-氧化剂混合物中传播是最简单的燃烧现象之一，在此火焰 中，化学动力学以及能量和组分扩散输运起重要作用。通过守恒方程和状态方程可以导出Rankine-Hugoniot 曲线。该曲线把在一维层流预混火焰中未燃气和已燃气状态联系起来。已燃气体位于Rankine-Hugoniot 曲线下分支（缓燃），并相应于未燃气体状态Rayleigh 线与具有适当反应热的Rankine-Hugoniot 曲线交点L ，如图6.1-6.2中所示。 图6.1 层流预混火焰坐标系 图6.2 一维燃烧波的Rankine-Hugoniot 曲线和Rayleigh 线 Rayleigh 线的斜率与相对于未燃气体的波的传播速度，即层流火焰速度有关。 22)()/(/u u u A m dv dP ρ?=?=& ==)(u u S u 层流火焰速度=)/()/1(dv dP u ρ? 由于缓燃Rayleigh 线斜率比爆震Rayleigh 线斜率小得多，所以缓燃速度比爆震速度小 得多。 虽然守恒方程和状态方程提供了缓燃的未燃气体和已燃状态之间的关系，但不能唯一确定层流火焰速度u S 。为了确定u S ，必须将守恒方程通过缓燃波积分。由于在第5章中推
导的方程是非线性耦合微分方程，其准确解只有通过数值积分才能获得。它需要很大的计算资源。为了考察层流火焰的某些特征（如火焰速度和厚度）以及这些特征与燃烧参数如燃料类型、化学配比、压力及未燃气体的温度的关系，对方程组进行了简化，以便能分析求解。要得到简化的模型，需要引入一系列的假设。我们从考察参考系建立在火焰上的层流火焰结构的某些方面入手。如前所述，这些计算是针对等压过程进行的。但是对一维缓燃的Rankine-Hugoniot 曲线，如图6.2所示，已燃气的压力小于未燃气的压力。现在我们需要考察压力减少的数值是否小到可以忽略的程度。如果能假设压力近似不变，则可以减少一个需要求解的方程数，动量方程将减少到P=常数。 对于稳态一维燃烧波，质量守恒方程变成： 常数=?=u dx u d ρρ0/)( 忽略粘性影响和体积力（浮力），动量方程可写成： 0)/(/=+dx du u dx dP ρ 应用以上两个方程估算通过火焰的压力降， [][] 1)/(1)/()()/(22 ?=??≈Δ?=Δ=ΔΔΔ?≈Δb u u u u b u u u b u u u u u u u u P u u u u u x x u u P ρρρρρρρ 由理想气体状态方程， )/(~)/)(/)(/(/u b u b u b b u b u T T T T R R P P =ρρ 由于反应物与产物的分子量近似相同，预期穿过火焰的压力降与温度增加相比是很小的，因此 []1)/(2 ??=Δu b u u T T u P ρ 碳氢燃料与空气混合物在大气条件下的层流火焰速度典型值在15-40cm/s 范围内。 u b T T /的典型值在5-7范围内，u ρ的典型值等于33/101cm g ?×。因此P Δ的典型值为： 2650.1~1/(10~10)P N m atm ???Δ= 因此，忽略通过火焰的压力降是很合理的。 §6-2 一维层流预混火焰模型 描述一维预混火焰的方程组是：压力为常数的条件下的质量守衡（=u ρ常数），能量守恒（H=常数），组分守恒以及理想气体状态方程。如附录G 中公式（G-40），（G-41）所示， 若利用Le=1近似，可以假设能量方程和组分方程为同一种形式： []RR dx dx d D u d =?/)/(ηρηρ 式中： [] )(/)()] /[('ref R sebsible T i i i i T H T h Y Δ?=?=ηυυη 只要求解能量方程和组分方程中的一个，加上适当的边界条件，就可以完全决定火焰的结构和速度。
第章预混燃烧模拟
第十五章预混燃烧模拟 FLUENT 有一个预混湍流燃烧模型，基于反应过程参数方法。有关这一模型的内容按以下节次给出： 15.1概述和限制 15.2 预混燃烧模型 15.3 使用预混燃烧模型 15.1 概述和限制 15.1.1概述在预混燃烧中，燃料和氧化剂在点火之前进行分子级别的混合。火焰前锋传入未燃烧的反应物产生燃烧。预混燃烧的例子有吸气式内燃机，稀薄燃气轮机的燃烧器，气体泄露爆炸。 预混燃烧比非预混燃烧更难以模拟。原因在于（亚音速）预混燃烧通常做为薄层火焰产生，并被湍流拉伸和扭曲。火焰传播的整体速率受层流火焰速度和湍流涡旋控制。层流火焰速度由物质和热量逆流扩散到反应物并燃烧的速率决定。为得到层流火焰速度，需要确定内部火焰结构以及详细的化学动力学和分子扩散过程。由于实际的层流火焰厚度只有微M 量级或更小，求解所需要的开销是不可承受的。 湍流的影响是使传播中的层流火焰层皱折、拉伸，增加了薄层的面积，并因此提高了火焰速度。大的湍流涡使火焰层皱折，而小的湍流涡，如果它们比层流火焰的厚度还小，将会穿过火焰层并改变层流火焰结构。 与之相比，非预混燃烧可以极大地简化为一个混合问题（例如，14.1节中介绍的混 合物组分方法）。预混燃烧模拟的要点在于捕获湍流火焰速度，它受层流火焰速度和湍流的影响。 在预混火焰中，燃料和氧化剂在进入燃烧设备之前已经紧密混合。反应在燃烧区发生，这一区域将未燃烧的反应物和燃烧产物隔开。部分预混火焰具有预混和扩散火焰两方面的性质。它们发生在有额外的氧化剂或燃料气流进入预混系统，或是当扩散火焰离开燃烧器以在燃烧前产生某些预混的情况。 预混和部分预混火焰FLUENT的有限速率公式（见13章濮拟。还可以参阅16章了解更多有关FLUENT部分预混燃烧模型方面的信息。如果火焰是完全预混合的，则只有一股具有单一混合比的气流进入燃烧器，可以使用预混燃烧模型。 15.1.2限制 在使用预混燃烧模型时有以下限制：必须使用非耦合求解器。预混燃烧模型在两种耦合求解器中都不能得到。预混燃烧模型只对湍流、亚音速模型有效。这一类型的火焰成为爆燃。在爆炸中，可燃混合物被冲击波后面的热量点燃，这一类型的燃烧可以使用非耦合和耦合求解器用有限速率模型模拟。有关限速率模型见13章。预混燃烧模型不能和污染物（如碳烟和 NOx ）模型一起使用。但完全预混系统可以用部分预混模型（见16 章）模拟。 不能用预混燃烧模型模拟反应的离散相粒子。只有惰性粒子可以使用预混燃烧模型。 15.2预混燃烧理论 湍流预混燃烧模型基于Zimont 等人的工作[275，276，278]，涉及求解一个关于反应过程变量的输运方程。这一方程的封闭基于湍流火焰速度的定义。 15.2.1 火焰前锋的传播 在许多工业预混系统中，燃烧发生在一个非常薄的火焰层中。当火焰前锋移动时，未燃的反应物燃烧，变为燃烧产物。因此预混燃烧模型用火焰层将反应的流场分为已燃物区和未燃物区。反应的传播等同于火焰前锋的传播。 火焰前锋传播的模拟通过求借一个关于标量c的输送方程，c为(Favre平均)反应进 程变量。
湍流燃烧及其数值模拟
湍流燃烧及其数值模拟研究 1. 湍流燃烧 1.1湍流燃烧基本概念 当流动雷诺数数较小时，由于流体粘性的作用，流体呈层流流态。当流动的特征雷诺数超过相应的临界值，流动从层流转捩到湍流。湍流燃烧是指湍流流动中可燃气的燃烧，在能源、动力、航空和航天等工程领域,经常遇到的实际燃烧过程几乎全部都是湍流燃烧过程。湍流燃烧实质是湍流，化学反应和传热传质等过程相耦合的结果。湍流对燃烧的影响与湍流强度和湍流涡旋尺度有关。小尺度湍流通过湍流扩散使火焰区内的输运效应增加，从而使化学反应速率增加。但气流脉动不会火焰面产生皱褶，只能把火焰变成波纹状。大尺度湍流对火焰内部结构没有影响，但使火焰阵面出现皱褶，增加其燃烧面积，造成火焰表现传播速度增加。当湍流强度及湍流尺度均较大时，火焰前沿不再连续而分裂成四分五裂。 燃烧对湍流的影响主要表现在燃烧释放的热流流团膨胀，影响气体的密度和运动速度，从而影响当地的涡旋，湍流强度和湍流结构。 1.2湍流燃烧分类 湍流燃烧按其燃料和氧化剂的初始混合状态可以分类为：湍流非预混燃烧、预混燃烧和部分预混燃烧。在湍流非预混燃烧燃料和氧化剂事先是分离的，燃料和氧化剂一边混合一边燃烧，燃烧速率主要受湍流混合过程控制，而在湍流预混燃烧中，燃料和氧化剂在进入核心燃烧区以前已经充分混合，化学反应的速率由火焰前缘从炽热的燃烧区向冷态无反应区的传播所控制。上面两种燃烧方式是湍流燃烧的两个极限情形，很多情况下两种燃烧模式是并存的，称为部分预混燃烧。部分预混燃烧可出现在下列情形中叫：(1)在一个完全以非预混燃烧为配置的燃烧装置发牛了局部熄火；(2)当预混火焰前缘穿过非均匀的混气时；(3)射流非预混火焰发生抬举，其根部是一个典型的部分预混火焰。这三种部分预混燃烧情形涉及了经常受到关注的燃烧研究话题如局部熄火、火焰稳定等，它们对研究湍流燃烧过程的机理有很大意义。 在湍流燃烧中,湍流流动过程和化学反应过程有强烈的相互关联和相互影响.湍流通过强化混合而影响着时平均化学反应速率,同时化学反应放热过程又影响着湍流,如何定量地来描述和确定这种相互作用是湍流燃烧研究的一个重要内容. 湍流是非常复杂的，它包括湍流问题，湍流与燃烧的相互作用，流动参数与化学动力参数之间的耦合机理等问题。因此湍流燃烧是工程科学中最复杂的领域之一。 湍流燃烧的研究已进行多年，研究的方法有试验研究，理论分析和数值模拟等。计算流体力学和计算机技术的发展，数值模拟由于它的廉价性和可操作性在国际上受到越来越多的重视，得到了广泛的应用。 2.湍流燃烧数值模拟 2.1湍流燃烧数值模拟简介 湍流燃烧数值模拟(Numerical Simulation of Turbulent Combustion)是指应用计算机为工
第14章 非预混燃烧模拟
第十四章非预混燃烧模拟 Chapter 14. Modeling Non-Premixed Combustion 在非预混燃烧中，燃料和氧化剂以相异流进入反应区。这与预混燃烧系统截然不同。在预混燃烧系统中，反应物在燃烧以前以分子水平混合。非预混燃烧的例子包括甲烷燃烧、粉煤炉和内部燃烧柴油（压缩）发动机。 在一定假设条件下，热化学可被减少成一个单一的参数：混合分数。混合分数，用f表示，是来自燃料流的质量分数。换句话说，混合分数就是在所有组分（CO2、H2O、O2等）里，燃烧和未燃烧燃料流元素（C、H等）的局部质量分数。因为化学反应中元素是守恒的，所以这种方法极好。反过来，质量分数是一个守恒的数量，因此其控制输运方程不含源项。燃烧被简化为一个混合问题，并且与近非线性平均反应率相关的困难可以避免。一经混合，即可用层流小火焰（laminar flamelet）模型将化学反应模拟成为化学平衡或近化学平衡。 模型包括以下几个部分： 14.1：平衡混合分数／PDF模型（Description of the Equilibrium Mixture Fraction/PDF Model）； 14.2：非预混平衡化学反应的模拟方法（Modeling Approaches for Non-Premixed Equilibrium Chemistry）； 14.3：非预混平衡模型的用户输入（User Inputs for the Non-Premixed Equilibrium Model）； 14.4：层流小火焰模型（The Laminar Flamelet Model）； 14.5：在prePDF数据库中添加新种类（Adding New Species to the prePDF Database）； 14.1：平衡混合分数／PDF模型 非预混模拟方法包括解一或两个守恒量（混合分数）的输运方程。不解单个组分方程。取而代之的是每个组分的浓度用预混分数场得到。热化学计算在prePDF中进行，并列成表以便于在FLUENT中查询。紊流和化学的相互作用考虑为一个概率（几率）密度函数（PDF）。 关于非预混混合分数／PDF模型的信息在以下分节中讲述： 14.1.1：非预混方法的优点和局限（Benefits and Limitations of the Non-Premixed Approach）； 14.1.2：非预混方法的细节（Details of the Non-Premixed Approach）； 14.1.3：非预混模拟的限制和特有案例（Restrictions and Special Cases for Non-Premixed Modeling）； 见14.2：模拟和解决顺序的回顾，以及14.3；应用模型指导。 14.1.1非预混方法的优点和局限 非预混方法的优点（Advantages of the Non-Premixed Approach）：非预混模拟方法已被明确用于模拟进行快速化学反应的紊态扩散火焰的研究。对这样的系统，该方法有许多点优于第十三章中描述的有限率公式。非预混
火焰法向传播速度
用本生灯法测定石油液化气体火焰法向传播速度的实验 一、实验所需的器材和原理 实验台由移动式不锈钢实验台、空气泵、流量计、压力计、本生灯、液化石油气瓶等组成，实验中我们通过本生灯发或称动力法进行测定。 首先，法向火焰传播速度我们可以理解为垂直于层流火焰前沿面方向上火焰前沿面相对于未燃混合气的运动速度。 我们将气泵产生的空气通过泻流阀、稳压罐、流量计、调压阀后进入本生灯，同时可燃气体经减压器、流量计、防回火器、调压阀后进入本生灯与空气预混合点燃，在理想的稳定的燃烧火焰中，混合气流的法向分速度与未燃混合气流的运动速度即法向火焰传播速度相平衡，这样才能保持燃烧前沿面在法线方向上的燃烧速度后通过测量内焰锥高度计算火焰的传播速度
图一 火焰传播速度测试原理 0sin s u u α=? 式中：u s －混合气体的流速（CM/S ）； α－火焰锥角之半。 或 0q u =式中：q v －混合气体的体积流量（L/S ）； h －火焰内锥高度（cm ）； r －喷口半径（cm ）。 在我们的实验中，可燃气体即液化石油天然气的体积流量v q 是用流量计分别测定燃气与空气的体积流量而得到的，内锥焰面底部圆的半径r 可取本生灯喷口半径；内焰锥高度h 可由测高尺测量。 二、实验的操作步骤 1，检查实验所用到的器具是否完好，进行气密性实验。 2，打开电源启动气泵，调节本生灯出口大小使流速稳定，用流量计测定空气流量。 3，打开燃气阀，点燃火焰，这时火焰会呈扩散式燃烧；慢慢开启空气泵调节阀，送入空气。当混合气体发生完全燃烧出现火焰内锥时（即呈现完全的蓝色火焰并且火焰面三角形清晰可见），用气量计测出燃气与空气的体积流量，测高仪测得火焰内锥高度（从火焰底部，即喷口出口断面处到火焰顶部间的距离）。 4，多次适当增加或减少空气量，即改变一次空气系数，测出相应的火焰内锥高度。将测得的数据填入测试表中。根据理想气体状态方程式(等温)，将燃气和空气测量流量换算成(当地大气压下)喷管内的流量值，然后计算出混合气的总流量，求出可燃混合气在管内的流速u s ，并求出燃气在混合气中的百分数。记录室温，计算出0u 值。 三、结论 通过用本生灯法测定石油液化气体火焰法向传播速度的实验，了解火焰传播的的形式。并得出可燃混合流量与火焰传播速度的关系。
层流预混火焰传播速度测定
层流预混火焰传播速度的测定 实验成员：徐俊卿 郑仁春 韩超 一、实验的理论基础 许多工业设备都应用预混气燃烧作为热和能量的生成方式。如火花点火发动机（汽油机），煤气炉内的燃烧，灾害性的火灾和爆炸都涉及到预混气的燃烧和火焰传播问题。 研究预混气燃烧的最重要参数是层流火焰传播速度。火焰速度是预混气的基本特性，是研究火焰稳定性以及湍流预混气燃烧的基础。 层流火焰速度定义为给定可燃预混气的一维平面预混火焰在没有热损失时相对于未燃气的移动速度。用S 0表示。该定义给出的火焰速度是预混气的单一的固有特性，而与外界流动条件无关，在某些精心设计的实验设备，如相向流火焰设备上，采用激光多普勒速度仪，可以精确测定S 0。普通的预混火焰设备很难完成满足上述定义中的所有条件。如采用本生灯测定火焰速度，由于火焰面呈锥形，不是一维火焰，顶端和底部火焰有弯曲。不可避免地有热损失。因此测到的是被测点当地的火焰速度或称局部火焰速度，用S 表示。S 除与可燃预混气的气/油比有关外，还受热损失，火焰拉伸等动力学因素影响。用其它的实验方法，如平面火焰法，火焰推进法，肥皂泡法，球弹法和圆管法都是只能测定局部火焰速度。 层流火焰理论指出，预混火焰的稳定位置总是位于预混气在火焰面的法向速度分量与火焰速度（总与火焰面垂直）大小相等，方向相反的地方。当这两个速度不相等时，火焰面就要移动，而扩散火焰总是驻定在燃料与氧化剂为化学计量值的位置上。在这一位置，燃料与氧化剂混合最均匀，反应率最快。偏离这一位置，不可能组织起燃烧，扩散火焰没有火焰传播速度的概念，这是预混火焰和扩散火焰最主要的区别之一。 二、实验原理 实验采用本生灯测定（局部）火焰传播速度，实验设备与实验二相同。 设计良好的本生灯火焰呈锥形，除顶端和底部火焰弯曲外，中间有较长一段的平直火焰，假定预混气速度沿出口截面分布均匀，火焰前沿各处的气流法向速度相等，把驻定在管口的火焰面简化为正锥形，如图3-1所示。预混气的速度为u 0，火焰面平直的上点P 的火焰速度为S 。S 在数值上等于u 0在P 点垂直于火焰面的法向分量，即 S u =0sin α （1） α为火焰锥的半顶角。可用测高仪测出火焰高度和底部直径后算得，也可用量角器直接量出。 u 0由浮子流量计测定流量后，根据出口直径计算（d 0为管口内径）： Q u A d u == 004 02 0π （2） S ∠ 图2-1本生火焰 三、实验设备 预混火焰装置 五、数据数据及处理
fluent中燃烧模型分类
FLUENT燃烧模型 化学反应 模拟方法 方法描述计算反应的选择 有限速率模型需要求解组分质量分数的 输运方程，化学反应机理 由用户自己定义。反应速 率在组分输运方程中作为 源项，并由阿累尼乌斯公 式计算。应用范围最广泛。 应用：模拟化学组分混合、 输运和反应的问题；壁面 或粒子表面反应问题 层流有限 速率模型 使用Arrhenius公式计算化学源项，忽略湍流脉动的影响。对于化学动力学控制的燃烧（如层流燃 烧），或化学反应相对缓慢的湍流燃烧是准确的。但对一般湍流火焰中Arrhenius化学动力学的高度 非线性一般不精确；对于化学反应相对缓慢、湍流脉动较小的燃烧（如超音速火焰）可能可以接受。 漩涡破碎模型 Eddy Dissipation 大部分燃料快速燃烧，整体反应速率由湍流混合控制。复杂且常是未知的化学反应动力学速率可以 完全的被忽略掉。化学反应速率由大尺度涡混合时间尺度k/ε控制。只要k/ε（湍流）出现，燃烧 即可进行，不需要点火源来启动燃烧。（缺点：未能考虑分子输运和化学动力学因素的影响） 适用条件：高雷诺数湍流预混燃烧过程。 EBU-Arrehenius 模型 EDC模型 假定化学反应都发生在小涡中（精细涡），反应时间由小涡生存时间和化学反应本身需要的时间共 同控制。EDC模型能够在湍流反应中考虑详细的化学反应机理。但是他们的数值积分计算开销很大。 使用条件：只有在快速化学反应假定无效的情况下才能使用这一模型（如快速熄灭火焰中缓慢的 CO烧尽、选择性非催化还原中的NO转化问题）。 非预混燃烧模型不求解每个组分的质量分数输运方程，求解混合分数输运方程和一个或两个守恒标量的方程，然后从预测的混合分数公布推导出每一个组分的浓度。通过概率密度函数或PDF来考虑湍流的影响。 应用：主要用于模拟湍流扩散火焰的反应系统。这个系统要求接近化学平衡，氧化物和燃料以两个或者三个进口进入计算域。 预混燃烧模型主要用于单一、完全预先混合好的燃烧系统。反应物和燃烧产物被火焰前沿分开。求解出反应发展变量来预测前沿的位置。湍流的影响通过湍流火焰速度计算。 部分预混燃烧模型描述非预混燃烧完全预混燃烧相结合的系统。结合混合分数方程和反应物发展变量来分别确定组分浓度和火焰前沿位置。适用于计算域内具有变化等值比率的预混火焰情况。通过求解混合分数方程和反应过程参数来确定火焰峰面的位置。 PDF输运方程模型结合CHEMKIN可以考虑详细的化学反应机理，高度的非线性化学反应项是精确模拟，无须封闭模型，可以合理的模拟湍流和详细化学反应动力学之间的相互作用，是模拟湍流燃烧的精确模拟方法。但计算量特别大。 优点：可以计算中间组分；考虑分裂影响；考虑湍流-化学反应之间的作用；无需求解组分输运方程 缺点：系统须满足（靠近）局部平衡；不能用于可压缩或非湍流流动；不能用于预混燃烧。
实验一 层流预混火焰稳定浓度界限测定
实验一层流预混火焰浓度稳定界限的测定 一、实验的理论基础 现代工业要求尽可能高效地利用燃料，这就要求在稳定燃烧条件下，使燃料与空气以适当的比例混合燃烧。因此对燃料燃烧特性及流体力学和热力学方面需要有详细了解。 预混可燃气燃烧时，如果预混气体的速度在火焰锋面上的法向分量大于火焰传播速度，火焰将向下游移动，最后完全熄灭，称为吹脱或吹熄。反之，如果预混气的法向速度小于火焰传播速度，火焰将逆流向上游移动，进入燃烧器内部，即出现回火现象。 在燃烧过程中，出现回火和熄灭都是不允许的。回火会引起爆炸。熄灭使动力机械停止工作，并向周围扩散有毒气体，有中毒和爆炸的危险。 回火现象只能出现在预混燃烧过程中。在扩散燃烧中，燃料和空气是分别送入燃烧室的，在燃烧器内两者并不接触，因此没有回火现象。熄灭或吹脱现象在预混燃烧和扩散燃烧中均有可能出现。 回火煤气 8吹脱甲烷 4回火 黄尖 2 10203040V F （l/min 图 1-1甲烷/空气煤气/空气稳定曲线 要维持正常的稳定燃烧，就需要避免出现回火或熄灭现象，因而要求知道燃料的稳定燃烧范围。这一稳定界限与燃料/空气比和环境的温度和压力有关。图1—1是甲烷/空气煤气/空气在常温常压下预混燃烧的火焰稳定范围。图中下部是富油燃烧区，稳定范围宽，不易回火和吹脱。但火眼呈黄色，发光。下面一根曲线为火焰尖顶出现黄色的界限。图中上方为火焰吹脱熄灭界限。当预混气处于该曲线上方时，火焰被吹脱。 左边回火界限呈斜山立形，预混气的状态处于山丘内时出现回火现象。回火界限与预混气的空气气/燃气比有关。当空气/燃气比接近化学计量值时，火焰传播速度快，稳定范围小，容易出现回火。但这时燃烧效率最高。 二、实验原理 本实验用一套预混火焰装置（本生灯）进行实验，既可以定性观察预混火焰的各种现象，又可以定量地测定火焰浓度稳定的界限，火焰传播速度等各种燃烧特性。 图1-2实验系统图示意图 图1—2为实验装置系统图，小型压气机供给压缩空气，通过一个稳压减压阀保持气压稳定。流量由浮子流量计读数。燃料可用丁烷气、城市煤气或石油液化气等。由浮子流量计控制燃料流量。空气和燃料分别送入混合块进行混合。为保证燃料空气充分混合，并保持实验