FLUENT燃烧简介FLUENT软件中包含多种燃烧模型、辐射模型及与燃烧相关的湍流模型，适用于各种复杂情况下的燃烧问题，包括固体火箭发动机和液体火箭发动机中的燃烧过程、燃气轮机中的燃烧室、民用锅炉、工业熔炉及加热器等。

1.1 FLUENT燃烧模拟方法概要燃烧模型是FLUENT软件优于其它CFD软件的最主要的特征之一。

FLUENT可以模拟宽广范围内的燃烧问题。

然而，需要注意的是：你必须保证你所使用的物理模型要适合你所研究的问题。

FLUENT在模拟燃烧中的应用可如下图所示：图 1 FLUENT模拟过程中所需的物理模型1.1.1 气相燃烧模型一般的有限速率形式(Magnussen模型)守恒标量的PDF模型(单或二组分混合分数)层流火焰面模型(Laminar flamelet model)Zimount 模型1.1.2 离散相模型煤燃烧与喷雾燃烧1.1.3 热辐射模型DTRM，P-1，Rosseland 和Discrete Ordinates 模型1.1.4 污染物模型NOx模型，烟(Smoot)模型2.1气相燃烧模型·在FLUENT中，针对不同的燃烧现象，采用了不同的化学动力学处理手段，以减少计算成本，如下：有限速率燃烧模型---预混、部分预混和扩散燃烧混合分数方法(平衡化学的PDF模型和非平衡化学的层流火焰面模型)---扩散燃烧反应进度方法(Zimont模型)---预混燃烧混合物分数和反应进度方法的结合---部分预混燃烧2.2.1 有限速率模型化学反应过程一般采用总包机理(即简化化学反应，如单步反应)进行描述。

求解积分的输运方程，得到每种组分的时均质量分数值，如下：-----(1)其中组分j的反应源项为所有反应K个反应中，组分j的净生成速率：-----(2)-----(3)计算所需参数包括：1、组分及其热力学参数值；2、反应及其速率常数值。

有限速率模型的有缺点：优点：适用于预混、部分预混和扩散燃烧，简单直观；缺点：当混合时间尺度和反应时间尺度相当时缺乏真实性，难以解决化学反应与湍流的耦合问题，难以预测反应的中间组分，模型常数具有不确定性。

这种模型求解反应物和生成物输运组分方程，并由用户来定义化学反应机理。

反应率作为源项在组分输运方程中通过阿累纽斯方程或涡耗散模型。

应用领域：该模型可以模拟大多数气相燃烧问题，在航空航天领域的燃烧计算中有广泛的应用。

2.2.2守恒标量的PDF模型守恒标量的PDF模型仅适用于扩散(非预混)燃烧问题，该方法假定了反应是受混合速率所控制，即反应已经达到化学平衡状态，每个单元内的组分及其性质是由燃料和氧化剂的湍流混合强度所控制，其中涉及的化学反应体系由化学平衡计算来处理(利用FLUENT的组件程序PrePDF)。

该方法通过求解混合物分数及其方差的输运方程获得组分和温度场，而不是直接求解组分和能量的输运方程。

-----(4)-----(5)其中-----(6)混合分数定义-----(7)其中Zk代表元素k的元素质量分数，下标F和O分别代表燃料和氧化剂的进口值。

对于简单的燃料/氧化剂体系，每一计算单元内的混合物分数代表了该单元内的燃料质量分数，由于混合物分数是守恒标量，因此在求解输运方程时不再考虑反应源项。

在该方法中，化学反应认为足够快，体系中的组分立刻达到平衡状态。

化学平衡组分在混合物空间的分布可示意如下：化学反应和湍流之间的相互作用采用概率密度函数(PDF)的方法处理：上图代表了概率密度函数p(V)的定义，因此在混合物分数空间，f标量的时均值可由下式计算：守恒标量PDF模型的优缺点：优点：可以预测中间组分的浓度，可以考虑流动中的耗散现象，可以考虑化学反应与湍流之间的相互作用，该方法不需要求解大量的组分和能量的输运方程，因此可以缩短计算时间。

缺点：研究的流动体系必须接近于局部化学平衡状态，且不能用于非湍流流动，同时亦不能处理预混燃烧问题。

该模型不要求用户显式地定义反应机理，而是通过火焰面方法(即混即燃模型)或化学平衡计算来处理，因此比有限速率模型有更多的优势。

应用领域：该模型应用于非预混燃烧（湍流扩散火焰），可以用来计算航空发动机的环形燃烧室中的燃烧问题及液体/固体火箭发动机中的复杂燃烧问题。

2.2.3 层流火焰面模型层流火焰面模型的基本思想是把湍流扩散火焰看作是层流对撞扩散火焰面的系统。

该方法可以看作是守恒标量PDF模型的一个扩展，它可用于处理非化学平衡状态的体系，即可以利用化学反应动力学的方法处理反应流。

不同于守恒标量的PDF模型，标量是混合物分数和标量耗散率的函数，而非混合物分数的函数：-----(8)指定混合物分数f的PDF符合β函数分布，标量耗散率χ的PDF符合狄拉克-δ函数分布，因此，时均标量值可以通过在f和χ空间求标量的统计平均来得到(即，考虑化学反应与湍流的相互作用)：-----(9)层流火焰面模型的计算过程如下：1、计算不同标量耗散率下，标量在混合物分数f空间的分布，即求解火焰面方程，组分方程为：-----(10)能量方程为：-----(11)从而得到标量随混合物分数和标量耗散率的变化关系，即式(8)，并以火焰面数据库文件的形式保存结果。

