fluent燃烧简介

1、百度文库-让每个人平等地提升自我FLUENT燃烧简介FLUENT软件中包含多种燃烧模型、辐射模型及与燃烧相关的湍流模型，适用于各种 复杂情况下的燃烧问题，包括固体火箭发动机和液体火箭发动机中的燃烧过程、燃气轮机中的燃烧室、民用锅炉、工业熔炉及加热器等。1.1 FLUENT燃烧模拟方法概要燃烧模型是FLUENT软件优于其它CFD软件的最主要白特征之一。FLUENT可以模拟宽广 范围内的燃烧问题。 然而，需要注意的是：你必须保证你所使用的物理模型要适合你所研究 的问题。FLUENT在模拟燃烧中的应用可如下图所示：图1 FLUENT模拟过程中所需的物理模型1.1.1 气相燃烧模型一般的有限速率形式（

2、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模型和非平衡化学的层流火焰面模型）-扩散燃烧百度文库-让每个人平

3、等地提升自我反应进度方法（Zimont模型）-预混燃烧混合物分数和反应进度方法的结合-部分预混燃烧2.2.1 有限速率模型化学反应过程一般采用总包机理（即简化化学反应，如单步反应）进行描述。求解积分的输运方程，得到每种组分的时均质量分数值，如下：0-战-（1）其中组分j的反应源项为所有反应K个反应中，组分j的净生成速率：&"-(2)计算所需参数包括：1、组分及其热力学参数值；2、反应及其速率常数值。有限速率模型的有缺点：优点：适用于预混、部分预混和扩散燃烧，简单直观；缺点：当混合时间尺度和反应时间尺度相当时缺乏真实性，难以解决化学反应与湍流的耦合问题，难以预测反应的中间组分，

4、模型常数具有不确定性。这种模型求解反应物和生成物输运组分方程，并由用户来定义化学反应机理。反应率作为源项在组分输运方程中通过阿累纽斯方程或涡耗散模型。应用领域：该模型可以模拟大多数气相燃烧问题，在航空航天领域的燃烧计算中有广泛 的应用。2.2.2 守恒标量的PDF模型守恒标量的PDF模型仅适用于扩散（非预混）燃烧问题，该方法假定了反应是受混合速率 所控制，即反应已经达到化学平衡状态，每个单元内的组分及其性质是由燃料和氧化剂的湍 流混合强度所控制， 其中涉及的化学反应体系由化学平衡计算来处理（利用FLUENT的组件程序 PrePDF）。该方法通过求解混合物分数及其方差的输运方程获得组分和温度场，

5、而不是直接求解组分和能量的输运方程。(4)其中(6)混合分数定义氯淖）产）=77山区了）Cg铲4 Siu«kr百度文库-让每个人平等地提升自我其中Zk代表元素k的元素质量分数，下标F和O分别代表燃料和氧化剂的进口值。对于简单的燃料/氧化剂体系，每一计算单元内的混合物分数代表了该单元内的燃料质量分数， 由于混合物分数是守恒标量，因此在求解输运方程时不再考虑反应源项。在该方法中，化学反应认为足够快， 体系中的组分立刻达到平衡状态。化学平衡组分在混合物空间的分布可示意如下:化学反应和湍流之间的相互彳用采用概率密度函数(PDF)的方法处理:上图代表了概率密度函数p(V)的定义，因此在混合物分

6、数空间,f标量的时均值可由下式计算:守恒标量PDF模型的优缺点:优点：可以预测中间组分的浓度,可以考虑流动中的耗散现象,可以考虑化学反应与湍流之间的相互作用，该方法不需要求解大量的组分和能量的输运方程，因此可以缩短计算时间。缺点：研究的流动体系必须接近于局部化学平衡状态,且不能用于非湍流流动，同时亦不能处理预混燃烧问题。该模型不要求用户显式地定义反应机理，而是通过火焰面方法(即混即燃模型)或化学平衡计算来处理，因此比有限速率模型有更多的优势。应用领域：该模型应用于非预混燃烧(湍流扩散火焰),可以用来计算航空发动机的环形燃烧室中的燃烧问题及液体 /固体火箭发动机中的复杂燃烧问题。百度文库-让每个

