Abstract

In this study, the porosities of C60 high-strength concrete after 0, 30, 60, and 90 freeze-thaw cycles determined via the water retention method are 1.30%, 3.65%, 5.14%, and 7.34%, respectively. Furthermore, a mathematical model of porosity varying with the number of freeze-thaw cycles is established. Using an artificial environment simulation experimental system and the natural diffusion method, the chloride diffusion law of C60 high-strength concrete after 0, 30, 60, and 90 freeze-thaw cycles is obtained. The corresponding diffusion coefficients are calculated based on the experimental results and Fick’s law, where 0.3431 × 10−12, 0.5288 × 10−12, and 0.6712 × 10−12, and 0.8930 × 10−12 m2/s are obtained, respectively, and a mathematical model of diffusion coefficient with freeze-thawing is established. Transport control equations comprising solution flow and solute migration control equations are established for chloride ions in concrete after freeze-thawing cycles. The equations consider the effects of freeze-thawing, solution pressure, solution concentration, solution density, convection, mechanical dispersion, and chemisorption on chloride ion transport in concrete. Using COMSOL numerical software, the transport control equations for chloride ions are solved using a real concrete numerical model, and the chloride ion corrosion process in concrete after freeze-thaw cycles is simulated. The simulation results are consistent with the experimental values.

1. Introduction

With the continuous increase in energy demand, mineral resources extraction has gradually shifted to deeper strata, and the development depth of vertical shafts has increased consequently, along with the increase in the cross-sectional size of the shaft [112]. For the development of deep and thick unstable aquifer mines, the artificial ground freezing method is widely used [1315]. When constructing using the freezing method, the freezing shaft lining is crucial, and its main function is to block groundwater and withstand the temporary load to ensure the safety and smooth operation of shaft sinking construction [1619]. The structural forms of freezing shaft lining primarily include single-layer reinforced concrete, double-layer reinforced concrete, and reinforced concrete composite shaft lining. Currently, concrete materials are essential in the construction of all freezing shaft lining. With the increase in mining depth, high-strength and high-performance concrete materials such as C60 and C70 have been used for freezing shaft lining, and many scholars have conducted relevant studies [2024].

As the freezing shaft lining is generally cast on site, pouring concrete will produce a large amount of heat of hydration, and the heat of hydration will contribute to part of the freezing lining temperature. After the heat of hydration dissipates, the freezing lining will refreeze, and the temperature will decrease; therefore, the freezing shaft lining will be in a freeze-thaw environment with a changing temperature field, posing a significant challenge to the strength and stability of the freezing lining. Hence, many scholars have investigated the effect of temperature on the rupture of vertical shaft lining [2032]. Reinforced concrete structures for freezing shaft lining not only are affected by temperature and stress fields but also are susceptible to the corrosive effects of the subsurface environment. For example, sulfate has a certain corrosive and destructive effect on concrete materials, and the intrusion of chloride ions accelerates the rusting of reinforcing steel bars, thereby affecting the stability of the structure; this is closely associated with the durability problem of concrete structures and materials [3337].

By understanding the corrosion law of chloride ions in C60 high-strength concrete after freeze-thawing, we aim to establish control equations for chloride ion transport in concrete considering freeze-thawing, solution pressure, solution concentration, solution density, convection, mechanical dispersion, and chemisorption and verify them by numerical calculations to provide some theoretical basis for the design, construction, and safe use evaluation of freezing shaft lining.

2. Porosity of Concrete after Freeze-Thawing Cycles

2.1. Materials and Mix Proportions

They are cement: 42.5 grade ordinary silicate cement supplied by Xuzhou Zhonglian Cement Plant, the composition of which is shown in Table 1; fine aggregate: river sand with a fineness modulus of 2.8, good gradation, apparent density of 2679 kg/m3, and silt content of 2.49%; and coarse aggregate: stone, particle size of 5–20 mm, apparent density of 2719 kg/m3, and silt content of 0.54%. The mix proportions of concrete were 453 kg/m3 cement, 740 kg/m3 sand, 1112 kg/m3 stone, and 145 kg/m3 water; the water-cement ratio was 0.32. The compressive strength of the sample was 57.0 MPa after 28 d of standard curing.

