Abstract
Portland cement is usually substituted with granulated blast furnace slag to make a slag-blended cement. A slag-incorporating cement has a more complicated hydration process than Portland cement due to the interactions between the slag reaction and the hydration of Portland cement in the cementitious systems. Understanding the effect of slag substitution on the hydration of cement is still challenging; hence, this research is aimed at predicting the hydration of slag-incorporating cement. To achieve this, the extended CEMHYD3D model was employed to predict the hydration of a slag-blended cement. The simulation was done with cement paste samples made with various water-to-cement ratios and different slag substitution levels, in which the interaction between the hydration of Portland cement and the reaction of slag was considered. The prediction model has been validated with experimental results and verified to be successful in predicting the hydration of slag-blended cement.
1. Introduction
Supplementary cementitious materials (SCMs) such as ground-granulated blast furnace slag, often called slag cement, fly ash (FA), and silica fume (SF) are widely used in the cementitious composites due to their durability, environmental friendliness, and minimal cost. When added to Portland cement, SCMs improve the characteristics of a cementitious material via either their pozzolanic activity and/or hydraulic activity. They reduce the energy used and the CO2 released while producing Portland cement [1], and they are widely used in cement-based composites. When compared to the composites made only with Portland cement, composites incorporating SCMs have better durability regarding several deterioration processes, such as chloride penetration and sulfate attack [2]. Slag is among the most frequently used SCMs and is made from iron blast-furnace slag during the reduction of iron ore with coke at a temperature of about 1,500°C. Slag shows high activity during hydration and comprises over 90% of glass because of quick water quenching [3]. Studies conducted by Choi et al. [4], Lee et al. [5], and Cahyani and Rusdianto [6] explained that the substitution of slag for Portland cement enhances the workability and hydration of the process.
At the microscale level, cement particles are typically made up of tricalcium silicate (C3S), dicalcium silicate (C2S), tricalcium aluminate (C3A), and tetra-calcium aluminate ferrite (C4AF) mineral phases. Cementitious materials require a supply of water for hydration, and moisture exchange can shift the water balance within the material. The volume of water in the component changes due to these two processes. Hydration is highly sensitive to the water-to-cement (w/c) ratio of the cementitious system immediately after it has been prepared. Once the cement and water are mixed, a number of chemical reactions can be performed: C3S and C2S produce calcium hydroxide (CH) and calcium–silicate–hydrates (C─S─H), while C3A and C4AF in the hydration process contribute to the formation of C3AH6 and FH3 [7].
The substitution of slag for Portland cement can alter the rate of hydration and the microstructural properties of the cementitious materials. The effect of slag on the hydration of cement has been experimentally studied by Kolani et al. [3], Binici et al. [8], and Zhao et al. [9]. In other studies, advanced characterization methods like scanning electron microscopy (SEM) [10] and X-ray diffraction (XRD) [11] were employed to determine the hydration and microstructure of various blended cements. But techniques like SEM and XRD analysis can also be time-consuming and costly. Computer models are helpful to simulate the hydration of cement. As long as water and space are available, the volume of cement and hydration change over time and continue to do so during the hydration process. Computer models simulate these changes using volume stoichiometries based on the material properties. On the basis of hydration models for cementitious systems that use solely Portland cement, a number of computer-based models such as CEMHYD3D [12], HYMOSTRUC3D [13, 14], μic [15], and DuCOM [16] were developed to predict the hydration of cement and how its microstructure changes over time. Several multiscale models for the hydration of Portland cement have been developed and are now in use throughout the last few decades [17–19]. The existing models of hydration were developed to predict the hydration of composites made of only Portland cement. The rate of hydration in a single cement system was predicted previously by Park [20], in which the rate of hydration that results from the free water reduction and reduced interaction area within the hydration products was considered. The hydration of cements incorporating SCMs like slag is more complex than that of cementitious systems involving only Portland cement owing to the mineral reactions involved in SCMs. The mineral reactions during the hydration of cements involving SCMs such as slag, FA, and limestone are considered distinct from those involving regular Portland cement, and some reactions are considered by looking at the concentration of calcium hydroxide and free water.
A general hydration model was developed by Maekawa et al. [16] for concrete containing FA and slag. A hydration model for limestone containing cement was developed and validated by Mohamed et al. [21]. As described by Swaddiwudhipong et al. [22], the study predicted the hydration of cements containing SF in concrete with a computer model. A multiphase hydration model was proposed by Merzouki et al. [19] and Meinhard and Lackner [23]. The heat generation and the various microstructural characteristics of slag-incorporating cement are all dependent on the hydration of the minerals involved in each cement in the composite. However, there are currently no models available to predict the rate of hydration of slag-incorporating cements. Furthermore, the prediction models that were developed in the previous studies did not consider how slag influences the rate of hydration of a cement system, considering the w/c ratio and the slag replacement level in the system. Therefore, there is a need to consistently predict the hydration of slag-incorporating cementitious systems with the computer models.
2. Hydration Modeling
A modeling software called CEMHYD3D was created for the first time by the National Institute of Standards and Technology to describe the hydration and microstructure of cement [12]. In this model, the microstructure is explained as a grid of volume pixels, each of which stands for a different type of hydration phase. CEMHYD3D is an effective tool for simulating cement hydration that can also involve SCMs [24]. Because the model simulates the cement particles based on actual cement pictures and particle size distributions, it can accurately model the compositions of cement particles and distributions of phases. So, the extended CEMHYD3D version 3.0 modeling package is proposed to predict the hydration of slag-incorporating cements in this study.
2.1. Hydration Model for Portland Cement
The hydration of a single cement particle can be calculated using the basic equation expressed in Equation (1) for each chemical composition in the cement [25]. The activation of chemical reaction and the diffusion process during hydration of a single cement are all taken into account in this model:where is the degree of hydration of cement, is the mass ratio of water to cement, is the physically bound water in the C─S─H gel, is the effective surface area of the cement particles that are in contact with water, is the total surface area, is the physically free water in the C─S─H gel, is the radii of an unhydrated cement particle, is the density of cement, and is the density of water, is the reaction coefficient in the dormant period, is the coefficient of the rate of reaction for each mineral, and is the effective diffusion coefficient of water in the hydration products for each mineral constituents.
The reaction coefficient, in Equation (2) is assumed to be proportional to the rate of hydration of cement, whereas the values of coefficients B and C govern the initial rates of shell formation and decay, respectively.
In regard to hydrates, the effective diffusion coefficient of water changes based on the radii of the gel pore in the cementitious system. The following Equation (3) presents an expression for this phenomenon in terms of the rate of hydration.
Moreover, when the hydration of the minerals in each cement continues, the free water in the capillary pores reduces. Therefore, it is important to consider that some amount of water can be trapped in the gel pores and may be used for further hydration at every stage of hydration [20]. Hence, as shown in Equation (4), the amount of water in the capillary pore varies with the rate of hydration.where and represent the initial percentages by mass of cement and water in the mixture.
Hydration modeling covers the chemical effects of the cement minerals in the hydration reactions. According to the theories put forth by Papadakis [26], Lu et al. [27], and Maciel et al. [28], the reactions that the minerals in Portland cement undergo during the hydration stage can be described with Equations (5)–(8) as follows:
Both CH and C─S─H are formed when C3S and C2S are hydrated with water. As long as C3A or C4AF aluminate clinker phases are present, aluminate is absorbed by C─S─H which produces a C─A─S─H product. If no remaining aluminate reactants are left, the C3S and C2S reactions further yield CSH without aluminum substitution. However, using Equations (5)–(8), the mass per volume of the CH and C─S─H products in Portland cement can be calculated.
2.2. Hydration Model for Slag-Incorporating Cement
There are five major hydration products that occur when slag is added to cement: C─S─H, C─S─A─H, hydrotalcite (), ettringite (), and AFm phase () are usually formed [29]. While both slag and Portland cement produce similar hydration products, part of the portlandite produced when cement clinker hydrates is actually consumed when slag is hydrated [30–32]. Nevertheless, a high rate of slag reaction requires a lower w/c ratio [26]. The reaction of slag during hydration is appropriately summarized by the approximate critical values shown in Table 1, as reported by Maekawa et al. [33], based on his study conducted on the blended cements.
In accordance with the stoichiometry and hydration model of the slag reaction [33], the mass of calcium hydroxide, the chemically bound water, and the gel water may all be computed with Equations (9)–(11) in which the hydration of slag-incorporating cements can be described as follows:where , , and are the mass of calcium hydroxide, capillary water, and chemically bound water, respectively. sg is the rate of the slag reaction, P is the mass of the slag in the mix, and are the mass of calcium hydroxide formed in the hydration of cement, and the slag reaction in that order, where as and are the mass of chemically bound water and the gel water in the slag reaction, respectively.
3. Experiments
An OPC as a reference cement confirmed to ASTM C 150 and a slag confirmed to ASTM C 989 standard requirements were used in this study. Before the simulation of the rate of hydration of slag-incorporating cement with the proposed model, the individual cements were first characterized according to their particle size distribution (PSD) and chemical compositions. The average PSD for OPC and slag was measured with a laser diffraction (LD) technique and plotted in Figure 1. The hydration model development needs the consideration of some reaction chemistry, for which standard abbreviations are used throughout the model: C stands for CaO, S for SiO2, A for Al2O3, F for Fe2O3, M for MgO, and H for H2O. The hydration equations of the four cement clinker phases (C3S, C2S, C3A, and C4AF) and slag reactions are considered in the model.

