跳到主要内容
Sandbox Physics
电化学与离子输运

L3 · 第 4 章

离子来不及时,电流怎样接近极限

选择性界面不断抽走正离子,附近的盐就会变稀。浓度越低,同样电流需要越大的电场。电中性近似给出一个扩散极限;继续加压时,空间电荷层会变厚,输运规律也随之改变。

从盐水库走到只通正离子的界面

沿 x∈[0,1]x\in[0,1] 从左侧盐水库走向右侧选择性界面。两种浓度按盐水库初值归一化,电势按热电压归一化。左端取 c+=c−=1c_+=c_-=1、ϕ=0\phi=0;右端取 c+=1c_+=1、J−=0J_-=0、ϕ=−V\phi=-V。

右端固定的正离子浓度,代表一个规定正离子活度的理想选择性边界。它与左端恒浓度盐水库共同维持稳态输运。这个半电池模型没有解析真实膜的厚度与固定电荷分布;电压包括靠近右端的界面电荷层。

对应阅读:MIT 10.626 · Rubinstein & Zaltzman

电中性体相里,两个通量条件足够推出浓度

远离薄界面层时,先近似 c+=c−=cc_+=c_-=c。负离子不能穿过右端,稳态下它在任何截面的净通量都为零。令负离子的 Nernst–Planck 通量为零,得到电势斜率等于浓度的对数斜率。

再把这个斜率放进正离子通量。扩散与电迁移各贡献一半,合起来得到 J+=−2c′J_+=-2c\prime。把恒定正离子通量记为 jj,从左端的浓度一积分,浓度便是一条下降直线。

ϕ=ln⁡c,j=−2∂xc,c(x)=1−j2x\phi=\ln c,\qquad j=-2\partial_xc,\qquad c(x)=1-\frac j2x
这里无量纲电流 j=J+−J−j=J_+-J_-,因 J−=0J_-=0 而等于正离子通量。对应有量纲电流密度尺度是 eDn0/LeDn_0/L。

对应阅读:MIT 10.626

把右端的薄层接上,总电压会多一段

记紧邻界面层外侧的体相浓度为 csc_s,体相电势为 ϕs=ln⁡cs\phi_s=\ln c_s。体相浓度一般不等于边界上规定的正离子浓度一;薄层负责把两者连接起来。

在薄层近似局部平衡时,正离子电化学势连续。体相一侧是 ln⁡cs+ϕs\ln c_s+\phi_s,边界一侧是 ln⁡1−V\ln1-V。两者相等,得 2ln⁡cs=−V2\ln c_s=-V。一半对数电压降在稀释体相,另一半在浓度分配产生的 Donnan 层。

cs=e−V/2,jEN(V)=2(1−e−V/2),jlim=2c_s=e^{-V/2},\qquad j_{\rm EN}(V)=2\left(1-e^{-V/2}\right),\qquad j_{\rm lim}=2
这是该边界约定的电中性外解。若忽略界面电势差,把右端的边界值直接套到体相,会得到不同且不适用这里的电流—电压关系。

对应阅读:MIT 10.626 · Rubinstein & Zaltzman

手算一条接近极限的曲线

令 V=2V=2,室温对应总电压约 51.4  mV51.4\;\mathrm{mV}。界面外侧的盐浓度降到盐水库的 e−1=0.368e^{-1}=0.368,电流为 1.2641.264,达到扩散极限的约百分之六十三。

把电压再加到 V=6V=6,盐浓度只剩 e−3=0.0498e^{-3}=0.0498,电流增到 1.9001.900。电压提高了三倍,电流只增约一半。下表完全由电中性解析式计算,展示平台形成的原因。

电中性、薄界面层解析参考;电流按同一扩散尺度归一化
VVcsc_sjENj_{\rm EN}jEN/jlimj_{\rm EN}/j_{\rm lim}
10.6070.7870.393
20.3681.2640.632
40.1351.7290.865
60.04981.9000.950

浓度接近零时,电中性近似先遇到困难

直线浓度解若想承载 j>2j>2,右端附近就会出现负浓度。问题出在近似适用范围:离子越少,屏蔽长度越长,电势的曲率和两种离子之差开始重要。保留 Poisson 方程,可以形成扩展空间电荷,承担电中性解无法描述的一部分电压。

有限 Debye 长度的一维 PNP 解因而可以显示超过电中性扩散极限的电流。这种解仍然没有横向流动。真实装置的超极限电流还可能来自表面导电、侧壁电渗或电对流,要结合几何、场强和观察到的流动判断。

模型比较应先在离子尚未严重耗尽、薄层分离良好的条件下检查两者靠近,再向高电压走。高驱动下的曲线分离,还要用网格加密检查电荷层是否解析充分,才能归因于物理。

对应阅读:Rubinstein & Zaltzman · Dydek et al.

动手检验预测

沿电流平台追踪离子耗尽

  1. 先用较小电压比较 PNP 与电中性参考,预测哪段浓度最先下降。
  2. 逐步加压,记录最低盐浓度、电流和空间电荷;在剖面中分清界面外侧体相与薄层。
  3. 减小 Debye 长度比,判断低驱动下两解是否更接近,再解释高驱动下为什么必须保留电荷。

练习与解释性答案

  1. 电流达到电中性扩散极限的九成时,界面外侧浓度是多少?室温总电压约多少?

    展开解释

    由直线浓度解,cs=0.1c_s=0.1。由薄层匹配,V=−2ln⁡0.1=4.605V=-2\ln0.1=4.605,室温总电压约 118  mV118\;\mathrm{mV}。这是模型内的参考,不含真实膜、接触或外电路的额外压降。

  2. 一个一维程序给出了超极限电流。能据此给电对流的起始电压吗?

    展开解释

    不能。电对流需要横向电场、空间电荷和流速扰动相互放大。一维静止流体模型没有这些自由度。它能给出静止输运背景,稳定性分析还要加入至少一个横向方向及流体方程。

来源与计算范围

  1. MIT 10.626, Ion concentration polarization, 2011 Lecture 23

    选择性输运、浓度梯度与扩散极限电流。课程另行定义盐水库至界面的正离子活度边界,并把 Donnan 电势差计入总电压。

  2. Rubinstein & Zaltzman, The driving factors of electro-convective instability in concentration polarization, arXiv:1602.00137v1 (2016)

    第 I—II 节、式 (3)—(9) 与图 1—3 区分选择性导电、流体耦合和横向稳定性计算。课程讲读这项已发表研究;一维 PNP 实验没有复现中性稳定曲线。

  3. Dydek et al., Overlimiting current in a microchannel, arXiv:1108.2325v1 (2011)

    图 1—3 及其中分别讨论的表面导电、侧壁电渗与失稳机制。论文的通道几何和电压约定与课程半电池不同;本课程 PNP 曲线没有复现这些图。

反馈这个页面有错误、看不明白,或者缺了什么?告诉我们
反馈类型

提交时会附带本页地址、Lab ID、网站语言和设备类型。请不要填写个人信息。反馈最多保存 12 个月。隐私政策