化学动能

化学动能涉及化学组分在分子级别的反应机制和反应率。

对反应组分输运进行建模时,为传输方程提供源项需要计算化学机制中所有反应的反应率(化学动能)。

任何基元化学反应均可用以下常规方程表示。
1. EQUATION_DISPLAY
υAA+υBB+...kbkfυC′′C+υD′′D+...
(3353)
其中,带上标 ' 的化学计量系数 υ 表示反应物,带 '' 的系数表示产物。净化学计量系数 υij 给出反应 j 产生或消耗的组分 i 的摩尔总数:
2. EQUATION_DISPLAY
υij=υiυi
(3354)
反应可以正向进行(即,反应物变为产物)或逆向进行(即,产物形成反应物)。 kf 为正反应率常数, kb 为逆反应率常数。 k 为同时考虑 kf kb 的总速率常数。
3. EQUATION_DISPLAY
k=kfkb
(3355)

反应率

基元反应的反应率计算如下。

以反应机制中的一个反应为例:
4. EQUATION_DISPLAY
2CO+O2k2CO2
(3356)

有两种反应物组分和一种产物组分。要使反应继续下去,反应物分子必须与等于或大于激活能 Ea 的能量碰撞。总体反应率 rk 还取决于速率常数 k 和两种反应物的浓度(是其化学计量系数的次幂)。

Eqn. (3356) 的速率方程写作:
5. EQUATION_DISPLAY
rk=kf[CO]2[O2]kb[CO2]2
(3357)
对于常规反应,Eqn. (3353),这表示为:
6. EQUATION_DISPLAY
rk=kf[A]υA[B]υB...kb[C]υC[D]υD...
(3358)

方括号中组分的化学式表示该组分在反应混合物中的浓度。

通过将反应率乘以净化学计量系数 Eqn. (3354),计算反应的每个组分的变化率。例如,Eqn. (3356) 中 CO 的反应率为:
7. EQUATION_DISPLAY
d[CO]dt=2(kf[CO]2[O2]+kb[CO2]2)
(3359)
对于常规一步反应,Eqn. (3353),这表示为:
8. EQUATION_DISPLAY
d[A]dt=υA[kf[A]υA[B]υB...kb[C]υC[D]υD...]
(3360)
其中
υA=υAυA
(3361)

组分 i 的生成(或消耗)率 ri 取决于生成和消耗该组分的反应的反应率以及参与这些反应的所有组分的浓度。因此,总净速率为组分 i 参与的所有反应的 Eqn. (3360) 之和。在以下方程中, j (= 1、2… Nj )表示反应, i 表示组分。

组分 i 的净反应率,表示为 ω˙i ,是所有反应( rj ,位于 Eqn. (3358))之和乘以化学计量系数。

9. EQUATION_DISPLAY
ω˙i=d[A]dt=j=1Njυijrj
(3362)
10. EQUATION_DISPLAY
rj=kf,jΠi[A]υj,ikb,jΠi[A]υj,i
(3363)

其中,正反应率系数 kf 和逆反应率系数 kb Eqn. (3365) 确定。

用户定义反应率的计算方法如下:

11. EQUATION_DISPLAY
rj=kjΠi[A]υj,i
(3364)
其中, kj 为反应 j 的用户反应系数。

速率常数

速率常数 k 通常使用阿雷尼乌斯方程确定,其中速率与温度非线性相关:
12. EQUATION_DISPLAY
k=ATβexp(EaRuT)
(3365)
其中, Ru 为通用气体常数,以下值为化学机制提供的经验因子:
  • 温度指数 β
  • 活化能 Ea
  • 频率因子或指前因子 A

化学放热率

化学放热率定义为:

13. EQUATION_DISPLAY
h˙=ω˙khko
(3366)

其中, ω˙k kth 组分的净反应率, hko kth 组分的形成焓。使用组分输运模型时,Eqn. (3366) 可作为名为 Chemical Heat Release Rate(化学放热率) 的场函数计算和提供。

对于小火焰生成流形 (FGM) 以及部分预混化学平衡 (CE) 和稳态层流小火焰 (SLF) 模型,化学放热率计算如下:

14. EQUATION_DISPLAY
h˙=c˙hR
(3367)

其中, c˙ 为进程变量源项, hR 为反应热。 hR 的计算需要假设燃料在一步反应中与氧气发生碳氢化合物反应:

15. EQUATION_DISPLAY
CnHm+(n+m4)O2nCO2+m2H2O
(3368)

反应热计算如下:

