电池中的迁移、扩散离子通量分析(仿真视角)

P2D 模型中电解质离子通量的物理含义与守恒性

在电化学 P2D(伪二维)模型中,电解质相的锂离子通量 $\mathbf{N}_+$ 通常由三个部分组成。这三部分共同构成了电解质中的总物质通量,并严格遵循质量守恒定律。

1. 离子通量的三项分解(修正版)

根据 Newman 的浓溶液理论,电解质中的锂离子通量 $\mathbf{N}_+$(单位:$\text{mol} \cdot \text{m}^{-2} \cdot \text{s}^{-1}$)的标准表达式如下:

\[\mathbf{N}_+ = \underbrace{-\frac{t_+}{F}\kappa_{\text{eff}} \nabla \phi_e}_{\text{1. 电场驱动项}} + \underbrace{\frac{2RT}{F^2} t_+(1-t_+)\left(1+\frac{d\ln f}{d\ln c_e}\right)\kappa_{\text{eff}} \nabla \ln c_e}_{\text{2. 浓差过电势驱动项}} - \underbrace{D_{\text{eff}} \nabla c_e}_{\text{3. Fick 扩散项}}\]

1.1 电场驱动项

\[\mathbf{N}_{\text{elec}} = -\frac{t_+}{F}\kappa_{\text{eff}} \nabla \phi_e\]
  • 物理意义
    • 这是由于电解液内部的 电势梯度(电场) 引起的离子迁移。
    • 根据欧姆定律,电流密度 $\mathbf{i} = -\kappa_{\text{eff}} \nabla \phi_e$。由于只有部分电流由锂离子承载(份额为迁移数 $t_+$),因此锂离子通量正比于 $t_+ \mathbf{i}$。
    • 符号解释:放电时,电势梯度 $\nabla \phi_e$ 通常为负(电势从负极到正极降低),因此 $-\nabla \phi_e$ 为正,表示离子从负极向正极迁移。

1.2 浓差过电势驱动项

\[\mathbf{N}_{\text{conc}} = \frac{2RT}{F^2} t_+(1-t_+)\left(1+\frac{d\ln f}{d\ln c_e}\right)\kappa_{\text{eff}} \nabla \ln c_e\]
  • 物理意义
    • 这一项本质上是 由浓度梯度引起的附加离子迁移
    • 当电解液中存在浓度梯度时,会产生浓差电势。虽然它表现为电势梯度,但其根源是浓度分布不均。
    • 该项来源于修正的欧姆定律中的浓差电势项,必须乘以 $t_+/F$ 转换为物质通量
    • 关键系数:$t_+(1-t_+)$ 反映了阴阳离子迁移数的差异。如果 $t_+ = 1$(锂离子承载所有电流),则此项为零,因为此时浓度梯度不会产生液接电势来驱动额外迁移。

1.3 Fick 扩散项

\[\mathbf{N}_{\text{diff}} = -D_{\text{eff}} \nabla c_e\]
  • 物理意义
    • 这是经典的分子扩散项,由浓度梯度直接驱动,与电场无关。
    • $D_{\text{eff}}$ 为考虑多孔电极孔隙率和曲折度后的有效扩散系数(通常通过 Bruggeman 关系修正:$D_{\text{eff}} = D \cdot \varepsilon^{1.5}$,或者通过孔隙迂曲度校正)。
    • 它的作用是“抹平”浓度梯度,使离子从高浓度向低浓度扩散。

2. 总通量与物质守恒

2.1 总物质通量

将上述三项相加,即得电解质中锂离子的总通量: \(\mathbf{N}_{\text{total}} = \mathbf{N}_{\text{elec}} + \mathbf{N}_{\text{conc}} + \mathbf{N}_{\text{diff}}\) 这完整地描述了单位时间内通过单位面积的锂离子总物质的量。

2.2 质量守恒

该通量表达式严格遵循质量守恒定律。

在 P2D 模型中,电解质相的质量守恒方程(连续性方程)为:

\[\frac{\partial (\varepsilon_e c_e)}{\partial t} = -\nabla \cdot \mathbf{N}_{\text{total}} + s_{\text{rxn}}\]

其中:

  • $\varepsilon_e$:电解质体积分数(孔隙率)。
  • $s_{\text{rxn}}$:电化学反应源项(单位体积的锂离子生成/消耗速率)。

守恒性的具体体现:

  1. 在隔膜区域
    • 没有电化学反应发生,$s_{\text{rxn}} = 0$。
    • 此时方程变为 $\nabla \cdot \mathbf{N}_{\text{total}} = 0$。这意味着通量是守恒的,即进入隔膜一侧的锂离子量等于流出另一侧的量(对于稳态)。
  2. 在多孔电极区域
    • 存在电化学反应(嵌入/脱出)。
    • 此时 $\nabla \cdot \mathbf{N}_{\text{total}} \neq 0$。通量的散度(流入与流出之差)正好等于电化学反应消耗或生成的锂离子量。这体现了源头守恒:反应消耗了多少锂离子,通量就减少多少

3. 总结

通量项 驱动力 符号特征 典型作用
电场驱动 电势梯度 $\nabla \phi_e$ $-\nabla \phi_e$ 放电时从负极指向正极,主要贡献项
浓差过电势 浓度梯度 $\nabla \ln c_e$ 系数含 $t_+(1-t_+)$ 高倍率下显著,修正因离子迁移数差异引起的扩散势
Fick 扩散 浓度梯度 $\nabla c_e$ $-\nabla c_e$ 从高浓度指向低浓度,缓解浓差极化