Abstract

The objective of this study is to examine the compressive strength and durability of very high-volume SCM concrete. The prepared 36 concrete specimens were classified into two groups according to their designed 28-day compressive strength. For the high-volume SCM, the FA level was fixed at a weight ratio of 0.4 and the GGBS level varied between the weight ratio of 0.3 and 0.5, which resulted in 70–90% replacement of OPC. To enhance the compressive strength of very high-volume SCM concrete at an early age, the unit water content was controlled to be less than 150 kg/m3, and a specially modified polycarboxylate-based water-reducing agent was added. Test results showed that as SCM ratio ( ) increased, the strength gain ratio at an early age relative to the 28-day strength tended to decrease, whereas that at a long-term age increased up to of 0.8, beyond which it decreased. In addition, the beneficial effect of SCMs on the freezing-and-thawing and chloride resistances of the concrete decreased at of 0.9. Hence, it is recommended that needs to be restricted to less than 0.8–0.85 in order to obtain a consistent positive influence on the compressive strength and durability of SCM concrete.

1. Introduction

Ordinary Portland cement (OPC), an essential construction material, has contributed substantially to building and infrastructure development. However, since the late 1990s the concrete industries have exerted considerable effort and made investments to minimize the use of OPC, partly because of serious worldwide issue to reduce greenhouse gas emissions. It is generally estimated that the production of one ton of OPC consumes approximately 2.8 tons of raw materials such as limestone and coal and that it releases about 0.7–0.95 tons of carbon dioxide (CO2) into the Earth’s atmosphere from the decarbonation of lime in the kiln and the combustion of fuels [1, 2]. Because of the high CO2 inventory of OPC, the annual emission of greenhouse gases from the worldwide production of OPC is estimated to be approximately 1.35 billion tons [3]. Furthermore, the average electricity consumption in cement manufacturing is given as 106 kWh/ton, which is equivalent to approximately 1.2 GJ/ton in primary energy [3]. For these reasons, a stronger effort is required for the development of an alternative practical concrete technology that ensures low CO2 emissions.

The use of high-volume supplementary cementitious materials (SCMs) as partial replacement for OPC in concrete has become increasingly attractive for the development of sustainable construction materials with low CO2 emissions. As a result, the practical application of by-products such as fly ash (FA) and ground granulated blast-furnace slag (GGBS) as SCMs has gradually increased in the construction industry because of their environmentally beneficial recycling effect and remarkably low CO2 inventory [4]. Furthermore, the appropriate addition of SCMs in place of OPC can improve concrete properties as follows [58]. The pozzolanic activity of SCMs is effective for forming a denser matrix, leading to higher strength (especially at a long-term age) and better durability of the concrete; namely, the pozzolanic activity improves the impermeability of the concrete through the formation of calcium silicate hydrate (CSH) and calcium aluminate hydrate (CAH) gels. FA with spherically shaped particles improves the workability of fresh concrete, which reduces the demand for water for targeted workability and leads to reduced bleeding and less shrinkage deformation of the concrete. The temperature increase during cement hydration is controlled, which helps reduce cracking in mass concrete at early ages. To attain these positive effects, the typical individual limitation for OPC replacement is commonly estimated to be 15–20% for FA and 40–50% for GGBS [9]. On the other hand, it has been commonly pointed out [10, 11] that a large amount of SCMs is not helpful in improving the workability of concrete to any considerable extent because of their low density. Furthermore, the longer curing time owing to their slow pozzolanic reaction which converts soluble alkali into a more stable CSH gel requires longer curing time needed to gain targeted concrete strength. This indicates that a relatively higher strength gain at an early age is one of the essential considerations for the practical use of high-volume SCM concrete.

Malhotra et al. [1012] did pioneering work on high-volume FA concrete and conducted extensive studies to establish and improve the characteristics of concrete containing large amounts of SCMs. Mahmoud et al. [13] showed that concrete mixes made with a ternary binder that incorporated both FA and GGBS have an advantage in terms of early strength development over concretes with FA alone. Huang et al. [14] confirmed the feasibility of using up to 80% Class F of FA as an OPC replacement in concrete if rational mixture proportions are provided. Chen et al. [6] proposed that the amount of cement paste and the water content need to be minimized in order to obtain good quality concrete containing a high volume of FA and GGBS. Lee and Wu [15] reported that FA with a high loss-on-ignition (LOI) value has an adverse influence on the strength and durability of concrete. Yazıcı [16] demonstrated that the chloride-ion penetration depth of concrete decreased with the increase in substitution level of FA up to 30%, beyond which it is marginally affected by the FA content. Overall, from a review of recent experimental observations, it can be concluded that the extent of improvement of the strength and durability of high-volume SCM concrete depends on the mixture proportions of each ingredient for concrete and the chemical composition and the quality of the SCMs. Moreover, the optimization of high-volume SCMs needs to be qualified for required specification in the intended application of concrete.

