Advances in Materials Science and Engineering

Volume 2017, Article ID 4090265, 11 pages

https://doi.org/10.1155/2017/4090265

## Kinetic Hydration Heat Modeling for High-Performance Concrete Containing Limestone Powder

College of Engineering, Department of Architectural Engineering, Kangwon National University, Chuncheon-si 200701, Republic of Korea

Correspondence should be addressed to Xiao-Yong Wang; rk.ca.nowgnak@evarbxw

Received 15 June 2017; Accepted 1 August 2017; Published 5 September 2017

Academic Editor: Doo-Yeol Yoo

Copyright © 2017 Xiao-Yong Wang. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Limestone powder is increasingly used in producing high-performance concrete in the modern concrete industry. Limestone powder blended concrete has many advantages, such as increasing the early-age strength, reducing the setting time, improving the workability, and reducing the heat of hydration. This study presents a kinetic model for modeling the hydration heat of limestone blended concrete. First, an improved hydration model is proposed which considers the dilution effect and nucleation effect due to limestone powder addition. A degree of hydration is calculated using this improved hydration model. Second, hydration heat is calculated using the degree of hydration. The effects of water to binder ratio and limestone replacement ratio on hydration heat are clarified. Third, the temperature history and temperature distribution of hardening limestone blended concrete are calculated by combining hydration model with finite element method. The analysis results generally agree with experimental results of high-performance concrete with various mixing proportions.

#### 1. Introduction

The use of limestone powder blended cement is a common practice in the modern concrete industry. The benefits from technical, economic, and ecological aspects can be achieved by using limestone blended concrete [1]. Technical benefits mainly refer to limestone powder which can increase the early-age performance of concrete. Economic benefits mean obtaining cement with a compressive strength similar to control concrete at low production costs. The ecological aspects are the reduction of greenhouse gas emission by using limestone.

Many studies in experimental or theoretical aspects have been done about early-age properties and durability of limestone blended concrete. Bonavetti et al. [2] found that the addition of limestone can increase the early-age strength of concrete. However, the late-age strength is impaired due to the dilution effect of limestone addition. Mohammadi and South [3] reported that limestone addition can reduce the bleeding of fresh concrete and increase the viscosity and cohesiveness of fresh concrete. Mohammadi and South [4] also reported that the concrete with various limestone contents up to 12% has a similar drying shrinkage and sulfate expansion resistance with control concrete. Chen and Kwan [5] measured heat generation of concrete with different limestone stone replacement ratios and binder contents. They found that the addition of limestone can significantly reduce the heat generation of concrete. Palm et al. [6] found that high-level limestone addition can increase the carbonation depth and chloride migration coefficient of concrete. Based on life cycle assessment, Palm et al. [6] also found high-level limestone additions can reduce the CO_{2} emission about 25% in comparison with average cement with the same performance.

Compared with abundant experimental studies, the theoretical models for limestone blended concrete are relatively limited. Lothenbach et al. [7] proposed a thermodynamic model for limestone blended concrete. The formation of monocarboaluminate and bulk compositions of hydrating cement is calculated by thermodynamic models. Bentz [8] proposed a hydration model which analyzed the dilution effect, nucleation effect, and chemical effect of limestone addition. Similarly, Mohamed et al. [9] also proposed a model to evaluate reaction degree of hydration for concrete with different limestone additions. However, Lothenbach et al. [7], Bentz [8], and Mohamed et al.’s [9] studies mainly focus on the degree of hydration in cement-limestone blends. The hydration heat is scarcely simulated in their studies. Poppe and Schutter [10] and Ye et al. [11] proposed models to analyze heat evolution of limestone blended self-compacting concrete. But their study mainly focuses on the isothermal condition or adiabatic temperature rise. For concrete structures in construction sites, the heat release from hydration and heat transfer to ambient environments occurs simultaneously. Poppe and Schutter [10] and Ye et al.’s [11] studies do not consider semiadiabatic temperature rise in real construction sites.

To overcome the shortcomings in current models [7–11], we proposed an integrated numerical procedure to analyze the temperature history and temperature distribution of hardening limestone blended concrete. The hydration model is combined with finite element method. The heat of hydration is calculated from the degree of hydration. The reduction of hydration heat due to limestone additions is clarified through analysis.