As a reference cement, an OPC with the chemical composition presented in Table 2 is used to determine the influence of the substitution of slag into cement. While simulating the hydration of the blended cement, the w/c ratio and the slag replacement level were considered.
The proposed CEMHYD3D model from open source offers the possibility of simulating the hydration of cement. Moreover, the code in the program allows us to adjust the model for extension as necessary.
To extend the CEMHYD3D model and enhance its predictive capabilities, a critical refinement has been introduced by incorporating the reaction of slag into the program. This significant modification is applied because the rate of hydration in cementitious systems is intricately linked to the mineral composition of each specific type of cement.
The compound composition of the slag used in the model is presented in Table 3. In a slag-blended cement system, the simplest and most direct approach is to assume that the slag reacts with CH in the hydration process and yields a single product with a mixed phase. This process provides a C─S─H product that has a lower ratio of Ca/Si than the typical C─S─H product, amorphous silica (AFm), and some hydrotalcite products [34, 35].
Various samples of cement pastes were prepared with a w/c ratio of 0.3, 0.4, 0.5, and 0.6 at a slag replacement level of 10%, 20%, 30%, and 40% by mass of cement. The pastes were made by mixing the raw materials in a vertical axis pot mixer with water first, then slowly adding cementitious materials while the blades were rotating; the subsequent mixing process included a 3-min low speed mixing, a half-minute pause, and a 2-min high speed mixing. After being mixed, fresh mixtures were cast into plastic tubes of φ15 mm, and the two ends of the tubes were sealed to prevent moisture exchange. The sealed specimens were cured in a saturated condition at a standard temperature of 20°C until the age of testing. The specimens were tested at specific ages of 1, 3, 7, 28, and 90 days. The contents of hydrate water and CH in the Portland cement and blended cement pastes were determined using thermogravimetric analysis (TGA). The specimens were removed from the plastic tubes, ground to a fine powder with particles of size <60 µm in a chamber, and flushed with methanol in a thistle tube connected to a vacuum to arrest the hydration. After drying, about 40 mg of the powder sample was analyzed by measuring the weight while heating the sample up from standard temperature to 550°C in TGA. Weight loss in the ranges of 450–550°C was used to compute the contents of hydrate water and CH, respectively [36–38]. The amount of hydrate water per gram of binder (wH) was calculated by normalizing the hydrate water content by the mass loss of the sample between 450–550°C divided by the mass of the ignited sample, corrected for the measured loss-on-ignition of the unhydrated cement (or of the unhydrated cement divided by the percentage of slag for slag-blended cement). The values of wH were converted to estimated degrees of hydration (α) based on the phase composition of the cement used in the study.
4. Simulation of the Influence of Slag on Hydration of Cement
Cement pastes with various values of w/c and replacement levels of OPC with slag were created using the extended CEMHYD3D version model to match the experimental mix. In all simulations of the proposed cement paste, a conversion factor of 0.00035 hr/cycle2 as described by Bentz [12] has been applied to make conversions between the model hydration cycles and the real time. As discussed in previous sections, it was assumed in this study that the slag in a blended cement reacts with CH and provides a single product of hydration in a mixed phase. With the proposed model, the hydration of slag-incorporating cement was simulated using cement paste samples made with a replacement level of 10%, 20%, 30%, and 40% by mass of slag and a w/c of 0.3, 0.4, 0.5, and 0.6. The prediction has been validated with measured values from experiments based on the standard test method of ASTM C186. Previously, experimental studies were conducted by Zhao et al. [9], Durdziński et al. [39], and Hlobil et al. [40]. In slag-incorporating cement, the amount of CH and the rate of reaction of mineral components determine the rate of hydration of the cement [3, 9]. This proves that the hydration product of slag usually develops in close proximity to the slag particles [41–43]. Using the approach suggested by Saeki and Monteiro [44], the equation for the reaction of slag can be expressed with Equations (12)–(14) as follows, considering the aforementioned factors:where is the mass of calcium hydroxide formed in the hydrated slag-blended mixes obtained from Equation (9), is the initial diffusion coefficient, is the density of slag, is the rate of reaction coefficient of slag in the dormant phase ( and are coefficients), and is the rate of reaction coefficient of slag.
5. Analysis and Discussion
Figures 2–5 shows the hydration of slag-incorporated cements, which were simulated with the extended CEMHYD3D model for different w/c values and various replacement levels. The simulations with the computer model were carried out under the same conditions as the experiments. The model provided a good prediction of the influence of slag on the hydration of cement in all simulations, which agreed with the experimental results.

