AAM Advances in Applied Mathematics 2324-7991 Scientific Research Publishing 10.12677/AAM.2021.108297 AAM-44736 AAM20210800000_22396806.pdf 数学与物理 混凝土中变扩散系数的时间分数阶氯离子扩散模型研究 Research on the Time-Fractional Diffusion Model of Chloride Ion with Variable Diffusion Coefficient in Concrete 王 一晗 2 1 白 羽 2 1 3 张 艳 2 1 3 巩 璇 4 1 北京建筑大学理学院，北京 北京建筑大学建筑结构与环境修复功能材料北京市重点实验室，北京；北京建筑大学土木与交通工程学院，北京 null 北京建筑大学理学院，北京；北京建筑大学建筑结构与环境修复功能材料北京市重点实验室，北京；北京建筑大学大数据建模理论与技术研究所，北京 03 08 2021 10 08 2853 2861 © Copyright 2014 by authors and Scientific Research Publishing Inc. 2014 This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/

研究氯离子在混凝土中的扩散机理可以有针对性地对混凝土进行修护，从而有效地延长混凝土的耐久性。对于氯离子在混凝土中扩散的一维问题，首先基于Caputo分数阶导数建立了变扩散系数的时间分数阶扩散模型，结合L1算法得到了数值求解该模型的有限差分格式。随后，针对一个实例，利用已有实验数据和最小二乘原理优化了模型中不同暴露时间下分数阶导数的阶数；通过数据拟合得到了时间分数阶导数的阶数与暴露时间之间的关系，进而建立了混凝土中变扩散系数的时间变分数阶的氯离子扩散模型。 Research on the diffusion mechanism of Chloride ion in concrete can be applied to repair concrete, which can effectively prolong the durability of concrete. This paper is focusing on one-dimensional diffusion of Chloride ion in concrete. Firstly, the time-fractional diffusion model of Chloride ion with variable diffusion coefficient in concrete is established based on Caputo fractional derivative, which is solved by the finite difference scheme combining with L1 algorithm. Next, the fractional derivative parameter on different exposure time is optimized by using the existing experimental data and the least square method in the model. At last, the relationship between the optimal fractional derivative parameter and the exposure time is fitted, which results in the time variable fractional diffusion model of Chloride ion with variable diffusion coefficient in concrete.

 混凝土，氯离子扩散，分数阶导数，有限差分方法, Concrete Chloride Ion Diffusion Fractional Derivative Finite Difference Method 摘要

研究氯离子在混凝土中的扩散机理可以有针对性地对混凝土进行修护，从而有效地延长混凝土的耐久性。对于氯离子在混凝土中扩散的一维问题，首先基于Caputo分数阶导数建立了变扩散系数的时间分数阶扩散模型，结合L1算法得到了数值求解该模型的有限差分格式。随后，针对一个实例，利用已有实验数据和最小二乘原理优化了模型中不同暴露时间下分数阶导数的阶数；通过数据拟合得到了时间分数阶导数的阶数与暴露时间之间的关系，进而建立了混凝土中变扩散系数的时间变分数阶的氯离子扩散模型。

 关键词

混凝土，氯离子扩散，分数阶导数，有限差分方法

 Research on the Time-Fractional Diffusion Model of Chloride Ion with Variable Diffusion Coefficient in Concrete<sup> </sup>

Yihan Wang1, Yu Bai1,2,3*, Yan Zhang1,2,3, Xuan Gong2,4

1School of Science, Beijing University of Civil Engineering and Architecture, Beijing

2Beijing Key Laboratory of Functional Materials for Building Structure and Environment Remediation, Beijing University of Civil Engineering and Architecture, Beijing

3Institute of Big Data Modelling and Technology, Beijing University of Civil Engineering and Architecture, Beijing

4School of Civil and Transportation Engineering, Beijing University of Civil Engineering and Architecture, Beijing

Received: Jul. 23rd, 2021; accepted: Aug. 15th, 2021; published: Aug. 24th, 2021

 ABSTRACT

