Next Article in Journal
Ice as a Green-Structure-Directing Agent in the Synthesis of Macroporous MWCNTs and Chondroitin Sulphate Composites
Next Article in Special Issue
The Results of 45 Years of Atmospheric Corrosion Study in the Czech Republic
Previous Article in Journal
The Fabrication of Porous Si with Interconnected Micro-Sized Dendrites and Tunable Morphology through the Dealloying of a Laser Remelted Al–Si Alloy
Previous Article in Special Issue
Analysis of Historic Copper Patinas. Influence of Inclusions on Patina Uniformity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Corrosion Prediction with Parallel Finite Element Modeling for Coupled Hygro-Chemo Transport into Concrete under Chloride-Rich Environment

1
R&D Division, Hyundai E&C, Yongin-si, Gyeonggi-do 16891, Korea
2
Computer Science, University of Colorado Boulder, Boulder, CO 80309, USA
3
Civil, Environmental, and Architectural Engineering, University of Colorado Boulder, Boulder, CO 80309, USA
*
Author to whom correspondence should be addressed.
Materials 2017, 10(4), 350; https://doi.org/10.3390/ma10040350
Submission received: 26 February 2017 / Revised: 19 March 2017 / Accepted: 23 March 2017 / Published: 28 March 2017
(This article belongs to the Special Issue Fundamental and Research Frontier of Atmospheric Corrosion)

Abstract

:
The prediction of the chloride-induced corrosion is very important because of the durable life of concrete structure. To simulate more realistic durability performance of concrete structures, complex scientific methods and more accurate material models are needed. In order to predict the robust results of corrosion initiation time and to describe the thin layer from concrete surface to reinforcement, a large number of fine meshes are also used. The purpose of this study is to suggest more realistic physical model regarding coupled hygro-chemo transport and to implement the model with parallel finite element algorithm. Furthermore, microclimate model with environmental humidity and seasonal temperature is adopted. As a result, the prediction model of chloride diffusion under unsaturated condition was developed with parallel algorithms and was applied to the existing bridge to validate the model with multi-boundary condition. As the number of processors increased, the computational time decreased until the number of processors became optimized. Then, the computational time increased because the communication time between the processors increased. The framework of present model can be extended to simulate the multi-species de-icing salts ingress into non-saturated concrete structures in future work.

1. Introduction

One of the main long-term durability problems of reinforced concrete structure is the corrosion of reinforcing bars (rebar) in concrete. The corrosion can be resulted from several necessary conditions, and one of them is a high concentration of chloride ions near rebars, where the chloride ions come from deicing chemicals used in the winter on roadways and parking structures or seawater for the off-shore structures. Once the chloride content on the surface of steel reinforcement reaches a threshold value and the moisture and oxygen are sufficiently provided, the corrosion of steel bar is initiated [1,2]. The corrosion initiation period is the time during which substances such as water, chloride ions, oxygen and carbon dioxide penetrate through the concrete cover. The length of period depends on the resistance of concrete to the transport processes and the severity of the environmental conditions they are exposed to [3,4,5,6]. In practical engineering, the chloride contaminated RC structures have many internal cracks due to rebar corrosion and the prediction of crack width and propagation is very important for safety and serviceability [7,8]. Recently, the severity of corrosion damage was demonstrated through the corrosion-induced modeling with a statistical approach on bridge structures [9,10,11]. In concrete bridges, there are various concentrations of moisture and chemicals on the top surface of bridge decks due to the scatter of de-icing chemicals when they are applied. Especially, moisture and chloride concentrations have spatial distributions and seasonal variations. The service life of concrete is affected by a combination of material properties and microclimate. The microclimate is a term for the climatic conditions at the concrete surface or very close to the surface. This condition on the surface has a more decisive effect on the conditions inside the concrete than most other parameters [12]. The model of microclimate such as environmental humidity and temperature was proposed by Bazant et al. [13].
To simulate more realistically the durability performance of reinforced concrete structures, sophisticated numerical methods for fully coupled moisture and chloride transport processes and reliable material models for the transport parameters are needed. In addition, in order to predict accurately the moisture and chloride distributions in a large reinforced concrete structure, a large number of fine finite element meshes are needed. For the use of the large meshes in numerical analysis, in this study, a parallel algorithm will be employed to a parallel processing system with up to 2000 processors. To consider more realistic boundary conditions on concrete structures, humidity model applied with random stochastic process will be demonstrated and approximately 1.5 million nodes and three million elements will be used. Furthermore, multi-boundary conditions to describe the actual chloride concentrations and humidity distribution will be applied instead of only constant boundary conditions.

2. Basic Diffusion Formulation of Unsaturated Concrete

2.1. Governing Equation

