第13章物质输送及有限速率化学反应.doc

第十三章 组分输运和有限速率化学反应
FLUENT可以通过求解描述每种组成组分的对流、扩散和反应源的守恒方程来模拟混合和输运，可以模拟多种同时发生的化学反应，反应可以是发生在大量相（容积反应）中，和/或是壁面、微粒的壁面。包括反应或不包括反应的组分输运模拟能力，以及当使用这一模型时的输入将在本章中叙述。
注意你可能还希望使用混合物成分的方法（对非预混系统，在14章介绍）、反应进程变量的方法（对预混系统，在15章介绍），或部分预混方法（在16章介绍）来模拟你的反应系统。见12章FLUENT中反应模拟方法的概述。
本章中的分为以下章节：
l 14.1 容积反应
l 14.2 壁面壁面反应和化学蒸汽沉积
l 14.3 微粒壁面反应
l 14.4 无反应组分输运

14．1 容积反应
与容积反应有关的组分输运和有限速率化学反应方面的信息在以下小节中给出：
l 14.1.1 理论
l 14.1.2 模拟组分输运和反应的用户输入概述
l 14.1.3 使能组分输运和反应，并选择混合物材料
l 14.1.4 混合物和构成组分的属性定义
l 14.1.5 定义组分的边界条件
l 14.1.6 定义化学组分的其他源项
l 14.1.7 化学混合和有限速率化学反应的求解过程
l 14.1.8 组分计算的后处理
l 14.1.9 从CHEMKIN导入一个化学反应机理

14.1.1 理论
组分输运方程
当你选择解化学组分的守恒方程时，FLUENT通过第种组分的对流扩散方程预估每种组分的质量分数，Yi。守恒方程采用以下的通用形式：
（14.1-1）
其中是化学反应的净产生速率（在本节稍后解释），为离散相及用户定义的源项导致的额外产生速率。在系统中出现N种组分时，需要解N-1个这种形式的方程。由于质量分数的和必须为1，第N种组分的分数通过1减去N-1个已解得的质量分数得到。为了使数值误差最小，第N种组分必须选择质量分数最大的组分，比如氧化物是空气时的N2。
层流中的质量扩散
在方程14.1-1中，是组分的扩散通量，由浓度梯度产生。缺省时，FLUENT使用稀释近似，这样扩散通量可记为：
（14.1.2）
这里是混合物中第种组分的扩散系数。
对于确定的层流流动，稀释近似可能是不能接受的，需要完整的多组分扩散。在这些例子中，可以解Maxwell-Stefan方程，详细情况见7.7.2节。
湍流中的质量扩散
在湍流中，FLUENT以如下形式计算质量扩散：
(14.1.3)
其中是湍流施密特数，（缺省设置值为0.7）。
注意，湍流扩散一般将淹没层流扩散，在湍流中指定详细的层流性质是得不到保证。
能量方程中的组分输运处理
在许多多组分混合流动中，组分扩散引起焓的传递。

这种扩散对于焓场有重要影响，不能被忽略。特别是，当所有组分的Lewis数
（14.1-4）
远离1时，忽略这一项会导致严重的误差。
FLUENT缺省地包含这一项。在方程14.1-4中，k为热导率。
进口处的扩散
在FLUENT的非耦合求解器中，入口的组分净输运量包括对流分量和扩散分量，对于耦合解算器，只包括对流分量。对流分量由用户指定的组分浓度确定。扩散分量取决于入口处计算所得的组分场。因此，扩散分量（从而使净输运量）不预先指定。见14.1.5节有关指定入口净输运量的信息。
反应建模的通用有限速率形式
反应速率作为源项在方程14.1-1中出现，在FLUENT中根据以下三种模型中的一个计算：
层流有限速率模型：忽略湍流脉动的影响，反应速率根据Arrhenius公式确定。
涡耗散模型：认为反应速率由湍流控制，因此避开了代价高昂的Arrhenius化学动力学计算。
涡耗散概念（EDC）模型：细致的Arrhenius化学动力学在湍流火焰中合并。注意详尽的化学动力学计算代价高昂。
通用有限速率对于范围很广的应用，包括层流或湍流反应系统，预混、非预混、部分预混燃烧系统都适用。
层流有限速率模型
层流有限速率模型使用Arrhenius公式计算化学源项，忽略湍流脉动的影响。这一模型对于层流火焰是准确的，但在湍流火焰中Arrhenius化学动力学的高度非线性，这一模型一般不精确。对于化学反应相对缓慢、湍流脉动较小的燃烧，如超音速火焰可能是可以接受的。
化学组分的化学反应净源项通过有其参加的NR个化学反应的Arrhenius反应源的和计算得到。