16. EQUATION_DISPLAY
hR=YkFhkonhCO2om2hH2Oo
(3369)

其中, YkF 为燃料流中的 kth 组分的质量分数。

对于非预混 SLF,化学放热率计算如下:

17. EQUATION_DISPLAY
h˙=χhR
(3370)

其中, χ 为标量耗散率。

注意,小火焰模型的化学放热率近似 Eqn. (3366) 中的精确放热率,场函数名为 Chemical Heat Release Rate Indicator(化学放热率指示器)

第三体反应

第三体效率
一些反应需要“第三体”才能进行,如下所示:
18. EQUATION_DISPLAY
H+O2+MHO2+M
(3371)
其中, M 表示第三体。第三体充当一种催化剂。通常,混合物中的任何组分均可作为第三体。第三体影响反应进行的速率。
需要第三体的反应的反应率的表达式如下:
19. EQUATION_DISPLAY
rj=(Σi=1N(αij)[Xi])(kf,jΠi=1N[Xi]υijkb,jΠi=1N[Xi]υij)
(3372)
对于组分 (i=1,...,N)
  • αij 表示第三体的碰撞效率
  • kf,j 为反应的正向速率常数 j
  • kb,j 为反应的反应反向速率常数 j
αij=1 时,混合物中的所有组分都将作为第三体同等发挥作用,混合物的总浓度为:
20. EQUATION_DISPLAY
[M]=Σi=1N[Xi]
(3373)
但是,当 αij1 时,反应中的有效第三体的浓度为:
21. EQUATION_DISPLAY
[M]=Σi=1N(αij)[Xi]
(3374)
其中, αij 是为反应机制中的第三体组分指定的值。
压力相关化学反应率
压力相关化学反应(脱落反应)在低压下为一阶比率,随着第三体组分的浓度增加而变为零阶。在此类反应中,速率受压力和温度影响,Simcenter STAR-CCM+ 使用反应机制提供的指前因子 A 计算速率常数:
22. EQUATION_DISPLAY
k=k(Pr1+Pr)F
(3375)
其中,折算压力 Pr 为:
23. EQUATION_DISPLAY
Pr=k0[M]k
(3376)
其中, [M] 为具有第三体效率的组分 M 的浓度。
低压限值 k0 为:
24. EQUATION_DISPLAY
k0=A0Tβ0exp(E0RuT)
(3377)
高压限值 k 为:
25. EQUATION_DISPLAY
k=ATβexp(ERuT)
(3378)
混合(脱落)函数 F 的计算方法有所不同,具体取决于指定压力相关参数所使用的方法。
  • 使用 Lindemann 方法时, F 计算为合一值, F=1
  • 使用 Troe 方法时, F 计算如下:
    26. EQUATION_DISPLAY
    logF=[1+[logPr+cnd(logPr+c)]2]1logFcent
    (3379)
    其中:
    27. EQUATION_DISPLAY
    Fcent=(1α)exp(TT***)+αexp(TT*)+exp(T**T)
    (3380)
    • c=0.40.67logFcent
    • n=0.751.27logFcent
    • d=0.14
    • α T*** T* T** 分别为在已导入的化学机制中指定的常数
  • 使用 (SRI) 方法时, F 计算如下:
    28. EQUATION_DISPLAY
    F=d[aexp(bT)+exp(Tc)]XTe
    (3381)
    其中,a、b、c、d 和 e 是在导入的化学机制中指定的常数。如果机制中未指定 d 和 e,则默认情况下会将它们分别设为 1 和 0。
    29. EQUATION_DISPLAY
    X=11+log2Pr
    (3382)
对于 Troe 和 SRI 方法,高压和低压限值的 F 都接近合一值。

化学活性双分子反应率

在化学活性双分子反应中,由于反应中间体的碰撞稳定,因此反应率随压力的增加而降低(脱落)。正反应率常数与低压速率常数成比例:
30. EQUATION_DISPLAY
k=k0(11+Pr)F
(3383)
混合函数 F 也使用 Troe、Lindemann 或 SRI 方法定义。

基于对数插值的压力依赖性 (PLOG)

使用 PLOG 方法描述反应率的压力依赖性时,可使用已针对不同压力进行修正的速率参数( β Ea A )通过阿雷尼乌斯表达式 Eqn. (3365) 计算反应率。