The objective of the present study is to examine the practical feasibility of producing very-high-volume SCM concrete (incorporating FA and GGBS) with relatively good strength gain at an early age. A total of 36 high-volume SCM concrete mixes with different mixture proportions were prepared according to designed concrete compressive strengths of 24 MPa (Group I) and 30 MPa (Group II). As a partial replacement for OPC, the weight ratio of FA was fixed at 0.4, whereas that of GGBS varied between 0.3 and 0.5; as a result, 70–90% of the OPC was replaced with FA and GGBS. Concrete mixes with of 0.25 and of 0.15 were also prepared as control specimens in each group. To achieve good strength gain, especially at an early age, the unit water content was controlled to be less than 150  , and a polycarboxylate-based water-reducing agent was added after being specially modified through the adjustment of the amount of polyethylene glycol alkyl ether and the addition of an amine. Simple equations to predict strength development of the very-high-volume SCM concrete samples are proposed based on the nonlinear multiple regression analysis of the measured results. Four very-high-volume SCM specimens, together with the companion control mixes, were selected in order to examine their durability under the following environments: repeated freezing and thawing, chloride penetration, and sulfate attack.

2. Experimental Program

2.1. Materials

OPC (ASTM Type I) was partially replaced with commercially available FA and GGBS powders, which produces a ternary-type binder. The chemical compositions of these materials were determined by X-ray fluorescence (XRF) analysis and the results are given in Table 1. The FA had low calcium oxide (CaO) and a silicon oxide (SiO2)-to-aluminum oxide (Al2O3) ratio by mass of 1.91, which belongs to Class F of ASTM C618 [17]. The LOI and 28-day activity coefficient of FA were 0.89% and 92%, respectively. GGBS conforming to ASTM C989 [17] had high CaO and a SiO2-to-Al2O3 ratio by mass of 2.29, which is very similar to that of OPC. The basicity of GGBS calculated from the chemical composition was 1.94. The specific gravity and specific surface area, respectively, were 3.15 and 3466 cm2/g for OPC, 2.23 and 3720 cm2/g for FA, and 2.91 and 4497 cm2/g for GGBS.

Locally available natural sand with a maximum particle size of 5 mm and crushed granite with a maximum particle size of 25 mm were used for fine aggregates and coarse aggregates, respectively. The specific gravity and water absorption were 2.61 and 1.16%, respectively, for fine aggregate and 2.62 and 1.78% for coarse aggregate, as given in Table 2. The moduli of fineness of the fine and coarse aggregates were 2.83 and 7.05, respectively.

To maintain good workability of the concrete at lower unit water content, the molecular structure of a polycarboxylate-based water-reducing agent was specially modified as follows. The degree of polymerization of the main chain in the acryl acid-type polycarboxylate polymer was reduced by a factor of 10 (see Figure 1). The molecular weight of a polyethylene glycol mono-alkyl ether monomer was increased to 2000 in order to increase the length of the graft chain in the polycarboxylate polymer. The decreased length of the main chain and increased length of the graft chain are effective for enhancing the dispersibility of polycarboxylate polymers in cement pastes. Furthermore, to obtain an increase in the strength of the concrete at an early age, an amine was added to the modified polycarboxylate polymer. It is known [18] that the addition of an amine is helpful in catalyzing the hydration reaction of cement at an early age because it accelerates the leaching rate of Ca2+ and OH ions from the mineral compositions of the cement. From previous tests [19], the optimum dosage of the amine was determined to be 3% of the modified polycarboxylate polymer weight.

2.2. Specimens and Mixture Proportions

