如何选择正确的电流分布接口?

2014年 2月 10日

在设计电化学电池时,我们需要考虑电解质和电极中的三类电流分布:一次分布二次分布三次分布。不久之前,我们介绍了电流分布的基本理论;本文则以线电极为例,详细解释不同的电流分布类型,帮助你在 COMSOL Multiphysics 中选择合适的电流分布接口,顺利执行电化学电池仿真。

三种电流分布接口

上一篇博客文章一样,我们可以借助线电极模型的案例比较三种电流分布接口。下方是线电极的几何结构示意图:

以电化学电池的几何模型为例介绍电流分布理论
使用 COMSOL Multiphysics 求解的线电极模型。电解质可以在导线和平坦表面之间的开放空间内流动。

本文提及的三种电流分布情形采用了相同的几何模型:在两个平坦的电极表面之间放置一个线电极结构。该电化学单元电池可视为一个体积更大的线网格电极中的一个单元,后者是许多大规模工业生产过程中常见的电化学电池装置。

回顾基本方程

以下是我们在上一篇博客中详细说明的基本方程式:

Nernst-Planck 方程:

(1)

\textbf{N}_i = -D_i\nabla c_i-z_i u_{m,i} F c_i\nabla \phi_l+c_i\textbf{u}

包含 Nernst-Planck 方程的电流密度表达式:

(2)

\textbf{i}_l = -F \left(\nabla \sum_i z_iD_i c_i\right)-F^2\nabla \phi_l \sum_i z^2_i u_{m,i} c_i+\textbf{u}F\sum_i z_ic_i

一般电解质电流守恒:

\nabla\cdot\mathbf{i}_l=Q_l

一次电流分布

一次电流分布仅适用于解释溶液电阻造成的损耗,而忽略了电极动力学和浓度依赖性效应造成的损耗。它假定电解液中的电荷转移遵守欧姆定律。在这里,我们做出了两个假设:第一,电解质是电中性的,它抵消了对流项在表达式(2)中对电流密度的贡献;第二,电解质的成分变化微不足道(即分布均匀),这抵消了扩散项在表达式(2)中对电流密度的贡献,允许我们将离子强度看作一个常数。由此一来,表达式(2)的剩余项即为欧姆定律用于描述电解质电流密度。

我们认为,电极-电解质界面上的电解反应发生得非常快,以至于可以忽略电极动力学的影响,因此,电极-电解质界面上的电势差基本不会偏离平衡值。换句话说,没有活化过电位,而且通过电解会产生任意的电流密度。由此可见,一次电流分布仅取决于正极和负极的几何构造。

COMSOL Multiphysics 的一次电流分布 接口定义了两个相关变量:一个是电解质电势(\phi_l\),另一个是电极电势(\phi_s\)。基于上述的一次电流分布假设,可得以下方程:

电极:\textbf{i}_s = -\sigma_s\nabla\phi_s\\nabla\cdot\textbf{i}_s = Q_s

电解质:\textbf{i}_l = -\sigma_l\nabla \phi_l\\nabla\cdot\textbf{i}_l = Q_l

电极与电解质界面:\phi_s-\phi_l = E_{\mathrm{eq},m}

其中,\sigma_l 表示电解质的电导率,根据上述假设它是一个常数。下标 s 代表电极,l 代表电解质。E_{\mathrm{eq},m} 表示反应 m 的平衡电位。

在下图中,我们展示了线电极示例的一次电流密度分布。如您所见,电流密度分布在直面负极板的导线的拐角处,在导线结构的中心区域则趋近于零,原因是从几何结构上该区域与负极隔开了。

一次电流分布
一次电流分布,Ecell = 1.65 V。正极上的电流密度分布(无量纲)。

何时使用一次电流分布 接口?

如果电池的电解质浓度(相对于电流密度)较高,或者经剧烈搅拌混合后,电解质浓度可以假设为均匀。并且,电化学反应速度非常快,以至于反应电阻的大小与电阻损耗(溶液电阻)相比可忽略不计。如此,便可以使用一次电流分布来模拟电池。

