4.4 化学与电化学上手：扩散-反应与电流分布

4.4 化学与电化学上手：扩散-反应与电流分布

本节摘要：化学物质传递是耦合性最强的物理场家族：浓度场自己不动，全靠流动带着走、靠反应消耗着变。本节用稀物质传递接口跑一个有解析解的扩散-反应算例，再给电化学建模画出"一次、二次、三次电流分布"的层级路线图，需要电化学模块的部分明确标注。

从扩散-反应起步

稀物质传递接口（td）在核心许可里就有，方程一眼可读：左手是浓度的积累与扩散，右手是反应源项。案例设计成有解析解的形式，方便对账：一条 2 mm 长、0.5 mm 见方的通道条，物质 A 从左端面以 1 mol/m³ 的浓度注入，在体内被一级反应消耗，消耗速率 R = −k·c，扩散系数 D = 1e-9 m²/s，速率常数 k = 0.01 1/s。动手前先算特征数：反应-扩散长度 λ = √(D/k) = √(1e-9/0.01) ≈ 0.32 mm——浓度大约每零点三毫米衰减一个特征尺度，两毫米的通道足够让它衰减殆尽；无 preceded 量纲的判据是 Damköhler 数 Da = k·L²/D = 0.01 × 4e-6/1e-9 = 40，远大于一，反应主导、扩散供料，这套系统属于"吃得多、送得慢"的 типичный 传质限制工况。

模型向导：三维、稀物质传递（td）、稳态。几何一个 2 × 0.5 × 0.5 mm 的长方体。传质属性节点里把扩散系数 D_c 填 1e-9（自定义，属 3.1 的空材料套路：属性树→传质属性→扩散系数）。边界两笔：右键 td → 浓度，选左端面，填 c0 = 1 mol/m³；右端面保持默认零通量。反应：右键 td → 反应，选择域全部，R 填 −k_react*c，k_react 在全局定义里建参数 0.01[1/s]——写成参数而非裸数字，第 7 章扫描 k 时这就是现成的旋钮。

网格这里有个讲究：衰减长度只有 0.32 mm，而域长 2 mm，浓度在左端头一小段里急剧变化，均匀网格会把这段"吞"过去。给整个域加"尺寸"子节点、最大单元 0.05[mm]，沿 x 向约四十个单元，衰减区被稳稳覆盖。计算，沿中轴线做一维线图：浓度从左端的 1 单调衰减，到右端只剩千分之几。

图：扩散-反应浓度剖面与电化学电流分布层级

对账与解读

稳态扩散-反应方程在两端边界下的解析解是双曲余弦形式：c(x) = c0·cosh((L−x)/λ)/cosh(L/λ)。仿真线图与这条曲线叠在一起，几乎分不出彼此；用派生值的"点求值"在 x = 1 mm 处读数，与公式代入值之差在小数点后第三位，账干净利落。解读时把两个特征数挂嘴边：衰减长度 λ 说的是"反应吃料的势力范围"，Damköhler 数说的是"消耗与供给谁的势力大"。Da 为 40 意味着注入面附近一小层就把反应消耗扛了大半，出口浓度趋零——若工艺目标是让反应物跑满全程（比如均匀镀层），这个工况就不合格，得调 k 或流速。仿真此刻给出的是方向性判断，不是装饰数字。

电化学：一张层级路线图

电化学把电流、电位、物质传递三件事拧在一起，接口随模块分层：电化学模块与电池模块提供一次、二次、三次电流分布接口。层级选择有明确逻辑：一次只解电解液与电极里的欧姆导电，电极反应按理想化处理，适合动力学极快、浓差可忽略的体系，比如电导率极高的强电解液粗算；二次加入 Butler-Volmer 动力学，交换电流密度与过电位的关系进入方程，能算出真实的极化曲线，是电镀、电解精炼的工程主力；三次再把物质传递与浓差纳入，大电流密度下反应物来不及扩散到位、浓度梯度显著时必须用它，电池全极片模型是它的旗舰应用。选择原则与湍流模型如出一辙：从最简单的层级开始，逐级加复杂度，每一级的判据都是与实验数据的对照。操作节奏不变：向导、几何、材料（电解液电导率、迁移参数）、边界（电极电位或电流）、网格、求解、对账——第五步的对账对象从公式换成了实验极化曲线，这也是电化学建模对实验依赖更重的原因。