7、人平等地提升自我2.2.3 层流火焰面模型层流火焰面模型的基本思想是把湍流扩散火焰看作是层流对撞扩散火焰面的系统。该方法可以看作是守恒标量 PDF模型的一个扩展，它可用于处理非化学平衡状态的体系，即可 以利用化学反应动力学的方法处理反应流。不同于守恒标量的 PDF模型，标量是混合物分数和标量耗散率的函数，而非混合物分数的函数：m"ue|（8）指定混合物分数f的PDF符合3函数分布，标量耗散率X的 PDF符合狄拉克-8函数分布， 因此，时均标量值可以通过在 f和x空间求标量的统计平均来得到 （即，考虑化学反应与湍 流的相互作用）：I * 依JJ嫉j/=打,/乃U以犷 s 口（9）层流火

8、焰面模型的计算过程如下：1、计算不同标量耗散率下，标量在混合物分数f空间的分布，即求解火焰面方程，组分方程为：能量方程为:犯 11丁。+&X焉+ ,为期第如（铲叫r-司 ）从而得到标量随混合物分数和标量耗散率的变化关系，即式 （8）,并以火焰面数据库文件的形式保存结果。2、火焰面数据库文件也可由其它软件生成，若得到的库文件为单标量耗散率，则需计算不同标量耗散库的库文件，最后将它们合并。3、利用（9）式计算火焰面的PDF库，从而得到时均标量随平均混合物分数和平均混合物方差 的变化关系。4、利用求解平均混合物分数及其平均方差的输运方程的方法，在流场中计算这两个量，然后再利用得到的 PDF库

9、查找时均标量值。2.2.4预混燃烧的Zimont模型湍流预混燃烧的化学反应采用反应进度（progress variable）,c,进行表征，如下:百度文库-让每个人平等地提升自我匚=£工4丫丁"-(12)其中Yp和Yp(ad)分别代表当前和完全绝热燃烧后燃烧产物的质量分数，其取值范围在 0到1之间，0代表未燃混合物，1代表已燃混合物。若用反应进度c代表其平均值，则其输运方程可表达如下:+ D ,(即)± (&V4 (13)上式中平均反应速率项如下求解:原一(14)P v代表未燃物密度，Ut代表湍流火焰传播速度。湍流预混燃烧的关键在于求解湍流火焰传播速度(位

10、于湍流火焰表面的法线方向),该速度受 两方面馆因素的影响：一是层流火焰传播速度，即决定于燃料和氧化剂的浓度、初始温度、组分的扩散特性以及化学反应动力学特性；二是有大涡褶皱和拉伸以及由小涡决定的火焰表面厚度。根据上述讨论，FLUENTZ中的湍流火焰传播速度可表达为：ut :*仰'尸/ %-.产-(15)式中，A模型常数，u'速度均方值，Ul层流火焰传播速度，& =k/p Cp未燃物的分子导热系数，I=CdU3/g湍流长度尺度，T t为湍流时间尺度，t c为化学反应时间尺度。为考虑火焰拉伸所导致的吹熄现象，在反应速率源项中可乘以一个拉伸因子G,它代表了拉伸所导致火焰不熄火的

11、概率。=#卡一层卜件)+2(16)2.2.5部分预混燃烧模型部分预混燃烧系统是指这样一种预混火焰，其燃料/氧化剂之比不唯一。FLUENT中的部分预混模型是非预混模型和预混模型的结合。预混燃烧的反应进度，c,决定了火焰前锋的位置，在火焰前锋的后面 (c=1),混合物已燃，使用守恒标量 PDF或层流火焰面模型的解； 在火焰前车t的前面(c=0)，组分质量分数，温度和密度由混合但未燃烧混合物分数来计算。在 火焰内部，未燃和已燃混合物的线性结合的方法被使用。部分预混模型求解平均反应进度c,平均混合物分数f和混合物分数方差f2的输运方程。平均标量可由如下的 f和c的PDF来计算：f f(/者女卜-(17

12、)在火焰很薄的假设下，由于存在未燃的反应物和已燃的产物，则平均标量可如下计算：百度文库-让每个人平等地提升自我。匚三，伽COM/）讶+0-司，1£/（力中具体煤粉燃烧算例1、建立求解模型：连续相（气体）只有在非耦合求解时非预混燃烧模型才是可用的。Definemodels-solver.2、打开RNG k-e湍流模型Define-models-viscous.百度文库-让每个人平等地提升自我3、打开非预混燃烧模型Define-models-species.a、在 Model 下选择 Non-Premixed Combustion 。当点OK时，FLUENT将打开一个对话框，要求输入在模