为相关压力范围内的一组(至少两个)不同离散压力给定了速率参数。

  • 当模拟中的压力 P 介于给定压力之间(例如,介于 Pi Pi+1 之间)时,可通过将 lnk 线性插值为 lnP 的函数来计算速率参数:
    31. EQUATION_DISPLAY
    lnk=lnki+(lnki+1lnki)lnPlnPilnPi+1lnPi
    (3384)
  • 当模拟中的压力高于或低于任何给定压力时,分别使用最高或最低压力的速率参数。

化学平衡

当正反应率和逆反应率相等时,反应被视为已达到平衡且未观察到净反应。只要压力(对于气相反应)、温度和组分浓度保持恒定,反应物和产物的比例即为常数。

可以使用平衡常数 K 确定在平衡状态下每种组分的数量。

对于常规反应,Eqn. (3353),平衡常数定义为:
32. EQUATION_DISPLAY
K=(aC)υC(aD)υD...(aA)υA(aB)υB...
(3385)
其中, a 为相应组分的活性。气相组分(例如,组分 A )的活性定义为:
33. EQUATION_DISPLAY
aA=[A]c0
(3386)

其中, c0 为参考浓度。

在化学平衡状态下:
34. EQUATION_DISPLAY
kf[A]υA[B]υB...=+kb[C]υC[D]υD...
(3387)
也可以表示为:
35. EQUATION_DISPLAY
[C]υC[D]υD...[A]υA[B]υB...=kfkb
(3388)
对于在平衡状态下已知组分浓度的反应,相对浓度 Kc 的平衡常数可由下式确定:
36. EQUATION_DISPLAY
Kc=[C]υC[D]υD[A]υA[B]υB
(3389)
这意味着:
37. EQUATION_DISPLAY
Kc=kfkb
(3390)
相对浓度的平衡常数也可以使用亥姆霍兹自由能 A 根据热力学数据确定:
38. EQUATION_DISPLAY
Kc=exp(ΔRA¯0RuT)
(3391)
或者,对于相对压力的平衡常数,使用吉布斯自由能 G 确定:
Kp=exp(ΔRG¯0RuT)
(3392)

表面化学

在表面化学反应期间,可能发生三个过程:
  • 气相分子与表面分子碰撞,然后被吸附。表面点位(表面的活性分子)成为块体的一部分,先前的气相分子在块体顶部形成新的表面。
  • 表面上的分子之间或表面上的一个分子和气体中的一个分子之间的表面反应会形成新的表面、气相或块体分子。
  • 分子在反应后进入气相,或者只需克服键能,即可解吸表面分子。

例如,在下图中,气体分子被吸附在先前吸附的表面组分旁边的表面上。气体分子分解为其构成原子。吸附的某些原子发生反应,形成离开表面的小气体分子。气体分子中的一个原子仍然作为表面组分被吸附,而先前吸附的表面组分则成为块体组分。



上图中的反应可以表示为:

SiH4(G)+Si(S)Si(S)+Si(B)+2H2(G)

如果先前覆盖块体组分的表面分子被解吸或参与反应,则块体组分不会被视为活性,但最终可以位于表面点位上。

使用表面气体相互作用模型必须定义三种类型的组分:
  • 气相中的组分, G
  • 气体表面交界面处表面上的组分, S 。表面反应可描述为使用开放点位形式或原子点位形式。
  • 固体层中的组分,块体组分, B

反应 j 的通用表面反应形式如下:

39. EQUATION_DISPLAY
i=1Ngg′i,jGi+i=1Nss′i,jSi+i=1Nbb′i,jBi=i=1Nggi,jGi+i=1Nssi,jSi+i=1Nbbi,jBi
(3393)

其中:

  • Gi Si Bi 为气体、点位和块体组分。
  • Ng Ns Nb 为气体、点位和块体组分的数量。
  • g′i,j s′i,j b′i,j 为气体、点位和块体反应物的化学计量系数。
  • gi,j si,j bi,j 为气体、点位和块体产物的化学计量系数。

j 个反应的速率为:

40. EQUATION_DISPLAY
rj=kf,jΠi=1Ng[Gi]wallpg(i,j) Πi=1Ns[Si]wallps(i,j)
(3394)

其中, pg ps 为反应物的比率指数, [*]wall 表示壁面上的摩尔浓度。反应 j 的正向速率常数 kf,j 由阿雷尼乌斯表达式给出。对于最简单的反应类型,该常数由 Eqn. (3365) 给出(其中, k kf,j )。

对于第 n 种点位类型,表面组分 [Si] 的浓度(以摩尔/平方米为单位)由点位密度 Γn 、第 i 个组分 ZSi,n 覆盖的点位分数和组分覆盖 σSi,n 的点位数量确定。
41. EQUATION_DISPLAY
[Si]=ΣnZSi,nΓnσSi,n
(3395)