2、火焰面数据库文件也可由其它软件生成，若得到的库文件为单标量耗散率，则需计算不同标量耗散库的库文件，最后将它们合并。

3、利用(9)式计算火焰面的PDF库，从而得到时均标量随平均混合物分数和平均混合物方差的变化关系。

4、利用求解平均混合物分数及其平均方差的输运方程的方法，在流场中计算这两个量，然后再利用得到的PDF库查找时均标量值。

2.2.4 预混燃烧的Zimont模型湍流预混燃烧的化学反应采用反应进度(progress variable),c,进行表征，如下：-----(12)其中Yp和Yp(ad)分别代表当前和完全绝热燃烧后燃烧产物的质量分数，其取值范围在0到1之间，0代表未燃混合物，1代表已燃混合物。

若用反应进度c代表其平均值，则其输运方程可表达如下：-----(13)上式中平均反应速率项如下求解：-----(14)ρv代表未燃物密度，Ut代表湍流火焰传播速度。

湍流预混燃烧的关键在于求解湍流火焰传播速度(位于湍流火焰表面的法线方向)，该速度受两方面馆因素的影响：一是层流火焰传播速度，即决定于燃料和氧化剂的浓度、初始温度、组分的扩散特性以及化学反应动力学特性；二是有大涡褶皱和拉伸以及由小涡决定的火焰表面厚度。

根据上述讨论，FLUENTZ中的湍流火焰传播速度可表达为：-----(15)式中，A模型常数，u`速度均方值，Ul层流火焰传播速度，α=k/ρCp未燃物的分子导热系数，I=CdU3/ε湍流长度尺度，τt为湍流时间尺度，τc为化学反应时间尺度。

为考虑火焰拉伸所导致的吹熄现象，在反应速率源项中可乘以一个拉伸因子G，它代表了拉伸所导致火焰不熄火的概率。

-----(16)2.2.5 部分预混燃烧模型部分预混燃烧系统是指这样一种预混火焰，其燃料/氧化剂之比不唯一。

FLUENT中的部分预混模型是非预混模型和预混模型的结合。

预混燃烧的反应进度，c，决定了火焰前锋的位置，在火焰前锋的后面(c=1),混合物已燃，使用守恒标量PDF或层流火焰面模型的解；在火焰前锋的前面(c=0),组分质量分数，温度和密度由混合但未燃烧混合物分数来计算。

在火焰内部，未燃和已燃混合物的线性结合的方法被使用。

部分预混模型求解平均反应进度c，平均混合物分数f和混合物分数方差f2的输运方程。

平均标量可由如下的f和c的PDF来计算：-----(17)在火焰很薄的假设下，由于存在未燃的反应物和已燃的产物，则平均标量可如下计算：-----(18) 具体煤粉燃烧算例1、建立求解模型：连续相(气体)只有在非耦合求解时非预混燃烧模型才是可用的。

Define---models---solver...2、打开RNG k-e湍流模型Define---models---viscous...3、打开非预混燃烧模型Define---models---species...a、在Model下选择Non-Premixed Combustion。

当点OK时，FLUENT将打开一个对话框，要求输入在模拟中要用到的PDF 文件。

b、在选择文件对话框中，选择并读入非绝热的PDF文件(coal.pdf)当FLUENT读入非绝热PDF文件时，它会自动激活能量求解方程，所以你可以不用打开能量面板激活传热方程。

3、选择P-1以激活辐射模型。

Define---Models---Radiation....P1模型是能求解气体和颗粒间辐射传热的模型之一。

2、建立求解模型：离散相FLIENT会用离散相模型来模拟煤粉的流动。

这模型会预示出单个煤粒的轨迹，每一个都代表煤的连续流，在交替计算离散相的轨迹和气相连续方程时，煤粒与气相间热量、动量、质量的传递都将包含其中。

1、耦合离散相与连续相流动预报。

Define---Models---Discrete Phase.....a、在Interaction 下，选中Interaction with Continuous Phase 选项。

这选项将激活耦合求解，在求解中，离散相的轨迹将会对气相产生影响，如果不选中这选项，你仍可以看到煤粒的轨迹，但上述参数对连续相的流动将没有任何影响。

b、定义耦合参数，设定Number of continous phase interations per DPM interation 为200；在一些有着高质量粒子和较大网格尺寸的问题中，应该给这参数设定高一点的值，这对低频率轨迹是很有好处的，为了更完全地聚合气相方程，应先对轨迹进行反复的计算。

c、在Tracking paramete 下，为MAX Number of steps 输入500000。

d、打开Specify Length Scale ，保持Lenth Scale的默认值为0.01m。

Length Scale 控制离散相轨迹综合中用到的每一次步数的大小。

这儿用到的值0.01m意味着10m长的一段轨迹要计算1000步左右。

e、在Options下，选择Particle Radiation Interation。

2、创建离散相煤的射入轨迹。

煤粉流用初始条件定义，这初始条件认为煤是进入到气体中。

在颗粒的运动方程的每一次综合中，FLUENT将用这些初始条件作为计算的开始点。

在这里，煤的总质量流的比率为2.4653kg/s，在10-160微米直径方向上，假设颗粒是服从Rosin-Rammler尺寸分布的。

其他的初始条件以及适当的输入程序将在下面做详细说明。

Define----Injections.....a、在Injections面板中点击Creat按钮。