The two fully-coupled partial differential equations governing the coupled chloride and moisture diffusion through non-saturated concrete are simply derived by employing the mass balance equations and Fick’s law [14,15].
First, the flux of chloride ions (J) through a unit area of porous media depends on the gradient of chloride ions as well as the gradient of moisture as. Thus, Fick’s law is modified in this study to include the coupling effects between moisture and chloride transfer. The governing equations are shown in Equations (1) and (2).
J c l = ( D c l C f + ε D H H )
J H = ( δ D c l C f + D H H )
where Dcl = chloride diffusion coefficient (cm2/day); Cf = free chloride concentration (in gram of free chloride per gram of concrete, g/g); DH = humidity diffusion coefficient; H = pore relative humidity; ε = humidity gradient coefficient, which represents the coupling effect of moisture diffusion on chloride penetration; and δ = chloride gradient coefficient, which represents the coupling effect of chloride ions on moisture diffusion.
When chloride ions ingress into the concrete, some of them are bound to the internal surface of the cement paste and aggregates, which are called bound chlorides and freely go through the concrete. Steel corrosion is related only to the free chloride content but not to the total chloride content [14].
The mass balance of chloride ions and moisture can be expressed using Fick’s second law as Equations (3) and (4),
C t t = C t C f C f t = d i v ( J C l ) = d i v ( D C l C f + ε D H H )
w t = w H H t = d i v ( J H ) = d i v ( δ D C l C f + D H H )
where
  • C t : total chloride concentration (in gram of free chloride per gram of concrete, g/g);
  • w : water content;
  • C t C f : chloride binding capacity; and
  • w H : moisture binding capacity.
Equations (3) and (4) can be rewritten as,
C t C f C f t = ( D C l C f + ε D H H ) = ( D C l C f + D ε H )
w H H t = d i v ( δ D C l C f + D H H ) = ( D δ C f + D H H )
where D ε is coupling parameter, ε D H , and D δ is coupling parameter, δ D C l .
The general boundary conditions are as below,
C f = C 0                   on    Γ 1
D c l C f n + J C l + D ε H n + α c l ( C f C f a ) = 0 on    Γ 2
H = H 0                   on    Γ 3
D H H n + J H + D δ C f n + α H ( H H a ) = 0    on    Γ 4
where α C l = convective chloride coefficient, α H = convective relative humidity coefficient, C f a = ambient chloride ions, and H a = ambient relative humidity.
Γ 1 and Γ 3 are the part of boundary with constant chloride ions and relative humidity and Γ 2 and Γ 4 are the part of boundary subjected to specified chloride ions and relative humidity flux, respectively. Γ 1 and Γ 2 form the complete boundary surface for the chloride diffusion problem, and Γ 3 and Γ 4 form the moisture diffusion problem [16].

2.2. Material Parameters

To numerically solve the chloride diffusion problem, many material parameters must be determined: the moisture capacity ( w / H ), chloride binding capacity ( C t / C f ), humidity diffusion coefficient ( D H ), and chloride diffusion coefficient ( D C l ) (in Equations (5) and (6)) [15].

2.2.1. Moisture Capacity

The moisture capacity of the concrete was developed based on the multiphase and multiscale model [17]. Assuming that the effect of the shrinkage of concrete can be evaluated simply by the average of the moisture capacities of the aggregate and the cement paste as shown in Equation (11),
w H = f a g g ( w H ) a g g + f c p ( w H ) c p
where f a g g and f c p = weight percentages of the aggregate and cement paste, and ( w H ) a g g and ( w H ) c p = moisture capacities of aggregate and cement paste, which can be calculated based on the model developed (Xi et al., 1994a, b) and Xi (1995a, b) [18,19,20,21].

2.2.2. Chloride Binding Capacity

A modified relationship between the bound chloride C b , and the free chloride C f , was established by Tang and Nilson (1993) based on Freundlich isotherm, and was proposed by Xi and Bazant (1999);
C b = β C S H 1000 ( C f 35.45 β s o l ) A 10 B
which is differentiated with respect to C f and then one would obtain the chloride binding capacity of concrete as defined by Xi and Bazant (1999) to be;
C f C t = 1 1 + A 10 B β C S H 35450 β s o l ( C f 35.45 β s o l ) A 1
where A and B = chloride adsorption related constants, 0.3788 and 1.14, βsol = ratio of pore solution volume to concrete weight, L/g, and β C - S - H = weight ratio of C-S-H gel to concrete (g/g) [14,22].
Based on the definition of βsol, the following Equation (14) can be easily derived;
β s o l = f c p n c p + f a g g n a g g ρ s o l
where f c p and f a g g = weight percentages of cement paste and aggregate in concrete mix (g/g), n c p = cement paste adsorption isotherm, n a g g = aggregate adsorption isotherm, and ρ s o l = density of pore solution measured in g/L.
Specific gravities of concrete and C-S-H are similar, therefore, Xi and Bazant (1999) assumed that the weight fraction of C-S-H in concrete, β C - S - H , is equal to the volume fraction of C-S-H in concrete, f C - S - H , then,
β C - S - H = f C - S - H = V t o t a l V 1 V c p V t o t a l = 1 f 1 f c p
where f 1 = volume fraction of anhydrous pores of cement particles, and f c p = volume fraction of capillary pores of cement paste [14].

