Chinese Journal of Ship Research
金属表面点蚀的稳态过程数值研究
崔静，杨帆，杨霆浩，杨广峰 300300中国民航大学 机场学院，天津
崔静，杨帆，杨霆浩，杨广峰
摘 要：［目的］腐蚀会严重影响舰船设备的使用寿命，了解腐蚀机理及对腐蚀进行预测能有效地减缓腐蚀对 舰船设备的危害。［方法］构建了格子Boltzmann腐蚀模型，该模型可描述包含多相多组分流动与扩散、电化学反 应和金属点蚀的稳态全过程。应用此模型，研究浸没于液体腐蚀环境中的金属表面单坑点蚀稳态过程；分析腐 蚀化学反应速率、腐蚀溶液扩散系数、腐蚀产物扩散系数对腐蚀程度的影响。通过数值模拟，获得金属表面点 蚀坑的形貌变化特征。［结果］结果显示，对于金属表面单坑点蚀的稳态过程，由于钝化膜与金属基体构成了“大 阴极，小阳极”的电化学腐蚀体系，使腐蚀稳态点蚀将持续向金属材料基体的纵深发展，并且初生蚀孔除自身发 生点蚀外，还将在蚀孔底部产生次生蚀孔。在不同影响因素下，腐蚀程度随腐蚀反应速率和反应物组分扩散 系数的增大而增大，而随腐蚀产物扩散系数的增大而减小。［结论］通过该腐蚀模型，可模拟出与真实金属相 近的腐蚀形貌变化特征。 关键词：格子Boltzmann方法；金属点蚀；电化学反应；钝化概率；结点体积法；点蚀形貌 中图分类号：U661.5 文献标志码：A DOI：10.19693/j.issn.1673-3185. 01189
Numerical study on steady pitting process of metal surface
Cui Jing，Yang Fan，Yang Tinghao，Yang Guangfeng Airport College，Civil Aviation University of China，Tianjin 300300 Abstract：［Objectives］The service life of the shipboard equipment will be seriously affected in the harsh corrosive environment of seawater. Therefore，it can effectively reduce the equipment damages due to corrosion by understanding corrosion mechanism and predicting corrosion.［Methods］Lattice Boltzmann （LBM） corrosion model can be used to describe multiphase multicomponent flow and transmission， electrochemical reaction and metal steady pitting process. Based on this model，a numerical study is conducted to simulate a single corrosion pit on a metal surface immersed in the liquid corrosive environment；the effects of corrosion reaction rate，corrosion solution diffusion coefficient and corrosion product diffusion coefficient on the degree of corrosion are analyzed.［Results］Morphological changes of the pit on the metal surface can be obtained by the numerical simulation. The numerical results show that， the pitting corrosion will be deeper because of the electrochemical corrosion system of large cathode and small anode formed by passivation film and metal matrix during the steady corrosion process of the single pit on the metal surface；the primary pitting hole will produce secondary pitting hole at the bottom. A conclusion obtained by changing the factors in the model is that，the degree of corrosion increases as the corrosion reaction rate increases， increases as the diffusion coefficient of the reactant components increases，and decreases as the corrosion product diffusion coefficient increases. ［Conclusions］ The similar corrosion morphology of the real metal can be simulated through this corrosion model. Key words： Lattice Boltzmann method；metal pitting；electrochemical reaction；passivation probability； volume of pixel；morphology of pitting corrosion
0引言
现代舰船装备结构的系统化和复杂化程度越来越高，在舰船上使用的新材料也越来越广泛，因此舰船的腐蚀防护工作面临着巨大挑战。金属腐蚀会导致舰船装备的损坏，造成经济损失以及腐蚀污染。预测腐蚀的损伤过程，对于保证舰船装备的安全和正常运行起着重要的作用。舰船所处的海洋环境复杂，运用实际测量和经验估算等方法对舰船材料表面腐蚀进行检测时，评估标准不统一，不能实时观测腐蚀发展过程，而采用计算机技术进行数值模拟，在一定程度上可弥补传统方法的不足。腐蚀过程是一个包含了多组分之间的质量传输、化学反应或电化学反应、界面动力学等多重Boltzmann因素影响的复杂过程，而格子 方法（Lattice Boltzmann Method，LBM）作为一种介观方法（Mesoscopic technique），近年来发展迅速，成功地应用于多种微观相互作用的流体系统问题中。LBM既能兼顾实际的计算成本，又能保证模型或LBM方法的合理性。国内外运用 进行腐蚀模拟的研究较少，更多学者采用的是元胞自动机模型（Cellular Automata Model，CAM）进行腐蚀模拟，采CAM用 进行腐蚀模拟对于多组分之间的质量传输和腐蚀化学反应只进行了简单的元胞转化，并He LBM等［1］首没有理论支撑。 次将化学反应与结合，解决了流体—壁面表面反应与扩散问题的Kang等［2-3］对溶解沉淀问题中单相流体和耦合。 LBM多相流体的传输及流动进行了研究，将 与结法（Volumn of Pixel，VOP）结合来模拟物质点体积的传输过程。Chen等［4］为了解决多相流动中的相LBM VOP变和溶解沉淀问题，将 和 方法在孔隙尺度下相结合，解决了多相流中的相分离、化学反应、质量传输和溶解沉淀等问题。张婷等［5］使用MRT-LBM来模拟腐蚀溶液的流动，并分析了不Liu D3Q19等［6同参数对腐蚀的影响。 ］使用 模型Zhang等［7］在腐蚀模型模拟了岩石的腐蚀情况。中加入了电场并与流场成功耦合，更好地模拟了Boltzmann（LB）腐电化学腐蚀。但是，现有的格子蚀模型中，并未对具体的金属材料进行模拟，更多是使用化学腐蚀的方法进行研究，仅考虑了阳极固体腐蚀，很少涉及阴极反应的发生。LB本文将针对现有 腐蚀模型的不足，采用阴阳极共同作用的方式研究电化学腐蚀，将使用改进后的模型对金属表面点蚀现象的稳态过程进行模拟，得到点蚀坑的形貌变化特征以及各组分浓 度的变化情况。通过数值模拟，分析腐蚀化学反应速率、腐蚀溶液扩散系数、腐蚀产物扩散系数对腐蚀程度的影响。
1 格子Boltzmann腐蚀模型
腐蚀过程包含多相多组分的流动与扩散、电化学反应、金属的溶解以及腐蚀产物沉淀。根据以上问题建立多个子模型，并将子模型耦合形成LB腐蚀模型。其中，子模型包括电化学反应模型、多组分模型、质量传输模型、固体边界迁移模型和结点钝化模型。
1.1 电化学反应模型
模拟浸没于液体腐蚀环境的金属（以不锈钢1为例）表面的点蚀现象，其腐蚀过程如图 所示。
为了便于数值建模，本文将复杂的电化学反应过程进行如下简化。阳极： Fe-2e- →Fe2+ （1）
阴极：
2H2O+O2 +4e- →4OH - （2）沉淀反应： Fe 2 + + 2OH - → Fe（OH） 2↓ （3）沉淀氧化： 4Fe(OH )2 + 2H2O + O2 → 4Fe（OH） 3↓ （4）亚铁离子氧化： Fe2+ -e- →Fe3+ （5）铁离子水解： Fe3+ +3H2O→Fe （OH） + 6H （6） + 3自催化反应： Fe+2H + →Fe2+ +H2↑ （7）假设腐蚀环境已达到蚀孔可发生稳态点蚀时（腐蚀深度达到一定值，腐蚀产物溶度达到其饱和60%~80%［8］），如图 1（a浓度的 ）所示。蚀孔内金属基体处于活化溶解状态，作为阳极，蚀孔外钝化膜作为阴极，孔内外形成了“大阴极，小阳极”的电化学腐蚀体系，阴阳极发生式（1）和式（2）的反应， 1（b）所示；由于孔外发生式（2）反应，孔口处如图OH-逐渐增多，孔内的金属离子向孔外扩散，发生
式（3）和式（4）的反应，在孔口附近产生沉淀；随着孔内金属不断腐蚀溶解，金属离子逐渐增多，继而5 H+发生式（ ）的水解，孔内 增多，呈酸化环境，进一步加快式（6）的腐蚀，如图1（c）所示。
1.2 多组分模型
Shan-Chen模型［2-4］来处理多组分介质间使用的流动，其演化方程为fα ( x + ceα Dt t +D t )- fα ( x t )= - 1 [ )] 8 fα ( x t )- fαeq ( x t （ ） τv式中： fα ( x t ) 为格子 x 处 t 时刻的密度分布函数；fαeq ( x t ) 为平衡分布函数；c =D x/Dt ，为格子声速，其中 Dx 为空间步长，Dt 为时间步长；eα 为离散速度； τv 为流场中的松弛时间。流场采用D2Q9 ［1-3］模型，其平衡分布函数的计算公式为eq 3 fα = ωα ρ 1+ ( eα·u )+9 ( eα·u)2 - 3 u2 （9） c2 2c4 2c2式中： ωα 为权系数， ω = 4/9 ，ω = 1/9 ， 0 1 -4 ω = 1/36 ；宏观量密度 ρ 和速度 u 可由密度分5 - 8布函数求得： （10） ρ =å fα （11） ρu =å ( fαeα) Shan-Chen在 模型中，为了降低误差和保证计算结果的稳定性，使用了分子作用力对速度进行修正，分子间的作用力计算如下： N (| )ψ F =- G ψ ( x )å w eα |2 ( x + eα )eα 12 （ ） σσ σσ σ σ α =1式中：F 为液体分子间的作用力；G 为作用强σσ σσ (| )度；ψσ为有效密度函数［9-12］ ；w e |2 为权系数，其α中 | e1 |= 1 ，| e5 |= 2 ，w ( 1)= 1/3 ，w ( 2)= 1/12 。- 4 - 8 (| ) N F =- w ψ ( x )å w eα |2 s ( x + eα )eα （13） σw σw σ α =1式中：F 为液体分子与固相分子之间的作用力； σw w 为液体与固体壁面的作用强度，可通过接触σw 度［13-15］，与角的变化进行调整；s为固体壁面密 通过分布函数算出的密度不是同一概念。将上述计算得到的分子间作用力代入下式进行速度修正：
τv Fσ （14） u = u′ + ρ 11式中：F 为该结点所受到的合力；u′为式（ ）所σ求的速度；式（14）计算的 需代入式（9）进行平衡u分布函数的计算。
1.3 质量传输模型
采用质量传输模型来解决溶液组分之间的相 互传输问题，本文涉及到的液体组分为腐蚀溶液LB组分和腐蚀产物组分。 质量传输模型的演化方程为gk ( x + ceα Dt t +D t )- gk ( x t )= α α - 1 [gk )] （15） ( x t )- gk ( x t eq α α τkg式中： gk ( x t ) 为第 k 种组分的浓度分布函数α （k=1，2）；g eq为平衡浓度分布函数； τkg为浓度k α场内第 k 组分的松弛时间。在浓度场中使用的是D2Q5 ［16-18］。模型平衡浓度分布函数计算如下： ae eq 1 16 gk = Ck Jk + eα·uö （ ） è 2 ø α α 由式（11）式中：Ck 为第 k 组分的宏观量浓度；u计算得到；系数 J为kα Jk0 α =0 （17） Jk = α (1 -J α = 1 2 34 0)/4 k 0~1其中，Jk 的取值范围在 之间。0
18 Ck =å gkα （ ）浓度场中的扩散系数 Dk 为
( )( （19） D = 1 1 - Jk τk - 0.5) k 0 g 2 0.9。式中，Jk 取值为0
1.4 固体边界迁移模型
由于发生腐蚀溶解反应，被腐蚀的金属基体表面结点的体积会随之减小，并且当新生成的腐蚀产物浓度达到其饱和浓度时，在固体表面沉淀结点的体积会随之增大，因此，固体边界在腐蚀与沉淀的影响下不断更新。为了能更好地实时追踪VOP新的固体边界，采用 方法予以解决。根据式（1）反应，金属基体表面结点的体积变化［19］如下： dV （20） Fe =- AVˉFe k CO 1 dt根据式（2）反应，新生成的沉淀结点在固体表面的体积变化如下： dV （21） Fe(OH) = AVˉFe(OH) k 2C dt Fe3式（20）~式（21）中：V Fe和V分别为金属Fe Fe(OH)3(s) Fe（OH）3与沉淀 结点的无量纲体积； VˉFe(s)和VˉFe(OH)3(s)分别为金属Fe Fe（OH）3与沉淀 的摩尔体O2 Fe3+的积； CO 和 C 分别为腐蚀边界上 和Fe3 为式（1）或式（6）浓度；A为固体结点的表面积；k1 4）反反应中的腐蚀溶解反应速率； k 为式（ 应的2
沉淀反应速率；通过式（23）~式（24）计算腐蚀边界O2 Fe3+的浓度。式（20）和式（21）的离散格组分 与式为： （22） V ( t +D t )= V ( t )- AVˉFe k 2C Dt Fe Fe O2 V ( t +D t )= V ( t )+ AVˉFe(OH)3 k 2C Dt（23） Fe(OH)3 Fe(OH)3 Fe3 22），当根据式（ V ( t )的值减小到一定的数Fe值，即该固体结点完全被腐蚀，固相结点转化为液体结点；当V ( t )的值增大到一定的数值，即该Fe(OH)3液体结点完全沉淀，液相结点转化为固体结点。式（1）反应的金属与腐蚀溶液在边界处发生腐蚀，腐蚀边界处各组分浓度变化［7］为： ¶CO （24） D 2(aq) = k C O ¶n 1 O 2(aq) 2(aq) ¶C +( （25） Fe3 D aq) =- k C Fe 3+ ¶n +( 1 Fe3 ( aq) aq) O2式中： CO 为腐蚀边界上 的消耗浓度； C Fe3 +( 2(aq) aq) Fe3+的生成浓度；n为腐蚀边界上 为反应固体表面O2的外法线方向； DO 为 扩散系数；D为Fe3 + 2(aq) ( aq) Fe3+扩散系数。
1.5 结点钝化模型
模拟过程中的腐蚀边界，每个结点都有一定的概率被钝化，导致该处的腐蚀速率变慢，使腐蚀坑的形状呈现不对称性，这也符合真实的腐蚀情VOP况。在 方法中，腐蚀边界处结点发生体积的变化，如果腐蚀边界某一结点发生钝化，该结点的体积减小速率会变慢。结点发生钝化的概率称为钝化概率P，钝化概率系数ε为0( R > P) （26） ε= 1( R P) 0~1式中，R为 之间的任意随机实数。因此，需对式（22）进行修正： （27） V ( t +D t )= V ( t )- εAVˉ k C Dt Fe Fe Fe 1 O VOP在 方法中引入钝化概率函数，可以避免腐蚀界面上的阳极金属受到的腐蚀作用相同，以使腐蚀坑的形貌随机变化且不对称，更加符合真实的腐蚀情况。
2 数值模拟 2.1 物理模型
LB基于上文构建的 腐蚀模型，对金属表面点100×160蚀稳态过程进行模拟。计算区域选取 的网格。整个计算域内，初始的初生蚀孔假设为一71，80），半 10，钝部分圆，其圆心位于（ 径为 化膜 5，如图1（a）所示。的厚度为
2.2 边界条件
对于液体四周的边界采用了非平衡外推格式进行处理，而溶液与金属之间、溶液与钝化膜以及溶液与腐蚀产物沉淀之间的流动碰撞采用反弹边界进行处理。
2.3 模拟工况
本文模拟金属表面点蚀稳态过程的形貌的演1化过程，参数如表 所示，表中 v 为运动粘度。采用变参数的数值模拟得到腐蚀程度的影响因素，各2所示［3，14］。所有扩散系数均由各浓工况参数如表度场的松弛时间决定，如式（14）所示。为了满足0.5。运行过程中的数值稳定性，选取数值需大于 上述工况中所有参数均为无量纲参数，均通1过表 中的转化关系得到。由于化学反应速率受温度、反应接触面积、催化剂等影响，且腐蚀介质扩散系数受浓度、温度等影响，这两种参数可选取值1广［4，15］，因此选取参数时，仅在表范围较 中的模拟数值附近选取。
3 数值模拟结果与分析 3.1 腐蚀演化过程分析
通过模拟得到不锈钢材料表面点蚀稳态过程的形貌变化以及各组分浓度变化情况。模拟工况=8.14 × 10-5，C =1 134.7， =0.028，为 k D 1 O O =0.027 5。图2 D 所示为不锈钢材料点蚀坑形Fe貌随时间变化的情况。从图中可以看到，腐蚀破损区域逐渐变大，并且逐渐向金属基体纵深发Fe（OH）3展。当蚀孔出现 沉淀，随着时间的推移2（c））；除了初生的沉淀增多，逐渐覆盖了蚀孔（图蚀孔发生腐蚀，在蚀孔底部出现次生蚀孔继续向2（d）~图 2（f））。金属基体深处腐蚀（图 3图 所示为含氧腐蚀溶液的浓度随时间变化 的情况。随着时间的延续，腐蚀的区域变大，并向金属基体纵深发展，腐蚀溶液向蚀坑内渗入，腐蚀溶液逐渐被消耗，在蚀孔处出现沉淀，阻碍了腐蚀溶液向蚀孔内扩散，导致蚀孔内腐蚀溶液的浓度逐渐降低。4图 所示为蚀孔内阳极金属发生腐蚀溶解Fe3 +的浓度变化情况。随着时间的时，腐蚀产物延续，腐蚀产物在腐蚀边界不断生成，并且向蚀坑Fe3 +扩散至蚀孔与孔外外以及液体环境扩散；当OH- Fe（OH）3，Fe（OH）3大量的 结合产生 达到其饱Fe3和浓度继而沉淀。蚀孔产生沉淀，阻碍了 +向Fe3+升高，导致Fe3+进一步发生孔外扩散，孔内的水解反应。5图 所示为蚀孔外钝化膜作为阴极及其表面OH- OH-产生的 浓度变化及 扩散的情况。随着时
间的延续，蚀孔内阳极金属发生失电子的腐蚀溶
解反应，电子通过金属基体转移至电位较高的钝化膜，O2在钝化膜表面得到电子，发生阴极反应，在OH- ；OH-表面产生大量的 ，并向液体环境扩散 扩Fe3+生成 Fe（OH）3散至蚀孔处与 发生沉淀反应，孔OH-向蚀孔内的扩散，导致孔口沉淀的增多阻碍了OH-浓度降低。内 6 H+图 所示为 浓度随时间的变化情况。由于Fe3 +向孔外扩散蚀孔处的沉淀不断增多，阻碍了OH- Fe3+聚集，继以及 向孔内扩散，导致孔内大量H+而发生水解反应，产生大量的 ，初生蚀孔底部H+浓度变大，进一步腐蚀金属基体，金属加速溶解，产生次生孔，推动点蚀的稳定发展，蚀孔逐步
向深处生长。
7图 所示为整个腐蚀过程中被腐蚀的金属体t=0~20 000积随时间的变化关系。当 时，金属腐蚀仅发生初生蚀孔的吸氧腐蚀，初生蚀孔被腐蚀的体积V逐渐变大，蚀孔的深度h逐渐增大；当t>20 000时，由于孔内闭塞及酸化的腐蚀环境，导致金属自催化腐蚀反应的发生，腐蚀速率大于先前单独发生的吸氧腐蚀的速率，并且V随 t的增加而增大，h也逐渐增大。
3.2 影响因素分析 3.2.1 腐蚀化学反应速率
模拟过程中，通过改变式（1）反应的化学反应速率，得到某一时刻（t=50 000）不同腐蚀化学反应8速率下点蚀坑的形貌（图 ）。从图中可以看出，在相同腐蚀时间下，随着腐蚀化学反应速率的降
低，腐蚀程度减小。8通过对图 中 V 和 h进行统计测量，得到t=50 000 9）。从图9可时V和h与 k1的关系曲线（图就越大。由于式（1）以看出，随着 k1 增加，V越大，h中化学反应速率的增加，导致了吸氧腐蚀加快，腐Fe3+浓度增大，加速了 Fe3+水蚀产物增多，蚀孔内 H+解过程，产生更多的 ，使孔内阳极金属进一步加速腐蚀，因此呈现出以上结果，并且该结果与真 实情况相符。 8与图9综合分析图可知，由于 k1 增加，单位时间内V变大，h增加。
3.2.2 腐蚀溶液扩散系数
通过改变模拟过程中浓度场含氧溶液的扩散
系数，研究腐蚀反应中介质间流动与扩散对腐蚀10 刻（t=50 000），不的影响。图 给出了某一时 同腐蚀溶液扩散系数下腐蚀坑的形貌。从图中可以看出，相同的腐蚀时间下，随着 D 的增大，腐蚀程O2
度加重。 10通过对图 中V 和 h进行统计测量，得到t=50 000 11时V和 h与 D 的关系曲线（图 ）。从O2 11图 可以看出，随着 D 的增大，V越大，h越大。O2由于 D 的增加，导致腐蚀溶液向蚀孔内扩散速O2率加快，金属腐蚀边界的浓度相对较高，加快了腐蚀反应，因此，呈现出以上模拟结果。10 11综合分析图 与图 可知，由于 D 增大， O2导致腐蚀溶液加速向蚀孔内扩散，致使腐蚀结点处的腐蚀溶液浓度相对增加，加快了腐蚀反应。由此得到的结论是：D 越大，V越多，h越大。O2
3.2.3 腐蚀产物扩散系数
Fe3+的扩散在模拟中改变浓度场的腐蚀产物系数，研究腐蚀产物扩散系数对腐蚀程度的影12所示为某一时刻（t=50 000）不同腐蚀产响。图
Fe3+物 扩散系数下腐蚀孔的形貌特征。从图中可以看出，在相同的腐蚀时间下，随着 D 的减Fe3
小，腐蚀程度加重。12通过对图 中V 和 h进行统计测量，得到t=50 000 13）。从时V和 h与 D 的关系曲线（图Fe3 13图 可以看出，随着 D的降低， V越大， h越Fe3 Fe3+的大量聚大。由于 D 的降低，导致蚀孔内Fe3 Fe3+集，加速了 的水解过程，产生了更多的氢离子，使孔内的酸度增加，金属基体的腐蚀溶解速率增加。因此，呈现出以上模拟结果。 12 13综合分析图 与图 可知，D增加，将加Fe3 Fe3+水解过程，使孔内酸度增加，加快腐速蚀孔内蚀反应。由此得到的结论是，D 越小，V越大，h Fe3
越大。
4结论
LB本文基于改进后的 腐蚀模型，将钝化概率VOP函数引入 方法中，对腐蚀化学反应速率进行LB修正，建立了 腐蚀模型。使用该模型，模拟了浸没于液体腐蚀环境中的金属表面单坑点蚀稳态过程，通过模拟得到以下结论： 1）LBM算法可以明确得到溶液中不同组分的流动与扩散，由此得到其浓度分布，通过流场和浓度场之间的耦合模拟真实的液体腐蚀环境；引入钝化概率函数使得模拟的腐蚀坑的形貌随机变化且不对称，更加符合真实的腐蚀坑形貌；阴、阳极共同作用使电化学腐蚀的模拟更加完善；VOP方法可以控制腐蚀的快慢，而不是简单的固、液相随机转换。2 ）腐蚀反应速率越大，在相同的腐蚀时间下，在金属表面点蚀的稳态过程中，金属被腐蚀的体积越大，腐蚀的深度就越大；腐蚀含氧溶液扩散系数越大，在相同的腐蚀时间下，金属被腐蚀的体 Fe3+扩散系积越大，腐蚀的深度就越大；腐蚀产物数越大，在相同的腐蚀时间内，金属被腐蚀的体积越小，腐蚀的深度就越小。
参考文献：