其中是第种组分的分子量，为第种组分在第r个反应中的产生/分解速率。反应可能发生在连续相反应的连续相之间，或是在壁面沉积的壁面处，或是发生在一种连续相组分的演化中。
考虑以如下形式写出的第r个反应：
（14.1-6）
其中——系统中化学组分数目；
——反应r中反应物i的化学计量系数;
——反应r中生成物i的化学计量系数;
——第i种组分的符号；
——反应r的正向速率常数；
——反应r的逆向速率常数；
方程14.1-6对于可逆和不可逆反应（FLUENT中缺省为不可逆）都适用。对于不可逆反应，逆向速率常数简单地被忽略。
方程14.1-6中的和是针对系统中的所有组分，但只有作为反应物或生成物出现的组分才有非零的化学计量系数。因此，不涉及到的组分将从方程中清除。
反应r中组分i的产生/分解摩尔速度以如下公式给出：
（14.1-7）
其中：——反应r 的化学组分数目；
——反应r中每种反应物或生成物j的摩尔浓度；
——反应r中每种反应物或生成物j的正向反应速度指数；
——反应r中每种反应物或生成物j的逆向反应速度指数；

见14.1.4节有关输入整体正向反应（不可逆）和单元反应（可逆）的化学计量系数和速率指数方面的内容。
表示第三体对反应速率的净影响。这一项由下式给出：

其中为第r个反应中第j种组分的第三体影响。在缺省状态，FLUENT在反应速率计算中不包括第三体影响。但是当你有它们的数据时，你可以选择包括第三体影响。
反应r的前向速率常数通过Arrhenius公式计算：

其中，——指数前因子（恒定单位）；
——温度指数（无量纲）；
——反应活化能();
——气体常数（）
你（或者数据库）可以在FLUENT的问题定义中提供，并可选择提供。
如果反应是可逆的，逆向反应常数可以根据以下关系从正向反应常数计算：

其中为平衡常数，从下式计算：

其中表示大气压力（101425Pa）。指数函数中的项表示Gibbs自由能的变化，其各部分按下式计算：

其中和是标准状态的熵和标准状态的焓（生成热）。这些值在FLUENT中作为混合物材料的属性指定。

压力决定反应
FLUENT可以用以下三种方法之一来表示压力决定反应（或压力下降）反应的速率表达式。“压力下降”反应是发生在Arrhenius高压和低压极限之间的反应，因而不再单独取决于温度。
有三种方法表示在“fall-off”区域的速率表达式，最简单的是Lindemann[140]形式。还有其它两种相关的方法，Troe方法[77]和SRI方法[230]，它们提供了更精确的”fall-off”区域表达形式。
Arrhenius速率参数对于高压和低压极限都是需要的。两个极限的速率系数融合以产生光滑的压力决定表达式。在Arrhenius形式中，高压极限和低压极限的参数如下：

在任意压力下，净反应速率常数为：

其中定义为：

为溶液气体的浓度，可以包括第三体效率。如果方程14.1-16函数F为1，则是Lindemann形式。FLUENT提供了两种其他形式来表述F，称为Troe方法和SRI方法。
在Troe方法中，F按下式给出：

其中，

参数做为输入确定。
在SRI方法中，混合函数F近似为

其中

除了低压极限表达式中的三个Arrhenius参数以外，你还需要提供F表达式中的a, b, c, d, e。
！化学动力学机理通常为一组高度非线性、且强烈耦合的方程。求解过程指导见14.1.7节。如果你有一个CHEMKIN形式的化学反应机理[112]，你可以将这一机理导入FLUENT，如14.1.9节。

涡耗散模型
大多数燃料快速燃烧，且总反应速率由湍流混合控制。在非预混火焰中，湍流缓慢地使燃料和氧化剂对流/混合，从而进入反应区，在反应区快速燃烧。在预混火焰中，湍流对流/混合冷反应物和热生成物进入反应区，从而迅速发生反应。在这些情况下，燃烧称为混合限制、复杂、未知的，通常化学反应动力学速率可以安全地忽略。
FLUENT提供了湍流-化学反应相互作用模型，基于Magnussen 和Hjertager[149]的工作，称为涡耗散模型。

