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

L4 · 第 9 章

贫盐层什么时候会自己流动起来

Sandbox Physics ·

膜前的贫盐层一旦流动，液体就能把远处的盐带近，原先的扩散限制也会改变。但电流超过扩散极限，还不足以判定流动已经出现。要解释电对流，必须追踪一个横向小扰动怎样改变电场、推动液体，再由液体改变离子分布。

先让一个小扰动有地方生长

把第四章的通道向侧面展开。坐标 XX 仍从盐水库指向膜，YY 沿膜面增加，两者以米计。原来的静止解只有沿 XX 的变化。现在假设膜前有一小片区域比旁边稍浓：这一点浓度差就引入了沿膜面的变化。

较浓区域更容易导电，周围电势随之重新分配，可能出现横向电场。横向场作用于电荷，可以带起一小股液流。若这股流动又把更多盐带到最初较浓的位置，扰动得到放大；若它把浓度差摊平，扰动反而衰减。电对流失稳问的正是这个反馈的正负。

一维空间电荷计算提供静止的导电背景，允许电场在耗尽区中变强。它没有横向位置，也就没有形成循环、把侧面盐水卷进来的自由度。下一步要让横向扰动、流体运动和离子守恒一起演化，观察盐怎样重新分布。

对应阅读：Rubinstein & Zaltzman

在离子通量中加入液体的搬运

流速 u=(uX,uY)\mathbf u=(u_X,u_Y) 用每秒米，数密度 n±n_\pm 用每立方米。离子除了扩散和电迁移，还随液体运动；数通量增加一项数密度乘流速。符号 ∇\nabla 表示同时取两个空间方向的变化率，点积则把各方向的输运加起来。

液体很慢时，惯性可以忽略。压力 pp、黏性和单位体积电力相互平衡；η\eta 是黏度，单位为帕秒，ρe=e(n+−n−)\rho_e=e(n_+-n_-) 是电荷密度。压力项包含可吸收为梯度的平衡贡献。它会抵消某些电力，所以“有电荷、有电场”本身还不能确定一股环流。

∂τn±=−∇⋅N±,N±=n±u−D(∇n±±en±kBT∇Φ)\partial_\tau n_\pm=-\nabla\cdot\mathbf N_\pm,\qquad\mathbf N_\pm=n_\pm\mathbf u-D\left(\nabla n_\pm\pm\frac{en_\pm}{k_BT}\nabla\Phi\right)
这是等扩散系数、一价稀理想电解质的二维输运。无量纲化前先保留单位，可以直接检查每一项都是粒子数通量。
−∇p+η∇2u+ρeE=0,∇⋅u=0,E=−∇Φ-\nabla p+\eta\nabla^2\mathbf u+\rho_e\mathbf E=0,\quad\nabla\cdot\mathbf u=0,\quad\mathbf E=-\nabla\Phi
第二式表示不可压缩。离子分布仍与 Poisson 方程联立；固体边界上的无滑移或薄层有效滑移，须与所解析的电荷层范围一致。

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

算一次搬运盐的时间

设贫盐区厚度 ℓ=100  μm\ell=100\;\mathrm{\mu m}，离子扩散系数 D=10−9  m2 s−1D=10^{-9}\;\mathrm{m^2\,s^{-1}}。若已经观察到速度 U=20  μm s−1U=20\;\mathrm{\mu m\,s^{-1}} 的循环，液体跨过这一长度约需五秒，扩散则需十秒。两者之比定义 Péclet 数，衡量给定流动搬运浓度的能力。

本例 Pe=2\mathrm{Pe}=2，平流已可能显著改变贫盐区。若速度只有 5  μm s−15\;\mathrm{\mu m\,s^{-1}}，平流时间变为二十秒，Pe=0.5\mathrm{Pe}=0.5。这些是教学数量级估计；速度是题目给定的观测条件，没有由一维电流反推。