13、拟中要用到的PDF文件。b、在选择文件对话框中，选择并读入非绝热的PDF文件（coal.pdf）当FLUENT读入非绝热PDF文件时，它会自动激活能量求解方程，所以你可以不用打开能 量面板激活传热方程。3、选择P-1以激活辐射模型。百度文库-让每个人平等地提升自我Define-ModelsRadiation.P1模型是能求解气体和颗粒间辐射传热的模型之一。2、建立求解模型：离散相FLIENT会用离散相模型来模拟煤粉的流动。这模型会预示出单个煤粒的轨迹，每一个都代表煤的连续流，在交替计算离散相的轨迹和气相连续方程时，煤粒与气相间热量、动量、质量的传递都将包含其中。1、耦合离散相与连续相流动预报。

14、Define-Models-Discrete Phase.a、在 Interaction 下，选中 Interaction with Continuous Phase 选项。这选项将激活耦合求解，在求解中，离散相的轨迹将会对气相产生影响，如果不选中这选项，你仍可以看到煤粒的轨迹，但上述参数对连续相的流动将没有任何影响。百度文库-让每个人平等地提升自我b、 定义耦合参数， 设定 Number of continous phase interations per DPM interation 为 200;在一些有着高质量粒子和较大网格尺寸的问题中，应该给这参数设定高一点的值， 这对低频 率轨迹是很

15、有好处的，为了更完全地聚合气相方程，应先对轨迹进行反复的计算。c、在 Tracking paramete 下，为 MAX Number of steps 输入 500000。d、打开 Specify Length Scale ,保持 Lenth Scale 的默认值为 0.01m。Length Scale控制离散相轨迹综合中用到的每一次步数的大小。这儿用到的值 0.01m意味着 10m长的一段轨迹要计算1000步左右。e、在 Options T,选择 Particle Radiation Interation 。2、创建离散相煤的射入轨迹。煤粉流用初始条件定义，这初始条件认为煤是进入到气体中。

16、在颗粒的运动方程的每一次综合中，FLUENT将用这些初始条件作为计算的开始点。在这里，煤的总质量流的比率为2.4653kg/s,在10-160微米直径方向上，假设颗粒是服从Rosin-Rammler尺寸分布的。其他的初始条件以及适当的输入程序将在下面做详细说明。Define-Injections.百度文库-让每个人平等地提升自我a、在Injections面板中点击 Creat按钮。b、Injection Type 向下列表中选择 Surface。c、在 Particle Type 下，选择 Combusting。通过Combustion ,可以激活煤液化作用和煤燃尽的子模型。类似地，如果选择D

17、roplet将会激活小液滴蒸发和沸腾的子模型。d、在Material向下列表中选择 Coal-mv。Material列表包含了 FLUENT数据库中的燃烧颗粒物质。你可以在列表中选择一种适当的10百度文库-让每个人平等地提升自我煤，然后在Materials面板中对其进行修改。e、在 Diameter Distribustion 向下列表中选择 rosin-rammler。煤粒的尺寸分布是不均匀的，其直径大小从10微米到160微米，这些尺寸分布是与rosin-rammler方程相适应的。f、在Oxidizing Species向下列表中选择 02。g、在Point Priperties下指定初始

18、条件。h、定义湍流差量。点击Turbulent Disoersion ,面板将示出相关的输入数据。在 Stochastic Tracking 下，打开 Stochastic Modek。随机跟踪可以模拟湍流颗粒轨迹的作用，在模拟实际颗粒的分散， 将随机跟踪包含其中是很 重要的。3、物质：连续相当需要用到非预混燃烧模型时，prePDF化学数据库的热力学数据包括密度，比热和成分的始将被用到，这些性质将作为pre-mixture物质被转换到FLUENT中去，这其中只有可以传输性质，如粘性和导热系数才需要定义。11百度文库-让每个人平等地提升自我DefineMaterials1、设定 Thermal

19、Conductivity 为 X。2、设定 Viscosity3、在Absorption Coefficient 的向下列表中选择 wsggm-cell-based。这个模型指定了一个随着 成分变化而变化的吸热系数。4、点击 Change/Create 按钮。4、物质：离散相。DefineMaterials12百度文库-让每个人平等地提升自我1、从 Materials Type 歹U表中选择 combustion-particle。2、在Combustion Particle Materials 列表中，保持当前的选择（coal-mv）。你可以点击Database.按钮来查看可用的 combu