2.2.3. Humidity Diffusion Coefficient

The humidity diffusion coefficient of concrete depends on the diffusion coefficients of aggregate and cement paste. Using the composite theory (Christensen, 1979), the effective diffusion coefficient of concrete can be evaluated as following Equation (16),
D H = D H c p ( 1 + g i / [ ( 1 g i ) 3 + 1 / ( D H a g g D H c p 1 ) ] )
where g i = aggregate volume fraction, D H c p = humidity diffusion coefficient of the cement paste, and D H a g g = humidity diffusion coefficient of the aggregates [23].
The humidity diffusion coefficient of aggregates in concrete is very small due to the fact that the pores in aggregates are discontinuous and enveloped by cement paste and so it can be neglected. The humidity diffusion coefficient of cement paste can be predicted by using the empirical formula described in Xi et al. (1994b) [19].

2.2.4. Chloride Diffusion Coefficient

Chloride diffusion coefficient in saturated concrete was studied by Xi and Bazant (1999) as following Equation (17),
D c l = f 1 ( w / c , t 0 ) f 2 ( g i ) f 3 ( H ) f 4 ( T ) f 5 ( C f )
In Equation (17), the functions of the chloride diffusivity consist of the concrete curing time ( t 0 ), gravel volume fraction ( g i ), relative humidity (H), temperature (T), and free chloride concentration ( C f ) [14].
First factor of chloride diffusion coefficient accounts for the effect of the water–cement ratio ( w / c ) and curing time ( t 0 ) in Equation (18),
f 1 ( w / c , t 0 ) = 28 t 0 62500 + ( 1 4 + ( 28 t 0 ) 300 ) ( w / c ) 6.55
The second influence factor is to consider the effect of composite action of the aggregate and the cement paste in Equation (19),
f 2 ( g i ) = D c p ( 1 + g i [ 1 g i ] / 3 + 1 / [ ( D a g g / D c p ) 1 ] )
where g i = volume fraction of aggregate in concrete, and D a g g and D c p = chloride diffusion coefficient of aggregate and cement paste.
The third factor, f3(H), in Equation (20) is to consider the effect of relative humidity level on the chloride diffusion coefficient. A model proposed by Bazant et al. (1972) can be used, which was developed initially for moisture diffusion [24].
f 3 ( H ) = [ 1 + ( 1 H 1 H C ) 4 ] 1
where H c = critical humidity level, 0.75.
Arrhenius law was used by Xi and Bazant (1999) to introduce the temperature effect of the forth factor in chloride diffusion coefficient as shown in Equation (21),
f 4 ( T ) = exp [ U R ( 1 T 0 1 T ) ]
where
  • T 0 and T = reference and current temperatures in Kelvin, T0 = 296 K;
  • R = gas constant, 8.314 J/mol·K; and
  • U = diffusion process activation energy, depending on w/c ratio [14].
The detail description of diffusion process activation energy, U can be found in the paper of Ababneh et al. (2003) [15].
The dependence of chloride diffusion coefficient on free chloride concentration C f for the fifth factor is presented in Equation (22);
f 5 ( C f ) = 1 k i o n ( C f ) m
where k i o n and m = 8.333 and 0.5, respectively, according to Xi and Bazant (1999) [14].

3. Finite Element Formulation