总和将在所有不同的点位类型中运行。

假设反应率与块体组分的浓度无关。

每个组分 i 每单位表面积的净摩尔变化率(结果或消耗)由以下公式给出:

42. EQUATION_DISPLAY
rˆi,gas=j=1Nreac(gi,j-g′i,j)rj
(3396)
43. EQUATION_DISPLAY
rˆi,site=j=1Nreac(si,j-s′i,j)rj
(3397)
44. EQUATION_DISPLAY
rˆi,bulk=j=1Nreac(bi,j-b′i,j)rj
(3398)

对于某些表面反应,反应率受所涉及组分的表面覆盖特征影响。可以修正阿雷尼乌斯速率表达式 Eqn. (3365) 来创建校正速率常数,如下所示:

  • 点位覆盖校正率表达式

    某些表面反应由组分 i 的表面点位覆盖率 zi 体现。

    45. EQUATION_DISPLAY
    kf,j=AjTβjexp(EaRuT)Πi=1Ns10ηi,jziμi,jexp(εi,jZiRuT)
    (3399)

    其中, ηi,j μi,j εi,j 是可为反应表面上的组分 i 和反应 j 指定的三个表面覆盖参数。在导入的反应机制文件中,表面覆盖修正反应由关键词 COV 表示。请参见反应机制格式

  • 粘附系数校正率表达式
    可以使用粘附系数或粘附概率(而非反应率系数)表示正向表面反应率。粘附系数表示组分粘附到表面的概率。反应可以包括一个气相组分和任意数量的表面组分,且反应率与所有涉及的表面组分的覆盖率成比例。粘附系数 sc 计算如下:
    46. EQUATION_DISPLAY
    sc=min1,AjTβjeEa,jRuT
    (3400)
    速率表达式修改如下:
    47. EQUATION_DISPLAY
    kf,j=scΠi=1Nsσiυi,j(Γtot)mRuT2πMwi
    (3401)

    Γtot 是总表面站点浓度、 m 是所有表面反应物化学计量系数之和。 Mwi 为气相组分 的分子量,而 为每个表面组分占据的点位数 - 为该组分的反应级数。 i σi υi,j

  • Motz-Wise 校正率表达式
    使用 Motz-Wise 校正因子的前向反应率表达式为:
    48. EQUATION_DISPLAY
    kf,j=(sc1sc/2)Πi=1Nsσiυi,j(Γtot)mRuT2πMwi
    (3402)
    其中, sc Eqn. (3400) 计算。
  • 博姆校正率表达式
    模拟包含正离子的反应时,可以修改速率常数来包含 Bohm 速度校正 RRuuTTeeMMwwii ,这考虑到了反应的发生概率:
    49. EQUATION_DISPLAY
    kf,j,Bohm=AjTβjeEa,j/RTΠi=1Nsσiυi,j(Γtot)mRuTeMwi
    (3403)

    其中 为电子温度,而 为正离子的分子量。 Te Mwi

  • Langmuir-Hinshelwood

    Langmuir-Hinshelwood 为全局表面反应,将组分吸附、表面反应和解吸合并到单个步骤中。Langmuir-Hinshelwood 全局速率表达式适用于假设吸附和解吸处于平衡状态,并且表面上吸附组分之间的反应决定速率的情况。假设表面上的吸附点位相互独立(单点位吸附),则这些点位是等效的,并且表面覆盖仅会减少可供吸附的点位数,但不会改变吸附/解吸的能量特性。

    通过吸附/解吸平衡,包含由各种组分阻碍的表面点位的效应。
    过程速率由以下公式给出:
    54. EQUATION_DISPLAY
    q=kΠi[Xi]li(1+ΣiKi[Xi]ni)m
    (3408)
    其中 为速率常数,由 Eqn. (3365) 求出, 为气相组分 的浓度,而 为吸附/解吸步骤的平衡常数,由以下公式给出: k [Xi] i Ki
    55. EQUATION_DISPLAY
    Ki=AiTβiexp(Hi/RT)
    (3409)
    上标指数 m ,当 m=2 时对应于 Langmuir-Hinshelwood 反应,当 m=1 时对应于 Eley-Rideal 反应。

颗粒化学

对于颗粒反应,速率常数 k 使用阿雷尼乌斯方程 Eqn. (3365) 确定。由于颗粒化学反应率考虑随动能产生的扩散,因此理论指南的颗粒反应一节介绍了这些计算。