2.2. Test Procedure

The test procedure, as shown in Figure 1, is as follows.

Three cylindrical concrete samples of diameter d 100 mm and height h 50 mm (standard curing for 28 d) were used. Their actual height was measured, and then they were subjected to 0, 30, 60, and 90 freeze-thaw cycles (based on the GBT50082-2009 standard, China, shown in Figure 1(a)) and porosity tests (shown in Figures 1(a)–1(d)).

The porosity test was performed using the water saturation method. The freeze-thaw treated samples were first dried at 40°C for 48 h (shown in Figure 1(b)) and then weighed using an electronic balance, placed in a water tank with a 16 mm diameter PVC pipe underneath, and filled with room-temperature water such that the water surface was 4-5 mm above the bottom of the sample (shown in Figure 1(c)). The samples were removed at 5, 10, 30, 60, 120, 720, and 1440 min, respectively, and then weighed (shown in Figure 1(d)) and placed back immediately until all intervals were completed. Finally, the samples were completely submerged, and their water-saturated mass was measured after another 48 h.

2.3. Porosity

Porosity is calculated using the following equation:where is the porosity, is the volume of water in the sample when saturated, and is the actual volume of the sample.

and are calculated using the two following equations, respectively:

In equation (2), is the mass of water in the sample when saturated, and is the density of water; in equation (3), the diameter of the sample  = 100 mm, and is the actual height of the sample.

is calculated using the following equation:where is the mass of the concrete sample when saturated and is the initial mass of the concrete sample after drying.

The porosity of concrete after 0, 30, 60, and 90 freeze-thaw cycles was calculated to be 1.30%, 3.65%, 5.14%, and 7.34%, respectively, and a model of porosity varying with the number of freeze-thaw cycles n was obtained, as shown in the following equation:where is the porosity of concrete after freeze-thawing; is the initial porosity of concrete before freeze-thawing; and are the material parameters, set at 1.61459 and 0.0143, respectively; is the number of freeze-thaw cycles. The experimental values were compared with the calculated values, as shown in Figure 2.

3. Corrosion Pattern of Chloride Ions in Concrete after Freeze-Thawing Cycles

3.1. Test Methodology

The sample measured 100 mm × 100 mm × 400 mm; the materials and mix proportions are shown in Section 2.1, and the test procedure is shown in Figure 3. In the test, four groups (three samples in each group) were formed and treated with 0, 30, 60, and 90 freeze-thaw cycles, respectively (based on GBT50082-2009, shown in Figure 3(a)); subsequently, a 100 mm × 400 mm rectangular surface of each sample was used as the corrosive surface, and the other surfaces were sealed with epoxy resin (shown in Figure 3(b)). The sealed samples were placed into an artificial marine climate simulation laboratory for salt spray corrosion test (shown in Figure 3(c)), and the humidity in the environmental chamber was controlled at 75%–80%. The salt spray corrosion settings were as follows: salt spray “on” for 12 h/d (hour/day) + salt spray “off” for 12 h/d and NaCl mass concentration of 5% in the salt spray solution. When the corrosion time reached 30, 110, 190, and 270 d, the samples were removed, concrete powder was obtained by drilling to depths of 5, 10, 15, 20, 25, 30, 35, and 40 mm, respectively (shown in Figure 3(d)), and the chloride ion concentration values were measured using a DY-2501A chloride ion concentration rapid tester (shown in Figure 3(e)). After each drilling of powder, the borehole was sealed with a sealant. Subsequently, the samples were placed back into the artificial marine climate simulation laboratory and continued to be corroded until the next testing time point (shown in Figure 3(f)), that is, until the end of 270 d of corrosion.

3.2. Chlorine Ion Concentration