上述条件的应用之一是在施加外加电流阴极保护系统的阳极上,与此同时,在一次电流密度分布接口中,我们可以将相当于氧还原反应中受质量传递限制的电流,以恒定电流的条件设为阴极的边界条件(参照案例)。这也是对反应速度相对较快的电化学过程的一个有效近似,例如氯碱工艺流程中氯离子的快速氧化。

由于一次电流分布 接口易于求解,且不涉及非线性动力学表达式,因此通常适合于计算基准近似结果,为处理较复杂的模型做好准备。

二次电流分布

二次电流分布用于解释电极动力学和溶液电阻的影响。它关于电解质的组成和特性的假设与一次电流分布是一致的,也就是说电解质电流同样遵守欧姆定律。一次和二次电流分布的区别在于对电解质和电极界面上的电化学反应的描述。

二次电流分布引入了电极动力学的影响;由于电解反应速率有限导致额外的阻抗,电位差可能偏离其平衡值。实际电位差和平衡电位差的差值即活化过电位(\eta)。因此,你会发现一次电流分布二次电流分布 接口的域方程完全相同,但是由于过电位的存在,电极-电解质界面的电位方程式存在差别:

电极:\textbf{i}_s = -\sigma_s\nabla\phi_s\\nabla\cdot\textbf{i}_s = Q_s

电解质:\textbf{i}_l = -\sigma_l\nabla \phi_l\\nabla\cdot\textbf{i}_l = Q_l

电极与电解质界面:\eta_m = \phi_s-\phi_l-E_{eq,m}

二次电流分布 接口中,由于电化学反应而产生的电流密度被描述为过电位的函数。物理场接口可以使用电流密度和过电位之间的任何关系,例如 COMSOL 的常用内置选项 Butler-Volmer 方程式(3)和 Tafel 方程式。

(3)

i_{loc,m} = i_{0,m}\left(e^\frac{\alpha_{a,m} F \eta_m}{RT}-e^\frac{-\alpha_{c,m} F \eta_m }{RT}\right)

上述 Butler-Volmer 方程针对反应 m 规定:i_{loc,m} 代表局部电荷转移电流密度,i_{0,m} 代表交换电流密度,\alpha_{a,m} 代表阳极电荷传递系数,\alpha_{c,m} 代表阴极电荷传递系数。R 是通用气体常数。这个方程式描述的情形是:单电子的电荷转移是净电荷转移反应的速率控制步。我们可以以描述均质化学反应的 Arrhenius 方程进行类比,从而推导出方程式,前提是假设带电物质的自由能会受到电位的影响。因此,活化能将随电极-电解质界面的电位差而变化。

根据下方公式,将全部电极的反应电流之和用作电极和电解质域之间边界上的电流密度条件:

-\textbf{i}_s \cdot \textbf{n} = \textbf{i}_l \cdot \textbf{n} = i_\mathrm{DL} + \sum_m{i_{loc,m}}

附加电容电流 i_\mathrm{DL} 产生于双电层的充电和放电过程。

一般情况下,通过活化过电位来解释电极动力学的影响,往往会使电流分布更加均匀。下图中的线电极示例说明了这一点。

与一次电流分布相比,二次电流分布更加平滑,最小值与最大值之间的差异更小。引入了活化过电位的影响后,局部高电流密度将在电极表面引发高的局部活化过电位,自然导致电流流向其他的路径。换句话说,你可以理解为电化学反应以有限的速率进行。在某些区域,反应动力学受到了限制,因此表面的电流密度分布会比不受限制的快速反应的情形更加平滑。

二次电流分布
二次电流分布,Ecell = 1.65 V。正极的电流密度分布(无量纲)。

何时使用二次电流分布 接口?

二次电流分布 是电化学工业应用建模的“主力”接口。你可以将二次电流分布用于电池建模,前提是电解质充分混合或浓度相对较高,因而可以忽略浓度过电位,若电极动力学引起的损耗与电阻损耗相比不可忽略,则不可行。众所周知,在工业应用中充分混合高浓度的电解质不是一个难题。你还可以将使用二次电流分布 接口作为电化学电池仿真的第一步,先计算出活化损耗,之后再引入浓度依赖反应动力学。