Research on the diffusion mechanism of Chloride ion in concrete can be applied to repair concrete, which can effectively prolong the durability of concrete. This paper is focusing on one-dimensional diffusion of Chloride ion in concrete. Firstly, the time-fractional diffusion model of Chloride ion with variable diffusion coefficient in concrete is established based on Caputo fractional derivative, which is solved by the finite difference scheme combining with L1 algorithm. Next, the fractional derivative parameter on different exposure time is optimized by using the existing experimental data and the least square method in the model. At last, the relationship between the optimal fractional derivative parameter and the exposure time is fitted, which results in the time variable fractional diffusion model of Chloride ion with variable diffusion coefficient in concrete.

Keywords:Concrete, Chloride Ion Diffusion, Fractional Derivative, Finite Difference Method

Copyright © 2021 by author(s) and beplay安卓登录

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

 1. 引言

在现代建筑物的建造中，混凝土由于价格低、建筑综合效益好，且具有良好的耐冲击性、可塑性和耐火性等特点，成为了最普遍的施工材料，从而其耐久性和服役寿命成为备受关注的问题。针对目前的土木工程，混凝土结构的耐久性不容乐观，建筑物承载力不足、开裂等问题屡见不鲜 [ 1 ]。混凝土中的氯离子主要通过扩散作用、渗透作用、毛细管作用和电化学作用等侵入到混凝土和钢筋表面。通常，氯离子的侵蚀是几种侵入方式的组合。但是，许多情况下认为扩散是最主要的侵入方式 [ 2 ]。

许多研究人员基于整数阶导数的Fick第二扩散定律来描述氯离子的扩散行为，Buenfeld等 [ 3 ] 假定混凝土中的氯离子传输过程是一个非稳态的扩散过程，并提出了此过程的假设条件；Mangat等 [ 4 ] 将时间因素考虑到氯离子的扩散系数中，并给出扩散系数与时间的函数表达式：

D = D i t − m (1)

式中Di为t等于1个时间单位时的扩散系数， t ≥ 0 ；m为时间依赖性常数。同年Thomas等 [ 5 ] 修正了该模型中Di不可测的缺点，提出氯离子的扩散系数是时变函数的观点，并给出下式：

D t = D 0 ( t 0 t ) m (2)

式中D0和Dt分别为扩散时间为t0和t时混凝土的氯离子扩散系数， t ≥ 0 。Thomas给出的(2)式说明，t0时测定混凝土的氯离子扩散系数为D0。

近年来，研究人员发现分数阶微积分可以更简洁准确地描述具有历史记忆和空间全域相关性等复杂力学与物理现象的过程，因而基于分数阶导数的扩散模型开始受到关注。宦俊蒙 [ 6 ] 使用拉普拉斯变换法及Crump反演方法求解分数阶扩散方程的半解析解。Chen等 [ 7 ] 考虑了氯离子和混凝土之间会发生相互作用，从而导致一定比例的氯离子与混凝土结合，所以使用了朗缪尔等温线描述氯离子浓度。Wei等 [ 8 ] 建立了氯离子整数阶扩散模型和Caputo型氯离子时间分数阶扩散模型，并通过均方位移的方法近似得到扩散系数Dα。

本文将基于Caputo型分数阶导数，考虑形如(2)式的变扩散系数，对混凝土中氯离子扩散时间分数阶的一维模型进行研究。

 2. 变扩散系数的时间分数阶扩散模型 2.1. 模型建立

基于Fick第二定律和Caputo型分数阶导数 [ 9 ]，考虑氯离子在混凝土中的一维扩散问题，建立变扩散系数的时间分数阶扩散模型，描述为：