反应r中组分的产生速率由下面两个表达式中较小的一个给出：

在方程14.1-25和14.1-26中，化学反应速率由大涡混合时间尺度控制，如同Splading[227]的涡破碎模型一样。只要湍流出现（），燃烧即可进行，不需要点火源来启动燃烧。通常对于非预混火焰这是可以接受的，但在预混火焰中，反应物一进入计算区域（火焰稳定器上游）就开始燃烧。为了修正这一点，FLUENT提供了有限速率/涡耗散模型，其中Arrhenius（方程14.1-7）和涡耗散(方程14.1-25和14.1-26)反应速率都进行计算。净反应速率取两个速率中较小的。实际上，Arrhenius反应速率作为一种动力学“开关”，阻止反应在火焰稳定器之前发生。一旦火焰被点燃，涡耗散速率通常会小于Arrhenius反应速率，并且反应是混合限制的。
！尽管FLUENT允许采用涡耗散模型和有限速率/涡耗散模型的多步反应机理（反应数>2），但可能会产生不正确的结果。原因是多步反应机理基于Arrhenius速率，每个反应的都不一样。在涡耗散模型中，每个反应都有同样的湍流速率，因而模型只能用于单步（反应物—产物）或是双步（反应物—中间产物，中间产物—产物）整体反应。模型不能预测化学动力学控制的组分，如活性组分。为合并湍流流动中的多步化学动力学机理，使用EDC模型（下面介绍）。
！涡耗散模型需要产物来启动反应（见方程14.1-26）。当你初始化定常流的解时，FLUENT设置所有产物的质量分数为用户指定的初始值，其最大值0.01，通常足够启动反应。但是，如果你最初收敛一个混合解，其中所有的产物的质量分数均为0，你可能必须在反应区域中补入产物以启动反应。详细内容见14.1.7节。

LES的涡耗散模型
当使用LES湍流模型时，湍流混合速率（方程14.1-25和14.1-26中的）被亚格子尺度混合速率替代。计算为：

涡-耗散-概念（EDC）模型
涡-耗散-概念（EDC）模型是涡耗散模型的扩展，包括湍流流动中详细的化学反应机理[148]。它假定反应发生在小湍流结构中，称为良好尺度。良好尺度的length fraction按下式建模[80]：

其中*表示良好尺度量，
——容积比率常数=2.1477；
——运动粘度

良好尺度的体积分数被计算为，species are assumed to react in the fine structures over a time scale.
其中为时间尺度常数，等于0.4082
在FLUENT中，良好尺度中的燃烧视为发生在定压反应器中，初始条件取单元中当前的组分和温度。反应经过时间尺度后开始进行，由方程14.1-7的Arrhenius速率控制，并且用ISAT法则进行数值积分[45]。ISAT法则可通过2到3个量加速化学计算，从而减少运行时间。ISAT功能强大，但要求用户必须仔细。
经过一个时间的反应后组分状态记为
组分的守恒方程14.1-1中的源项计算公式为：

EDC模型能将详细的化学反应机理并入湍流反应流动中。但是，典型的机理总具有刚性，其数值积分计算代价昂贵。因而，只有在快速化学反应这一假设无效的情况下才使用该模型，例如在快速熄灭火焰中缓慢的CO烧尽、在选择性非催化还原中的NO转化。
推荐使用双精度求解器以避免刚性机理中固有的大指数前因子和活化能产生的舍入误差。见14.1.7节获得使用EDC模型求解的指导。