Table 3 shows the main mixture parameters for concrete specimens using FA and GGBS to achieve the targeted properties. All concrete mixes were classified into two groups according to the designed 28-day compressive strength of 24 MPa (Group I) and 30 MPa (Group II). The selected test parameters in each group were as follows. Two levels of the water content were used, 140  and 150  . The unit binder content for each water content level was varied as 310  , 330  , and 350  for Group I and 370  , 390  , and 410  for Group II; as a result, the water-to-binder ratios for the Group I mixes were calculated to be 45.2%, 42.4%, and 40.0%, respectively, for of 140  and 48.4%, 45.5%, and 42.9% for of 150  , while those in the Group II mixes were 37.8%, 35.9%, and 34.2% for of 140  and 40.5%, 38.5%, and 36.6% for of 150  . SCM level as a partial replacement for OPC was varied as 0.7, 0.8, and 0.9. At each , was fixed to be 0.4, whereas varied as 0.3, 0.4, and 0.5. The addition of FA as a partial replacement of OPC is favorable to the reduction of hydration heat of concrete but unfavorable to the strength development of concrete at an early age. Considering this fact, the present study selected to be 0.4. The volumetric fine aggregate-to-total aggregate ratio was designed to be 48% for of 140  and 46% for of 150  . For comparison, a control mix with a typical ( of 0.25 and of 0.15) was also prepared for each group. Considering the demand increase trend on the use of SCM, the typical SCM concrete was selected for the control mix, instead of OPC concrete. From the practical mixture proportions of ready-mixed concrete batches, the unit water and binder contents determined for the control mixes were 184  and 342  , respectively, for Group I and 165  and 400  for Group II. The targeted air content and initial slump of all concrete mixes were 4.5 ± 1.5% and 210 ± 25 mm, respectively. To meet the designed initial air content and slump , an air entraining agent and the specially modified polycarboxylate-based high-range water-reducing agent were added, as given in Table 4. The state of moisture in aggregates was measured before the mix of concrete, and the surface water on aggregates was then reflected through the correction of the unit water content.

For easy recognition of test parameters, the concrete specimens were notated sequentially using the targeted compressive strength, water content, binder content, and SCM level as a partial replacement for OPC. For example, specimen I-140-310-0.7 indicates a concrete with of 24 MPa produced from the following mixture proportions: of 140  , of 310  , and of 0.7 ( of 0.4 and of 0.3). Concrete specimens denoted by I-C and II-C indicate the control concrete with a typical value in each group.

2.3. Casting, Curing, and Testing

All concrete specimens were mixed using a twin forced mixing-type mixer with 0.35  capacity. The initial slump and air content of fresh concrete were measured in accordance with the ASTM C143 and C231 provisions, respectively [17]. All specimens were cured under water with temperature of 23 ± 2°C until testing at a specified age. All steel molds were removed after aging for 36 h.

The compressive strength of the concrete was measured using cylindrical specimens of 100 mm in diameter and 200 mm high at ages of 3, 7, 28, 56, and 91 days in accordance with ASTM C39 [17]. The durability properties (freezing-and-thawing, chloride ion penetration, and sulfate resistances) were examined for the four selected very-high-volume SCM concrete mixes and two control mixes. All specimens used to measure the durability were demolded at an age of 1 day. The resistance to the freezing-and-thawing cycle of concrete was determined using 100 × 75 × 400 mm prisms in accordance with procedure A specified in ASTM C666 [17]. Prior to the rapid freezing-and-thawing test, the prism specimens were cured for 14 days and saturated in lime water for 48 h. With the start of tests, the relative dynamic modulus of elasticity was recorded at intervals of 30 cycles of freezing-and-thawing up to a maximum of 300 cycles. The resistance to chloride penetration was measured at ages of 28 and 91 days in accordance with a nonsteady-state migration test described in NT Build 492 [20]. Concrete cylinders (100 mm in diameter and 200 mm long) were sawn into disks with 50 mm thick. After vacuum saturation of the cylindrical test specimens in a Ca(OH)2 solution (4 g/L), an external electrical potential was applied axially across the specimen, forcing the chloride ions outside to migrate into the specimen. The catholyte solution was a 10% NaCl solution, whereas the anolyte solution was a 0.3 N NaOH solution. The penetration depth, measured from the visible white silver chloride precipitation at saturation ages of 28 and 91 days, was then converted into the chloride migration coefficient according to the procedure specified in NT Build 492. The sulfate resistance of the concrete was evaluated from the variations of compressive strength of the specimens saturated in a curing tank containing 5% sulfuric acid solution for 28 days.

