电化学 · 选择性界面前的盐逐渐耗尽

浓差极化与空间电荷

膜吸走离子的速度快过扩散补给时，电流会怎样变化？ 选择性膜只让一种离子通过，膜前便逐渐形成贫盐层。电中性近似给出一个扩散限制电流；耗尽足够强时，局部净电荷不可忽略，完整输运模型就会偏离这个极限。

让阳离子流向理想的选择性膜，观察贫盐区逐渐形成。把最终剖面与电中性近似并排比较，再增大电压，看看空间电荷怎样使电流越过经典输运极限。

加压后的时间 $0.120\,\mathrm{ms}$

电流密度大小 $0.0200\,j_0$

最低盐浓度 $0.120\,c_0$

电流与经典极限之比 $0.500$

一次加压，一组演化中的解

让膜吸走离子

施加电压

V/(k_{\rm B}T/e)

4

屏蔽长度与间距之比

\lambda_{\rm D}/L

0.06

拖动池底的发光标记，回看任一时刻。彩色粒子呈现计算出的离子密度。

电压阶跃正在改变离子分布。

单位与数值精度

L=2\,\mu\mathrm m

D=10^{-9}\,\mathrm{m^2\,s^{-1}}

T=298\,\mathrm K

t_0=L^2/D=4.00\,\mathrm{ms}

V_T=k_{\rm B}T/e

V_T=25.7\,\mathrm{mV}

改变屏蔽长度相当于改变初始盐浓度；溶剂按相对介电常数为 78.5 的水计算。有限体积解采用 112 个单元，每步以保持正浓度的隐式离子通量和 Poisson 方程联立迭代。

\Delta j_{\rm tot}<10^{-5}\,j_0

c_0=6.44\,\mu\mathrm{mol\,L^{-1}}

N_+/N_{+,0}=1.00

N_-/N_{-,0}=1.00

保存文件含物理单位、各时刻剖面、边界通量、极板电荷、离子总量与总电流的空间差异。理想电压切换瞬间的脉冲不计入记录。

模型与近似范围

储液池固定两种离子的浓度。膜把阳离子浓度保持为储液池值，阻挡阴离子，电势为施加电压的负值。电中性虚线是晚期外层参考，不含膜边电荷层。

稀溶液点离子模型省略分子堆积与关联。电压更高时，可能需要有限离子体积、电极反应或浓溶液模型；这里不计算流体运动。

模型假设

一维模型采用稀理想电解质和相同的离子扩散系数。左侧储液池保持浓度与电势；右侧界面固定阳离子浓度和电势、阻挡阴离子。施加电压包含界面 Donnan 电势差。两微米间距与放大的屏蔽长度用于教学。模型不含流动、电对流失稳、表面导电、水解离和真实膜微结构。

模拟输运记录

把离子分布和电流连起来看

穿过电池的离子浓度

穿过电池的电势

加压后的电流

阳离子阴离子电势／总电流虚线：电中性稳态参考

电中性参考渐近于电流密度单位的两倍。完整模型包括带电界面层及其瞬态；图中均为模拟或解析模型结果。

方程与原始来源

J_\pm=-\partial_xc_\pm\mp c_\pm\partial_x\phi,\quad -2\delta^2\partial_x^2\phi=c_+-c_-,\quad j_{\rm tot}=J_+-J_-+2\delta^2\partial_tE

Bazant、Thornton 与 Ajdari（2004），阻挡电池的充电动力学

Rubinstein 与 Zaltzman（2016），浓差极化与电对流

继续学习扩散、电荷与流动的联系

先抓重点

核心问题: 膜吸走离子的速度快过扩散补给时，电流会怎样变化？
一句话直觉: 贫盐浓度与空间电荷电场需要一起求解。一维电流超过电中性扩散极限，并不能用来判定流体什么时候失稳。

核心数学模型

扩散与电场漂移

J_\pm=-\partial_xc_\pm\mp c_\pm\partial_x\phi,\qquad \partial_tc_\pm=-\partial_xJ_\pm

距离以完整间距为单位，浓度以储液池浓度为单位，电势以热电压为单位，时间以间距平方除以扩散系数为单位。通量向右为正，两种离子的扩散系数相同。

电荷改变电势

-2\delta^2\partial_x^2\phi=c_+-c_-,\qquad \delta=\lambda_{\rm D}/L

两种浓度之差给出这套单位中的净电荷。屏蔽长度相对间距越小，电荷越容易集中在窄层里。数值网格要解析这些窄层，收敛性另作检查。

电中性稳态参考

c(x)=1-\frac{j x}{2},\quad \phi(x)=\ln c(x),\quad j=2(1-e^{-V/2})

左侧储液池浓度为一、电势为零。右侧膜固定阳离子浓度为一、电势为外加电压的负值，并阻挡阴离子。电压包含 Donnan 电势差，因此电中性外层的电流渐近于二。

常见难点

电中性极限有适用范围

容易误解经典极限电流就是所有电解质体系的绝对上限。

正确理解电中性扩散极限假设局部近似中性，且没有其他输运机制。空间电荷会改变电场漂移；流动与表面输运则是另外的拓展。

动手实验

01
先看扩散足够快时
选“扩散来得及补充”，等演化接近稳定。先猜正负离子浓度会在哪里几乎重合，再与虚线的电中性参考比较。
你会看到： 离开膜边的电荷层后，两种离子浓度几乎重合，并朝膜方向下降。参考线是稳态外层解，早期瞬态剖面不必与它相同。
02
追踪贫盐区
选“贫盐层形成”，拖动池底的时间标记，比较早期和晚期。随后增大“施加电压”。
你会看到： 膜前的盐逐渐耗尽，电势降也越来越集中在这里。离子通量和电流都由这同一组浓度、电势数组计算。
03
查看空间电荷
选“空间电荷扩展”，把完整电流与经典极限相比，再查看正负离子浓度曲线分开的位置。
你会看到： 带净电荷的区域让完整的一维模型可以承载超过电中性极限的电流。若要描述流动、表面导电或水解离，还需要加入相应物理过程。

浓差极化与空间电荷

把离子分布和电流连起来看

穿过电池的离子浓度

穿过电池的电势

加压后的电流