14.1.2 组分传输及反应建模的用户输入回顾
设定涉及组分输运和反应问题的基本步骤如下，每一步的详细执行过程见14.1.3-14.1.5节。有关设定和求解的附加信息在14.1.6-14.1.8节中提供。
1． 选定组分输运和容积反应，指定混合物材料。见14.1.3节（混合物材料概念在下面解释）。
2． 如果你还要模拟壁面或颗粒壁面反应，则要打开壁面和/或微粒壁面反应。细节见14.2和14.3节。
3． 检查和/或定义混合物的属性。（见14.1.4节）。混合物属性包括：
l 混合物中的组分
l 反应
l 其他物理属性（如粘度、比热）
4． 检查和/或设置混合物中单个组分的属性（见14.1.4节）
5． 设置组分边界条件（见14.1.5节）
在很多情况下，当你选择混合物材料是，求解器从材料数据库中得到组分性质、反应等，因而你将不需要修改任何物理属性。但有一些性质可能在数据库中没有定义。如果有任何性质需要设置时，你将被警告，这样你可以指定这些性质的适当值。你还可能希望检查数据库中这些性质的值，以确定它们对你的应用是否正确。修改已存在的混合物材料或从最开始创建一个新的材料的详细内容见14.1.4节。混合物材料的修改包括以下方面：
l 组分的添加和删除
l 改变化学反应
l 修改混合物的其他材料属性
l 修改混合物本构组分的材料属性
如果求解一个反应流问题，你常常希望将混合物的比热定义为组分的函数，将每种组分的比热定义为温度的函数。你还可能对其它一些性质希望也做这样的定义。缺省时，将使用常量属性，但对一些组分的性质，在数据库中存在一个温度的分段多项式函数可供使用。如果你知道更多适合于你的问题的函数的话，你还可以选择指定一个不同的温度依赖函数。

混合物材料
在FLUENT中提出混合物材料的概念，从而便于设置组分输送和反应流动。混合物材料被以认为是一组组分和一列控制它们相互作用的规则。混合物材料包括以下性质：
l 一列本构组分，涉及“流体”组分
l 一列混合定律，指示如果希望得到组分依赖的属性，混合属性（密度，粘度，比热等）如何从单个组分得到
l 如果希望属性不依赖组分，直接指定混合物属性
l 其它与单个组分无关的材料属性（如吸收和辐射系数）
l 一组反应，包括反应类型（有限速率，涡耗散等）和化学计量和速率常数
混合物材料和流体材料都储存在FLUENT的材料数据库中。包括许多常见的混合物材料（如甲烷-空气，丙烷-空气）。通常，数据库中定义了一步/两步反应机理和大量混合物及其构成组分的属性。当指定了你希望使用哪种混合物材料后，适当的混合物材料，流体材料和属性将被装载到求解器中。如果缺少任何所选材料（或构成流体材料）必须的属性，求解器将通知你需要指定它。另外，你可以选择修改任何预定义的属性。见7.1.2节了解有关FLUENT数据库属性数据源的信息。
例如，如果你计划模拟一种甲烷-空气的燃烧，你不需要明确指定反应中涉及的组分和反应本身。只需要简单地选择甲烷-空气作为使用的混合物材料，相关的组分（CH4，O2，CO2，H2O和N2）和反应数据将从数据库装入求解器。然后你可以检查组分、反应和其它属性，并定义其它缺少的属性，和/或修改任何你希望使用不同值或函数的属性。通常你希望定义一个与组分、温度相关的比热，还可能希望将其它属性定义为温度和/或组分的函数。
混合物材料的使用具有灵活性，可以使用大量预定义混合物中的一种，修改这些混合物，或是创建你自己的混合物材料。自定义混合物材料在Materials面板中进行，在14.1.4节中讲述。

14.1.3 选定组分输运和反应，并选择混合物材料
组分输运和容积反应的问题设置总组分模型（Species Model）面板开始（图14.1.1）。