{ ∂ α C ( x , t ) ∂ t α = D 0 ( t 0 t ) m ∂ 2 C ( x , t ) ∂ x 2 , 0 < x < + ∞ , t > 0 C ( x , 0 ) = C 0 , 0 ≤ x ≤ + ∞ C ( 0 , t ) = C s , t > 0 C ( + ∞ , t ) = C 0 , t > 0 (3)

式中∂α/∂tα为Caputo型分数阶算子；x为侵蚀深度， x ∈ [ 0 , + ∞ ) ；t为混凝土结构中氯离子暴露于环境的时间， t ≥ 0 ；C(x, t)表示t时刻x深处下的氯离子浓度；D0和Dt分别为扩散时间为t0和t时的扩散系数；α为分数阶导数的阶数， 0 < α ≤ 1 ；Cs为混凝土表面氯离子浓度；C0为氯离子内部初始浓度；C0和Cs是常数；m为时间依赖性常数。

 2.2. 模型求解的有限差分格式

这里结合L1算法 [ 9 ]：

D t α f ( t k ) = Δ t − α Γ ( 2 − α ) ∑ s = 0 k − 1 α s [ f ( t k − s ) − f ( t k − s − 1 ) ] + O ( Δ t 2 − α ) = Δ t − α Γ ( 2 − α ) [ f ( t k ) − α k − 1 f ( t 0 ) − ∑ s = 1 k − 1 ( α s − 1 − α s ) f ( t k − s ) ] + O ( Δ t 2 − α ) (4)

其中 α s = ( s + 1 ) 1 − α − s 1 − α , s = 0 , 1 , 2 , ⋯ ，建立求解模型(3)的有限差分格式。

首先，选取足够大的空间长度L和暴露时间T。定义 x = i h   ( i = 0 , 1 , ⋯ , M ) ， t = k τ   ( k = 0 , 1 , ⋯ , N ) ，其中h = L/M为空间步长，τ=T/N为时间步长。这样在划分求解区域后， C i k 表示氯离子浓度在点(xi, tk)处的近似值。

接下来，将整数阶导数项在 ( x = x i , t = t k − 1 2 ) 处离散：

C | ( x = x i , t = t k − 1 2 ) = C ( x i , t k ) + C ( x i , t k − 1 ) 2 + O ( Δ t 2 ) (5)

∂ C ∂ x | ( x = x i , t = t k − 1 2 ) = C ( x i , t k ) − C ( x i − 1 , t k ) + C ( x i , t k − 1 ) − C ( x i − 1 , t k − 1 ) 2 Δ x + O ( Δ x ) ， (6)

∂ 2 C ∂ x 2 | ( x = x i , t = t k − 1 2 ) = C ( x i + 1 , t k ) − 2 C ( x i , t k ) + C ( x i − 1 , t k ) + C ( x i + 1 , t k − 1 ) − 2 C ( x i , t k − 1 ) + C ( x i − 1 , t k − 1 ) 2 Δ x 2 + O ( Δ x 2 + Δ t 2 ) (7)

然后，将(4)~(7)代入扩散模型(3)，并用 C i k 近似C(xi, tk)，得到有限差分格式：

r 1 [ C i k + C i k − 1 2 − α k − 1 C i 1 + C i 0 2 − ∑ s = 1 k − 1 ( α s − 1 − α s ) C i k − s + C i k − s − 1 2 ] = r 2 C i + 1 k − 2 C i k + C i − 1 k + C i + 1 k − 1 − 2 C i k − 1 + C i − 1 k − 1 2 ( i = 1 , 2 , ⋯ , M − 1 ; k = 1 , 2 , ⋯ , N ) (8)

对应的初始条件和边界条件为：

C i 0 = C 0 ( i = 0 , 1 , ⋯ , M ) (9)

C 0 k = C s , C M k = C 0 ( k = 1 , 2 , ⋯ , N ) (10)

其中 r 1 = h 2 Γ ( 2 − α ) ， r 2 = D 0 [ t 0 ( k − 1 ) τ ] m τ α 。代入参数值，通过求解有限差分格式(8)~(10)，就可以得到变扩散系数的时间分数阶扩散模型的数值解。

 3. 实例

Thomas等 [ 5 ] 进行了一系列氯离子扩散实验，所得的PC (混凝土预制件)中氯离子含量观测值如表1所示，各常数参数值如表2所示。

Chloride concentration of Precast Concrete on different exposure time [5
种类 深度(m) 氯离子浓度(%)
6个月 1年 2年 3年 6年 8年
PC 0.005 0.267 0.493 0.313 0.370 0.257 0.288
0.015 0.140 0.193 0.273 0.273 0.237 0.264
0.025 0.050 0.043 0.143 0.190 0.197 0.203
0.035 0.010 0.017 0.093 0.117 0.153 0.183
0.045 0.000 0.000 0.053 0.073 0.133 0.159

表1. PC不同暴露时间下的氯离子浓度值 [ 5 ]

Initial concentration, surface concentration and values of parameters [<xref ref-type="bibr" rid="hanspub.44736-ref5">5</xref>] [8
种类 C0(%) Cs(%) D0(m2/s) m
PC 0 0.35 8 × 10−12 0.1

表2. 初始浓度、表面浓度和参数值 [ 5 ] [ 8 ]

3.1. 固定分数阶导数的阶数

利用有限差分格式(8)~(10)，将表1和表2所提供的数据和参数值代入计算，就可以得到变扩散系数的时间分数阶扩散模型(3)在不同暴露时间下的氯离子浓度曲线。当取分数阶导数的阶数α = 0.9时，PC混凝土块在6个月、1~3年、6年和8年时的氯离子浓度曲线如图1(a)~(f)所示。图中绿色的数据点为实验数据 [ 5 ]，蓝色曲线为扩散模型(3)的数值解，数值解在一定程度上与实验数据一致，验证了模型的有效性。

 3.2. 优化分数阶导数的阶数

为了更准确地反映氯离子的扩散行为，采用最小二乘法对分数阶导数的阶数α进行优化，也就是寻找α使得扩散模型的数值解与实验数据的误差的平方和Q达到最小值，即

min α Q = min α ∑ i = 1 5 ( C i | 扩 散 模 型 数 值 解 − C i | 实 验 数 据 ) 2 (11)

在不同暴露时间下，以精度0.0001的要求求解问题(11)，得到的结果如表3所示。可见，随着暴露时间的变化，最优的分数阶导数的阶数也在改变。利用表3中的结果，重新求解扩散模型(3)，得到了不同暴露时间下使用优化的分数阶导数的阶数α的氯离子浓度曲线，如图2(a)~(f)所示。

图1. 不同暴露时间下α = 0.9时氯离子浓度曲线

图2. 不同暴露时间下优化的分数阶导数的阶数α的氯离子浓度曲线

Optimal fractional derivative parameter α and the corresponding error square sum Q on different exposure tim
暴露时间 6个月 1年 2年 3年 6年 8年
最优阶数α 0.975 0.956 0.950 0.941 0.875 0.890
误差平方和Q 9.3430 × 10−5 0.0466 0.0027 0.0047 0.0043 0.0019

表3. 不同暴露时间下最优阶数α和对应的误差平方和Q

 3.3. 变扩散系数的时间变分数阶扩散模型

根据表3的结果，去除暴露时间为1年时的异常点后，采用三阶多项式拟合最优分数阶导数的阶数α与暴露时间t之间的关系，得到：

α ( t ) = 0.0008 t 3 − 0.0082 t 2 + 0.0048 t + 0.9731 (12)

拟合效果如图3所示。

图3. 最优分数阶导数的阶数α与暴露时间t之间的关系

用(12)中的函数α(t)替换扩散模型(3)式中的分数阶导数的阶数α，就可以得到PC混凝土块的变扩散系数的氯离子时间变分数阶扩散模型：

{ ∂ α ( t ) C ( x , t ) ∂ t α ( t ) = D 0 ( t 0 t ) m ∂ 2 C ( x , t ) ∂ x 2 , 0 < x < + ∞ , t > 0 C ( x , 0 ) = C 0 ⋅ δ ( x ) , 0 ≤ x ≤ + ∞ C ( 0 , t ) = C s , t > 0 C ( + ∞ , t ) = C 0 , t > 0 (13)

式中 0 < α ( t ) ≤ 1 满足(12)。

 3.4. 氯离子扩散浓度预测

当暴露时间分别为7年和9年时，对于扩散模型(3)选取α = 0.9计算氯离子浓度分布，如图4(a)和图4(b)中蓝色曲线所示。此时扩散模型(13)中的α(7) = 0.8748和α(9) = 0.9304，计算氯离子浓度如图4(a)和图4(b)中红色曲线所示。图4(a)中黑色星号为6年实验数据，图4(b)中绿色圆圈为8年实验数据。从图中不难看出，扩散模型(3)和扩散模型(13)都可以在一定程度上预测氯离子浓度的扩散情况，相比较而言扩散模型(13)更准确一些。

图4. 暴露时间分别为7年和9年时氯离子浓度的预测值

 4. 结论

· 变扩散系数的时间分数阶氯离子扩散模型可以有效地预测氯离子在混凝土中的扩散情况。

· 结合L1算法建立的有限差分格式实现了分数阶扩散模型的数值求解。

· 利用已有实验数据和最小二乘法优化了分数阶导数的阶数，进而建立了混凝土中变扩散系数的时间变分数阶氯离子扩散模型。

 基金项目

北京市属高校基本科研业务费专项资金项目(X20142, X21027, X21030, X21031)，北京建筑大学教育科学研究项目(Y2015)，北京建筑大学研究生教学质量提升项目(J2020004)。

 文章引用

王一晗,白 羽,张 艳,巩 璇. 混凝土中变扩散系数的时间分数阶氯离子扩散模型研究Research on the Time-Fractional Diffusion Model of Chloride Ion with Variable Diffusion Coefficient in Concrete[J]. 应用数学进展, 2021, 10(08): 2853-2861. https://doi.org/10.12677/AAM.2021.108297

 参考文献 References 夏瑜. 混凝土结构耐久性及加固技术探究[J]. 四川水泥, 2021(5): 11-12. Saetta, A.V., Scotta, R.V. and Vitaliani, R.V. (1993) Analysis of Chloride Diffusion into Partially Saturated Concrete. ACI Materials Journal, 90, 441-451.
https://doi.org/10.14359/3874 Buenfeld, N.R., Glass, G.K., Hassanein, A.M., et al. (1998) Chloride Transport in Concrete Subjected to Electric Field. Journal of Materials in Civil Engineering, 10, 220-228.
https://doi.org/10.1061/(ASCE)0899-1561(1998)10:4(220) Mangat, P.S. and Limbachiya, M.C. (1999) Effect of Initial Curing on Chloride Diffusion in Concrete Repair Materials. Cement and Concrete Research, 29, 1475-1485.
https://doi.org/10.1016/S0008-8846(99)00130-1 Thomas, M. and Bamforth, P.B. (1999) Modelling Chloride Diffusion in Concrete: Effect of Fly Ash and Slag. Cement and Concrete Research, 29, 487-495.
https://doi.org/10.1016/S0008-8846(98)00192-6 宦俊蒙. 混凝土中氯离子传输机理的理论研究[D] : [硕士学位论文]. 舟山: 浙江海洋大学, 2020. Chen, W., Zhang, J. and Zhang, J. (2013) A Variable-Order Time-Fractional Derivative Model for Chloride Ions Sub-Diffusion in Concrete Structures. Fractional Calculus & Applied Analysis, 16, 76-92.
https://doi.org/10.2478/s13540-013-0006-y Wei, S., Chen, W. and Zhang, J. (2015) Time-Fractional Derivative Model for Chloride Ions Sub-Diffusion in Reinforced Concrete. European Journal of Environmental and Civil Engineering, 21, 1-13.
https://doi.org/10.1080/19648189.2015.1116467 刘发旺, 庄平辉, 刘青霞. 分数阶偏微分方程数值方法及其应用[M]. 北京: 科学出版社, 2015: 107-144.
Baidu
map