Research Article  Open Access
Spatial Distribution of Voids in Insulating Concrete Analyzed by MicroCT Images and Probability Functions
Abstract
Insulating concrete is a multiphase material designed for reduced thermal conductivity, and the void distribution in concrete strongly affects its physical properties such as mechanical response and heat conduction. Therefore, it is essential to develop a method for identifying the spatial distribution of voids. To examine the voids of insulating concrete specimens, microCT (computed tomography) images can be effectively used. The microCT images are binarized to visualize the void distribution and stacked to generate 3D specimen images. From the obtained images, the spatial distribution of the voids and the microscopic constituents inside the insulating concrete specimens can be identified. The void distribution in the material can be characterized using loworder probability functions such as twopoint correlation, linealpath, and twopoint cluster functions. It is confirmed that microCT images and loworder probability functions are effective in describing the relative degree of void clustering and void connectivity in insulating concrete.
1. Introduction
The phase distribution of multiphase materials strongly affects their properties [1–3]. The properties of concrete, a random heterogeneous multiphase material that is widely used in infrastructure, are affected by the spatial distribution of each constituent, namely, cement, sand, and aggregate. For example, the strength [4, 5] and percolation [6] of a concrete specimen are affected by the spatial distribution of voids within it; therefore, it is essential to investigate the spatial distribution of the constituents of the specimen to clarify its properties.
Insulating concrete is a type of concrete that is designed to have reduced thermal conductivity. Insulating concrete, when cast, contains numerous voids owing to the use of airentraining admixtures or lightweight aggregates. These voids are critical in reducing heat conduction. When the void ratio of the material is increased, certain properties of concrete behave contrarily; with a higher void ratio of the specimen, the insulation is increased, but the strength is decreased. To overcome this drawback, the introduction of anisotropic voids is a promising approach; several studies have examined this topic [5–7]. Therefore, proper methodology to examine the spatial distribution of voids in concrete specimen is necessary for investigating anisotropy of the void distribution in insulating concrete.
Various experimental techniques such as scanning electron microscopy (SEM), optical microscopy (OM), and mercury intrusion porosimetry (MIP) have been used for identifying the inner structure of the void distribution within insulating concrete specimens. SEM and OM can be used to identify the component distribution of cement paste; however, these methods can only provide information obtainable from crosssectional images. Additional preparation processes, such as polishing and drying, that may damage the specimen are required to obtain more information. MIP is considered effective for investigating the spatial distribution of voids [8, 9]; however, it cannot describe the overall void distribution within the material, and the specimen may be damaged during drying before MIP testing and pressure loading [10].
In this study, computed tomography (CT), a nondestructive and noninvasive method, is used to investigate the spatial distribution of the constituents of insulating concrete. In this method, a series of crosssectional images of insulating concrete are generated using Xrays without damaging the specimen. These images are expressed by a pixel, which serves as a unit for building twodimensional (2D) images. A CT imaging process that converts images from an 8bit image to a binary image is used to examine the specific constituents of insulating concrete. Subsequently, threedimensional (3D) images of the insulating concrete can be generated by stacking the 2D images.
Numerous studies have reported on the use of CT imaging for investigating material characteristics. Lu et al. [11] used μCT to examine the pore connectivity and its correlation with chloride migration in order to evaluate the permeability of concrete. Masad et al. [12] used CT to evaluate the permeability and anisotropy of asphalt by simulating the microstructure. Gallucci et al. [13] used μCT to investigate the degree of the pore network in concrete. CT images have also been used for studying other properties of concrete materials. Burlion et al. [14] used CT images to analyze the microstructure of cementitious materials during leaching. Chotard et al. [15] used CT to evaluate the early stage of the hydration of calcium aluminate cement. Helfen et al. [16] used CT to examine the evolution of the cement microstructure during hydration. Chung et al. [17] used CT images to evaluate the correlation between the anisotropy of the void distribution and the stiffness of lightweight aggregates. Rougelot et al. [18] used CT images to examine microcracking in cement composites. Zhang et al. [10] used μCT to investigate mass diffusivity in cement paste. These studies have adopted CT imaging methods to investigate the characteristics and properties of concrete materials.
In addition to CT imaging, probabilitybased methods are employed to describe the void distribution in insulating concrete. Here, the void distribution of insulating concrete was characterized using loworder probability functions such as twopoint correlation [3, 19], linealpath [1, 20], and twopoint cluster [21, 22]. The degree of phase clustering can be identified using a twopoint correlation function [1, 3]. The linealpath function can be used to obtain information about the phase connectivity along a specific direction. The twopoint cluster function is effective for describing the continuity of phaseclusteringrelated material characteristics such as percolation and diffusivity [1, 21].
In this study, the spatial distribution of voids in concrete specimens with different constituents is investigated. In insulating concrete, hollow glass beads (3M) with diameters of the order of a few micrometers are added to increase the insulating effect. Therefore, a μCT image expressed by a micrometer order pixel size is required to identify the voids induced by the glass beads in the insulating concrete specimen. To obtain highresolution CT images, the highenergy microCT facility in the 6D Xray imaging beamline at Pohang Accelerator Laboratory (PAL, Korea) is used. The spatial distribution of glass beads as well as voids in the insulating concrete is examined using μCT images and loworder probability functions. Using loworder probability functions, the characteristics of the void distribution can be described using only a small amount of data. Anisotropy of the void distribution inside the specimen is also examined using those methods. In addition, the difference between the void distribution of concrete specimens with and without glass beads is compared to investigate the effect of the glass beads on the cementitious material.
2. Insulating Concrete Specimen and μCT Imaging
2.1. Insulating Concrete Specimens
Insulating concrete is widely used as an exterior and interior structural material to improve energy efficiency by reducing heat loss through its insulating effect [23–25]. Studies have extensively focused on using lightweight aggregates to produce numerous voids in concrete and thus to improve the insulation [26, 27]. In this study, hollow glass beads (3M) were added to insulating concrete to enhance its insulating effect. The insulating concrete specimen is prepared using the mix design shown in Table 1. The material properties of the glass bead used are listed in Table 2. These beads have diameters of 10–65 m, and they are used to produce a spherical void within the concrete. The beads are expected to be more effective for insulation than lightweight aggregates because they can ensure regularly dispersed voids within the concrete.