Figures 4(a)4(d) present the variation patterns of chloride ion concentration value c with depth l after 30, 110, 190, and 270 d of chloride ion corrosion, respectively, the analysis of which shows that c of the samples decreased with increasing l. In this experiment, the salt spray cycle corrosion mechanism was used, and the driving force for chloride ion intrusion into concrete was primarily provided by capillary action and concentration gradient. The chloride ion concentration gradient in the surface layer of concrete was relatively large, whereas the chloride ion concentration gradient in the deep layer was relatively small; hence, more chloride ions were accumulated in the surface layer of concrete. As the corrosion time progressed, the total chloride ion corrosion and corrosion depth increased. As the number of freeze-thaw cycles increased, the porosity of the concrete increased, the water and air permeability of the samples increased, chloride ions were more easily accessible and stored, and the corrosion concentration of chloride ions increased (at the same depth).

3.3. Chloride Ion Diffusion Coefficient

Fick’s second law was used to describe the corrosion pattern of chloride ions in concrete, as shown in the following equation:

C0 is the initial chloride ion mass concentration at any section of concrete, Cs is the chloride ion mass concentration at the concrete surface at any moment t, D is the chloride ion diffusion coefficient, t is the time, l is the distance from the concrete surface, and erf is the error function.

The inverse solution of equation (6) yields the equation for the chloride ion diffusion coefficient as follows:

The mass concentration of chloride ions on the concrete surface Cs is calculated as follows:where is the porosity of concrete; is the density of sodium chloride solution, that is, 1034 kg/m3; is the density of concrete, that is, 2400 kg/m3; is the mass concentration of sodium chloride in the solution, that is, 5%.

Based on equations (7) and (8), the chloride ion diffusion coefficient D at l = 5 mm was calculated using MATLAB, and the results are shown in Table 2.

A mathematical model correlating the chloride ion diffusion coefficient Dn with the number of freeze-thaw cycles n was developed, as shown in the following equation:where is the diffusion coefficient of concrete after freeze-thawing; is the diffusion coefficient of concrete before freeze-thawing; and are the material parameters, which were set to 1.02151 and 0.00916, respectively; is the number of freeze-thaw cycles. The experimental and calculated values were compared, as shown in Figure 5.

4. Control Equations for Chloride Ion Transport in Concrete after Freeze-Thawing Cycles

4.1. Basic Hypotheses

To derive the control equations for chloride ion transport in freeze-thaw corroded concrete, the following hypotheses are proposed: the liquid in concrete is an ideal liquid; the effect of temperature on liquid flow is not considered; and only a single solute (chloride ion) is considered. The transport of chloride ions in concrete can be classified into two main processes: solution flow and solute migration, in which the effects of freeze-thawing on both primarily affect the porosity and diffusion coefficient.

4.2. Control Equation for Solution Flow
4.2.1. Mass Continuity Equation

where is the density of the solution, is the porosity, t is the time, , and .

4.2.2. Porosity

Introducing the parameter , pore deformation due to pore pressure changes can be characterized.where denotes the pore volume.

The volume of each component phase of the porous medium is expressed as follows:where is the total volume of the porous medium, is the volume of the solid skeleton of the porous medium, and is the porosity of the porous medium. Based on the above hypothesis that the solid particles do not undergo deformation, is a constant; that is, .

Subsequently, equation (12) was used to obtain a model of porosity varying with pore pressure as follows:where and is as shown in equation (5).

4.2.3. Solution Flow Rate

It is assumed that the solution flow in the porous medium conforms to Darcy’s law and that the effect of gravity is not considered.where is the permeability of the solution and is the dynamic viscosity of the solution.

The relationship between permeability and porosity is described using the Kotyakhov [38] model, as shown in the following equation:where is the effective diameter of the porous medium particles. Based on equation (15), the following can be derived:where is the initial permeability of the solution. Because and , equation (16) can be simplified as follows:

Therefore, based on equations (13) and (17), the solution velocity can be expressed as

4.2.4. Solution Density

As chloride ions are transported through concrete, the pore pressure and solution concentration change, resulting in changes in the solution density. Hence, the density of solution is a function of the pore pressure P and concentration c; that is,

Taking the derivative of both sides of equation (19) yields