三次电流分布

三次电流分布用于解释电解质组成和离子强度的变化对电化学过程的影响,以及溶液电阻和电极动力学的影响。为此,它完整求解了 Nernst-Planck 方程(1),借此描述每一种化学物质通过扩散、迁移和对流的质量传输。此外,物质的浓度需要满足电中性假设。电化学反应的动力学表达式解释了活化过电位和浓度过电位,这意味着在电极-电解质界面,反应物的耗尽会导致电解反应的速率受到传输的限制。这也暗示了,电解液中的全部离子和全部活性物质都必须包含在模型中。

与一次和二次电流分布不同的是,三次电流分布中的电解质电流密度不再假定遵循欧姆定律。根据方程式(2),电中性的施加仍然意味着对流通量对电流密度没有贡献,但现在不能忽视电解质浓度变化带来的影响。因此,方程式(2)中的扩散项可能是非零项。

在电极-电解质的界面上,电荷转移反应的电流密度不仅被表示为过电位的函数,也是界面处电活性物质浓度的函数。相对于取决于单电子电荷转移步骤的反应速率,反应动力学方程采用了计算电荷传递电流密度 i_{loc,m} 的 Butler-Volmer 表达式(与方程式(3)对比),在三次电流分布情况下表达式还包含了浓度依赖性。

COMSOL 软件中的三次电流分布 接口求解了电解质电位(\phi_l\)、电极电位(\phi_s\)和一组物质浓度 c_i。根据上述假设,可得到以下方程:

电极:\textbf{i}_s = -\sigma_s\nabla\phi_s\\nabla\cdot\textbf{i}_s = Q_s

电解质:\textbf{i}_l = F\sum_{i=1}^n z_i (-D_i\nabla c_i-z_i u_{m,i} F c_i\nabla \phi_l)\nabla\cdot\textbf{i}_l = Q_l

电解质电中性:\sum_i z_ic_i = 0

电极与电解质界面:\eta_m = \phi_s-\phi_l-E_{eq,m}

典型的电流密度表达式:i_{loc,m} = i_{0,m}\left(\frac{c_\mathrm{Red}}{c_\mathrm{ref}} e^\frac{\alpha_{a,m} F \eta_m}{RT}-\frac{c_\mathrm{Ox}}{c_\mathrm{ref}}e^\frac{-\alpha_{c,m} F \eta_m }{RT}\right)

对于所有参与反应的物质而言,参考浓度 c_\mathrm{ref} 必须是相同的,这确保了过电位在零电流密度(平衡)状态下会遵守热力学 Nernst 方程。

下图展示了线电极示例的三次电流分布情况。由于具有浓度依赖性,三次电流分布会受电解质流动的影响,因此会受到质量传递控制的反应物的影响。当导线之间的流速很小时,为了获得法拉第电流而消耗的电解质不能得到有效补充,从而使电池内的这些部分形成了反应物的耗尽区。这大大降低了局部电流密度,这一现象可描述为“质量传递限制”,它进一步导致更多电流产生于线电极的外边缘。由于电流受到质量传递的限制,因此可以观察到电压降相应地增加,这就是所谓的“浓度过电位”。

三次电流分布
三次电流分布,Ecell = 1.65 V。正极上的电流密度分布(无量纲)。

何时使用三次电流分布 接口?

当电解质混合不均匀或者(与净电流密度相比)浓度相对较低时,可以使用三次电流分布对电池进行模拟,原因是电池内的电解质的组成变化很大,电阻损耗不遵守欧姆定律。在三次电流分布中,利用电流密度和电中性来求解所有 物质浓度的 Nernst-Planck 方程,会使方程组变得非线性和复杂化,致使仿真时间拉长,内存需求增加。在对三次电流分布进行模拟之前,最好利用二次电流分布来预测和判断电化学电池的可能性能。

还有哪些其他选项?