本节要点回顾

• 特征数先行：衰减长度 λ = √(D/k) 与 Damköhler 数 Da = k·L²/D 两个手算决定网格与解读方向。
• 对账有解析解：双曲余弦剖面叠图比对，点求值误差在小数点后第三位。
• 网格跟着 λ 走：衰减区单元尺寸压到特征长度的零头，均匀网格会吞掉关键梯度。
• 电化学三级递进：一次解欧姆、二次加动力学、三次加传质，从简到繁逐级与实验对照。
• td 接口核心许可可用，电化学深度案例需模块；流动与传质的耦合预告了下一章的多场世界。

补电化学仿真的入门路径。COMSOL 的电化学模块覆盖电池、腐蚀、电解等场景，新手入门建议从"一次电流分布"接口开始——它只考虑欧姆降，忽略所有动力学与传质效应，方程最简，适合验证几何与边界条件。掌握一次分布后升级到"二次电流分布"（加入Butler-Volmer方程描述电极动力学），最后再到"三次电流分布"（加入传质——扩散与对流）。三步升级的每一步都增加一个物理因素：第一步看几何，第二步加动力学，第三步加传质。新手常犯的错误是一步到位选最全的接口——方程多了、参数多了、不收敛了就不知道该调哪个。正确的做法是逐级增加物理复杂度，每增加一级都对比上一级的结果看变化在哪里——如果二次分布和一次分布的结果几乎一样，说明动力学不是你体系的瓶颈，可以放心用简单模型继续。这种"由简入繁"的方法论是所有多物理场仿真（不只电化学）的最佳实践。

补电化学仿真中电极动力学的入门要点。Butler-Volmer 方程是电化学的核心方程，描述电极表面的电流密度与过电位（实际电位与平衡电位之差）的指数关系。在 COMSOL 中使用时需要指定：交换电流密度 i₀（决定反应速率的快慢）、传递系数 α（决定氧化与还原的方向偏好，通常取 0.5）、以及反应物的浓度依赖。新手常犯的错误是把交换电流密度设得太大——理论上无穷大电流密度在物理上意味着"反应瞬间完成"，会导致数值不收敛。实际操作：i₀ 的典型值在 10⁻⁶ 到 10⁻³ A/cm² 之间，从文献取值或从实验的 Tafel 图反推。另外一个关键参数是平衡电位 E_eq——不同材料的平衡电位不同（锂约 -3V vs 氢 0V），多电极体系中谁的电位高谁就是阴极。掌握这两个参数的物理含义，你就理解了电化学仿真的核心——剩下的只是几何和边界条件的变化。

一个电化学仿真的入门案例

以铜的电解沉积为例走一遍电化学仿真。场景：两个铜电极浸入硫酸铜溶液，外加电压后铜从阳极溶解、在阴极沉积。设置：物理场选二次电流分布（含 Butler-Volmer 动力学）；几何是二维矩形表示电解槽截面；材料参数：交换电流密度 i₀ = 0.001 A/cm²、传递系数 α = 0.5、平衡电位 E_eq = 0.34 V（铜/硫酸铜标准电位）；电解液电导率约 5 S/m。边界条件：阳极面加恒电位 0.4 V、阴极面加 0.3 V，其余面绝缘。求解后关注两个量：阴极面的电流密度分布（均匀性决定沉积质量）和总沉积速率（法拉第定律把电流换算为沉积厚度）。结果的物理直觉校验：电流密度在阴极边缘应该略高（电场集中），但不应出现极端的局部峰值——如果峰值高出平均值数倍，说明电解槽的几何设计有问题（阳极阴极距离不够或对齐偏差）。这个案例展示了电化学仿真从"一个氧化还原反应"到"数值可预测的沉积过程"的完整入门路径。