3. Test Results and Discussion

3.1. Initial Slump and Air Content

The ratios of the modified polycarboxylate-based water-reducing agent and air entraining agent to the total binder by weight used to achieve the target and are given in Table 4. In general, a greater amount of was required for the very-high-volume SCM concrete mixes than for the companion control mixes, regardless of and values. The value of was between 0.028% and 0.042% for Group I mixes and between 0.032% and 0.045% for Group II mixes, indicating that of fresh concrete without the air-entraining agent was commonly lower in Group II mixes than in Group I mixes. To achieve the target compressive strength, a greater amount of binder was needed for the Group II mixes than for the Group I mixes at the same water content. This implies that increasing at the same water content is accompanied by a decrease in the number of macrocapillaries and artificial air pores [9]. The specially modified polycarboxylate-based water-reducing agent was commonly added in the amount of 0.7–1.0% of the binder weight for the concrete mixes tested. The value of added to meet the targeted was slightly higher for the Group II mixes than for the Group I mixes. This is attributed to the fact that of the Group II mixes was lower than that of the Group I mixes. On the other hand, the value of tended to be independent of , indicating that the GGBS content has little influence on the workability of concrete [9].

3.2. Compressive Strength at 28 Days

Most concrete mixes with of 140  met the targeted 28-day compressive strength , as given in Table 4. However, some specimens with of 150  failed to achieve , in particular, for the concrete with of 0.9, and for the Group I concrete with of 48.4% and the Group II concrete with of 40.5%. As expected, the measured 28-day compressive strength ( ) decreased with increasing and . The ratio of of the very-high-volume SCM concretes relative to that of the control concrete is shown in Figure 2. The relative 28-day strength commonly decreased with increasing , indicating that the rate of the decrease was greater for Group II mixes than that for Group I mixes. All concrete mixes with more than 0.8 developed lower than the control concrete. Furthermore, of the concrete with of 150  was commonly lower by approximately 10% than that of the control concrete with of 140  , even at the same , indicating that of very-high-volume SCM concrete is somewhat affected by . Overall, to obtain a value of equivalent to that of a conventional concrete with a typical , very-high-volume SCM concrete should have and .

In general, is taken to be inversely proportional to and [9]. Considering this fact, Yang [21] proposed an empirical model to predict the value of of concrete with various SCMs based on a nonlinear multiple regression (NLMR) analysis using an extensive amount of test data collected from the available literature. In the database for the regression analysis, the primarily ranges of the main parameters are as follows: = 0.25–0.6, = 0.1–0.4, and = 0.2–0.4. The number of ternary-type-binders using OPC, FA, and GGBS in the database is small and is mostly within 0.5. Overall, the following equation, proposed by Yang, is thought to be suitable for concrete with a typical not exceeding 0.5: where (=10 MPa) is the reference value for the 28-day compressive strength of concrete and is the silica fume level as a partial replacement for OPC.

Table 4 clearly shows that of high-volume SCM concrete is somewhat sensitive to , though sensitivity depends on the type and level of SCMs. Furthermore, to obtain the same of OPC concrete or concrete with a typical SCM level, a lower is required for high-volume SCM concrete as compared to OPC concrete or typical SCM concrete. Considering these experimental observations, (1) was modified using the current test data to predict the of high-volume SCM concrete (see Figure 3). Consider where (=100  ) is the reference value for the unit water content.

Comparisons of the measured 28-day compressive strength and predictions obtained from the Yang’s model (1) and the current model (2) are plotted in Figure 4. The current model gives lower values of than the Yang’s model. The mean and standard deviation of the ratios between the experimental results and the predicted results are 0.89 and 0.103, respectively, for the Yang’s model and 0.99 and 0.062 for the current model. This indicates that the Yang’s model based on concrete mixes with typical SCM levels is likely to overestimate the 28-day compressive strength of high-volume SCM concrete.

3.3. Compressive Strength Development