In order to solve a time-dependent coupled moisture-chloride diffusion problem, a large linear system equation was derived with finite element method. The finite element formulation will be briefly introduced in this section. The continuous variables in the coupled chloride and moisture diffusion equations, free chloride ( C f ) and relative humidity ( H m ) are spatially discretized over the space domain, Ω . The domain discretization can be described as shown in Equation (23).
Ω = e = 1 n e l Ω e
in which nel is the total number of elements in space domain and Ω e is an element. It is also defined Ω as the boundary of computational domain and Ω e the boundary of subdomain.
The unknown variables in Equations (24) and (25) are defined in terms of nodal values, { C f } and { H m } ,
C f N { C f }
H m N { H m }
where N is the triangle element shape function. The notations and { } are row and column vectors, respectively. The element shape functions are expressed as following Equation (26),
N = N 1 N 2 . . . N n
in which Ni is the shape function for node i and n is the total numbers of nodes in an element. The unknown vectors of free chloride { C f } in Equation (27) and relative humidity { H m } in Equation (28) can be defined as,
{ C f } { C 1 , C 2 , C 3 , , C n }
{ H m } { H 1 , H 2 , H 3 , , H n }
The nodal free chloride concentrations and relative humidity are solved by substituting the approximated values of Equations (24) and (25) into governing equations of Equations (5) and (6), and applying the Galerkin procedure to the weak forms, then the finite element matrix can be obtained as shown in Equation (29):
d d t ( [ C e ( ϕ ) ] { ϕ } ) = [ K e ( ϕ ) ] { ϕ }
where the element matrices and vector are as following Equations (30)–(32),
[ C e ] = [ C c 0 0 C h ]
[ K e ] = [ K c c K c h K h c K h h ]
{ ϕ } = C f H m
In detail, the components in element matrices are as Equations (33)–(38),
[ K c c ] = Ω e N c T D C f N c d Ω + Ω e N c T D C f N c d Γ
[ K c h ] = Ω e N c T D ε N h d Ω + Ω e N c T D ε N h d Γ
[ K h c ] = Ω e N h T D δ N c d Ω + Ω e N h T D δ N c d Γ
[ K h h ] = Ω e N h T D H m N h d Ω + Ω e N h T D H m N h d Γ
[ C C ] = Ω e N C T C C N C d Ω
[ C h ] = Ω e N h T C h N h d Ω
Finally, Equation (39) is also discretized in time space with time interval t = t ξ + 1 t ξ as following,
( [ C e ( ϕ ) θ Δ t K e ( ϕ ) ] { ϕ } ) ξ + 1 = ( [ C e ( ϕ ) ( 1 θ ) Δ t K e ( ϕ ) ] { ϕ } ) ξ
The value of parameter θ is related to the solution method adopted in the program. Typical values of θ are 0, 1/2 and 1 correspond to fully explicit, semi-implicit and fully implicit methods, respectively. The semi-implicit method called Crank–Nicholson method is used in this study.
Equation (39) is simplified as linear system equation.
[ A ] ξ + 1 { ϕ } ξ + 1 = { b } ξ
where
[ A ] ξ + 1 = [ C e ( ϕ ) θ Δ t K e ( ϕ ) ] ξ + 1
{ b } ξ = ( [ C e ( ϕ ) ( 1 θ ) Δ t K e ( ϕ ) ] { ϕ } ) ξ

4. Implementation of Parallel Finite Element Method

4.1. Various Programs Adapted in Parallel Finite Element Program

In order to implement the parallel finite element program for the coupled moisture and chloride problem, various programs were employed such as Triangle for mesh generation, ParMETIS, PETSc (Portable, Extensible Toolkit for Scientific Computation), and MPI (Message Passing Interface) [25,26,27,28,29,30,31].
In this study, Triangle was used for the mesh generation of triangle element, which was created at Carnegie Mellon University. Triangle generates exact delaunay triangulations, and are suitable for finite element analysis [25].
PETSc (3.0.0 p8) is a large and versatile package integrating distributed vectors, distributed matrices in several sparse storage formats, Krylov subspace methods, preconditioners, and Newton-like nonlinear methods with built-in trust region or linesearch strategies and continuation for robustness. It is designed to provide the numerical infrastructure for application codes involving the implicit numerical solution of PDEs, and it uses on MPI for portability to most parallel machines. The PETSc library is written in C, but may be accessed from user codes written in C, Fortran, and C++ [32].
MPI (Message Passing Interface) is a standardized and portable message-passing system designed to function on a wide variety of parallel computers. The standard defines the syntax and semantics of library routines and allows users to write portable programs in the main scientific programming languages (Fortran, C, or C++) [27,28,29].
ParMETIS extends the functionality of METIS and includes routines based on a parallel graph-partitioning algorithm that are especially suited for parallel computations and large-scale numerical simulations involving unstructured meshes. In typical FEM computations, ParMETIS dramatically reduces the time spent in interprocess communication by computing mesh decompositions such that the number of interface nodes/elements is minimized [30].

4.2. Overlapping Domain Decomposition Method with Additive Schwarz Preconditioner

A parallel program is typically developed by dividing the program into multiple fragments that can execute simultaneously, each on its own processor. In the finite element analysis, this can be accomplished by applying a domain decomposition method. Domain decomposition method is the method usually used for solving large scale system equations and it is also suitable for parallel programming because of data locality [15]. There are two types of domain decomposition methods, overlapping and non-overlapping methods. In this study, the overlapping method was employed to solve the linear sparse matrix. That is because the advantage of overlapping domain decomposition is easier to setup in algebraic approach and faster convergence than non-overlapping domain decomposition. Furthermore, the boundaries of extended subdomains are smoother than non-overlapping subdomains.
Iterative solver must be used in the iterative domain decomposition method. In this study, for the iterative solver, GMRES (Generalized Minimal Residual method) was chosen for both global and local matrix. GMRES is mainly chosen because of its ability to solve non-symmetric linear system as in the case of our problem. To improve the convergence of this problem, the additive Schwarz method preconditioner was applied. Figure 1 shows the flow chart of parallel pre-process and FE solver.

5. Numerical Results

5.1. Applied Bridge Overview

The Castlewood Canyon bridge is located on Highway 83 in the Black Forest of central Colorado. The original two-lane reinforced concrete arch bridge was built in 1946. The arches are 1.93 m wide by 1.78 m deep at the base with the depth tapering down to 1.0 m at the highest point. This bridge was severely dilapidated and was in need of repairing, enlarging, and strengthening. Parts of the concrete had spalled off of the deck, columns, and arches, and the steel rebar under concrete were severely rusted. In 2003, the original arch was repaired, the spandrel columns and decks were replaced and widened from about 0.9 m to 1.0 m including railings, and the overall length of the bridge was increased from 114 m to 123 m as shown in Figure 2.