这也解释了一个常见误用：平流和扩散一样快，不等于此时必定失稳。Péclet 数使用的是已经存在的流速；要预测静止液体会不会自己动起来，还要判断小扰动所驱动的流是否反过来增强它。

τD=ℓ2D=10.0  s,τA=ℓU=5.00  s,Pe=UℓD=τDτA\tau_D=\frac{\ell^2}{D}=10.0\;\mathrm s,\qquad\tau_A=\frac\ell U=5.00\;\mathrm s,\qquad\mathrm{Pe}=\frac{U\ell}{D}=\frac{\tau_D}{\tau_A}
选择同一个长度比较平流与扩散。失稳起点则需由扰动增长率确定。
给定流速下的教学时间尺度估计
U/(μm s−1)U/(\mathrm{\mu m\,s^{-1}})τA/s\tau_A/\mathrm sPe\mathrm{Pe}
5.0020.00.500
10.010.01.00
20.05.002.00

从扰动的增长率读失稳曲线

设某个量原来只有静止剖面 f0(X)f_0(X)。在上面叠加很小的横向波纹，把它写成随时间变化的振幅与沿膜面重复的波形。用复指数只是把正弦和余弦一起记下来；最后取实部，得到可观测的实数。增长率 ss 的实部为正时，波纹越长越大。波数 kk 决定重复间距，单位为每米。

把这个形式代入输运和流动方程，只保留扰动振幅的一次项，就得到线性稳定性问题。对每个外加电压、每个波数，求出增长率；最快的模式刚好从衰减变成增长时，静止解失去稳定性。非线性演化再决定最后形成稳态涡旋、摆动还是更复杂的流动。

打开 Rubinstein 与 Zaltzman 的 arXiv:1602.00137v1 图 2。横轴是沿界面的波数，纵轴是电压，均用论文自己的无量纲尺度；每条曲线上方为不稳定。实线保留空间电荷，虚线取电中性极限。该图取论文的层宽比 L=1L=1、Debye 长度参数 ε=10−4\varepsilon=10^{-4}，并分别改变膜固定电荷参数 NN。比较曲线最低点前，先读清几何和边界。

f=f0(X)+ℜ ⁣[f^(X)esτ+ikY],λY=2πk,ℜ(s)=0f=f_0(X)+\Re\!\left[\widehat f(X)e^{s\tau+ikY}\right],\qquad\lambda_Y=\frac{2\pi}{k},\qquad\Re(s)=0
最后一个条件定义中性稳定曲线。例如 λY=50.0  μm\lambda_Y=50.0\;\mathrm{\mu m} 对应 k=1.26×105  m−1k=1.26\times10^5\;\mathrm{m^{-1}}。选 ℜ(s)=0.200  s−1\Re(s)=0.200\;\mathrm{s^{-1}} 作教学示例，五秒后小扰动振幅约增为 e1=2.72e^1=2.72 倍。实际增长率要由选定边界下的线性稳定性计算给出。

对应阅读：Rubinstein & Zaltzman

同样的超极限电流，可以由不同过程承载

先看带负电的窄通道壁。它周围需要额外阳离子屏蔽壁面电荷，即使中间液体已经很稀，这些阳离子仍能沿壁导电。通道变窄时，单位液体体积对应的壁面面积增大，表面贡献就更值得检查。这个过程可以增加电流，却不要求出现涡旋。

再把通道加宽，带电侧壁上的电渗流可能把盐水送到膜前，并在封闭端附近回流。更宽的通道还可能容许膜端贫盐层的横向扰动长成电渗失稳。Dydek 等人的 arXiv:1108.2325v1 图 1 把表面导电、侧壁电渗、膜端失稳分成三幅物理图。它们都是特定几何和壁面条件下的机制，不宜只凭一条电流曲线命名。

实验可以同时记录电流和示踪流速，改变通道高度、壁面电荷或盐浓度，看哪一种预测随之变化。若怀疑有水解离，还需要检测生成物或酸碱变化。把这些观测与电流曲线相互对照，就能逐项检验额外输运的来源。

