FLUENT中的燃烧模拟.docx
- 文档编号:23736112
- 上传时间:2023-05-20
- 格式:DOCX
- 页数:26
- 大小:572.39KB
FLUENT中的燃烧模拟.docx
《FLUENT中的燃烧模拟.docx》由会员分享,可在线阅读,更多相关《FLUENT中的燃烧模拟.docx(26页珍藏版)》请在冰豆网上搜索。
FLUENT中的燃烧模拟
第六章,FLUENT中的燃烧模拟
6.1燃烧模拟的重要性
•而向实际装置(如锅炉、内遇机、火箭发动机、火灾等)
•而向实际现象(如点火、熄火、燃烧污染物生成等)
6.2FLUENT燃烧模拟方法概要
•FLUENT可以模拟宽广范国内的燃烧(反应流)问题。
然而,需要注意的是:
你必须保证你所使用的物理模型要适合你所研究的问题。
FLUENT在燃烧模拟中的应用可如下图所示:
•气相燃饶模型
一般的有限速率形式(Magniisscn模型)守恒标量的PDF模型(单或二组分混合物分数)层流火焰面模型(Laminarflameletmodel)Zimont模型
•离散相模型
煤燃烧与喷雾燃饶
•热辐射模型
DTRM.P-1.Rosseland和DiscreteOrdinates模型
•污染物模型
NOX模型,烟(Soot)模型
6.3气相燃烧模型6.3.1燃烧的化学动力学模拟
实际中的燃烧过程是湍流和化学反应相互作用的结果,燃烧的化学反应速率是强非线性和强刚性的。
通常的化学反应机理包含了几十种组分和几百个基元反应,而且这些组分之间的反应时间尺度相差很大(10勺〜IO?
秒),因此在实际问题的求解过程中计算量和存储屋极大,目前应用尚不现实。
在FLUENT中,针对不同的燃烧现象,采用了不同的化学动力学处理手段,以减少计算成本,如下:
•有限速率燃烧模型一一>预混、部分预混和扩散燃烧
•混合物分数方法(平衡化学的PDF模型和非平衡化学的层流火焰而模型)一一>扩散燃烧
•反应进度方法(Zimont模型)一一>预混燃烧
•混合物分数和反应进度方法的结合一一>部分预混燃烧
6.3.2一般的有限速率模型
•化学反应过程一般采用总包机理(即简化化学反应,如单步反应)进行描述
•求解组分的输运方程,得到每种组分的时均质量分数值,如下:
Q—
毎s耳)+厂・(松)=—▽•必+用+&61
其中组分j的反应源项为所有K个反应中,组分j的净生成速率:
k
6-2式中,反应k中的组分j的反应速率可按照Arrhenius公式、混合(mixing)速率或"eddybreakups速率的方法求解。
在混合(mixing)速率方法中,混合速率和涡的时间尺度,k/&有关,其物理意义为化学反应受限于湍流导致的组分和热量的混合速率。
Ji表达如下:
朕十如+註回6_3
•计算所需参数包括:
(i)组分及英热力学参数值:
(ii)反应及貝速率常数值。
其中,FLUENT提供了一个混合物组分的数据库可供查找选用,另外也提供了一个化学反应机理以及组分热力学性质的数据库可供查找选用。
•有限速率模型的优缺点:
优点:
适用于预混、部分预混和扩散燃烧;简单直观
缺点:
当混合时间尺度和反应时间尺度相当时(即Da»l)缺乏真实性;难以解决化学反应与湍流的耦合问题:
难以预测反应的中间组分:
模型常数具有不确左性
6.3.3守恒标量的PDF模型
•守恒标量的PDF模型仅适用于扩散(非预混)燃烧问题,该方法假左了反应是受混合速率所控制,即反应已达到化学平衡状态,每个单元内的组分及其性质由燃料和氧化剂的湍流混合强度所控制。
英中涉及的化学反应体系由化学平衡计算来处理(利用FLUENT的组件程序PrePDF)o
•该方法通过求解混合物分数及英方差的输运方程获得组分和温度场,而不是直接求解组分和能量的输运方程。
|(p7)+v-(^/)=v.(牛呵+s他+s喚
篇(P严)+V•(丽)R(訴严)+Cg也(⑺刀-吠严+S喚
其中
•混合物分数的立义:
yrPQfm.qQFlcwhtlnc.
其中以代表元素£的元素质量分数,下标尸和o分别代表燃料和氧化剂的进口值。
对于简单的燃料/氧化剂体系,每一计算单元内的混合物分数代表了该单元内的燃料
质量分数。
由于混合物分数是守恒标量,因此在求解输运方程时不再考虑反应源项匚
•在该方法中,化学反应认为足够快,体系中的组分立刻达到平衡状态。
化学平衡组分在混合物空间的分布可示意如下:
02
C02
H20
N2
C0
H2
CHEMICALEQUILIBRIUM
INSTANTANEOUSSPECIESCOMPOSITION
•化学反应和湍流之间的相互作用采用几率密度函数(PDF)的方法进行处理:
p(V)AV=
上图代表了几率密度函数P(V)的左义,因此在混合物分数空间,f.标量的时均值可由下式计算:
0
英中时均标疑迓可以是组分浓度、密度或温度。
•守恒标量PDF模型的优缺点:
优点:
可以预测中间组分的浓度,可以考虑流动中的耗散现象,可以考虑化学反应与湍流之间的相互作用。
该方法不需求解大量的组分和能量的输运方程,因而可以缩短计算时间。
缺点:
研究的流动体系必须接近于局部化学平衡状态,且不能用于非湍流流动,同时亦不能处理预混燃烧问题。
6.3.4层流火焰面模型
•层流火焰面模型的基本思想是把湍流扩散火焰看作是层流对撞扩散火焰而的系综。
该方法可以看作是守恒标量PDF模型的一个扩展,它可用于处理非化学平衡状态的体系,即可以利用化学反应动力学的方法处理反应流。
•不同于守恒标量PDF模型,标量是混合物分数和标量耗散率的函数,而非混合物分数的函数:
0=0(/,力)<6_4)
指泄混合物分数./•的PDF符合8函数分布,标量耗散率刀的PDF符合狄拉克一0函数分布,因此,时均标量值可以通过在f和刀空间求标量的统计平均来得到(即,考虑化学反应与湍流的相互作用):
10C
恥JM(/,力)勺(/必⑵d河(6-5)
00
•层流火焰面模型的计算过程如下:
(1)计算不同标疑耗散率下,标疑在混合物分数f空间的分布,即求解火焰而方程,组分方程为:
1
px^~df
dYi11房匕「]
Q而=尹适忑萨+艮25/
脇訴€)(等+P单評)]
能量方程为:
从而得到标量随混合物分数和标量耗散率的变化关系,即式(6-4),并以火焰面数据库文件的形式保存结果。
(2)火焰面数拯库文件也可由苴它软件(如CHEMKIN的OPPDIF生成),若得到的库文件为单标疑耗散率,则需计算不同标量耗散率的库文件,最后将它们合并。
(3)利用式(6-5)计算火焰而的PDF库.从而得到时均标量随平均混合物分数和平均混合物方差的变化关系。
(4)利用3.3节的求解平均混合物分数及其平均方差的输运方程的方法,在流场中计算这两个量,然后再利用得到的PDF库查找时均标量值。
6.3.5预混燃烧的Zimont模型
•湍流预混燃烧的化学反应采用反应进度(progressvariable),c进行表征,如下:
pP
其中s齐N分别代表当前和完全绝热燃烧后燃烧产物的质量分数,其取值范帀在o到1之间,0代表未然混合物,1代表已燃混合物。
•若用反应进度C代表苴平均值,则其输运方程可表达如下:
a
詁pc)+V•(河©=V•
上式中平均反应速率项如下求解:
psc=puUt\Vc\
內代表未燃物密度,U代表湍流火焰传播速度。
•湍流预混燃烧的关键在于求解湍流火焰传播速度(位于湍流火焰表而的法线方向),该速度受两方而因素的影响:
一是层流火焰传播速度,即决泄于燃料和氧化剂的浓度、初始温度,组分的扩散特性以及化学反应动力学特性:
二是有大涡决左的火焰皱褶和拉伸以及由小涡决左的火焰表而厚度。
根据上述讨论,FLUENT中的湍流火焰传播速度可表达为:
U,=A(M,)3,4(//1,2 T.. 式中,A模型常数,/速度均方值,U层流火焰传播速度,a=k/pcp未燃物的分子导热系数,/,=Cdu'3/£湍流长度尺度,T,=lthi'湍流时间尺度,rc=a/Ut2化学反应时间尺度。 为考虑火焰拉伸所导致的吹熄(熄火)现象,在反应速率源项中可乘以一个拉伸因子,G,它代表了拉伸所导致火焰不熄火的槪率: 式中: €cr= BUf 以上冬式中出现的一些常数值在FLUENT默认条件下为: A=O.52»Ci)=0.31,如心0・26,B=0.5 •湍流预混火焰温度的计算: 对绝热火焰,有: ©=(1—c)Tu4- 7;代表未燃混合物的温度,7R代表绝热条件下燃烧产物的温度。 对非绝热火焰,解如下的能量方程: 訥)+y•(剛"•(于%)+Sy+Sg •湍流预混火焰组分密度和浓度的计算: FLUENT利用理想气体状态方程求解密度,对绝热火焰,假左压力变化可以忽略,混合物平均分子量为常数,有 Pb^b—PuTj 英中下标方代表已燃物,下标“代表未燃物,求解疑为处,已知量为D,久,T,对非绝热火焰,有 pT=puTu 求解量为”,已知量为丁(由上面的能量方程计算),內,几。 对于组分浓度,则按下式计算: Y=Y,.(l-c) 6.3.6部分预混燃烧模型 •部分预混燃饶系统是指这样一种预混火焰,其燃料/氧化剂之比不唯一。 FLUENT中的部分预混模型是非预混模型(3.3和3.4节)和预混模型(3.5*i*j)的结合。 预混燃烧的反应进度,c,决迫了火焰前锋的位豊,在火焰前锋的后而(el),混合物已燃,使用守恒标量PDF或层流火焰而模型的解: 在火焰前锋的前而((=0),组分质疑分数,温度和密度由混合但未燃烧混合物分数来汁算。 在火焰内部(0 •部分预混模型求解平均反应进度忑(以决左火焰前锋的位巻)、平均混合物分数了和混合物分数方差户彳的输运方程。 在火焰前而(c=0),燃料和氧化剂混合但未燃烧,在后而(c=l),混合物已燃。 平均标量(如平均温度、密度和组分质量分数)可由如下的./•和c的PDF来计算: JoJo 在火焰很薄的假设下,由于存在未燃的反应物和已燃的产物,则平均标量可如下计算: 不=凯^(/)p(/)df+卩一©/0"⑴P⑺附 6.3.7PrePDF4.0使用方法(以FLUENT6.0为例) (1)开始PrePDF 在WINDOWS2000或WINDOWSNT环境下,用鼠标点击开始一>程序—>FlucntInc—>PrePDF4.02(对FLUENT6.0版),则进入PrePDF环境。 (2)AllocateMemory 首先为计算的问题设置数组空间和分配内存等,如下: Setup—>MemoryAllocation MaximumNumberofSpecies 为PDF表中最大物质数,默认值为20,最大值为65。 MaximumNumberoff-meanPoints 为PDF表中混合物分数离散点最大个数,默认值为45,最大值为100。 MaximumNumberoff-varPoints 为PDF表中混合物方差离散点最大个数,默认值为22,最大值为30。 MaximumNumberofEnthalpyPoints 为PDF表中焰离散点最大个数,默认值为45,最大值为100。 MaxiinuinNumberofScalarDissipationPointsinAdiabaticFlameletPDFTable为PDF表中标量耗散率离散点最大个数,默认值为45,最大值为100。 MaximumNumberofFlamelets 为火焰而模型中火焰面最大个数,默认值为20,最大值为30。 (3)设定问题 根据所要计算的问题进行设置 Setup—>Case DefineCase 」SecondaryStream 」PwtkdlyPremisedModel ”AdiabaticyNon-Adiabatic HeAttransferoptionsChemistrymodeh AEquilibriumChemistry«Stcichiom&tricReactionvLaminarFlamelets EmpiricallyDefinedStre^ins JFuelstream J艮理伽参 PDFmodel? yDeltaPDF eBetaPDF FlameletOptions *时辭 ApplyClge|Help 注: 若在Chemistrymodels中选择LaminarFlamelets,则参见从(13)开始的步骤 (4)定义要考虑的化学组分 (5)左义燃料和氧化剂的初始组分 Setup—>Species—>Composition Lp I6498 (6)左义OperatingConditionsSetup—>OperatingConditions OperatingConditions AbsolutePressure(Pa) 鬻|衣.了乞窮妙e颅&rv{K} Max.了饬驴g撚檢re{K} InletTemperature(IQ Apply|Close| (7)左义化学反应计量系数 SetupTSpeciesTStoichiometry 例如对反应: CH4+2O2TCC>2+2H2O CH4,O2,CO2和H? 0分别定义为1,2,一1,-2(如上图所示)。 (8)定义创建的PDF表中所用的参数 SetupTSolutionParameters (9)将输入保存到文件中 T—>Input (10)计算并生成PDF表 CalculateTPDFTabic (11)将PDF表保存到文件中 TTPDF (12)计算结果的显示 DisplayT MixtwerationF CHEMICALEQUILIBRIUM INSTANTANEOUSSPECIESCOMPOSITION B.IIOII 0.00E*fl0 0“】 kOOE-QJ «.ooE-a) MOE-Q] MixtureTxactionF TErnperlature MONO CHEMICALEQUILIBRIUM INSTANTANEOUSTEMPERATURE(KELVIN) 上面两图分別给出了由平衡化学计算得出的混合物分数空间标量(组分和温度)的变化。 上图给出了绝热体系中固左标量耗散率时的PDF表中温度随平均混合物分数和混合物分数方差的变化关系。 (13)若使用层流火焰而模型(即使用化学动力学反应而非化学平衡),则重新设置计算的问题,如下图。 SetupTCasc (14)左义火焰而 实际上,该过程是立义一个层流对撞扩散火焰的构型及燃料和氧化剂的初始状态。 Han)e<9tGeneration CbcmblryFile Apply Ctewislryfiles... dose Help /uso^/Fluon匕.lno/prsf^l^."W/hyA二wonig.uho 上而的面板包括了如下信息的输入(a〉提供化学反应动力学机理文件的输入,扩展名为*心代(FLUENT6.0也提供了部分燃料氧化反应的动力学机理供使用),(b)初始燃料和氧化剂的组分浓度(质量分数和摩尔分数),(C)混合物分数空间的离散点数,标量耗散率的最大最小值及貝离散点数,解的控制参数,辐射选项等。 (15)左义火焰而的操作条件 Setup-^OperatingConditions 主要是指定燃料和氧化剂来流的温度和压力(16)左义PDF表的参数 SetupTSolutionParameters(同(8)) (17)左义火焰而模型参数 SetupTFlameletParameters (17)存输入文件I叩utFile 输入文件包括火焰而的左义和其它指左参数。 T—>Input... (18)计算火焰而 CalculateTFlamclct (19)存PDF文件 ->~>PDF... (关于PDF文件的内容可以在PrePDF中以图形的方式查看,这里不再讨论)得到PDF文件后,在FLUENT的计算中即可调用。 6.6污染物模型 6.6.1NOX模型 •NOx主要由NO(nitricoxide)所组成,此外还包括少量的NO? 和20。 其中,NO是造成光化学烟雾、酸雨以及臭氧空洞的主要元凶。 FLUENT中提供了NO’产生的如下反应机理: (1)ThermalNOX-Zeldovich机理(大气中N的氧化产物),在高温条件下,该机理很重要, (2)PromptNOX-Fenimore机理,该机理对NOx的生成贡献相对较小,但在富燃料区缺很重要, (3)FuelNO「经DeSoete,Williams等人发展的经验机理(燃料中含N组分的氧化产物),在富含N的煤燃烧过程中占主要地位(燃烧温度通常较低), (4)NO、reburnchemistry-通过在富燃料区使NO与碳氢化介物发生反.应以减少NOx的化学反应过程。 •在FLUENT中将求解NCK的输运方程,对ThermalNOX机理和PromptNO、机理,仅仅求解NO的输运方程: 8 毎(“K 对FuelNO、机理,除求解上述方程外,还要求解与NO生成密切相关的组分NFB和HCN的输运方程: g^(plHCN)+V•(awYhcn)=V•(pDYhcn)+Shcn 8 ^(plNHs)+V•(awVnh3)=▽•(p^Vnhs)+Snh$ •湍流流动中NO的生成中还要考虑湍流和化学反应之间的相互作用,为此,FLUENT求解密度加权的时均NS方程以求解温度,速度和组分(或平均混合物分数)。 为计算NO的生成速率,须利用平均流场的信息求讣算区域内每点的时均NO生成速率,可以采用PDF的方法: 对单变呈: : Sno=/Sno(VI)丹(Vl)dVi 对两变量: 5no=jfSno(旳,旳,)P(旳,旳)d! 4d旳 相同的处理也适用与NHs和HCN“ 662烟(soot)模型 •FLUENT提供了烟生成的两种模型,即Khan和Greeves的单步模型(One-stepmodel))以及Tesner的两步模型(Two-Stepmodel)a •对于单步模型,FLUENT求解烟质量分数的输运方程: 爲S)+V•(出+碣 x^soot/(6-6) 上式中,Kg代表烟净生成速率源项,即: 其中,烟形成速率由如下的经验关系式来获得: 烟燃烧速率Rsooucomb则如下计算: ^sootjcomb=min|[冗冗2] 忌=ApY^oot- 隔_彳(吒X)/^oot"soot)£ ^soot^soot+^fuel^fuel/k以上各■式中的符号含义如下: C: 一模型常数,匕吹一燃料分压,①一当虽比,r-当量比指数,E/R—活化温度,A—模型常数,Yux,丫冋一氧化剂和燃料的质量分数,joe,vluel一烟和燃料的计量系数。 •对于两步模型,首先预测核(radicalnuclei)的形成,然后再计算在这些核上烟的生成。 因此,FLUENT将汁算上述两个标量的输运方程,即式(6-6)和下式: ^(Ac)+v-«uc)=v.(总叽)+砾 其中,bj"•为正交的核浓度,RJuc为正交的核的净生成速率。 Rs®仍遵循(6-7)式的形式,其中 对RSuc则有, 其中: ^nWjform=3+(丁一9)Yg-的C羸Moot 事—E[ 乃0=GaCfuelE 上式中出现的各种参数和系数在FLUENT中均设有默认值,具体的改动和参见FLUENT6.0使用手册第17章。 ♦此外,烟模型还可考虑烟与辐射之间的相互作用,此部分内容将在热辐射中介绍。 6.7FLUENT中燃烧模拟计算的步骤和原则 ♦Startin2Dor3D (1)确定物理模型的应用范围, (2)划分讣算网格(必要时应根据初步il•算结果调整网格疏密), (3)确定求解量和计算收敛判据。 ♦Boundaryconditions (4)燃烧问题通常对进口边界条件十分敏感,利用已知的(或合理的)速度和标量分布作为边界条件是必要的, (5)壁面传热对于整个计算也是很重要的,若已知,应指立壁而温度,而非指定边界条件中的内部对流、辐射等。 (6)尽管稳态问题的解不依赖于初始条件,但很差的初始条件会导致问题不能收敛(由于输运方程的数量和非线性), (7)对一些燃烧问题,可先求解冷态问题,以此为初始条件求气相燃烧问题,再求解离散相问题,再求解有辎射的问题, (8)对强旋流,应逐渐增加英涡旋度。 ♦UnderrelaxationFactors (9)松弛的效果是针对高度非线性问题的, ■DecreasethedivergingresidualURFinincrementsof0」 ・使用混合物分数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-smodel开始计算, (22)再转化到RNGRealizablek-s或雷诺应力模型(RSM)以获得与实验数拯较好
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- FLUENT 中的 燃烧 模拟