第十三章 物质输送和有限速率化学反应
FLUENT可以通过求解描述每种组成物质的对流、扩散和反应源的守恒方程来模拟混合和输运，可以模拟多种同时发生的化学反应，反应可以是发生在大量相（容积反应）中，和/或是壁面、微粒的表面。包括反应或不包括反应的物质输运模拟能力，以及当使用这一模型时的输入将在本章中叙述。
注意你可能还希望使用混合物成分的方法（对非预混系统，在14章介绍）、反应进程变量的方法（对预混系统，在15章介绍），或部分预混方法（在16章介绍）来模拟你的反应系统。见12章FLUENT中反应模拟方法的概述。
本章中的分为以下章节：
13.1 容积反应 13.2 壁面表面反应和化学蒸汽沉积 13.3 微粒表面反应 13.4 无反应物质输运
13．1 容积反应
与容积反应有关的物质输运和有限速率化学反应方面的信息在以下小节中给出：
13.1.1 理论
13.1.2 模拟物质输运和反应的用户输入概述 13.1.3 使能物质输运和反应，并选择混合物材料 13.1.4 混合物和构成物质的属性定义 13.1.5 定义物质的边界条件 13.1.6 定义化学物质的其他源项
13.1.7 化学混合和有限速率化学反应的求解过程 13.1.8 物质计算的后处理
13.1.9 从CHEMKIN导入一个化学反应机理
13.1.1 理论 物质输运方程
当你选择解化学物质的守恒方程时，FLUENT通过第i种物质的对流扩散方程预估每种物质的质量分数，Yi。守恒方程采用以下的通用形式：
YivYiJiRiSi （13.1-1） t其中Ri是化学反应的净产生速率（在本节稍后解释），Si为离散相及用户定义的源项导致的额外产生速率。在系统中出现N种物质时，需要解N-1个这种形式的方程。由于质量分数的和必须为1，第N种物质的分数通过1减去N-1个已解得的质量分数得到。为了使数值误差最小，第N种物质必须选择质量分数最大的物质，比如氧化物是空气时的N2。 层流中的质量扩散
在方程13.1-1中，Ji是物质i的扩散通量，由浓度梯度产生。缺省时，FLUENT使用稀释近似，这样扩散通量可记为：
JiDi,mYi （13.1.2）
WORD完美格式 下载可编辑
专业资料
这里Di,m是混合物中第i种物质的扩散系数。
对于确定的层流流动，稀释近似可能是不能接受的，需要完整的多组分扩散。在这些例子中，可以解Maxwell-Stefan方程，详细情况见7.7.2节。 湍流中的质量扩散
在湍流中，FLUENT以如下形式计算质量扩散：
tJiDi,mScYi (13.1.3)
t其中Sct是湍流施密特数，
t（缺省设置值为0.7）。 Dt注意，湍流扩散一般淹没层流扩散，在湍流中指定详细的层流性质是不允许的。 能量方程中的物质输送处理
在许多多组分混合流动中，物质扩散导致了焓的传递。
nhiJi i1这种扩散对于焓场有重要影响，不能被忽略。特别是，当所有物质的Lewis数
Leik （13.1-4）
cpDi,m远离1时，忽略这一项会导致严重的误差。
FLUENT缺省地包含这一项。在方程13.1-4中，k为热导率。 进口处的扩散
在FLUENT的非耦合求解器中，入口的物质净输送量由对流量和扩散量组成，对耦合解算器，只包括对流部分。对流部分由你指定的物质浓度确定。扩散部分依赖于计算得到的物质浓度场。因此，扩散部分（从而使净输送量）不预先指定。见13.1.5节有关指定入口净输送量的信息。 反应建模的一般有限速率形式
反应速率作为源项在方程13.1-1中出现，在FLUENT中根据以下三种模型中的一个计算： 层流有限速率模型：忽略湍流脉动的影响，反应速率根据Arrhenius公式确定。
涡耗散模型：认为反应速率由湍流控制，因此避开了代价高昂的Arrhenius化学动力学计算。
涡耗散概念（EDC）模型：细致的Arrhenius化学动力学在湍流火焰中合并。注意详尽的化学动力学计算代价高昂。
通用有限速率对于范围很广的应用，包括层流或湍流反应系统，预混、非预混、部分预混燃烧系统都适用。
层流有限速率模型
层流有限速率模型使用Arrhenius公式计算化学源项，忽略湍流脉动的影响。这一模型对于层流火焰是准确的，但在湍流火焰中Arrhenius化学动力学的高度非线性，这一模型一般不精确。对于化学反应相对缓慢、湍流脉动较小的燃烧，如超音速火焰可能是可以接受的。
化学物质i的化学反应净源项通过有其参加的NR个化学反应的Arrhenius反应源的和计算得到。
ˆ RiMw,iRi,ri1Nr WORD完美格式 下载可编辑
专业资料
ˆ为第i种物质在第r个反应中的产生/分解速率。反应可能发生其中Mw,i是第i种物质的分子量，Ri,r在连续相反应的连续相之间，或是在表面沉积的壁面处，或是发生在一种连续相物质的演化中。
考虑以如下形式写出的第r个反应：
i1N'i,rMii'',rMi （13.1-6） kb,ri1N其中N——系统中化学物质数目；
i',r——反应r中反应物i的化学计量系数; i'',r——反应r中生成物i的化学计量系数;
Mi——第i种物质的符号；
kf,r——反应r的正向速率常数； kb,r——反应r的逆向速率常数；
方程13.1-6对于可逆和不可逆反应（FLUENT中缺省为不可逆）都适用。对于不可逆反应，逆向速率常数kb,r简单地被忽略。
方程13.1-6中的和是针对系统中的所有物质，但只有作为反应物或生成物出现的物质才有非零的化学
计量系数。因此，不涉及到的物质将从方程中清除。
反应r中物质i的产生/分解摩尔速度以如下公式给出：
ˆ'''Ri,ri,ri,rj,rj,rNrNrkf,rCj,rkb,rCj,r （13.1-7）
j1j1'''其中：Nr——反应r 的化学物质数目；
Cj,r——反应r中每种反应物或生成物j的摩尔浓度；
'j,r——反应r中每种反应物或生成物j的正向反应速度指数； 'j',r——反应r中每种反应物或生成物j的逆向反应速度指数；
见13.1.4节有关输入整体正向反应（不可逆）和单元反应（可逆）的化学计量系数和速率指数方面的内容。
表示第三体对反应速率的净影响。这一项由下式给出：
其中j,r为第r个反应中第j种物质的第三体影响。在缺省状态，FLUENT在反应速率计算中不包
WORD完美格式 下载可编辑
专业资料
括第三体影响。但是当你有它们的数据时，你可以选择包括第三体影响。
反应r的前向速率常数kf,r通过Arrhenius公式计算：
其中，Ar——指数前因子（恒定单位）；
r——温度指数（无量纲）；
Er——反应活化能(J/kmol);
R——气体常数（J/kmolK）
你（或者数据库）可以在FLUENT的问题定义中提供i,r,i,r,j,r,j,r,r,Ar,Er，并可选择提
供j,r。
如果反应是可逆的，逆向反应常数kb,r可以根据以下关系从正向反应常数计算：
其中Kr为平衡常数，从下式计算：
其中patm表示大气压力（101325Pa）。指数函数中的项表示Gibbs自由能的变化，其各部分按下式计算：
其中Si0和hi0是标准状态的熵和标准状态的焓（生成热）。这些值在FLUENT中作为混合物材料的
属性指定。
压力独立反应
FLUENT可以用以下三种方法之一来表示压力独立反应（或压力下降）反应的速率表达式。“压力下降”
WORD完美格式 下载可编辑
专业资料
反应是发生在Arrhenius 高压和低压限制之间的反应，因而不仅仅依赖于温度。
有三种方法表示在“fall-off”区域的速率表达式，最简单的是Lindemann[140]形式。还有其它良种相关的方法，Troe方法[77]和SRI方法[230]，它们提供了更精确的”fall-off”区域表达形式。
Arrhenius速率参数对于高压和低压限制都是需要的。两个限制的速率系数融合以产生光滑的压力独立表达式。在Arrhenius形式中，高压限制k和低压限制klow的参数如下：
其中pr定义为：
在任意压力下，净反应速率常数为：
M为溶液气体的浓度，可以包括第三体效率。如果方程13.1-16函数F为1，则是Lindemann形式。
FLUENT提供了两种其他形式来表述F，称为Troe方法和SRI方法。
在Troe方法中，F按下式给出：
其中，
参数,T3,T2,T1做为输入确定。 在SRI方法中，缝合函数F近似为
除了低压限制表达式中的三个Arrhenius参数以外，你还需要提供F表达式中的a, b, c, d, e。 ！化学动力学机理中有很高的非线性并且形成了一组强烈耦合的方程。求解过程指导见13.1.7
WORD完美格式 下载可编辑
专业资料
节。如果你有一个CHEMKIN形式的化学反应机理[112]，你可以将这一机理导入FLUENT，如13.1.9节。
涡耗散模型
大部分燃料快速燃烧。整体反应速率由湍流混合控制。在非预混火焰中，湍流缓慢地通过对流/混合燃料和氧化剂进入反应区，在反应区它们快速地燃烧。在预混火焰中，湍流对流/混合冷的反应物和热的生成物进入反应区，在反应区迅速地发生反应。在这些情况下，燃烧称为混合限制的，复杂，常常是未知的化学反应动力学速率可以安全地忽略掉。
FLUENT提供了湍流-化学反应相互作用模型，基于Magnussen 和Hjertager[149]的工作，称为涡耗散模型。
反应r中物质i的产生速率Ri,r由下面两个表达式中较小的一个给出：
在方程13.1-25和13.1-26中，化学反应速率由大涡混合时间尺度k/控制，如同Splading[227]的涡破碎模型一样。只要湍流出现（k/0），燃烧即可进行，不需要点火源来启动燃烧。这通常对于非预混火焰是可接受的，但在预混火焰中，反应物一进入计算区域（火焰稳定器上游）就开始燃烧。为了修正这一点，FLUENT提供了有限速率/涡耗散模型，其中Arrhenius（方程13.1-7）和涡耗散(方程13.1-25和13.1-26)反应速率都进行计算。净反应速率取两个速率中较小的。实际上，Arrhenius反应速率作为一种动力学开关，阻止反应在火焰稳定器之前发生。一旦火焰被点燃，涡耗散速率通常会小于Arrhenius反应速率，并且反应是混合限制的。
！尽管FLUENT允许采用涡耗散模型和有限速率/涡耗散模型的多步反应机理（反应数>2），但可能会产生不正确的结果。原因是多步反应机理基于Arrhenius速率，每个反应的都不一样。在涡耗散模型中，每个反应都有同样的湍流速率，因而模型只能用于单步（反应物—产物）或是双步（反应物—中间产物，中间产物—产物）整体反应。模型不能预测化学动力学控制的物质，如活性物质。为合并湍流流动中的多步化学动力学机理，使用EDC模型（下面介绍）。
！涡耗散模型需要产物来启动反应（见方程13.1-26）。当你初始化求解的时候，FLUENT设置产物的质量比率为0.01，通常足够启动反应。但是，如果你首先聚合一个混合解，其中所有的产物质量比率都为0，你可能必须在反应区域中补入产物以启动反应。详细内容见13.1.7节。
LES的涡耗散模型
当使用LES湍流模型时，湍流混合速率（方程13.1-25和13.1-26中的/k）被亚网格尺度混合速率替代。计算为：
WORD完美格式 下载可编辑
专业资料
涡-耗散-概念（EDC）模型
涡-耗散-概念（EDC）模型是涡耗散模型的扩展，以在湍流流动中包括详细的化学反应机理[148]。它假定反应发生在小的湍流结构中，称为良好尺度。良好尺度的容积比率按下式模拟[80]：
其中*表示良好尺度数量，
C——容积比率常数=2.1377；
——运动粘度
认为物质在好的结构中，经过一个时间尺度
后开始反应。
其中C为时间尺度常数，等于0.4082
在FLUENT中，良好尺度中的燃烧视为发生在定压反应器中，初始条件取为单元中当前的物质和温度。反应经过时间尺度后开始进行，由方程13.1-7的Arrhenius速率控制，并且用普通微分方程求解器CVODE进行数值积分[45]。经过一个时间的反应后物质状态记为Yi 物质i的守恒方程13.1-1中的源项计算公式为：
**
EDC模型能在湍流反应流动中合并详细的化学反应机理。但是，典型的机理具有不同的刚性，它们的数值积分计算开销很大。因而，只有在快速化学反应假定无效的情况下才能使用这一模型，例如在快速熄灭火焰中缓慢的CO烧尽、在选择性非催化还原中的NO转化。
推荐使用双精度求解器以避免刚性机理中固有的大指数前因子和活化能产生的舍入误差。见13.1.7节获得使用EDC模型求解的指导。
13.1.2 物质输送和反应模拟输入概览
设定涉及物质输送和反应问题的基本步骤如下，每一步的详细执行过程见13.1.3-13.1.5节。有关设定和求解的附加信息在13.1.6-13.1.8节中提供。
1． 选定物质输送和容积反应，指定混合物材料。见13.1.3节（混合物材料概念在下面解释）。 2． 如果你还要模拟壁面或微粒表面反应，则要打开壁面和/或微粒表面反应。细节见13.2和13.3
节。
3． 检查和/或定义混合物的属性。（见13.1.4节）。混合物属性包括：
混合物中的物质 反应
其他物理属性（如粘度、比热）
WORD完美格式 下载可编辑
专业资料
4． 检查和/或设置混合物中单个物质的属性（见13.1.4节） 5． 设置物质边界条件（见13.1.5节）
在很多情况下，当你选择混合物材料是，求解器从材料数据库中得到物质性质、反应等，因而你将不需要修改任何物理属性。但有一些性质可能在数据库中没有定义。如果有任何性质需要设置时，你将被警告，这样你可以指定这些性质的适当值。你还可能希望检查数据库中这些性质的值，以确定它们对你的应用是否正确。修改已存在的混合物材料或从最开始创建一个新的材料的详细内容见13.1.4节。混合物材料的修改包括以下方面：
物质的添加和删除 改变化学反应
修改混合物的其他材料属性 修改混合物本构物质的材料属性
如果你在求解一个反应流问题，你常常希望将混合物的比热定义为组成的函数，将每种物质的比热定义为温度的函数。你还可能对其它一些性质希望也做这样的定义。缺省状态下，将使用恒定的属性，但对一些物质的性质，在数据库中存在一个温度的分段多项式函数可供你使用。如果你知道更多适合于你的问题的函数的话，你还可以选择指定一个不同的温度依赖函数。
混合物材料
在FLUENT中提出混合物材料的概念以方便物质输送和反应流动的设置。混合物材料可以认为是一组物质和一列控制它们相互作用的规律。混合物材料带有以下性质：
一列本构物质，相对于“流体”物质
一列混合定律，指示如果希望得到组分依赖的属性，混合属性（密度，粘度，比热等）如何从单
个物质得到
如果希望属性不依赖组分，直接指定混合物属性
其它与耽搁物质无关的材料属性（如吸收和辐射系数）
一组反应，包括反应类型（有限速率，涡耗散等）和化学计量和速率常数
混合物材料和流体材料都储存在FLUENT的材料数据库中。包括许多常见的混合物材料（如甲烷-空气，丙烷-空气）。通常，在数据库中定义了一步/两步反应机理和大量混合物及其构成物质的属性。当你指定了你希望使用哪种混合物材料后，适当的混合物材料，流体材料和属性将被装载到求解器中。如果缺少任何所选材料（或构成流体材料）必须的属性，求解器将通知你需要指定它。另外，你可以选择修改任何预定义的属性。见7.1.2节了解有关FLUENT数据库属性数据源的信息。
例如，如果你计划模拟一种甲烷-空气的燃烧，你不需要明确指定反应中涉及的物质和反应本身。只需要简单地选择甲烷-空气作为使用的混合物材料，相关的物质（CH4，O2，CO2，H2O和N2）和反应数据将从数据库装入求解器。然后你可以检查物质、反应和其它属性并定义其它任何缺少的属性，和/或修改任何你希望使用不同值或函数的属性。通常你希望定义一个与组分、温度相关的比热，还可能希望将其它属性定义为温度和/或组分的函数。
混合物材料的使用给你提供了一种灵活性，可以使用大量预定义混合物中的一种，修改这些混合物，或是创建你自己的混合物材料。自定义混合物材料在Materials面板中进行，在13.1.4节中讲述。
13.1.3 选定物质输送和反应，并选择混合物材料
物质输送和容积反应的问题设置总物质模型（Species Model）面板开始（图13.1.1）。
WORD完美格式 下载可编辑
专业资料
1． 在Model下，选择Species Transport.
2． 在Reaction下，选择Volumetric reactions
3． 在Mixture Properties下的Mixture Material下拉列表中选择在你的问题中希望使用
的混合物材料
下拉列表中将包括所有在当前数据库中定义的混合物。为检查一种混合物材料的属性，选择它，并点击View…按纽。如果你所希望使用的混合物不在列表中，选择混合物模板（mixture-template）材料，并参阅13.1.4节了解设置你自己的混合物属性的详细内容。如果有一种混合物材料和你你所希望使用的混合物相似，可以选择这一材料并参阅13.1.4节修改已存在材料性质的详细信息。
当选择Mixture Material时，混合物中的Number of Volumetric Species将在面板中显示，表达你的信息。
!注意如果你在已经选定物质输送后，重新打开Species Model面板时，只有你的案例中可得到的混合物材料才会显示在列表中。你可以通过从数据库中拷贝在你的案例中增加更多混合物材料，如7.1.2节所述，或是创建一个新的混合物，如7.1.2节和13.1.4节所述。
正如在13.1.2节中提到的，物质输送的模拟参数和反应（如果有关）将自动从数据库中装入。如果缺少任何信息，当你点击Species Model面板中的OK按纽后将被告知缺少什么。如果你希望检查或修改混合物材料的任何属性，你将使用Materials面板，如13.1.4节所述。
4． 选择湍流-化学反应相互作用模型，可以使用四种模型： 层流有限速率：只计算Arrhenius速率（见方程13.1-7），并忽略湍流-化学反应相互作用。 涡耗散模型（针对湍流流动）： 只计算混合速率（方程13.1-25和方程13.1-26）。 有限速率/涡耗散模型（针对湍流流动）：计算Arrhenius速率和混合速率，并使用其中较小的一个。 EDC模型（湍流流动）：使用详细的化学反应机理模拟湍流-化学反应相互作用（见方程13.1-25和13.1-26）。
5． 如果你选择EDC模型，你可以选择修改容积比率常数和时间尺度常数（方程13.1-28中的C和方程13.1-29中的C，尽管通常推荐缺省值。此外，为减少化学反应计算的开销，你可以增加每次化学反应更新的流动迭代（Flow Iteration Per Chemistry Update）次数。缺省时，FLUENT每十次流动迭代更新化学反应一次。
WORD完美格式 下载可编辑
专业资料
6． （可选）如果你希望模拟完整的多组分扩散或热扩散，打开完整多组分扩散或热扩散Full Multicomponent Diffusion 或Thermal Diffusion选项。
13.1.4 定义混合物的属性和构成物质
如13.1.2节所讨论的，如果你使用来自数据库的混合物材料，大部分混合物和物质属性已经定义了。你可以跟随这一节的过程检查当前的属性、修改某些属性或是设定一个你从头开始定义的全新的混合物材料的所有属性。
记住你将需要定义混合物材料和其构成物质的属性。由于物质属性输入可能依赖于你定义混合物属性的方法，在设定构成物质的属性之前定义混合物属性非常重要。建议按照如下的属性输入顺序：
1． 定义混合物物质和反应，定义混合物物理属性。记住在设定混合物材料属性时点击
Change/Create按纽。
2． 定义混合物中物质的物理属性。记住在设定了每种物质的属性后点击Change/Create按纽。 所有这些步骤都在Materials面板中进行，在本节中将详细叙述。
定义混合物中的物质
如果你使用数据库中的混合物材料，混合物中的物质已经为你定义了。如果你创建你自己的材料或是修改已存在材料中的物质，你将需要自己定义它们。 在Materials面板中（图13.1.2），检查材料类型Material Type是否已经设置为混合物，并且你的混合物是否已经在混合物材料列表Mixture Materials list中选定。点击Mixture Species右边的Edit…按纽打开Species面板（图13.1.3）。
Species面板概览
在Species面板中，已选物质Selected Species列表显示所有混合物中的流体相物质。如果你模拟壁面或微粒表面反应，已选物质Selected Species列表将显示所有混合物中的表面物质。表面物质是那些从壁面边界或是离散相微粒(如Si(s))产生或散发出来的，以及在流体相物质中不存在的物质。表面物质和壁面反应将在13.2节中叙述，微粒表面反应的有关内容见13.3节。 ！已选物质Selected Species列表中物质的顺序非常重要。FLUENT认为列表中最后的物质是大量的物质。因此，当你从混合物材料中增加或是删除物质时，必须小心将最丰富（按质量）的物质作为最后一个物质。
可获得的物质Available Material列表显示可获得，但不在混合物中的材料。通常你可以在列表中看到空气air，因为缺省时，空气通常是可获得的。
WORD完美格式 下载可编辑
专业资料
在混合物中增加物质
如果你从头创建了一个混合物，或者从一种已存在的混合物开始并增加一些缺少的物质，你首先需要从数据库中导入希望的物质（或是创建它们，如果它们不在数据库中出现的话），以便在求解器中可以得到，增加物质的过程在下面列出。（在开始前你需要关闭Species面板，因为它是一个模式面板，不允许在它打开的时候做别的事）
1． 在Materials面板中，点击Datebase…按纽以打开Database Materials面板，并拷贝希望的
WORD完美格式 下载可编辑
专业资料
物质，如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)列表的末尾。
4． 对所有希望的物质重复上述步骤。结束后点击OK按纽。
！ 增加一种物质到列表中将改变物质的顺序。你必须确保列表中的最后一个物质是大量的物质，并且你需要检查你以前所设定的所有边界条件，欠松弛因子或其它求解参数，如详细叙述的那样。
在混合物中去除物质
为从混合物中去除一种物质，只需要简单地在Selected Species列表（或是Selected Surface Species列表）中选定它，并点击列表下的Remove按纽。这种物质将从列表中去除并增加到Available Materials列表中。
！ 去除列表中的一种物质将改变物质的顺序。你必须确保列表中的最后一个物质是大量的物质，并且你需要检查你以前所设定的所有边界条件，欠松弛因子或其它求解参数，如详细叙述的那样。
重排物质
如果你发现Selected Species列表中最后一种物质不是最丰富的物质（应该是），你需要重排物质以得到正确的顺序。
1． 将最大量的物质从Selected Species列表中去除。它将出现在Available Materials列表中。 2． 再次增添这种物质，它将自动放在列表的末尾。
物质的命名和顺序
如上面讨论的，你必须在增加或去除物质时保持最丰富的物质作为Selected Species列表中的最后一种物质。在增加或去除物质时还需要知道的一些考虑事项在这里给出。
一种物质有三个特征在供求解器识别：名字，化学分子式和在Species面板中物质列表中的位置。改变这些特征将会产生以下效果：
你可以改变一种物质的名字（使用Materials面板，如7.1.2节所述），而不产生任何影响。 你不能改变一种物质的给定的化学分子式。
如果你增加或去除物质，将改变物质列表的顺序。这时，所有的边界条件，求解器参数和物质的
求解数据将被从设为缺省值。（其它流动变量的求解数据，边界条件，求解器参数将不受影响）因此，如果你增加或去除物质，你需要注意这一新定义问题的物质边界条件和求解参数。另外，你必须认识到基于原来物质顺序给出的物质浓度或是存储在数据文件中的浓度将会与新定义的问题不匹配。你可以使用数据文件作为初始猜测值，但你必须知道数据文件中的物质浓度将可能对新定义的模型提供一种不好的初始猜测。
定义反应
如果你的FLUENT模型中涉及化学反应，你可以接着定义参与的已定义物质的反应。这只有在你从头开始创建一种混合物、修改了物质或是出于某些其他原因希望重定义反应时才是必须的。
WORD完美格式 下载可编辑
专业资料
在Materials面板的Reaction下拉列表中显示适当的反应机理，依赖于你在Species Model面板中选择的湍流-化学反应相互作用模型（见13.1.3节）。如果你使用层流有限速率或EDC模型，反应机理将是有限速率的，如果你使用涡耗散模型，反应机理将是涡耗散的；如果使用有限速率/涡耗散模型，反应机理将是有限速率/涡耗散的。
反应定义的输入
为定义反应，点击Reaction右侧的Edit…按纽。 将打开Reaction面板（图13.1.4）。 定义反应的步骤如下： 1． 在Total Number of Reaction区域中设定反应数目（容积反应，壁面反应和微粒表面反应）
（使用箭头改变数值，或是键入值并按回车键。
注意如果你的模型包括离散相的燃烧微粒，只有在你计划使用表面燃烧的多表面反应模型时，才必须在反应数目中包括部分表面反应（s）（如碳的燃烧，多样碳粒氧化）
2. 设定你希望定义的反应的Reaction ID
2. 如果是流体相反应，保持缺省选项Volumetric作为反应类型。如果是壁面反应（在13.2节
中描述）或者颗粒表面反应（13.3节描述），选择Wall Surface或Particle reaction 作为反应类型。有关定义颗粒表面反应的进一步信息见13.3.2节。
3． 通过增加Number of Reactants和Number of Products的值指定反应中涉及的反应物和生成物的
WORD完美格式 下载可编辑
专业资料
数量。在Species下拉列表中选择每一种反应物或生成物，然后在Stoich. Coefficient和Rate Exponent区域中设定它的化学计量系数和速率指数。（化学计量系数是方程13.1-6中的常数i,r和
i,r，速率指数是方程13.1-7中的常数j,r和j,r。共有两种普通类型的反应可以在Reactions
面板中处理。因此正确输入每种反应的参数非常重要。反应的类型如下：
整体正向反应（无逆向反应）：产物一般不影响正向速率，因此所有产物的速率指数（j,r）
应该为0。对于反应物，设定速率指数（j,r）为期望的值。如果某种反应不是基元反应，速率指数一般不等于这种物质的化学计量系数。整体正向反应的一个实例是甲烷的燃烧：
图13.1.4表示了甲烷燃烧的系数输入（还可在DataBase Materials面板中查阅甲烷/空气混合物材料）。 注意：在某些情况下，你可能希望模拟产物影响正向速率的反应。对于这些情况，设定产物速率指数（j,r）为期望的值。这种反应的一个例子是气体转换反应（见DataBase Materials面板中的CO/空气混合物材料），其中水的存在对反应速率有影响：
在气体转换反应中，速率表达式可以定义为：
可逆反应：假定每种物质的化学计量系数等于速率指数，这种反应的一个例子是SO2氧化为SO3：
参见下面的第6步如何选定可逆反应。
4． 如果你使用 层流/有限速率，有限速率/涡耗散或是EDC模型模拟湍流-化学反应的相互作用，在
Arrhenius Rate标题下输入Arrhenius速率的以下参数：
ˆ的单位，指数前因子（方程13.1-9中的常数Ar）。Ar的单位必须指定为方程13.1-5中的反应速率Ri,r即mol/vol*time（如kgmol/m-s）和方程13.1-5中容积反应速率的单位，即质量/容积-时间（如kg/
3
m-s）.
!注意如果你选择了英制单位系统，Arrhenius因子还是会按公制单位输入。这是因为当你使用英制单位时， FLUENT没有对你的Ar输入使用转换因子，正确的转换因子依赖于你的i,r，r等输入。
WORD完美格式 下载可编辑
3
专业资料
活化能（正向速率常数表达式中的常数Er，方程13.1-9）。 温度指数（方程13.1-9中的常数r）。 第三体效率（方程13.1-8中的j,r。如果你有这一效率的精确值，并且希望在反应速率中包括这种影
响（即在方程13.1-7中包括），打开Third Body Efficiencies选项，并点击Specify…按纽以打开Reaction Parameters面板（图13.1.5）。对于面板中的每种物质，指定Third Body Efficiencies。 ！包括第三体效率不是必须的。你不需要选定Third Body Efficiencies选项，除非你有这些参数的精确值。
压力依赖反应（如果相关）如果你使用层流/有限速率，有限速率/涡耗散或是EDC模型模拟湍流-化学反应的相互作用，并且反应是压力下降反应（见13.1.1节），打开对于Arrhenius Rate的Pressure Dependent Reaction选项，并点击Specify…按纽以打开Pressure Dependent Reaction面板（图13.1.6）。 欠松弛参数，选择适当的反应类型Reaction Type（Lindemann,Troe,Sri）。这三种方法的细节见13.1.1节。然后，如果要将混合物的Bath Gas Concentration定义为一种混合物构成物质的浓度的话，你需要指定这一点，通过在下拉列表中选择适当的项。
在Reactions面板中Arrhenius Rate下指定的参数表示高压Arrhenius参数，但是你可以对Low Pressure Ahhrenius Rate下的以下参数指定值：
In（指数前因子）（方程13.1-15中的Alow），指数前因子Alow常常是一个非常大的数，因此你需要这一项的自然对数值。
活化能（方程13.1-15中的Elow）
WORD完美格式 下载可编辑
专业资料
温度指数（方程13.1-15中的low）
如果你选择Troe作为反应类型，你可以在Troe Parameter下指定Alpha,T1,T2,T3的值（方程13.1-22中的，T1，T2和T3）。如果你选择SRI反应类型, 你可以在SRI Parameter下指定a,b,c,d,e的值（方程13.1-22中的a，b，c，d和e）。
6．如果你使用层流/有限速率或是EDC模型模拟湍流-化学反应的相互作用，且反应是可逆的，则打开对于Arrhenius Rate的Include Backward Reaction选项。当选定这一选项时，你将不能编辑产物的Rate Exponent，这些值将被设定为与相应的Stoich.系数相等。如果你不希望使用FLUENT的缺省值，或者你在定义你自己的反应，你将还需要指定标准状态觞和标准状态焓，以在逆向反应速率常数计算中使用（方程13.1-10）。注意可逆反应选项对于涡耗散或有限速率/涡耗散湍流-化学反应相互作用模型是不可获得的。 7．如果你使用湍流-化学反应相互作用的涡耗散或有限速率/涡耗散模型，你可以在Mixing Rate标题下输入A和B的值。但是注意除非你有可靠的数据，不要改变这些值/在大多数情况下，你只需要简单地使用缺省值。
A是湍流混合速率的常数A（方程13.1-25和13.1-26），当一种物质作为反应物在反应中出现时用于这
WORD完美格式 下载可编辑
专业资料
种物质。缺省值为4.0，根据Magnussen等人给出的经验值[149]。 B是湍流混合速率的常数B（方程13.1-25和13.1-26），当一种物质作为产物在反应中出现时用于这种物质。缺省值为0.5，根据Magnussen等人给出的经验值[149]。
8．对于每一种你需要定义的反应重复步骤2-7。完成所有反应后，点OK。
定义燃料混合物的物质和反应 经常会遇到这种情况，燃烧系统中的燃料不能用一种纯物质（例如CH4或C2H6）来描述）。复杂的烃类，包括燃料油乃至木材片），很难用这种纯物质来定义。但是，如果你已经得到了这种燃料的热值和最终分析（单元组成），你可以定义一种等价的燃料物质和等价的热公式。例如，考虑一种含有50%C，6%H和44%O（按重量）的燃料。除以原子质量后，你可以得到一种“燃料”物质，其分子式为C4.17H6O2.75。你可以从一种类似的，已存在的物质开始，或者从头开始创建一种物质，并分配给它一个分子量100.04(4.1712612.7516)。化学反应可以认为是：
你将需要对这一反应设定适当的stoichiometric系数。
燃料物质的燃烧热（或标准焓）可以从已知的热值H计算得到，由于
其中，hi0是1mol的标准状态焓。注意方程13.1-31的符号约定：当反应为放热时，H为负。
定义混合物的物理属性
当你的FLUENT模型包括化学物质时，需要由你或数据库定义混合物材料的以下物理属性： 密度，可以用气体定律或作为组分比容的函数来定义 粘度，可以定义为组分的函数
热导率和比热（在涉及求解能量方程的问题中），你可以定义为组分的函数 标准状态焓，如果你模拟可逆反应 这些属性输入的详细介绍在第7章中提供。
！包括一步和两步的总包反应机理不可避免地忽略中间物质。在高温火焰中，忽略这些分离物质可能会导致温度的过高预测。通过增加每种物质的比热容可以得到更理想的温度场。Rose和Cooper[252]提供了一组关于温度的比热多项式函数。每种物质的比热容可按下式计算：
表13.1.1为修正的Cp多项式系数
WORD完美格式 下载可编辑
专业资料
13.1.5 定义物质边界条件
在你的模拟中，需要指定入口处每种物质的质量分数。另外，对于压力出口，你需要指定出口处的物质质量分数以在回流情况中使用。在壁面上，FLUENT将对所有物质使用0梯度（0通量）边界条件，除非你已经在壁面上定义了表面反应（见13.2节）或是你选择指定壁面上的物质质量分数。边界条件的输入在第6章论述。
！注意你只需要明确指定前N-1种物质的质量分数。求解器通过用1减去指定物质质量分数的和来计算最后一种物质的质量分数，如果你需要明确指定最后一种物质的质量分数，你必须在列表中（Materials面板）记录这种物质，如13.1.4节所述。 进口处的扩散，使用非耦合求解器
正如13.1.1节中所提到的，当使用非耦合求解器时，没有指定入口处的物质扩散部分（因此也没有净入口输送量）。在某些情况下，你可能希望通过你的计算区域入口的只有物质的对流输送。你可以通过取消进口物质扩散做到这一点。在缺省状态下，FLUENT在入口包括物质的扩散通量。为关闭入口扩散，使用define/models/species-transport/inlet-diffusion? Text命令。 13.1.6 定义化学物质的其他源项
你可以通过在Fluid面板中定义一个源项来在计算区域中定义一个化学物质的源或是容器。当你的问题中存在物质源，但你又不希望通过化学反应机理来模拟它的时候，可以选择这一方法。6.27节叙述了在你的FLUENT模型中定义物质源所需要遵从的步骤。如果源项不是常数，你可以使用用户定义函数。用户定义函数的细节见单独的UDF手册。
13.1.7 化学混合和有限速率反应的求解步骤
尽管许多涉及化学物质的模拟在求解过程中不需要特殊的步骤，你可能发现本节中提供的一种或多种求解技术会对加速收敛或提高更复杂模拟的稳定性有所帮助。如果你的问题涉及许多物质和/或化学反应，尤其是模拟燃烧流动时，以下列出的技术可能特别重要。 反应流中的稳定性和收敛
在反应流中获得收敛解非常困难，有很多原因。首先，化学反应对基本流型的影响可能非常强烈，导致模型中质量/动量平衡和物质输运方程的强烈耦合。在燃烧中，反应导致大的热量释放和相应的密度变化以及流动中很大的加速度，上述耦合尤其明显。但是，当流动属性依赖于物质浓度时，所有的反应系统都具有一定程度的耦合。处理这些耦合问题的最好方法是使用下面介绍的两步求解过程及使用22.9节中
WORD完美格式 下载可编辑
专业资料
介绍的欠松弛方法。
反应流中的第二个收敛问题涉及反应源项的强度。当你的FLUENT模型涉及非常快的反应速率（即比对流和扩散速率快得多），物质输运方程的求解在数值上非常困难。这种系统称为“刚性”系统，当你定义涉及非常快的动力速度的模型，尤其是这些速度描述可逆反应或竞争反应，这种系统得以创建。在涡耗散模型中，较慢的湍流速率去除了非常快的反应速率。对非预混系统，反应速率从模型中去除。对于层流化学反应的刚性系统，推荐使用耦合求解器代替非耦合求解器。对湍流有限速率机理（可能是刚性的），推荐使用EDC模型，这一模型对化学反应使用一个刚性的ODE积分器。求解刚性化学反应系统的其它指南见下述内容。
两步求解过程（冷流动模拟）
将一个反应流动作为两步过程求解对于获得你的FLUENT问题的稳定收敛解是一个实用的方法。在这一过程中，你从求解不带反应的流动、能量和物质方程（“冷流动”，或无反应流动）开始。当建立基本的流型后，你可以再选择反应，并重新开始计算。冷流动求解提供了燃烧系统计算的初始解。这种燃烧模拟的两步方法可以采用以下步骤完成：
1． 设定包括所有感兴趣物质和反应的问题
2． 通过关闭Species Model面板中的Volumetric Reactions选项暂时不选择反应计算
3． 关闭Solution Controls面板中的产物计算
4． 计算初始解（冷流动）。（注意通常获得完全收敛的冷流动解没有实际价值，除非你对无反应解也有兴趣）
5． 打开Species Model面板中的Volumetric Reactions选项使能反应计算
6． 打开所有反应。如果你使用层流有限速率、有限速率/涡耗散，或是EDC模型模拟湍流-化学反应相互作用，你可能需要增添一个点火源（如下所述）。 密度欠松弛
燃烧模拟难以收敛的一个主要原因是温度的剧烈变化引起密度的剧烈变化，从而导师流动求解的不稳定性。当你使用非耦合求解器时，FLUENT允许你欠松弛密度的这种变化以降低收敛的困难。密度欠松弛因子的缺省值为1，如果你遇到收敛问题，你可以将这个值减少到0.5到1之间（在Solution Controls面板中）。
燃烧模拟的点火
如果你将燃料引入氧化剂，自发的点火不会发生，除非混合物的温度超过了维持燃烧所需要的活化能阈值。这一物理问题在FLUENT中也会出现。如果你使用层流有限速率、有限速率/涡耗散或EDC模型模拟湍流-化学反应相互作用，你将不得不提供一个点火源以启动燃烧。这个点火源可以是加热的表面或温度超过点火温度的入口质量流。但是，这常常等同于一个火花：一个初始求解状态使得燃烧可以进行。你可以通过在FLUENT模型中一个包含有足够燃料/空气混合物以使点火能发生的区域给一个高的温度，来提供这个初始火花。
根据模型的不同，你可能需要提供温度和燃料/氧化剂/产物浓度以在你的模型中产生点火。点火。这种补缀对于最终的稳态解没有影响——不超过火柴的位置对它点燃的火炬最终流型的影响。缀入初始值的细节见22.13.2节。
刚性层流化学反应系统的求解
当使用层流有限速率模型模拟层流反应系统时，你可能需要在反应机理是刚性的时候使用耦合求解器。
WORD完美格式 下载可编辑
专业资料
（反应流中的求解收敛 问题见以上讨论。注意你可以对湍流火焰使用层流有限速率模型，这表示忽略湍流-化学反应的相互作用。）另外，你可以通过使用stiff-chemistry 文本命令为耦合求解器提供进一步的求解稳定性。
这一选项对于大Courant(CFL)数特例也适合，尽管需要额外的运算来计算化学雅可比行列式的特征值[258]。当你使能刚性化学反应求解器时，你将会被询问指定下列参数：
温度的正性速率限制：以这个因子乘旧的温度限制新的温度改变。缺省值为0.2。 温度时间步长减少因子：当温度变化过快时限制当地CFL数。缺省值为0.25。 最大允许时间步长/化学反应时间尺度的比值：当化学反应时间尺度（化学雅可比行列式的特征值）
过大时限制当地CFL数，以保持好条件矩阵。缺省值为0.9。 这些缺省值在大多数实例中都是可用的。
注意stiff-chemistry选项对非耦合求解器是没有的；它只能用于耦合求解器（隐式的或显式的）。 EDC模型求解步骤
如果你使用EDC模型，推荐使用双精度求解器（见1.5节），以避免刚性机理中固有的大指数前因子和活化能产生的截断误差。
由于EDC模型需要很大的计算开销，建议你采用以下步骤，以用非耦合求解器得到解： 1． 用涡耗散模型和简单的单步或两步放热机理计算一个初始解。
2． 用适当的物质使能EDC化学反应机理。如果你有一个CHEMKIN[112]格式的机理，参见13.1.9如何
将它导入。
3． 如果物质的数目和反应顺序改变，你将需要改变物质边界条件。
4． 通过关闭Species Model面板中的Volumetric Reaction选项暂时取消反应计算。
5．在Solution Controls面板中只使能物质方程的求解。
6．对物质混合场计算一个解。
7． 打开Species Model面板中的Volumetric Reaction选项，选定反应计算，并在
Turbulence-Chemistry Interaction下选择EDC模型。 8． 在Solution Controls面板中使能Energy方程的求解。
9． 对复合了物质和温度的场计算一个解。如果火焰吹熄，你可能还需要补缀一个高温区域。 10． 打开所有方程。 11． 计算最终解。
13.1.8 物质计算的后处理
FLUENT可以报告化学物质的质量分数、摩尔分数和摩尔浓度。你还可以显示层流和有效质量扩散系数。物质输运和反应模拟的后处理可以得到以下变量：
物质n的质量分数 物质n的摩尔分数 物质n的浓度
物质n的的层流扩散系数
WORD完美格式 下载可编辑
专业资料
物质n的有效扩散系数
物质n的的焓（仅对非耦合求解器） 物质n的源项（仅对非耦合求解器）
物质n的表面沉积速率（只对颗粒表面反应计算） 相对湿度
时间步长尺度（只对刚性化学反应求解器）
反应n的速率
反应n的Arrhenius速率 反应n的湍流速率
这些变量包含在后处理面板的变量选择下拉列表中的Species…，Temperature…，和Reactions…栏中。完整的流动变量、场函数和它们的定义见27章，25、26章解释了如何产生图形显示和数据报告。
平均物质浓度
使用Surface Integrals面板，可以得到模型中的入口、出口和选择面（即用Surface菜单项创建的表面）上的平均浓度。
在Field Variable 下拉列表中选择适当物质的Concentration of species-n。
13.1.9 从CHEMKIN导入一个化学反应机理
如果你有一个CHEMKIN格式的气相化学反应机理，你可以使用Chemkin Import面板（图13.1.7）将机理文件导入FLUENT。
WORD完美格式 下载可编辑
专业资料
在Chemkin Import面板中，在Chemkin Mechanism Flie下输入CHEMKIN文件的路径（如路径/file.che），并指定Thermodynamic Data File(THERMO.DB)的位置。有关数据库的更多信息见14.3.1节。
当你指定了两个文件的正确路径后，在Material Name下输入化学反应机理的名字，并点击Import按纽。FLUENT将以这个指定的名字创建一个种材料，它将包含物质和反应的CHEMKIN数据，并且将加到Material面板中的可得到Mixture Materials中去。
13.2 壁面表面反应和化学蒸汽沉积
对于气相反应，反应速率是在容积反应的基础上定义的，化学物质的形成和摧毁成为物质守恒方程中的一个源项。沉积的速率由化学反应动力和流体到表面的扩散速率控制。壁面表面反应因此在丰富相中创建了化学物质的源（和容器），并决定了表面物质的沉积速率。 有关表面反应和化学蒸气沉积的信息按以下小节给出：
13.2.1 概述和限制 13.2.2 理论 13.2.3 壁面表面反应的用户输入 13.2.4 壁面表面反应的求解过程 13.2.5 壁面表面反应的后处理
13.2.1 概述和限制
表面物质和壁面表面反应的概述
FLUENT把沉积在表面的化学物质与气体中的相同化学物质分开处理。类似地，涉及沉积的表面反应定义为单独的表面反应，因而其处理也与涉及相同化学物质的丰富相反应不同。举例说明，考虑从硅烷分解的如下硅沉积机理：
反应1（表面）：SiH4(g)Si(s)2H2(g)； 反应2（气体）：SiH4(g)SiH2(g)H2(g)； 反应3（气体）：SiH2(g)Si(g)H2(g)； 反应4（表面）：SiH2(g)Si(s)H2(g)；
WORD完美格式 下载可编辑
专业资料
反应5（表面）：Si(g)Si(s)；
在这一反应体系中，Si(s)和Si(g)作为两中不同的物质处理，你将需要分别定义。与气相物质相比，表面物质（如Si(s)）在流体相中不存在，不求解它们的质量守恒方程。FLUENT将预测每种表面物质的质量沉积速率。类似地，反应3和反应4也定义为两种不同的反应。你将定义每种反应为发生在流体相（反应3）或是发生在选择的与流体相邻的表面（反应4）。注意表面反应可能会被限制发生在某些壁面边界（而其它壁面边界保持不发生表面反应）。最后，表面反应是以单位面积来定义和计算的，而流体相反应基于单位容积。 限制
表面反应采用的连续方法在高Knudsen数（非常低压力下的流动）下不适用。
13.2.2 理论
壁面表面反应的Arrhenius反应速率
考虑如下的写为通用形式的第r个壁面表面反应：
其中，N-系统中总的化学物质数目
i,r-反应r中反应物I的化学计量系数 i,r-反应r中产物I的化学计量系数 Mi-物质I的记号
kf,r-反应r的正向速率常数
方程13.2-1中的和是针对系统中的所有化学物质，但只是那些具有非零化学计量系数的反应物和产物。因此，不涉及的物质将在方程中去除。
ˆ）由下式控制： 反应r中物质I的产生/摧毁摩尔速率（方程13.1-5中的Ri,r
其中，Nr-反应r中的化学物质数目
Cj,r-反应r中每个反应物和产物j的摩尔浓度（kgmol/m）
3
j,r-反应r中每个反应物和产物j的正向反应指数 反应r的正向反应常数kj,r按Arrhenius公式计算：
WORD完美格式 下载可编辑
专业资料
其中，Ar-指数前因子（恒定单位）
r -温度指数（无维）
Er -反应活化能（j/kmol）
R-通用气体常数（j/koml-K）
你（或数据库）需要提供i,r，i,r，j,r，r，Ar和Er的值。
壁面表面反应边界条件
对于壁面表面边界条件，反应表面的物质浓度计算基于进入（或离开）表面的每种物质的对流和扩散平衡，以及它在表面消耗（或产生）的速率。物质I的这种通量平衡可以记为：
其中，n为垂直于表面的单位矢量；
Jt为物质i的扩散通量；
Ri为由于表面反应的物质i产生速率；
dep为总的质量沉积速率； m Yi,wall为壁面上物质i的质量分数
用方程13.2-4可以得到壁面处物质i的质量分数和单位面积物质i的净产生速率的表达式。这些表达式
在FLUENT中用来计算反应表面处的气相物质浓度，采用点对点耦合刚性求解器。 在连续性中包括进入表面的质量输运
在以上所述的表面反应边界条件中，壁面法向速度的影响或输运到壁面的丰富相质量没有包括在物质输运的计算之中。离开表面的净质量通量的动量也忽略掉了，由于这一动量和邻近表面单元里的流动动量相比通常很小。但是，你可以通过激活Species Model面板中的Mass Deposition Source选项，在连续性方程中包括表面质量输运的影响。 能量方程中的壁面表面质量输运影响
能量方程中由于壁面表面反应产生的物质扩散影响包含在13.1.1节所述的正常物质扩散项中。
如果你使用非耦合求解器，可以通过关闭Species Model面板中的Diffusion Energy Source 选项忽略掉这一项。（对耦合求解器，通常是包括这一项的；你不能关闭）忽略物质扩散项意味着在涉及带有显著热容差异的物质混合，尤其是对带有Lewis数远大于1的物质的问题中，可能会对温度预测带来误差。当Le=1时物质扩散的影响应为零，你可能会发现由于物质和能量方程中数值积分的差异而产生微妙的影响。
模拟壁面表面反应产生的放热
缺省时，FLUENT忽略壁面表面反应产生的放热。但你可以通过激活Species Model面板中的Heat of Surface Reactions选项，并设定Materials面板中的生成焓来选择包含表面反应的热量。
13.2.3 壁面反应的用户输入
WORD完美格式 下载可编辑
专业资料
设置一个涉及壁面反应问题的基本步骤与13.1.2节中所述的设置一个只有气相反应的问题大致相同，但有一些额外的设置：
1． 在Species Model面板中：
（a） 使能Species Transport，选择Reactions下的Volumetric和Wall Surface，并指定
Mixture Material。这一过程的细节见13.1.3节及13.1.2节关于混合物材料概念的解释。
（b） （可选）如果希望模拟壁面反应的放热，打开Heat of Surface Reactions选项。
（c） （可选）如果希望在连续性方程中包括表面质量输运的影响，打开Mass Deposition
Source选项。
（d） （可选）如果使用非耦合求解器，并且不希望在能量方程中包括物质扩散的影响，关闭
Diffusion Energy Source选项。细节见13.2.2节。
（e） （可选，但对CVD推荐）如果希望模拟完整的多组分扩散或热扩散，打开Full
Multicomponent Diffusion 或Thermal Diffusion选项。细节见7.7.2节。
2．检查和/或定义混合物属性（见13.1.4节）：
混合物属性包括：
混合物中的物质 反应
其他物理属性（如粘度，比热）
！你将在Fluid Materials列表中找到所有物质（包括表面物质）。对于象Si一类的物质，你可以在流体材料类型的材料列表中找到Si(g)和Si(s)。如果你在模拟13.2.1节开始处实例中的硅沉积反应，你将需要在混合物中包括两中Si元素（气体和固体）。
！注意如果你的模型中包括稀释混合物中的物质，Selectecd Species列表中的最终气相物质必须是载体气体。（因为FLUENT不会求解最后物质的输运方程。）还需要注意的是任何物质的重排、增减都必须小心处理，如13.1.4节所述。
3．检查和/或设定混合物中单独物质的属性（见13.1.4节）。注意如果你模拟表面反应的放热，你必须检查（或定义）每种物质的生成焓。 4．设置物质边界条件
除了13.1.5节中所述的边界条件之外，你还需要指定表面反应对每个壁面是否有效，并考虑热边界条件的选择。为使能一个壁面上的表面反应影响，在Wall面板的Species区域打开Surface Reactions选项。壁面的边界条件输入细节见6.13.1节。
！当在一个给定壁面使能表面反应后，这一壁面上对混合物材料定义的所有表面反应都被激活。
13.2.4 壁面反应的求解过程
正如所有的CFD模拟一样，如果模拟从一个简单的问题描述开始，在求解向前推进时增加复杂性，可能会使你的表面反应模拟工作更加成功。对于表面反应，你可以遵从13.1.7节中对流体相反应给出的指南。 另外，如果你模拟表面反应的放热，而且遇到了收敛性方面的麻烦，你可以尝试暂时关闭Species Model
WORD完美格式 下载可编辑
专业资料
面板中的Surface reactions和Mass Deposition Source选项4。
13.2.5 表面反应的后处理
对壁面反应，除了13.1.8节中列出的变量之外，你还可以显示/报告沉积在一个表面上的固体物质的沉积速率。在变量选择下拉列表中的Species…栏中选择物质n的Surface Deposition Rate 。
13.3 颗粒表面反应
如19.3.6节所述的，可以定义多颗粒表面反应以模拟离散相颗粒的表面燃烧。颗粒表面反应的信息按以下小节给出：
13.3.1 理论
13.3.2 颗粒表面反应的用户输入
13.3.3 使用离散相颗粒燃烧的多表面反应
13.3.1 理论 一般描述
Smith提出了计算碳粒燃烧速率的关系，并进行了详细讨论[218]。颗粒反应速率R可以表示为：
其中，D0-bulk扩散系数；
； Cg大量物质中的平均反应气体物质浓度（kg/m3）
Cs颗粒表面的平均反应气体物质浓度（kg/m3） Rc化学反应速率系数；
N-显式反应级数（无维）；
在方程13.3-1中，颗粒表面处的浓度Cs是未知的，因此需要消掉，表达式改写为如下形式：
这一方程需要通过一个迭代过程求解，除去N1或N0的特例。当N1时，方程13.3-2可以写
为：
在N0情况下，如果在颗粒表面具有有限的反应物浓度，固体损耗速度等于化学反应的速度。如果在表面没有反应物，固体损耗速度根据扩散控制速率突然变化。在这种情况下，处于稳定性的原因，FLUENT采用化学反应速率。
FLUENT模型形式
WORD完美格式 下载可编辑
专业资料
图13.3.1表示了一个正在气相中进行放热反应的颗粒。Tp和T是方程19.3-78中的温度。
基于以上分析，FLUENT用以下方程表示在气相物质n中颗粒表面物质i的第r个反应的速率。在这中情况下，反应r的化学计量表达式描述为
反应的速率为：
其中，Rj,r-颗粒表面物质的消耗速率（kg/s）； ； Ap-颗粒表面积（m2）
Yj-颗粒状态的表面物质j的质量分数；
； r -效率因子（无维）
； Rj,r-单位面积的颗粒表面物质反应速率（kg/m2-s）； pn-气相物质的分压力（Pa）
D0,r-反应r的扩散速率系数
； Rkin,r-反应r的动力学速率（单位变化）
Nr-反应r的显式级数
WORD完美格式 下载可编辑
专业资料
效率因子r与表面积有关，可以用于多反应例子中的每一个反应。D0,r由下式给出：
反应r的动力学速率定义为：
反应级数Nr1的颗粒表面物质消耗速率由下式给出：
反应级数Nr0时，
多气相反应物的化学表达式扩展
当有超过一种气相反应物参与到反应中时，反应的表达式需要扩展：
为表示有nmax种气相物质出现的反应中某种颗粒表面物质j的第r个反应的速率，需要对每种固体颗粒反应定义扩散限制的物质，即在载体和颗粒表面之间浓度梯度最大的物质。对其它物质，表面和载体浓度认为相等。扩散限制物质的浓度在图13.3.1中一Cd,b和Cd,s表示，其它物质的浓度记为Ck。对多气相反应物的化学表达式，方程13.3-4和13.3-8中的载体分压力pn就是反应r中扩散限制物质的大量分压力
pr,d。
反应r的动力学速率定义为：
其中，pn-气体物质n的大量分压力
Nr,n-物质n的反应级数。
当选定这一模型后，常数C1,r（方程13.3-6）和 效率因子r（方程13.3-4）在Reactions面板中输入，如13.3.2节所述。
13.3.2 颗粒表面反应的用户输入
颗粒表面反应的输入过程只需要在13.1.2节-13.1.6节介绍的容积反应过程的基础上增加一些输入即
WORD完美格式 下载可编辑
专业资料
可。这些附加的输入如下：
在Species Model面板，打开Reactions下的Particle Surface选项。
当你指定颗粒表面反应中涉及的物质时，确保确定表面物质，如13.1.4节所述。 ！ 你将在Fluid Materials列表中找到所有的物质（包括表面物质）。例如，如果你模拟煤气化，你可以在流体材料类型的材料列表中找到固态碳C(s)。
对每一种颗粒表面反应，在Reactions面板中选择Reaction Type为Particle Surface，并
指定以下参数（除了13.1.4节所述的之外）：
扩散限制物质：当有超过一种的气相反应物参与颗粒表面反应时，扩散限制物质是在载体和颗粒表面反应之间浓度梯度最大的物质。这一概念的描述见13.3.1节。在大多数情况下，只有一种气相反应物，扩散限制物质没有必要定义。
扩散速率常数：方程13.3-6中的C1,r 效率因子：方程13.3-4中的r
13.3.3 对离散相颗粒燃烧使用多表面反应模型
当你使用多表面反应模型时，设定一个涉及离散相问题的过程与19.6节概括的过程稍有不同。修正的过程如下：
1． 使能所有相关的离散相模拟选项，如19.7节所述。 2． 指定初始条件，如19.9节所述。 3． 定义边界条件，如19.10节所述。 4． 定义材料属性，如19.11节所述。
！在进入下一步之前，先需要在Materials面板中的Combustion Model下拉列表中选择multi-surface-reactions。
5． 如果你已经定义了超过一种颗粒表面物质，如碳（Cs）和硫（Ss），你需要返回到Set Injection Properties面板（或Set Multiple Injection Properties面板），指定燃烧颗粒中每种颗粒表面物质的质量分数。点击Multiple Reactions表格，并输入Species Mass Fractions。这些质量分数指的是燃烧颗粒中的可燃组分，和必须为1。如果混合物材料中只有一种表面物质，这种物质的质量分数将被设定为1，你不需要在Multiple Surface Reactions下指定任何东西。 6． 求解参数并求解问题，如19.12节所述。 7． 检查结果，如19.13节所述。
！ 当使用多表面反应模型时，不能进行非稳态粒子追踪。
13.4 无反应的物质输运
除了上述的容积和表面反应之外，你还可以用FLUENT求解无反应的物质混合问题。FLUENT要求解的物质输运方程在13.1.1节中已经介绍，设定无反应物质输运问题的过程与13.1.2-13.1.6节中所述的过程相同，某些地方有所简化。 基本步骤如下所列：
1．在Species Model面板中选定Species Transport，并选择适当的混合物材料。
WORD完美格式 下载可编辑
专业资料
有关混合物材料的概念见13.1.2节，使用Species Model面板的更多细节见13.1.3节。 2．（可选）如果你希望模拟完全多组分扩散或热扩散，打开Full Multicomponent或Thermal Diffusion选项。细节见7.7.2节。
3．检查和/或定义混合物的属性及其构成物质。
混合物属性包括： 混合物中的物质
其他物理属性（如粘度、比热） 细节见13.1.4节。
4．设定边界条件，如13.1.5节所述。
无反应的物质输运一般不需要特殊的求解过程。计算完成后，你可以显示或报告以下量： 物质n的质量分数 物质n的摩尔分数 物质n的浓度
物质n的的层流扩散系数 物质n的有效扩散系数
物质n的的焓（仅对非耦合求解器） 相对湿度
这些变量包含在后处理面板中出现的变量选择下拉列表中的Species…栏。流动变量、场函数及其定义的完整列表见27章。25章和26章解释如何产生图形显示和数据报告。
WORD完美格式 下载可编辑
因篇幅问题不能全部显示，请点此查看更多更全内容