To effectively characterize the density as a function of pore pressure and solution concentration, two parameters and were introduced, which are the pressure compression coefficient and concentration compression coefficient, respectively.

Parameters and represent the ratio of the change in density due to a unit change in pressure to the initial density and the ratio of the change in density due to a unit change in concentration to the initial density, respectively.

Substituting equation (21) into equation (20) yields

Integrating both sides of equation (22) yieldswhere denotes the initial density, denotes the initial pore pressure, and denotes the initial concentration.

4.2.5. Control Equation for Solution Flow

Using equations (10), (13), (18), and (23), the flow control equation for chloride ion solution in concrete porous medium was obtained as follows:

4.3. Control Equation for Solute Migration

Chloride migration in concrete is governed by convection, molecular diffusion, mechanical dispersion, precipitation dissolution, complexation, and chemisorption. However, only the effects of convection, mechanical dispersion, and chemisorption on chloride migration were considered in this study.

4.3.1. Convective Effects

The convective migration of chloride ions is described by convective flux, as shown in the following equation: where is the convective flow vector; is the solution flow velocity vector.

4.3.2. Mechanical Dispersion

The mechanical dispersion migration flux of chloride ions is expressed by Fick’s first law, as shown in the following equation:where is the dispersion flux vector; represents the mechanical dispersion coefficients in the x-, y-, and z-directions, respectively. If we assume that the porous medium is isotropic, then and ; therefore, .

4.3.3. Chemisorption

The Langmuir adsorption theory was used to describe the chemisorption of chloride ions in concrete transport, as shown in the following equation:where is the dry density of the skeleton material, is the isothermal adsorption concentration, is the Langmuir isothermal nonlinear adsorption partition coefficient, and is the maximum adsorption capacity.

4.3.4. Control Equation for Solute Migration

Equation (28) is the chloride ion convection-diffusion equation, where [D] denotes the chloride ion diffusion coefficient tensor. It is assumed that the diffusion coefficient is uniform everywhere in the concrete; therefore, the diffusion coefficient can be expressed as shown in equation (9). Considering adsorption, the mass of chloride ions in the representative elementary volume is expressed as , and the convection-diffusion equation for chloride ions in concrete considering adsorption can be obtained by replacing in equation (28) with the above equation; that is,

Substituting equation (13) into equation (29) yields the convection-diffusion equation for chloride ions expressed in terms of pressure P and concentration c, as in the following equation:

4.4. Control Equations for Chloride Ion Transport in Concrete

The control equations for chloride ion transport in concrete considering the freeze-thaw action comprise two parts, that is, the solution flow control equation and the solute migration control equation, as shown in 31a and 31b, respectively, which contain the effects of convection, diffusion, and chemisorption, reflecting the coupling between the seepage and solute migration fields, where and D are shown in equations (13) and (9), respectively.

The initial and boundary conditions are solved as follows.

The solution flow and solute migration can be described using the Dirichlet and Neumann boundary conditions as follows:where and are the fluid pressure and fluid flow velocity on the solution domain boundary and are the chloride ion concentration and chloride ion diffusion flux on the solution domain boundary , respectively. The initial conditions for fluid flow and chloride ion diffusion can be expressed aswhere and are the initial fluid pressure and chloride ion concentration distribution over the solution domain , respectively.

5. Model Calculation and Validation

5.1. Calculation Scheme
5.1.1. Numerical Model

Equations (31a) and (31b) were solved using COMSOL numerical calculation software to simulate the chloride ion corrosion in concrete after freeze-thawing. The numerical model measured 100 mm × 100 mm, the left boundary was set as the chloride ion corrosion boundary, and the other boundaries were no-flow boundaries. Point A is a point on the midline at a distance 5 mm from the left boundary, as shown in Figure 6. To study the effect of freeze-thaw concrete chloride ion corrosion, numerical simulation scenarios involving 0, 30, 60, and 90 freeze-thaw cycles and 5% chloride salt concentration were investigated. Because of the complex coupled model and boundary conditions, the solution was complex; therefore, the model parameters were partially simplified in the calculation; for example, the porosity ϕ, solution density ρw, chloride ion diffusion coefficient D, and other parameters were set as constants in the calculation.