#### 2. Hydration Heat Model of Limestone Powder Blended Concrete

##### 2.1. Hydration Model of Portland Cement

Wang and Lee [12] proposed a kinetic hydration model for Portland cement which takes into account the effects of a water to binder ratio (), compound compositions of cement, fineness of cement, and capillary water contents on the hydration of cement. The kinetic hydration model analyzes the involved kinetic processes of cement hydration, such as initial dormant process, phase boundary reaction-controlled process, and diffusion-controlled process. The kinetic hydration model is valid for concrete with various types of Portland cement, various mixing proportions, and various curing temperatures. The equation for the kinetic hydration model is shown as follows:In (1), is hydration degree of cement, is the hydration rate coefficient in the initial dormant period, is the hydration rate coefficient of phase boundary reaction-controlled process, is the hydration rate coefficient in the diffusion-controlled stage, denotes the amount of capillary water at the exterior of hydration products, denotes the effective contacting surface area between the cement particles and capillary water, and denotes the total surface area if hydration products develop unconstrained [12]. Equation (1) also considers chemical and physical aspects of cement hydration. In (1), is the mass of chemically bound water for one-gram hydrated cement (=0.25), is the mass of physically bound water for one-gram hydrated cement (=0.15), denotes the density of water, denotes the density of the cement, denotes the radius of unhydrated cement particles.

The determinations of reaction coefficients , , and capillary water content are shown in (2)–(4), respectively.Equation (2) can be used to determine reaction coefficient . In (2),* B* is the rate of the initial impermeable layer formation, and* C* is the rate of the initial impermeable layer decay.

Equation (3) can be used to determine reaction coefficient . In (3), is the initial diffusion coefficient, and decreases as cement hydration proceeds.

Equation (4) can be used to determine capillary water content . In (4), is the cement content in mixing proportion, is the water content in the mix proportion, and* r* () is an empirical parameter that considers the accessibility of water into an inner anhydrous part through an outer hard shell of cement particles. For high-strength concrete with low ratio at late ages, has a significant influence on the rate of hydration.

The influences of temperature on reaction coefficients can be described by using Arrhenius’s law [12] as follows:where , , , and denote the activation energies of , , , and , respectively. , , , and denote the values of reaction coefficients , , , and at 293 K, respectively.

Based on the degree of hydration of concrete with various types of Portland cement and various curing temperatures, Wang [13] proposed that the reaction coefficients of hydration model, such as , , , and , can be determined from compound compositions of cement. The temperature sensitivity coefficients can be approximately regarded as constants for different types of cement [13]. These relationships are shown as follows:

Summarily, the kinetic hydration model is composed of four rate determining coefficients, that is, the rate of formation of the initial impermeable layer (*B*), the rate of destruction of initial impermeable layer (*C*), the rate of phase boundary reaction-controlled process (), and the rate of diffusion-controlled process (). By using compound compositions of cement, the reaction coefficients of kinetic hydration model can be determined. Furthermore, the degree of hydration can be calculated by using (1). The proposed hydration model is valid for Portland cement concrete with various materials properties and curing conditions.

##### 2.2. Effect of Limestone Addition on Cement Hydration

Wang [13] reported that the addition of limestone presents dilution effect, nucleation effect, and chemical effect on cement hydration. Dilution effect is when cement is partially replaced by limestone, the content of cement is reduced and water to cement ratio increases correspondingly. Nucleation effect is that limestone can serve as nucleation sites of hydrating cement particles. Hydration of cement can accelerate due to nucleation effect. Chemical effect is the formation of monocarboaluminate due to limestone reaction in preference to a monosulfoaluminate.

On the other hand, because the reactivity of limestone is very weak compared with other supplementary cementitious materials, limestone can be approximately regarded as chemical inert filler [1, 11]. Hence in this study, the chemical effect of limestone is not considered.

