Abstract

In order to further understand the carbonation process of concrete structures, the time- and temperature-dependent diffusion process of CO2 in concrete is simulated based on the law of the CO2 mass conservation, and a two-dimensional mass transfer equation is established for the CO2 diffusion in concrete. The concrete block is discretized into triangular elements, and the CO2 concentrations at different positions are calculated based on finite element method. A computational algorithm is programed through the Matlab platform. The time- and temperature-dependent property and difference of the CO2 concentration at different positions of the structure are considered in the proposed model. Then, an accelerated carbonation experiment is carried out using concrete blocks with different mix proportions to investigate the influence of the water-cement ratio and temperature on the concrete carbonation. The experimental results effectively verify the correctness of the finite element model, and the proposed finite element method reasonably simulates the concrete carbonation through calculating the carbonation in practical engineering compared with other methods in references. An experimental-numerical correlation has been performed. The ratio of carbonation depth at the corner of the concrete members to the other positions is about 1.35. The carbonation depth is increased about 1.9 times when the temperature changes from 20°C to 40°C.

1. Introduction

CO2 concentration is an important factor affecting carbonation in the atmospheric environment. When the CO2 in the air permeates into concrete, the acid-base neutralization reaction happens between the CO2 and the alkaline materials, which leads to the decrease of the concrete alkalinity and the change of chemical composition. When the carbonation front reaches the surface of steel bars, the passivation film on the surface of steel rebar will be destroyed, and then, the corrosion of steel rebars occurs. Overtime, along with the increase of the corrosion products, it is well recognized that the volume of corrosion product is 3∼8 times larger than that of the steel material. This will result in cover cracking and even spalling and also lead to a reduction of the structural load carrying capacity. These phenomena affect the working performance and durability of concrete structures.

With the development of global economy and the sharp increase of population, the CO2 concentration also increases rapidly. Before the second industrial revolution, the CO2 concentration was about from 265 ppm to 290 ppm in the global atmosphere [1], was about 365 ppm in 2000 [2] (1 ppm = 0.00188 × 10−3 kg/m3), and reached 400 ppm in 2015, which broke the historical extreme. According to the prediction of International Panel on Climate Change (IPCC), the CO2 concentration will exceed 1000 ppm by 2100 [2]. The CO2 concentrations under various different emission strategies are demonstrated in Figure 1 (A1F1, B1, and Best mitigation) [3]. A1F1 describes a world of very rapid economic growth, a global population that peaks in midcentury, and rapid introduction of new and more efficient technologies that remain fossil intensive, and this is the worst CO2 emission scenario. B1 describes a world with intermediate population and economic growth, emphasizing local solutions to economic, social, and environmental sustainability, and this is the lowest CO2 emission scenario. Best mitigation describes that the CO2 concentration is kept stable at 2010 levels (386 ppm) due to stabilization and reduction of CO2 emissions. The climate change leads to the increase of CO2 concentration, which accelerates the carbonation process of concrete structures, and it is seriously harmful to the structural durability and safety. In present experimental and theoretical studies, the effects of the CO2 concentration on the carbonation of concrete have already been explored by Visser [4], Cui et al. [5], and Castellote et al. [6]; however, the above scholars do not consider the time-dependent character of the CO2 concentration, which will affect the precision of the prediction of the carbonation depth, so it is necessary to consider the time-dependent character of CO2 concentration resulting from climate change when investigating the carbonation progress.

The greenhouse effect has become globalization problem in recent years. According to the prediction of atmospheric temperature from China’s National Meteorological Center, the change of temperature in the Beijing district is shown in Figure 2. The average value, lower value and upper value of temperature are presented in Figure 2. Jianxin et al. [7] stated that the diffusion coefficient of concrete carbonation will improve with the increase of temperature, so the temperature change should be considered when predicting the carbonation progress. However, in existing researches [811], the carbonation progress rarely considers the change of temperature. This leads to underestimate concrete carbonation behavior. Therefore, the influence of the temperature change on the carbonation progress should be considered in theoretical and experimental analysis.

Except the CO2 concentration and temperature, the effects of other environmental factors and the concrete material compositions on concrete carbonation have been investigated in-depth in existing researches. The environmental factor includes relative humidity [1114], and the compositions include water-cement ratio [15], various supplementary cementing materials [1618], aggregates [19, 20], and mix proportions [21, 22]. In these studies, there are few researches about the carbonation situation on the carbonation surface. However, the carbonation degree at the different place on the carbonation face is different, which is meaningful for the research of spatial reliability [23], so it is necessary to carry out the related research.

