Advanced Cementitious Materials: Mechanical Behavior, Durability, and Volume StabilityView this Special Issue
Research Article | Open Access
Fractal Modeling of Pore Structure and Ionic Diffusivity for Cement Paste
Pore structure in cement based composites is of paramount importance to ionic diffusivity. In this paper, pore structure in cement paste is modeled by means of the recently proposed solid mass fractal model. Moreover, an enhanced Maxwell homogenization method that incorporates the solid mass fractal model is proposed to determine the associated ionic diffusivity. Experiments are performed to validate the modeling, that is, mercury intrusion porosimetry and rapid chloride migration. Results indicate that modeling agrees well with those obtained from experiments.
Ionic (i.e., chloride ion) diffusivity is one of the critical parameters in service life design of cement based composites for coastal infrastructures . In recent years, much attention has been drawn to the direct assessment of ionic diffusivity on a basis of pore structure which could be obtained from experiment or modeling . However, a practical challenge still exists as a result of the intrinsic complexity of cement based composites. In particular, the pore size distribution may range from nanometer (C-S-H gel pore) to micrometer (capillary pore) of several orders . The multiscale approach is therefore an option, which deals with C-S-H gel and capillary pore in terms of the packing of basic globules and hydration products, respectively [4, 5]. Nevertheless, as it has to assign a large number of parameters and assumptions, the multiscale approach often leads to substantial disadvantages in efficiency and viability. More robust approaches would be beneficial to describe the complex pore structure as well as ionic diffusivity for cement based composites.
The fractal characters have been recognized among natural and artificial porous materials such as limes, soils, rocks, and ceramics [6–10]. In essence, fractal describes a natural phenomenon or object that exhibits similar pattern at different scales . Usually, if one-dimensional length of fractal is magnified, the occupied space of fractal is then magnified by a noninteger power. This power is called the fractal dimension. For a fractal phenomenon or object, the character of similarity at different scales facilitates property characterization.
In an early work, Winslow reported the fractal nature of internal surface for cement paste . Thereafter, the fractal characters in cement paste have been discussed a lot during past decades [13–20]. It was revealed that three types of fractal might be present in cement paste, that is, the pore mass fractal, the pore surface fractal, and the solid mass fractal . Based on the fractal geometry theory, the complex pore structure in cement paste could be modeled via simple iterations, as done in literatures [9, 10, 17, 21]. Moreover, the fractal geometry theory was able to account for various transport properties for porous materials such as permeability [22, 23] and electrical or thermal conductivity [24–28]. In this regard, the fractal modeling shows a great potential to well describe the complex pore structure as well as ionic diffusivity for cement based composites.
In this paper, pore structure in cement paste is modeled by means of the recently proposed solid mass fractal model . Moreover, an enhanced Maxwell homogenization method that incorporates the solid mass fractal model is proposed to determine the associated ionic diffusivity. Experiments are performed to validate the modeling, that is, mercury intrusion porosimetry (MIP) and rapid chloride migration (RCM).
2. Solid Mass Fractal Model
As shown in Figure 1, for the solid mass fractal to model porous medium, the generator refers to two distinct phases, that is, the pore phase and the iterating phase. At the initial step, a blank region with linear length is defined in E-dimension Euclidean space which can further be divided into identical subregions. Let be the number of subregions in each dimension such that . Moreover, let and b, respectively, be the proportion and number of iterating phase in the generator such that . Upon successive iterations, the pattern of generator is repeated for the iterating phase at different scales. The generated pore structure can be analyzed as follows.
Step 1. The number of iterating phase and pore phase with the size reads and , respectively.
Step 2. The number of iterating phase and pore phase with the size reads and , respectively.
Step i. The number of iterating phase and pore phase with the size reads and , respectively.
For the generated porous medium, the solid fraction or density is defined as the ratio of solid volume to the total volume thatwhere is the total volume and is the resultant volume of pore phase. As a result, the solid fraction at Step i of iteration holds that . If a constant is defined, the solid fraction can then be rewritten as follows with respect to and :Equation (2) suggests that the generated porous medium is a solid mass fractal with the dimension of . In a log-log diagram, the linearity shall be detected for the solid fraction against the pore size . In addition, the pore size distribution in terms of versus , that is, the cumulative porosity f that counts from the large to the small pores, can be expressed as follows:
3. Modeling of Pore Structure
The MIP test was performed to validate the solid mass fractal model . Figure 2 shows the pore size distribution of cement paste with water cement (w/c) of 0.4 and age of 28 days, which was cured at the condition of relative humidity and 20°C ± 1°C temperature. For the solid fractal mass model, solid fraction is plotted against pore diameter in the log-log diagram, as shown in Figure 3. It can be noted that the sound piecewise linearity exists. Moreover, the slope of linearity varies within two pore ranges. Let and be the lower and upper limit of pore diameter, and let be the threshold of pore diameter for the two pore ranges. The values of , , and can be determined as μm, μm, and μm. In other words, the two pore ranges, that is, the large and the small pore range, are located as and , that is, μm and μm. The large pore range corresponds to capillary pores, while the small pore range corresponds to C-S-H gel pores.
With the determined pore ranges (, ), the least square method is then applied to fit the associated slopes (, ) as well as the fractal dimensions (, ) since , , and . Meanwhile, it is possible to specify the partial porosity of large and small pore range (, ) with the total porosity f holding that . In rough comparison, the fractal dimension of small pore range is smaller than that of large pore range, that is, . The partial porosity of small pore range is larger than that of large pore range, that is, . The fractal modeling of pore structure in cement paste refers to two specific structures, that is, two fractal generators for the large and the small pore range.
For the large pore range, let two integers and be the step of iterations and the number of subregions such that . The solid mass fractal model requires , , and . In this paper, the symbol denotes the power operator. Without complicated mathematical analysis, it leads to and . The value of is obtained via . The number of iterating phase holds that . The intrinsic porosity of generated pore structure holds that . Details of analysis for the example are shown as follows: The large pore range: , The proportion of iterating phase: The number of iterating phases: The intrinsic porosity:
For the small pore range, let two integers and be the step of iterations and the number of subregions, which holds that . The solid mass fractal model requires , , and . It leads to and . The value of is obtained via . The number of iterating phase holds that . The intrinsic porosity of generated pore structure holds that . Details of analysis for the example are shown as follows: The small pore range: , The proportion of iterating phase: The number of iterating phases: The intrinsic porosity:
The build-up of pore structure in cement paste follows from the large to the small pore range. In particular, after steps of iteration to generate the large pore range, the size and number of iterating phase is and , respectively. Note that during the process to generate the small pore range, not all the iterating phases with size of are to be further iterated but a proportion of γ . In other words, the proportion of iterating phases with size of shall not be performed with further iterations. With respect to the chemical composition of cement paste, such and proportion of iterating phases are introduced as C-S-H gel and solid phases, respectively. The solid phases may be anhydrous cement, calcium hydroxide crystal, Aft, and gypsum. It is necessary to determine to fulfill the build-up of pore structure in cement paste. Herewith, consider the total porosity satisfying . Then, , , , , and lead to . Figure 4 shows one realization of modeled pore structure in cement paste. In particular, it refers to two specific structures, that is, the capillary structure and the C-S-H gel structure. The capillary structure of large pore range consists of three phases, that is, the (capillary) pore phase, the solid phase (anhydrous cement, calcium hydroxide crystal, AFt, and gypsum), and the C-S-H gel phase. The C-S-H gel structure of small pore range consists of two phases, that is, the (gel) pore phase and the solid phase (molecule layers within C-S-H gel). The pore size distribution can be written in terms of the piecewise form: for , , , ; for , , , .In particular, for example, it holds that for μm, , , ; for μm, , , .
4. Enhanced Maxwell Homogenization Method
The calculation of effective properties of heterogeneous materials that contain inclusions of diverse shape and/or properties withstands a long history. Numerous homogenization schemes have been developed so far [30–32]. Among them, the Maxwell homogenization method is probably the oldest and also most famous one . It considers heterogeneous materials with unknown effective ionic diffusivity () consisting of bulk phase with ionic diffusivity () and inclusion phase with ionic diffusivity (). Mathematically, the Maxwell homogenization method can be described as follows:where is the volume fraction of inclusions.
An enhanced Maxwell homogenization method is proposed that incorporates the solid mass fractal model to calculate ionic diffusivity of cement paste. In particular, cement paste is viewed as the composite material consisting of two phases. One refers to the low-diffusivity C-S-H gel structure being bulk. The other refers to the high-diffusivity capillary structure being inclusions. The treatment of two distinct phases was also applied in the general effective media theory [34–36]. For bulk (C-S-H gel structure) and inclusion (capillary structure), the associated ionic diffusivity can be approximated in a power form of intrinsic porosity [37, 38]:where and are constants, that is, ionic diffusivities in water corresponding to C-S-H gel and capillary structures. and are the cementation exponent in Archie’s law for C-S-H gel and capillary structures. Other parameters, that is, , , and are defined in the solid mass fractal model. In addition, the volume fraction of inclusions is derived from
As presented above, the Maxwell homogenization method is enhanced with critical parameters determined from the solid mass fractal model and Archie’s law. Then, the enhanced Maxwell homogenization method can be applied to predict effective ionic diffusivity of cement paste.
5. Rapid Chloride Migration Test
The RCM test was performed to measure ionic (chloride) diffusivity of cement paste [39, 40]. Cement paste samples of 50 mm in thickness and 100 mm in diameter were prepared. Samples were then subjected to the vacuum saturation treatment. The saturated samples were treated to be surface-dry. Three samples of replicates were tested at the same time. Power sources with constant voltage outputs (adjustable in the range of V, accuracy of 0.05 V) were applied. The used catholyte and anolyte were 10% NaCl solution and 0.3 M NaOH solution, respectively. The electrolytes were refreshed after each series of tests. After the migration test, the samples were split and sprayed with a 0.1 M AgNO3 solution in order to determine the penetration depth of chlorides. The non-steady-state migration coefficient , that is, ionic (chloride) diffusivity, was calculated according to the NT Build 492 as follows:where is the non-steady-state migration coefficient, ×10−12 m2/s; is absolute value of the applied voltage, V; is average value of the initial and final temperatures in the anolyte solution, °C; is thickness of the specimen, mm; is average value of the penetration depth, mm; and is the test duration, h.
6. Results and Discussions
Cement pastes with w/c ratios of 0.4, 0.5, and 0.6 were compared, as shown in Figure 5. All samples were cured at the same condition (28 days, 95% ± 10% relative humidity, 20°C ± 1°C temperature). Figure 6 shows the parametric analysis with respect to the solid mass fractal model. Results are summarized in Tables 1 and 2. Table 1 suggests that as w/c ratio increases, the partial porosities increase, while the fractal dimensions decrease. Table 2 suggests that some parameters, that is, , , , and , do not change with w/c ratio, while some parameters, that is, , , and , do. That means the three parameters are critical to characterize the effect of w/c ratio on geometry of capillary and C-S-H gel structure. Figure 7 shows the comparison of modeled and measured pore size distributions of cement pastes.
With the enhanced Maxwell homogenization method, ionic diffusivity of cement paste can be obtained. The parameters listed in Tables 1 and 2 are used, that is, () and (, , , , ). Besides, values of the parameters are suggested in Table 3. In particular, m2/s is empirical. It leads to the ionic (chloride) diffusivity of C-S-H gel around m2/s, which is comparable to the value suggested in . m2/s is taken after . The cementation exponents of and for C-S-H gel and capillary structures are suggested from extensive tests, which differs from the usual range of 1.8–2.0 for consolidated sandstones .
For the RCM test, the measured parameters are listed in Table 4, which leads to an easy calculation of the non-steady-state migration coefficient . As shown in Figure 8, results of the fractal modeling agree well with those obtained from experiment (). In this regard, the present fractal modeling provides an optional description of pore structure as well as ionic diffusivity for cement paste.
7. Concluding Remarks
In this paper, modeling of pore structure as well as ionic diffusivity for cement paste is performed based on the fractal geometry theory. The solid mass fractal model is applied to model pore structure, and an enhanced Maxwell homogenization method that incorporates the solid mass fractal model is proposed to calculate the associated ionic diffusivity. Some general conclusions are drawn as follows:(1)The solid mass fractal model well describes pore structure in cement paste in terms of two pore ranges, that is, the large (capillary) and the small (C-S-H gel) pore range. The fractal generator varies with the pore range.(2)The enhanced Maxwell homogenization method considers cement paste consisting of low-diffusivity C-S-H gel structure and high-diffusivity capillary structure. Archie’s law is applicable to approximate ionic diffusivity for C-S-H gel and capillary structures.
|:||Linear length of modeling space|
|:||Dimension of Euclidean space|
|:||Total number of identical subregions|
|:||Number of identical subregions in one dimension|
|:||Proportion of iterating phase in generator|
|:||Number of iterating phases in generator|
|:||Step of iteration|
|:||Size of iterating phase|
|:||Solid fraction or density|
|:||Total volume of porous medium|
|:||Volume of pore phase|
|:||Cumulative porosity f that counts from large to small pore|
|:||Equivalent pore diameter|
|:||Lower, threshold, and upper limit of equivalent pore diameter|
|:||Slope of large and small pore range|
|:||Fractal dimension of large and small pore range|
|:||Partial porosity of large and small pore range|
|:||Step of iteration to generate large and small pore range|
|:||One-dimension number of subregions for large and small pore range|
|:||Proportion of iterating phase in generator for large and small pore range|
|:||Number of iterating phases in generator for large and small pore range|
|:||Intrinsic porosity of large and small pore range|
|:||Proportion of iterating phase being C-S-H gel|
|:||Effective ionic diffusivity of porous medium|
|:||Ionic diffusivity of bulk and inclusion phase|
|:||Volume fraction of inclusions|
|:||Ionic diffusivity in water corresponding to gel and capillary pore|
|:||Cementation exponent in Archie’s law|
|:||Non-steady-state migration coefficient|
|:||Absolute value of the applied voltage|
|:||Average value of the initial and final temperatures in anolyte solution|
|:||Thickness of the specimen|
|:||Average value of the penetration depth|
The authors declare that they have no competing interests.
The authors gratefully acknowledge the financial support from Key Basic Research Project (Grant No. 2015CB655105), State Key Laboratory for Track Technology of High-speed Railway (Grant no. 2014G004-F), and Jiangsu Natural Science Foundation (Grant no. BK20150621).
- L. Tang, “Electrically accelerated methods for determining chloride diffusivity in concrete—current development,” Magazine of Concrete Research, vol. 48, no. 176, pp. 173–179, 1996.
- E. J. Garboczi and D. P. Bentz, “Multiscale analytical/numerical theory of the diffusivity of concrete,” Advanced Cement Based Materials, vol. 8, no. 2, pp. 77–88, 1998.
- H. M. Jennings, J. W. Bullard, J. J. Thomas, J. E. Andrade, J. J. Chen, and G. W. Scherer, “Characterization and modeling of pores and surfaces in cement paste: correlations to processing and properties,” Journal of Advanced Concrete Technology, vol. 6, no. 1, pp. 5–29, 2008.
- K. Maekawa, T. Ishida, and T. Kishi, “Multiscale modeling of concrete performance integrated material and structural mechanics,” Journal of Advanced Concrete Technology, vol. 1, no. 2, pp. 91–126, 2003.
- G. W. Sun, W. Sun, Y. S. Zhang, and Z. Y. Liu, “Prediction of the effective diffusion coefficient of chloride ions in cement-based composite materials,” Journal of Materials in Civil Engineering, vol. 24, no. 9, pp. 1245–1253, 2012.
- D. Giménez, E. Perfect, W. J. Rawls, and Y. Pachepsky, “Fractal models for predicting soil hydraulic properties: a review,” Engineering Geology, vol. 48, no. 3-4, pp. 161–183, 1997.
- J. L. Pérez Bernal and M. A. Bello, “Fractal geometry and mercury porosimetry. Comparison and application of proposed models on building stones,” Applied Surface Science, vol. 185, no. 1-2, pp. 99–107, 2001.
- M. Arandigoyen, J. L. Perez Bernal, M. A. Bello López, and J. I. Alvarez, “Lime-pastes with different kneading water: pore structure and capillary porosity,” Applied Surface Science, vol. 252, no. 5, pp. 1449–1459, 2005.
- C. Atzeni, G. Pia, and U. Sanna, “Fractal modelling of medium-high porosity SiC ceramics,” Journal of the European Ceramic Society, vol. 28, no. 14, pp. 2809–2814, 2008.
- C. Atzeni, G. Pia, U. Sanna, and N. Spanu, “A fractal model of the porous microstructure of earth-based materials,” Construction and Building Materials, vol. 22, no. 8, pp. 1607–1613, 2008.
- B. B. Mandelbrot, P. Pfeifer, O. Biham, O. Malcai, D. A. Lidar, and D. Avnir, “Is nature fractal?” Science, vol. 279, no. 5352, pp. 783–786, 1998.
- D. N. Winslow, “The fractal nature of the surface of cement paste,” Cement and Concrete Research, vol. 15, no. 5, pp. 817–824, 1985.
- D. Pearson and A. J. Allen, “A study of ultrafine porosity in hydrated cements using small angle neutron scattering,” Journal of Materials Science, vol. 20, no. 1, pp. 303–315, 1985.
- M. Kriechbaum, G. Degovics, J. Tritthart, and P. Laggner, “Fractal structure of portland cement paste during age hardening analyzed by small-angle X-ray scattering,” Progress in Colloid and Polymer Science, vol. 79, pp. 101–105, 1989.
- D. A. Lange, H. M. Jennings, and S. P. Shah, “Image analysis techniques for characterization of pore structure of cement-based materials,” Cement and Concrete Research, vol. 24, no. 5, pp. 841–853, 1994.
- D. N. Winslow, J. M. Bukowski, and J. F. Young, “The fractal arrangement of hydrated cement paste,” Cement and Concrete Research, vol. 25, no. 1, pp. 147–156, 1995.
- X. Ji, S. Y. N. Chan, and N. Feng, “Fractal model for simulating the space-filling process of cement hydrates and fractal dimensions of pore structure of cement-based materials,” Cement and Concrete Research, vol. 27, no. 11, pp. 1691–1699, 1997.
- A. Heinemann, H. Hermann, and F. Häussler, “SANS analysis of fractal microstructures in hydrating cement paste,” Physica B: Condensed Matter, vol. 276–278, pp. 892–893, 2000.
- Y. Wang and S. Diamond, “A fractal study of the fracture surfaces of cement pastes and mortars using a stereoscopic SEM method,” Cement and Concrete Research, vol. 31, no. 10, pp. 1385–1392, 2001.
- Q. Zeng, K. F. Li, T. Fen-Chong, and P. Dangla, “Surface fractal analysis of pore structure of high-volume fly-ash cement pastes,” Applied Surface Science, vol. 257, no. 3, pp. 762–768, 2010.
- Y. Gao, K. Wu, and J. Jiang, “Examination and modeling of fractality for pore-solid structure in cement paste: starting from the mercury intrusion porosimetry test,” Construction and Building Materials, vol. 124, pp. 237–243, 2016.
- T. J. Miao, B. Yu, Y. G. Duan, and Q. T. Fang, “A fractal analysis of permeability for fractured rocks,” International Journal of Heat and Mass Transfer, vol. 81, pp. 75–80, 2015.
- S. Wang, T. Wu, H. Qi, Q. Zheng, and Q. Zheng, “A permeability model for power-law fluids in fractal porous media composed of arbitrary cross-section capillaries,” Physica A, vol. 437, pp. 12–20, 2015.
- G. Pia and U. Sanna, “Intermingled fractal units model and electrical equivalence fractal approach for prediction of thermal conductivity of porous materials,” Applied Thermal Engineering, vol. 61, no. 2-3, pp. 186–192, 2013.
- W. Wei, J. Cai, X. Hu, and Q. Han, “An electrical conductivity model for fractal porous media,” Geophysical Research Letters, vol. 42, no. 12, pp. 4833–4840, 2015.
- H.-Q. Jin, X.-L. Yao, L.-W. Fan, X. Xu, and Z.-T. Yu, “Experimental determination and fractal modeling of the effective thermal conductivity of autoclaved aerated concrete: effects of moisture content,” International Journal of Heat and Mass Transfer, vol. 92, pp. 589–602, 2016.
- S. Negrelli, R. P. Cardoso, and C. J. Hermes, “A finite-volume diffusion-limited aggregation model for predicting the effective thermal conductivity of frost,” International Journal of Heat and Mass Transfer, vol. 101, pp. 1263–1272, 2016.
- G. Pia, “High porous yttria-stabilized zirconia with aligned pore channels: morphology directionality influence on heat transfer,” Ceramics International, vol. 42, no. 10, pp. 11674–11681, 2016.
- R. A. Cook and K. C. Hover, “Mercury porosimetry of hardened cement pastes,” Cement and Concrete Research, vol. 29, no. 6, pp. 933–943, 1999.
- R. W. Zimmerman, “Elastic moduli of a solid containing spherical inclusions,” Mechanics of Materials, vol. 12, no. 1, pp. 17–24, 1991.
- M. Kachanov, I. Tsukrov, and B. Shafiro, “Effective moduli of solids with cavities of various shapes,” Applied Mechanics Reviews, vol. 47, no. 1, pp. 151–174, 1994.
- S. K. Kanaun and D. Jeulin, “Elastic properties of hybrid composites by the effective field approach,” Journal of the Mechanics and Physics of Solids, vol. 49, no. 10, pp. 2339–2367, 2001.
- I. Sevostianov and A. Giraud, “Generalization of Maxwell homogenization scheme for elastic material containing inhomogeneities of diverse shape,” International Journal of Engineering Science, vol. 64, pp. 23–36, 2013.
- L. Cui and J. H. Cahyadi, “Permeability and pore structure of OPC paste,” Cement and Concrete Research, vol. 31, no. 2, pp. 277–282, 2001.
- J. Zhang and Z. J. Li, “Application of GEM equation in microstructure characterization of cement-based materials,” Journal of Materials in Civil Engineering, vol. 21, no. 11, pp. 648–656, 2009.
- J. Zhao, J.-J. Zheng, and G.-F. Peng, “An effective medium approach for predicting the intrinsic permeability of heated cement paste,” Advances in Cement Research, vol. 27, no. 4, pp. 240–246, 2015.
- G. E. Archie, “Introduction to petrophysics of reservoir rocks,” American Association of Petroleum Geologists Bulletin, vol. 34, no. 5, pp. 943–961, 1950.
- H. S. Salem and G. V. Chilingarian, “The cementation factor of Archie's equation for shaly sandstone reservoirs,” Journal of Petroleum Science and Engineering, vol. 23, no. 2, pp. 83–93, 1999.
- L. Tang and L. Nilsson, “Rapid determination of the 317 chloride diffusivity in concrete by applying an electric field,” ACI Materials Journal, vol. 89, no. 1, pp. 49–53, 1992.
- Z. Yu and G. Ye, “New perspective of service life prediction of fly ash concrete,” Construction and Building Materials, vol. 48, pp. 764–771, 2013.
- E. J. Garboczi and D. P. Bentz, “Computer simulation of the diffusivity of cement-based materials,” Journal of Materials Science, vol. 27, no. 8, pp. 2083–2092, 1992.
- S. Kerisit and C. X. Liu, “Molecular simulations of water and ion diffusion in nanosized mineral fractures,” Environmental Science and Technology, vol. 43, no. 3, pp. 777–782, 2009.
Copyright © 2016 Yun Gao et al. 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.