湍流预混火焰模型(1)
第六讲 预混合气体火焰-1
f
f
即为Hugoniot曲线上任意点的切线斜率。
（3）Raleigh-Hugoniot方程
• 火焰锋面后方的熵变：
∵
dq du pdV
T f dS f de f Pf dV
dS f Ps Pf 1 1 1 dPf Tf ( ) 1 1 1 1 d ( ) 2 s f d f s f f
• 或将u f 代入式（2）
f ( Ps Pf ) u s ( s f )
2 s
s2us2
s f ( Ps Pf ) ( s f )
6.1.2 爆燃与缓燃-基本方程
• 同理： •
Ps Pf u u 1 1 f s
2 2 s s 2 2 f f
过点S的直线（且 C <0）在第一、三象限不复存在。 凡满足瑞利方程的均为过点S的直线簇。
C 0 ， s Pf ，等压线， P
m 0 。
等密线，m 。
（ C ，Pf
s f ），
（2）Hugoniot（雨果尼特）方程
将式（1）与能量方程耦合
递的热量，提高自身温度，火焰逐渐向未燃区传播。 燃烧前后反应物浓度、温度变化情况见上图示意。 描述燃烧波中这些量得变化称为火焰结构。 预混火焰可以是固定的，也可以是以一定速度传播 的。
6.1.2 爆燃与缓燃
• 一根很长的等截面水平管内充满可燃混气，在左端 点燃，燃烧波将以某恒定的速度向右传播。
6.1.2 爆燃与缓燃
s f
( Ps
C
s
)
C
f
(4)
• 瑞利方程另一种形式。
湍流燃烧火焰面模式理论及应用(孙明波,白雪松,王振国著)PPT模板
0 5
2.5湍流预混 燃烧算例验证
0 6
2.6带自点火 特性的预混火 焰传播模型
第2章湍流预混燃 烧
参考文献
第2章湍流预混燃烧
2.1层流预混火焰
2.1.1层流 预混火焰结 构
2.1.2层流 预混火焰温 度
第2章湍流预混燃烧
2.2湍流预混火焰
0 1 2.2.1湍流预混火焰的基本性质
02
2.2.2湍流脉动与火焰的相互作 用
第1章湍流燃烧及其数值模拟概述
1.1湍流燃烧基本特性
1.1.1湍流 的基本特 性
1.1.2湍流 燃烧的特 点
第1章湍流燃烧及 其数值模拟概述
1.2化学反应流的数学描 述
1
1.2.1化学反应流控制方程
2
1.2.2化学反应机理及反应速率
第1章湍流燃烧及 其数值模拟概述
1.3湍流燃烧模拟的一般方 法
2.4.4G方程 和C方程比较
第2章湍流预混燃烧
2.5湍流预混燃烧算例验证
1
2.5.1均匀各向同性湍流中的火 焰核增长
2
2.5.2三角棱柱火焰稳定器的燃 烧模拟
3
2.5.3低旋流燃烧器的火焰稳定
4
2.5.4本生灯的火焰形状
第2章湍流预混燃烧
2.6带自点火特性的预混火焰传播模型
2.6.1预混 火焰自点火 耦合模型
n解和化学 平衡解
04
03
3.2.4火焰面结构的 渐近解
3.2.3详细化学反应 机理对层流扩散火 焰的影响
第3章扩散燃烧
3.3湍流扩散燃烧火焰面模型
01 3 .3 .1 扩散火焰 面模 02 3 .3 .2 火焰面模 型方
型合理性验证
浙大高等燃烧学_湍流燃烧理论模型_程乐鸣_2013_9
率决定于末燃气团在湍流作用下破碎成更小气团的速
率，而破碎速率与湍流脉动动能的衰变速率成正比。
R fu ~ / k
湍流燃烧速率
对比用k - ε模型和混合长度模型计算湍流粘度的公式
t C k 1/ 2 C C D k 2 /
假定 k 1/ 2 正比于混合长度与均流速度梯度绝对值的乘
对于层流火焰，在一定条件下，火焰传播速度与试验装 置无关。
在研究湍流燃烧时，针对湍流火焰，同样期望确定其传 播速度时，不要与装置本身有关，以带有共性，仅与料量比： λ、μ、D等量数有关。 事实证明这是不可能的。
在某些化学当量比下，湍流中有效热扩散系数要比层流 中分子的热扩散系数大 100倍，因此，湍流火焰的理论概念 不象层流火焰那样容易定义。
分析湍流火焰时，不仅要考虑湍流的 输运特性，还必须考虑湍流的脉动特性。 建立湍流燃烧模型中，要把混合过程 的控制作用和湍流脉动的影响有机地统一 起来。 基于此，Spalding提出了k-ε-g模型
几率分布函数
几率分布函数，即：一个用于描述湍流燃烧系 统中的因变量。 对于某个量我们关心的是它取某个值的几率。 无量纲混合分数的几率分布函数定义如下： P(f)df=f(t)处于(f,f+df)范围内的那段时间间隔t的 时间分数，即几率。 式中，P(f)称为瞬态混合分数f的几率分布密度 PDF。
F Sl w0 FL
F ST w0 FT
湍流火焰锥外 表面面积
研究湍流火焰过程中发展起来的方法
一类为经典的湍流火焰传播理论，包括皱折层流火焰的 表面燃烧理论与微扩散的容积燃烧理论。 另一类是湍流燃烧模型方法，是以计算湍流燃烧速率为 目标的湍流扩散燃烧和预混燃烧的物理模型，包括几率 分布函数输运方程模型和ESCIMO湍流燃烧理论。
湍流预混火焰模型概要
(3-56)
流场较均匀的区域，合理地估算层流火焰传播速度是 正确运用拉切滑模型的关键之一 。 层流火焰传播速度是可燃气的物理化学性质，它取决 于混合物的热力学状态(如压力和温度)，对温度尤为 敏感。 丙烷和空气当量比混合物的火焰传播速度 S 0.113(T / 298)2 0.186(T / 298) 0.02 (m/s) (3-57) 求S的问题转化为求T。
反应度τ的脉动均方值 gτ
定义
g 2 ( ) 2 (m fu m fu ) 2 /(m fu ,b m fu ,u ) 2 m /(m fu ,b m fu ,u )
2 fu 2
(3-60)
2 gτ和 mfu 应当遵守同一类型的微分方程。
反应度τ的平均值和脉动均方值 gτ的确定
对比用k - ε模型和混合长度模型计算湍流粘度 t C k 1/ 2 C CD k 2 / 的公式 2) 假定 k 1/ 2正比于混合长度与均流速度梯度绝对 值的乘积 则ε/k正比于均流速度梯度的绝对值 3) 燃烧速率一定与燃料浓度有关 二维边界层问题湍流燃烧速率 u (3-42) R c m
5 平面管道内火焰稳定器后面的燃烧场
Spalding et al 结果优于只用 阿伦纽斯类型 的公式(3-48)得 到的结果，与 实验数据的趋 势符合
6 对旋涡破碎模型的评价
功绩在于正确地突出了流动因素对燃烧速率的 控制作用，给出了简单的计算公式，为湍流燃 烧过程的数学模拟开辟了道路。 不足：该模型未能考虑分子输运和化学动力学 因素的作用 适用范围：一股说来，EBU模型只适用于高雷 诺数的湍流预混燃烧过程。
湍流燃烧模型-PDF
PDF 模型概率密度函数PDF方法以随机的观点来对待湍流问题，对解决湍流化学反应流的问题具有很强的优势。
在湍流燃烧中存在一些非输运量( 如反应速率, 密度, 温度及气相体积分数等) 的湍流封闭问题。
尽管这些量没有输运方程, 但它们常常是输运变量的已知函数。
平均或者过滤高度非线性的化学反应源项会引起方程的封闭问题。
因此，用PDF的方法来解决这些非输运量的湍流封闭问题显然是一个既简单又直接的途径。
PDF方法是一种较为流行的湍流燃烧模型，能够较为精确的模拟任何详细的化学动力学过程, 适用于预混、非预混和部分预混的任何燃烧问题。
目前, 确定输运变量脉动概率密度函数的方法有输运方程和简化假定两种, 分别称之为输运方程的PDF和简化的PDF。
前者建立输运变量脉动的概率密度输运方程，通过求解该方程来获得输运变量脉动的概率分布。
后者假定输运变量脉动的概率密度函数的具体形式, 通过确定其中的一些待定参数来获得输运变量脉动的概率分布。
湍流燃烧中, 后者应用最为普遍和广泛。
在简化的PDF 中, 输运变量脉动的概率密度函数常常采用双 D 分布、截尾高斯分布和B 函数分布等形式。
PDF在理论上可以精确考虑任意详细的化学反应机理，但是其具体求解时需借助其它的模型和算法，而且计算量相对较大。
PDF的方程是由N-S方程推导而来，其中的化学反应源项是封闭的，但压力脉动梯度项以及分子粘性和分子扩散引起的PDF的分子输运项是不封闭的，需要引入模型加以封闭。
例如，在速度- 标量-湍流频率PDF中，必须采用小尺度混合模型、随机速度模型和湍流频率模型加以封闭。
模化后的输运方程难以用有限容积、有限差分和有限元等方法来求解，比较可行的一种方法是蒙特卡洛(MonteCarlo)方法,在该方法中输运方程被转化为拉格朗日(Lagrangian)方程,流体由大量遵循Lagrang ian方程的随机粒子的系统来描述, 最后对粒子作统计平均得到流场物理量和各阶统计矩。
第五章fluent预混燃烧模型ppt课件
1、预混模型理论
• 火焰前锋的传播：预混燃烧时，火焰发生在一个 非常薄的火焰层中，火焰前锋移动时，未燃反应 物燃烧变为产物，火焰层将反应的流场分为已燃 物区和未燃物区，反应的传播等同于火焰前锋的 传播
预混燃烧--Zimont模型
•
反应进程变量c：c
Yp
/
Y
ad p
p
p
Yp：当前产物的质量分数；
Ypad ：完全绝热燃烧后产物的质量分数；
预பைடு நூலகம்模型总结
• 适用条件
湍流 快速化学反应 只有预混合
• 限制条件
不能模拟运动学细节中的实际现象 （如点燃、熄灭和低Da数）。
实例演练四：预混燃烧
混合燃料入口2
混合燃料入口1
烟气出口
湍流长度尺寸常数CD 湍流火焰速度常数A
拉伸系数 湍流施密特数Sct
拉伸系数
• 为了考虑火焰面拉伸所导致的吹熄现象，在 反应源项中乘以一个拉伸因子 G，即GSC ：
其中：
• 以上各式中出现的一些常数值在FLUENT默认条件下为：
A=0.52，CD=0.37，μstr=0.26， Sct =0.7
温度的计算
• 关键：捕获湍流火焰速度，受层流火焰速度和湍流的 影响。
预混模型使用限制
• 必须使用非耦合求解器； （define-models-solver: Pressure based） • 只对湍流、亚音速模型有效； • 不能和污染物模型（如NOx）一起使用； • 不能模拟离散相粒子的反应，只有惰性粒
子才能与预混模型一起使用。
的情况。
3、FLUENT相关设置
1、选择预混模型
2、确定绝热或非绝热
（如果有fluent材料库 中的模型，可以首先选 择一种）
流动及燃烧的模型及计算
第一章绪论1.1计算燃烧学的研究对象和目标(1)研究对象：对流体流动、传热传质和燃烧过程进行计算机模拟的基本方程（连续方程、动量方程、能量方程、组分方程）、理论模型（湍流输运、湍流燃烧、辐射换热、多组分化学反应和多项问题）、数值方法（研究体系的网格化、控制方程的离散化和求解方法）计算机程序（计算程序、绘图程序、仿真程序）(2)计算燃烧学的研究目标：构造、检验和发展基本方程及理论模型，提高他们的可靠性、准确性和实用性；改进数值方法，在保证计算精度的同时提高计算速度和经济性；改善绘图及仿真软件，通告速度和直观性；提供“使用方便性”强的计算、绘图及仿真软件，，方便使用。
1.2 计算燃烧学的意义(1)使燃烧上升到系统理论(2)是设计、科研和教学的手段(3)有助于学科发展和开拓新领域1.3计算燃烧学的发展简史燃烧的定义：燃烧室一种带有剧烈放热化学反应的流动现象，它包含着流动、传热、传质和化学对流体流动、传热传质和燃烧过程进行计算机模拟的基本方程（连续方程、动量方程、能量方程、组分方程）理论模型（湍流输运、湍流燃烧、辐射换热、多组分化学反应和多项问题）数值方法（研究体系的网格化、控制方程的离散化和求解方法）计算机程序（计算程序、绘图程序、仿真程序）在科研、工程和教学中的应用反应以及它们之间的相互作用。
实际燃烧过程几乎都是湍流过程。
计算燃烧学特点：(1)兼顾研究理论模型和数值方法(2)及研究理论模型和数值方法，有中式计算机程序的编制和更新。
(3)既重视通用模型、通用方法和通用程序的研究，又重视与实际应用得密切结合。
1.4 内容梗概(1)构造基本方程和理论模型燃烧过程涉及的基本定律：物质不灭定律、牛顿第二定律、能量转换和守恒定律、组分转换和平衡定律等。
控制燃烧过程的基本方程组：连续方程、动量方程、能量方程、组分方程需要模化的分过程有：湍流输运、燃烧、辐射转换、多项流动和燃烧。
(2)数值方法数值方法主要包括：研究体系的网格化、控制方程的离散化、离散化方程的求解方法。
燃烧学-4.预混合气燃烧及火焰传播
uz uzuz
流场中某一点总湍流强度：
u
1 3
ux2
uy2
uz2
湍流尺度(turbulent scale)
与湍流涡团大小及其变化过程有直接关系。 物理意义：涡团在无规则运动中，保持自由前进而不与其 他涡团碰撞的距离。或者说，流体涡团在运动过程中消失 前运动的距离，也可以认为是涡团的一种平均自由程。
第四章 预混合气燃烧及火焰传播
层流火焰概念、结构特征、传播机理、传播速度计算，层 流火焰传播速度影响因素 ，湍流火焰概念 ，湍流火焰传播理论 与传播速度，爆震燃烧理论。
层流火焰结构、传播机理，湍流火焰传播两种理论
层流火焰传播的数学模型建立与推导，湍流火焰传播理论
概述
一、预混合燃烧概念
定义
燃料和氧（或空气）预先混合成均匀的混 合气，此可燃混合气称为预混合气，预混合气 在燃烧器内进行着火、燃烧的过程称为预混合 燃烧（premixed combustion）。
五、层流火焰传播速度影响因素分析 S lC p 2 2 2T n f!W T Q n 2 E E T f2 n 1 T T f ne x p E R T 1 f T 1
Lagrange湍流尺度： L1u0RduT
Euler湍流尺度：
L2 0 Rrdr
湍流扩散系数(turbulent diffusion coefficient)
反映了湍流传质和传热的特性。 依照分子扩散与湍流扩散的相似性，可以认为， 流体层中体积基元（涡团）的无规则运动与分子的无 规则运动相似，用下式定义湍流扩散系数：
二、层流火焰传播速度
6-湍流预混火焰讲解
湍流火焰的特点
均匀、各向同性的湍流流场，可以用两 个特征量表示湍流特征：湍流强度和湍 流尺度
湍流尺度：
（1）流动特征尺度（与管径、绕流物体尺度有关） （2）积分尺度（湍流宏观尺度，大涡尺度） （3）泰勒微尺度（与平均应变率有关） （4）柯尔莫戈洛夫尺度（最小尺度，与旋涡耗散有关）
湍流火焰的特点
小尺度湍流预混火焰传播速度确定
湍流火焰传播速度和层流火焰传播速度之比等 于二者传输率之比的平方根
ut un
T n
1/ 2
T n
/ 0cp / 0cp
1/ 2
λt表示湍流热传导系数 λl表示层流热传导系数 根据相似性原理，分子导温系数α= λn/(ρ0cp), 故 湍流导温系数αt= λT/(ρ0cp)。在湍流中湍流导温 系数取决于湍流尺度和脉动速度乘积，即
a）小尺度湍流火焰（2300<Re<6000） 条件： l<δL
现象：能够保持规则的火焰锋面，火焰前 沿仍然平滑，只是增加了厚度，火焰锋面 不发生皱折，湍流火焰面厚度δT> δn
特点：小尺度湍流只是由于湍流增强了物 质的输运特性，从而使热量和活性粒子的 传输增加，使湍流火焰传播速度比层流火 焰传播速度快，而在其它方面没有什么影 响
预混湍流火焰结构的分形特征
第33卷 第2期 1999年2月 西 安 交 通 大 学 学 报JO URNAL OF XI′AN JIAOTONG UN IVERSITY Vol.33 №2 Feb.1999预混湍流火焰结构的分形特征蒋德明,马凡华,杨 迪(西安交通大学,710049,西安)摘要:使用分形几何学的方法,对由高速纹影摄影法获得的容弹内预混湍流火焰照片进行了处理.结果表明,容弹内预混湍流火焰结构具有分形特征,可以用分数维数来定量描述其折皱、扭曲程度及与湍流强度的关系.关键词:火焰结构;分形几何;燃烧学中国图书资料分类法分类号:TK411The Fractal Nature of Tu rbulent PremixedFlame StructureJiang Dem ing,Ma Fanhua,Yang Di(Xi′an Jiaotong University,Xi′an710049,C hina)A bstract:The fractal geometry method is used to analyze the hig h speed schlieren film of turbulentpremixed combustion for a co nstant volume combustion bom b.The fractal dimension can be used to describe the ex tent of w rinkling deg ree quantitatively and to exhibits a relation for the root mean square turbulent velocity.Keywords:flame structure;fractal geometry;combustion theory 大量研究表明,在实际火花点火发动机中的预混湍流燃烧均属于高度折皱的层流火焰[1].但长期以来,对高度折皱层流火焰的结构一直缺乏定量的描述和研究.1987年,F.C.Gouldin首先将分形的概念引入预流湍流燃烧研究之中[2].本文的目的是给出求取容弹内预混湍流火焰分数维数的具体方法,以及分数维数随u′/S L(u′和S L分别为湍流强度和层流燃烧速度)的变化规律,并导出相应的结论.1 从火焰照片求取分数维数的方法根据分形几何学理论,对各向均匀的分形,只需求出相应曲线的分数维数D2,即可用下式求得曲面的分数维数D3D3=D2+1(1)此外,有收稿日期:1998-03-02.作者简介:蒋德明,男,1934年4月生,能源与动力工程学院汽车工程系,教授,博士生导师. 基金项目:国家“八五”攀登计划资助项目(85-25-2-1).-D 2=lg N (ε)lg ε(2)1-D 2=lg L (ε)lg ε(3)1.1 直接法用不同的测量尺度ε沿曲线逐段测量,即可得到不同测量尺度时的总测量次数N (ε)以及曲线的总长度L (ε),L (ε)=ε·N (ε),然后,将测量结果以对数形式绘成坐标,即可得到其分数维数及内、外界.由于N (ε)或L (ε)是直接测得的(相对下面而言),故作者称之为“直接法”.1.2 面积法画一个半径为ε的圆,该圆中心位于曲线上;然后,沿曲线移动该圆(该圆的中心始终在曲线上),即形成宽度为2ε的“带子”,如图1所示;通过测量得到这一“带子”的面积A (ε)后,由A (ε)/2ε,即可求得相应的曲线长度(带长)L ;改变测量尺度ε,并重复上述过程,即可得到不同ε时的曲线长度L (ε),这样得到L (ε)与ε的对数坐标图,从而根据对数坐标图上的直线段,即可求得曲线的分数维数和相应的内、外界.由于曲线长度是通过测量“带子”的面积而得到的,故称之为“面积法”.图1 求取曲线分数维数的面积法示意图对火焰前锋外廓线,一般按上述方法用计算机进行处理.火焰前锋外廓线在计算机中是以离散的像素点来表示的,在运用面积法时,将所有到火焰外廓线(由一系列离散的像素点组成)的距离小于或等于ε的像素点累加起来,而每个像素点都对应于一定的面积,由此即可求得“带子”的面积,进而可求得其分数维数与内、外界.1.3 计算结果运用直接法和面积法,对本实验拍摄的火焰纹影胶片进行了计算处理,所得结果分别如图2和图3所示,图中坐标均为对数值.由图可知,这2个图中分别存在相应的直线段,其斜率所对应的分数维数D 2分别为1.095和1.101(相应地,其D 3值分别为2.095和2.101),这表明该湍流火焰具有分形特征,其分数维数D 2取为1.10.由图2可以看出,用直接法得不到明显的内、外界值εi 和εo ,而用面积法能够得到相应的εi 、εo 值.由图3可知,其外界εo 大约为8.3mm (100.92=8.3mm ),而内界εi 为1.7mm (100.24=1.7mm ).应该指出,由于图形空间分辨率的限制(大约在0.5mm 左右),由图2所得到的内界是不准确的.采用直接法求取分数维数时,在实际计算过程中只考虑了彼此连结的像素点的情况,因而所得到的曲线长度数值有较明显的波动,所以得不到明显的内、外界值;而采用面积法时,将所有不相连接的像素点都考虑进去,因而所得曲线长度值相当准确,且得到的分数维数比用直接法所得到的值略大一点,这与R .D .M atthew s 等人的结果是一致的[3].图2 用“直接法”处理火焰外廓线的结果图3 用“面积法”处理火焰外廓线的结果2 试验结果及分析图4为点火时刻在不同的湍流条件下,火焰半径R 大约为18m m 时容弹中的湍流火焰纹影照23第2期 蒋德明,等:预混湍流火焰结构的分形特征 (a )u ′i =0.2m /s ;L i =3.0mm (b )u ′i =0.8m /s ;L i =3.2mm (c )u ′i =1.5m /s ;L i =4.0mm 甲烷-空气; =0.9; p i =101.3kPa ; T i =300K图4 点火时刻不同湍流强度时的火焰纹影照片片.由图可以看出,点火时刻的湍流强度增大,则火焰前锋折皱程度加大,火焰结构也变得更为精细.图5为容弹内湍流火焰的分数维数D 3值随u ′/S L 的变化曲线.随着u ′/S L 的增大,在相同火 图5 容弹中湍流火焰的分数维数D 3值随u ′/S L 的变化曲线焰半径R 处,火焰前锋面的分数维数也增大.这表明其火焰前锋面的折皱、扭曲(粗糙)程度加大,从而增加了火焰前锋面的面积,因此燃烧反应在更大的火焰前锋面内进行,其瞬时质量燃烧率也随之增加,火焰燃烧速度加快.这一结果揭示出,在中、小程度的湍流强度条件下,湍流燃烧速度随湍流强度增大而增加的真正原因. 表1为4个实验工况下对应不同火焰半径处火焰前锋面的分数维数D 3值.由此可知,在同一实验工况下,当火焰半径增大时,其分数维数D 3是增加的.这表明,在容弹内预混湍流燃烧过程中,随着火焰半径增加,火焰前锋的折皱程度增加,因此其火焰燃烧速度也相应地增大.同时也进一步证实,本研究容弹内的预混湍流燃烧属非充分发展的湍流火焰,湍流对燃烧影响(促进作用)是一个渐进的过程,随着火焰的不断发展,这一促进作用愈加明显和强烈.表1 4个实验工况下的分数维数D 3值点火时刻的湍流参数u ′i =1.5m /s L i =4.0mmu ′i =0.8m /s L i =3.2mmu ′i =1.0m /s L i =3.5mmu ′i =0.7m /s L i =3.0mm当量燃空比0.90.90.70.7D 3(R =13mm )2.082.062.092.07D 3(R =25mm )2.122.092.132.093 结 论(1)火花点火发动机内的预混湍流火焰具有分形特征.(2)在中、小程度的湍流强度下,分数维数随湍流强度的增加而增加,湍流燃烧速度增大.(3)火花点火发动机中的预混湍流火焰属于非充分发展湍流火焰.参考文献:[1] 蒋德明.火花点火发动机的燃烧.西安:西安交通大学出版社,1992.[2] Gouldin F C .A n application of fractals to modelling pre -mix ed turbulent flame .Combustion and F lame ,1987,68:249～266.[3] Hall M J ,M atthew s R D .Fractal analysis of turbulentpremixed flame images from SI engines .SAE P aper 922242,1992.(编辑 蒋慧姝)24西 安 交 通 大 学 学 报 第33卷。
预混燃烧的燃烧模型
预混燃烧的燃烧模型摘要为了达到抑制污染物排放，实现燃料的清洁燃烧的目的，人们采取了很多办法。
“节能减排”促使燃烧系统采用贫燃燃烧技术，它具有降低NOx、CO等污染物，提高燃烧效率的作用。
但这种燃烧方式的燃烧极限范围很窄，而且火焰稳定性差，容易诱发燃烧系统的不稳定性，如火焰的热声耦合振荡，这种不稳定性会造成更大的污染和浪费。
新型燃烧器的设计必须克服这些缺点，以达到“节能减排”的目的。
首先本文以FLUENT软件为平台，构建了合理的数学物理模型，对甲烷-空气预混燃烧过程进行了数值模拟，实验证明，贫燃料燃烧及贫氧燃烧都可以起到降低污染物排放的目的。
并利用数值模拟的方法针对不同燃烧模型的情况下甲烷的预混燃烧的特性进行分析，观察其NO）的分布情况，发现预混燃烧的相关规律，寻求燃烧速度场、温度场、以及污染物（X的最佳工况。
其次本文了解不同燃烧模型对流场结构、燃烧结构的影响，与实验结果比较，探讨如何改进数值模拟，提高设计精度，同时找出预混火焰稳定性规律，探讨抑制燃烧不稳定性的策略。
本文通过数值计算，得到了在不同燃烧模型下柱状燃烧室内甲烷燃烧的数值模拟结果，分析发现，燃烧模型的不同对甲烷燃烧特性的影响也不同。
通过对燃烧速度分布图，火焰温度分布云图，燃烧的污染物NO的云图进行分析研究，得出结论。
关键词预混燃烧数值模拟FLUENT 部分预混燃烧Title Pre-mixing combustion combustion modelAbstractIn order to achieve inhibit pollutants, realize fuel clean burning purpose, people taken a lot of measures. "Energy conservation and emission reduction" prompted combustion system using poor fuel combustion technology, it has to reduce pollutants such as NOx, CO, increase the combustion efficiency role. But this kind of combustion way combustion limit range is very narrow, and flame stability is poor, and likely to cause combustion system instability, such as flame of thermoacoustic oscillation, the coupling instability will cause more pollution and waste. New burner's design must overcome these shortcomings, to achieve "the purpose of energy saving and emission reduction".Firstly this paper with FLUENT software for the platform, and constructs the reasonable mathematical physics model of methane - air pre-mixing combustion process was simulated, the experiment proof, the poor fuel combustion and poor oxygen burning can reducing pollutant purpose. And using the method of numerical simulation of combustion model for different under the condition of pre-mixing combustion characteristics of methane areanalyzed, observe its velocity field and temperature field, and the distribution of pollutants (), found the relevant law pre-mixing combustion, seeking the best condition burning. Then this paper to understand different combustion model convection field structure, the influence of combustion structure, compared with the experimental results, this paper discusses how to improve the design accuracy numerical simulation, and at the same time, improve the stability pre-mixed flame out rules and explore the inhibiting combustion instability strategy.This article through numerical calculation, obtained in different combustion model columnar combustion chamber under the numerical simulation results of methane combustion, analysis, we found that the different combustion model for the influence of methane combustion characteristic of different also. Through the burning rate distribution, the flame temperature distribution of convective, the combustion pollutants analysis of NO cloud, draws the conclusion.Keywords：Pre-mixing combustion Numerical simulation FLUENT Part pre-mixing combustion绪论课题的研究背景及意义燃烧室作为燃气轮机中最重要的部件，是利用燃料的燃烧，提高进入涡轮的气流温度的装置。
3.2火焰面模型
1
2
3
混合物分数及组分、能量方程
Z Z Z u1 ( D ) t x1 x1 x1
（3-87） （3-88） （3-89）
Y Y Y u1 ( D ) W t x1 x1 x1
T T T C p u1C p ( ) h W t x1 x1 x1 1
G u G S L | G | t
对皱褶的层流火焰面和波纹板式的燃烧机制
G 0 u G S L | G | DL K | G | LS | G | t
0 L
其中 S 是不受拉伸的平面火焰的传播速度，K是火焰面的 曲率，S是流场应变率，L和层流火焰面厚度成正比，且 数量级相同，D 是一个扩散系数。
（3-101） （3-102）
Z 2 Z 2 Z Z 2 2 u ( Z 2u 2 D( ) j ) 2 Z u j t x j x j x j x j
梯度输运模型
Z u j Dt Z x j ， Z 2u j Dt Z 2 x j
（3-91） （3-92）
标量的瞬时耗散率
Pitsch等提出
2 D(
Z 2 ) x1
（3-93）
（3-94）
st f (Z )
Z 2 ln Z ， f (Z ) 2 Z st ln Z st
在湍流燃烧中，标量耗散率的平均值
Z Z st
P( )d f ( Z ) P( Z )dZ
f (G; xi , t ) ，
G( xi , t )
G2 ( xi , t )
湍流预混火焰模型和几率密度函数的输运方程
应在这两种气团的交界面上发生，认为平均化学反应
率决定于末燃气团在湍流作用下破碎成更小气团的速
率，而破碎速率与湍流脉动动能的衰变速率成正比。
R fu ~ / k
湍流燃烧速率
对比用k - ε模型和混合长度模型计算湍流粘度的公式
t C k 1/ 2 C C D k 2 /
得
R fu ,T 0.4m fu
u y
对比EBU模型的燃烧速率表达式 R fu ,T
u cEBU m fu y
易知EBU模型是拉切滑模型在一定条件下的简化形式
进一步分析公式
R fu ,T
u y u ln 1 u / S y (m fu ,u m fu ,b ) u
可见在不均匀性很强的流场中，湍流燃烧速率主要
取决于流体应变率。
在较均匀的流场中（
），公式可简化为
R fu ,T (m fu m fu ,b ) S / u
可见此时湍流燃烧速率受层流火焰传播速度的影响
较大，故合理地估算层流火焰传播速度是正确运用拉 切滑模型的关键之一 。
层流火焰传播速度是可燃气的物理化学性质，它取决于 混合物的热力学状态(如压力和温度)，对温度尤为敏感。
模化方法
是概率分布函数和脉动速度的二阶关联项，按 照“梯度准则”进行模拟，在物理上表示几率分布函 数由于湍流而引起的输运特性。
所以：
P ( )ui c3 P ( ) xi k
2
模化方法
2 因子 ( / xi ) 主要与微尺寸的小脉动有关，而因子
P(φ )主要受大尺度的大脉动控制，不妨认为这两者不
下标u和b分别表示末燃状态和已燃状态
第八章 预混湍流火焰
图 8.4 用 Rayleigh（瑞利）散射方法测得的瞬时气体密度
上述数据表明，湍流预混火焰传播可以描述为在湍流流动中传播的薄的层流预混火焰。 湍流流动使火焰发生了扭曲，扭曲的程度则取决于当地的湍流程度。这一观察表明，除了
火焰厚度 δ f 以外，还存在另一与湍流速度脉动有关的长度尺度和速度尺度。 湍流预混火焰速度可以用流经火焰的可燃预混气的体积流量 Q& 与湍流火焰的表观面积
142
长度缩短，并且在直接长时间曝光照片中，可见到火焰区域扩大了，出现了皱褶；在火花 纹影照片中，我们可以看到火焰前沿出现了不规则的皱褶。图 8.2 是一系列湍流火焰边界纹 影照片的叠加。
观察到的利用钝体或值班火焰稳定的预混火焰具有相似的结构特点。 另一种方法是在未燃混气中加入小粒子利用米氏散射来显示预混湍流火焰结构。在这 种方法中，用激光片来照射流动区域，火焰前锋位于亮和暗的分界线上，亮的地方是激光 片被粒子散射的结果，而暗的地方则是由于燃烧放热，气体膨胀，从而使粒子的密度减小， 减少了对激光的散射所致。图 8.3 是在活塞发动机中、不同的活塞运动速度下，用该法拍摄 的一些预混湍流火焰传播的图象。我们再次看到，火焰前锋出现了褶皱，并且褶皱的程度 与活塞运动的速度有关。
于 100％（在高的剪切、低的均速流中）。
145
湍流强度在描述湍流流动特性时有用，但是它没有反映湍流流动中参数的脉动尺度或 频率。湍流流动通常有一个宽广的脉动频率（或者波长，振幅）范围，包括从高频（小尺 度）的能量耗散涡流微团到低频（大尺度）的涡流微团，这些低频微团宏观输运着质量、 动量和能量。通常，用单一的尺度是无法全面描述湍流流动的所有特征的，对尺度参数的 选择取决于我们感兴趣的那种湍流现象。在湍流燃烧中，影响湍流火焰传播的因素很多， 我们可以考虑定义湍流的尺度的几种方法。
6-湍流预混火焰讲解
大尺度弱湍流传播速度确定——小 火焰模型（表面理论）
设薄层焰锋的传播速度仍然是un，那么单位时间内焰锋锋 面烧掉的混合气是Acun，它应与湍流火焰传播速度ut和湍 流焰锋的平均面积Ap的乘积相等，即：
Acun=Aput 或 ut=Acun/Ap 因为Ac>Ap,故ut>un, 若把湍流气团设想成凹凸不平的很多 小的焰锋，则ut>un, 等于这些小的椎体表面积和底面积之 比。 ——小火焰模型，亦称湍流火焰传播的表面理论
燃烧学
6-湍流预混火焰
湍流预混火焰传播 湍流预混火焰传播图域 湍流预混火焰传播速度确定 湍流火焰传播速度影响因素
第一节 湍流预混火焰传播
研究湍流火焰的目的
（1）工程中的燃烧装置多为湍流燃烧 （2）确定湍流特性对火焰传播的影响 雷诺数Re=ρvL/μ 直管段中： Re<2300时，层流； Re>3200时，湍流 此时火焰为湍流火焰
（5）混气浓度 化学恰当比或偏富时速度最大
St
m
Au
湍流预混火焰传播速度要比层流预混火焰传 播速度快
湍流火焰比层流火焰传播快的原因
（1）湍流流动使火焰变形，火焰表面积增加，因而增大了 反应区； （2）湍流加速了热量和活性中间产物的传输，使反应速率 增加，即燃烧速率增加； （3）湍流加快了新鲜混合气和燃气之间的混合，缩短了混 合时间，提高了燃烧速度。 湍流火焰理论基于上述概念发展起来的。 湍流火焰传播理论主要有两种： （1）表面褶皱理论（邓克尔和谢尔金） （2）容积燃烧理论（萨默菲尔德和谢京科夫）
第二节 湍流预混火焰传播图域
湍流预混火焰的性质既依赖于预混层流火焰的特性（如SL和
第五章 湍流燃烧模型
grad mox
gradvj mfu mox
K
f
(mox
Sm fu
)mfu mox
Smox mf2u
m fu
mo2x
(mo x
Smfu
)mfu
mo x
式中D1表层流交换系数；S表示化学当量比；
（131）
K f B exp(E / RT )
方程（131）中第三、五、六项需进行模化，才能使其封闭。 其方法如下：
与能量的湍流扩散和颗粒的经历效应 Euler坐标系中处理气相；Lagrange坐标系中描述颗粒相 有反应颗粒相的连续介质 -- 轨道模型和考虑颗粒经历效应的多
流体模型
连续介质―轨道模型的基本方法
用多流体模型求解Euler坐标系中颗粒相的连续与动量方程， 求出颗粒速度与浓度分布，同时沿着Euler坐标系中计算得到的 轨道或流线追踪因反应和传热引起的颗粒质量和温度的变化， 使用常微分方程和代数式
一般地，组分分布需要通过求解两个以平均化学反应率为源项 和耦合的二阶非线性偏微分方程
在简单化学反应系统的假设下，通过引入如下定义的质量分
数 f 简化
f m fu mox / S
（107）
式中（m fu 及 mox 分别为燃料及氧化剂的质量分数的时均值）
就可以将确定组分质量分布转变为只需求解一个有源方程和一
如果 f 0 ，则 m fu mox 0 。
时均值 f 及其脉动均方值 g f 2 的输送方程形式分别为
(f )
t
xj
(
vj
f
)
xj
( f
f )
xj
（109）
t
(g)
xj
( vjg)
第五章1,燃烧模型
预热区 O(1) CH4
氧化层 O(ε)
T
T0
103
sL
内层 O(δ ) O2 CO
涡旋尺度小于火焰内层尺度
破碎反应区
Ka 1
H2
x / lF
10 2
l
/ sL
10
Re 1
涡旋尺度小于火焰尺度
薄反应区
F 涡旋尺度大小于火焰尺度
l
Ka 1
1
波纹板式火焰
涡旋速度与火焰传播速度比较
皱褶层流火焰面燃烧机制
AT
m
A
sL
sT
m
在横截面积为A的湍流火焰中，瞬时层流火焰面的形状 是不断变化的。根据质量守衡，通过该湍流火焰面AL 的流量和某一时刻通过瞬时皱褶层流面面积的流量是 相等的 m u SL AL u ST A
湍流火焰传播速度（皱褶）
• 湍流火焰传播速度ST和层流火焰传播速度SL之 ST AL 比为
L
湍流火焰传播速度（皱褶）
在薄反应区燃烧机制中，湍流涡旋深入 预混区内部，使扩散速度增加。我们用 湍流扩散系数 t u0l 代替层流扩散系 数 D ，即得湍流预混火焰的传播速度
( D W )1 2
SL (DW )1 2
u0 l 1 2 ST ( ) SL S L
0
尺度涡旋不能进 入火焰的内层。
10 4
层流火焰
0.1
皱褶层流火焰
10 10 2 103
0.1
1
l / lF
波纹板式火焰面燃烧机制和薄 反应区机制之间的关系
在大Re数条件下，小尺度湍流是间歇性的，耗 散率满足一定的统计分布。 在对湍流的Kolmogorov理论进行修正时，提出 了间歇性或流场中湍流的作用是不规则的这样 一些概念。反映到燃烧问题中，大Re数条件下， 湍流火焰可以有被小尺度涡旋产生的很强的局 部混合区（如在薄反应区机制中），也可以有 光滑的火焰面区。这时就产生了波纹板式的火 焰面燃烧机制。 上述讨论的两种机制可以在一个湍流火焰中同 时存在。
层流火焰与湍流火焰结构的比较
层流火焰与湍流火焰的结构一、层流火焰分析层流火焰速度和火焰结构一维层流火焰在预混燃料-氧化剂混合物中传播是最简单的燃烧现象之一，在此火焰中，化学动力学以及能量和组分扩散输运起重要作用。
通过守恒方程和状态方程可以导出Rankine-Hugoniot 曲线。
该曲线把在一维层流预混火焰中未燃气和已燃气状态联系起来。
已燃气体位于Rankine-Hugoniot 曲线下分支（缓燃），并相应于未燃气体状态Rayleigh 线与具有适当反应热的Rankine-Hugoniot 曲线交点L，如图1-2中所示。
图1 层流预混火焰坐标系图 2 一维燃烧波的Rankine-Hugoniot 曲线和Rayleigh 线Rayleigh 线的斜率与相对于未燃气体的波的传播速度，即层流火焰速度有关。
/ ( / )2 ( )2 u u dP dv = −m& A = −ρu(= ) = u u u S 层流火焰速度= (1/ ) (dP / dv) u −ρ由于缓燃Rayleigh 线斜率比爆震Rayleigh 线斜率小得多，所以缓燃速度比爆震速度小得多。
虽然守恒方程和状态方程提供了缓燃的未燃气体和已燃状态之间的关系，但不能唯一确定层流火焰速度Su 。
为了确定Su ，必须将守恒方程通过缓燃波积分。
由于 96方程是非线性耦合微分方程，其准确解只有通过数值积分才能获得。
它需要很大的计算资源。
为了考察层流火焰的某些特征（如火焰速度和厚度）以及这些特征与燃烧参数如燃料类型、化学配比、压力及未燃气体的温度的关系，对方程组进行了简化，以便能分析求解。
要得到简化的模型，需要引入一系列的假设。
我们从考察参考系建立在火焰上的层流火焰结构的某些方面入手。
如前所述，这些计算是针对等压过程进行的。
但是对一维缓燃的Rankine-Hugoniot 曲线，如图2 所示，已燃气的压力小于未燃气的压力。
现在我们需要考察压力减少的数值是否小到可以忽略的程度。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
单变量概率分布函数输运方程
D Dt
P( )
xi
P( )ui
P( )S ( )
P( )
xi
2
(3-76)
S(φ)是变量φ的源或汇
为使方程封闭，必须对有关的项进行模化。
模化方法
概率分布函数和脉动速度的二阶关联项
按照“梯度准则”进行模拟，在物理上表示概率分 布函数的湍流输运
或把式(3-81)写成
Rfu B 2mfumox exp(E / RT )[1 F]
(3-82)
F 概括了湍流脉动对平均化学反应率的影响，是
对燃烧速率进行模化的困难所在。
对燃料和氧化剂质量分数脉动值 的二阶关联项 mfumox 的控制方程
Borghi等人为简化模拟过程，略去温度脉动的影响，提
出了在F中的 mfumox 的控制方程
xi
ui mfumox
1
Dl grad mfu mox 2
2Dl gradmfu 3
gradmox
2t f
gradm fu 4
gradmox 5
graduimfu mox
K f [(mox Smfu )mfumox Smox mfu2 m fu mo2x (mox Smfu )6mfumox ]
评价
在简单的湍流火焰计算中获得与实验基本符合 的结果，仍需改进和完善。
建立双变量（混合分数和反应度）的联合概率 分布函数的输运方程(Pope)
§3.4 平均反应速率的 输运方程模型
湍流流动模型：模拟雷诺应力，建立了雷诺应 力的输运方程模型，在某些情况下获得了优于 应用湍流粘性系数模型得到的结果
(3-83)
式中Dl表示层流交换系数，S表示化学当量比， Kf = B exp(E/RT)
mfumox 控制方程中各项的意义
第一项
xi
ui mfumox
表示因变量的对流速率
第二项 Dlgradmfumox 表示该因变量的层流输运
Dl gradmfu gradmox 表示分子扩散引起的脉动耗散 第不四均项匀2性ft g引rad起m的fu g因ra变dm量ox 的表净示增由长浓率度均值的空间
2Dl gradmfu
gradmox
2
k
mfu mox
(3-84)
mfumox 控制方程中各项的模拟-2
第五项―湍流输运项的模拟：可以参照常用的
“梯度模拟准则”
xi
ui mfu mox
xi
ut
f
xi
mfu mox
(3-84)
第六项：可以用简单地略去脉动的三阶关联项 来模拟
模拟后 mfumox 控制方程
§3.3 概率分布函数 的输运方程模型
回顾
湍流脉动在湍流燃烧过程中起着重要作用 湍流脉动在某种程度上具有随机量的特征
在湍流燃烧的数学模型研究中引用概率分布
➢ k-ε-g 模型中的混合分数瞬时值的城墙式分布，相应 于双δ函数形的概率分布
➢ 截尾正态分布概率分布函数 ➢ Beta函数分布等形式的概率分布函数
第五项 graduimfumox 表示因变量的湍流输运 第六项(其余)表示化学反应引起的因变量的变化
mfumox 控制方程中各项的模拟-1
为求得方程的封闭，必须对方程的第三、五和 六项进行模化。
第三项：在湍流场中因变量的耗散速率可以考
虑成与该变量本身、密度和湍流耗散速率成正
比（适用于Schmit数接近于1的情况）
不同模型算得的气 体燃料锅炉内轴上 温度分布与实验数 据的对比。
用上述两个模型和 旋涡破碎模型算出 的结果大致相同， 然而本节的两个模 型却要多解好几个 微分方程，这个代 价没有换取到应有 的成功。
P( )ui
c3
k2
xi
P( )
(3-77)
模化方法
因子 ( / xi )2主要与微尺寸的小脉动有关，而 因子P(φ)主要受大尺度的大脉动控制，不妨认 为这两者不相互关联
2
P( )
xi
P( )
xi
P(
)c4
k
2
(3-78)
2 表示φ的脉动均方值
单变量PDF输运方程
D Dt
P( )
在湍流燃烧模型的研究中人为模拟平均化学反 应速率，是否可以建立平均反应率的输运方程 模型呢？
定压燃烧过程瞬时反应速率
设反应机理用单步不可逆反应来表征，瞬时反 应速率遵守双分子碰撞模型的Arrehnius公式
Rfu B 2mfumox exp(E / RT )
(3-80)
经过对因变量mfu，mox和T的雷诺分解以及对 公式的雷诺平均，并利用泰勒级救展开，略去
脉动值的三阶以上关联量，便可得到平均化学
反应速率的表达式
平均化学反应速率
R fu
B 2m fu mox
exp(E / RT )[1
a0
mfu mox m fu mox
Hale Waihona Puke a1T 2 T2
a2
T mox Tmox
T mfu Tm fu
]
(3-81)
式中，a0＝1，a1=1/2(E/R)2-(E/R)，a2= E/R。
xi
ui
mfu mox
xi
eff f
mfu mox xi
c1eff
m fu xi
mox xi
2
k
mfu mox
Bmfumox exp(E / RT )
(mox
sm
fu
)
mfu m fu
mox mox
mo2x mox
mfu2 m fu
(3-85)
为使该方程封闭，还需给出以 mfu2 和 mox2 为因变 量的微分方程
评价
物理思想和推理较清楚 由于需要模化的量的数目增加和完全略去了温
度脉动对燃烧速率的影响，没有能够达到提高 计算精度的目的。 在该模型的基础上进一步考虑温度脉动对燃烧 速率的影响，建立、模化并求解以 m0xT， mfuT 和T 2 为因变量的微分方程，但改进 甚微。
气体燃料锅炉内轴上温度分布
PDF→受输运方程控制的因变量
对燃烧现象的更深入的实验研究发现：
➢ 在不同的燃烧过程中，因变量的概率分布形式不同； ➢ 即使对于同一个燃烧过程，在不同的区域，各个因
变量的概率分布函数也不尽相同。
启示：概率分布函数是否也是一个受输运方程 控制的因变量，与火焰特性、流动图象、湍流 特性及边界条件有关。
xi
c3
k2
xi
P( )
P(
)S
(
)
c4
k
2
P(
)
(3-79)
空间位置和变量φ的函数 辅之以适当的定解条件
二维管道火焰稳定器后面的湍 流预混火焰
Khalil 不同轴向位置的横截
面上轴向速度
计算结果与实验符合 得较好
用RBU模型和PDF输 运方程进行计算的结 果大抵相同，但是后 者付出的计算时间和 计算机的贮存量的代 价比前者大得多