电镀模拟
电镀的理论基础
电镀模拟概述
什么是电镀模拟？它预测什么？
简单地说，电镀模拟预测镀膜厚度的均匀性。工件表面各处应该厚度相同，但实际上由于电流分布不均，某些地方厚、某些地方薄。通过提前用数值计算来预见这一点，就是电镀模拟的作用。
为什么厚度会不均匀？
电流"喜欢"走最近的路。例如给汽车保险杠镀铬时，阳极靠近的凸起部分电流密度大，导致镀层厚；而凹陷处和孔的深处电流密度小，镀层就薄。这个"缓冲能力"问题是电镀模拟诞生的最大动力。
实际应用中什么地方用到？
几个典型例子——
- PCB通孔内镀层：直径0.1〜0.3 mm的孔要在内壁均匀沉积Cu。孔的深宽比越大越困难
- 半导体大马士革布线：宽度10 nm以下的沟槽要用Cu完全填充（超级填充）
- 汽车装饰镀层：复杂3D形状的保险杠和格栅要均匀镀铬
- 航空零部件硬铬镀：耐磨性要求高的零件需严格膜厚管理
- 镀具与治具设计：辅助阴极和遮蔽板的配置优化来提高均匀性
这些产品都离我们很近呢！那理论上怎么计算？
电镀模拟按精度分为3个阶段：一次电流分布、二次电流分布、三次电流分布。精度越高，物理描述越完整，但计算成本越大。我们逐个来看。
一次电流分布（拉普拉斯方程）
一次电流分布是什么？
一次电流分布完全忽略电极反应的阻力，只从电解液的电位分布求电流。假设电解液导电率 $\kappa$ 均匀且无化学反应影响，电位 $\phi$ 就满足拉普拉斯方程：
电流密度由电位梯度给出（欧姆定律的微分形式）：
就这么简单？
方程式确实简单，但一次分布给出最不均匀的分布，是"最坏情形"。因为把电极面当作等电位，电流纯粹按几何最短路走。凸起部分电流集中。实际的电镀因为电极反应有阻力，结果会比一次分布均匀些。所以一次分布可作为"均匀性的下界"。
二次电流分布（Butler-Volmer方程）
二次分布增加了什么？
二次分布考虑电极表面电荷转移反应的速率论。金属离子接收电子变成金属原子这个反应有一定的"阻力"。Butler-Volmer方程描述这个关系：
这么多符号，一个一个解释好吗？
| 符号 | 含义 | 典型值·说明 |
|---|---|---|
| $i$ | 电极面电流密度 | [A/m²] |
| $i_0$ | 交换电流密度 | 取决于反应系（Cu：约1〜10 A/m²，Ni：约0.1〜1 A/m²） |
| $\alpha_a, \alpha_c$ | 阳极侧·阴极侧的转移系数 | 一般 $\alpha_a + \alpha_c = 1$（单电子反应时） |
| $F$ | 法拉第常数 | 96,485 C/mol |
| $\eta$ | 过电压（活性化过电压） | $\eta = \phi_m - \phi_s - E^{\text{eq}}$（金属电位 - 溶液电位 - 平衡电位） |
| $R$ | 气体常数 | 8.314 J/(mol·K) |
| $T$ | 绝对温度 | [K] |
过电压 $\eta$ 是关键参数吧？大的话反应快？
没错。过电压 $\eta$ 是从平衡状态的偏离程度，$|\eta|$ 越大电流越强。关键是这个式子是非线性边界条件。电解液内部还是拉普拉斯方程，但电极面上不是把 $\phi$ 固定，而是用Butler-Volmer式把 $i$ 和 $\eta$（即 $\phi$）连接。这个非线性性使得原本集中到凸起的电流被"反应阻力的墙"挡住，结果电流分布就变均匀了。
反应慢的（$i_0$小的）镀液反而会均匀？
眼光敏锐！正是如此。$i_0$ 越小，电极反应越成为限速步骤，电流分布趋向均匀。青化铜镀液的均匀性优异，就是因为反应速度论效应很强。
三次电流分布（Nernst-Planck方程）
三次分布还要加什么？
三次分布纳入物质输运（质量转移）。镀着进行中，电极附近的金属离子浓度下降（浓度分极化），补给跟不上。离子的移动由Nernst-Planck方程支配：
右边3项分别是：
- $-D_i \nabla c_i$：扩散（浓度梯度驱动，斐克定律）
- $-\frac{z_i F D_i c_i}{RT} \nabla \phi$：泳动（电场对荷电离子的拖拽）
- $c_i \mathbf{v}$：对流（液体流动的运输）
其中 $D_i$ 是扩散系数，$c_i$ 是物质 $i$ 的浓度，$z_i$ 是电荷数，$\mathbf{v}$ 是流速矢量。
对流也要算，那需要流体力学的计算？
严格做的话是这样。用Navier-Stokes方程求出流场，代入Nernst-Planck方程的对流项。计算量相当大。但实务中，常用Nernst扩散层模型的近似法：只在电极附近厚度 $\delta_N$（典型值10〜500 μm）的扩散层内做计算，把这个参数给定就行。
比如PCB的通孔，孔内液体不动，浓度容易降低？
完全正确。孔的深处几乎停滞，只能靠扩散供应离子。孔的深宽比越大，浓度分极化越严重，膜就越薄。这正是三次分布发挥威力的地方。
Wagner数与均匀电着性
到底该用一次还是二次还是三次？选择标准是什么？
计算代价相差很大，所以选择合适精度很关键。判断的工具是Wagner数（瓦格纳数）Wa：
$\kappa$ 是电解液导电率，$L$ 是工件的代表尺寸，$\partial\eta/\partial i$ 是分极曲线的斜率（电极反应的微分阻抗）。物理上是"电极反应阻抗 / 电解液阻抗"的比值。
- Wa ≪ 1：电解液阻抗占主导 → 接近一次分布（不均匀）
- Wa ≫ 1：电极反应阻抗占主导 → 电流分布均匀
- Wa ≈ 1：两者影响相当 → 必须做二次分布计算
Wa数看一下就知道要用几次分布了。实际中的值是多少？
| 镀液类型 | 典型Wa值 | 均匀性倾向 |
|---|---|---|
| 硫酸铜镀（酸性） | 0.1〜1 | 低〜中 |
| 青化铜镀 | 5〜50 | 高（均匀性优秀） |
| 瓦特镍镀 | 0.5〜5 | 中 |
| 硬质铬镀 | 0.01〜0.1 | 极低 |
| 锡镀 | 1〜10 | 中〜高 |
硬质铬镀Wa才0.01？会非常不均匀呢。
没错。所以铬镀中要用辅助阴极（罗宾）和遮蔽板来调节电流分布。这些治具配置的最优化，正是模拟的强处。
膜厚成长模型
电流分布求出来后，膜厚怎么算？
用法拉第电解定律，从局部电流密度算膜的成长速度：
- $\delta$：膜厚 [m]
- $M$：沉积金属的摩尔质量 [kg/mol]（Cu: 0.0636, Ni: 0.0587, Cr: 0.0520）
- $i_n$：电极面法向电流密度 [A/m²]
- $n$：电子转移数（Cu²⁺→Cu: $n=2$, Cr⁶⁺→Cr: $n=6$）
- $\rho$：沉积金属的密度 [kg/m³]
- $\eta_{\text{CE}}$：电流效率（扣除副反应。铬镀极低，仅10〜25%）
例如硫酸铜镀，$i = 300$ A/m²、电流效率100%时，成长速度约3.7 μm/min。
手机处理器里的镀层模拟
你的智能手机处理器里有数十亿条宽度10 nm以下的铜布线，每一条都是用电镀（大马士革工艺）做出来的。最难的是在细深沟槽的底部从下往上长出铜——"超级填充"。这需要PEG（聚乙二醇）等抑制剂、SPS（双(3-磺基丙基)二硫化物）等加速剂和JGB（甲烷绿B）等均一化剂三种有机添加剂互相竞争，在电极表面吸附、扩散、消耗。用物质转移模型计算这种复杂竞争，找到从沟槽底部填充的工艺条件。一块芯片需要的镀层模拟计算量是一辆汽车保险杠的数千倍，据说。
电镀的数值计算方法
FEM和BEM的应用选择
镀层模拟用什么数值方法？有限元法吗？
主要两种方法。
- 有限元法（FEM）：整个电解液区域分割成单元。对二次三次分布和非线性材料（导电率浓度依赖等）有优势。COMSOL和Ansys采用
- 边界元法（BEM）：只需对电极和绝缘壁表面网格化。一次分布（拉普拉斯方程）用BEM体积网格省了，网格生成轻松。3D复杂镀具设计的优势明显。ElSyCA PlatingMaster代表
BEM听着简便，缺点呢？
BEM得到密矩阵（full matrix），内存随单元数平方增长。数万单元超过就计算变重。还有，电解液导电率空间变化（浓度依赖等）时BEM不好用。实务中，一次二次分布用BEM，三次分布需要用FEM，这样分工。
移动边界的处理
镀着进行中形状变了。那怎么计算？
镀沉堆积导致电极形状时刻变化——这是移动边界问题。主要方法3种：
- ALE法（任意拉格朗日欧拉法）：网格随膜成长变形。COMSOL的Deformed Geometry模块用。网格歪斜过大时需重新网格
- 水平集方法：隐式追踪界面。对拓扑变化（孔填满等）有利
- 准定态近似：膜厚变化缓慢的情况，在该时刻形状下解电流分布，更新膜厚，进行下一步。很多实际情形精度足够
准定态近似用得到的时候怎么判断？
膜厚变化量相对工件代表尺寸充分小就行。比如膜厚30 μm、工件凹凸数mm以上，形状变化小于0.1%，准定态成立。但半导体沟槽填充（孔宽100 nm、膜厚50 nm）形状剧变，非用ALE法或水平集法不可。
非线性求解器和收敛策略
Butler-Volmer式非线性，求解不稳定吗？
其实就是。Butler-Volmer式的指数函数对过电压变化剧烈，初值稍差就牛顿法发散。实务的诀窍总结一下：
- 初值用一次分布的解：拉普拉斯方程的解做热启动
- 带阻尼的牛顿法：迭代更新量乘0.3〜0.7的缓和系数防发散
- 电流分阶段增加：最终电流不一次用上，0.1倍→0.5倍→1.0倍逐步加（伪时间步进）
- 从泰费尔近似切换：先用泰费尔近似（对数函数，非线性温和）找初值，再切到Butler-Volmer
一步到位目标条件解的、段階的。
对。电化学的非线性方程需要"热身"。COMSOL有Auxiliary Sweep分阶段改参数。Ansys Fluent的Pseudo Transient让计算稳定。
电镀的实务应用
分析工作流程
实际做镀层模拟，怎么进行？
标准流程是这样：
- CAD模型准备：镀槽、阳极、阴极（工件）、治具（辅助电极、遮蔽板）的3D几何。电解液实际充满的区域提取
- 电解液域定义：导电率 $\kappa$、反应物的扩散系数 $D_i$、体积浓度 $c_\infty$ 设定
- 电极反应参数设定：交换电流密度 $i_0$、转移系数 $\alpha$、平衡电位 $E^{\text{eq}}$、电子数 $n$
- 网格划分：电极附近细，内部粗
- 求解器设置和运行：应用上述非线性收敛策略
- 后处理：电流密度分布、电位等值线、膜厚分布可视化，评估均匀性
网格划分的要点
镀层模拟特有的网格注意事项？
结构分析和很不同。重点排序：
- 电极边缘细分：阴极的棱、孔口等电流汇聚处。从棱向外几mm范围内单元尺寸0.1〜0.5 mm
- 扩散层解像：三次分布要在电极附近数十μm的扩散层内3〜5层单元。用棱柱层网格（prism mesh）有效
- 阳阴间隙方向：两电极距离方向至少10单元。这里的解像度直接关系电解液的电位降
- 内部粗化OK：电极远处电位梯度小，网格10倍粗也结果不变的多
结构分析違応力集中部電極。
物理同、"拉普拉斯方程的解特异的场所"网格集中。90度角的場合理論上 $r^{-1/3}$（$r$ 距離）的特異性。结构分析的応力特異点本質同数学。
边界条件设置
边界条件设定的一般的？
| 境界 | 一次分布 | 二次分布 | 三次分布 |
|---|---|---|---|
| 阴极面 | $\phi = \phi_c$（等电位） | Butler-Volmer式 | Butler-Volmer式 + 浓度依赖 |
| 阳极面 | $\phi = \phi_a$（等电位） | Butler-Volmer式或全电流指定 | 同左 + 溶解反应速率 |
| 绝缘壁 | $\partial\phi/\partial n = 0$ | 同左 | 同左 + $\partial c/\partial n = 0$ |
| 液面（开放面） | $\partial\phi/\partial n = 0$ | 同左 | $c = c_\infty$（体液浓度） |
实务上的要点：阳极分溶解型（Cu、Ni）和不溶性（Pb-Sn、Ti-Pt、DSA）两类，不溶性阳极还要含放氧反应的速率论。
验证和有效性确认
模拟结果的正确认方法？
镀层模拟的V&V（验证与有效性）的定番方法3段階：
- 分析解的比较（Verification）：同心円筒电极（的分析解）（Hull Cell）的理論解一次分布验证。全電流的積分值入力電流一致必確认
- 実験的比较（Validation）：267 mL標準化試験槽、阳极阴极中的距离直线的変。実测的膜厚分布仿真結果比较
- 実機试镀的比较：最后实际镀液试作品镀、蛍光X線（XRF）膜厚实测、比较
実験仿真的検证使！
。「一次電流分布的分析解知簡単形状」、仿真代码的基本的正确认的最適。ASTM B571规定标準试験法。
电镀的软件对比
支持工具与功能对比
镀层模拟用什么软件？
| 软件 | 解法 | 电流分布 | 移动边界 | 添加剂模型 | 优势 |
|---|---|---|---|---|---|
| COMSOL Multiphysics （Electrodeposition模块） | FEM | 一次〜三次 | ALE / 水平集 | 自定义PDE | 多物理耦合的灵活性。电化学领域实绩最多 |
| Ansys Fluent （电化学模型） | FVM | 二次、三次 | 动态网格 | UDF | 流动耦合优势。大规模并行计算 |
| STAR-CCM+ （电化学） | FVM | 二次、三次 | 网格变形 | Java API | 多面体网格自动生成 |
| ElSyCA PlatingMaster | BEM | 一次、二次 | 表面网格更新 | 有限 | 3D复杂形状镀具设计专用。高速 |
| ELSYCA CuPlate | BEM+FEM | 一次〜三次 | 支持 | 支持 | PCB镀专用。孔填充强 |
结局使？
用途次第。
- 研究论文：COMSOL压倒的。自由自定義方程式書、引用多
- 汽车装饰镀的治具设计：PlatingMaster的BEM固有的高速定評
- 流的影响大（喷流镀镀）：Ansys FluentSTAR-CCM+的CFD连成强
- 半导体/PCB孔埋：CuPlateCOMSOL的水平集法
会社Ansys持、Fluent的电化学模型始的初期投资现実的。
电镀的前沿研究
添加剂建模
実際的镀液添加剂入？ 仿真？
。特半导体Cu镀添加剤的必須。3種的添加剤竞争的電極表面吸着析出速度制御：
- 抑制剤（Suppressor）：PEG。電極吸着析出遅。凸部多吸着
- 加速剤（Accelerator）：SPS。表面吸着析出速。凹部的底蓄积
- 均一化剂（Leveler）：JGB。对流依存的吸着、凸部的加速剤効果打消
表面被覆率 $\theta$ 追跡表面扩散吸着模型、中的輸送方程式连立。COMSOL的Surface Reactions实装。
脉冲镀模拟
脉冲镀的。仿真難？
脉冲镀ON/OFF（/）周期的繰返。DC镀均匀性结晶性向上。仿真上：
- ON期間：通常的電流分布计算（秒〜秒）
- OFF期間：電流拡散層内的濃度回復（扩散方程式解）
- ：逆電流表面的凸部優先的溶解
時間秒秒大異、適応的重要。ON/OFF切替直後 $\Delta t$ 極小（μs）、定常近大。
数字孪生与工艺优化
最近聞、镀？
先進工場镀的实用化進。具体的：
- 濃度：液分析的数据仿真模型入力、膜厚分布予測
- 治具形状的自動优化：遮蔽板的形状位置参数、膜厚的標準偏差最小化优化
- 模型：高精度FEM大量走結果機械学習模型構築、推論時間秒単位短縮
- 过程可视化：電流密度温度液流速的組合、品質OK/NG領域明示
特自動車业、新车种的复雑外装部品的镀具试作设计流程広。
电镀的问题排查
常见错误与对策
镀层仿真教！
| 症状 | 原因 | 対策 |
|---|---|---|
| 牛顿法发散 | Butler-Volmer初值不良 | 一次分布的解初值使用。阻尼系数0.3〜0.5开始 |
| 付近電流密度異常大 | 一次电流分布的特异性（幾何学的效果） | 细分化（要素0.1mm以下）。丸付（实形状忠实） |
| 膜厚合計則合 | 電流効率未设定、网格的電流积分误差 | 電流効率 $\eta_{\text{CE}}$ 正设定。境界積分的精度网格収束确認 |
| 三次分布負的濃度出 | 拡散層的要素粗 | 境界層网格追加（最低3層）。上流化適用 |
| ALE法网格潰 | 膜厚成長対网格頻度不足 | 网格间隔短。変形量要素的30%超网格 |
| 计算時間異常長 | 脉冲镀時間小 | ON/OFF切替後的微小。定常部大自動切替 |
的特異性、模型的角丸解消？
実物的完全、微小丸（ R=0.1〜0.5 mm）付的物理的正。仿真上特异性緩和网格収束性大幅改善。结构分析応力集中部付的同考方。
的一次分布自分解。COMSOL的例題。手動的一番的近道。
错误