2. 电化学基础与COMSOL建模准备

在COMSOL中建立电化学模型之前，必须先厘清几个贯穿始终的基础概念：电极反应如何描述、电流与物质变化如何换算、电位与浓度如何耦合、以及为什么实际电极电位会偏离平衡值。这些概念决定了后续选择哪个物理场接口、需要输入哪些参数、边界条件如何设置。本节先回顾核心理论，再给出在COMSOL中搭建电化学模型的基本步骤与参数设置方法。

一、电极反应与法拉第定律

电化学的核心是发生在电极/电解质界面的氧化还原反应。以常见的铜沉积为例，阴极反应可写为 Cu²⁺ + 2e⁻ → Cu，阳极则为 Cu → Cu²⁺ + 2e⁻。反应式中的电子数 n 是后续所有计算的关键参数。法拉第定律把通过的电量与物质变化量联系起来：

m = (M · I · t) / (n · F)
其中：m 为沉积/溶解质量(g)，M 为摩尔质量(g/mol)，
I 为电流(A)，t 为时间(s)，n 为电子转移数，
F = 96485 C/mol 为法拉第常数。

在COMSOL中，法拉第定律并不需要手动输入，而是由“电分析”或“二次电流分布”接口自动处理。用户真正需要设置的是反应式、电子数 n、以及交换电流密度等动力学参数。理解这一点，才能明白为什么模型里要填“参与电子数”而不是直接填质量。

二、能斯特方程与平衡电位

平衡状态下，电极电位由能斯特方程决定：

E_eq = E⁰ + (RT / nF) · ln(a_Ox / a_Red)

其中 E⁰ 为标准电极电位，R 为气体常数 8.314 J/(mol·K)，T 为温度(K)，a 为活度（稀溶液可近似用浓度代替）。以 Cu²⁺/Cu 体系为例，E⁰ = +0.34 V（相对标准氢电极），若 Cu²⁺ 浓度为 0.1 mol/L、温度 298 K，则 E_eq ≈ 0.34 + (0.0257/2)·ln(0.1) ≈ 0.31 V。这个数值就是COMSOL中“平衡电位”边界条件的输入依据。若体系涉及多个反应，需要为每个反应分别设定平衡电位。

三、过电位与极化曲线

实际电极反应发生时，电极电位 E 会偏离平衡电位 E_eq，这个差值就是过电位：η = E − E_eq。过电位是驱动反应进行的“额外代价”，它由三部分组成：活化过电位（电荷转移阻力）、浓差过电位（传质限制）和欧姆过电位（溶液电阻）。在COMSOL中，活化过电位通常用Butler-Volmer方程描述：

i_loc = i₀ · [exp(α_a·F·η / RT) − exp(−α_c·F·η / RT)]

其中 i₀ 为交换电流密度，α_a、α_c 为阳极/阴极传递系数。当 η 较大时，方程退化为Tafel形式，这也是实验中常通过Tafel斜率提取动力学参数的依据。极化曲线（E 对 log|i| 作图）能直观反映反应从活化控制到传质控制的转变，是验证模型是否合理的重要参照。

四、在COMSOL中建立电化学模型的基本步骤

打开COMSOL后，第一步是选择物理场接口。对于电极动力学主导的体系，选择“电化学”>“电分析”或“二次电流分布”；若需考虑浓度极化，则选“三次电流分布”并耦合“稀物质传递”。接口选定后，在“模型开发器”中依次完成以下操作。

几何建模：根据实际电池或电解池尺寸绘制二维或三维几何。例如一个简单的平行板电解池，可画两个矩形电极（宽 20 mm、高 30 mm）和中间的电解质域（间距 10 mm）。几何单位统一用 mm，避免后续参数换算出错。

参数设置：在“全局定义”>“参数”中集中输入所有物理量，便于统一修改。典型参数包括：温度 T = 298.15 K，法拉第常数 F = 96485 C/mol，气体常数 R = 8.314 J/(mol·K)，Cu²⁺ 初始浓度 c0 = 100 mol/m³，交换电流密度 i0 = 1 A/m²，阳极传递系数 α_a = 0.5，阴极传递系数 α_c = 0.5，电子数 n = 2，标准电位 E⁰ = 0.34 V。把这些量写成带单位的参数，COMSOL会自动处理量纲。

电极反应设置：在“电分析”接口下，右键“电极”添加“电极反应”。选择反应类型为“阴极”或“阳极”，输入化学计量系数和电子数。对于 Cu²⁺ + 2e⁻ → Cu，在阴极反应中设置 Cu²⁺ 的化学计量系数为 1，电子数为 2。若使用Butler-Volmer动力学，需指定 i0、α_a、α_c 和平衡电位 E_eq。E_eq 可以直接填能斯特方程计算值，也可以勾选“由能斯特方程计算”，让COMSOL根据浓度自动更新。

边界条件：电解质域的外边界通常设为“绝缘”，电极表面设为“电极反应”边界。若模拟恒电流过程，在阳极边界设置“总电流”或“电流密度”；若模拟恒电位过程，则设置“电位”并指定相对于参比电极的电压值。注意COMSOL中的电位是相对于“地”节点的，需要明确参比电极的电位偏移。

网格与求解：电极附近的浓度和电位梯度较大，需要在电极表面设置边界层网格。对于二维模型，电极表面第一层网格厚度可设为 0.01 mm，增长率 1.2，层数 8。求解器选择“稳态”用于极化曲线，选择“瞬态”用于恒电流阶跃或循环伏安。瞬态求解时，时间步长建议先粗后细，例如 0 到 1 s 步长 0.01 s，1 到 100 s 步长 0.5 s。

五、参数设置中的常见陷阱

初学者最容易出错的地方有三处。第一，单位不统一：浓度用 mol/L 而COMSOL默认 mol/m³，1 mol/L = 1000 mol/m³，忘记换算会导致结果偏差三个数量级。第二，平衡电位符号：能斯特方程中的对数项在阴极反应时通常取负号，因为还原态浓度在分母。第三，交换电流密度 i0 的量级：i0 从 10⁻⁶ 到 10² A/m² 不等，填错量级会让极化曲线完全偏离预期。建议在“参数”中把 i0 写成 10^logi0 的形式，方便扫描不同动力学快慢的影响。

完成上述设置后，模型即可求解。后处理中重点查看三个量：电极表面的局部电流密度分布（判断边缘效应）、电解质中的电位等值线（判断欧姆降）、以及浓度场随时间的变化（判断传质限制）。将模拟得到的极化曲线与实验Tafel斜率对比，若斜率偏差在 10% 以内，通常认为动力学参数设置合理。这些判断依据将贯穿后续所有电化学仿真章节。