3. 一次电流分布仿真

一次电流分布（Primary Current Distribution）是电化学仿真中最基础的电位场求解模型。它假设电解质中离子浓度均匀、无浓度极化，电极反应动力学无限快（即过电位可忽略），因此整个体系的电位分布完全由电解质欧姆定律和电极表面的电荷守恒决定。理解一次电流分布，是后续引入二次分布（含动力学过电位）和三次分布（含传质）的起点。

物理本质与适用边界

在一次电流分布中，电解质被视为具有均匀电导率 κ 的欧姆导体，电位 φ 满足拉普拉斯方程：

∇·(κ∇φ) = 0    （电解质域内）
i = -κ∇φ        （电流密度矢量）
n·i = i_n       （电极边界法向电流密度）

由于不考虑电极反应动力学，电极边界上的局部电流密度 i_n 与局部过电位无关，仅由几何形状和电导率决定。这意味着一次电流分布对电极边缘、尖角等几何突变极为敏感——在理想情况下，尖锐边缘处的电流密度理论上趋于无穷大。因此，一次电流分布最适合用于以下场景：电解质电导率较高、电极反应极快（如高浓度支持电解质中的简单氧化还原）、或需要快速评估几何均匀性的初步分析。

COMSOL 建模流程与关键设置

在 COMSOL Multiphysics 中，一次电流分布通过“电化学模块”下的“一次电流分布”接口实现。建模的核心步骤包括几何构建、电解质域设置、电极边界定义和求解器配置。

几何构建：以平行板电极电解池为例，假设两块 20 mm × 20 mm 的矩形电极平行放置，间距 10 mm，电解质充满其间。在二维轴对称或三维中建立对应几何。若研究边缘效应，需在电极边缘附近加密网格。

电解质电导率：在“电解质”节点中输入电导率 κ。对于 0.5 mol/L 硫酸钠水溶液，25°C 时 κ ≈ 5.5 S/m；对于 0.1 mol/L KCl，κ ≈ 1.3 S/m。电导率可以是常数，也可以是温度或浓度的函数，但在一次分布中通常设为常数。

电极边界：在“电极”节点中指定电极表面。对于一次电流分布，电极反应动力学被简化为“无限快”，因此只需指定总电流或平均电流密度。例如，设置阴极总电流为 0.1 A，或平均电流密度为 250 A/m²（0.1 A / 0.0004 m²）。COMSOL 会自动将总电流分配到各电极边界，分配比例由电位场决定。

边界条件：电解质外边界通常设为绝缘（n·i = 0），除非有电流从侧面流出。对称面可设为对称条件以减小计算量。

网格与求解：一次电流分布是线性问题，使用“稳态”研究即可。由于电流密度在电极边缘可能急剧变化，建议在电极边缘使用边界层网格或局部加密。求解后，可得到电解质电位 φ、电流密度矢量 i、电极表面局部电流密度 i_n 等物理量。

结果解读与几何/电导率影响分析

求解完成后，在“结果”中绘制电解质电位 φ 的等值线图和电流密度矢量的箭头图。对于平行板电极，中心区域的电流线基本平行，电位等值线近似为平面；但在电极边缘附近，电流线向边缘弯曲，局部电流密度急剧升高。这种“边缘效应”是一次电流分布的典型特征。

几何影响：电极间距越小，边缘效应越显著；电极面积越大，中心均匀区域占比越高。若将电极边缘倒圆角（如半径 1 mm），可有效降低边缘电流密度峰值。在 COMSOL 中，可通过参数化扫描圆角半径，观察最大电流密度与平均电流密度的比值变化。

电导率影响：电导率 κ 均匀时，电位分布形状与 κ 无关（因为拉普拉斯方程是线性的，κ 只影响电流密度的绝对值）。但若电解质中存在电导率梯度（如不同浓度区域），电流线会发生折射，高电导率区域电流密度更大。在一次电流分布中，通常假设 κ 为常数，因此电导率主要影响槽电压：U = I·R_Ω，其中 R_Ω 与 1/κ 成正比。

网格与求解器注意事项

一次电流分布是线性稳态问题，COMSOL 默认使用直接求解器（如 MUMPS）即可快速求解。但需注意：电极边缘的电流密度奇异性会导致网格收敛性变差。建议在电极边缘使用“边界层”网格，第一层厚度约为电极间距的 1/100 至 1/50。若使用自适应网格细化，可设置基于电流密度梯度的误差估计，让求解器自动加密边缘区域。

此外，若模型包含多个电极对或复杂几何，应检查总电流守恒：所有电极边界上的积分电流之和应为零（阳极流出等于阴极流入）。在 COMSOL 中可通过“派生值”下的“积分”功能验证。

从一次分布到二次分布的过渡

一次电流分布忽略了电极反应动力学，因此当电极反应速率有限时，实际电流分布会比一次分布更均匀。二次电流分布通过引入 Butler-Volmer 方程，将局部电流密度与过电位关联，从而在电极表面形成“动力学阻力”，抑制边缘电流集中。在实际仿真中，可先运行一次电流分布获得快速估算，再切换到二次电流分布进行更精确的分析。若电解质浓度变化显著（如低浓度或高电流密度），则需进一步使用三次电流分布，耦合物质传递方程。

掌握一次电流分布的设置与解读，不仅有助于理解电位场的基本行为，也为后续更复杂的电化学仿真奠定基础。在实际操作中，建议从简单平行板几何入手，逐步调整电导率、电极尺寸和间距，观察电流密度分布的变化规律，建立对一次电流分布的直观认识。