(a)

(b)

(c)

(d)

(a)

(b)

(c)

(d)

(a)

(b)

(c)

(d)

(a)

(b)

(c)

(d)
The hydration of slag-incorporated cement is simulated in Figure 2 separately as a, b, c, and d in volume fractions of each cement. The model predicted the hydration of blended cement pastes without and with 10%, 20%, 30%, and 40% slag substitution made with a w/c of 0.3 for all mixes under sealed conditions at a temperature of 20°C. It was observed that the cement clinker phase reacted much faster than the reaction of slag in the early ages. Slag can significantly reduce the heat needed for hydration when substituted for cement [45]. However, the slag continues to react at a later age while the clicker phase reaction decreases. In all mixes, the slag reaction increases as the substitution percentage of slag increases, and it was observed that the rate of hydration at a later age increases when a lower w/c of about 0.3 is used in the mix.
As shown in Figure 3 below, the model predicted the influence of slag on the rate of hydration of cement for pastes without and with 10%, 20%, 30%, and 40% slag substitution for a w/c of 0.4 for all mixes under sealed conditions at 20°C. The prediction indicated that no substantial effect of slag on the rate of hydration of cement was observed when up to 10% slag was substituted for a w/c of 0.4. Even though the reaction of clicker phases in the early age is faster than the reaction of slag, the simulation illustrates that the hydration of the blend increases as the substitution level of slag into cement increases in the later age, when a w/c of about 0.4 is used for the paste.
The prediction in a, b, c, and d of Figure 4 illustrates that the rate of reaction of slag at a replacement level of 10%, 20%, 30%, and 40% is similar in most cases with the paste with 0% slag cement when a w/c of 0.5 is utilized for the mix. From the model prediction, no significant effect was observed between the hydration of Portland cement and a slag-incorporating cement in the early ages, when about 10% slag was substituted. In contrast, using a w/c of 0.5 accelerates the rate of hydration in the later age at about 90 days and increases as the substitution percentages of slag in the blend increase.
The predicted values for hydration of a slag-blended cement made with a w/c of 0.6 and varying slag substitution levels are also plotted in a, b, c, and d of Figure 5. The rate of hydration of Portland cement and up to 10% by mass slag substituted cement paste is similar when w/c is around 0.6. When compared with the clicker phase, the reaction of slag slightly increases as the slag substitution level increases, even at the later age of 90 days when the same w/c is used in both mixes.
The hydration of a slag-incorporating cement was accurately predicted using the extended CEMHYD3D model. The hydration of cement pastes made with a w/c ratio of 0.3, 0.4, 0.5, and 0.6 and with a slag replacement level of 10%, 20%, 30%, and 40% was simulated for the age of 90 days cured in a saturated condition at a standard temperature of 20°C. The model prediction and the measured values from experiments were in good agreement in all mixes. The simulation shows that the slag-incorporating cement was predicted to have greater hydration than the Portland cement in the later ages. The rate of hydration of the blend with slag is found to be higher for a w/c of 0.3 and negligible when a w/c of about 0.6 is used. Hence, a greater slag reaction for hydration is preferred at a higher w/c. This makes sense given that there is more space and water in the paste available for reactions and the formation of products. Moreover, there is more water accessible for cement hydration at the paste level. The simulation clearly shows that the rate of reaction of the slag in the blend increases proportionally as the slag substitution level increases. According to Ahmad et al. [46], when the fraction of slag is increased up to 60% in cement, it reduces the amount of time required to reach a maximal rate of hydration. However, no substantial change was observed in the rate of hydration for a slag substitution up to 10% in Portland cement.
6. Conclusions
The paper presented the simulation of the influence of slag substitution on the hydration of cement using the extended CEMHYD3D version 3.0 model. The model predictions agreed with the experimental values regarding the w/c ratio and substitution percentage of slag in Portland cement. The possible influence of slag on the hydration of cement is higher when a lower w/c of about 0.3 is used in the mixture in the later ages as the replacement level of slag increases. The increase in w/c ratio to about 0.6 can accelerate the early hydration of slag-blended cement systems but it has a negligible effect in the later ages as the replacement level of slag increases. Thus, slag substitution in Portland cement is recommended with a higher w/c ratio to achieve an optimal hydration result. The rate of slag reaction in a blended cement system is variable and dependent on the substitution level of the slag. The consideration of the chemical composition of slag improves the predictability of hydration in cement, and the same method was applied in this model. Both in the early age and up to the age of 90 days, the hydration of slag-incorporating cement can be predicted with the proposed model.
Data Availability
The data supporting 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
This study was supported by Ministry of Education of the Federal Democratic Republic of Ethiopia.