In recent years, with the rapid development of computer science, many scholars have made numerical simulations of concrete carbonation. Saetta et al. [24, 25] proposed a model to simulate all the phenomena based on carbonation process considering the multidimensional transport of moisture, heat, and carbon dioxide through concrete and taking in account the reduction of the porosity due to carbonation. Burkan Isgor and Razaqpur [10] built a nonlinear finite element approach for tracing the spatial and temporal advancement of the carbonation front in concrete structures with and without cracks. However, the above two models do not consider the time- and temperature-independence of the carbonation coefficient. Saetta et al. [11, 12] established the differential equation of the carbonation process considering the factors such as water, CO2 diffusion, and the variation of temperature and applied to practical engineering for predicting the carbonation depth of a prestressed RC beam, while the detail numerical simulation progress was not carried out for the control equation of the carbonation reflection. Peter et al. [26] proposed the system of reaction–diffusion equations independent of the space dimension, while their numerical simulations only deal with the situations of one-dimensional diffusion. Wang et al. [27] developed the numerical simulation model of the concrete carbonation and determined the value range of the factors of the CO2 diffusion coefficient and the reaction rate of carbonation. However, only two factors were considered in the model, and the model was realized by the secondary development of software. As a result, its applicability and versatility were limited. Zha et al. [28] proposed a mathematical model to simulate the physical–chemical coupling process of supercritical carbonation in cement-based materials. This model took into account the rate of chemical reaction, mass conservation for gas-liquid two-phase flow, diffusion and dispersion of CO2 in water, energy conservation for porous medium, and the solubility of CO2 in water, while this model did not consider the influence of the time and temperature variation.

Based on the existing carbonation test and numerical simulation, this research firstly builds a finite element model to predict the carbonation depth considering the time (t) and temperature (T) dependent carbonation diffusion coefficient and carbonation position difference for concrete structures. An experiment is also carried out to investigate the influence of multiple factors such as water-cement ratio, temperature, humidity, and CO2 concentration on the concrete carbonation. The accuracy of the finite element model is also validated by the experimental results and the practical measuring value in engineering.

2. Numerical Simulation of the Concrete Carbonation considering Time and Temperature Effects

2.1. Mass Transfer Diffusion Equation of CO2 in Concrete

Carbonation is a series of physical and chemical reactions between acid gases such as CO2 and hydration products of the concrete. The acid gas enters the concrete through pores in concrete. In order to obtain more accurate prediction values of the carbonation depth, it is necessary to simulate the real diffusion process of CO2 in concrete and establish the related mass transfer equation to obtain the value of CO2 mass concentration in concrete.

A control body (microunit with a constant volume) was selected as the research object in the structure. The microunit is cuboid. The side length of the microunit was assumed as dx, dy, and dz, respectively. According to the principle of mass conservation, a balance equation for CO2 gas is

According to Fick’s diffusion law, the diffusion is the main transferring behavior for CO2 in concrete. Considering only the diffusion flux generated by the concentration gradient in various directions and ignoring the impact of the flowing of CO2 gas, the rate of CO2 gas influx in mass can be formulated aswhere , and are the diffusion fluxes in the x-, y-, and z-axis direction, respectively. Their units are all kg·s−1·m−2. , , and are the values of , , and in the x-, y-, and z-axis direction, respectively, as shown in Figure 3.

The corresponding quality rate of CO2 gas efflux in mass in the microunit can be expressed as

The rate of cumulative CO2 in mass iswhere is the mass concentration of CO2.

If the rate of CO2 mass consumption per unit volume is , the rate of CO2 gas mass consumption in the microunit is

Substituting (3)–(5) into (2), the equilibrium equation can be given as

Taking the limit values for , , and as being closer to zero, (6) can be expressed as

Assuming the materials in the concrete structure are isotropic, the diffusion fluxes of CO2 in concrete can be calculated by Fick’s first law aswhere is the diffusion coefficient and is the concentration gradient. The minus sign means the diffusion direction is contrary to the concentration gradient.

The CO2 diffusion process in concrete is unsteady, and its concentration varies with the distance x and the diffusion time t. Combining (8) with (7), we get

Most existing studies only consider the carbonation value of concrete under a certain apparent diffusion coefficient. In fact, with the development of carbonation and hydration, the porosity of the carbonation region becomes smaller, and the diffusion rate decreases. According to the research results from the Yoon et al. [29] and Jianxin et al. [7], the diffusion coefficient is time and temperature dependent. Therefore, this paper considers the time- and temperature-dependent diffusion coefficient, and the calculation model iswhere is the CO2 diffusion coefficient after a year, is the age of CO2 diffusion coefficient, is the activation energy during the diffusion progress (40 kJ/mol), and is the gas constant (0.008134 kJ/mol·K).