对应阅读：Dydek et al.

前沿模型怎样接受检验

Rubinstein 与 Zaltzman 的工作已发表于 2016 年 Physical Review Fluids；这里按其固定的预印本版本读图。它用给定膜选择性和边界的线性稳定性计算比较不同驱动项。Dydek 等人的微通道研究发表于 2011 年 Physical Review Letters，结合简化模型和二维数值结果讨论几何怎样选择输运机制。两者都应按各自的计算问题阅读。

一个下一步计算可以从已验证的一维导电解开始：逐个横向波长加入很小扰动，测量它早期的增长率，再分别加密电荷层网格与时间步。改变扰动幅度时，线性阶段的增长率应保持一致；改变侧面边界时，也要重新判断允许哪些波长。这样才能把静止导电、失稳起点和后来的流动形态分开检验。

强驱动还可能让界面附近的离子拥挤。Kilic、Bazant 与 Ajdari 的有限离子体积模型说明，点离子理论对电容和盐吸附的外推会改变。若问题的关键已经是拥挤、浓溶液活度或多孔电极内的反应分布，就应另立相应模型。当前的一维输运和弱场圆柱电渗实验，提供的是这条研究链中可独立复算的起点。

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

动手检验预测

确定静止导电计算已经回答了什么

  1. 预测增加电压会先改变浓度、电场还是横向流速；说明各个量能否由一维模型给出。
  2. 比较“扩散来得及补充”和“空间电荷扩展”，保存同一边界条件下的轨迹，记录晚期电流与最低盐浓度。
  3. 为一项后续二维计算写出还需增加的量：横向波长、流速场、侧面与膜面的流体边界，再说明应如何测增长率。

区分给定的电渗流和自发失稳

  1. 观察固定圆柱周围的稳态流线，先预测反转外场会怎样改变它。
  2. 把外场反向并比较流动方向，再旋转外场改变四瓣结构的朝向。
  3. 解释这个由固定表面和给定电场确定的解析流场，与贫盐层中小扰动自行增长的问题有何不同。

练习与解释性答案

  1. 一条电流曲线越过经典扩散极限，但显微镜没有分辨出涡旋。还应考虑什么？

    展开解释

    可以先检查一维空间电荷输运、沿带电壁的表面导电以及化学生成新离子等过程。未分辨出流动也可能受成像灵敏度限制。比较几何与壁电荷依赖、提高流速测量灵敏度，并按需要测酸碱或反应产物，才能逐项区分。

  2. 保持流速与扩散系数不变，把输运长度减半。扩散时间、平流时间和 Péclet 数各怎样改变？这能否确定失稳电压？

    展开解释

    扩散时间变为四分之一，平流时间减半，Péclet 数也减半。扩散相对更有机会摊平浓度差。但失稳电压还取决于扰动如何产生流动、膜选择性和流体边界；这些量不能只用给定流速的时间尺度代替。

  3. 某扰动的振幅在四秒内从一个单位增长到两个单位。若它处在线性指数增长阶段，增长率是多少？为什么还要换一个初始振幅再做一次？

    展开解释

    增长率为 ℜ(s)=ln⁡2/(4  s)=0.173  s−1\Re(s)=\ln2/(4\;\mathrm s)=0.173\;\mathrm{s^{-1}}。换更小的初始振幅应得到相同早期增长率；若增长率随振幅明显改变，可能已经进入非线性阶段，或数值误差尚未受控。

来源与计算范围

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

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

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

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

  3. Kilic, Bazant & Ajdari, Steric effects in the dynamics of electrolytes at large applied voltages I, arXiv:physics/0611030v1 (2006)

    第 I—II 节解释点离子稀溶液理论为何可能高估表面浓度，并介绍有限尺寸的替代模型。这些扩展是进阶阅读，尚未加入浏览器的理想 PNP 求解器。

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

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