电解模拟
理论:气泡降低电解液电导率
概述
电解仿真中除了电极反应,还需要计算什么?
气泡。水电解中阴极产生氢、阳极产生氧,以气泡形式出现。气泡不导电,电极周围的电解液混入气泡后,离子能通过的截面减少,电压损失增大。而且气泡向上运动,电极下部和上部气泡量不同,电压损失也随位置变化。损失大的地方电流难以流过,电流偏向电极下部,下部又产生更多气泡。也就是说电化学反应、电解液与气泡的流动、电流分布相互影响。把这些耦合求解是电解仿真的核心。
气泡分数与有效电导率
j_g和j_l为气体和液体的表观流速,u_d为气泡比液体上升快的速度(0.20 m/s),V_m为80℃、1个大气压下的气体摩尔体积,n为每摩尔气体的电子数(氢2、氧4),d为电极与隔膜的间隙(3 mm)。κ为无气泡电解液的电导率(1.2 S/cm,示例值)。
间隙电压损失
“电解”这个词
法拉第在1833〜1834年证明了电解产生的物质量与电量成正比。当时他与学者休厄尔商量,确定了电极、阳极、阴极、电解质、离子等至今仍在使用的术语,电解(electrolysis)也是其中之一。用电分离物质的思想随着术语一起得到整理,推动了此后电化学工业的发展。
算例
例1:电极高度与气泡、电压损失(0.3 A/cm²,高1 m)
| 距下端高度 | 阴极侧气泡 | 阳极侧气泡 | 间隙损失 | 气泡引起的增量 |
|---|---|---|---|---|
| 0 m | 0% | 0% | 0.150 V | 0 |
| 0.25 m | 11% | 6% | 0.172 V | +0.022 V |
| 0.5 m | 20% | 11% | 0.195 V | +0.045 V |
| 1.0 m | 34% | 20% | 0.244 V | +0.094 V |
沿高度平均,间隙损失为0.196 V(比无气泡多0.046 V)。电压1 V相当于每kg氢26.6 kWh,所以整个间隙为5.2 kWh/kg,其中气泡部分为1.2 kWh/kg。这里假设电流均匀。
例2:有间隙与零间隙
| 结构 | 间隙损失(平均) | 每kg氢 |
|---|---|---|
| 3 mm间隙·无气泡(假想) | 0.150 V | 4.0 kWh |
| 3 mm间隙·有气泡 | 0.196 V | 5.2 kWh |
| 零间隙 | 约0(隔膜电阻仍在) | 约0 |
做成零间隙气泡问题就消失了吗?
间隙电解液的损失没有了,但气泡问题换个地方仍然存在。零间隙中气泡穿过带孔电极逃到背面,但气泡附着在电极表面或孔内时,可反应面积减少,电流集中到其余部分。背面流道积聚的气泡难以排出时,电解液难以到达电极。所以电极孔的大小形状、背面流道宽度、电解液流速都会改变气泡的排出,研究这些就是仿真的任务。提高电流密度气泡量也增加,要确认电极上部和出口附近气泡不会过多。气泡分数高还会影响氢和氧穿过隔膜混合的量,从安全角度也需要掌握气泡分布。
建模流程
- 把电解液电导率作为温度和浓度的函数给出。
- 用极化曲线表示电极反应,求解电解液中的电位。
- 由电流分布求产气量,求解气泡与液体的流动(两相流)。
- 按气泡分数更新电导率,重新求解电流分布。
- 由发热求温度,反映到电导率和反应中直到收敛。
“只有上半部发黑的电极”
拆开某碱性水电解装置,发现只有电极上半部变色,催化层变薄。把气泡与电流耦合分析表明,上部气泡积聚、电解液难以到达,局部电压长时间偏高。增加电解液循环、加宽电极背面流道使气泡易于排出后,上下气泡差变小,此后检查未再发现变色。
常见错误
错误与对策
| 错误 | 后果 | 对策 |
|---|---|---|
| 忽略气泡 | 电压低估 | 按气泡分数降低电导率 |
| 气泡与液体同速 | 气泡分数高估 | 计入气泡上升速度 |
| 电流均匀 | 漏掉上部劣化 | 求解电流分布 |
| 氢氧气泡等量 | 阴极侧低估 | 体积为2比1 |
| 温度取常数 | 电导率出错 | 与热耦合 |
我想了解相关内容。
相关页面参考水电解(质子交换膜型)、电解加工、电镀、腐蚀仿真和燃料电池建模。
详细
报告