In practical engineering, the length of concrete structure is generally large. For convenience, the CO2 diffusion problem can be simplified from three-dimensional to two-dimensional, which sets dz = 1 in the control body and only considers diffusion along the x and y directions. When the mass transfer of CO2 gas in concrete is dominated by diffusion, the reaction process ends in the carbonization region, and the rate of CO2 mass consumption is zero (). Therefore, (9) can be simplified as

Equation (11) is the diffusion equation of CO2 gas in concrete, which demonstrates the diffusion law of molecules in the unsteady state. is the Laplace operator. The solution of this equation is the CO2 concentration function at time . The concentration values at any time can be obtained by solving this equation. Finally, the concentration distribution of CO2 gas in the concrete can be obtained. The calculation process can be achieved by using the Matlab program.

2.2. Finite Element Solution of Diffusion and Mass Transfer Equations

The main chemical reaction of concrete carbonation is that CO2 penetrates into the pore of concrete and interacts with CH, C3S, C2S, and other basic hydration products produced by the hydration process of cement. Finally, calcium carbonate and other substances are formed.

Following conditions are generally known during the reaction and presented as(i)Initial condition: when t = 0, is a certain value.(ii)Boundary condition: when x = y = 0,  =  (i.e., the concentration of CO2 is equal to the initial mass concentration). When x =  and y = ( and tends to infinity),  = 0.(iii)According to Yoon’s model [29], the carbonation front equals carbonation depth and the concentration of CO2 is zero at the end of the carbonation depth, and this setting is applied in this paper.

The following differential equation is established considering boundary conditions:where is the boundary (x = 0, y = 0) and is the constant. This equation is a parabolic partial differential equation.

According to the variation principle and Euler equation, the corresponding function of (12) is

Hence, solving the partial differential solution (13) can be converted to obtain the extreme value of function in the carbonation region , and the mass concentration of on the boundary of the region is known.

According to the function form, the region is divided into several small triangular meshes, and a typical unit is taken out for analysis, which is the six-node triangular element [30] (Figure 4), and the node (a, b, c) on the edge of the triangular element is located at the middle point of the corresponding edge. Equation (13) can be expressed in the area of the unit as

The mass concentration in the unit i is discretized into a concentration function which is only related to the concentration at the node. The column vector of the mass concentration of node is

A quadratic concentration function of a six-node triangular element is set as

By using the area coordinate Li (i, j, and m), the element concentration function can be expressed by shape functions as

The matrix form of is

The elements in the shape function matrix can be expressed aswhere is the area of the triangular element.

Assuming that is only the function of the position (x, y), the functional variation at a particular time and the variation of time parameter t are considered. The variation calculation of (14) is

The derivative of (17) with respect to each mass concentration of node is

Considering the unsteady diffusion process of CO2 in the concrete, is a function of time parameter t. Therefore, we have

By combining (20) to (22) and term by term integration of (20), the following equations can be obtained:

It can be obtained from the above equation as

The coefficient matrices and are calculated as

Other factors can be recursive based on the above equations.

Based on the forward difference method, the concentration of CO2 can be given as

Substituting (26) into (27) gives

According to the finite element method, the carbonation area is discretized into finite triangular elements and n nodes, and the mass concentration of the corresponding elements is distributed to each node. Therefore, it is necessary to solve the concentration values at each node and then solve the mass concentration distribution in the area . is the function of triangular elements, and the function of the entire region is . The function is a multivariate function to express the mass concentration of the unknown nodes. Combining with the variation principle, solving functional extremum can be converted to solve the extremum value of multivariate functions. The extreme condition for the function is

The coefficient matrixes are synthesized as and , respectively, and the recursion equation can be obtained by combining (27) aswhere is the initial mass concentration of CO2, and the concentration of the CO2 value at any time interval () can be obtained through setting the initial concentration of CO2 and solving (29).

2.3. Modeling of Carbonation Process

