Abstract
Reconstruction of random heterogeneous media has been an increasingly popular theme in materials science, as such media is primarily found in nature and manufactured materials. An ideal reconstruction includes every microstructural feature of the reference image and allows simulations of physical quantities that agree with experimental data. In this work, three 2D reconstructions and three 3D reconstructions are produced. All six are based on a single planar section. Its stereological and metallographic features are compared to those of the reference. 2D and 3D reconstructions from single planar sections were carried out with simulated annealing and a sampling method. The first reconstruction used the co-occurrence correlation function (CCF) and orthogonal sampling. The second employed the two-point correlation function () and the lattice point algorithm (LPA) as the sampling method. Finally, the third reconstruction used the together with the two-point cluster function () and the LPA as the sampling method. One of two sampling techniques (orthogonal and LPA) were used. The reconstructions that were done using and and LPA sampling provided the best results both in 2D and 3D cases, combining realistic morphology and good compatibility with the reference stereological measures.
Keywords:
Optimization; Simulated annealing; Heterogeneous media reconstruction; Microstructure; Characterization
1. Introduction
It is well-known that microstructure is a crucial issue in materials science. That is so because, to a considerable extent, macroscopic material properties, such as elasticity, yield strength, tensile strength, thermal and electrical conductivity, and fracture toughness, depend on the microstructure11 Oberwinkler B. Modeling the fatigue crack growth behavior of Ti-6Al-4V by considering grain size and stress ratio. Mater Sci Eng A. 2011;528(18):5983-92. http://doi.org/10.1016/j.msea.2011.04.046.
http://doi.org/10.1016/j.msea.2011.04.04...
2 Ghosh A, Sivaprasad S, Bhattacharjee A, Kar SK. Microstructure-fracture toughness correlation in an aircraft structural component alloy Ti-5Al-5V-5Mo-3Cr. Mater Sci Eng A. 2013;568:61-7. http://doi.org/10.1016/j.msea.2013.01.017.
http://doi.org/10.1016/j.msea.2013.01.01...
3 Xu J, Zeng W, Zhao Y, Jia Z. Effect of microstructure evolution of the lamellar alpha on impact toughness in a two-phase titanium alloy. Mater Sci Eng A. 2016;676:434-40. http://doi.org/10.1016/j.msea.2016.09.017.
http://doi.org/10.1016/j.msea.2016.09.01...
-44 Shao H, Zhao Y, Ge P, Zeng W. Crack initiation and mechanical properties of TC21 titanium alloy with equiaxed microstructure. Mater Sci Eng A. 2013;586:215-22. http://doi.org/10.1016/j.msea.2013.08.012.
http://doi.org/10.1016/j.msea.2013.08.01...
. Therefore, microstructural characterization and modeling are of particular interest.
Traditionally, one observes the microstructure on a planar section, whereas one typically seeks to determine three-dimensional microstructural quantities. A branch of mathematics, stereology, gives sound relationships between measures carried out on a planar section and certain three-dimensional features, such as volume fraction and interfacial area per unit of volume55 Russ JC, Dehoff RT. Practical stereology. 2nd ed. New York: Kluwer Academic/Plenum Publishers; 2000.. When full 3D reconstruction is required, a helpful technique is serial sectioning66 Núñez C, Domingo S. Statistical considerations on uniform grain size. Metall Trans, A, Phys Metall Mater Sci. 1988;19(12):2937-44. http://doi.org/10.1007/BF02647720.
http://doi.org/10.1007/BF02647720...
7 Krauss G, Marder AR. The morphology of martensite in iron alloys. Metall Trans. 1971;2(9):2343-57. http://doi.org/10.1007/BF02814873.
http://doi.org/10.1007/BF02814873...
8 Velichko A, Mücklich F. Quantitative 3D characterisation of graphite morphology in cast iron: correlation between processing, microstructure and properties. Int J Mater Res. 2009;100(8):1031-7. http://doi.org/10.3139/146.110148.
http://doi.org/10.3139/146.110148...
9 Kral MV, Mangan MA, Spanos G, Rosenberg RO. Three-dimensional analysis of microstructures. Mater Charact. 2000;45(1):17-23. http://doi.org/10.1016/S1044-5803(00)00046-2.
http://doi.org/10.1016/S1044-5803(00)000...
-1010 Baldissera MR, Rios PR, Hein LRO, Sandim HRZ. Three-dimensional characterization of pores in Ti-6Al-4V alloy. Mater Res. 2011;14(1):102-6. http://doi.org/10.1590/S1516-14392011005000014.
http://doi.org/10.1590/S1516-14392011005...
, which measures consecutive planar sections with an approximately constant spacing. Earlier works used "manual" calculation66 Núñez C, Domingo S. Statistical considerations on uniform grain size. Metall Trans, A, Phys Metall Mater Sci. 1988;19(12):2937-44. http://doi.org/10.1007/BF02647720.
http://doi.org/10.1007/BF02647720...
,77 Krauss G, Marder AR. The morphology of martensite in iron alloys. Metall Trans. 1971;2(9):2343-57. http://doi.org/10.1007/BF02814873.
http://doi.org/10.1007/BF02814873...
, but later works greatly benefitted from computational reconstruction methods using the planar sections as input. Nonetheless, it would be highly desirable if one could reconstruct the 3D microstructure from measurements on a single planar section. We investigate this possibility and its limitations in the present paper.
The microstructure can be identified as a "random set," "random heterogeneous material," or "random media"1111 Rios PR, Villa E. Transformation kinetics for inhomogeneous nucleation. Acta Mater. 2009;57(4):1199-208. http://doi.org/10.1016/j.actamat.2008.11.003.
http://doi.org/10.1016/j.actamat.2008.11...
,1212 Jiao Y, Stillinger FH, Torquato S. Modeling heterogeneous materials via two-point correlation functions: I. Basic principles. Phys Rev E Stat Nonlin Soft Matter Phys. 2007;76:031110. http://doi.org/10.1103/PhysRevE.76.031110.
http://doi.org/10.1103/PhysRevE.76.03111...
. By random media, one means a geometric interpretation in a scale far greater than the molecular scale but much smaller than the size of the specimen. The generation of realizations, either 2D or 3D, of "random media" is a problem of considerable current interest in materials science and engineering and various areas of knowledge.
Specifically, in metallurgy and materials science, macroscopic properties such as elastic modulus, elasticity, yield, tensile strength, fracture toughness, and electrical and thermal conductivity are microstructure-dependent. For example, specific properties such as the elastic modulus are less sensitive to the defects in a single crystal. Still, the apparent elastic modulus may depend on the microstructure when two microstructural components with distinct elastic moduli are present11 Oberwinkler B. Modeling the fatigue crack growth behavior of Ti-6Al-4V by considering grain size and stress ratio. Mater Sci Eng A. 2011;528(18):5983-92. http://doi.org/10.1016/j.msea.2011.04.046.
http://doi.org/10.1016/j.msea.2011.04.04...
2 Ghosh A, Sivaprasad S, Bhattacharjee A, Kar SK. Microstructure-fracture toughness correlation in an aircraft structural component alloy Ti-5Al-5V-5Mo-3Cr. Mater Sci Eng A. 2013;568:61-7. http://doi.org/10.1016/j.msea.2013.01.017.
http://doi.org/10.1016/j.msea.2013.01.01...
3 Xu J, Zeng W, Zhao Y, Jia Z. Effect of microstructure evolution of the lamellar alpha on impact toughness in a two-phase titanium alloy. Mater Sci Eng A. 2016;676:434-40. http://doi.org/10.1016/j.msea.2016.09.017.
http://doi.org/10.1016/j.msea.2016.09.01...
-44 Shao H, Zhao Y, Ge P, Zeng W. Crack initiation and mechanical properties of TC21 titanium alloy with equiaxed microstructure. Mater Sci Eng A. 2013;586:215-22. http://doi.org/10.1016/j.msea.2013.08.012.
http://doi.org/10.1016/j.msea.2013.08.01...
. Therefore, the reconstruction of a random media must be sufficiently accurate to be the starting point of computer simulations to obtain the properties listed above.
Early reconstruction papers1313 Quiblier JA. New three-dimensional modeling technique for studying porous media. J Colloid Interface Sci. 1984;98(1):84-102. http://doi.org/10.1016/0021-9797(84)90481-8.
http://doi.org/10.1016/0021-9797(84)9048...
,1414 Adler PM, Jacquin CG, Quiblier JA. Flow in simulated porous media. Int J Multiph Flow. 1990;16(4):691-712. http://doi.org/10.1016/0301-9322(90)90025-E.
http://doi.org/10.1016/0301-9322(90)9002...
used autocorrelation, probability density, and Gauss filters to obtain a binary matrix. Oolitic calcarium and Fontainebleau arenite, chosen for their simple microstructure, are the first reconstructed random media. These methods presume isotropy and can only use the most straightforward correlation functions. Their extension to multiphase systems and anisotropic systems is complicated.
Rintoul and Torquato1515 Rintoul MD, Torquato S. Reconstruction of the structure of dispersions. J Colloid Interface Sci. 1997;186(2):467-76. http://doi.org/10.1006/jcis.1996.4675.
http://doi.org/10.1006/jcis.1996.4675...
provided a better reconstruction methodology based on simulated annealing (SA). Those authors reported promising results for systems with low aggregation of pixels. However, Rintoul and Torquato concluded that even when there is good agreement between the reference correlation function and the correlation function obtained from the reconstructed microstructure, this does not necessarily imply that the reconstruction is faithful to the reference. Consequently, it is paramount to understand which correlation function is more satisfactory for each class of random media one wishes to reconstruct.
Torquato and Yeong1616 Yeong CLY, Torquato S. Reconstructing random media. Phys Rev E Stat Phys Plasmas Fluids Relat Interdiscip Topics. 1998;57(1):495-506. http://doi.org/10.1103/PhysRevE.57.495.
http://doi.org/10.1103/PhysRevE.57.495...
extended previous work to heterogeneous random media in a seminal paper. According to Torquato and Yeong, the main advantages of their new formulation are: (a) easy implementation, (b) the possibility of using diverse correlation functions, and (c) adaptability to multiphase and anisotropic systems. Their algorithm is briefly described in the Methodology section below. Yeong-Torquato's method's main characteristic is its flexibility. Torquato's work and Chen et al.1717 Chen D, Teng Q, He X, Xu Z, Li Z. Stable-phase method for hierarchical annealing in the reconstruction of porous media images. Phys Rev E Stat Nonlin Soft Matter Phys. 2014;89(1):013305. http://doi.org/10.1103/PhysRevE.89.013305.
http://doi.org/10.1103/PhysRevE.89.01330...
represent a class of reconstruction methods that employ optimization.
Other methods include multipoint statistics1818 Tahmasebi P, Sahimi M. Reconstruction of three-dimensional porous media using a single thin section. Phys Rev E - Stat Nonlinear. Soft Matter Phys. 2012;85:066709.
19 Hajizadeh A, Safekordi A, Farhadpour FA. A multiple-point statistics algorithm for 3D pore space reconstruction from 2D images. Adv Water Resour. 2011;34(10):1256-67. http://doi.org/10.1016/j.advwatres.2011.06.003.
http://doi.org/10.1016/j.advwatres.2011....
-2020 Okabe H, Blunt MJ. Prediction of permeability for porous media reconstructed using multiple-point statistics. Phys Rev E Stat Nonlin Soft Matter Phys. 2004;70:066135. that aim at decreasing computational costs. More recently, neural networks2121 Feng J, Teng Q, Li B, He X, Chen H, Li Y. An end-to-end three-dimensional reconstruction framework of porous media from a single two-dimensional image based on deep learning. Comput Methods Appl Mech Eng. 2020;368:113043. http://doi.org/10.1016/j.cma.2020.113043.
http://doi.org/10.1016/j.cma.2020.113043...
22 Mosser L, Dubrule O, Blunt MJ. Reconstruction of three-dimensional porous media using generative adversarial neural networks. Phys Rev E. 2017;96(4):043309. http://doi.org/10.1103/PhysRevE.96.043309.
http://doi.org/10.1103/PhysRevE.96.04330...
23 Bostanabad R, Bui AT, Xie W, Apley DW, Chen W. Stochastic microstructure characterization and reconstruction via supervised learning. Acta Mater. 2016;103:89-102. http://doi.org/10.1016/j.actamat.2015.09.044.
http://doi.org/10.1016/j.actamat.2015.09...
-2424 Zhang F, He X, Teng Q, Wu X, Dong X. 3D-PMRNN: reconstructing three-dimensional porous media from the two-dimensional image with recurrent neural network. J Petrol Sci Eng. 2022;208:109652. http://doi.org/10.1016/j.petrol.2021.109652.
http://doi.org/10.1016/j.petrol.2021.109...
have been employed. However, neural network models may take a long time to train. Nonetheless, they have the advantage that these models can quickly generate new microstructure realizations after training.
Several researchers have attempted to improve Yeong-Torquato's method2525 Tang T, Teng Q, He X, Luo D. A pixel selection rule based on the number of different-phase neighbours for the simulated annealing reconstruction of sandstone microstructure. J Microsc. 2009;234(3):262-8. http://doi.org/10.1111/j.1365-2818.2009.03173.x.
http://doi.org/10.1111/j.1365-2818.2009....
26 Ju Y, Huang Y, Zheng J, Qian X, Xie H, Zhao X. Multi-thread parallel algorithm for reconstructing 3D large-scale porous structures. Comput Geosci. 2017;101:10-20. http://doi.org/10.1016/j.cageo.2017.01.003.
http://doi.org/10.1016/j.cageo.2017.01.0...
-2727 Feng J, Teng Q, He X, Qing L, Li Y. Reconstruction of three-dimensional heterogeneous media from a single two-dimensional section via co-occurrence correlation function. Comput Mater Sci. 2018;144:181-92. http://doi.org/10.1016/j.commatsci.2017.11.030.
http://doi.org/10.1016/j.commatsci.2017....
. Feng et al.'s2727 Feng J, Teng Q, He X, Qing L, Li Y. Reconstruction of three-dimensional heterogeneous media from a single two-dimensional section via co-occurrence correlation function. Comput Mater Sci. 2018;144:181-92. http://doi.org/10.1016/j.commatsci.2017.11.030.
http://doi.org/10.1016/j.commatsci.2017....
co-occurrence correlation function (CCF) described in the Background section below is one such attempt.
A critical point in microstructure reconstruction is determining the correlation function, for example, the two-point correlation function (). A natural way of doing this is using the square or cubic lattice. For this reason, sampling in two orthogonal directions was initially employed. Torquato et al.2828 Jiao Y, Stillinger FH, Torquato S. A superior descriptor of random textures and its predictive capacity. Proc Natl Acad Sci USA. 2009;106(42):17634-9. http://doi.org/10.1073/pnas.0905919106.
http://doi.org/10.1073/pnas.0905919106...
suggested a new way of sampling random media. They proposed a "lattice point algorithm" (LPA). This formalism can be used to determine the two-point correlation function but can be extended for the two-point cluster correlation function ()2929 Torquato S, Beasley JD, Chiew YC. Two-point cluster function for continuum percolation. J Chem Phys. 1988;88(10):6540-7. http://doi.org/10.1063/1.454440.
http://doi.org/10.1063/1.454440...
.
In the present work, we employ the Yeong and Torquato1616 Yeong CLY, Torquato S. Reconstructing random media. Phys Rev E Stat Phys Plasmas Fluids Relat Interdiscip Topics. 1998;57(1):495-506. http://doi.org/10.1103/PhysRevE.57.495.
http://doi.org/10.1103/PhysRevE.57.495...
and the Feng et al.2727 Feng J, Teng Q, He X, Qing L, Li Y. Reconstruction of three-dimensional heterogeneous media from a single two-dimensional section via co-occurrence correlation function. Comput Mater Sci. 2018;144:181-92. http://doi.org/10.1016/j.commatsci.2017.11.030.
http://doi.org/10.1016/j.commatsci.2017....
methods to reconstruct 2D and 3D microstructures of nodular cast iron. Feng et al.2727 Feng J, Teng Q, He X, Qing L, Li Y. Reconstruction of three-dimensional heterogeneous media from a single two-dimensional section via co-occurrence correlation function. Comput Mater Sci. 2018;144:181-92. http://doi.org/10.1016/j.commatsci.2017.11.030.
http://doi.org/10.1016/j.commatsci.2017....
method was chosen because it requires low computational resources compared to other methods. The reconstructions were carried out from a single reference 2D planar section of the nodular cast iron. Nodular cast iron is widely used in many applications, and its microstructure is relatively simple.
2. Background
2.1. Two-point correlation function
Let be the indicator function of the phase and the position of a matrix pixel, and are the endpoints of :
The two-point correlation function is defined by the probability, , of finding, and in phases and . Specifically, the two-point correlation function for phase , , is the probability of finding and within the same phase :
where the angular brackets are symbols for mean value.
There are four possible two-point correlation functions, namely, , , and . Only one of these four is independent.
If the microstructure is statistically homogenous, the correlation function depends only on the relative positions of and . Moreover, if the microstructure is statistically isotropic, then the correlation function depends only on the distance between and :
where .
The relevant property of is that it is an even function. Therefore, for a square with pixels (assuming periodic boundary conditions), one needs to calculate the function up to pixels. Figure 1a illustrates which endpoints belong to the two-point correlation function.
Examples of correlation functions. (a) Two-point correlation function of one phase. (b) Two-point cluster correlation function of one phase. (c) Two-point lineal path correlation function of one phase.
2.2. Two-point cluster correlation function
Let be the indicator function of the phase and the position of a matrix pixel, and are the endpoints of :
The two-point cluster correlation function of phase can be defined as:
For homogenous and isotropic microstructures:
Figure 1b illustrates which endpoints belong to the two-point cluster correlation function.
2.3. Lineal path function
Let be the indicator function of the phase and the position of a matrix pixel, and are the endpoints of :
The two-point lineal path correlation function of phase can be defined as:
For homogenous and isotropic microstructures:
Figure 1c illustrates which endpoints belong to the two-point lineal path correlation function.
2.4. Co-occurrence correlation function
Feng et al.'s2727 Feng J, Teng Q, He X, Qing L, Li Y. Reconstruction of three-dimensional heterogeneous media from a single two-dimensional section via co-occurrence correlation function. Comput Mater Sci. 2018;144:181-92. http://doi.org/10.1016/j.commatsci.2017.11.030.
http://doi.org/10.1016/j.commatsci.2017....
co-occurrence correlation function can be defined as the probability, , of finding and in phases and :
For a homogeneous and isotropic microstructure and one has a
Thus utilizes the four two-point correlation functions even though only one is independent. Results from Feng et al.2727 Feng J, Teng Q, He X, Qing L, Li Y. Reconstruction of three-dimensional heterogeneous media from a single two-dimensional section via co-occurrence correlation function. Comput Mater Sci. 2018;144:181-92. http://doi.org/10.1016/j.commatsci.2017.11.030.
http://doi.org/10.1016/j.commatsci.2017....
suggest that, among other advantages, represents connectivity information better than using a single function.
3. Methodology
The reconstructions are based on experimental results from Freitas3030 Freitas AFP. Metalografia quantitativa tridimensional dos ferros fundidos nodular e cinzento [dissertation]. Volta Redonda: Universidade Federal Fluminense; 2008. on nodular cast iron. One of Freitas's micrographs from a planar section is our reference micrograph. The micrograph was binarized. A matrix of 300×300 square grid represented the micrograph for the 2D reconstruction. A matrix of 122×122 square grid represented the binarized image for 3D reconstruction.
For the reconstruction, one considers that the nodular cast iron is a two-phase isotropic microstructure. A function, , is determined for the reference matrix, for example, the two-point correlation function. A matrix with random black and white cells is the starting point. Black represents the second phase, and white represents the "parent-matrix" phase. The volume fraction of the second phase in the reference micrograph is equal to the volume fraction of the second phase, black cells, in the matrix. The same function is applied to the starting matrix that evolves from its starting point to the reference matrix. At a certain time-step, , this evolution of the function is calculated. A quantity that may be considered the energy is determined by:
For our matrix to evolve towards the reference matrix, thus minimizing , one interchanges two randomly selected cells of distinct phases. Of course, this interchange does not change the volume fraction. After the interchange is conducted, the energy is determined. The difference in energy between these two steps is .
The Metropolis algorithm is used to accept the interchange or not:
where is a "Monte Carlo temperature."
The objective of the reconstruction is that reaches a value such that the energy is as close as possible to zero. Specific details of this method applied to the present case can be found in Ferreira3131 Ferreira CCA. Reconstrução microestrutural 2D e 3D de ferro fundido nodular pelo método de Yeong-Torquato [dissertation] Volta Redonda: Universidade Federal Fluminense; 2021..
The simulated annealing was carried out here for the three functions defined in the Background section: the two-point correlation function (); a combination of the and the two-point cluster correlation function (); and the co-occurrence correlation function (CCF).
Two different sampling methods were used. The orthogonal sampling method and periodic boundary conditions were employed for the CCF2727 Feng J, Teng Q, He X, Qing L, Li Y. Reconstruction of three-dimensional heterogeneous media from a single two-dimensional section via co-occurrence correlation function. Comput Mater Sci. 2018;144:181-92. http://doi.org/10.1016/j.commatsci.2017.11.030.
http://doi.org/10.1016/j.commatsci.2017....
. For the and for the combination of the lattice point algorithm (LPA)2828 Jiao Y, Stillinger FH, Torquato S. A superior descriptor of random textures and its predictive capacity. Proc Natl Acad Sci USA. 2009;106(42):17634-9. http://doi.org/10.1073/pnas.0905919106.
http://doi.org/10.1073/pnas.0905919106...
was used. When LPA was employed, one did not use periodic boundary conditions. We cropped the reconstructed matrix border by 5% to minimize border effects. Full details of the methodology can be found in Ferreira3131 Ferreira CCA. Reconstrução microestrutural 2D e 3D de ferro fundido nodular pelo método de Yeong-Torquato [dissertation] Volta Redonda: Universidade Federal Fluminense; 2021..
4. Generation of 2D Microstructures from the Experimental Determination on a Single Plane
Figure 2 shows the microstructures generated by the 2D reconstructions from an experimental micrograph of nodular cast iron. In what follows, qualitative observations resulting from a visual inspection are presented.
2D reconstruction from a single experimental section, microstructures cut to final size 270×270 pixels. (a) Experimental-reference. (b) reconstruction using CCF. (c) reconstruction using . (d) hybrid reconstruction using .
Using the CCF method, the reconstructed microstructure reproduces well the smaller nodules of the reference microstructure, shown in Figure 2a. The spatial distribution of the nodules also looks good. By contrast, the nodule's morphology in Figure 2b differs from those in Figure 2a. Figure 2b depicts elongated instead of equiaxed nodules. This behavior suggests that CCF performs well regarding the spatial distribution of nodules. Nonetheless, CCF falls short as far as connectivity/morphology information is concerned.
Figure 2c presents the reconstruction done by the method using the LPA sampling. It is apparent in Figure 2c that LPA sampling results in a rounded morphology of the nodules much better than the CCF method. Nevertheless, one can see some elongated nodules analogous to those of CCF, Figure 2b. This result suggests that, like the CCF method, the method cannot preserve connectivity/morphology information well.
Figure 2d presents the microstructure reconstructed using a combination of and . Joint use of and results in a microstructure that closely approaches the reference microstructure, Figure 2a. The nodules have sizes and shapes similar to the experimental result. Few elongated nodules are present. All in all, the method conveys a result better than the CCF and methods. The area per unit of volume, , or equivalently, the mean intercept length are virtually the same, see Table 1.
2D reconstruction from a single experimental section: comparison of stereological measurements.
Table 1 presents a quantitative summary of the above results. The volume fraction, , is about the same for all three methods. It should be the same as it is kept constant in all three methods. The small difference is due to the 5% crop carried out in and methods.
The number of nodules per unit of area, , and Feret diameter, , using the and methods are close to those of the experimental microstructure. In contrast, and obtained in the microstructure generated by the CCF method differ from the experimental microstructure. This confirms the result in Figure 2, which suggests that the CCF method gives the worst reconstruction. It is fair to say that the CCF method has the advantage of being much less computer-intensive than the and methods besides it employs a worse sampling method, orthogonal sampling.
5. 3D Reconstruction of the Microstructure from a Single Experimental Section
Figure 3 shows the microstructures generated by 3D reconstructions from an experimental micrograph of nodular cast iron. Qualitative observations resulting from the visual inspection are presented in what follows. Figures 3b-3d are 2D sections of the 3D reconstructions. Those 2D sections are like those obtained by optical microscopy. They are compared with the reference microstructure in Figure 3a.
2D sections of 3D reconstructions from a single experimental section. (a) Experimental-reference. (b) reconstruction using CCF. (c) reconstruction using . (d) hybrid reconstruction using .
A comparison of the CCF reconstruction, Figure 3b, with the reference, Figure 3a, demonstrates that the CCF method cannot yield a nodule morphology resembling the experimental equiaxed nodule morphology. CCF nodules are not rounded, as is a typical graphite nodule. By contrast, the spatial distribution of the nodules looks reasonable. The reconstruction with the method results in a better morphology. An even better morphology can be obtained using the method. Likewise, Figure 2 and Figure 3 suggest that the method is the best of the three methods employed here. As mentioned above, the method utilizing considerably improves the connectivity representation. The inclusion of causes a better reconstruction of the nodule's shape. Only an occasional elongated nodule is present compared to the CCF and methods.
Table 2 displays a quantitative summary comparing the three reconstruction methods. All three reconstruction methods gave comparable values. In other words, the planar section of the 3D reconstructions resulted in similar stereological quantities. In contrast to the reconstruction of the 2D sections shown above, Table 2 shows that the stereological measurements on the 2D sections of 3D reconstructions did not reflect the discrepancy observed in Figure 3 among the reconstruction methods.
3D reconstruction from a single experimental section: comparison of stereological measurements.
Figure 4 displays a 3D opaque view of the reconstructions. Figures 4a-c exhibit a similar trend. Thus, the visually best reconstruction is obtained by method. The following best is the method using followed by the CCF method, which presented the worst result of the three methods.
3D opaque views of the reconstructions. (a) reconstruction using CCF. (b) reconstruction using S2. (c) hybrid reconstruction using .
The trend observed in Figure 4 becomes clearer using 3D views of the reconstructions. Figures 5a-c show transparent views of the reconstructions. The transparent views reinforce that the reconstruction is worse with CCF improves with and becomes much better with the method. The reconstruction is inferior to the CCF method. The shape of the nodules is far from the approximately spherical shape of the graphite nodules. The shape improves somewhat, but it still is not satisfactory with the method. Finally, the method results in a significant improvement.
3D transparent views of the reconstructions. (a) reconstruction using CCF. (b) reconstruction using . (c) hybrid reconstruction using .
6. Discussion
The CCF method is closely related to the two-point correlation functions of the microstructure2727 Feng J, Teng Q, He X, Qing L, Li Y. Reconstruction of three-dimensional heterogeneous media from a single two-dimensional section via co-occurrence correlation function. Comput Mater Sci. 2018;144:181-92. http://doi.org/10.1016/j.commatsci.2017.11.030.
http://doi.org/10.1016/j.commatsci.2017....
. According to Feng et al.2727 Feng J, Teng Q, He X, Qing L, Li Y. Reconstruction of three-dimensional heterogeneous media from a single two-dimensional section via co-occurrence correlation function. Comput Mater Sci. 2018;144:181-92. http://doi.org/10.1016/j.commatsci.2017.11.030.
http://doi.org/10.1016/j.commatsci.2017....
, the CCF method gives better reconstruction results than a method that uses the function only, as done here. Nonetheless, the results obtained here appear to indicate otherwise. The results from the CCF used here do not contradict Feng et al.'s2727 Feng J, Teng Q, He X, Qing L, Li Y. Reconstruction of three-dimensional heterogeneous media from a single two-dimensional section via co-occurrence correlation function. Comput Mater Sci. 2018;144:181-92. http://doi.org/10.1016/j.commatsci.2017.11.030.
http://doi.org/10.1016/j.commatsci.2017....
statement. This is because the sampling method used for CCF and methods were different. As described in the Methodology section, for CCF, one used orthogonal sampling, whereas for the LPA was used. This algorithm compensated a likely worse reconstruction when one uses the function only.
Still, a better sampling algorithm was not enough to cause a satisfactory result using solely the function. Further microstructural information is required. For this reason, we added the function to the function to get the method. The inclusion of the function improved connectivity information, and the overall results were better than those obtained by the CCF and the methods.
Figure 6 supports this reasoning for the 2D reconstructions. Figure 6 depicts the lineal path correlation function of the reconstructions by the three methods. This graph distinctly shows that the method gives the best results.
For the 3D reconstructions, Figure 7 corroborates method better representation of connectivity. In Figure 7, likewise in Figure 6, the connectivity, the function approaches that of the reference as the method changes from CCF, to the method.
Nevertheless, Figure 7 demonstrates that the method does possess a that is closer to the reference than the other methods but remains some distance from the reference .
Therefore, our results suggest that even if the method produced a reasonable reconstruction, the combination of and is unable to thoroughly represent all microstructural information from the nodular cast iron reference micrograph.
7. Conclusions
The Yeong-Torquato1616 Yeong CLY, Torquato S. Reconstructing random media. Phys Rev E Stat Phys Plasmas Fluids Relat Interdiscip Topics. 1998;57(1):495-506. http://doi.org/10.1103/PhysRevE.57.495.
http://doi.org/10.1103/PhysRevE.57.495...
method applied to a nodular cast iron's 2D and 3D microstructural reconstruction showed promising results. The quality of results varies according to the correlation functions and sampling method used in each case.
Stereological measurements typical of 2D reconstructions using and were compatible with the experimental-reference microstructure.
The 3D reconstructions using CCF, and presented interesting results. The best reconstruction was obtained by method, followed by the method, and the CCF method, which presented the worst result of the three methods. The same trend was obtained in the 2D reconstructions. However, our results suggest that although the methods produced a reasonable reconstruction, they could still not reproduce the complete microstructural information from the nodular cast iron reference micrograph. Despite that, even with some degree of uncertainty, they may be helpful to the researcher to get information and inferences from just one experimental planar section, with a relatively easy implementation and low computational cost.
Our results also demonstrate that the reconstruction technique proposed by Yeong and Torquato1616 Yeong CLY, Torquato S. Reconstructing random media. Phys Rev E Stat Phys Plasmas Fluids Relat Interdiscip Topics. 1998;57(1):495-506. http://doi.org/10.1103/PhysRevE.57.495.
http://doi.org/10.1103/PhysRevE.57.495...
shows a wide possibility of use in metallurgical and materials engineering. Using this method associated with other techniques may bring even better results in the future as new correlation functions and sampling techniques can be developed.
8. Acknowledgments
This study was financed in part by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior - Brasil (CAPES) - Finance Code 001. The authors are also grateful to Conselho Nacional de Desenvolvimento Científico e Tecnológico, CNPQ, and Fundação Carlos Chagas Filho de Amparo à Pesquisa do Estado do Rio de Janeiro, FAPERJ, for the financial support. The authors gratefully acknowledge the Center for Computational Materials Science, Institute for Materials Research, Tohoku University for the use of MASAMUNE-IMR (Materials science Supercomputing system for Advanced Multi-scale simulations towards Next-generation - Institute for Materials Research) supercomputer. In particular, the authors are grateful to Professors Yayoi Terada and Tetsuo Mohri for their collaboration through the Global Institute for Materials Research Tohoku (GIMRT) program.
9. References
-
1Oberwinkler B. Modeling the fatigue crack growth behavior of Ti-6Al-4V by considering grain size and stress ratio. Mater Sci Eng A. 2011;528(18):5983-92. http://doi.org/10.1016/j.msea.2011.04.046
» http://doi.org/10.1016/j.msea.2011.04.046 -
2Ghosh A, Sivaprasad S, Bhattacharjee A, Kar SK. Microstructure-fracture toughness correlation in an aircraft structural component alloy Ti-5Al-5V-5Mo-3Cr. Mater Sci Eng A. 2013;568:61-7. http://doi.org/10.1016/j.msea.2013.01.017
» http://doi.org/10.1016/j.msea.2013.01.017 -
3Xu J, Zeng W, Zhao Y, Jia Z. Effect of microstructure evolution of the lamellar alpha on impact toughness in a two-phase titanium alloy. Mater Sci Eng A. 2016;676:434-40. http://doi.org/10.1016/j.msea.2016.09.017
» http://doi.org/10.1016/j.msea.2016.09.017 -
4Shao H, Zhao Y, Ge P, Zeng W. Crack initiation and mechanical properties of TC21 titanium alloy with equiaxed microstructure. Mater Sci Eng A. 2013;586:215-22. http://doi.org/10.1016/j.msea.2013.08.012
» http://doi.org/10.1016/j.msea.2013.08.012 -
5Russ JC, Dehoff RT. Practical stereology. 2nd ed. New York: Kluwer Academic/Plenum Publishers; 2000.
-
6Núñez C, Domingo S. Statistical considerations on uniform grain size. Metall Trans, A, Phys Metall Mater Sci. 1988;19(12):2937-44. http://doi.org/10.1007/BF02647720
» http://doi.org/10.1007/BF02647720 -
7Krauss G, Marder AR. The morphology of martensite in iron alloys. Metall Trans. 1971;2(9):2343-57. http://doi.org/10.1007/BF02814873
» http://doi.org/10.1007/BF02814873 -
8Velichko A, Mücklich F. Quantitative 3D characterisation of graphite morphology in cast iron: correlation between processing, microstructure and properties. Int J Mater Res. 2009;100(8):1031-7. http://doi.org/10.3139/146.110148
» http://doi.org/10.3139/146.110148 -
9Kral MV, Mangan MA, Spanos G, Rosenberg RO. Three-dimensional analysis of microstructures. Mater Charact. 2000;45(1):17-23. http://doi.org/10.1016/S1044-5803(00)00046-2
» http://doi.org/10.1016/S1044-5803(00)00046-2 -
10Baldissera MR, Rios PR, Hein LRO, Sandim HRZ. Three-dimensional characterization of pores in Ti-6Al-4V alloy. Mater Res. 2011;14(1):102-6. http://doi.org/10.1590/S1516-14392011005000014
» http://doi.org/10.1590/S1516-14392011005000014 -
11Rios PR, Villa E. Transformation kinetics for inhomogeneous nucleation. Acta Mater. 2009;57(4):1199-208. http://doi.org/10.1016/j.actamat.2008.11.003
» http://doi.org/10.1016/j.actamat.2008.11.003 -
12Jiao Y, Stillinger FH, Torquato S. Modeling heterogeneous materials via two-point correlation functions: I. Basic principles. Phys Rev E Stat Nonlin Soft Matter Phys. 2007;76:031110. http://doi.org/10.1103/PhysRevE.76.031110
» http://doi.org/10.1103/PhysRevE.76.031110 -
13Quiblier JA. New three-dimensional modeling technique for studying porous media. J Colloid Interface Sci. 1984;98(1):84-102. http://doi.org/10.1016/0021-9797(84)90481-8
» http://doi.org/10.1016/0021-9797(84)90481-8 -
14Adler PM, Jacquin CG, Quiblier JA. Flow in simulated porous media. Int J Multiph Flow. 1990;16(4):691-712. http://doi.org/10.1016/0301-9322(90)90025-E
» http://doi.org/10.1016/0301-9322(90)90025-E -
15Rintoul MD, Torquato S. Reconstruction of the structure of dispersions. J Colloid Interface Sci. 1997;186(2):467-76. http://doi.org/10.1006/jcis.1996.4675
» http://doi.org/10.1006/jcis.1996.4675 -
16Yeong CLY, Torquato S. Reconstructing random media. Phys Rev E Stat Phys Plasmas Fluids Relat Interdiscip Topics. 1998;57(1):495-506. http://doi.org/10.1103/PhysRevE.57.495
» http://doi.org/10.1103/PhysRevE.57.495 -
17Chen D, Teng Q, He X, Xu Z, Li Z. Stable-phase method for hierarchical annealing in the reconstruction of porous media images. Phys Rev E Stat Nonlin Soft Matter Phys. 2014;89(1):013305. http://doi.org/10.1103/PhysRevE.89.013305
» http://doi.org/10.1103/PhysRevE.89.013305 -
18Tahmasebi P, Sahimi M. Reconstruction of three-dimensional porous media using a single thin section. Phys Rev E - Stat Nonlinear. Soft Matter Phys. 2012;85:066709.
-
19Hajizadeh A, Safekordi A, Farhadpour FA. A multiple-point statistics algorithm for 3D pore space reconstruction from 2D images. Adv Water Resour. 2011;34(10):1256-67. http://doi.org/10.1016/j.advwatres.2011.06.003
» http://doi.org/10.1016/j.advwatres.2011.06.003 -
20Okabe H, Blunt MJ. Prediction of permeability for porous media reconstructed using multiple-point statistics. Phys Rev E Stat Nonlin Soft Matter Phys. 2004;70:066135.
-
21Feng J, Teng Q, Li B, He X, Chen H, Li Y. An end-to-end three-dimensional reconstruction framework of porous media from a single two-dimensional image based on deep learning. Comput Methods Appl Mech Eng. 2020;368:113043. http://doi.org/10.1016/j.cma.2020.113043
» http://doi.org/10.1016/j.cma.2020.113043 -
22Mosser L, Dubrule O, Blunt MJ. Reconstruction of three-dimensional porous media using generative adversarial neural networks. Phys Rev E. 2017;96(4):043309. http://doi.org/10.1103/PhysRevE.96.043309
» http://doi.org/10.1103/PhysRevE.96.043309 -
23Bostanabad R, Bui AT, Xie W, Apley DW, Chen W. Stochastic microstructure characterization and reconstruction via supervised learning. Acta Mater. 2016;103:89-102. http://doi.org/10.1016/j.actamat.2015.09.044
» http://doi.org/10.1016/j.actamat.2015.09.044 -
24Zhang F, He X, Teng Q, Wu X, Dong X. 3D-PMRNN: reconstructing three-dimensional porous media from the two-dimensional image with recurrent neural network. J Petrol Sci Eng. 2022;208:109652. http://doi.org/10.1016/j.petrol.2021.109652
» http://doi.org/10.1016/j.petrol.2021.109652 -
25Tang T, Teng Q, He X, Luo D. A pixel selection rule based on the number of different-phase neighbours for the simulated annealing reconstruction of sandstone microstructure. J Microsc. 2009;234(3):262-8. http://doi.org/10.1111/j.1365-2818.2009.03173.x
» http://doi.org/10.1111/j.1365-2818.2009.03173.x -
26Ju Y, Huang Y, Zheng J, Qian X, Xie H, Zhao X. Multi-thread parallel algorithm for reconstructing 3D large-scale porous structures. Comput Geosci. 2017;101:10-20. http://doi.org/10.1016/j.cageo.2017.01.003
» http://doi.org/10.1016/j.cageo.2017.01.003 -
27Feng J, Teng Q, He X, Qing L, Li Y. Reconstruction of three-dimensional heterogeneous media from a single two-dimensional section via co-occurrence correlation function. Comput Mater Sci. 2018;144:181-92. http://doi.org/10.1016/j.commatsci.2017.11.030
» http://doi.org/10.1016/j.commatsci.2017.11.030 -
28Jiao Y, Stillinger FH, Torquato S. A superior descriptor of random textures and its predictive capacity. Proc Natl Acad Sci USA. 2009;106(42):17634-9. http://doi.org/10.1073/pnas.0905919106
» http://doi.org/10.1073/pnas.0905919106 -
29Torquato S, Beasley JD, Chiew YC. Two-point cluster function for continuum percolation. J Chem Phys. 1988;88(10):6540-7. http://doi.org/10.1063/1.454440
» http://doi.org/10.1063/1.454440 -
30Freitas AFP. Metalografia quantitativa tridimensional dos ferros fundidos nodular e cinzento [dissertation]. Volta Redonda: Universidade Federal Fluminense; 2008.
-
31Ferreira CCA. Reconstrução microestrutural 2D e 3D de ferro fundido nodular pelo método de Yeong-Torquato [dissertation] Volta Redonda: Universidade Federal Fluminense; 2021.
Publication Dates
-
Publication in this collection
27 May 2024 -
Date of issue
2024
History
-
Received
15 Sept 2023 -
Reviewed
05 Mar 2024 -
Accepted
28 Mar 2024