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

L4 · 第 8 章

从原论文重算圆柱周围的诱导电荷流

Sandbox Physics ·

打开 Squires 与 Bazant 的《Induced-charge electro-osmosis》，选 arXiv:physics/0304090v2。图 3(a) 的四瓣流动，可以从一个圆柱外的电势一路算出来。先分清图上的电荷属于金属还是液体，再检查滑移方向，最后用表 1 的速度场重画流线。做完这三步，就能判断自己的图是否真的满足这篇论文的模型。

先把两个标题和两个版本放对

JFM 长文《Induced-charge electro-osmosis》的选定版本是 physics/0304090v2，提交修订日期为 2004 年 2 月 9 日，后来发表于 Journal of Fluid Mechanics 第 509 卷，217–252 页。这里使用 PDF 中第 3.1 节的公式编号和图 3(a)；部分 HTML 页面会重新连续编号，找不到式号时，用公式内容和章节定位。

短文《Induced-charge Electrokinetic Phenomena: Theory and Microfluidic Applications》对应 physics/0306100v2，修订日期为 2004 年 2 月 8 日，发表于 Physical Review Letters 92，066101。它的图 1(c) 画出同一类中性圆柱流动。两篇文章可以互相帮助阅读，复算时仍要写清用的是哪一版、哪条公式。

这次只重算稀电解质、薄双电层、弱外场下的固定中性圆柱直流稳态。论文还讨论交流、充电瞬态、带净电荷的圆柱和表面涂层；这些内容各有自己的条件，不随这一小段复算一起完成。

对应阅读：Squires & Bazant · Bazant & Squires

图 3(a) 先读电荷，再读箭头

把图中的外场方向当成极角 θ=0\theta=0。半径 aa 是圆心到金属表面的距离，径向坐标 rr 向液体外侧增加。金属固定、总电荷为零；“中性”指一周积分后的净电荷为零，表面各处仍能有正负电荷。

向着外场的一极，金属表面电势相对附近体相较高，金属带正电，旁边的弥散离子层带负电。另一极的符号相反。图中的正负号如果没有先分清属于哪一侧，就很容易把流动方向画反。

紧接着看表面附近的箭头：流体从两极朝赤道移动，外层则在两极方向流入、在赤道方向流出。图没有频率轴、时间轴或实验误差棒；它展示理论模型的稳态流动结构。定量比较还要检查流速、边界和守恒。

对应阅读：Squires & Bazant

沿同一个电势检查三个符号

PDF 式 (3.4) 给出充电后的外部电势。取金属电势零，远处外场沿正轴，极角逆时针增加。先把 r=ar=a 代入外部电势，再用金属电势减去这个值，得到式 (3.5) 的 ζ\zeta。沿圆周求导则得到式 (3.6) 的切向场。

在 θ=45∘\theta=45^\circ，正的 zeta 电势对应负的弥散电荷，切向电场也是负的。电场给负离子的切向力为正，所以液体沿角度增加的方向流动，朝向赤道。把结果代入滑移关系，应该得到正的切向速度。

ϕ=−E0(r+a2r)cos⁡θ,ζ=2E0acos⁡θ,Eθ(a)=−2E0sin⁡θ\phi=-E_0\left(r+\frac{a^2}{r}\right)\cos\theta,\quad\zeta=2E_0a\cos\theta,\quad E_\theta(a)=-2E_0\sin\theta
三项来自同一个外部电势。外部体相近似电中性，圆柱外缘的法向电流为零；金属内部仍是等势体。
qd=−εζλD,us=−εζEθη=2U0sin⁡(2θ),U0=εaE02ηq_d=-\frac{\varepsilon\zeta}{\lambda_D},\quad u_s=-\frac{\varepsilon\zeta E_\theta}{\eta}=2U_0\sin(2\theta),\quad U_0=\frac{\varepsilon aE_0^2}{\eta}
这里 qdq_d 是弥散离子电荷面密度，ε\varepsilon 是液体介电常数，η\eta 是动力黏度。PDF 式 (3.7)–(3.8) 给出滑移及速度尺度。

对应阅读：Squires & Bazant

先算三个数，再画一条流线

用 a=10.0 μma=10.0\,\mu\mathrm m、E0=0.450 V mm−1E_0=0.450\,\mathrm{V\,mm^{-1}}，以及 ε=6.95×10−10 F m−1\varepsilon=6.95\times10^{-10}\,\mathrm{F\,m^{-1}}、η=10−3 Pa s\eta=10^{-3}\,\mathrm{Pa\,s}。先把场强换成 450 V m−1450\,\mathrm{V\,m^{-1}}、半径换成米，得到 U0=1.407375 μm s−1U_0=1.407375\,\mu\mathrm{m\,s^{-1}}。

表面四十五度处的速度应为 2U02U_0。再移到离圆心两倍半径的两点，用 PDF 式 (3.9) 或表 1 计算。下表用 U0U_0 归一化，数值不依赖所选水参数，适合先手算再检查程序。