The finite element method is implemented in the Matlab platform, and the flow chart of the program is shown in Figure 5. The carbonation surface of the test specimen is simplified (75 mm × 75 mm), and a simulation process at a carbonation age is shown in Figure 6. The direction of arrow represents the diffusion direction of CO2 gas after penetrating into the concrete and along the straight line OA (Figure 6). Areas from red to blue represent the transition from the carbonated zone to the partially carbonated zone and to the noncarbonated zone, respectively. The L-shaped curves represent the distribution curve of CO2 concentration in the concrete, and the values of curves are , 0.9, … , 0.1, respectively. Carbonation depth values at the corner and other positions are calculated for comparison, and the carbonation depth values at the corner (Xj) and other positions (Xq) for C25, C30, and C35 concrete are obtained under the condition of three kinds of CO2 emissions strategy (A1F1, B1, and Best mitigation) [3] for the next 100 years from 2000 as shown in Figure 7. The corner refers the line OA, and other positions are the upper triangular and lower triangular positions except the corner. It can be found that the ratio of the carbonation depth value at the corner (Xj) to that at the other positions (Xq) is between 1.3 and 1.4. Therefore, when predicting the carbonation depth of concrete structures, it is necessary to consider the problem of different carbonation depths at different positions.

3. Experimental Test Programs

3.1. Material Properties

The sizes of standard concrete blocks in this test were 150 mm × 150 mm × 150 mm. Three different concrete water-cement ratios (0.65, 0.55, and 0.45) were designed. The designed concrete compressive strengths in this experiment were 25 MPa, 30 MPa, and 35 MPa, respectively. The mix compositions for concrete blocks are shown in Table 1. The blocks were taken out after the 28 days standard curing period and baked for 48 hours at 60 ± 2°C and then put into a closed carbonation chamber (Figure 8). In addition, the test also used 1% phenolphthalein ethanol solution to determine the carbonation depth.

3.2. Test Instruments

The test was carried out in both the high- and low-temperature carbonation chambers. The CO2 gas supply device (including gas cylinders and pressure gauge) was used to maintain the CO2 concentration in the chamber. The test block was broken by the universal testing machine. The carbonation depths were measured by the digital carbonation depth detector (HC-TH01).

3.3. Test Parameters

Based on the different concrete strengths, water-cement ratio (), temperature, and test time, nine groups of the concrete samples were set up in the test. There were totally fifteen samples in each group that is divided into five subgroups. In each subgroup, three concrete samples were tested for every test time. Table 2 presents the detailed setting of each parameter for each group in the test.

3.4. Test Progress

One side of the block was tested for carbonation, while the remaining five surfaces were protected with epoxy resin, and the space between components was larger than 5 cm as shown in Figure 8. These samples were placed in a chamber containing 20 ± 3% concentration of CO2 at 20 (30 or 40) ± 2°C and 70 ± 5% relative humidity. During the test, the temperature and humidity were measured and recorded every four hours, and the CO2 mass concentration on the surface of the block in the chamber was kept constant. These blocks were taken out after designed time (Table 2) and were broken by a universal testing machine from the center of the uniaxial carbonation surface to measure the carbonation depth. After the test, the ear syringe was used to blow away the ash and slag on the concrete surface, and then the phenolphthalein ethanol solution was dropped on the concrete surface for about 30 seconds. The carbonation depth value of each measurement point whose space was set as 1 cm on the surface of concrete block was measured by the carbonation depth detector. The detailed carbonation schematic diagram is shown in Figure 9(a), the uniaxial carbonation surface is divided into corner areas (Xj) and other areas except corner areas (Xq) as shown in Figure 9(b).

When the coarse aggregate is just located at the dividing line of a certain measured point, the average carbonation depth of the tested block on the left and right sides of the coarse aggregate is taken as the depth of the point xi. The formula for the average carbonation depth of concrete at different ages is calculated aswhere is the average carbonation depth of the test blocks after t days and n is the total point numbers on two sides surface.

4. Test Results and Discussion

4.1. Influence of Carbonation Position

According to the experimental data, Figure 10 presents the carbonation depth curve in the corner areas (Xj) and other areas (Xq) for concrete with three different water-cement ratios under the conditions of different temperatures. Due to the continuous reaction of the alkaline hydration product and carbon dioxide in the pore of concrete, which generates the difficult soluble substances CaCO3 and blocking gas channel, the carbonation recreation is hindered. This leads to a larger slope for the initial curve, while a smaller slope for the later carbonation depth curve. Under the condition of different water-cement ratios, the carbonation depth of concrete at the corner position is about 1.35 times larger than that at the other positions, which is in good agreement with the finite element results and verifies the correctness of calculation results. The reason of the difference between corner and other area is that the carbonation rate of each direction is different in the closed carbonation chamber, and the corner position is subjected to bidirectional carbonation corrosion on both sides of the corner position.

4.2. Influence of the Ratios