1． 在Model下，选择Species Transport.
2． 在Reaction下，选择Volumetric reactions
3． 在Mixture Properties下的Mixture Material下拉列表中选择在你的问题中希望使用的混合物材料
下拉列表中将包括所有在当前数据库中定义的混合物。为检查一种混合物材料的属性，选择它，并点击View…按纽。如果你所希望使用的混合物不在列表中，选择混合物模板（mixture-template）材料，并参阅14.1.4节了解设置你自己的混合物属性的详细内容。如果有一种混合物材料和你你所希望使用的混合物相似，可以选择这一材料并参阅14.1.4节修改已存在材料性质的详细信息。
当选择Mixture Material时，混合物中的Number of Volumetric Species将在面板中显示，表达你的信息。
!注意如果你在已经选定组分输运后，重新打开Species Model面板时，只有你的案例中可得到的混合物材料才会显示在列表中。你可以通过从数据库中拷贝在你的案例中增加更多混合物材料，如7.1.2节所述，或是创建一个新的混合物，如7.1.2节和14.1.4节所述。
正如在14.1.2节中提到的，组分输运的模拟参数和反应（如果有关）将自动从数据库中装入。如果缺少任何信息，当你点击Species Model面板中的OK按纽后将被告知缺少什么。如果你希望检查或修改混合物材料的任何属性，你将使用Materials面板，如14.1.4节所述。
4． 选择湍流-化学反应相互作用模型，可以使用四种模型：
层流有限速率：只计算Arrhenius速率（见方程14.1-7），并忽略湍流-化学反应相互作用。
涡耗散模型（针对湍流流动）： 只计算混合速率（方程14.1-25和方程14.1-26）。
有限速率/涡耗散模型（针对湍流流动）：计算Arrhenius速率和混合速率，并使用其中较小的一个。
EDC模型（湍流流动）：使用详细的化学反应机理模拟湍流-化学反应相互作用（见方程14.1-25和14.1-26）。
5． 如果你选择EDC模型，你可以选择修改容积比率常数和时间尺度常数（方程14.1-28中的和方程14.1-29中的，尽管通常推荐缺省值。此外，为减少化学反应计算的开销，你可以增加每次化学反应更新的流动迭代（Flow Iteration Per Chemistry Update）次数。缺省时，FLUENT每十次流动迭代更新化学反应一次。
6． （可选）如果你希望模拟完整的多组分扩散或热扩散，打开完整多组分扩散或热扩散Full Multicomponent Diffusion 或Thermal Diffusion选项。
7. 对于层流反应，激活KINetics from Reaction Design，允许你使用反应速率，并利用Reaction Design的求解规则，该规则基于chemkin技术。对于EDC湍流－化学相互作用和组分PDF输运模型，激活KINetics from Reaction Design，将允许你使用Reaction Design中KINetics模块的反应速率，取代FLUENT默认的反应速率。FLUENT的ISAT规则用于积分这些速率请参考Fluent manual中的KINetics。

14.1.4 定义混合物的属性和构成组分
如14.1.2节所讨论的，如果你使用来自数据库的混合物材料，大部分混合物和组分属性已经定义了。你可以跟随这一节的过程检查当前的属性、修改某些属性或设定一个全新混合物材料的所有属性。
记住你将需要定义混合物材料和其构成组分的属性。由于组分属性输入可能依赖于你定义混合物属性的方法，在设定构成组分的属性之前定义混合物属性非常重要。建议按照如下的属性输入顺序：
1． 定义混合物组分和反应，定义混合物物理属性。记住在设定混合物材料属性时点击Change/Create按纽。
2． 定义混合物中组分的物理属性。记住在设定了每种组分的属性后点击Change/Create按纽。
所有这些步骤都在Materials面板中进行，在本节中将详细叙述。

定义混合物中的组分
如果你使用数据库中的混合物材料，混合物中的组分已经为你定义。如果你创建你自己的材料或修改已存在材料中的组分，你将需要自己定义它们。
在Materials面板中（图14.1.2），检查材料类型Material Type是否已经设置为混合物，并且你的混合物是否已经在混合物材料列表Mixture Materials list中选定。点击Mixture Species右边的Edit…按纽打开Species面板（图14.1.3）。

Species面板概览
在Species面板中，已选组分Selected Species列表显示所有混合物中的流体相组分。如果你模拟壁面或微粒壁面反应，已选组分Selected Species列表将显示所有混合物中的壁面组分。壁面组分是那些从壁面边界或是离散相微粒(如Si(s))产生或散发出来的，以及在流体相组分中不存在的组分。壁面组分和壁面反应将在14.2节中叙述，微粒壁面反应的有关内容见14.3节。
！已选组分Selected Species列表中组分的顺序非常重要。FLUENT认为列表中最后的组分是大量的组分。因此，当你从混合物材料中增加或是删除组分时，必须小心将最丰富（按质量）的组分作为最后一个组分。
可获得的组分Available Material列表显示可获得，但不在混合物中的材料。通常你可以在列表中看到空气air，因为缺省时，空气通常是可获得的。