20、stion-particle物质。3、为物质coal-mv设置下列常量的参数值。13百度文库-让每个人平等地提升自我Dtiisity1300 kg m *Cp1000 Jkg-KThgi tital CoudiK tii in0.0454 uiii-kLMeot Heat0Vaporization400 KVolatile Component Fraction (° &)128Bluary DifTusi'vity5e-4 m ®/sPm hcle Lmissivit'09Pni tkk Staturmg Factor65ISn elliiior|z

21、Bui nout Stokhioineh k Kiuio2.67Combustible Fraction (%)64FLUENT将用到下列输入值：Density将影响颗粒的惯性和质量力。Cp决定了颗粒温度变化所需要的热量。Latent Heat是挥发份蒸发所需要的热量。当在煤燃烧中用非预混燃烧模型的时侯，这值通常设为0。如果为了节省燃料的热量，而选中挥发份成份的话，那么潜热也会包含在其中(但是，如果水被作为气相的水蒸汽而含在挥发分中的话，你就应该用非0的潜热值)。Vaporization Temperature是煤开始液化时的温度。它应该与prePDF中用到的燃料入口温度相等。Volatile

22、 Component Fraction 确定了液化的每一煤粒的质量。Binary Diffusivity颗粒表面的氧化剂的扩散率，用在有扩散率限制的煤燃烬 率方面。Particle Emissivity颗粒的辐射系数。在计算辐射到颗粒上的热量的时侯会 用到。Particle Scattering Factor 颗粒的散射系数。Swelling Coefficient确定了煤在液化的过程中直径的变化。膨胀系数为 2就 意味着在挥发份的挥发过程中颗粒尺寸会变为原来的两倍。Burnout Stoichiometric Ratio用在有扩散控制的燃烬率的计算方面。另外， 在用非预混燃烧模型时，这参数将

23、没有用。当采用有限速率化学的时候，化学 计量比定义了需要单位质量碳的氧化剂的质量。默认值为将C(s)氧化成CO2的所需量。Combustible Fraction是煤粒中碳的质量分数。它决定了被碳燃烬子模型消 耗的每一煤粒的质量。4 .为 Devolatilization Model选择 Single Rate Devolatilization Model 。(a)在 Devolatilization Model向下列表中选择 single-rate 选项。14百度文库-让每个人平等地提升自我这将打开 Single Rate Devolatilization Model 面板。(b)保持面板中

24、默认的液化模型参数。5.为 Combustion Model .选择 kinetics/diffusion-limited。(a)在 Combustion Model向下列表中选择 kinetic/diffusion-limited 选项。这将打开 Kinetics/Diffusion Limited Combustion Model 面板。(b)保持默认值不变。6.点击 Change/Create ,并关闭 Materials 面板。5、边界条件Define Boundary conditionsBojndaiy CcriditionsZonea-h-db a-h-mla-h-m2in le

25、t”nt3-h-m3a-h-pz1 a-h-p ah-zb a-h-zd a-h-2C? a-h-?jt a-h-zj2intake-fan interface m a s&f low-inlet outflow outlet-vent prcssurc"fa-field pressure-inlet pressure-outlet svmmclryCloseinletHelpwall1、速度边界条件15百度文库-让每个人平等地提升自我Zone Namea-h-ni3Momentum |Thermal Radiation Species DPNI Multiphase UD5

26、Velocity Specification Method Magnitude, Normal to BoundaryReference FrameAbsolutecofista ntVelocity Magnitude (m/s 12Th.3TurbulenceSpecification Methodpntenslty ®nd Hydraulic DismeierTurbulent Intensity qgHydraultc Diameter (m) 6.377OK I Ctmcd Help 注意：在这儿将根据其强度和水力直径来定义湍流参数。对于助燃空气流来说10%的湍流强度是很具

27、有象征性的。水力直径的值为计算值。在非预混燃烧计算中，需要定义进口的Mean Mixture Fraction 和MixtureFraction V ariance。在煤燃烧中，所有的燃料都是离散相的，所以进口处的混合分数的值为0。所以，保持默认的 0值不变。2、为wall区域设定条件（炉壁）ZoneTypeinte-shadowinicl-vent4new_outintake-fannew_wsllinterface久w?tHmass-flow-in letwall.2S5outflowwatlloutlet-ventpressure-far-fteMwalllGpressure-inletwalll 7pressure-outletwatllBsymmetrywalllSvelocity-inletwall20wallIDSet i Copy”. Close Help我们可以将炉壁当成是等温边界，温度为 700K16百度文库-让每个人平等地提升自我tfdllMnmtniiiim Thermal Plndiiatinin! | Spntira | 01MMniNiphnBc UDSThenii的 I Co