Image Anal Stereol 2004;23:45-52 Original Research Paper ANALYSIS OF POROUS STRUCTURE IN PLASMA-SPRAYED COATING R adka Jüzkovä1, P avel C tibor2 and V iktor B eneš3 1 Private College of Economic Studies, Ltd., Lindnerova 575/1, 18000 Prague 8-Liben, Czech Republic, Institute of Plasma Physics ASCR, Za Slovankou 3, 182 21 Prague 8-Liben, Czech Republic, 3Faculty of Mathematics and Physics, Charles University, Sokolovska´ 83, 18600 Praha 8-Karl´in, Czech Republic e-mail: radka.juzkova@svses.cz, ctibor@ipp.cas.cz, benesv@karlin.mff.cuni.cz (Accepted February 23, 2004) ABSTRACT The paper is devoted to the quantitative analysis of the microstructure of a ceramic plasma-sprayed coating. The aim is to describe the size and spatial distribution of approximately convex pores. A whole 3D specimen was subsequently grinded to obtain serial sections of small distance. Using image analysis techniques the centroids and size of individual pores are detected and the obtained data analysed by means of summary spatial statistics. The tests of complete spatial randomness are performed using various 3D characteristics and important practical conclusions for the material fabrication procedure are obtained. Keywords: ceramic coating, morphological, porosity, second order characteristics, spatial distribution, volume reconstruction. INTRODUCTION The analysis of the microstructure of ceramic plasma-sprayed coatings is of great importance (Ctibor et al, 2002; 2003). A typical specimen, sprayed in the Institute of Plasma Physics ASCR, Prague, Czech Republic, is investigated in this paper. Its microstructure forms pores of various shapes. In the analysis approximately convex pores are detected and measured, their size and spatial distribution is a subject of interest. For hard materials, usually a single section is available that is studied by means of stereological methods (Ohser and Mu¨cklich, 2000). To avoid ill-posed stereological problems, with great laboratory effort, we obtained serial sections along the whole specimen. This approach enhances the quantitative analysis considerably, since three-dimensional (3D) distance and second order characteristics can be directly estimated. However preliminary image analysis of serial sections requires careful and tedious work. Also, a computer program for 3D pore reconstruction has been prepared. The choice of estimators of contact distribution, nearest neighbour distance distribution, K- and L-function is discussed in the literature (Ohser and Mu¨cklich, 2000; Baddeley and Gill, 1997; Moller and Waagepetersen, 2003) and we contribute with our practical results. Standard statistical Monte Carlo tests are applied for complete spatial randomness and the alternative of mild clustering is recognised. Porosity affects important functional properties of any thermal spray coating like corrosion resistance, cohesion under induced mechanical stress/strain, electric resistivity and permittivity and others. The factors governing the porosity level and also some spatial porosity features (such as clustering, inhomogeneity) are given by plasma spray gun parameters combination. The detection of pores clustering means, for us, that more attention has to be paid to powder feeding (especially the diameter of the feeder nozzle vs. plasma stream diameter, feeding point position adjustment). MATERIAL Titanates, in general, form a wide and important group of electronic ceramics. A mixture of MgTiO3 and CaTiO3, with a Mg:Ca ratio equal to 94:6, was plasma sprayed at the Institute of Plasma Physics ASCR, Prague, Czech Republic (Ctibor et al, 2002). It was proven that the microstructure of a selected material is representative enough to serve as a pilot material for studying the microstructure of ceramic plasma deposits in general. A sample was manufactured using a water-stabilised plasma spray system WSP® PAL 160. The feedstock powder size was 63-125 um, final coat thickness was about 1 mm. Volume studied in present paper has dimensions 450 x 350 x 240 um and is situated far enough from a free surface and from the substrate to prevent material “edge effects” (Fig 1). Lamellas (usually called splats) are results of spreading and solidification of droplets accelerated against the substrate. The lamellas are principal structural units, basic “bricks” of a microstructural architecture of the coatings, in general more 45 JÜzkovÄ R et al: Analysis of porous structure in plasma-sprayed coating important for coating properties than crystallites local orientation, grain boundaries or similar structural features. The steel substrate covered by plasma-sprayed ceramic coating was cut by a low speed saw (Leco, USA) equipped with a diamond wheel to prevent structural modifications. A cross section of the substrate-coating system was mounted in a resin. Then, grinding was completed with a rough grain diamond paste. Subsequent polishing of individual planes for large area section imaging was performed manually using diamond polishing paste (Struers, Denmark). Fig. 1. The investigation of 3D specimen of plasma-sprayed ceramic coating. Eighty serial sections are obtained perpendicular to the grinding direction. Each section is a montage of 5x5 images observed by an optical microscope. The spraying direction is perpendicular to the substrate surface. The creation of the large-area section involved a sequence of steps, namely: 1. Identification of a suitable large region on the first cross section plane covering 25 fields of view of the optical microscope, Neophot 32 (Zeiss, Germany), and RGB camera, TK-1070E (JVC, Japan). 2. Creation of a seamless “montage” of images covering all fields in this region of the first plane. 3. Physical removal of a small amount of material by polishing and identification of the same region on the second plane. 4. Creation of a seamless “montage” of all microstructural fields of view of the same region on the second plane. 5. Measurement of the distance between the subsequent planes by measuring the decrease in the sizes of pyramidal markers. The resulting “montage” (cf. Tewari and Gokhale, 2000), was essentially a microstructural image of a large area (25 fields of view at 800 x magnification in this case) having a high resolution (~0.13 um/pixel). By repeating of the steps 1-5 the stack of large area images is prepared. A combination of the resolution in one plane and the distance between planes 3 um was used with respect to the main structural features of plasma-sprayed coatings. Together, 80 aligned planes, each consisting of 25 high-resolution “montage” images were documented and all 2000 images stored in the computer to undergo the image analysis. One large section after montage represents 450 x 350 x 240 um, i.e. 3460 x 2690 pixels x 240 um (the z axis represents 80 times removed 3 um by polishing). METHODS As described in the previous paragraph the whole three-dimensional (3D) specimen was cut in serial sections using grinding. The image of a section is two-dimensional (2D). The analysis of microstructural images was performed using automatic image analysis softwares Microimage (Olympus, 1999), and Lucia G3.52. The image analysis performed on each single image consists of the following steps. Segmentation was done using the same threshold level for all pictures. It leads to the isolation of pores of all types (including both cracks and approximately convex pores). Using morphological opening (Serra, 1982) with spherical structural element other cracks than pores are deleted (Fig. 2). Another way based on erosion/reconstruction as described in Soille (1999), was not used since we believe that the selection induced by opening does not influence our statistics substantially. The reason is that pores do not contain much sharp corners of boundaries. Then, the parameters of pore sections are measured (centroid coordinates, area, perimeter, angle and length of maximum diameter). In the second step, a computer program was prepared for the reconstruction of 3D pores from the section measurements. The centroid coordinates, number of intersected sections, minimum and maximum diameter and volume were estimated for each detected pore. The aim of stereological methods used is to reconstruct 3D parameters of the microstructure from the measured parameters in 2D sections. Classical stereological methods (Ohser and Mu¨cklich, 2000), are not needed because the stack of serial 46 Image Anal Stereol 2004;23:45-52 sections in fact yields 3D information directly. The spatial distribution is described by 3D contact distribution, nearest neighbour distribution, K—, L— and the pair correlation function. Estimators of these characteristics are taken from Baddeley and Gill (1997), Moller and Waagepetersen (2003) and Ohser and Mu¨cklich (2000). Special care should be given to the correction of edge effects. Comparing these estimated functions with those corresponding to a Poisson process the deviations from spatial randomness are classified, which is important for the description of the material studied. Consider a convex 3D set W, a subset of the Euclidean space R3. In our application, the centroids of pores in a specimen W, which is a rectangular parallelepiped (see Fig. 1), are observed. Let us denote the centroids X1,...,Xn, they are considered as a realisation of a homogeneous isotropic point process as defined in Moller and Waagepetersen (2003). By (B), the (random) number of points a set B is denoted. The intensity X of such a process is the mean number ofpoints in a unit volume, i.e. X =E([0, 13). An unbiased estimator of intensity is X = n/V(W), where V(W) is the volume of W and n is a number of observed points. Another estimator for X2 (Ohser and Mu¨cklich, 2000), is The L-function is just a normed version of K—function: A2 n(n— 1) VW)2 ¦ The distance between points Xi,Xj is denoted d ij = 11 X — X j 11 , where 11.11 is the Euclidean norm. The window function cW is defined on R3 as cW(h) = V(W 0 {o,h}), i.e. by morphological erosion of W with a two point set {o,h}, o being the origin of coordinate system. Its isotropised version on R1 is c— W{r) 4 1 S V(WG{o,(r,u)})du, where S2 is the unit sphere in R3 and (r,u) polar coordinates (r size, u orientation) of a vector h€ R3. The K—function of is defined as XK(r) = E<š>(b(o,r) | o G «I»), where on the right hand side there is the mean number of points in a ball b(o,r) = {x € R3; \\x\\ < r} under the condition that o is the point of the process . Since this condition has zero probability the Palm theory (Soille, 1999) is used for a rigorous definition. The relationship between K—function and the pair correlation function g is for r>0 K{r) = 4n ( x2g{x)dx. (1) L(r)=[—K(r)) ,r>0 (2) For the Poisson point process, it holds L(r) = r, r > 0. An unbiased estimator of K—function is used from Ohser and Mu¨cklich (2000): For the rectangular parallelepiped W with side lengths a < b < c, we have the window function . . , 2(ab + bc + ac) c— W{r) = abc----------------------r + 2 „3 2{a + b + c)r r 3tt 4n , 0 to its neighbour. Its distribution function D(r) = P(D < r) can be written as (Xi is any point of ) D{r) = 1-P[^{b{Xi,r)) = 1), r>0. (4) Spherical contact distance H can be defined as a distance of the nearest point of from the origin. Its distribution function HS is HS(r) = 1-P(0. (5) The point process is homogeneous Poisson if and only if it holds 4% D(r) = H S(r) = 1-exp{-X3 r 3},r>0. (6) The deviations from this theoretical form are used to recognise the interactions within the point process, either attractive or repulsive, possibly on different scales. The clustering of the point process is observed in the graph for the distribution function of nearest neighbour distances when the curve is above the curve (Eq. 6) and in the graph for contact distribution function when the curve is below the curve (Eq. 6). 47 JÜzkovÄ R et al: Analysis of porous structure in plasma-sprayed coating In the opposite case, for distribution function of nearest neighbour distances below and for the contact distribution function above the curve (Eq. 6), we observe regularity in the point process. Estimation of D(r) and Hs(r) is a problem because of edge effects, since for the points near to the boundary of W the nearest neighbour need not be observed. Several methods for edge effect corrections are described in Baddeley and Gill (1997), Ohser and Mu¨cklich (2000), Stoyan et al. (1995). We used the so called Kaplan-Meier estimators, see Baddeley and Gill (1997). Let d j be the nearest neighbour distance of Xj g W and p(Xj, dW) be the distance of Xj to the boundary of the window W for each j = 1,... n. Then, the Kaplan-Meier estimator of D(r) is D(r) = 1-n 1- N{j:dj = s,djs,p{XhdW)>s} )' (7) where r > 0, setting 0/0 = 0. N denotes the number of points (cardinality) of the set. To estimate the contact distribution function, a random or systematic point test set is used. We used a binomial point pattern {X{, ...,X'} of p random test points in W and denote d'j the contact distances of nearest point of the process in W to the test point, p(XpdW) denotes the distance of Xj to the boundary of the window W for each j = 1,...n. Then, the Kaplan-Meier estimator Hs used is r U{ N{j:d)>s,p(Xtj,dW)>s} J (8) where r > 0, setting 0/0 = 0. N denotes the number of points (cardinality) of the set. One of the main aims of the statistical analysis of measured data was to test the hypothesis of complete spatial randomness of convex pores within the specimen. Monte Carlo tests (Ohser and Mu¨cklich, 2000) present a standard tool for this purpose. Nineteen independent realisations of the Poisson point process with the same intensity as that estimated from data are simulated. For each realisation the K—, L—, Hs—, D—functions are estimated (using Eqs. 3, 2, 8, 7, respectively) and their pointwise maximum and minimum is plotted. Using more realisations the bounds might be smoother, but as will be seen in the numerical results, the amount 19 is enough for our analysis. If the estimator of the same function obtained from data gets outside the plotted bounds for some r, we reject the null hypothesis of complete spatial randomness on a 95% confidence level at distance r. NUMERICAL RESULTS The microstructure of a ceramic plasma-sprayed deposit specimen of size 450 x 350 x 240 um was monitored on 80 parallel sections considered as equidistant with distance 3 um. A typical result of morphological operations is in Fig. 2. ¦ Ji W ¦ ~ rj * : '.'j>r^<;;ry ^ * 100 pm Fig. 2. The binary image of the 2D microstructure in a section including pores and cracks (left) and the same microstructure after morphological opening (right). The important pore volume histogram (Fig.3) is the starting point for further analysis, since from a physical point of view it is desirable to study pores of defined size classes. The pores with maximum diameter between 3-40 um were detected. The mode of the distribution is about 4-6 um. 100 150 200 250 300 350 400 450 500 550 BOO Volume of pores (pm3) Fig. 3. Pore volume histogram; the volume of pores was estimated using Cavalieri method, i.e. V = tLAi, where t is the distance between sections and Ai is the area corresponding to convex pores in 2D images to 3D pores. The vertical lines in histogram show limits of size categories (lines corresponding to diameters equal to 4,6,10\\,m)for a spherical shape model. 48 Image Anal Stereol 2004;23:45-52 The size categories were selected with respect to lamellae thickness, which is for this spraying technique and material typically given 3 to 6 um. Interlamellar globular pores are usually smaller than 4 um. Also, the smallest size class used (4-6 um) is supposed to correspond dominantly to pores between lamellas. Larger size class (6-10 um) is “intermediate”, while the largest one (10-20 um) is supposed to correspond dominantly to voids in proximity of not well melted (i.e. not well flattened) lamellas. Also the different size classes could exhibit different behaviour. But, given the present state of the art, the above expressed opinion is not more than a hypothesis. In Table 1, important size classes are defined with corresponding pore numbers in the sample corresponding to studied size classes. Table 1. Number of pores in corresponding size categories Diameter of pores Number of pores any 2116 4-20 um 1601 4-6 um 598 6-10 um 548 10-20 um 455 The experiment of plasma spraying, as described in the Materials section, may lead to hypoteses that there are some anisotropies or gradients in the spatial distribution of pores. A simple method to detect them is to plot projections of all pores onto coordinate planes. In Figs. 4-6 these plots, moreover, include information on pore size and no apparent phenomena of that kind are observed. A finer analysis is based on summary statistics. In Figs. 7-9 the results are presented for different size classes, as suggested in Table 1. In each Figure the estimated L-, HS-, and D-functions are analysed. The argument r of the functions corresponds to the distance. The starting value of r in the graph is equal to the lower bound of diameters of pores in the studied size class because the distance between centroids of pores cannot be lower than the sum of their half diameters (e.g. the starting value of r for the class 4-6 um is 4). The value of L-function at the starting value of r is equal to this value for the Poisson process and, therefore, the vertical axis of the graph starts at this r. Fig. 4. Projection on the planes XY, XZ and YZ (see Fig. 1)for the pores of diameter equal to 4-6 \im. The rejection of the hypotesis of complete spatial randomnes is mostly apparent for L-function (in all three size classes) and for D-function (for smaller pores). Since these functions lie above the bounds for the Poisson process in some intervals, rejection tends towards clustering. For HS- function, this corresponds to values below the bounds, which is apparent for the smallest size class and 27 < r < 40 um. Finally, if the whole set of events between 4-20 um is evaluated, then since the smaller pores prevail, the null hypothesis is again rejected towards the alternative of clustering. JÜzkovÄ R et al: Analysis of porous structure in plasma-sprayed coating Projection on XY warnte % Tii t < S3&Äfc^?1 240 200j M' 130 J BO -C' 0 Projection ¦ nT h*" i"** SO 100 150 200 250 300 35C Fig. 5. Projection on the planes XY, XZ and YZ (see Fig. 1)for the pores of diameter equal to 6-10 \im. Fig. 6. Projection on the planes XY, XZ and YZ (see Fig. 1)for the pores of diameter equal to 10-20 \im. CONCLUSIONS Globular pores represent one category of voids in plasma-sprayed coatings. These pores are distributed isotropically and their potential to deteriorate coating properties is proportional to their size. Moreover, all size categories of globular pores are prospective for the connection between the knowledge of the structure and the knowledge of behaviour of powder in the plasma. Spherical pores arise as defects in powder partly solidified before impacting on the substrate and as defects during spreading of impacting particles. Plasma spraying is a very rapid process — in the setup used for the sample manufacturing studied approximately 6x105 powder particles are in 1 cm3 of plasma jet and travel together with them towards the substrate surface by velocity of hundreds of m/s. It means that each 0.001 s new 1.8x105 particles are in the same volume. Particles are mixed with plasma which burns in the air. From the thermal point of view plasma gives to the particles energy for complete melting before access of the substrate. From the kinetic point of view plasma gives energy for spreading of particles on the substrate. 50 Image Anal Stereol 2004;23:45-52 Fig. 7. The diagrams indicate values of L—function L(-), the nearest neighbour distribution function D(-) and the spherical contact distribution function HS(-) for pores of diameter equal to 4-6 /j,m. The starting value of r corresponds to minimum of the theoretical distance between the centroids of the smallest pores (in this case, 4). It must be noted that the relative velocity of particles versus the surrounding environment (plasma + air) is not constant. The whole process is very dynamic and fluctuations of plasma parameters are at a time scale high enough to have an influence on plasma volume occupied in the same moment by tens of particles or more. The microstructure of the final coating contains structural features inherited from the “history” of individual particles in the plasma jet. Fig. 8. The diagrams indicate values of L—function L(-), the nearest neighbour distribution function D(-) and the spherical contact distribution function HS(-) for pores of diameter equal to 6-10 /j,m. The starting value of r corresponds to minimum of the theoretical distance between the centroids of the smallest pores (in this case, 6). The present paper is based on data gained from microstructure on a high level of characterisation ability. Precise alignment of the in-flight behaviour to our microstructural data is behind the possibility of the present state of the art. But, as was documented earlier (Ctibor et al., 2003), there exists clustering of special structural features in the coating. Mild clustering of globular pores is supported by the results of the present paper. JÜzkovÄ R et al: Analysis of porous structure in plasma-sprayed coating 5 so -sola -30-20 ¦ 10-c ^^ ^—^^ ^^ -**L> •**>- Data Minimum *tr -rh^, Maximum 1 O C 3 O O 3 _ *T ID 0.6 - J^V Q 04- ¦ Vzord< — 0 - R M B 5 if S T 1 -D.a ¦ D.6 -D.4-Q.2 ¦ .* ¦r^- i'S' Data L* c 1 \Ci ü en r 3 1Q o id T ID ID Fig. 9. The diagrams indicate values of L—function L(-), the nearest neighbour distribution function D(-) and the spherical contact distribution function HS(-) for pores of diameter equal to 10-20 \im. The starting value of r corresponds to minimum of the theoretical distance between the centroids of the smallest pores (in this case, 10). ACKNOWLEDGEMENTS The research was supported by grants GACR 201/03/0946 and MSM113200008. The autors wish to thank to their referees for helpful suggestions and to O. Roussel, A. Tricoire (E.N.S.I.L. Limoges, France) for their help with metallography and microscopy. REFERENCES Baddeley AJ, Gill RD (1997). Kaplan-Meier estimators of distance distributions for spatial point processes. Ann Stat 25:263-92. Ctibor P, Chra´ska P, Roussel O, Tricoire A (2002). 3D reconstruction of structure of plasma spray deposits. Proceedings from the 1st International Surface Engineering Congress and the 13th IFHTSE Congress, Columbus, Ohio, 357-60. Ctibor P, Roussel O, Tricoire A (2003). Unmelted particles in plasma sprayed coatings. J Eur Ceram Soc 23:2993-9. Moller J, Waagepetersen RP (2003). Statistical Inference and Simulation for Spatial Point Processes. Boca Raton: Chapman & Hall/CRC. Ohser J, Mu¨cklich F (2000). Statistical Analysis of Microstructures in Materials Science. New York: Wiley. Olympus MicroImage 2, Reference Guide (1999). Serra J (1982). Image analysis and mathematical morphology. London: Academic Press. Soille P (1999). Morphological Image Analysis. Berlin: Springer. Stoyan D, Kendall WS, Mecke J (1995). Stochastic geometry and its applications (second edition). New York: Wiley. Tewari A, Gokhale AM (2000). Efficient estimation of number density in opaque material microstructures: the large area disector. J Microsc 200:277-83. 52