The typical compressive strength development rate of high-volume SCM concrete is shown in Figure 5. On the same figure, predictions determined from the ACI 209 equation [22] are plotted for comparison. It was difficult to determine the effect of on the strength development rate. As increased, the strength gain ratio at an early age relative to the 28-day strength tended to decrease, whereas that at a long-term age increased up to of 0.8, beyond which it decreased somewhat. A slightly higher ratio at an early age and a slightly lower ratio at a long-term age were observed for Group II mixes as compared to Group I mixes, indicating that the strength development rate is affected by . Relative to the 28-day strength of high-volume SCM concrete, the strength gain ratio at an age of 3 days ranged between 0.2 and 0.28 for Group I mixes and between 0.27 and 0.33 for Group II mixes, whereas that at age of 91 days ranged between 1.33 and 1.46 for Group I mixes and between 1.27 and 1.43 for Group II mixes. As compared with the predictions from the ACI equation, those values are lower by approximately 27–50% at 3 days and higher by approximately 14–31% at 91 days. This indicates that, by the ACI 209 equation, the compressive strength of very-high-volume SCM concrete is likely to be slightly overestimated at an early age or, conversely, underestimated at a long-term age. Although the specially modified polycarboxylate-based water-reducing agent was added to enhance the early strength of very high-volume SCM concrete, a strength gain lower than that found using the ACI 209 equation was measured at the ages of 3 and 7 days. However, it can be estimated that these low gains at an early age are not detrimental because the early strength gain of concrete with typical is frequently found to be 10–40% lower than that of OPC concrete or the values predicted using the ACI 209 equation [7, 9]. Hence, the specially modified polycarboxylate-based water-reducing agent is expected to contribute to the early strength gain of very-high-volume SCM concrete.

The ACI 209 provision [22] empirically recommends the following parabolic strength development equation based on test results of OPC concrete: where is the compressive strength according to age (in days). The strength development rate at early and long-term ages is determined by the variation of the constants and . In general, a lower value of leads to a higher compressive strength gain at an early age. For OPC concrete cured by air drying, it is recommended that the values of and are 4.0 and 0.85, respectively. However, these values need to be modified for very-high-volume SCM concrete in order to minimize the error observed in Figure 5. To fit the strength development characteristics of high-volume SCM concrete, the values of both constants were determined using test data (see Table 4). All specimens had a high correlation coefficient of more than 0.93, as listed in Table 4. With the increase of the , , and , the constant tends to increase, whereas decreases. The determined values of the constants appear to be more significantly affected by than by , whereas they are independent of . Based on regression analysis using these influencing parameters, the two constants and in (3) were proposed by the following linear equations (Figure 6):

Comparisons of the measured and predicted compressive strengths at various ages are shown in Figure 7. Note that in (3) is determined using (2). The mean , standard deviation , and coefficient of variation of the ratios between the experimental and predicted results are also given in the same figure. Compressive strengths at different ages predicted using (2)–(5) are mostly within ±12.5% of the measured values, giving values of and that range between 0.952 and 1.059 and between 0.061 and 0.097, respectively. The values of and for all tested ages were calculated to be 1.0 and 0.082, respectively. The proposed equations describe well the compressive strength development of very high-volume SCM concrete according to age.

3.4. Durability

The relative dynamic modulus of elasticity recorded every 30 cycles of freezing-and-thawing is shown in Figure 8. The control mixes maintained values of above 98% throughout the 300 freezing-and-thawing cycles. The high-volume concrete with of 0.8 showed the same behavior as the control mixes. For the very-high-volume concrete with of 0.9, the value of remained at 98% until the 210th freezing-and-thawing cycle, beyond which it gradually decreased to 90% until the end of the tests (300 cycles). This indicates that the freezing-and-thawing resistance of the selected high-volume SCM concrete mixes is comparable to that of the control mixes with the typical .

Figure 9 presents the nonsteady-state chloride migration coefficients of concrete specimens at ages of 28 and 91 days, which were calculated from the measured chloride penetration depth, in accordance with the procedure specified in NT Build 492 [20]. As expected, the concrete with a designed strength of 30 MPa (Group II) had lower values of than that with designed strength of 24 MPa (Group I). Furthermore, the value of tended to decrease with increasing age. The ratios of values between ages of 91 and 28 days were calculated to be 0.81 and 0.61 for the I-C and II-C control specimens, respectively, and 0.44 and 0.33 for the I-140-330-0.8 and II-140-390-0.8 specimens, indicating that the decrease of with age is higher in the very-high-volume SCM concrete than in the control concrete. At an age of 28 days, a slightly higher was calculated for the very-high-volume SCM concrete than for the control concrete, regardless of concrete compressive strength. However, at an age of 91 days, tended to decrease with increasing up to 0.8, beyond which it started to increase.