在混合物中增加组分
如果你从头创建了一个混合物，或者从一种已存在的混合物开始并增加一些缺少的组分，你首先需要从数据库中导入希望的组分（或是创建它们，如果它们不在数据库中出现的话），以便在求解器中可以得到，增加组分的过程在下面列出。（在开始前你需要关闭Species面板，因为它是一个模式面板，不允许在它打开的时候做别的事）
1． 在Materials面板中，点击Datebase…按纽以打开Database Materials面板，并拷贝希望的组分，如7.1.2节所述。记住混合物的构成组分是流体材料，因而你需要在Database Materials面板中选择Fulid Material Type来看到正确的选择列表。注意可获得的固体组分（对壁面反应）同样也在流体列表中。
！ 如果你在数据库中没有见到你搜寻的组分，可以创建对这种组分创建一个新的流体材料，依造7.1.2节中的指导，然后继续以下的第二步。
2. 再次打开Species面板。将看到你从数据库拷贝（或创建）的流体材料在Available Materials列表中列出。
3. 为了在混合物中增加一种组分，在Available Materials列表中选择它，并点击Selected Species列表下面（或是Selected Surface Species列表下面，以定义一种壁面组分）的Add按纽。这种组分将从Available Materials列表移动到Selected Species(或Selected Surface Species)列表的末尾。
l 对所有希望的组分重复上述步骤。结束后点击OK按纽。
！ 增加一种组分到列表中将改变组分的顺序。你必须确保列表中的最后一个组分是大量的组分，并且需要检查以前所设定的所有边界条件，欠松弛因子或其它求解参数

在混合物中去除组分
为从混合物中去除一种组分，只需在Selected Species列表（或Selected Surface Species列表）中选定它，并点击列表下的Remove按纽。这种组分将从列表中去除并增加到Available Materials列表中。
！ 去除列表中的一种组分将改变组分的顺序。你必须确保列表中的最后一个组分是大量的组分，并且你需要检查你以前所设定的所有边界条件，欠松弛因子或其它求解参数，如详细叙述的那样。

重排组分
如果你发现Selected Species列表中最后一种组分不是最丰富的组分（应该是），你需要重排组分以得到正确的顺序。
1． 将最大量的组分从Selected Species列表中去除。它将出现在Available Materials列表中。
2． 再次增添这种组分，它将自动放在列表的末尾。

组分的命名和顺序
如上所述，你必须在增加或去除组分时保持最丰富的组分作为Selected Species列表中的最后一种组分。在增加或去除组分时还需注意以下事项。
一种组分有三个特征在供求解器识别：名字，化学分子式和在Species面板中组分列表中的位置。改变这些特征将会产生以下效果：
l 你可以改变一种组分的名字（使用Materials面板，如7.1.2节所述），而不产生任何影响。
l 你不能改变一种组分的给定的化学分子式。
l 如果你增加或去除组分，将改变组分列表的顺序。这时，所有的边界条件，求解器参数和组分的求解数据将被从设为缺省值。（其它流动变量的求解数据，边界条件，求解器参数将不受影响）因此，如果你增加或去除组分，你需要注意这一新定义问题的组分边界条件和求解参数。另外，你必须认识到基于原来组分顺序给出的组分浓度或是存储在数据文件中的浓度将会与新定义的问题不匹配。你可以使用数据文件作为初始猜测值，但你必须知道数据文件中的组分浓度将可能对新定义的模型提供一种不好的初始猜测。
定义反应
如果你的FLUENT模型中涉及化学反应，可以接着定义参与的已定义组分的反应。这只有在你从头开始创建一种混合物、修改了组分或是出于其他原因希望重定义反应时才是必须的。
在Materials面板的Reaction下拉列表中显示适当的反应机理，取决于你在Species Model面板中选择的湍流-化学反应相互作用模型（见14.1.3节）。如果你使用层流有限速率或EDC模型，反应机理将是有限速率的，如果你使用涡耗散模型，反应机理将是涡耗散的；如果使用有限速率/涡耗散模型，反应机理将是有限速率/涡耗散的。

反应定义的输入
为定义反应，点击Reaction右侧的Edit…按纽。
将打开Reaction面板（图14.1.4）。
定义反应的步骤如下：
1． 在Total Number of Reaction区域中设定反应数目（容积反应，壁面反应和微粒壁面反应）（使用箭头改变数值，或是键入值并按回车键。
注意如果你的模型包括离散相的燃烧微粒，只有在你计划使用壁面燃烧的多壁面反应模型时，才必须在反应数目中包括部分壁面反应（s）（如碳的燃烧，多样碳粒氧化）

