第六章FLUENT中的燃烧模拟.docx
《第六章FLUENT中的燃烧模拟.docx》由会员分享,可在线阅读,更多相关《第六章FLUENT中的燃烧模拟.docx(21页珍藏版)》请在冰豆网上搜索。
![第六章FLUENT中的燃烧模拟.docx](https://file1.bdocx.com/fileroot1/2022-12/13/eefe668c-05d0-42c7-be9a-1854cac7c3ec/eefe668c-05d0-42c7-be9a-1854cac7c3ec1.gif)
第六章FLUENT中的燃烧模拟
第六章,FLUENT中的燃烧模拟
6.1燃烧模拟的重要性
面向实际装置(如锅炉、内燃机、火箭发动机、火灾等)面向实际现象(如点火、熄火、燃烧污染物生成等)
6.2FLUENT燃烧模拟方法概要
FLUENT可以模拟宽广范围内的燃烧(反应流)问题。
然而,需要注意的是:
你必须
保证你所使用的物理模型要适合你所研究的问题。
FLUENT在燃烧模拟中的应用可如
下图所示:
气相燃烧模型
一般的有限速率形式(Magnussen模型)守恒标量的PDF模型(单或二组分混合物分数)
层流火焰面模型(Laminarflameletmodel)
Zimont模型
离散相模型
煤燃烧与喷雾燃烧
热辐射模型
DTRM,P-1,Rosseland和DiscreteOrdinates模型
污染物模型
NOx模型,烟(Soot)模型
6.3气相燃烧模型
6.3.1燃烧的化学动力学模拟
实际中的燃烧过程是湍流和化学反应相互作用的结果,燃烧的化学反应速率是强非线性
和强刚性的。
通常的化学反应机理包含了几十种组分和几百个基元反应,而且这些组分之间
—92
的反应时间尺度相差很大(10〜102秒),因此在实际问题的求解过程中计算量和存储量极
大,目前应用尚不现实。
在FLUENT中,针对不同的燃烧现象,采用了不同的化学动力学处理手段,以减少计算成本,如下:
有限速率燃烧模型一一>预混、部分预混和扩散燃烧
混合物分数方法(平衡化学的PDF模型和非平衡化学的层流火焰面模型)一一>扩散燃
烧
反应进度方法(Zimont模型)>预混燃烧
混合物分数和反应进度方法的结合一一>部分预混燃烧
6.3.2一般的有限速率模型
化学反应过程一般采用总包机理(即简化化学反应,如单步反应)进行描述求解组分的输运方程,得到每种组分的时均质量分数值,如下:
細)+0(刚)一"+皿+禺
其中组分j的反应源项为所有K个反应中,组分j的净生成速率:
Rj=Rjk
k
式中,反应k中的组分j的反应速率可按照Arrhenius公式、混合(mixing)速率或“eddybreakup"速率的方法求解。
在混合(mixing)速率方法中,混合速率和涡的时间尺度,k/.有关,其物理意义为化学反应受限于湍流导致的组分和热量的混合速率。
Ji表达如下:
宀十%+瓷冲
计算所需参数包括:
(i)组分及其热力学参数值;(ii)反应及其速率常数值。
其中,FLUENT提供了一个混合物组分的数据库可供查找选用,另外也提供了一个化学反应机理以及组分热力学性质的数据库可供查找选用。
有限速率模型的优缺点:
优点:
适用于预混、部分预混和扩散燃烧;简单直观
缺点:
当混合时间尺度和反应时间尺度相当时(即Da>>1)缺乏真实性;难以解决化学
反应与湍流的耦合问题;难以预测反应的中间组分;模型常数具有不确定性
6.3.3守恒标量的PDF模型
守恒标量的PDF模型仅适用于扩散(非预混)燃烧问题,该方法假定了反应是受混合速率所控制,即反应已达到化学平衡状态,每个单元内的组分及其性质由燃料和氧化剂
的湍流混合强度所控制。
其中涉及的化学反应体系由化学平衡计算来处理(利用
FLUENT的组件程序PrePDF)。
该方法通过求解混合物分数及其方差的输运方程获得组分和温度场,而不是直接求解组
分和能量的输运方程。
鲁的+口(佛N(2町)+亦+s叶
备(卩严)+▽•—W(牛Vf")+CgM-Cdp*f2+Suwr
其中
/=/-/
混合物分数的定义:
Zk-Zk,O
f
Zk,F-Zk,o
其中Zk代表元素k的元素质量分数,下标F和0分别代表燃料和氧化剂的进口值。
对于简单的燃料/氧化剂体系,每一计算单元内的混合物分数代表了该单元内的燃料质量分数。
由于混合物分数是守恒标量,因此在求解输运方程时不再考虑反应源项。
在该方法中,化学反应认为足够快,体系中的组分立刻达到平衡状态。
化学平衡组分在
混合物空间的分布可示意如下:
p(V)V=lim1'j
FTi
1
:
「p(f)i(f)df
0
其中时均标量i可以是组分浓度、密度或温度。
守恒标量PDF模型的优缺点:
优点:
可以预测中间组分的浓度,可以考虑流动中的耗散现象,可以考虑化学反应与湍流之间的相互作用。
该方法不需求解大量的组分和能量的输运方程,因而可以缩短计算时间。
缺点:
研究的流动体系必须接近于局部化学平衡状态,且不能用于非湍流流动,同时亦不能处理预混燃烧问题。
6.3.4层流火焰面模型
层流火焰面模型的基本思想是把湍流扩散火焰看作是层流对撞扩散火焰面的系综。
该方
不同于守恒标量的函数:
法可以看作是守恒标量PDF模型的一个扩展,它可用于处理非化学平衡状态的体系,即可以利用化学反应动力学的方法处理反应流。
PDF模型,标量是混合物分数和标量耗散率的函数,而非混合物分数
广i(f,)
指定混合物分数f的PDF符合B函数分布,标量耗散率x的PDF符合狄拉克—3函数分
布,因此,时均标量值可以通过在f和x空间求标量的统计平均来得到(即,考虑化学反应与湍流的相互作用):
1:
:
ii(f,)Pf(f)P()ddf
00
层流火焰面模型的计算过程如下:
(1)计算不同标量耗散率下,标量在混合物分数f空间的分布,即求解火焰面方程,
组分方程为:
从而得到标量随混合物分数和标量耗散率的变化关系,即式(6-4),并以火焰面数据
库文件的形式保存结果。
(2)火焰面数据库文件也可由其它软件(如CHEMKIN的OPPDIF生成),若得到的库文件为单标量耗散率,则需计算不同标量耗散率的库文件,最后将它们合并。
(3)利用式(6-5)计算火焰面的PDF库,从而得到时均标量随平均混合物分数和平均混合物方差的变化关系。
(4)利用3.3节的求解平均混合物分数及其平均方差的输运方程的方法,在流场中计算
这两个量,然后再利用得到的PDF库查找时均标量值。
6.3.5预混燃烧的Zimont模型
湍流预混燃烧的化学反应采用反应进度(progressvariable),c,进行表征,如下:
c八f丫『
其中Yp、Ypad分别代表当前和完全绝热燃烧后燃烧产物的质量分数,其取值范围在0
到1之间,0代表未然混合物,1代表已燃混合物。
若用反应进度c代表其平均值,则其输运方程可表达如下:
^(pc)+V*(pflc)=V-Vc)&
上式中平均反应速率项如下求解:
P代表未燃物密度,Ut代表湍流火焰传播速度。
湍流预混燃烧的关键在于求解湍流火焰传播速度(位于湍流火焰表面的法线方向),该
速度受两方面因素的影响:
一是层流火焰传播速度,即决定于燃料和氧化剂的浓度、初
始温度,组分的扩散特性以及化学反应动力学特性;二是有大涡决定的火焰皱褶和拉伸
以及由小涡决定的火焰表面厚度。
根据上述讨论,FLUENT中的湍流火焰传播速度可
表达为:
二Au'
(二)1/4
c
3/41/21/4.1/4
式中,A模型常数,u'速度均方值,
Ut=A(u')U|口一It
5层流火焰传播速度,:
•=k/「cp未燃物的分子
导热系数,|t=CDu'3/;湍流长度尺度,t=lt/u'湍流时间尺度,/U|2化学
反应时间尺度。
为考虑火焰拉伸所导致的吹熄(熄火)现象,在反应速率源项中可乘以一个拉伸因子,G,它代表了拉伸所导致火焰不熄火的概率:
G=rrfc{-保卜(?
)+1]}
式中:
以上各式中出现的一些常数值在FLUENT默认条件下为:
A=0.52,Cd=0.37,pstr=0.26,B=0.5
湍流预混火焰温度的计算:
对绝热火焰,有:
T二(1.-C)Ttt+C?
ad
Tu代表未燃混合物的温度,Tad代表绝热条件下燃烧产物的温度。
对非绝热火焰,解如下的能量方程:
帥)+V-(pvA)=V-(詈吶)++Sg
湍流预混火焰组分密度和浓度的计算:
FLUENT利用理想气体状态方程求解密度,对绝热火焰,假定压力变化可以忽略,混合物平均分子量为常数,有
P皿=內A
其中下标b代表已燃物,下标u代表未燃物,求解量为p,已知量为Tb,p,Tu。
对非绝热火焰,有
求解量为p,已知量为T(由上面的能量方程计算),Pu,Tu。
对于组分浓度,则按下式计算:
Y=Yu(1-c)
6.3.6部分预混燃烧模型
部分预混燃烧系统是指这样一种预混火焰,其燃料/氧化剂之比不唯一。
FLUENT中的
部分预混模型是非预混模型(3.3和3.4节)和预混模型(3.5节)的结合。
预混燃烧的反应进度,c,决定了火焰前锋的位置,在火焰前锋的后面(c=1),混合物已燃,使用
守恒标量PDF或层流火焰面模型的解;在火焰前锋的前面(c=0),组分质量分数,温
度和密度由混合但未燃烧混合物分数来计算。
在火焰内部(0物的线性结合的方法被使用。
部分预混模型求解平均反应进度(以决定火焰前锋的位置)、平均混合物分数f和混
合物分数方差f'2的输运方程。
在火焰前面(c=0),燃料和氧化剂混合但未燃烧,在后面(c=1),混合物已燃。
平均标量(如平均温度、密度和组分质量分数)可由如下的f和c的PDF来计算:
在火焰很薄的假设下,由于存在未燃的反应物和已燃的产物,则平均标量可如下计算:
6.3.7PrePDF4.0使用方法(以FLUENT6.0为例)
(1)开始PrePDF
在WINDOWS2000或WINDOWSNT环境下,用鼠标点击开始一>程序一>FluentInc—>PrePDF4.02(对FLUENT6.0版),贝U进入PrePDF环境。
(2)AllocateMemory
首先为计算的问题设置数组空间和分配内存等,如下:
Setup—>MemoryAllocation
MaximumNumberofSpecies
65。
为PDF表中最大物质数,默认值为20,最大值为
MaximumNumberoff-meanPoints
为PDF表中混合物分数离散点最大个数,默认值为45,最大值为100。
MaximumNumberoff-varPoints
为PDF表中混合物方差离散点最大个数,默认值为22,最大值为30。
MaximumNumberofEnthalpyPoints
为PDF表中焓离散点最大个数,默认值为45,最大值为100。
MaximumNumberofScalarDissipationPointsinAdiabaticFlameletPDFTable为PDF表中标量耗散率离散点最大个数,默认值为45,最大值为100。
MaximumNumberofFlamelets
为火焰面模型中火焰面最大个数,默认值为20,最大值为30。
(3)设定问题
根据所要计算的问题进行设置
Setup—>Case
DefineCase
」SecondaryStream
」PartiaJivPremixed
HeatIranjTerapitions
Model
Chemhtr^models
*Adiabatic
Non-Adiabatic
EquilibriumChEinistr^
StoichiometricReactionyLaminarFlamelets
PDFmodels
7DeltaPDF
*BetaPDF
EmpiricallyDefinedStreams
_lFielstream
riameletOptions
V丈逢強輕或*
Apply
CloseHelp
注:
若在Chemistrymodels中选择LaminarFlamelets,则参见从(13)开始的步骤
(4)定义要考虑的化学组分
(5)定义燃料和氧化剂的初始组分
SpecifyCompositionIn...
Idol#fractions
MassFractioils
SpeciesFrajction
褲馆弗丹械£鮒!
i
ApplyClo»e
(6)定义OperatingConditions
SetuprOperatingConditions
(7)定义化学反应计量系数
例如对反应:
"1:
|-11L1
CH4,O2,CO2和H2O分别定义为1,2,-1,—2(如上图所示)
(8)定义创建的PDF表中所用的参数
SetuprSolutionParameters
(9)将输入保存到文件中
FileTWriteInput
(10)计算并生成PDF表
CalculatePDFTable
(11)将PDF表保存到文件中
FileTWritePDF
(12)计算结果的显示
Display……
02
C02
H20
N2
CO
oie-ki
0
H2
tLllEIhl
±.NE-M
I.IIOII
0-09E4C
«-tOE-01
^OOE-Cl
MiJiUneFnctmF
CHEMICALEQUILIBRIUM
INSTANTANEOUSSPECIESCOMPOSFION
T
J.tfDII
l.HE>II
i.ii&ii
0-09E4C
主mm】
«-tOE-Oi
^OOE-CI
FnctmF
“MF
FlUECktllK.
CHEMICALE口UILIBRILJM
INSTANTANEOUSTEMPERATURE(KELVIN)
上面两图分别给出了由平衡化学计算得出的混合物分数空间标量(组分和温度)的变化。
E
ME*P
ME*D3
±1JE*|J
7.DE*Q2
aJ0E*Q2
JLIIE-IJ
F-MEAN
PDFTABLE-CHEMICALEQUILIBRIUM
MEANFLAMETEMPERATURE
FluEntlnc.
;3-UE-ll
F-lIII!
liftE.
上图给出了绝热体系中固定标量耗散率时的PDF表中温度随平均混合物分数和混合物
分数方差的变化关系。
(13)若使用层流火焰面模型(即使用化学动力学反应而非化学平衡),则重新设置计算的
问题,如下图。
SetupCase
Heattran^reroptions
Adiabatic
Non-Adiabatic
DefineCase
JPartidivPremisedM»d«l
Cherniatrvmodels
EquilibriumChemistryStoichiometricReactionLaminarFlaineleta
PDFmodels
v眈哦為卩於
*旣违u
EmpiricallyDefinedStreams
」fjieJAhftAfti
上面的面板包括了如下信息的输入(a)提供化学反应动力学机理文件的输入,扩展名
为*.che(FLUENT6.0也提供了部分燃料氧化反应的动力学机理供使用),(b)初始燃料
和氧化剂的组分浓度(质量分数和摩尔分数),(c)混合物分数空间的离散点数,标量
耗散率的最大最小值及其离散点数,解的控制参数,辐射选项等。
(15)定义火焰面的操作条件
SetupOperatingConditions
主要是指定燃料和氧化剂来流的温度和压力
(16)定义PDF表的参数
SetupSolutionParameters(同(8))
(17)定义火焰面模型参数
SetupFlameletParameters
(17)存输入文件InputFile
输入文件包括火焰面的定义和其它指定参数。
FileWriteInput...
(18)计算火焰面
CalculateFlamelet
(19)存PDF文件
FileWritePDF…
(关于PDF文件的内容可以在PrePDF中以图形的方式查看,这里不再讨论)得到PDF文件后,在FLUENT的计算中即可调用。
6.6污染物模型6.6.1NOx模型
NOx主要由NO(nitricoxide)所组成,此外还包括少量的NO2和N20。
其中,NO是造成光化学烟雾、酸雨以及臭氧空洞的主要元凶。
FLUENT中提供了NOx产生的如下反应机理:
(1)ThermalNOx-Zeldovich机理(大气中N的氧化产物),在高温条件下,该机理很重要,
(2)PromptNOx-Fenimore机理,该机理对NOx的生成贡献相对较小,但在富燃料区缺很重要,
(3)FuelNOx-经DeSoete,Williams等人发展的经验机理(燃料中含N组分的氧化产物),在富含N的煤燃烧过程中占主要地位(燃烧温度通常较低),
(4)NOxreburnchemistry-通过在富燃料区使NO与碳氢化合物发生反应以减少NOx的化学反应过程。
在FLUENT中将求解NOx的输运方程,对ThermalNOx机理和PromptNOx机理,仅仅求解NO的输运方程:
Q
毎(phio)十V-(puhio)=V-(pPVYno)+
对FuelNOx机理,除求解上述方程外,还要求解与NO生成密切相关的组分NH3和HCN
的输运方程:
£(卩旳cn)+V-(pO%cN)=V-SQFrcn)+5hcn
詁函血+v-(问伽)my1⑷丫呗)+
湍流流动中NO的生成中还要考虑湍流和化学反应之间的相互作用,为此,FLUENT
求解密度加权的时均NS方程以求解温度,速度和组分(或平均混合物分数)。
为计算
NO的生成速率,须利用平均流场的信息求计算区域内每点的时均NO生成速率,可以
采用PDF的方法:
对单变量:
Sno=ISNO(Vi)Pi(Vi)dVi
对两变量:
Sno=IfSN0(Vi^)P{Vi.V2)dVidVi
相同的处理也适用与NH3和HCN。
6.6.2烟(soot)模型
FLUENT提供了烟生成的两种模型,即Khan和Greeves的单步模型(One-stepmodel))以及Tesner的两步模型(Two-Stepmodel)。
对于单步模型,FLUENT求解烟质量分数的输运方程:
爲畑)+W=w(dv心+j
(6-6)
上式中,Rsoot代表烟净生成速率源项,即:
其中,烟形成速率Rsootform由如下的经验关系式来获得:
凡QOtjfuni=GPhiel化E/旳
烟燃烧速率Rsoot,comb则如下计算:
^sootjComb=血。
忆1,弘]
^BDOt^BOOl
以上各式中的符号含义如下:
Cs—模型常数,Pfuel-燃料分压,①一当量比,r—当量比
指数,E/R—活化温度,A—模型常数,Yox,Yfuel—氧化剂和燃料的质量分数,—烟和燃料的计量系数。
对于两步模型,首先预测核(radicalnuclei)的形成,然后再计算在这些核上烟的生成。
因此,FLUENT将计算上述两个标量的输运方程,即式(6-6)和下式:
其中,b*nuc为正交的核浓度,R*nuc为正交的核的净生成速率。
Rsoot仍遵循(6-7)式的形式,其中
^sootjform=-©Moot)®"
对R*nuc则有,
其中:
Rji叫form=用+(/_
^nvCjCQuib=p>mbg"yBOOt
上式中出现的各种参数和系数在FLUEN冲均设有默认值,具体的改动和参见FLUENT6.0
使用手册第17章。
此外,烟模型还可考虑烟与辐射之间的相互作用,此部分内容将在热辐射中介绍。
6.7FLUENT中燃烧模拟计算的步骤和原则
Startin2Dor3D
(1)确定物理模型的应用范围,
(2)划分计算网格(必要时应根据初步计算结果调整网格疏密),
(3)确定求解量和计算收敛判据。
Boundaryconditions
(4)燃烧问题通常对进口边界条件十分敏感,利用已知的(或合理的)速度和标量分布作为边界条件是必要的,
(5)壁面传热对于整个计算也是很重要的,若已知,应指定壁面温度,而非指定边界条件中的内部对流、辐射等。
Initialconditions
(6)尽管稳态问题的解不依赖于初始条件,但很差的初始条件会导致问题不能收敛(由于输运方程的数量和非线性),
(7)对一些燃烧问题,可先求解冷态问题,以此为初始条件求气相燃烧问题,再求解离散相问题,再求解有辐射的问题,
(8)对强旋流,应逐渐增加其涡旋度。
UnderrelaxationFactors
(9)松弛的效果是针对高度非线性问题的,
DecreasethedivergingresidualURFinincrementsof0.1
使用混合物分数PDF模型时应松弛密度(0.5)
对高浮力流应松弛速度
对高速流动应松弛压力
(10)一旦获得稳定解,应尝试增加所有量的松弛因子以尽可能地接近默认值。
Discretization
(11)首先以一阶精度的方法离散控制方程,收敛后再以二阶精度离散以提高计算结果的精度,
(12)对三角形或四边形网格,二阶离散是尤为必要的。
DiscretePhaseModel
为增强计算的稳定性,应
(13)Increasenumberofstochastictracks(oruseparticlecloudmodel)
(14)DecreaseDPMURFandincreasenumberofgasphaseiterationsperDPM
Magnussenmodel
(15)为有限速率/涡耗散方法(Arrhenius/Magnussen)的默认方法,
对非预混(扩散)火焰,应关闭有限速率方法选项预混火焰需要Arrhenius项,因此反应物早期不燃烧
(16)可能需要高温初始化/补丁(initialization/patch),
(17)使用依赖于温度的等压比热Cp以减少高温时的不合理性。
MixturefractionPDFmodel
(18)适合于所研究的问题符合该模型的假设,
(19)在PDF表中使用足够的离散点以保证插值的精确(在不增加计算时间的前提下),
(20)使用betaPDF的形状。
Turbulence
(21)应首先由standardk-;model开始计算,
(22)再转化到RNGk-;,Re