The variation of compressive strength of concrete measured from specimens saturated in 5% sulfuric acid solution for 28 days is shown in Figure 10(a). The appearance of those specimens is presented in Figure 10(b). The deterioration ratio of owing to the saturation in sulfuric acid solution was between 30 and 32% for the control concrete, whereas it decreased to 8–20% for the very-high-volume SCM concrete; in other words, the ratio of after saturation in a sulfuric acid solution for 28 days relative to the concrete cured at room temperature was measured to be 69%, 81%, and 92% for specimens I-C, I-140-330-0.8, and I-140-350-0.9, respectively. This indicates that the deterioration of owing to sulfate attack decreased with increasing . This trend was similarly observed in terms of the damage to the specimens; that is, the presence of damaged chips and flaws decreased with increasing . Hence, it can be proposed that the very-high-volume SCM concrete has superior sulfate resistance as compared to conventional concrete.

The beneficial effect of SCMs on the durability of concrete can be explained by improvement in both the impermeability and diffusion taking place in water-filled pores or by capillary suction. Gruyaert et al. [23] showed that the value of in concrete mixes with varying from 0 to 0.85 recorded at an age of 91 days decreases with increasing . However, the addition of SCM exceeding a certain limit would result in decreased impermeability, as demonstrated in freezing-and-thawing and chloride resistances. Hence, needs to be restricted to less than 0.8–0.85 in order to maintain a positive influence on the durability of concrete.

4. Conclusions

The present investigation needs to be further extended to examine the carbonation resistance and inelastic deformation of very-high-volume SCM concrete in order to improve the compressive strength. From the experimental observations on the compressive strength and durability in the current study, the following conclusions may be drawn.(1)The compressive strength of the concrete with of 150  was commonly lower than that of the companion concrete with of 140  by approximately 10%, even at the same , showing that of high-volume SCM concrete is somewhat sensitive to .(2)To achieve a value of equivalent to that of conventional concrete with typical , and in very-high-volume SCM concrete need to be restricted to less than 40% and to 0.7, respectively.(3)As increased, the strength gain ratio at an early age relative to the 28-day strength tended to decrease, whereas that at a long-term age increased until reaching of 0.8, beyond which it decreased somewhat.(4)Unlike the ACI 209 equation which overestimates the early strength of high-volume SCM concrete and underestimates the strength at a long-term age, the proposed equations describe well the compressive strength development of very-high-volume SCM concrete; the mean and standard deviations of the ratios between the experimental and predicted results were 1.0 and 0.082, respectively.(5)In general, the freezing-and-thawing, chloride, and sulfate resistances of the high-volume SCM concrete mixes were comparable to those of the control mixes with the typical . However, the beneficial effect of SCMs on the freezing-and-thawing and chloride resistances of concrete decreased at of 0.9.

Notations

: Initial air content of fresh concrete
:Unit binder content
:Nonsteady state chloride migration coefficient
:Relative dynamic modulus of elasticity
:Unit coarse aggregate content
:Measured concrete compressive strength at an age of 28 days
:Reference concrete compressive strength (=10 MPa)
:Concrete compressive strength at age (in days)
:Designed 28-day compressive strength of concrete
:Ratio of air entraining agent to binder by weight
:Ratio of fly ash to binder by weight
:Ratio of granulated ground blast-furnace slag (GGBS) to binder by weight
:Ratio of silica fume (SF) to binder by weight
:Ratio of supplementary cementitious materials (SCMs) to binder by weight
:Ratio of the modified polycarboxylate-based water-reducing agent to binder by weight
:Unit fine aggregate content
:Fine aggregate-to-total aggregate ratio by volume
:Initial slump of fresh concrete
:Unit water content
:Reference value for the unit water content (=100  )
:Water-to-binder ratio
:Mean of the ratios between experiments and predicted compressive strengths
:Standard deviation of
:Coefficient of variation of .

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

This work was supported by the Nuclear Power Core Technology Development Program of the Korea Institute of Energy Technology Evaluation and Planning (KETEP) with a Grant from the Ministry of Trade, Industry & Energy, Republic of Korea (no. 20131520100750).