第六章-FLUENT中的燃烧模拟.doc

第六章，FLUENT中的燃烧模拟

6.1燃烧模拟的重要性

面向实际装置〔如锅炉、内燃机、火箭发动机、火灾等〕

面向实际现象〔如点火、熄火、燃烧污染物生成等〕

6.2FLUENT燃烧模拟方法概要

FLUENT可以模拟宽广范围内的燃烧〔反响流〕问题。然而，需要注意的是：你必须保证你所使用的物理模型要适合你所研究的问题。FLUENT在燃烧模拟中的应用可如下列图所示：

气相燃烧模型

气相燃烧模型

离散相模型

离散相模型

预混燃烧

局部预混燃烧

扩散燃烧

输运方程

输运方程

连续性

连续性

动量〔湍流〕

能量

化学组分

热辐射和传热模型

热辐射和传热模型

污染物模型

气相燃烧模型

一般的有限速率形式〔Magnussen模型〕

守恒标量的PDF模型〔单或二组分混合物分数〕

层流火焰面模型〔Laminarflameletmodel〕

Zimont模型

离散相模型

煤燃烧与喷雾燃烧

热辐射模型

DTRM,P-1,Rosseland和DiscreteOrdinates模型

污染物模型

NOx模型，烟〔Soot〕模型

6.3气相燃烧模型

6.3.1燃烧的化学动力学模拟

实际中的燃烧过程是湍流和化学反响相互作用的结果，燃烧的化学反响速率是强非线性和强刚性的。通常的化学反响机理包含了几十种组分和几百个基元反响，而且这些组分之间的反响时间尺度相差很大〔10－9～102秒〕，因此在实际问题的求解过程中计算量和存储量极大，目前应用尚不现实。

在FLUENT中，针对不同的燃烧现象，采用了不同的化学动力学处理手段，以减少计算本钱，如下：

有限速率燃烧模型——预混、局部预混和扩散燃烧

混合物分数方法〔平衡化学的PDF模型和非平衡化学的层流火焰面模型〕——扩散燃烧

反响进度方法〔Zimont模型〕——预混燃烧

混合物分数和反响进度方法的结合——局部预混燃烧

6.3.2一般的有限速率模型

化学反响过程一般采用总包机理〔即简化化学反响，如单步反响〕进行描述

求解组分的输运方程，得到每种组分的时均质量分数值，如下：

6－1

其中组分j的反响源项为所有K个反响中，组分j的净生成速率：

6－2

式中，反响k中的组分j的反响速率可按照Arrhenius公式、混合〔mixing〕速率或“eddybreakup”速率的方法求解。在混合〔mixing〕速率方法中，混合速率和涡的时间尺度,k/?.有关，其物理意义为化学反响受限于湍流导致的组分和热量的混合速率。Ji表达如下：

6－3

计算所需参数包括：〔i〕组分及其热力学参数值；〔ii〕反响及其速率常数值。其中，FLUENT提供了一个混合物组分的数据库可供查找选用，另外也提供了一个化学反响机理以及组分热力学性质的数据库可供查找选用。

有限速率模型的优缺点：

优点：适用于预混、局部预混和扩散燃烧；简单直观

缺点

：当混合时间尺度和反响时间尺度相当时〔即Da1〕缺乏真实性；难以解决化学反响与湍流的耦合问题；难以预测反响的中间组分；模型常数具有不确定性

6.3.3守恒标量的PDF模型

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

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

其中

混合物分数的定义：

其中Zk代表元素k的元素质量分数，下标F和O分别代表燃料和氧化剂的进口值。

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

在该方法中，化学反响认为足够快，体系中的组分立刻到达平衡状态。化学平衡组分在混合物空间的分布可示意如下：

化学反响和湍流之间的相互作用采用几率密度函数〔PDF〕的方法进行处理：

上图代表了几率密度函数p〔V〕的定义，因此在混合物分数空间，f，标量的时均值可由下式计算：

其中时均标量可以是组分浓度、密度或温度。

守恒标量PDF模型的优缺点：

优点：可以预测中间组分的浓度，可以考虑流动中的耗散现象，可以考虑化学反响与湍流之间的相互作用。该方法不需求解大量的组分和能量的输运方程，因而可以缩短计算时间。

缺点：研究的流动体系必须接近于局部化学平衡状态，且不