An insulating concrete specimen composed Type I Portland cement, fine aggregates, and glass beads. The specimen is produced using a plastic tube with a length and diameter of 70 mm and 4 mm, respectively, as shown in Figure 1. The specimen volume is the maximum size that can be analyzed effectively using the Xray microCT device at PAL. The concrete specimen is cured under standard conditions (C) for 28 days. In this study, the fine aggregates used are washed ones having a fineness modulus (FM) of 2.64; aggregates smaller than 1.2 mm in size are only used in consideration of the specimen size.
2.2. CT Imaging for Constituents Analysis
It is difficult to investigate the constituent distribution of concrete without damaging the specimen. In this study, CT imaging is used to examine the spatial distribution of constituents as well as voids within insulating concrete. Crosssectional images of insulating concrete are obtained using the microCT facility at Pohang Accelerator Laboratory (PAL) in Korea. It is essential to use the synchrotron microCT device because the images from a general CT device (e.g., [6]) cannot visualize glass beads having diameters of the order of a few micrometers. At PAL, an Xray radiograph and a phase difference image with micrometer resolution can be obtained using the maximum 55 keV Xray source. Figure 2(a) shows the insulating concrete specimen setup and Figure 2(b) the synchrotron micro CT equipment. A crosssectional CT image of the insulating concrete specimen is generated when Xrays are injected from the device shown in Figure 2(c) and transmitted to the specimen.
Figure 3 shows sample 8bit microCT crosssectional images of the insulating concrete specimen, in which a pixel is represented by a value ranging from 0 to 255 (total of 256 values) depending on the relative density, where 0 is black and 255 is white in gray scale. A pixel value of 0 indicates that the relative density of the material for the pixel is minimum. Each microCT image is composed of 1633 × 1633 () pixels with 3.4 μm pixel size, and 1,500 images are obtained along the length direction over a 5.1 mm range. In the figure, the thin circular border is a plastic tube that encloses the insulating concrete specimen, and the blurred region outside the border is noise which is a byproduct of image reconstruction.
(a)
(b)
(c)
(d)
Assuming that cement is not completely hydrated and a unit quantity of water is fully evaporated, the constituents in the plastic tube are voids, glass beads ( = 0.125 kg/), mortar ( = 1.43 kg/), and the fine aggregates ( = 2.58 kg/) remaining after the hydration reaction, where is the material density. Therefore, the darkest region in Figure 3 can be considered as glass beads or voids and the brightest region as fine aggregate with relatively large density. The gray region observed across most of the specimen is considered to be mortar. As shown in Figure 3, the spatial distribution of the constituents in the insulating concrete specimen can be determined effectively and nondestructively using μCT images obtained from the synchrotron facility. To quantify the spatial distribution, probabilistic functions can be used.
3. LowOrder Probability Functions for Material Characterization
A suitable method is required for determining the spatial distribution of the constituents of the insulating concrete specimen. In this study, loworder probability functions, originally used to characterize the phase distribution of random heterogeneous materials [1], are used for this purpose.
3.1. TwoPoint Correlation Function
The twopoint correlation function () is used to characterize the phase distribution in this research. is the probability that any two points are located in the phase and phase [1, 3, 17, 19]. is the distance between two points; is the angle between a test line and the axis; and is the angle between the projection of a test line on the plane and the axis, as shown in Figure 4.
The twopoint correlation function only includes phase information at the end points of the test line; therefore, any additional information between two points is not included in this function. The general limits for the twopoint correlation function arewhere is the volume fraction of phase . In (1), when the distance between two points approaches zero, the probability of finding the same phase at both points converges to the phase volume fraction; otherwise, the probability of finding different phases at a point becomes zero. Equation (2) indicates that converges to the product of the phase volume fraction as increases. Gokhale et al. [19] suggested a simple form of the twopoint correlation function based on statistical information [28]In (3), is the number of intersections between a test line and the phase interface per unit test line length. In the case of twophase materials, the probability that both phases are located in the same () phase, , can be obtained from (3) and
3.2. LinealPath Function
To further investigate the constituent distribution, a linealpath function is used to compensate for the twopoint correlation function. The linealpath function, , is the probability that a randomly placed line of length is entirely in a phase [6, 29, 30]. Using the linealpath function, the continuous connectivity of the same phases for a specific direction can be evaluated. The general limits for the linealpath function areIn this study, Coker and Torquato’s method [20] is adopted to obtain the linealpath function for the insulating concrete specimen. To calculate the linealpath function, a random point inside the specimen is selected, and the phase information for the specific direction is analyzed by increasing the length of the test line. This process is repeated until a test line reaches a different phase from that of the starting point. From the linealpath function, the connectivity of the specific constituent can be evaluated because it considers the entire line segment between two end points of the test line.
3.3. TwoPoint Cluster Function
Another form of the probability function—the twopoint cluster function—is used to examine the connectivity of both the phase distribution and the phase clusters. The twopoint cluster function, , is defined as the probability of finding any two points in the same cluster of phase [1, 21]. Using the twopoint cluster function, the connected cluster of a phase can be identified even when the phase connectivity along the specific direction is interrupted by another phase. The general limits for the twopoint cluster function areTo calculate the twopoint cluster function, Lee and Torquato’s method [22] is used. In this method, two points in the material are arbitrarily selected, and it is checked whether these two points are located in the same phase and cluster. Then, the distance between these two points is calculated, and this process is repeated along the specific direction to obtain the twopoint cluster function with the change of . The twopoint cluster function can be used to obtain detailed information on the phase distribution because it provides information about the connectivity for both the phase and the cluster. Thus, the twopoint cluster function can be used for a material that needs to secure the channel pass, such as insulating and porous concrete.
4. Analysis of Constituents of Insulating Concrete Specimen
In this study, we determine the void distribution within insulating concrete using μCT images and analyze the spatial distribution of voids quantitatively using probability functions. The void distribution of insulating concrete strongly affects material properties such as thermal conductivity and strength. Using μCT images obtained from a synchrotron facility, voids within the specimen as well as visible ones on the surface are described. In addition, the spatial distribution of voids is statistically characterized using lowprobability functions.
4.1. CT Imaging of Insulating Concrete Specimen
The spatial distribution of the constituents of insulating concrete can be roughly identified from Figure 3. For example, regions 1–5 in Figure 5(a) are samples of circular voids. In this image, a circular void having a diameter of at most 20 pixels can be considered a glass bead because the maximum diameter of a glass bead is 65 μm (20 pixels based on a pixel size of 3.4 μm). Therefore, among the voids in Figure 5(a), regions 4 and 5 can be considered glass beads. Figure 5(b) shows a 3D image of a glass bead as obtained subsequently by stacking 8bit crosssectional images around region 5 in Figure 5(a); because a smooth sphere is identified, region 5 can be considered a candidate for a glass bead. As shown in Figure 5, the constituents of insulating concrete can be investigated using μCT images.
(a)
(b)
However, with the 8bit μCT image, it is difficult to detect minute components such as sparsely dispersed voids inside the specimen. In addition, the objective of the use of glass beads is to secure the void inside the specimen; the void phase, including the voids induced by glass beads, is the primary element of interest in this study. To describe the void distribution inside the insulating concrete effectively, a binary image is generated from the 8bit crosssectional image, as shown in Figure 6. To analyze the void distribution inside the insulating concrete specimen, the region of interest (Figure 6(b)) is extracted from the original μCT image shown in Figure 6(a). In the 8bit image in Figure 6(b), a set of pixels smaller than the threshold value is considered as voids. Here, a threshold value is selected using the Otsu method [31] from the graylevel histogram of the image, and the median filter is also used for noise reduction. By using an imaging process to separate voids and solids, a binary image is generated from the 8bit crosssectional image, as shown in Figure 6(c). For binary imaging, the image processing toolbox in MATLAB [32] is used. The binary image in Figure 6(c) is composed of 700 700 () pixels with a pixel size of 3.4 μm. In this figure, the white region is solid, and the black region is a void within the insulating concrete specimen. To describe the void distribution of the insulating concrete specimen in 3D, a series of binary images are subsequently stacked. The 3D image in Figure 7(a) is generated by stacking 700 binary images along the direction, and it is composed of 700 700 700 () voxels. The concrete model in Figure 7 is the region of interest (ROI) selected from the insulating concrete specimen. Figure 7(b) shows the void image of the specimen. From this image, both large void clusters and minute sparsely dispersed voids inside the specimen can effectively be detected; it can be denoted that the image resolution used in this study is sufficient to describe capillary pores including the void induced by glass beads. Using the image in Figure 7, the void ratio as well as the spatial distribution of voids inside the insulating concrete specimen can be identified. To quantify the void distribution in the specimen, the probabilistic methods discussed in Section 3 can be used.
(a)
(b)
4.2. Quantitative Description of Void Distribution Using Probability Functions
In this study, to quantitatively describe the spatial distribution of voids, three loworder probability functions—twopoint correlation (), linealpath (), and twopoint cluster () functions—are used. Here, the subscript indicates the void phase. Figure 8 shows probability functions for the ROI of the insulating concrete specimen for the , , and directions. In each function, is the distance between two points, and is the length of the specimen in the axial direction. The twopoint correlation function for voids, , is shown in Figure 8(a), and it is found to be similar for the , , and directions. The relative degree of void clustering for different directions can be identified using ; large values for a specific direction indicates that more voids are clustered in that direction. Therefore, Figure 8(a) indicates that the degree of void clustering in the specimen is almost isotropic or statistically identical for all axial directions.
(a)
(b)
(c)
Figure 8(b) shows for the , , and directions; it is similar in all three directions, indicating that the straightline void connectivity is isotropic in all three directions, even though a small difference is found in Figure 8(b). Similarly, the twopoint cluster function () in Figure 8(c) shows an isotropic void distribution when is 0–0.2. However, anisotropic fluctuation of is observed when is 0.2–0.6. This indicates the existence of anisotropic continuous void clusters with different sizes in the specimen. These results indicate that the relative size of the void clustering within a specimen can be quantitatively identified using the probability functions; Figure 8 indicates that a specific region of the specimen contains an anisotropic continuous void cluster. Continuous void clusters passing through the specimen that affect material properties, such as permeability, can also be identified using the twopoint cluster function. The existence of a void cluster passing through the entire specimen can be identified from the value when is 1; the function values are not equal to zero when the percolating void clusters exist in the specimen. In addition, the applicability of the specimen as a representative volume element (RVE, the smallest subset of material volume that describes representative properties) can be identified from the probability function; when the specimen is a RVE, each probability function for all directions should be the same because a RVE of concrete is isotropic. From the twopoint correlation and linealpath functions, the concrete sample investigated in this study can be considered as a RVE of the specimen because their values are almost the same for all axial directions. However, the result of the twopoint correlation function in Figure 8(c) demonstrates that this sample contains relatively large and anisotropic void clusters passing through the specimen, which confirms that this sample might not be considered as a RVE of the specimen; for more detailed investigation about the applicability of a RVE, the use of a larger sample will be needed.
To evaluate the effect of glass beads on the concrete specimen, CT images of the specimens with/without glass beads are compared. Figures 9(a) and 9(b), respectively, show crosssectional images of the specimen with and without glass beads. The two cases differ owing to the existence of glass beads in the specimen. The void ratio of the specimen with glass beads (Figure 9(a)) is 10.93%, and that of the specimen without glass beads (Figure 9(b)) is 6.84%. As shown in the binary images in Figure 9, the voids distributed between the fine aggregate (white region) are more distinctive in Figure 9(a) than in Figure 9(b). Voids induced by glass beads are found to restrict the chemical reaction between aggregates. Therefore, more voids are dispersed between aggregates in the specimen that contains glass beads, and this phenomenon can be inferred from the CT images in Figure 9.
(a)
(b)
The results of this study show that the void distribution within insulating concrete can be characterized using μCT imaging and probabilistic functions. In addition, the effect of specific constituents (e.g., glass beads) on the concrete material can also be examined from μCT images.
5. Conclusions
This study investigated methods for characterizing the spatial distribution of the constituents of insulating concrete. For insulating concrete, a specimen that contained hollow glass beads to increase the insulating effect was used. In general, glass beads having a diameter of a few micrometers can be utilized to secure regularly distributed voids within the specimen. To visualize the void distribution in the concrete specimen without damaging it, the highenergy microCT facility at the PAL synchrotron was used. The synchrotron device was necessary to obtain a highresolution image because the constituent, for example, glass bead, can only be described in the CT image with micrometer pixel size.
The specific constituents in the insulating concrete specimen were identified using μCT images. In particular, the void distribution within the ROI of the insulating concrete was described using a series of μCT images and CT image processing. The spatial distribution of voids in the specimen could be visualized using the 3D digital insulating concrete image generated by the subsequent stacking of 2D binarized images. To quantify the void distribution of the specimen, probabilistic methods such as twopoint correlation, linealpath, and twopoint cluster functions were adopted. From these loworder probability functions, the degree of void clustering and void connectivity can be characterized. The isotropy of the void distribution in the specimen was also demonstrated from the probability functions.
The results confirm that μCT imaging with a highresolution image and a loworder probability function can be effectively used to examine the spatial distribution of the voids in insulating concrete. Anisotropy of the void distribution was effectively examined by using these methods effectively, and the applicability of the specimen as a RVE was also identified. In particular, the heterogeneity and anisotropy of the void distribution were identified by using these methods, which can be used to design an optimized insulating concrete model with appropriate void size and distribution. Further simulation studies of heat conduction and mechanical behavior for insulating concrete samples might be necessary to evaluate the relationship between void distribution and material properties, such as thermal conductivity and strength of insulating concrete, in detail.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This research was supported by a Korea Research Foundation grant funded by the Korean Government (KRF20110029212 and KRF2012R1A1A2006629). This work was also supported by the Industrial Strategic Technology Development Program (10041589) funded by the Ministry of Knowledge Economy (MKE, Korea). The authors acknowledge Pohang Accelerator Laboratory (PAL) for beamline use (Grant 20121st6D004). In addition, the authors express their appreciation to Dr. Kwang Soo Youm at GS E&C Research Institute in the Republic of Korea, who provided insulating concrete specimens and physical property data.
References
 S. Torquato, Random Heterogeneous Materials, Springer, New York, NY, USA, 2002. View at: Publisher Site  MathSciNet
 T.S. Han and P. R. Dawson, “Representation of anisotropic phase morphology,” Modelling and Simulation in Materials Science and Engineering, vol. 13, no. 2, pp. 203–223, 2005. View at: Publisher Site  Google Scholar
 S.Y. Chung and T.S. Han, “Reconstruction of random twophase polycrystalline solids using loworder probability functions and evaluation of mechanical behavior,” Computational Materials Science, vol. 49, no. 4, pp. 705–719, 2010. View at: Publisher Site  Google Scholar
 D. N. Winslow, M. D. Cohen, D. P. Bentz, K. A. Snyder, and E. J. Garboczi, “Percolation and pore structure in mortars and concrete,” Cement and Concrete Research, vol. 24, no. 1, pp. 25–37, 1994. View at: Publisher Site  Google Scholar
 R. A. Dorey, J. A. Yeomans, and P. A. Smith, “Effect of pore clustering on the mechanical properties of ceramics,” Journal of the European Ceramic Society, vol. 22, no. 4, pp. 403–409, 2002. View at: Publisher Site  Google Scholar
 S.Y. Chung and T.S. Han, “Correlation between loworder probability distribution functions and percolation of porous concrete,” Magazine of Concrete Research, vol. 65, no. 7, pp. 448–460, 2013. View at: Publisher Site  Google Scholar
 R. Cabrillac, B. Fiorio, A.L. Beaucour, H. Dumontet, and S. Ortola, “Experimental study of the mechanical anisotropy of aerated concretes and of the adjustment parameters of the introduced porosity,” Construction and Building Materials, vol. 20, no. 5, pp. 286–295, 2006. View at: Publisher Site  Google Scholar
 S. Diamond, “Mercury porosimetry: an inappropriate method for the measurement of pore size distributions in cementbased materials,” Cement and Concrete Research, vol. 30, no. 10, pp. 1517–1525, 2000. View at: Publisher Site  Google Scholar
 R. Kumar and B. Bhattacharjee, “Study on some factors affecting the results in the use of MIP method in concrete research,” Cement and Concrete Research, vol. 33, no. 3, pp. 417–424, 2003. View at: Publisher Site  Google Scholar
 M. Zhang, Y. He, G. Ye, D. A. Lange, and K. V. Breugel, “Computational investigation on mass diffusivity in Portland cement paste based on Xray computed microtomography (μCT) image,” Construction and Building Materials, vol. 27, no. 1, pp. 472–481, 2012. View at: Publisher Site  Google Scholar
 S. Lu, E. N. Landis, and D. T. Keane, “Xray microtomographic studies of pore structure and permeability in Portland cement concrete,” Materials and Structures, vol. 39, no. 6, pp. 611–620, 2006. View at: Publisher Site  Google Scholar
 E. Masad, A. Al Omari, and H.C. Chen, “Computations of permeability tensor coefficients and anisotropy of asphalt concrete based on microstructure simulation of fluid flow,” Computational Materials Science, vol. 40, no. 4, pp. 449–459, 2007. View at: Publisher Site  Google Scholar
 E. Gallucci, K. Scrivener, A. Groso, M. Stampanoni, and G. Margaritondo, “3D experimental investigation of the microstructure of cement pastes using synchrotron Xray microtomography (μCT),” Cement and Concrete Research, vol. 37, no. 3, pp. 360–368, 2007. View at: Publisher Site  Google Scholar
 N. Burlion, D. Bernard, and D. Chen, “Xray microtomography: application to microstructure analysis of a cementitious material during leaching process,” Cement and Concrete Research, vol. 36, no. 2, pp. 346–357, 2006. View at: Publisher Site  Google Scholar
 T. J. Chotard, M. P. BoncoeurMartel, A. Smith, J. P. Dupuy, and C. Gault, “Application of Xray computed tomography to characterise the early hydration of calcium aluminate cement,” Cement and Concrete Composites, vol. 25, no. 1, pp. 145–152, 2003. View at: Publisher Site  Google Scholar
 L. Helfen, F. Dehn, P. Mikulik, and T. Baumbach, “Threedimensional imaging of cement microstructure evolution during hydration,” Advances in Cement Research, vol. 17, no. 3, pp. 103–111, 2007. View at: Publisher Site  Google Scholar
 S.Y. Chung, T.S. Han, T. S. Yun, and K. S. Youm, “Evaluation of the anisotropy of the void distribution and the stiffness of lightweight aggregates using CT imaging,” Construction and Building Materials, vol. 48, pp. 998–1008, 2013. View at: Publisher Site  Google Scholar
 T. Rougelot, N. Burlion, D. Bernard, and F. Skoczylas, “About microcracking due to leaching in cementitious composites: Xray microtomography description and numerical approach,” Cement and Concrete Research, vol. 40, no. 2, pp. 271–283, 2010. View at: Publisher Site  Google Scholar
 A. M. Gokhale, A. Tewari, and H. Garmestani, “Constraints on microstructural twopoint correlation functions,” Scripta Materialia, vol. 53, no. 8, pp. 989–993, 2005. View at: Publisher Site  Google Scholar
 D. A. Coker and S. Torquato, “Extraction of morphological quantities from a digitized medium,” Journal of Applied Physics, vol. 77, no. 12, pp. 6087–6099, 1995. View at: Publisher Site  Google Scholar
 S. Torquato, J. D. Beasley, and Y. C. Chiew, “Twopoint cluster function for continuum percolation,” The Journal of Chemical Physics, vol. 88, no. 10, pp. 6540–6547, 1988. View at: Publisher Site  Google Scholar  MathSciNet
 S. B. Lee and S. Torquato, “Measure of clustering in continuum percolation: computersimulation of the twopoint cluster function,” The Journal of Chemical Physics, vol. 91, no. 2, pp. 1173–1178, 1989. View at: Publisher Site  Google Scholar
 K.Y. Shin, S.B. Kim, J.H. Kim, M. Chung, and P.S. Jung, “Thermophysical properties and transient heat transfer of concrete at elevated temperatures,” Nuclear Engineering and Design, vol. 212, no. 1–3, pp. 233–241, 2002. View at: Publisher Site  Google Scholar
 W. Feist, J. Schnieders, V. Dorer, and A. Haas, “Reinventing air heating: convenient and comfortable within the frame of the Passive House concept,” Energy and Buildings, vol. 37, no. 11, pp. 1186–1203, 2005. View at: Publisher Site  Google Scholar
 S. Grieu, O. Faugeroux, A. Traoré, B. Claudet, and J.L. Bodnar, “Artificial intelligence tools and inverse methods for estimating the thermal diffusivity of building materials,” Energy and Buildings, vol. 43, no. 23, pp. 543–554, 2011. View at: Publisher Site  Google Scholar
 S.C. Ng and K.S. Low, “Thermal conductivity of newspaper sandwiched aerated lightweight concrete panel,” Energy and Buildings, vol. 42, no. 12, pp. 2452–2456, 2010. View at: Publisher Site  Google Scholar
 A. Sales, F. R. de Souza, W. N. Dos Santos, A. M. Zimer, and F. D. C. R. Almeida, “Lightweight composite concrete produced with water treatment sludge and sawdust: thermal properties and potential application,” Construction and Building Materials, vol. 24, no. 12, pp. 2446–2453, 2010. View at: Publisher Site  Google Scholar
 E. Underwood, Quantitative Stereology, AddisonWesley, Reading, Mass, USA, 1970.
 B. Lu and S. Torquato, “Linealpath function for random heterogeneous materials,” Physical Review A, vol. 45, no. 2, pp. 922–929, 1992. View at: Publisher Site  Google Scholar
 H. Singh, A. M. Gokhale, S. I. Lieberman, and S. Tamirisakandala, “Image based computations of lineal path probability distributions for microstructure representation,” Materials Science and Engineering A, vol. 474, no. 12, pp. 104–111, 2008. View at: Publisher Site  Google Scholar
 N. Otsu, “A threshold selection method from graylevel histograms,” IEEE Transactions on Systems, Man and Cybernetics, vol. 9, no. 1, pp. 62–66, 1979. View at: Publisher Site  Google Scholar
 MATLAB, Version 7.10 (R2010a), The MathWorks, Natick, Mass, USA, 2010.
Copyright
Copyright © 2015 SangYeop Chung 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.