2. 设定你反应的Reaction Name
3．设定你反应的ID
4. 如果是流体相反应，保持缺省选项Volumetric作为反应类型。如果是壁面反应（在14.2节中描述）或者颗粒壁面反应（14.3节描述），选择Wall Surface或Particle reaction 作为反应类型。有关定义颗粒壁面反应的进一步信息见14.3.2节。
3． 通过增加Number of Reactants和Number of Products的值，指定反应中反应物和生成物的数量。在Species下拉列表中选择每一种反应物或生成物，然后在Stoich. Coefficient和Rate Exponent中设定化学计量系数和速率指数。（化学计量系数是方程14.1-6中的常数和，速率指数是方程14.1-7中的常数和。共有两种普通类型的反应可以在Reactions面板中处理。因此正确输入每种反应的参数非常重要。反应的类型如下：
l 整体正向反应（无逆向反应）：产物一般不影响正向速率，因此所有产物的速率指数（）应该为0。对于反应物，设定速率指数（）为期望的值。如果某种反应不是基元反应，速率指数一般不等于这种组分的化学计量系数。整体正向反应的一个实例是甲烷的燃烧：

图14.1.4表示了甲烷燃烧的系数输入（还可在DataBase Materials面板中查阅甲烷/空气混合物材料）。
注意：在某些情况下，你可能希望模拟产物影响正向速率的反应。对于这些情况，设定产物速率指数（）为期望的值。这种反应的一个例子是气体转换反应（见DataBase Materials面板中的CO/空气混合物材料），其中水的存在对反应速率有影响：

在气体转换反应中，速率表达式可以定义为：

l 可逆反应：假定每种组分的化学计量系数等于速率指数，这种反应的一个例子是SO2氧化为SO3：

参见下面的第6步如何选定可逆反应。
4． 如果你使用 层流/有限速率，有限速率/涡耗散或是EDC模型模拟湍流-化学反应的相互作用，在Arrhenius Rate标题下输入Arrhenius速率的以下参数：
指数前因子（方程14.1-9中的常数）。的单位必须指定为方程14.1-5中的反应速率的单位，即mol/vol*time（如kgmol/m3-s）和方程14.1-5中容积反应速率的单位，即质量/容积-时间（如kg/ m3-s）.
!注意如果你选择了英制单位系统，Arrhenius因子还是会按公制单位输入。这是因为当你使用英制单位时， FLUENT没有对你的输入使用转换因子，正确的转换因子依赖于你的，等输入。
活化能（正向速率常数表达式中的常数，方程14.1-9）。
温度指数（方程14.1-9中的常数）。
第三体效率（方程14.1-8中的。如果你有这一效率的精确值，并且希望在反应速率中包括这种影响（即在方程14.1-7中包括），打开Third Body Efficiencies选项，并点击Specify…按纽以打开Reaction Parameters面板（图14.1.5）。对于面板中的每种组分，指定Third Body Efficiencies。
！包括第三体效率不是必须的。你不需要选定Third Body Efficiencies选项，除非你有这些参数的精确值。
压力依赖反应（如果相关）如果你使用层流/有限速率，有限速率/涡耗散或是EDC模型模拟湍流-化学反应的相互作用，并且反应是压力下降反应（见14.1.1节），打开对于Arrhenius Rate的Pressure Dependent Reaction选项，并点击Specify…按纽以打开Pressure Dependent Reaction面板（图14.1.6）。
欠松弛参数，选择适当的反应类型Reaction Type（Lindemann,Troe,Sri）。这三种方法的细节见14.1.1节。然后，如果要将混合物的Bath Gas Concentration定义为一种混合物构成组分的浓度的话，你需要指定这一点，通过在下拉列表中选择适当的项。

在Reactions面板中Arrhenius Rate下指定的参数表示高压Arrhenius参数，但是你可以对Low Pressure Ahhrenius Rate下的以下参数指定值：
In（指数前因子）（方程14.1-15中的），指数前因子常常是一个非常大的数，因此你需要这一项的自然对数值。
活化能（方程14.1-15中的）
温度指数（方程14.1-15中的）
如果你选择Troe作为反应类型，你可以在Troe Parameter下指定Alpha,T1,T2,T3的值（方程14.1-22中的，，和）。如果你选择SRI反应类型, 你可以在SRI Parameter下指定a,b,c,d,e的值（方程14.1-22中的，，，和）。