The dilution effect of limestone powder can be considered by using (4). Wang [13] and Kishi and Saruul [14] proposed that the nucleation effect of limestone relates to the ratio of surface area of cement particles to that of limestone powder. The nucleation effect of limestone powder is significant in phase boundary reaction-controlled process and diffusion-controlled process. The nucleation effect of limestone powder can be considered as follows:

In (7), is the limestone nucleation effect indicator, is the mass of limestone in concrete mixing proportions, and is the Blaine surface area of limestone powder.

In (8), is the updated phase boundary reaction coefficient, and is enhanced coefficients of . In (9), is the updated diffusion coefficient , and is enhanced coefficients of . Based on analysis shown later (Section 3.1), the values of and are set as 0.6.

Summarily, for a cement-limestone blend, the dilution effect is considered through capillary water concentration. The nucleation effect is considered by nucleation effect indicator which considers binder proportions and surface area of binders. Furthermore, by using updated reaction coefficients, the reaction degree of cement in cement-limestone blends can be determined.

##### 2.3. Temperature History Model of Hardening Concrete

Hydrate heat of hydrating concrete is dependent on both cement content and degree of hydration. The relation heat from hydration of concrete can be determined as follows [15–17]:where is hydration heat and is released hydration heat from a unit mass of cement. can be determined by using compound compositions and hydration heat of individual components of cement [15–17].

For hardening concrete, the temperature distribution is in a dynamic heat balance between the hydration heat generation inside the concrete and heat loss to the ambient. The heat generation comes from hydration reactions of the cement. The temperature distribution of hardening concrete is determined as follows [18]:where is the heat capacity of hydrating concrete and can be calculated as the sum of the individual components of concrete, is the thermal conductivity of concrete, is the concrete temperature, is time, and can be calculated based on the degrees of hydration of the cement (10).

For hardening concrete in construction sites, the boundary condition can be described as follows:where is the coefficient of heat convection between concrete surface and ambient; is the temperature on the concrete surface; and is ambient temperature.

Equation (11) is numerically solved by using finite element method. Three-dimensional (3D) eight-node isoparametric element is used to mesh the geometry model of hardening concrete. The iteration scheme of (11) is shown as follows [15–17]:where is the global mass matrix, is global stiffness matrix, is the global vector of temperature, and is integration parameter. The global stiffness matrix , global mass matrix , and global temperature vector can be obtained by assembling of element mass matrix, element stiffness matrix, and element temperature vector, respectively. Generally, to guarantee the stability of the numerical integration in a time domain, the value of integration parameter should be higher than 0.5. In this paper, according to the Galerkin method, the value of the parameter is used as 2/3 [12].

##### 2.4. Summary of the Proposed Numerical Procedure

The numerical procedure consists of a kinetic hydration model and a finite element model. The kinetic hydration model considers the dilution effect and nucleation effect from limestone additions. The heat of hydration of hydrating concrete is calculated by using the degree of hydration and cement content. The calculation results of the heat of hydration are used as a source term in finite element model. By using Galerkin method, the parabolic partial differential equation about temperature distribution of hardening concrete is solved. Temperature history of semiadiabatic temperature rise is calculated considering both concrete materials properties and ambient conditions. The proposed numerical procedure is valuable for thermal cracking analysis of hardening concrete and construction plan design and materials design of concrete structures.

#### 3. Verification of Proposed Model

##### 3.1. Heat of Hydration

Experimental results about isothermal heat evolution shown in [14] are used to verify the proposed limestone blended hydration model. Kishi and Saruul [14] measured isothermal hydration heat of limestone blended cement paste. The mixing proportions of paste specimens are shown in Table 1. In Kishi and Saruul’s study [14], to demonstrate the nucleation effect of limestone, the water to cement is the same for all specimens. Limestone is added as an additional binder, not replacing partial cement, which is different from the general applications of limestone. The limestone ranges within 10% and 30% of total binder, and the Blaine surface area of limestone ranges from 3500 cm^{2}/g to 8900 cm^{2}/g. The curing temperature is 20°C. The used cement is moderate heat Portland cement. By using Portland cement hydration model shown in Section 2.1, the hydration rate of control cement paste without limestone is calculated and shown in Figure 1(a). The -axis of Figure 1(a) represents the heat release for cement portions in cement-limestone blends.