5.2. Parallel Finite Method of Large-Scale Concrete Structure

5.2.1. Modeling of Castle Wood Canyon Bridge with Large Meshes

In order to analyze the penetration of chloride and humidity in a concrete bridge, a large number of meshes are needed to capture the diffusion phenomenon within thin layer from concrete top surface to steel rebar. For depicting in detail concrete cover depth on rebar, the information of nodes and elements are created by mesh generator. Triangle as mesh generator is specialized for creating two-dimensional finite element meshes.
For castle wood canyon bridge, approximately 1.5 million nodes and three million elements were generated for input files as shown in Figure 3. To visualize and check the mesh size and shape, ParaView was employed as illustrated in Figure 4 [33]. ParaView is an open-source, multi-platform data analysis and visualization application. It is developed to analyze extremely large datasets using distributed memory computing resources and can be run on supercomputers to analyze datasets of terascale as well as on laptops for smaller data. For large scale information, VTK file format is adopted because this format is easy to read and write by hand or programmatically. This file format is automatically created when the program runs. For partitioning origin meshes, Parmetis was used to be embedded in the parallel program to automatically partition. Parmetis is an MPI-based parallel library that implements a variety of algorithms for partitioning unstructured graphs, meshes, and for computing fill-reducing orderings of sparse matrices. ParMETIS extends the functionality provided by METIS and includes routines that are especially suited for parallel AMR computations and large scale numerical simulations. The algorithms implemented in ParMETIS are based on the parallel multilevel k-way graph-partitioning, adaptive repartitioning, and parallel multi-constrained partitioning schemes developed.

5.2.2. Environmental Humidity Model

Environmental humidity model plays an important role of diffusion of chloride in concrete. The humidity model used in this study was proposed by Bazant and Xi (1993) [13]. Even though the moisture and temperature transport in concrete were coupled, humidity fluctuation was considered and temperature was assumed to be constant during analysis time. To establish a practicable and realistic humidity input model, real climate record in Chicago was employed as shown in Figure 5.
The environmental humidity model consists of three different components as shown in Equation (43) and demonstrated in Figure 6. First, H 1 is the mean value, standing for the stable horizontal trend. Second, H 2 is a random phase process corresponding to the harmonic variation of humidity. Third, H 3 is a random normal distribution with a specific mean and variance [13].
H = H 1 + H 2 + H 3
where H 1 = H ¯ , mean value of environmental model, H 2 = A 1   cos ( 2 π   t   + φ 1 ) +   A 365   cos ( 2 π   t   + φ 365 ) , the random phase process of a one day period with A 1 = 0.1 and the random phase process of a one year period with A 365 = 0.08, φ 1 and φ 365 have uniform distributions with constant density 1 2 π , H 3 = 1 2 π σ 2 e [ ( x μ ) 2 2 σ 2 ] , normally distributed random numbers are generated with specific mean and variance, 0.0046 [13].
For the chloride concentration on the top surface of concrete, various boundary conditions were used to actualize the realistic boundary conditions as illustrated in Figure 7. One can see that the interference between the interfaces of different concentration.

5.2.3. Speed-Up for Parallel Algorithm

Speed-up is measured for the performance of the parallel implementation of finite element program. The definition of speed-up is a process for increasing in performance between two systems processing the same problem. With regard to speed-up in this study, 8–2048 processors were used in this study and the meshes with the number of about 1.5 million nodes were analyzed. When using up to eight processors, the program was not operated due to memory capacity when input data can be assigned on each memory which means the problem was too huge to solve with a couple of processors.
As the number of processors used in the analysis increased, speed-up was improved up to 512 processors because speed-up was better than ideal condition. After more than 512 processors were used, speed-up gradually increased because the communication time between the processors increased. As shown in Figure 8, one can see the optimum point of number of processors.

5.2.4. Effect of Boundary Condition