6．如果你使用层流/有限速率或是EDC模型模拟湍流-化学反应的相互作用，且反应是可逆的，则打开对于Arrhenius Rate的Include Backward Reaction选项。当选定这一选项时，你将不能编辑产物的Rate Exponent，这些值将被设定为与相应的Stoich.系数相等。如果你不希望使用FLUENT的缺省值，或者你在定义你自己的反应，你将还需要指定标准状态觞和标准状态焓，以在逆向反应速率常数计算中使用（方程14.1-10）。注意可逆反应选项对于涡耗散或有限速率/涡耗散湍流-化学反应相互作用模型是不可获得的。
7．如果你使用湍流-化学反应相互作用的涡耗散或有限速率/涡耗散模型，你可以在Mixing Rate标题下输入A和B的值。但是注意除非你有可靠的数据，不要改变这些值/在大多数情况下，你只需要简单地使用缺省值。
A是湍流混合速率的常数A（方程14.1-25和14.1-26），当一种组分作为反应物在反应中出现时用于这种组分。缺省值为4.0，根据Magnussen等人给出的经验值[149]。
B是湍流混合速率的常数B（方程14.1-25和14.1-26），当一种组分作为产物在反应中出现时用于这种组分。缺省值为0.5，根据Magnussen等人给出的经验值[149]。
8．对于每一种你需要定义的反应重复步骤2-7。完成所有反应后，点OK。

定义燃料混合物的组分和反应
经常会遇到这种情况，燃烧系统中的燃料不能用一种纯组分（例如CH4或C2H6）来描述）。复杂的烃类，包括燃料油乃至木材片），很难用这种纯组分来定义。但是，如果你已经得到了这种燃料的热值和最终分析（单元组成），你可以定义一种等价的燃料组分和等价的热公式。例如，考虑一种含有50%C，6%H和44%O（按重量）的燃料。除以原子质量后，你可以得到一种“燃料”组分，其分子式为C4.17H6O2.75。你可以从一种类似的，已存在的组分开始，或者从头开始创建一种组分，并分配给它一个分子量100.04()。化学反应可以认为是：

你将需要对这一反应设定适当的stoichiometric系数。
燃料组分的燃烧热（或标准焓）可以从已知的热值计算得到，由于

其中，是1mol的标准状态焓。注意方程14.1-31的符号约定：当反应为放热时，为负。

定义混合物的物理属性
当你的FLUENT模型包括化学组分时，需要由你或数据库定义混合物材料的以下物理属性：
l 密度，可以用气体定律或作为组分比容的函数来定义
l 粘度，可以定义为组分的函数
l 热导率和比热（在涉及求解能量方程的问题中），你可以定义为组分的函数
l 标准状态焓，如果你模拟可逆反应
这些属性输入的详细介绍在第7章中提供。
！包括一步和两步的总包反应机理不可避免地忽略中间组分。在高温火焰中，忽略这些分离组分可能会导致温度的过高预测。通过增加每种组分的比热容可以得到更理想的温度场。Rose和Cooper[252]提供了一组关于温度的比热多项式函数。每种组分的比热容可按下式计算：

表14.1.1为修正的Cp多项式系数

14.1.5 定义组分边界条件
在你的模拟中，需要指定入口处每种组分的质量分数。另外，对于压力出口，你需要指定出口处的组分质量分数以在回流情况中使用。在壁面上，FLUENT将对所有组分使用0梯度（0通量）边界条件，除非你已经在壁面上定义了壁面反应（见14.2节）或是你选择指定壁面上的组分质量分数。边界条件的输入在第6章论述。
！注意你只需要明确指定前N-1种组分的质量分数。求解器通过用1减去指定组分质量分数的和来计算最后一种组分的质量分数，如果你需要明确指定最后一种组分的质量分数，你必须在列表中（Materials面板）记录这种组分，如14.1.4节所述。
进口处的扩散，使用非耦合求解器
正如14.1.1节中所提到的，当使用非耦合求解器时，没有指定入口处的组分扩散部分（因此也没有净入口输送量）。在某些情况下，你可能希望通过你的计算区域入口的只有组分的对流输送。你