第十三章钢铁在防护层剥离情况下的缝隙腐蚀F.M.SONG,N.SRIDHAR.doc

第十三章 钢铁在防护层剥离情况下的缝隙腐蚀 F. M. SONG, N. SRIDHAR Mechanical and Materials Engineering, Southwest Research Institute, 6220, Culebra Road, San Antonio, Texas 78238, USA E-mail: fsong@ 缝隙腐蚀是降低地下管道寿命的主要腐蚀形式。腐蚀机理目前尚不明确，有待进一步的研究。这里通过在二维缝隙结构上耦合两个拉普拉斯方程来模拟这一问题，用到了有限元计算软件FEMLAB，问题得到了很好的求解。使用导出、积分耦合变量和试算方法来寻求没有阴极保护（CP）情况下的电势。这里计算了缝隙中溶解的氧气浓度，钢铁电势和电流。 1．背景 1.1 研究意义 粗略估计，2004年由于腐蚀给美国带来2760亿美元的损失，其中有70亿美元来自油气管道的腐蚀[1]。埋藏在地下的管道通常通过聚合物涂层和阴极保护（CP）来避免腐蚀。阴极保护是通过与地下水接触来降低管道金属化学电势的电化学防腐措施。即使管道保护起来，存在的腐蚀仍然很严重。尤其是当聚合物涂层和管道钢铁剥离时更为严重，这时地下水会在缝隙中聚集，导致阴极保护不在管用。本章提供了一个例子说明如何通过数学建模研究腐蚀机理。关于建模的更多信息可以参考文献[2]-[5]。 1.2 模型的独特性和挑战性 求解该问题的主要难点包括钢铁表面复杂的非线性边界条件和模拟真实管道涂层腐蚀带来的复杂几何结构。在求解这个问题中用到了导出、积分耦合变量和其他一些技术。为了高效计算，在不同边界条件处使用了不同的网格尺寸。 1.3 模型描述 图1中的模型几何结构示意性的描述了缝隙腐蚀现象。它是左右对称的。B1边界左侧的另一半没有显示。漏涂缺口和剥离涂层（B4）将整个几何结构分为两部分，上面是土壤溶液，下面是发生钢铁腐蚀的缝隙。假设边界（B5）为测定管道钢铁电势的参考电极。沿着该边界外侧，参考电极的电势和氧气浓度都设单位1。模型的坐标系为x和y。几何结构尺寸见表1。 图1 管道缝隙腐蚀结构示意图。B1是对称边界；B2是钢铁表面，有腐蚀流量；B3和B4是不渗水的涂层，没有腐蚀流量；B5是测量边界，具有单位电势和单位氧气浓度 表1 几何尺寸 变量 值 漏涂缺口半长 0.001m 涂层和钢铁表面见的间隙 0.0005m 缝隙长度 0.1m 上部区域高度 0.0995m 上部区域宽度 0.01m 溶解在土壤溶液中的氧气可以通过涂层缺口或直接穿过涂层进入缝隙中。由于氧气穿过涂层的渗透率很小，只考虑穿过涂层缺口的氧气。阴极保护电流只能通过涂层缺口进入缝隙。 在稳态情况下，模型的控制方程是以氧气浓度（cO2）表示的溶解氧的质量守恒方程： (1) 同时，假设缝隙内外侧为单位离子浓度，电荷守恒就可以用任意位置处相对于参考电极的金属电势来表示： (2) 前面已经介绍过B5的边界条件。由于对称性，B1边界的氧气浓度和电势都为零。B3和B4边界的通量为零，因为涂层隔绝了电流和氧气扩散。B2通量边界定义如下。 在管道钢铁表面，阳极和阴极反应同时存在。阳极的半电池反应是离子氧化或腐蚀反应： (3) 腐蚀电流密度用Tafel方程表示： (4) 这里是钢铁表面的。和分别表示交换电流密度和相同情况下方程(3)的平衡电势。是离子氧化的Tafel斜率。 对于非酸性土壤溶液，最主要的阴极半电池反应是水的还原反应： (5) 和氧气的还原： (6) 它们的Tafel方程分别为： (7) (8) 这里和是交换电流密度和在相同参考情况下方程(5)的平衡电势。是水还原的Tafel斜率。和是当氧气浓度为时，交换电流密度和在相同参考情况下方程(5)的平衡电势。是钢铁表面的氧浓度。是氧气还原的Tafel斜率。注意到方程(3)，(5)和(6)的参考情况是彼此独立。 在边界B2，氧气扩散通量和阴极还原速率相等： (9) 这里是地下水中的氧气扩散率。是基元反应(6)中转移的电子数，F是法拉第常数。根据欧姆定律可以得到以下方程： (10) 这里是地下水的电阻率，。 2．FEMLAB建模过程和求解 2.1 变量和常数 表2给出了该模型的应用模式和几何结构绘制说明。表3给出了模型所必需的常数项。在Option菜单中按照以下步骤建立常数、表达式和耦合变量。 表2模型建立——应用模式和绘制几何模型 运行FEMLAB 在模型导航栏中选择空间维数为二维 在FEMLAB的应用模式中，选择PDE General Form，其中从PDE Modes