5.1.2. Initial Pore Size Distribution

As shown in Figure 7, the processed numerical plots can accurately reflect the pore-fissure distribution in concrete, with smaller porosity corresponding to coarse aggregates and larger porosity corresponding to cracks. When the numbers of freeze-thaw cycles were 0, 30, and 60, no cracks were observed in the concrete; therefore, the pores were relatively dispersed, and large pores dominated the mortar, as shown in Figures 7(a)7(c); when the number of freeze-thaw cycles was 90, wider diagonal cracks appeared in the concrete, and the distribution of cracks can be observed well, as shown in Figure 7(d).

5.2. Calculation Results and Analysis
5.2.1. Concentration Distribution

Figure 8 shows the cloud plots of the chloride ion concentration distribution after 30, 110, 190, and 270 d of chloride ion corrosion of concrete at a corrosive surface concentration of 854 mol/m3 (5% NaCl solution) when the numbers of freeze-thaw cycles n were 0, 30, 60, and 90, respectively. As shown, the chloride ions eroded nonuniformly from the surface to the interior; with the increase in the number of freeze-thaw cycles, chloride ion corrosion accelerated significantly, primarily due to the increased porosity of the concrete. Meanwhile, the chloride ion diffusion and permeation of the concrete after 90 freeze-thaw cycles were significantly faster than those of the concrete after 60 freeze-thaw cycles, primarily due to the increased porosity of the concrete and the generation of microcracks.

5.2.2. Chloride Ion Diffusion

Figure 9 shows the numerical simulation curves of chloride ion corrosion concentration c in concrete with depth l (horizontal model midline) after different numbers of freeze-thaw cycles compared with the experimental values, where the scatter points represent the experimental values, and the curves represent the simulated value. As shown, the simulated value of chloride ion concentration with depth was not a smooth decreasing curve owing to the introduction of the real concrete numerical model. Furthermore, the results reflected the effect of the coarse aggregates, as shown in Figure 9(d). The simulated curves deviated slightly from the experimental values, particularly in the smaller numbers of freeze-thaw cycles, as shown in Figures 9(a) and 9(b).

6. Conclusion

The porosities of C60 high-strength concrete after 0, 30, 60, and 90 freeze-thaw cycles determined using the water retention method were 1.30%, 3.65%, 5.14%, and 7.34%, respectively, and a mathematical model of porosity varying with the number of freeze-thaw cycles was developed.

The chloride ion diffusion patterns of C60 high-strength concrete after 0, 30, 60, and 90 freeze-thaw cycles were obtained using an artificial environment simulation experimental system and the natural diffusion method, and the corresponding diffusion coefficients were calculated to be 0.3431 × 10−12, 0.5288 × 10−12, 0.6712 × 10−12, and 0.8930 × 10−12 m2/s, respectively. Furthermore, a mathematical model of diffusion coefficients varying with the number of freeze-thaw cycles was developed.

Transport control equations for chloride ions in concrete after freeze-thawing comprising a solution flow control equation and a solute migration control equation were developed, in which the effects of freeze-thawing, solution pressure, solution concentration, solution density, convective action, mechanical dispersion, and chemisorption on the transport of chloride ions in concrete were considered.

The chloride ion transport control equation was solved using a real concrete numerical calculation model in COMSOL numerical software to simulate the chloride ion corrosion process after different numbers of freeze-thaw cycles, and the simulated values of chloride ion concentration agreed well with the experimental values.

Data Availability

Some or all data, models, or code generated or used during the study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The authors gratefully acknowledge the financial support for this study provided by the and 51974296), research project of “333 project” in Jiangsu Province of China (BRA2019236), Natural Science Research Project of Higher Education of Jiangsu (General Programs 16KJB560017, 17KJA560004, and 18KJB560019), Research Program of Xuzhou Institute of Technology (XKY2018131), Xuzhou Key Research and Development Program (KC18090 and KC20176), and Science and Technology Project of Xuzhou City (KC18241).