ur=2U0a(a2−r2)r3cos⁡(2θ),uθ=2U0a3r3sin⁡(2θ)u_r=2U_0\frac{a(a^2-r^2)}{r^3}\cos(2\theta),\qquad u_\theta=2U_0\frac{a^3}{r^3}\sin(2\theta)
径向与切向分量共同决定一条流线。只把表面滑移向外复制，会漏掉不可压缩流动所要求的径向分量。
中性圆柱速度场的三个特殊位置；数值是解析模型的重算参考。
位置径向速度 ur/U0u_r/U_0切向速度 uθ/U0u_\theta/U_0
r=a, θ=45∘r=a,\ \theta=45^\circ0022
r=2a, θ=0r=2a,\ \theta=0−3/4-3/400
r=2a, θ=45∘r=2a,\ \theta=45^\circ001/41/4

对应阅读：Squires & Bazant

把复算保存下来，检查边界和守恒

计算参考文件 methods-iceo-reference.json 保存了 181181 个表面角度和 16811681 个速度网格点。表面数据列包含角度、zeta 电势、切向电场、弥散面电荷和滑移；网格列中的 x,yx,y 用微米，u,vu,v 用微米每秒。标记 insideMetal 的点属于固体，不用于外部流动积分。

文件中的参数、单位和来源版本应当跟着数据一起保留。它是我们从解析模型重算的数值，论文里的理论图是阅读目标；两者都不含金属丝实验的测量噪声。若要与实验比较，还需要独立的装置尺寸、表面状态和采集条件。

先检查圆柱表面法向速度为零，切向速度等于电势推得的滑移；再对圆周积分，确认净体积通量与总弥散电荷接近零。把速度作有限差分，检查一个小格的净流出是否为零，这就是散度检查。再计算描述局部旋转的涡量，检查它在两个方向的二阶导数之和是否为零；这检验消去压强后的稳态 Stokes 方程。最后才积分示踪轨迹，重画图 3(a) 的四瓣结构。

默认 λD=0.100 μm\lambda_D=0.100\,\mu\mathrm m 时，薄层比为 0.01000.0100，最大诱导电势为 9.00 mV9.00\,\mathrm{mV}。这些尺度支持有限的弱场薄层演示。继续增大场强时，先检查电势与热电压的比值；平方关系向强场的外推，需要重新检验电荷模型。

对应阅读：Squires & Bazant · Sandbox Physics analytic ICEO reference

下一篇论文也可以这样读

先固定版本与图号，写明哪些曲线是理论、模拟或测量。接着从图注找几何和坐标，从正文找边界与近似。挑出两三个可以手算的极限，等它们通过，再画完整曲线。保存参数与中间量，下一次比较就能追溯到具体假设。

如果新论文加入电极反应、壁面或浓度变化，先找出原来哪条边界或守恒关系被改写，再决定哪些公式还能沿用。一个复算片段的价值，在于让你知道下一次该重新检查哪一步。

动手检验预测

用三个特殊位置检验四瓣流动

  1. 在“四瓣电渗流”下下载计算剖面，按单位读取表面角度和速度网格。
  2. 先检查四十五度的表面滑移，再检查两倍半径处的径向和切向速度。
  3. 把电场反向：比较电荷符号、速度数值和流线，说明哪些项反号、哪些项保持相同。

练习与解释性答案

  1. 有人把图 3(a) 画成了四瓣结构，但圆柱表面存在向内的径向速度。外观相似，是否足以接受复算？

    展开解释

    不够。固定、不可穿透的圆柱要求表面法向速度为零。该图违反了装置的边界条件；应先检查极坐标分量、半径归一化和坐标转换，再比较流线外观。

  2. 把场强加倍、半径减半，最大诱导电势和速度尺度怎样变？这能否自动保证模型仍然适用？

    展开解释

    最大诱导电势与场强、半径的乘积成正比，因此不变；速度尺度与半径和场强平方成正比，因此加倍。若 Debye 长度不变，薄层长度比也加倍，薄层近似变差。电势条件相同，还需要重新检查几何尺度条件。

  3. 计算文件没有误差棒，能否据此说原论文的实验误差为零？

    展开解释

    不能。这份文件保存解析模型的计算值，选定的论文图也是理论流动图；其中没有实验随机误差的估计。数值精度要靠边界、守恒和差分误差检查，真实实验的不确定度需要另有采集数据。

来源与计算范围

  1. Squires & Bazant, Induced-charge electro-osmosis, arXiv:physics/0304090v2 (9 February 2004)

    固定 Journal of Fluid Mechanics 509, 217–252 的预印本版本：第 3.1 节、PDF 式 (3.4)—(3.9)、图 3(a) 与表 1。限定复算对象是固定中性圆柱的稳态电势、电荷、滑移和外部 Stokes 流。HTML 转换版采用连续公式编号。

  2. Bazant & Squires, Induced-charge Electrokinetic Phenomena: Theory and Microfluidic Applications, arXiv:physics/0306100v2 (2004)

    Physical Review Letters 92, 066101，主要对应式 (1) 与图 1(c)。这篇短文提供物理图景，与较长的 JFM 论文是两篇不同文献。原图仅链接，不重新分发。

  3. Sandbox Physics analytic ICEO reference, version 1 (original calculated values)

    由上述圆柱解析解计算的原创表面与速度采样，附单位和参数。它们是理论计算值，不是实验测量，也没有从论文图片取点。

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

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