《Materials Today Communications》：Modeling pitting corrosion in stainless steel: A comprehensive review
编辑推荐：
这篇综述系统梳理了不锈钢中点蚀腐蚀的建模进展。文章沿时间轴追溯了从早期随机模型到近期多物理场及无网格框架的演化历程，重点剖析了确定性模型（如相场法PFM）、随机模型及混合模型的核心公式、物理假设和预测能力，为理解局部腐蚀模拟的技术路径和发展方向提供了全面参考。
点蚀腐蚀：不锈钢的隐匿威胁
点蚀腐蚀是一种高度局域化且极具侵蚀性的材料降解过程，对不锈钢在役期间的长期可靠性构成严重威胁。它始于局部钝化膜的破损，通常在富含氯离子等侵蚀性阴离子的环境中被触发。一旦形成闭塞环境，内部会酸化并富集氯离子，导致自催化反应，使点蚀迅速向材料内部扩展，往往引发突然的灾难性失效。全球每年因腐蚀造成的经济损失高达约2.5万亿美元，而点蚀由于难以检测和预测，是其最具危险性的形式之一。准确预测点蚀行为对于航空航天、海洋工程和能源等关键领域至关重要。
点蚀生长动力学与反应机制
点蚀的扩展由发生在蚀坑-电解质界面上的电化学反应主导。核心是金属在闭塞坑内的阳极溶解反应（如Fe → Fe2++ 2e-），释放的电子被外部金属表面的阴极反应（如氧还原）消耗，形成一个局部的电化学电池。蚀坑的生长速率取决于阳极电流密度，在早期活化控制阶段，它遵循Butler-Volmer方程，与过电位呈指数关系。随着金属离子在坑内积累并发生水解，环境变得更为苛刻，体系会过渡到扩散控制阶段，此时金属溶解速率受限于离子通过坑内电解质的传输，达到一个极限电流密度。这一转变对于点蚀的稳定和特征性的碟形或叶瓣状形貌演化至关重要。
阴极反应虽然常被简化处理，但对于维持电流连续性和影响局部电位分布同样重要，环境因素（如pH、温度、氧浓度）会显著影响其效率。临界点蚀温度（CPT）是评估不锈钢点蚀稳定性的一个关键参数，它标志着从亚稳态点蚀向稳定生长转变的温度阈值。
离子传输与电荷平衡
准确模拟点蚀生长需要耦合描述离子传输（扩散、电迁移、对流）的Nernst-Planck方程、质量守恒定律以及电荷守恒条件（如电中性假设或Poisson方程）。在多数工程模型中，宏观电中性假设简化了计算，但在界面附近电荷分离显著的狭窄区域，则需要求解完整的Poisson-Nernst-Planck（PNP）方程组。这些控制方程构成了有限元、相场法等数值方法模拟点蚀的基础。
环境与力学耦合效应
在流动环境中，电解质运动会影响局部传质，从而改变点蚀行为。将腐蚀模型与Navier-Stokes方程耦合，可以模拟流速、剪切应力以及自然对流（通过Boussinesq近似）对点蚀的影响。机械载荷同样重要，它通过改变金属的化学势和界面反应驱动力来影响腐蚀。例如，拉应力会增加局部阳极溶解速率，而塑性应变会通过产生缺陷促进钝化膜破裂，加速应力腐蚀开裂（SCC）的萌生。氢脆（HE）机制也可能在某些环境中与SCC共同作用，加剧材料的脆性断裂倾向。
微结构与界面的影响
材料的微观异质性（如晶粒取向、晶界、非金属夹杂物、δ-铁素体条纹）是点蚀形核和扩展的关键控制因素。这些特征会改变局部成分和钝化膜的力学状态，形成薄弱点。例如，含钇不锈钢中，规则形状的Y2O3夹杂物相对惰性，而不规则夹杂物则容易在其边缘诱发快速点蚀。一旦点蚀形成，周围的微观结构（如晶界、第二相粒子）将继续影响其扩展路径和速率，导致局部溶解速率的空间差异。此外，金属-电解质界面的静电条件（如双电层结构）和电位分布也会影响侵蚀性离子（如Cl-）的局部分布和钝化膜的稳定性。
点蚀腐蚀的建模方法概览
模拟方法主要分为确定性、随机性和混合型三类。确定性模型基于物理定律（如扩散方程、电动力学），能够精确模拟点蚀生长动力学、形貌演化和传输过程，主要包括有限元法（FEM）、任意拉格朗日-欧拉法（ALE）、相场法（PFM）、水平集法（LSM）和近场动力学（PD）等。这些方法通常不处理形核的随机性。随机模型则专注于描述点蚀形核、亚稳态点蚀行为及再钝化等具有统计波动性的过程，常基于泊松过程、马尔可夫链或蒙特卡洛模拟。混合模型结合了上述两者的优势，例如用随机过程描述形核，再用确定性模型模拟后续的生长，为复杂点蚀场景提供了更全面的描述。
确定性模型详解
确定性模型通过求解控制方程来模拟点蚀扩展，提供了高分辨率的结果，但计算成本也相应较高。
- •
相场法（PFM）是一种强大的介观尺度模拟工具，它通过引入一个连续的序参量（相场变量）来区分金属和电解质相，界面由该变量的平滑过渡区域表示。通过求解相场方程和耦合的物理场（如浓度场、电势场）方程，PFM可以自然而优雅地处理复杂的界面形貌演化、拓扑变化（如合并、分裂）以及多物理场耦合问题，非常适合模拟点蚀的不规则生长、裂纹萌生以及微结构敏感腐蚀。
- •
水平集法（LSM）和任意拉格朗日-欧拉法（ALE）是显式界面追踪方法。LSM用水平集函数的零等值面表示界面，通过求解输运方程来更新界面位置，能很好地处理大变形和拓扑变化。ALE方法则通过动态调整网格节点来贴合移动界面，在界面附近保持高精度，但复杂形貌下可能需要频繁的网格重构。
- •
近场动力学（PD）是一种非局部连续介质理论，用积分方程代替偏微分方程，天然适用于模拟不连续问题（如腐蚀导致的材料缺损、裂纹）。在点蚀建模中，PD通过键的断裂来模拟材料的溶解，能够捕捉点蚀的萌生、扩展以及多蚀坑相互作用。
- •
元胞自动机（CA）是一种基于离散格子和简单演化规则的模型。在点蚀模拟中，每个元胞根据其邻居状态和局部规则（代表电化学反应和传质）更新其状态（如金属、电解质、钝化膜）。CA计算高效，易于并行，适合模拟大规模的微结构敏感腐蚀和形貌演化，但物理机制的描述可能不如连续介质方法精细。
随机与混合模型
随机模型将点蚀形核视为一个随机事件。极值统计模型（如Gumbel分布）常用于分析最大点蚀深度，评估最坏情况下的失效风险。泊松过程和马尔可夫链则用于模拟点蚀的形核率、亚稳态点蚀的存活与死亡过程。混合模型整合了随机形核与确定性生长，例如，用泊松过程在材料表面随机生成点蚀核，然后用基于扩散-反应方程的确定性模型模拟每个核的生长，从而更真实地再现实际腐蚀表面的统计特征。
总结与展望
这篇综述通过时间线和主题区块，系统回顾了不锈钢点蚀腐蚀建模的发展历程，从早期的静态传输模型到如今能够耦合电化学、力学、流体及微结构的多物理场高级框架。当前，基于有限元等方法的连续介质扩散-反应模型在局部腐蚀模拟中仍占主导地位，而相场法、近场动力学、元胞自动机和随机框架则被更选择性地应用于特定问题，如点蚀-裂纹转变、盖层形貌演化以及微结构敏感的生长模拟。未来的研究方向可能包括开发更高保真度的跨尺度模型、整合机器学习和高通量计算以加速材料筛选、以及建立更完善的模型验证与实验标定协议，以进一步提升对不锈钢点蚀腐蚀这一复杂现象的预测能力和理解深度。