(1) Constant boundary condition
The physical and material models of the castle wood canyon bridge used in this paper were mentioned above. The validation of these models was secured by Ababneh et al. and Suwito et al. [14,15]. Ababneh et al. proved the model with alternating-direction implicit (ADI) finite-difference method and the numerical solutions were compared with the experimental results obtained by the 90-day ponding test (AASHTO T 259-80) [14]. Suwito et al. verified the accuracy of the model implemented with finite element method and numerical results agreed well with the experimental data [15].
It was assumed that the concrete bridge contained initially no chloride ions and has 60% relative humidity (RH). The top surface of concrete bridge was exposed to 3% NaCl and 100% RH. To analyze the diffusion of coupled chloride and humidity, 8 to 2048 processors and approximately 1.5 million nodes were used with parallel finite element method.
Obviously, this example is not to simulate the concrete behavior under service condition. The triangle element meshes were employed and the whole domain was partitioned into 16 sub-domains to the same as the number of processors shown in Figure 3 and Figure 4. The sub-domains were divided with ParMetis so that each sub-mesh had almost the same number of elements.
Figure 9 and Figure 10 show the contour of chloride and humidity distribution into the entire domain. The total and free chloride concentration at 50 mm below from the top surface of concrete is Figure 11 and the variation of relative humidity is Figure 12 over time up to 550 days. As a result of analysis, the concentration of chloride and humidity is continuously accumulated inside concrete. Specially, the variation of humidity inside concrete might not be increased in reality due to the seasonally change of the humidity distribution, so that periodic humidity model should be applied as the boundary condition.
(2) Periodic boundary condition
In order to simulate the real diffusion phenomenon in concrete, the periodic boundary condition should be required. Figure 13 and Figure 14 describe the result of chloride concentration and relative humidity at different depth over time. As shown in Figure 13, total concentration is gradually increasing until almost one year. That means the chloride concentration can be affected by the periodic humidity boundary condition. After one year it suddenly increases because of the influence of the constant chloride initial condition. However, the total chloride concentration at 50 mm is reduced up to 15% rather than the concentration using constant boundary. In order to predict and demonstrate the realistic and seasonal chloride ingress through concrete cover, seasonal change information of chloride concentration on the top surface of concrete deck.
In Figure 13, the change of humidity distribution is influenced from the humidity model inputted as boundary condition. The short period random noise, one of random noise terms in environmental humidity model, can only reach a shallow portion near the top surface of concrete deck. As the concrete depth increased, the effect of short period random noise decreases. However, the seasonal variation of humidity concentration has an influence on 50 mm depth of concrete deck. At 50 mm from the top surface of concrete deck, the maximum value of humidity is about 77%, which is decreased up to 20% comparing with constant boundary condition. In addition, the result of about 500 days shows that the humidity trend is reversed between 10 mm and 50 mm, which means that evaporation of humidity just near concrete surface is easier than evaporation at 50 mm. Therefore, this prediction model is reasonable and effective to describe the diffusion phenomenon of chloride and humidity.

6. Conclusions

(1)
The parallel finite element program was developed based on the robust mathematical material model. The program can be used to simulate the coupled moisture and chloride penetration into non-saturated concrete structures. Chloride ion is one of sources causing the steel corrosion in reinforced concrete structures. The material parameters related to chloride and moisture diffusion in concrete are taken into account. These parameters include chloride and moisture diffusion coefficients, moisture capacity, chloride binding capacity, and the coupling parameters reflecting the coupling effects between moisture and chloride transfer in concrete.
(2)
For the implementation of parallel FE analysis, Triangle for mesh generation, ParMetis, PETSc, and MPI are employed. This program also used the overlapping domain decomposition method with additive Schwarz preconditioner. As the result of simulation, the computation time decreased until the number of processors became optimized when the number of processors increased. Then, the computational time increased because the communication between the processors increased.
(3)
The present model can be used to simulate the unsaturated concrete structures subjected to other aggressive chemicals from de-icing salt. The framework of present model can be extended to simulate the multi-species de-icing salts ingress into non-saturated concrete structures in future work.

Acknowledgments

The authors wish to acknowledge the partial support by the US National Science Foundation under Grant CNS- 0722023 to University of Colorado at Boulder. Opinions expressed in this paper are those of the authors and do not necessarily reflect those of the sponsor.

Author Contributions

