-
本工作利用腐蚀电化学计算腐蚀界面能, 构建了锆合金腐蚀过程的相场模型. 首先, 利用所建立的模型模拟了Zr-2.5Sn合金均匀腐蚀过程, 模拟结果显示该合金的腐蚀动力学曲线符合立方规律, 与实验结果一致. 分析发现, 在氧化层生成的初期, 氧化层生长速率很高, 但是受温度的影响不明显; 随着氧化层厚度的增长, 温度对氧化层生长曲线的影响变大, 温度越高腐蚀速率越快. 多晶Zr-2.5Sn合金腐蚀行为的模拟结果表明, 在锆合金基体晶界处由于具有更大的氧扩散速率, 氧化速率加快, 并在金属-氧化层界面朝向金属基体一侧形成了沿晶界的具有更高浓度的O2–带, 且对氧化腐蚀速率的影响主要表现在氧化初期, 相场模拟获得的腐蚀动力学曲线与实验结果符合非常好.
-
关键词:
- 相场模拟 /
- Zr-2.5Sn合金 /
- 腐蚀动力学 /
- 晶界腐蚀
Due to the small neutron absorption cross section and excellent thermal creep performance, zirconium alloy is one of the most important cladding materials for fuel rods in commercial fission reactors. However, quantitative analysis of the effects of temperature and grain boundaries on the corrosion microstructure evolution of zirconium alloys is still needed. The establishing of a phase field simulation for the corrosion process of polycrystalline zirconium alloy and the systematical investigating of the thermodynamic influence are both very important. In this study, the phase field model of the corrosion process in zirconium alloys is developed by combining corrosion electrochemistry through calculating the interfacial energy at the metal-oxide and oxide-fluid boundaries. Then the model is used to investigate the uniform corrosion behavior on the surface of Zr-2.5Sn alloy, which demonstrates that the corrosion kinetic curve follows a cubic rule. Subsequently, the influence of temperature on the corrosion thickening curve of zirconium alloy is examined, and good agreement between simulation and experimental results is achieved. It is observed that during early stage of oxide layer formation, there is a high growth rate with minimal temperature dependence; however, as the oxide layer thickness increases, temperature becomes a significant factor affecting its growth rate, with higher temperatures resulting in faster corrosion rates. Furthermore, the effect of polycrystalline zirconium alloy matrices on corrosion rate is investigated, revealing that the grain boundaries accelerate oxide layer thickening due to enhanced oxygen diffusion rates. At metal-oxide interface, O2– bands are formed in areas with higher O2– concentration along these grain boundaries towards the metal matrix, which mainly influences oxidation-corrosion rate during the initial oxidation stage.-
Keywords:
- phase-field simulation /
- Zr-Sn alloy /
- corrosion kinetics /
- grain boundary corrosion
1. 引 言
锆合金具有中子吸收截面小、热蠕变性能好等优点, 是目前商用裂变反应堆燃料棒最常采用的包壳材料[1]. 锆合金包壳长期暴露于高温、高压(300—400 ℃, 16 MPa)和强中子辐照的反应堆一回路水环境中, 因此在包壳外壁经常发生水侧氧化腐蚀, 严重限制其使用寿命[2]. 锆合金腐蚀由O2–在氧化层中的扩散控制, 其耐腐蚀性能取决于ZrO2氧化层对O2–和金属基体的隔离能力[3,4], 腐蚀增厚曲线表现为抛物线规律或立方规律[5]. 近年来, 国内外关于锆合金抗腐蚀性能的研究主要集中于ZrO2氧化层失效的机理性分析[6,7]. 然而, 从实验角度再现锆合金服役工况并对温度和辐照等条件进行对比研究的成本很高[8], 而且锆合金腐蚀周期很长, 实验成本较高. 因此, 对锆合金腐蚀行为展开理论模拟研究对揭示其腐蚀机理具有重要意义.
研究表明, 温度是影响锆合金腐蚀行为的重要因素[9–12]. 实验发现, 随着温度升高(50—100 ℃的升温), 常见Zr合金的腐蚀速率会增大, 造成这一现象的原因是O2–扩散系数随温度升高而提升[9,10]. 锆及其合金含氧体系的热力学自由能随温度的变化非常剧烈[11,12], 而目前的模拟工作将自由能参数简化为常数[13], 无法准确揭示温度对锆合金腐蚀行为的影响. 另一方面, 晶界是抵抗材料腐蚀的薄弱区域. 在高温高压和强辐照等条件下, 金属材料往往会出现晶间腐蚀甚至应力腐蚀开裂的严重问题[14]. 目前, 大部分关于锆合金腐蚀的模拟工作采用了均匀化假设, 无法模拟锆合金基体多晶微结构影响下的腐蚀形貌演化.
介观尺度的相场方法是模拟材料组织转变的重要工具, 近年来已经被成功应用于材料腐蚀形貌演化的研究[13,15–19]. Asle和El Kadiri [13]将腐蚀电化学融入相场模型的自由能势垒, 采用单一序参量描述了腐蚀层的一维生长, 对Zr在900 ℃条件下暴露在常压空气中的腐蚀行为进行了模拟; Mai和Soghrati [18]将相场模型的自由能势垒与腐蚀电化学相关联, 模拟了不锈钢的均匀腐蚀过程, 进而研究了应力对腐蚀速率的影响. Fang等[19]通过在腐蚀界面上引入感应微电流密度, 建立钢筋混凝土结构失效的腐蚀-扩散相场模型. 然而, 目前已有的锆合金腐蚀相场模型依然缺少对温度的定量分析, 也无法确定晶界等微结构对锆合金腐蚀行为的影响. 此外, Zr-Sn合金中合金元素Sn的添加对腐蚀行为的影响机理有待深入研究[10], 然而目前针对锆合金腐蚀过程的相场模型中, 关于合金元素对锆合金自由能参数的影响大多做了简化处理[13]. 因此, 构建多晶锆合金腐蚀过程的相场模拟, 系统研究腐蚀温度等热力学因素对腐蚀动力学的影响至关重要.
本工作结合锆合金热力学数据, 构建了描述锆合金腐蚀行为的相场模型, 引入锆合金腐蚀电化学能优化了相场模型. 之后研究了温度对Zr-2.5Sn合金腐蚀的影响, 绘制了对应温度下的腐蚀增厚曲线. 最后研究了引入锆合金金属基体晶界微结构后, 温度对腐蚀行为的影响. 模拟结果将有助于建立更精确地锆合金腐蚀动力学模型, 能够应用于锆合金腐蚀行为的机理性分析.
2. 物理模型
2.1 相场模型
图1所示为金属腐蚀组织的示意图及本工作的相场变量的示意图. 相场模型中, 将模拟区域分为金属相、氧化层相和腐蚀介质相三相. 氧化层中氧从腐蚀介质中定向朝向金属流动, 以及阳离子从金属内部朝外定向流动. 从图1(a)可以看出, 本工作采用非保守型相变序参量η1, η2和η3分别描述金属相(M相)、氧化层相(O相)和腐蚀介质相(F相). 3个相场变量的示意图如图1所示, η1 = 1时表示金属相, η1 = 0时则表示在氧化层相或腐蚀介质相, η1的取值在0—1之间连续变化. 对氧化层相η2和腐蚀介质相η3做同样处理, 即在所描述的物相内取值为1, 而在其他相的取值为0, 并且取值在0—1之间连续变化. 由于在阳极氧化物Zr4+在氧化层中的移动能力非常低[20], 因此本模型以氧离子O2–浓度cO作为保守型相场变量.
相场模型的总自由能方程为
F(cO,ηi)=∫V[fbulk(cO,ηi)+fgrad(ηi)+Eelas(cO,ηi)]dV, (1) 其中, fbulk为体自由能密度函数; fgrad 为梯度能密度函数; cO为O2–相对浓度; Eelas为弹性能, 由计算应力应变平衡方程后根据弹性模量和应力应变分布导出.
体自由能密度函数fbulk为
fbulk(cO,η)=3∑i=1higi(cO)+W12η21η22+W23η22η23, (2) 其中, W12及W23为针对氧化层-金属基体(O-M)界面及氧化层-腐蚀介质(O-F)界面的腐蚀电化学能势垒高度, 分别作用于O-M界面及O-F界面位置处; h1, h2, h3分别为对应相成分的插值函数, g1, g2, g3为表示金属、氧化层和腐蚀介质相的能量密度函数.
梯度能密度函数fgrad为
fgrad(ηi)=κ12∇η1⋅∇η2+κ23∇η2⋅∇η3, (3) 其中κ12与κ23分别为O-M界面与O-F界面的界面能系数, 其取值与界面能量和界面宽度相关[21].
本工作所使用的插值函数hi为 [15]
hi=η3i(6η2i−15ηi+10). (4) 该函数满足ηi相在稳态点ηi = 1和ηi = 0位置处同时满足hi(1) = 1, hi(1) = 0, h′i(1)=h′i(0)=0, 使得相成分在相界面处具有光滑且一阶导数连续的数学特性.
本工作以二元锆合金Zr-2.5Sn[22]作为模拟对象, 其自由能密度函数为
g1=3∑a=1caEa+3∑a=13∑b>aAabcacb, (5) 其中, a, b = 1, 2, 3分别表示O, Zr, Sn元素, ca表示成分a的相对含量(%); Ea为纯组分a的Gibbs自由能; Aab为二元组分a与b的混合系数.
Dinsdale[23]编写的SGTE元素热力学数据库报告提供了Zr, Sn, O的纯金属自由能密度函数参数与温度的关系. 在本工作的模拟温度范围内, Zr-Sn相图位于α-Zr相(HCP)[24], Zr-O相图位于α-Zr相(HCP)[25], Sn-O相图位于SnO2相[26], 根据热力学数据库引入了其二元组分混合系数[24–27].
Asle等[28]的工作中通过模拟O2–偏离饱和氧密度的差值来描述腐蚀介质相和氧化层相的自由能密度. 含O2–的氧化层、腐蚀介质相的相自由能密度函数g2和g3分别为
g2(cO)=g02+12k2(cO−cs2)2, (6) g3(cO)=g03+12k3(cO−cs3)2, (7) 其中k2和k3分别为O2–偏离饱和密度带来的自由能密度系数; cs2为氧化层相中饱和氧密度, 可由Zr-O的相图得到[25]. cs3为O-F界面腐蚀介质相侧的氧平衡浓度, 取决于一回路水的氧溶解度. 参考Zr-Sn-O相图[23], 本工作中cs2 = 0.667. g02为氧化层中氧含量为饱和氧密度cs2时, 纯ZrO2的自由能密度 [24]. g03为高温高压水中氧含量为饱和氧密度cs3时, O2的自由能密度. 本工作中应用到的Zr-Sn合金热力学参数如表1所示.
表 1 Zr-Sn合金自由能参数Table 1. Free energy parameters of Zr-Sn alloys.Parameter Value Reference EZr/(J·m–3) −7827.6+125.65T−24.1618TlnT−4.38e−3T3+34971T−1 [23] ESn/(J·m–3) 2524.7+4.00T−8.26T\lnT−16.81e−3T2+2.62e−6T3−1.08e6T−1 [23] EO/(J·m–3) −3480.9−25.50T−11.13TlnT−5.10e−3T2+0.66e−6T3−38365T−1 [23] AZrSn −148022.5+19.41T+(173681.9−22T)(cZr−cSn)+104271.96(cZr−cSn)2 [24] AZrO −37876.66+17.2915T−4471.4(cZr−cO) [25] ASnO 140878 – 23.9326T [26, 35] g02 −1117869+420.3T−69.6TlnT−0.003766T2+702910T−1+4.59e−21T7 [25] g03 0 k2/(J·m–3) 2.5 × 105 [13] k3/(J·m–3) 2.5 × 105 [13] 对非保守型相场序参量ηi 采用Allen-Cahn方程进行求解[29,30]:
∂ηi∂t=−LδF(cO,η)δηi=−L×[δf(cO,η)δηi+δEelas(cO,η)δηi−3∑j≠iκij∇ηi⋅∇ηj], (8) 其中L为非保守场变量ηi的相变迁移率.
对保守型相场变量O2–浓度序参量采用Cahn-Hilliard方程求解[31]:
∂cO∂t=∇⋅[Mc(∇δF(cO,η)δcc)]−P∂η2∂t=∇⋅[Mc∇(δf(cO,η)δcc+δEelas(cO,η)δcc)]−P∂η2∂t, (9) 其中P为单位M相转变为O相所消耗的O2–. 由于电荷守恒, 在O-M界面消耗O2–产生氧化层相后, O-F界面产生等量的O2–; Mc为O2–的迁移率, 其值由相变序参量插值函数h1, h2, h3对各相中O2–迁移率M1, M2, M3插值得到. 此外, 晶界扩散也会影响O2–在基体中的扩散, 因此需要考虑晶界迁移率MGB. 在锆合金腐蚀行为中, Zr, Sn的扩散速率大大小于O的扩散速率[32], 因此可忽略计算Zr, Sn元素的相场Cahn-Hilliard扩散方程. 本工作的迁移率系数为
Mc=h1M1+h2M2+h3M3+16h1⋅(1−η1)MGB. (10) 本工作中采用到的不同温度下O2–的扩散系数如表2所示[33,34]. 各相中O2–迁移率通过Einstein关系式由扩散系数推出, 推导关系式为Mi=Di(T)/(RT), 其中R为热力学常数, T为温度, Di(T)表示温度为T时i相中氧扩散系数.
Temperature/K D1/(m2·s–1) D2/(m2·s–1) D3/(m2·s–1) DGB/(m2·s–1) 633 3.99×10–17 1.61×10–17 2.98×10–17 1.46×10–16 653 1.26×10–16 5.18×10–17 8.77×10–17 8.11×10–16 673 3.09×10–16 1.11×10–16 2.25×10–16 3.76×10–15 693 1.04×10–15 3.69×10–16 7.64×10–16 2.82×10–14 Eelas(c,η)=12∫V∑pqrsCpqrs(η)(εpq−ε0pq)×(εrs−ε0rs)dV=12σ:ε, (11) 其中εpq与εrs为不同方向的总应变; ε0pq与ε0rs为不同方向的本征应变. 总应力矩阵σ和总应变矩阵ε遵循张量形式Hooker定律σ = C(η)ε. 总应变ε和总应力σ为3×3的矩阵形式, 弹性模量的张量形式C则为3×3×3×3的张量形式, 其分量为Cpqrs. 3×3×3×3的张量形式C可整理为6×6的矩阵形式, 其分量为Cmn. 通过计算力学本构, 实现了将弹性应变能引入腐蚀相场模型. 本征应变的取值体现包括非弹性应变等要素在内的组织受力的基本特性, 本工作目前的模拟中本征应变取0. 沿pq方向的总应变εpq为
εpq=12(∂up∂rq+∂uq∂rp)−εc, (12) 其中εc为由于氧成分导致的成分应变. 参考Asle和El Kadiri[13]的工作, 氧成分导致的成分应变εc为:
εc=h1e1(cO−cs1)+h2e2(cO−cs2)+∂η1∂rRPB(cO−cs1), (13) 其中e1, e2分别为由于金属相和氧化层相的氧平衡浓度偏离而导致的成分应变率, 本工作中e1, e2取0.002 [35]; cs1为金属相中饱和固溶氧浓度, 参考Zr-Sn-O相图[36], 本工作中cs1 = 0.281; RPB = 1.56为锆合金氧化体积膨胀比(PB比); r为氧化层朝向金属方向, (13)式右侧第3项表示由于锆氧化为ZrO2导致体积膨胀而带来的应变.
本工作的弹性模量C(ηi)为
C(ηi)=C1(η1)h1+C2(η2)h2+C3(η3)h3, (14) 其中C1(η1), C2(η2), C3(η3)分别为金属相、氧化层相、腐蚀介质相的弹性模量张量形式, 其6×6的矩阵形式分量为Cimn.
(1)式中弹性能Eelas对相场序参量和成分变量的导数分别为
∂Eelas∂ηi=12∂hi∂ηi∑pqrsCpqrsεpqεrs, (15) ∂Eelas∂c=12∑pqrsCpqrs(h1e1+h2e2+∂η1∂rR). (16) 应力应变的分布由应力-应变平衡方程给出:
∇⋅Cpqrs(εpq−ε0pq)=0. (17) 本工作中采用的金属基体及氧化层中弹性模量Cmn的取值及参考文献列于表3中[37–39].
2.2 数值方法与参数
为了便于数值模拟, 对相场方程中的变量和常数进行无量纲化处理. 无量纲化过程包括: r∗=r/Δl, t∗=tD/(Δl)2,M∗=MRT/D,κ∗12=κ12/[RT(Δl)2] , κ∗23=κ23/[RT(Δl)2] 其中Δl为空间离散步长, Δl = 10 nm. 得到(9)式和(10)式的无量纲化形式分别如下:
∂ηi(r∗,t∗)∂t∗=−LδF∗(cO(r∗,t∗),η(r∗,t∗))δηi(r∗,t∗), (18) ∂cO(r∗,t∗)∂t∗=∇∗⋅[M*c(r∗,t∗)(∇∗δF∗(cO(r∗,t∗),η(r∗,t∗))δcO(r∗,t∗))]−P∂η2(r∗,t∗)∂t∗. (19) 3. 结果与讨论
3.1 表面均匀腐蚀与腐蚀增厚曲线
图2为673 K下Zr-2.5Sn合金表面均匀腐蚀25 d后的氧化层的分布. 图中蓝色、绿色、黄色区域分别表示腐蚀介质相、氧化层相, 金属相. 模拟区域为1024×128网格, 网格宽度为空间离散步长Δl. 垂直于氧化层生长的方向取周期性边界, 左右方向取零通量边界. 在模拟初始时刻, 设置厚度为一个网格的氧化层隔离腐蚀介质相与基体相. 腐蚀2.5 d之后, 氧化层厚度为0.761 μm, 平均腐蚀速率约为0.304 μm/d. 随着腐蚀的进行, 氧需要从腐蚀介质相中穿过氧化层抵达O-M界面处. 腐蚀25 d之后, 氧化层厚度为1.902 μm, 平均腐蚀速率约为0.077 μm/d. 可见, 随着腐蚀的进行, 平均腐蚀速率的增速有所降低, 主要原因是: 随着氧化层的增厚, O2–扩散至O-M界面所需要的时间逐渐延长, 成为限制氧化层生长的主要因素. 模拟结果得到的腐蚀速率与Jiang等[9]统计的400 ℃下Zr-2.5Sn增厚速率接近.
图3为腐蚀层随时间的增厚曲线. 可以看出, 相场模拟计算得到的Zr-2.5Sn的腐蚀增厚曲线的与立方规律拟合较好, 拟合的R-Square为0.99325. 堆内实验表明[9], 常见的锆合金在腐蚀发生一个月后陆续发生腐蚀转折; 在腐蚀转折发生前, 锆合金的腐蚀行为呈现抛物线或立方规律[5]. O2–扩散到O-M界面的效率是决定锆合金腐蚀增厚曲线形状的关键因素; 在扩散控制下, 腐蚀增厚曲线表现为抛物线或立方规律[5]. 对比增厚曲线与拟合曲线的差异, 发现在腐蚀5 d以后的拟合结果相对5 d内的更好, 这是由于锆合金的腐蚀主要是受O2–扩散控制[5]. 在本模型的模拟初期ZrO2氧化层较薄, O2–扩散到O-M界面较快, 故在腐蚀初期腐蚀增厚曲线与立方规律拟合存在偏差[31].
3.2 温度对腐蚀增厚曲线的影响
图4为不同温度条件下腐蚀25 d的氧化层的分布图. 图中蓝色、绿色、黄色区域分别表示腐蚀介质相、氧化层相和金属相. 模拟区域为2048×256网格, 网格宽度为空间离散步长Δl. 垂直于氧化层生长的方向取周期性边界, 左右方向取零通量边界. 可以看出, 相同的腐蚀时间下, 腐蚀温度越高, 氧化层越厚. 当腐蚀温度从633 K上升到693 K, 氧化层的厚度从1.631 μm上升到2.147 μm, 氧化层的生长速率也从0.065 μm/d上升到0.086 μm/d, 工作温度的提升将加快ZrO2氧化层厚度的生长. Liu[40]的研究表明, 温度升高, 腐蚀增厚大幅增加, 锆合金包壳元件表面工作温度需要限制在400 ℃以下.
图5为不同温度下的腐蚀增厚曲线. 在氧化层生成的初期, 存在着较为明显的高速率氧化层生成过程, 且该过程的氧化层生长速率受温度影响不明显, 因而当腐蚀时间小于2 d时的增厚曲线几乎是重叠的. 氧化层的厚度生长到一定程度后, 由于O2–扩散到O-M界面的时间受到温度的影响, 导致温度越高扩散越快氧化层生长越快, 生长速率受到温度差异的影响. 温度的升高对O-F界面的移动影响并不明显, 在633—693 K范围内, O-F界面几乎处于平齐位置, 基本不受温度影响.
图6 为不同温度下O2–浓度沿氧化层生长方向的分布. O2–的扩散对锆合金氧化层的增厚有很大影响, 而温度是影响氧化层中O2–扩散的关键因素. 氧化层中O2–浓度始终大于氧化层相中饱和氧浓度cs2, 且越接近O-M界面, 其值越接近cs2. 随温度的升高, 整个氧化层中含有更多的O2–, 升温加快了O2–进入氧化层的速率. 此外, 随着温度的升高, O2–在氧化层中的扩散速率加快, 在O-M界面形成更高的O2–浓度, 相场自由能提升, 金属相向氧化层相的转化加快. 因此, 温度升高导致腐蚀速率加快是扩散系数加快和自由能势垒增大两方面因素共同导致的结果.
3.3 温度对多晶锆合金腐蚀的影响
图7为模拟多晶Zr-2.5Sn合金腐蚀过程的初始组织. 采用晶粒长大相场模型[41]的模拟结果得到的多晶结构作为锆合金基体M相的初始微结构, 组织形貌及每个晶粒的取向如图7所示. 模拟区域为256×1024网格, 每个网格取网格宽度为10Δl, 垂直于氧化层生长的方向(图7中左右方向)取周期性边界, 上下方向则取零通量边界. 所选组织中选用了平均晶粒直径为10.82 μm, 与工程用锆锡合金的平均晶粒尺寸基本一致[42].
图8为多晶Zr-2.5Sn合金在673 K腐蚀不同时间后的相场模拟结果. 可以看出, 随着腐蚀时间的延长, 腐蚀层的厚度不断增大. 晶界对氧化层的形貌有重要影响. 如图8(a3)中蓝色箭头部分所示位置, 晶粒内相对沿晶界方向的腐蚀层生长更慢. 由于O2–沿晶界的扩散较快, 形成了一段沿着晶界的O2–带, 深度为~1.2 μm. O2–沿晶界扩散深入到基体中, 是晶界导致氧化层生长加快的关键原因. 在金属材料的腐蚀中也观测到了元素扩散沿晶界加快[43]和沿着晶界的扩散元素浓度带[44]. O2–的含量在腐蚀初期相对更高, 在腐蚀5 d后的图8(b1)中, O2–带具有更高的亮度, 而腐蚀达到50 d后, O2–带的亮度变暗. 在腐蚀5 d后, 晶界处O-M界面的移动速率加快已经形成向下的“凸起”. 然而晶界处的“凸起”深度并未随着氧化时间加长而更加明显, 腐蚀时间达到20 d后“凸起”深度变化不大. 可以发现, 晶界导致的锆合金腐蚀速率加快主要作用于腐蚀发生的最早期阶段. Liu等[45]对Zr-4合金的腐蚀实验发现, 腐蚀层厚度达到2.6 μm时, 在O-M界面处发现了沿晶界朝向金属基体方向的突起, 突起平均长度小于1 μm. 这与本工作的模拟结果一致.
图9是不同温度下腐蚀25 d后多晶Zr-2.5Sn合金的相场模拟结果. 可以看出, 温度越高, 腐蚀速率越快, 晶界处的腐蚀速率增加明显大于晶粒内腐蚀速率, 且在O-M界面沿晶界形成的“凸起”也更加明显. 如图9(b)中红色平行线所示, 在温度653, 673, 693 K下, O-M界面的“凸起”纵向深度分别约为0.37, 0.51, 0.62 μm. 此外, 随着温度的升高, O-M界面偏向金属基体内的O2–带也越长. 温度升高后, 在晶界处的O2–将更容易扩散到金属基体中去, 从而在该位置形成更明显的O-M界面“凸起”. 综合图9(a), (b), 可以认为温度升高将会扩大锆合金基体中晶界对腐蚀的影响. 温度升高将导致O2–在晶界上的扩散系数增大更明显, 形成更深入的O-M界面“凸起”.
图10为不同温度下的腐蚀增厚曲线及其与实验结果的对比. 图中离散点为Jiang等[9]的实验结果. 可以看出, 增厚速率随着温度的升高而升高. 与锆锡合金单晶腐蚀的模拟结果(图5)相比, 引入多晶结构后的增厚曲线增厚了0.2—0.3 μm. 可见, 晶界加速了腐蚀行为. 在633, 653, 673和693 K下, 腐蚀25 d后的增厚分别为1.778, 2.017, 2.163和2.328 μm. 模拟结果与Jiang等[9]统计的关于Zr-2.5Sn的633 K和673 K下腐蚀实验结果匹配较好. Liu等[40]的研究也表明, 在腐蚀约一个月后表面氧化层的厚度达到约2 μm, 与本文的模拟结果也一致.
分析发现, 不同温度下氧化层增厚速率的差异主要表现在腐蚀初期. 在腐蚀5 d之后, 3条曲线之间的差距变小. 在腐蚀初期, 氧化层中O2–扩散受氧化层厚度限制较小, 温度升高对晶界扩散带来的影响相较于晶粒内扩散更大, 腐蚀速率随温度提升得更明显. 腐蚀层达到一定厚度后, O2–扩散至O-M界面需要的时间更长, 此时扩散成为限制氧化层进一步生长的主要原因, 晶界处扩散速率的提升已不再导致O-M界面的移动速率的明显加快. 因此, 晶界对腐蚀速率的影响主要表现在氧化层较薄的腐蚀初期, 而温度的提升会加剧晶界对腐蚀速率的影响.
4. 结 论
1)建立了Zr-2.5Sn二元合金的腐蚀过程的相场模型, 将腐蚀电化学能与相场界面参数关联起来, 模拟了表面均匀腐蚀行为. 腐蚀动力学曲线显示, Zr-2.5Sn合金的均匀腐蚀符合立方规律.
2)研究了温度对锆合金腐蚀增厚的影响, 发现温度越高腐蚀速率越快. 氧化层厚度的增长受O2–扩散控制, 在氧化层生成的后期, 温度对氧化层生长曲线的影响更大.
3)研究了多晶Zr-2.5Sn合金的腐蚀形貌演变, 并获得了腐蚀动力学曲线, 腐蚀增厚曲线与实验结果符合得非常好. 晶界上的O2–扩散速率更快, 因此晶界加速了腐蚀行为, 在腐蚀初期晶界扩散对氧化膜的生长的影响非常大.
[1] Kai J J, Huang W I, Chou H Y 1990 J. Nucl. Mater. 170 193
Google Scholar
[2] Zhao W Z, Quan Q W, Zhang S M, Yang X Y, Wen H Y, Wu Y C, Liu X B, Cao X Z, Wang B Y 2023 Radiat. Phys. Chem. 209 110986
Google Scholar
[3] Jones R H, Simonen E P 1994 Mater. Sci. Eng. 176 211
Google Scholar
[4] Hu J, Liu J L, Lozano P S, Grovenor C R M, Christensen M, Wolf W, Wimmer E, Mader E V 2019 Acta Mater. 180 105
Google Scholar
[5] Cann C D, So C B, Styles R C, Coleman C E 1993 J. Nucl. Mater. 205 267
Google Scholar
[6] Wei T G, Dai X, Long C S, Sun C, Long S J, Zheng J Y, Wang P F, Jia Y Z, Zhang J S 2021 Corros. Sci. 192 109808
Google Scholar
[7] Tian Z, Peng J C, Lin X D, Hu Y Y, Yao M Y, Xie Y P, Liang X, Zhou B X 2024 Corros. Sci. 202 111937
Google Scholar
[8] Ferreirós P A, Polack E C S, Lanzani L A, Alonso P R, Quirós P D, Mieza J I, Rubiolo G H 2021 J. Nucl. Mater. 553 153039
Google Scholar
[9] Jiang G Y, Xu D H, Yang W P, Liu L, Zhi Y W, Yang J Q 2022 Prog. Nucl. Energy 154 104490
Google Scholar
[10] Yuan R, Xie Y P, Li T, Xu C H, Yao M Y, Xu J X, Guo H B, Zhou B X 2021 Acta Mater. 209 116804
Google Scholar
[11] Mcdeavitt S M, Billings G W, Indacochea J E 2002 J. Mater. Sci. 37 3765
Google Scholar
[12] Jerlerud P R, Toffolon M C, Joubert J M, Sundman B 2008 CALPHAD 32 593
Google Scholar
[13] Asle Zaeem M, El Kadiri H 2014 Comput. Mater. Sci. 89 122
Google Scholar
[14] 张更, 王巧, 沙立婷, 李亚捷, 王达, 施思齐 2020 物理学报 69 226401
Google Scholar
Zhang G, Wang Q, Sha L T, Li Y J, Wang D, Shi S Q 2020 Acta Phys. Sin 69 226401
Google Scholar
[15] Chen L Q, Zhao Y 2022 Prog. Mater. Sci. 124 100868
Google Scholar
[16] Chen L Q 2002 Annu. Rev. Mater. Res. 32 113
Google Scholar
[17] 刘续希, 柳文波, 李博岩, 贺新福, 杨朝曦, 恽迪 2022 金属学报 58 943
Liu X X, Liu W B, Li B Y, He X F, Yang Z X, Yun D 2022 Acta Metal. Sin 58 943
[18] Mai W, Soghrati S 2017 Corros. Sci. 125 87
Google Scholar
[19] Fang X R, Pan Z C, Ma R J, Chen A R 2023 Comuput. Methods Appl. Mech. Engrg. 414 116196
Google Scholar
[20] 冯力, 王智平, 路阳, 朱昌盛 2008 物理学报 57 1084
Google Scholar
Feng L, Wang Z P, Lu Y, Zhu C S 2008 Acta Phys. Sin 57 1084
Google Scholar
[21] Yang C, Huang H B, Liu W B, Wang J S, Wang J, Jafri H M, Liu Y, Han G M, Song H F, Chen L Q 2021 Adv. Theor. Simul. 4 2000162
Google Scholar
[22] Zhou F Y, Qiu K J, Bian D, Zheng Y F, Lin J P 2014 J. Mater. Sci. Technol. 30 299
Google Scholar
[23] Dinsdale A T 1991 CALPHAD 15 317
Google Scholar
[24] Aricó S F, Gribaudo L M, 1999 Scripta Mater. 41 159
Google Scholar
[25] Wang C, Zinkevich M, Aldinger F 2004 CALPHAD 28 281
Google Scholar
[26] Yin T, Lee J, Moosavi K E, Jung I H 2021 Ceram. Int. 47 29267
Google Scholar
[27] Isomäki I, Hämäläinen M, Gierlotka W, Onderka B, Fitzner K 2006 J. Alloys Compd. 422 173
Google Scholar
[28] Fang X R, Pan Z C, Chen A R, Tian H, Ma R J 2023 Eng. Fract. Mech. 281 109131
Google Scholar
[29] Allen S M, Cahn J W 1979 Acta Metall. 27 1085
Google Scholar
[30] Larché F C, Cahn J W 1985 Acta Metall. 33 331
Google Scholar
[31] Cahn J W 1961 Acta Metall. 9 795
Google Scholar
[32] Cox B, Pemsler J P 1968 J. Nucl. Mater. 28 73
Google Scholar
[33] Ritchie I G, Atrens A 1977 J. Nucl. Mater. 67 254
Google Scholar
[34] Keneshea F J, Douglass D L 1971 Oxid. Met. 3 1
Google Scholar
[35] 姜彦博, 柳文波, 孙志鹏, 喇永孝, 恽迪 2022 物理学报 71 026103
Google Scholar
Jiang Y B, Liu W B, Sun Z P, La Y X, Yun D 2022 Acta Phys. Sin 71 026103
Google Scholar
[36] Qiu K, Wang R, Peng C, Lu X, Wang N 2015 CALPHAD 48 175
Google Scholar
[37] Fisher E, Renken C 1964 Phys. Rev. 135 A482
Google Scholar
[38] Fogaing E Y, Lorgouilloux Y, Huger M, Gault C P 2006 J. Mater. Sci. 41 7663
Google Scholar
[39] Mirgorodsky A, Smirnov M, Quintard P 1997 Phys. Rev. B 55 19
Google Scholar
[40] 刘建章 2007 核结构材料 (北京: 化学工业出版社) 第89页
Liu J Z 2007 Nuclear Structure Material (Beijing: Chemical Industry Press) p89
[41] Fan D, Chen L Q 1997 Acta Mater. 45 611
Google Scholar
[42] Bi Y B, Zhang X L, Lu L, Xu Z F, Xie Z G, Chen B B, Liang Z X, Sun Z Q, Luo Z 2023 J. Mater. Res. Technol. 26 5888
Google Scholar
[43] Gwinner B, Bataillon C, Chelagemdib R, Gruet N, Lorentz V, Puga B 2023 Electrochim. Acta 470 143334
Google Scholar
[44] Zhang M H, Liu S D, Jiang J Y, Wei W C 2023 T Nonferr Metal. Soc. 33 1963
Google Scholar
[45] Liu J L, Yu H B, Karamched P, Hu J, He G Z, Goran D, Hughes G M, Wilkinson A J, Sergio L P, Grovenor C R 2019 Acta Mater. 179 328
Google Scholar
-
表 1 Zr-Sn合金自由能参数
Table 1. Free energy parameters of Zr-Sn alloys.
Parameter Value Reference EZr/(J·m–3) −7827.6+125.65T−24.1618TlnT−4.38e−3T3+34971T−1 [23] ESn/(J·m–3) 2524.7+4.00T−8.26T\lnT−16.81e−3T2+2.62e−6T3−1.08e6T−1 [23] EO/(J·m–3) −3480.9−25.50T−11.13TlnT−5.10e−3T2+0.66e−6T3−38365T−1 [23] AZrSn −148022.5+19.41T+(173681.9−22T)(cZr−cSn)+104271.96(cZr−cSn)2 [24] AZrO −37876.66+17.2915T−4471.4(cZr−cO) [25] ASnO 140878 – 23.9326T [26, 35] g02 −1117869+420.3T−69.6TlnT−0.003766T2+702910T−1+4.59e−21T7 [25] g03 0 k2/(J·m–3) 2.5 × 105 [13] k3/(J·m–3) 2.5 × 105 [13] Temperature/K D1/(m2·s–1) D2/(m2·s–1) D3/(m2·s–1) DGB/(m2·s–1) 633 3.99×10–17 1.61×10–17 2.98×10–17 1.46×10–16 653 1.26×10–16 5.18×10–17 8.77×10–17 8.11×10–16 673 3.09×10–16 1.11×10–16 2.25×10–16 3.76×10–15 693 1.04×10–15 3.69×10–16 7.64×10–16 2.82×10–14 -
[1] Kai J J, Huang W I, Chou H Y 1990 J. Nucl. Mater. 170 193
Google Scholar
[2] Zhao W Z, Quan Q W, Zhang S M, Yang X Y, Wen H Y, Wu Y C, Liu X B, Cao X Z, Wang B Y 2023 Radiat. Phys. Chem. 209 110986
Google Scholar
[3] Jones R H, Simonen E P 1994 Mater. Sci. Eng. 176 211
Google Scholar
[4] Hu J, Liu J L, Lozano P S, Grovenor C R M, Christensen M, Wolf W, Wimmer E, Mader E V 2019 Acta Mater. 180 105
Google Scholar
[5] Cann C D, So C B, Styles R C, Coleman C E 1993 J. Nucl. Mater. 205 267
Google Scholar
[6] Wei T G, Dai X, Long C S, Sun C, Long S J, Zheng J Y, Wang P F, Jia Y Z, Zhang J S 2021 Corros. Sci. 192 109808
Google Scholar
[7] Tian Z, Peng J C, Lin X D, Hu Y Y, Yao M Y, Xie Y P, Liang X, Zhou B X 2024 Corros. Sci. 202 111937
Google Scholar
[8] Ferreirós P A, Polack E C S, Lanzani L A, Alonso P R, Quirós P D, Mieza J I, Rubiolo G H 2021 J. Nucl. Mater. 553 153039
Google Scholar
[9] Jiang G Y, Xu D H, Yang W P, Liu L, Zhi Y W, Yang J Q 2022 Prog. Nucl. Energy 154 104490
Google Scholar
[10] Yuan R, Xie Y P, Li T, Xu C H, Yao M Y, Xu J X, Guo H B, Zhou B X 2021 Acta Mater. 209 116804
Google Scholar
[11] Mcdeavitt S M, Billings G W, Indacochea J E 2002 J. Mater. Sci. 37 3765
Google Scholar
[12] Jerlerud P R, Toffolon M C, Joubert J M, Sundman B 2008 CALPHAD 32 593
Google Scholar
[13] Asle Zaeem M, El Kadiri H 2014 Comput. Mater. Sci. 89 122
Google Scholar
[14] 张更, 王巧, 沙立婷, 李亚捷, 王达, 施思齐 2020 物理学报 69 226401
Google Scholar
Zhang G, Wang Q, Sha L T, Li Y J, Wang D, Shi S Q 2020 Acta Phys. Sin 69 226401
Google Scholar
[15] Chen L Q, Zhao Y 2022 Prog. Mater. Sci. 124 100868
Google Scholar
[16] Chen L Q 2002 Annu. Rev. Mater. Res. 32 113
Google Scholar
[17] 刘续希, 柳文波, 李博岩, 贺新福, 杨朝曦, 恽迪 2022 金属学报 58 943
Liu X X, Liu W B, Li B Y, He X F, Yang Z X, Yun D 2022 Acta Metal. Sin 58 943
[18] Mai W, Soghrati S 2017 Corros. Sci. 125 87
Google Scholar
[19] Fang X R, Pan Z C, Ma R J, Chen A R 2023 Comuput. Methods Appl. Mech. Engrg. 414 116196
Google Scholar
[20] 冯力, 王智平, 路阳, 朱昌盛 2008 物理学报 57 1084
Google Scholar
Feng L, Wang Z P, Lu Y, Zhu C S 2008 Acta Phys. Sin 57 1084
Google Scholar
[21] Yang C, Huang H B, Liu W B, Wang J S, Wang J, Jafri H M, Liu Y, Han G M, Song H F, Chen L Q 2021 Adv. Theor. Simul. 4 2000162
Google Scholar
[22] Zhou F Y, Qiu K J, Bian D, Zheng Y F, Lin J P 2014 J. Mater. Sci. Technol. 30 299
Google Scholar
[23] Dinsdale A T 1991 CALPHAD 15 317
Google Scholar
[24] Aricó S F, Gribaudo L M, 1999 Scripta Mater. 41 159
Google Scholar
[25] Wang C, Zinkevich M, Aldinger F 2004 CALPHAD 28 281
Google Scholar
[26] Yin T, Lee J, Moosavi K E, Jung I H 2021 Ceram. Int. 47 29267
Google Scholar
[27] Isomäki I, Hämäläinen M, Gierlotka W, Onderka B, Fitzner K 2006 J. Alloys Compd. 422 173
Google Scholar
[28] Fang X R, Pan Z C, Chen A R, Tian H, Ma R J 2023 Eng. Fract. Mech. 281 109131
Google Scholar
[29] Allen S M, Cahn J W 1979 Acta Metall. 27 1085
Google Scholar
[30] Larché F C, Cahn J W 1985 Acta Metall. 33 331
Google Scholar
[31] Cahn J W 1961 Acta Metall. 9 795
Google Scholar
[32] Cox B, Pemsler J P 1968 J. Nucl. Mater. 28 73
Google Scholar
[33] Ritchie I G, Atrens A 1977 J. Nucl. Mater. 67 254
Google Scholar
[34] Keneshea F J, Douglass D L 1971 Oxid. Met. 3 1
Google Scholar
[35] 姜彦博, 柳文波, 孙志鹏, 喇永孝, 恽迪 2022 物理学报 71 026103
Google Scholar
Jiang Y B, Liu W B, Sun Z P, La Y X, Yun D 2022 Acta Phys. Sin 71 026103
Google Scholar
[36] Qiu K, Wang R, Peng C, Lu X, Wang N 2015 CALPHAD 48 175
Google Scholar
[37] Fisher E, Renken C 1964 Phys. Rev. 135 A482
Google Scholar
[38] Fogaing E Y, Lorgouilloux Y, Huger M, Gault C P 2006 J. Mater. Sci. 41 7663
Google Scholar
[39] Mirgorodsky A, Smirnov M, Quintard P 1997 Phys. Rev. B 55 19
Google Scholar
[40] 刘建章 2007 核结构材料 (北京: 化学工业出版社) 第89页
Liu J Z 2007 Nuclear Structure Material (Beijing: Chemical Industry Press) p89
[41] Fan D, Chen L Q 1997 Acta Mater. 45 611
Google Scholar
[42] Bi Y B, Zhang X L, Lu L, Xu Z F, Xie Z G, Chen B B, Liang Z X, Sun Z Q, Luo Z 2023 J. Mater. Res. Technol. 26 5888
Google Scholar
[43] Gwinner B, Bataillon C, Chelagemdib R, Gruet N, Lorentz V, Puga B 2023 Electrochim. Acta 470 143334
Google Scholar
[44] Zhang M H, Liu S D, Jiang J Y, Wei W C 2023 T Nonferr Metal. Soc. 33 1963
Google Scholar
[45] Liu J L, Yu H B, Karamched P, Hu J, He G Z, Goran D, Hughes G M, Wilkinson A J, Sergio L P, Grovenor C R 2019 Acta Mater. 179 328
Google Scholar
期刊类型引用(1)
1. 廖宇轩,申文龙,吴学志,喇永孝,柳文波. 陶瓷型复合燃料烧结过程的相场模拟研究. 物理学报. 2024(21): 7-17 . 百度学术
其他类型引用(0)
计量
- 文章访问数: 2774
- PDF下载量: 66
- 被引次数: 1