Image Anal Stereol 2007;27:79-85 Original Research Paper LOCAL STUDY OF DEFECTS DURING SINTERING OF UO2: IMAGE PROCESSING AND QUANTITATIVE ANALYSIS TOOLS Eric Girard1,2, Jean-Marc Chaix1, François Valdivieso2, Patrice Goeuriot2 and Jacques Lechelle3 1SIMAP –INPGrenoble-CNRS-UJF B.P.75, F-38402 Saint Martin d'Hères Cedex, France; 2Ecole Nationale Supérieure des Mines de Saint-Etienne, Centre SMS, CNRS UMR 5146 PECM, 158 Cours Fauriel, F-42023 Saint-Etienne cedex 02, France; 3Centre d'Etudes de Cadarache, CEA, F-13108 Saint Paul lez Durance cedex, France e-mail: jean-marc.chaix@ltpcm.inpg.fr; francois.valdivieso@emse.fr; patrice.goeuriot@emse.fr; jacques.lechelle@cea.fr (Accepted. May 6, 2008) ABSTRACT This paper describes the image analysis tools developed and used to quantify the local effects, of heterogeneities during sintering of ceramic materials applied in nuclear fuels. Specific materials, containing a controlled dispersion of well defined heterogeneities (dense or porous aggregates) in the ceramics matrix have been prepared and sintered. In order to characterize the materials in the vicinity of these likely isolated heterogeneities, large SEM images are first acquired around heterogeneities. Image processing tools are used to segment the relevant elements of the microstructure: ceramic grains and porosity (intra- and inter-granular). An analysis procedure is then developed to analyse the evolution of the main microstructure features (pore volume fraction, pore size, grain size) during sintering, with the distance to the matrix-heterogeneity interface as reference parameter. An evaluation of the anisotropy of repartition of pores and grains around the heterogeneities is also proposed. Some illustrative results are presented, which evidence a rather short range of mean interactions, and a significant rather long range of the anisotropy effect related to the matrix-heterogeneities interactions during the undidirectional compaction. Keywords: ceramics, image analysis, materials science, metallography, microstructure, sintering, sintering heterogeneities. INTRODUCTION Heterogeneities are always present in green compacts and evolve during sintering (Lange, 1984; Ciftcioglu et al., 1987; Lim et al. 2000). These heterogeneities (or “defects”), such as aggregates or large pores, are the source of microstructural defects which have a serious and detrimental influence on mechanical and thermo-mechanical properties of the material. The effects of these heterogeneities have been observed and in some cases interpreted (Evans, 1982; Uematsu et al., 1995; Shinohara et al., 2000), but most of the quantifications are macroscopic: although quantitative image analysis has been used by many authors in the field of sintering (Chaix, 2000) experimental quantitative analysis of local microstructures evolution during sintering with hetero- geneities are still lacking. The present paper is part of a study aiming at a better understanding of the effects of different defects on the sintering of their surrounding microstructure by means of a quantitative analysis of the microstructure characteristics (Girard 2004; Girard et al., 2004; 2005). To characterize these matrix-defect interactions, it was decided to perform a local analysis, i.e., an analysis of local variations of the main microstructure characteristics (pore volume fraction, pore size, grain size) in the vicinity of the defect-matrix boundary, and to apply this to a series of sintering experiments in different conditions: nature of the heterogeneities, sintering time and temperature. The present paper is centered on the image processing and image analysis tools developped, adapted and used for this materials science purpose. 79 Girard E et al: Image analysis tools for defects in UO2 sintering MATERIALS AND METHODS EXPERIMENTS The ceramic matrix was an UO2 powder, and the defects were made from the same powder, to avoid uncontrolled chemical effects of the heterogeneities. 20 p m (a) (b) (c) (d) Fig. 1. Local microstructure for different defects and sintering conditions: (a) porous aggregate, 1500 °C, 1minute. (b) porous aggregate, 1700 °C, 1 h. (c) dense aggregate, 1500°C, 1minute. (d) dense aggregate, 1700 °C, 1 h. Two types of defects have been used (Fig. 1): porous aggregates obtained from the natural aggregates of the powder, and dense aggregates obtained from powders compressed at 200MPa and broken by ball milling. Both kinds of aggregates were calibrated in the range 40-60 µm by sieving, then sintered in different conditions to get different levels of consolidation. Green compacts containing 3 vol% of each specific kind of aggregates were prepared and sintered at different temperatures: 1500°C, 1600°C, 1700 °C, 1700 °C 20 min or 1700 °C 60 min. This small volume fraction of defects dispersed in the matrix was chosen in order to have mostly independant defects, far enough from one another while keeping a high enough probability to observe them on cross sections. The cross sections of the samples were systematically realized in plans containing the comp action axis, and the images orientation carefully chosen so that the “y axis” corresponds to the compaction direction. Thermally etched polished sections of the samples were investigated in the vicinity of likely isolated heterogeneities. Only 40-50 µm defect sections were selected, to ensure that the analysed defects were cut in the vicinity of their diametral plane, so that the bias on the evaluation of the distance to the boundary was limited. A fine mapping of each heterogeneity and its surrounding matrix was needed. The examples on Fig. 1 illustrate the constraints on the image size and resolution: elementary objects (pores and grains) range from typically 0.1 to 10 µm, and the effects of interactions must a priori be characterized at distances up to several defect diameters. A compromise was adopted to get images with sufficient resolution and large enough fields. This first needed the careful acquisition of several elementary images (typically 7 × 9 images with 10% overlapping) for each defect, to cover the whole field of interest. These high resolution SEM images (Jeol JSM 6500 FEG-SEM, 15 kV acceleration voltage) were first processed to reduce the noise (Overton and Weymouth, 1979), limit the initial resolution, and homogenize the grey levels between the different elementary images. They were in a second step assembled in a single image using the “MIAC” (Multiple Image Alignment) module of the Analysis software (Soft Imaging System). The assembled images were typically 3500 x 3500 pixels, with effective pixel size of 45 nm, covering about 150 x 150 µm2. SEGMENTATION OF IMAGES: PORES AND GRAIN BOUNDARIES. An essential point for image analysis measurements is to have reliable segmented images in which the objects of interest, i.e., grains and pores, are clearly identified. In the present case, there is a strong dispersion in the characteristics of images (quantities and sizes): grain boundaries are differently revealed by the thermal etching (varying grey level and apparent width); grain sizes strongly differ from one sintering process to another; in some cases small and large grains are present in the same areas; contrast and luminosity differ from one image to another... 80 Image Anal Stereol 2007;27:79-85 (a) (b) (c) (d) Fig. 2. Segmentation of pores. (a) grey level image. (b) entropy threshold. (c) maximum contrast threshold, (d) geodesic dilation of (c) inside (b). This classic practical problem, i.e., that each microstructure has almost to be treated with its own specificities (Wojnar, 1999) has lead here to develop up to 5 adapted procedures, corresponding to the 5 sintering cycles. However, the methodology is based on similar processes, corresponding to similar basic image characteristics, with only some variations in the nature, size, and order of use of the filtering operations. The segmentation of pores (Fig. 2) is based on the combination of two automatic thresholding methods. Entropy thresholding (Kapur et al., 1985) leads to correct pores boundaries, but also detects dark objects which are not porosity (parts of grain boundaries and of dark grains). Maximum contrast thresholding (Kohler, 1981) does not lead to these undesired objects, but underestimates the pores. A correct reconstruction of pores is then obtained by performing geodesic dilation, based on classic mathematical morphology (Serra, 1982; Coster and Chermant, 1989), of the second image inside the first one. In some cases some enlarged grain boundaries are at this step considered as pores. This parasitic porosity is due to the excess of thermal etching which leads to a local enlargement of some grain boundaries; it is eliminated if necessary by morphologic tools. (a) (b) (c) (d) Fig. 3. Segmentation of grain boundaries and final extraction of information from the grey level initial image of Fig. 2a. (a) rough determination of grain boundaries. (b) grain boundaries + pores. (c) individual identification of grains. (d) separation of intra (green) and inter(red)-granular pores. The segmentation of grains (Fig. 3) involves 3 main steps: i. prefiltering of grey level image: classic noise reduction filters have been performed when necessary. To facilitate the final checking of segmentation, contrast enhancement has also been used: JASC PaintShop Pro 6 and Adobe Photoshop softwares proved to be suitable for this “eye-oriented” purpose. ii. extraction of a rough image of grain boundaries by over-thresholding: the obtained image (Fig. 3a) shows enlarged boundaries and undesired objects, but the grain boundaries are detected as an almost continuous network. iii. identification of grain boundaries using mathematical morphology transformations. The basic process involves small objects elimination (openings), skeleton transform which reduce the network to one pixel but leaves undesired objects and some disconnected grain boundaries, chain closure to reconnect them, and clipping of pending chains to eliminate separated objects. A chain closure method was specifically developed 81 Girard E et al: Image analysis tools for defects in UO2 sintering for some types of images (Fig. 4). The obtained one pixel wide network is then added to the previous pore image (Fig. 2d) to get the image of separated grains (Fig. 3b). ADCIS Aphelion Software was used here. Fig. 4. Chain closure procedure: for a given chain endpoint, a search area is defined from direction of the last two points (black line) and an angular parameter (here ± p / 4), and a maximum distance. The closest point obeing this constraints if any is used to close the chain. Finally, grains are individually identified (Fig. 3c). Combining the grains image (Fig. 3b) and the pores one (Fig. 2d), the pore population can be separated into intra- and inter- granular pores, located inside grains and at grain boundaries respectively (Fig. 3d). QUALITY CONTROL The efficiency of these procedures is of course not full, so that some excess or lacking grain boundaries always exist (Fig. 5). It is of great importance to evaluate the consequences of these defaults on the measurements. In the present work, some series of images have been manually corrected, and the measurements on uncorrected and corrected images have been compared. Some parameters are almost insensitive to these corrections (e.g., the total pore volume fraction or the pore sizes, as far as most corrections are concerned with grain boundaries), some are more sensitive (the intra-granular porosity for instance, for which less than 1% error in absolute value represents up to 20%, in relative value). In most cases, the relative difference remains acceptable for important parameters, such as in the case of grain size (Fig. 6). (1) > , 0 ** . ». '¦'%. • • 'i ,. v'" ..- % *.l • • j 5 |jm (2) (3) Fig. 5. Examples of manual corrections on images in a porous aggregate (top) and in a dense aggregate (bottom). (1) grey level image – (2) automatic segmentation – (3) hand corrected image. 4 3.5 3 2.5 2 1.5 1 1600° r 1500° 1700°C 20' _ ->5*** Q. '*jir------ 1700°C 0); erosions are used for inside distances (r < 0). The results have first been analysed (Girard et al, 2004) on the basis of the evolution of profiles to get informations on the local evolution during sintering. For instance, in the case of dense aggregates (1600 °C) (Fig. 9), the influence of distance r from the defect edge is limited to a relatively short range, smaller than the defect radius, while a strong influence (“constrained sintering”) is observed inside the defect. -15 -10 -5 0 5 10 15 20 25 30 35 40 distance r to interface (µm) Fig. 9. Evolution of the mean grain diameter as a function of distance r to the defect-matrix interface, for different sintering temperatures, in the case of dense aggregates consolidated at 1600 °C. Dotted lines indicate the mean value in the matrix. 83 Girard E et al: Image analysis tools for defects in UO2 sintering ANISOTROPY Angular measurements were also realized in order to analyse the effects of compaction. For that purpose, directional measurements in y (compaction axis) and x (perpendicular axis) directions, by measuring angular sectors of aperture a in both directions (Fig. 10). In practice, a = p/3 was used, corresponding to 3 equivalent angular directions : x, y, and intermediate. Fig. 10. Principle of anisotropy measurements. The parameters are measured at each distance r from the defect boundary in a given angular sector a in the considered direction. An example of results in the case of porous aggregates is presented on Fig. 11. The compaction induced anisotropy in the pore repartition is clearly evidenced, with a rather long typical range of action of about 30 µm, i.e., the order of the defect radius. Similar results are obtained for dense aggregates. 25 20 M 5 D- - ¦ 1600°C x ¦òr ¦ ¦ 1600°C y — 1600 average 1700°C x 1700°C y 1700 average M bd "'^<*%fi^ 10 20 30 40 distance from interface (pm) 50 Fig. 11. Anisotropic repartition of porosity (VVP) for porous aggregates at two temperatures (a = p/3). The lines indicate large distance average values in all directions. An index of anisotropy for any parameter p was defined: \{px-py dr \pdr (1 where px, py are the values in the x and y directions respectively, and p the value averaged on all directions. In practice, the integration is limited to a maximum value of 40 µm, above which the differences are only irrelevant noise errors. The evolution of the anisotropy index for porosity in the case of porous aggregates (Fig. 12) shows that anisotropy first increases during the sintering process, then tends to be reduced. 0.6 0.5 0.4 0.3 0.2 0.1 0.0 agi ag2 ag3 1500°C 1600°C 1700°C 1700°C 1700°C 20' 1h Sintering Cycle Fig. 12. Evolution of the anisotropy of pore volume fraction along the sintering cycle for porous aggregates. Anisotropy is almost constant for ag1 samples, which correspond to aggregates with small cohesion (presintered at 1600 °C); it is initially smaller but increases during the sintering in the case of defects with increasing strength ag2 and ag3, presintered at 1700 °C for 15 and 120 min respectively. CONCLUSION The methods developed in this work allow to get useful local quantitative microstructural data to study the effects of heterogeneities in the sintering of ceramics. The work is based on classic image processing and measuring tools, which have been adapted and in some cases improved to fit the specific needs of the studied materials: local and extended analysis of the microstructure evolution around structural localized defects. Some examples of results have been presented. We can also mention that the large number of local measurements realized in this work have lead to a set of microstructural data which I p 0 0 84 Image Anal Stereol 2007;27:79-85 form a database of microstructure evolution. A first attempt of statistical analysis of this database has been realized (Girard et al., 2005) REFERENCES Chaix JM (2000). Image analysis, microstructure evolution and sintering mechanisms. In: Wojnar L, Rozniatowski K, eds. STERMAT 2000. Proceedings of the 6th International Conference on Stereology and Image Analysis in Materials Science. 2000 Sep 20-23; Krakow, Poland, 69-78. Ciftcioglu M, Akinc M, Burkhart L (1987). Effect of agglomerate strength on sintered density for Yttria powders containing agglomerates of monosize spheres. J Am Ceram Soc 70:C329-34. Coster M, Chermant JL (1989). Précis d'analyse d'images. Paris: Editions du CNRS. De Hoff RT, Rhines FN (1968). Quantitative microscopy. New York: McGraw Hill. Evans AG (1982). Considerations of inhomogeneity effects in sintering. J Am Ceram Soc 65:497-501. Girard E (2004). Frittage du dioxyde d’uranium en présence d’hétérogeneites - caractérisation de l’évolution microstructurale. PhD Thesis, ENSMSE, Saint Etienne, France. Girard E, Chaix JM, Carry C, Valdivieso F, Goeuriot P, Léchelle J (2005). Microstructure evolution of ceramics during sintering: an analysis based on local image analysis measurements in the vicinity of controlled defects. In: 4th Int. Conf. on Science, Technology and Applications of Sintering, 2005 sep 1-5; Grenoble, France, 436-9. Girard E, Valdivieso F, Goeuriot P, Chaix JM, Léchelle J, Chiarelli G (2004). Sintering with Heterogeneities Characterization and Modelling of the Microstructural Evolution. Key Engineering Materials 264-268:197-200. Kapur JN, Sahoo PK, Wong AKC (1985). A new method for grey level picture thresholding using the entropy of the histogram. Comput Vision Graphics and Image Proc 29:273-85. Kohler R (1981). A segmentation system based on thresholding. Comput Vision Graphics and Image Proc 15:319-38. Lange FF (1984). Sinterability of agglomerated powders. J Am Ceram Soc 67(2):83-9. Lim LC, Wong P.M, Jan M (2000). Microstructural evolution during sintering of near-monosized agglomerate-free submicron Alumina powder compacts. Acta Materialia 48:2263-75. Overton KJ, Weymouth TE (1979). A noise reducing preprocessing algorithm. In: Proc. IEEE Computer Society Conf. on Pattern Recognition and Image Processing; Chicago, Illinois, 498-507. Serra J (1982). Image analysis and mathematical morphology. New-York: Academic Press. Shinohara N, Okumiya M, Hotta T, Nakahira K, Naito M, Uematsu K (2000). Morphological changes in process-related large pores of granular compacted and sintered Alumina. J Am Ceram Soc 83:1633-40. Uematsu K, Zhang Y, Uchida N (1995). Characteristics and formation mechanisms of defects in Alumina powder compacts. Ceram Trans 51:263-70. Wojnar L (1999). Image analysis: Applications in material engineering. New-York: CRC Press. 85