Figure 11 presents the relationship between the carbonation depth and water-cement ratio at the age of 28 and 84 days, respectively. Compared with different water-cement ratios of concrete, the carbonation depth of C25 concrete is less than that of the rest of the two other kinds of concrete under the same temperature and age. This is because the greater concrete water-cement ratio will result in the greater porosity which leads to the faster carbonation rate. Through the data fitting for the tested result, the relationship between the carbonation rate coefficient k and the water-cement ratio of each group is linear, as shown in Figure 12. This is consistent with the experimental results reported by Jin et al. [31].

4.3. Influence of Temperature

Figure 13 shows the relationship between the carbonation depth and temperature at the ages of 7 and 56 days, respectively. It can be found from Figure 6 that the carbonation rate increases as the temperature increases. Under the condition of the same humidity and CO2 concentration, the carbonation depth would increase about 1.9 times with an increasing temperature from 20°C to 40°C for the same concrete grade, which is consistent with the data measured by Li et al. [32]. According to the formula of temperature influence coefficient () proposed by Uomoto and Takada [33], using Origin software to fit the formula , the values of fitting coefficient q is given in Table 3.

It can be observed from Table 3 that the fitting coefficient q is basically consistent for the same concrete whether the carbonation depth is at the corner areas (Xj) or other areas (Xq). The values of the coefficient of carbonation rate for standard specimens at different temperatures can be calculated by the fitting formula, as shown in Table 4.

It can be obtained from Table 4 that the coefficient of carbonation rate is proportional to the ambient temperature. In the range from 10°C to 50°C, the coefficient increases as the temperature increases. The concrete with a smaller water-cement ratio is more sensitive to the temperature.

5. Comparison and Verification

The simulation results of concrete carbonation of C25, C30, and C35 three kinds of concrete are compared with that of the tested data, as shown in Figure 13. The tested results are also compared with the predicted value (Xq) obtained by Fick’s second law, as presented in Figure 14.

Figures 14 and 15 show the carbonation depth value at the corner (Xj) and other positions (Xq), respectively, and the simulation results are basically in good agreement with the tested results. The maximum error is about 5%. This demonstrates that the developed numerical simulation method for predicting carbonation depth in this paper is valid for concrete structures under the atmospheric environment. Moreover, as compared with the predicted results obtained based on Fick’s second law, the simulation results of the carbonation depth in this paper are more close to the tested results, which proves the correctness and reliability of the proposed method.

In order to verify the practicability of the finite element model in engineering, some carbonation data under the atmospheric environment from several references are used to compare with the finite element results. The environmental factors are listed in Table 5. The comparison results are shown in Table 6.

It can be found from Table 6 that the error between the finite element value and practical measuring value is from −4.8% to 13.4%, and the mean absolute percent error (MAPE) is 5.73%. At the same time, the absolute error of the finite element value is lower than the reference values. In addition, the MAPE between reference values and practical measuring values is 7.67%, which is larger than finite element value. So the finite element model in this paper can more precisely predict the carbonation progress of concrete under the atmospheric environment compared with references.

6. Conclusions

This paper firstly builds a 2-D finite element model to predict the carbonation depth considering the time (t) and temperature (T) dependent carbonation coefficients and carbonation position difference for concrete structures. A two-dimensional mass transfer equation on concrete cross section is proposed, which is solved by the finite element method. The time dependence of CO2 mass concentration in different positions of the concrete block is obtained. Secondly, a test program to investigate the diffusion process of CO2 in concrete is carried out. Finally, the simulation results are compared with the tested results. Several conclusions are obtained as follows:(1)The tested results show that the influence factors such as temperature, time, and water-cement ratio have an obvious effect on the concrete carbonation depth. It was also found from the tested results that under the same humidity and CO2 concentration, the carbonation depth increased about 1.9 times as the environmental temperature increased from 20°C to 40°C.(2)The carbonation depth curves obtained from the finite element method are consistent with the experimental results and the practical measuring value in engineering. This demonstrates that the numerical method can be applied to predict concrete carbonation of real concrete structures with a satisfactory accuracy.(3)The tested results and numerical results quantify the ratio of the carbonation depth value at the corner (Xj) to that at the other positions (Xq), which is from 1.3 to 1.4. It is necessary to consider the time dependence of CO2 concentration and the concentration differences at different structural positions when predicting the carbonation depth of concrete structures.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The financial supports of the National Basic Research Program of China (973 Program, Grant no. 2015CB057701), the National Natural Science Foundation of China (Grant no. 51378081), the Natural Science Foundation of Hunan Province (2018JJ2438), and the Natural Science Foundation of Guangdong Province (2015A0310141) are gratefully acknowledged.