在分析电化学电池的电流-电压关系时,一次、二次和三次电流分布的近似程度依次增加。除此之外,还有其它的建模方法可以最大限度获取有关电池性能的信息,同时尽可能地减少模型的复杂度。

化学物质质量传递和二次电流分布

如果电流密度受到电活性物质质量传递的限制,但电解质组成基本不变,则可能不需要完整求解三次电流分布。相反,离子强度恒定不变意味着我们可以假设电解质溶液服从欧姆定律,具有恒定的电导率,因此可以利用二次电流分布 来求解电解质电位。然而,将动力学速率定律耦合到一个求解化学物质扩散问题的传递模型(必要时也可以求解迁移和对流问题)后,动力学速率定律变得与浓度相关。

事实上,这种耦合方法被用于线电极的三次电流分布案例,这是因为在此案例中是反应物的消耗过程,而非大量电解质物质起到了决定作用。这种耦合方法的另一个应用示例是橙子电池模型。此外,在分析电池和燃料电池时,将电荷传递与质量传递部分耦合是一种普遍的方法。

电分析

当惰性(支持)电解质数量远远超过反应(电活性)物质数量时,就会发生一种特殊情况。因此,与法拉第电流密度相比,溶液的离子强度很大。在这种情况下,电场非常小,因此电解质电位几乎处处相等——电解质溶液电阻对电化学电池的特性没有显著的影响。

当电极动力学(活化)和电活性物质的质量传递的影响远大于电解质溶液电阻的影响时,你可以使用电分析 接口。这个化学物质的传递接口可求解质量传递的对流扩散方程,它包含了电极动力学边界条件,以驱使电极-电解质界面的化学物质通量随局部过电位的变化而变化。

溶液电阻为零的电分析近似适用于电化学技术的标准实验装置,例如循环伏安法、计时安培分析法和电化学阻抗谱。在“案例下载”库中,你可以查看循环伏安法模型的案例,此案例便采用了这种近似方法。

选择模拟方法的流程图

总结思考

本篇博文介绍了 COMSOL Multiphysics 中四个电化学相关模块中都可用的三个电流分布接口,以及何时、为什么要使用它们。COMSOL Multiphysics 软件的强大之处在于,它能够模拟所有电流分布类型(一次、二次和三次),因此在分析电化学电池时,你可以灵活地逐步引入并掌控理论模型的复杂度。

如果你有兴趣使用 COMSOL Multiphysics 来模拟电化学电池设计,或者遇到了本文未提及的问题,请与我们联系

浏览电化学加载模块


评论 (4)

正在加载...
WYZFJ
WYZFJ
2024-05-19

I have a mystery I haven’t figured out. In general, the definition of electrode potential in textbooks or literature is more complex, but it is emphasized that the electrode potential is not equal to the potential difference between the electrode and the solution. However, I found that COMSOL’s definition of electrode potential is this. Does this definition cause the simulation results to be inconsistent with the objective physical facts? Thank you!

chong shi
chong shi
2024-05-24 COMSOL 员工

The electrode potential is the difference in potential between the solid-liquid interface, and it is not directly defined in COMSOL. Generally, only the solid-phase potential is defined, while the liquid-phase potential and lectrode potential are calculated based on interface dynamics and solution properties.

硕 佟
硕 佟
2024-06-24

请问二次电流分布+稀物质传递和三次电流分布有什么不同,另外如果电极上电流的分布与溶液电阻相关(电极由一个整流柜供电,整流柜到电极之间还有复杂回路),是不是要手动把电解质电导率写成离子浓度的方程才可以算,有没有其他的方法

chong shi
chong shi
2024-06-25 COMSOL 员工

二次电流分布+稀物质传递等价于“支持电解质-三次电流分布”,从方程形式上两者完全相同,不过前者需要手动做耦合,后者耦合已经自动完成。适用于,电解液电导率对离子浓度变化不敏感的场景,比如加入了大量支持电解质。如果电解质电导率与离子浓度有关,建议采用“电中性-三次电流分布”,它默认考虑电解质离子电导率与浓度有关,确切将与离子通量有关。

浏览 COMSOL 博客