Okpin Na, he did all of the numerical implements, analyzed the numerical data, and wrote the manuscript. Xiao-Chuan Cai advised Okpin Na to code the parallel algorithm with C++ and PETSc and did the analysis of the numerical results. Yunping Xi directed and advised the theoretical research work, and revised the whole manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ali, H.A. Numerical simulation of chloride diffusion in RC structures and the implications of chloride binding capacities and concrete mix. IJCEE 2010, 10, 22–35. [Google Scholar]
  2. Masi, M.; Colella, D.; Radaelli, G.; Bertolini, L. Simulation of chloride penetration in cement-based materials. Cem. Concr. Res. 1997, 27, 1591–1601. [Google Scholar] [CrossRef]
  3. Conciatori, D.; Sadouki, H.; Bruhwiler, E. Capillary suction and diffusion model for chloride ingress into concrete. Cem. Concr. Res. 2008, 38, 1401–1408. [Google Scholar] [CrossRef]
  4. Conciatori, D.; Laferriere, F.; Bruhwiler, E. Comprehensive modeling of chloride ion and water ingress into concrete considering thermal and carbonation state for real climate. Cem. Concr. Res. 2010, 40, 109–118. [Google Scholar] [CrossRef]
  5. Isgor, O.B.; Razaqpur, A.G. Finite element modelling of coupled heat transfer, moisture transport and carbonation processes in concrete structures. Cem. Concr. Compos. 2004, 26, 57–73. [Google Scholar] [CrossRef]
  6. Lin, G.; Liu, Y.H.; Xiang, Z.H. Numerical modeling for predicting service life of reinforced concrete structures exposed to chloride environments. Cem. Concr. Compos. 2010, 32, 571–579. [Google Scholar] [CrossRef]
  7. Chernin, L.; Val, D.V.; Stewart, M.G. Prediction of cover crack propagation in RC structures caused by corrosion. Mag. Concr. Res. 2012, 64, 95–111. [Google Scholar] [CrossRef]
  8. Zhang, J.; Ling, X.; Guan, Z. Finite element modeling of concrete cover crack propagation due to non-uniform corrosion of reinforcement. Constr. Build. Mater. 2017, 132, 487–499. [Google Scholar] [CrossRef]
  9. Lim, S.; Akiyama, M.; Frangopol, D.M. Assessment of the structural performance of corrosion-affected RC members based on experimental study and probabilistic modeling. Eng. Struct. 2016, 127, 189–205. [Google Scholar] [CrossRef]
  10. Šomodíková, M.; Lehký, D.; Doležel, J.; Novák, D. Modeling of degradation processes in concrete: Probabilistic lifetime and load-bearing capacity assessment of existing reinforced concrete bridges. Eng. Struct. 2016, 119, 49–60. [Google Scholar] [CrossRef]
  11. Hu, N.; Burgueño, R.; Haider, S.W.; Sun, Y. Framework for Estimating Bridge-Deck Chloride-Induced Degradation from Local Modeling to Global Asset Assessment. J. Bridge Eng. 2016, 21, 06016005. [Google Scholar] [CrossRef]
  12. Tang, L.; Nilson, L.O. Prediction of Chloride Penetration into Concrete by Using the Computer Program CLINCON, Proceedings of the Second International Conference on Concrete under Severe Conditions 2-Environment and Loading, Tromsø, Norway, 21–24 June 1998; Gjørv, O.E., Sakai, K., Banthia, N., Eds.; E & FN Spon: New York, NY, USA, 1998; pp. 625–634. [Google Scholar]
  13. Bazant, Z.P.; Xi, Y. Stochastic Drying and Creep Effects in concrete structures. J. Struct. Eng. 1993, 119, 301–322. [Google Scholar] [CrossRef]
  14. Xi, Y.; Bazant, Z.P. Modeling chloride penetration in saturated concrete. J. Mater. Civ. Eng. 1999, 11, 58–64. [Google Scholar] [CrossRef]
  15. Suwito, A.; Cai, X.; Xi, Y. Parallel finite element method for coupled chloride moisture diffusion in concrete. Int. Numer. Anal. Model. 2006, 3, 481–503. [Google Scholar]
  16. Ababneh, A.; Benboudjema, F.; Xi, Y. Chloride penetration in Nonsaturated concrete. J. Mater. Civ. Eng. 2003, 15, 183–191. [Google Scholar] [CrossRef]
  17. Xi, Y.; Willam, K.; Frangopol, D.M. Multiscale modeling of interactive diffusion processes in concrete. J. Eng. Mech. 2000, 126, 258–265. [Google Scholar] [CrossRef]
  18. Xi, Y.; Bazant, Z.P.; Molina, L.; Jennings, H.M. Moisture diffusion in cementitious materials: Adsorption isotherm. Adv. Cem. Based Mater. 1994, 1, 248–257. [Google Scholar] [CrossRef]
  19. Xi, Y.; Bazant, Z.P.; Molina, L.; Jennings, H.M. Moisture diffusion in cementitious materials: Moisture capacity and diffusivity. Adv. Cem. Based Mater. 1994, 1, 258–266. [Google Scholar] [CrossRef]
  20. Xi, Y. A model for moisture capacities of composite materials. I: Formulation. Comput. Mater. Sci. 1995, 4, 65–77. [Google Scholar]
  21. Xi, Y. A model for moisture capacities of composite materials. II: Application to concrete. Comput. Mater. Sci. 1995, 4, 78–92. [Google Scholar] [CrossRef]
  22. Tang, L.; Nilson, L.O. Chloride binding capacity and binding isotherms of OPC pastes and mortars. Cem. Concr. Res. 1993, 23, 247–253. [Google Scholar]
  23. Christensen, R.M. Mechanics of Composite Materials, 2nd ed.; Wiley: New York, NY, USA, 1979. [Google Scholar]
  24. Bazant, Z.P.; Najjar, L.J. Nonlinear water diffusion of nonsaturated concrete. Mater. Constr. 1972, 5, 3–20. [Google Scholar] [CrossRef]
  25. Triangle. Available online: https://www.cs.cmu.edu/~quake/triangle.html (accessed on 27 March 2017).
  26. Balay, S.; Buschelman, K.; Eijkhout, V.; Gropp, W.D.; Kaushik, D.; Knepley, M.G.; McInnes, L.C.; Smith, B.F.; Zhang, H. PETSc-users manual. In Technical Report ANL-95/11–Revision 2.1.5; Argonne National Laboratory: Argonne, IL, USA, 2008. [Google Scholar]
  27. MPI Forum. MPI: A message passing interface standard. Int. J. Supercomput. Appl. 1994, 8, 159–416. [Google Scholar]
  28. MPI Forum. MPI2: A message passing interface standard. High Perform. Comput. Appl. 1998, 12, 1–299. [Google Scholar]
  29. MPICH Team. MPICH: A Portable Implementation of MPI. 1996–2005. Available online: http://www.mpich.org/documentation/guides/ (accessed on 27 March 2017).
  30. Karypis, G. ParMETIS: Parallel Graph Partitioning and Sparse Matrix Ordering. 1996–2005. Available online: http://glaros.dtc.umn.edu/gkhome/metis/parmetis/download (accessed on 27 March 2017).
  31. Dalcín, L.D.; Paz, R.R.; Anca, A.A.; Storti, M.A.; D’Elía, J. Parallel FEM application development in Python. Mec. Comput. 2005, 24, 1823–1838. [Google Scholar]
  32. PETSc. Available online: http://www.mcs.anl.gov/petsc/documentation/index.html (accessed on 27 March 2017).
  33. Paraview. Available online: http://www.paraview.org/ (accessed on 27 March 2017).
Figure 1. Framework of parallel FE method based on PETSc: (a) flow chart for parallel pre-process; and (b) flow chart for parallel FE solver.
Figure 1. Framework of parallel FE method based on PETSc: (a) flow chart for parallel pre-process; and (b) flow chart for parallel FE solver.
Materials 10 00350 g001
Figure 2. Castlewood Canyon bridge, Franktown, Colorado.
Figure 2. Castlewood Canyon bridge, Franktown, Colorado.
Materials 10 00350 g002
Figure 3. Modelling of Castlewood Canyon bridge.
Figure 3. Modelling of Castlewood Canyon bridge.
Materials 10 00350 g003
Figure 4. Partitioning of Castlewood Canyon bridge with 16 processors.
Figure 4. Partitioning of Castlewood Canyon bridge with 16 processors.
Materials 10 00350 g004
Figure 5. Environmental humidity Record, midway station, Chicago (Bazant et al., 1993) [13].
Figure 5. Environmental humidity Record, midway station, Chicago (Bazant et al., 1993) [13].
Materials 10 00350 g005
Figure 6. Environmental humidity model: (a) humidity model without random noise; and (b) humidity model with three components (including random noise).
Figure 6. Environmental humidity model: (a) humidity model without random noise; and (b) humidity model with three components (including random noise).
Materials 10 00350 g006
Figure 7. Comparison of single and various boundary conditions: (a) single boundary condition; and (b) two boundary conditions.
Figure 7. Comparison of single and various boundary conditions: (a) single boundary condition; and (b) two boundary conditions.
Materials 10 00350 g007
Figure 8. Speed up over number of processors.
Figure 8. Speed up over number of processors.
Materials 10 00350 g008
Figure 9. Contour of chloride distribution.
Figure 9. Contour of chloride distribution.
Materials 10 00350 g009
Figure 10. Contour of humidity distribution.
Figure 10. Contour of humidity distribution.
Materials 10 00350 g010
Figure 11. Chloride concentration over time at 50 mm depth.
Figure 11. Chloride concentration over time at 50 mm depth.
Materials 10 00350 g011
Figure 12. Humidity over time at 50 mm depth.
Figure 12. Humidity over time at 50 mm depth.
Materials 10 00350 g012
Figure 13. Chloride concentration over time at 50 mm depth.
Figure 13. Chloride concentration over time at 50 mm depth.
Materials 10 00350 g013
Figure 14. Humidity over time according to depth.
Figure 14. Humidity over time according to depth.
Materials 10 00350 g014

Share and Cite

MDPI and ACS Style

Na, O.; Cai, X.-C.; Xi, Y. Corrosion Prediction with Parallel Finite Element Modeling for Coupled Hygro-Chemo Transport into Concrete under Chloride-Rich Environment. Materials 2017, 10, 350. https://doi.org/10.3390/ma10040350

AMA Style

Na O, Cai X-C, Xi Y. Corrosion Prediction with Parallel Finite Element Modeling for Coupled Hygro-Chemo Transport into Concrete under Chloride-Rich Environment. Materials. 2017; 10(4):350. https://doi.org/10.3390/ma10040350

Chicago/Turabian Style

Na, Okpin, Xiao-Chuan Cai, and Yunping Xi. 2017. "Corrosion Prediction with Parallel Finite Element Modeling for Coupled Hygro-Chemo Transport into Concrete under Chloride-Rich Environment" Materials 10, no. 4: 350. https://doi.org/10.3390/ma10040350

APA Style

Na, O., Cai, X. -C., & Xi, Y. (2017). Corrosion Prediction with Parallel Finite Element Modeling for Coupled Hygro-Chemo Transport into Concrete under Chloride-Rich Environment. Materials, 10(4), 350. https://doi.org/10.3390/ma10040350

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop