Numerical Investigation into the Effect of Splats and Pores on the Thermal Fracture of Air Plasma-Sprayed Thermal Barrier Coatings

The effect of splat interfaces on the fracture behavior of air plasma-sprayed thermal barrier coatings (APS-TBC) is analyzed using finite element modeling involving cohesive elements. A multiscale approach is adopted in which the explicitly resolved top coat microstructural features are embedded in a larger domain. Within the computational cell, splat interfaces are modeled as being located on a sinusoidal interface in combination with a random distribution of pores. Parametric studies are conducted for different splat interface waviness, spacing, pore volume fraction and fracture properties of the splat interface. The results are quantified in terms of crack nucleation temperature and total microcrack length. It is found that the amount of cracking in TBCs actually decreases with increased porosity up to a critical volume fraction. In contrast, the presence of splats is always detrimental to the TBC performance. This detrimental effect is reduced for the splat interfaces with high waviness and spacing compared to those with low waviness and spacing. The crack initiation temperature was found to be linearly dependent on the normal fracture properties of the splat interface. Insights derived from the numerical results aid in engineering the microstructure of practical TBC systems for improved resistance against thermal fracture.


Introduction
Thermal barrier coatings (TBC) are heat insulation layers applied on the high-temperature regions of aircraft jet engines in order to protect the crucial structural components against overheating and consequently to extend the lifetime of these components. They also improve the combustion efficiency of the engine by allowing higher operational temperature. The most common TBC system is a multilayered system. It consists of two layers; an outer ceramic layer that protects the substrate from high-temperature gases and an intermediate metallic layer which provides adhesion between the substrate and the ceramic layer. The metallic layer also protects the substrate from high-temperature oxidation and corrosion. The outer ceramic layer is called top coat (TC) and is often made of yttria stabilized zirconia (YSZ). The metallic layer is called bond coat (BC) and is usually made of MCrAlY (M refers to Ni or Al or a mixture of these two. The coating system is subjected to a thermal cycle during operation where the temperature ranges from ambient to its operating value (around 1200 C). During operation, a thin thermally grown oxide (TGO) layer consisting of alumina (Al 2 O 3 ) is formed at the TC/BC interface due to oxidation of the bond coat. During each thermal cycle, thermal stresses are generated in the coating due to a mismatch in coefficients of thermal expansion (CTE) of the TBC layers. These thermal stresses along with the growth of the TGO layer lead to the nucleation and evolution of microcracks in the coating as illustrated in Fig. 1. With the increase in thermal cycles, these microcracks coalesce to form a larger crack which results in final failure of the coating known as spallation.
Conventional TBCs are predominantly manufactured either by an electron beam physical vapor deposition (EB-PVD) process or the air plasma spray (APS) technique (Ref 1). In the EB-PVD technique, the coatings are deposited on the substrate in vacuum using a high-energy electron beam (Ref 2). This technique produces a uniform coating with a fine columnar grain structure. These coatings offer superior strain tolerance, thermal shock and oxidation resistance which ultimately leads to prolonged service life. In the APS method, ceramic particles are melted and sprayed onto the substrate at atmospheric pressure using plasma jet (Ref 3). APS TBCs have a shorter lifetime compared to EB-PVD TBCs due to its layered microstructure and higher number of defects. However, APS remains the preferred manufacturing method in the fields of aerospace and power generation gas turbines because of its low cost and high production efficiency (Ref 4). A schematic diagram of the microstructure of the APS TBC and its failure behavior is presented in Fig. 1.
The distinctive microstructural features of APS TBCs are splats in addition to other features such as pores, microcracks and interface roughness. Splats are building blocks of the APS TBCs which are formed during the solidification process of the sprayed ceramic particles (Ref 5). The interfaces between the splats are usually weak. The characteristics of the splat interfaces are determined by the spraying parameters (Ref [6][7][8]. Literature studies on TBC microstructure show that the presence of splats along with microstructural pores significantly influences the material properties and failure behavior of the TBC (Ref [9][10][11][12][13]. Multiple efforts have also been undertaken in the literature to model and predict the influence of splat interfaces on thermo-elastic properties of TBCs (Ref [14][15][16][17]. In the work ( Ref 14), real microstructural images are used to predict the effective thermal conductivity and stiffness using finite element analysis. The stiffness reduction due to splats is obtained indirectly by calculating the difference between the stiffness values of finite element prediction of as-sprayed (with splats) and experimental measurements of thermal-cycled (no splats) samples, respectively. Similar type of work is also carried out in (Ref 15) using imagebased extended finite element modeling. In addition, the evolution of the stress intensity factor and the quenching stress as a function of intra-splat cracks and substrate temperature is also analyzed in their study. Thermal conductivity of splat interfaces is obtained in (Ref 16) through iterative finite element computations and by mapping the calculated effective thermal conductivity with the experimental measurements. A finite element model is generated using real splat interface distribution obtained from fracture image analysis. A finite element analysis on idealized YSZ model system to analyze the effect of splat interfaces on thermal conductivity of TBC is carried out in (Ref 17).
Literature studies to predict the influence of all aspects of the TBC microstructure on its failure behavior are limited. Some studies are focussed only on the effect on TC/  Ref 11,21,22). For instance, in addition to the microstructural pores, the stacking configuration and the bonding between the splats are also important for the fracture resistance of coatings (Ref 21). It is, therefore, the objective of this research to carry out a systematic numerical study to investigate the effect of interface topology and its fracture strength on the TBC fracture behavior in the presence of microstructural pores. Furthermore, concurrent multiscale modeling setup is adopted to include the free edge effect which is relevant in terms of interpreting experimental results of coated samples. The effect of splat interfaces is investigated using a parametric TBC model with realistic pore distribution corresponding to various splat interface features, namely waviness, spacing and also the strength. The scope of this research is relevant to the TBC community for microstructure optimization to enable them to choose the right set of processing and spraying parameters to get the desired splat interface characteristics.

Multiscale Finite Element Model of TBC Microstructure
Overall Geometry The modeling of the TBC is carried out in a two-dimensional domain under plane strain assumptions. A two-scale concurrent multiscale modeling approach is adopted in which the explicitly resolved TC microstructural features are embedded in a much larger domain of the TBC with an homogenized TC layer, as shown in Fig. 2. TBC microstructural features such as pores, interface roughness and splat interfaces are modeled explicitly in the TC layer of the embedded computational cell as depicted in Fig. 2 The dimensions of the TBC geometry considered in this work are derived from the disk-shaped APS TBC samples developed within the research consortium (Ref 23) for microstructure characterization and thermal cycling behavior. The total thickness and the radius of the diskshaped sample are 6 mm and 15 mm, respectively. The width of the embedded computational cell is considered as 480 lm ( Ref 24). More details about the geometry and the modeling setup can be found in ( Ref 18).
TBC computational cell is enriched with pore features that are representative of TC porosity measurements obtained as part of the collaborative work (Ref 23). For conciseness, only the statistical distribution of the porosity measurements is provided in Table 1. The pores are modeled as an ellipsoidal quantity with an aspect ratio (AR), orientation and size given in Table 1 for modeling convenience. The pores are distributed randomly within the TBC computational cell.
For simplicity, the interface waviness between TC and BC and also the horizontal splat interfaces are represented as a sinusoidal curve. The vertical coordinate y n of the n-th splat interface is given by as a function of the horizontal coordinates x as y n ðxÞ ¼ A 0 C n cosðx 2p=kÞ þ nh ðEq 1Þ where A 0 and k corresponds to the reference amplitude and wavelength of the sinusoidal curve, respectively. The subscript n refers to the types of interfaces in the TC layer (n = 0 corresponds to the TC/BC interface and n [ 0 corresponds to the splat interfaces), h and C denote the spacing between the two adjacent interface and interface waviness, respectively. The waviness parameter C ranges from C ¼ 0 (for a flat interface) to C ¼ 1 (for a sinusoidal interface with amplitude A 0 . The wavelength (k) and initial amplitude (A 0 ) of the interface are assumed to be constant with a specified value of 60 lm and 10 lm, respectively (Ref 24). The vertical splat interfaces are modeled as straight lines, and they are distributed randomly between the horizontal splat interfaces as shown in Fig. 2(c). The number of vertical splat interfaces between any two successive horizontal interfaces is set to a value of 10 over the length of the explicitly modeled domain.

Simulation Setup
The microstructure is meshed with three-noded plane strain triangular elements (CPE3) using GMSH to model the bulk response of the individual TBC layers. The fracture behavior of the TBC is modeled using a cohesive zone model which is implemented in the finite element framework using the zero thickness four-noded cohesive element (COH2D4). It is also important to mention that the fracture process is modeled only in the computational cell and the outer regions of the TBC only can respond as a homogenized featureless elastic medium. Also, the elastic response of the splat interfaces themselves are not considered in the analysis since the interfaces are very thin compared to the thickness of the splats. In other words, the constitutive behavior of the splat interfaces is governed only by its fracture properties, defined through zero thickness cohesive element. To enable arbitrary crack initiation and propagation, the cohesive elements are inserted at every bulk element interface of the TBC computational cell using a MATLAB script ( Ref 25). However, the cracks are allowed to propagate only along inter-element boundaries which lead to a mesh dependency effect. Hence, to obtain a converged fracture pattern, a fine mesh size of 1 lm is used in the computational cell region with a random mesh pattern that mitigates mesh dependency and provides meanvalue convergence. The outer elastic regions of the TBC are coarsely meshed with an average element size of 150 lm to reduce the computational time. The loading condition considered in this work corresponds to a single thermal cycle of as-deposited TBC system. A typical thermal cycle of a TBC consists of three phases namely the heating, dwell and cooling phase. In the heating and dwell phases of the thermal cycle, the TBC is assumed to be stress free as these coatings are deposited at high temperature. Hence, only the cooling phase of the thermal cycle is considered in the analysis, where the temperature is decreased from its typical operating value of 1100 to 30°C. Simulations are carried out under uniform temperature distribution which corresponds to the testing of samples in an oven (uniform temperature). The case of non-uniform temperature distributions has been studied experimentally in (Ref 26). Stresses appear in the TBC system during cooling due to the thermal contraction and the mismatch in coefficients of thermal expansion as shown in Fig. 1. The cooling rate does not play a role as the simulations are based on a quasistatic analysis. The nonlinear fracture simulations are carried out in the FEA package Abaqus using an implicit Newton-Raphson solver.

Constitutive Models and Material Properties
The bulk response of the different layers of TBC is assumed to be linearly elastic and isotropic. The fracture behavior of the TBC is modeled using a bilinear cohesive law (Ref 27) which is governed by the fracture strength and fracture energy. Material parameters used for individual TBC layers are summarized in Table 2. Elastic and thermal properties of the TBC layers are considered to be the same values used in the earlier work ( Ref 18). The fracture properties such as normal strength and toughness of the TC, TGO and BC layers are in correspondence with the values reported in (Ref [28][29][30]. The mode I fracture energies (G IC ) given in Table 2 are calculated from the toughness values (K IC ) under plane strain assumption. The ratio of shear strength to normal strength (c) for the TC layer, TGO layer and splat interfaces is assumed arbitrarily to have a value of 4 to suppress the failure of TBC in planes that are not parallel to the top surface in accordance with experimental observations. The mode I fracture energies G IC reported in Table 2 are calculated from the fracture toughness values. For convenience, the same value of c is also used to define the ratio between mode II and mode I fracture toughness. The effective (reduced) elastic modulus of 80 GPa is used for the outer (homogenized) TC layer. The fracture properties of TC/BC interface are assumed to have the same values as the BC layer. The cohesive stiffness parameter is assigned with a sufficiently large value of 10 10 GPa to eliminate the influence of artificial compliance introduced through the cohesive law.

Results and Discussions
Parametric simulations are carried out to analyze the influence of the TBC microstructure such as splat interfaces and pores on TBC failure behavior. In addition, the effect of splat interface fracture properties is also considered in this study. The results are discussed in terms of crack initiation temperature and total length of all microcracks formed in the computational domain at the end of the cooling cycle for each parametric case. The crack initiation temperature is defined based on a predefined initiation length to avoid a mesh dependency effect. Predefined initiation length is defined as the minimum total length of consecutively failed cohesive elements required to form a measurable microcrack. Based on the convergence analysis carried out in ( Ref 18,24), a predefined length of 3 lm is used to define the crack initiation temperature for all the parametric simulations. The failure of the cohesive element is defined based on energy dissipation and the element is considered to have failed when it dissipates 95% of its total fracture energy. The TC layer is assumed to have completely failed if it becomes detached from the BC layer.

Effect of TBC Microstructural Features
In this section, the effect of individual microstructural features on fracture behavior of the TBC is investigated. Three different TBC models with explicitly represented microstructures: (1) only pores, (2) only splat interfaces and (3) combined splat interfaces and pores are considered for this analysis along with the dense (defect free) TBC. The pores are distributed randomly in the TC layer with features given in Table 1. As discussed in the modeling setup, the horizontal splat interfaces are modeled using Eq 1 with the splat interface waviness and spacing of 0.8 and 12.5 lm, respectively. To include the effect of randomness of the pore distribution and vertical splat interfaces in the analysis, five different realizations are considered. To illustrate the effect of splat interfaces in the TBC failure process, stress distributions (r yy ) along with the crack pattern at the end of the thermal loading step for two different microstructural configurations (namely TC layer with pores and combined splat interfaces and pores) are shown in Fig. 3. The overall stress distribution (r yy ) is in correspondence with the studies carried out in the literature ( Ref 19,31). As observed in Fig. 3, the peaks and valleys of the TC/BC interface lead to tensile and compressive states of stress (r yy ), respectively. The CTE mismatch between the layers and the presence of pores result in a complex stress field especially at the vicinity of pores as shown in Fig. 3. Also, for both configurations, the cracks tend to nucleate and propagate from the edge due to tensile stress (r yy ) generated out of the free edge accompanied by local stress concentration near the pores. However, as shown in Fig. 3(a), (b), the presence of splat interfaces has a significant effect on the TBC cracking behavior both in terms of crack nucleation and propagation. For the TBC with combined splat interfaces and pores, the cracks propagate mostly along the splat interfaces. In case of TBC with only pores, the crack propagation is governed by the direction of local stress concentration generated by the pores. It can also be observed that the effect of the cracking remains relatively local since the stress distribution in the TC sufficiently far away from the cracked regions remains identical for both configurations.
The results are presented in terms of crack initiation temperature and total crack length in Table 3 for the four distinct TBC microstructural configurations as discussed.
It can be observed from Table 3 that the variations in TBC microstructural configuration have a significant influence on its failure behavior. For the microstructural configurations without splat interfaces, the presence of pores did not have a substantial effect on crack initiation temperature compared to the dense case whereas a positive influence of pores is observed for the configurations with splat interfaces (i.e., a lower crack initiation temperature compared to the splat interfaces only case indicates an increased resistance against crack nucleation during cooling). In terms of the total crack length, the presence of pores has a positive effect for both the configurations (with and without splat interfaces). However, this effect is more pronounced for the microstructure without splat interfaces. For instance, the total crack length of the dense TBC is approximately two times higher than the TBC with pores. This is due to the fact that the presence of pores increases the TBC compliance which reduces the strain energy contribution to the crack driving force (Ref 32). The positive influence of the pores on TBC lifetime is reported in the literature through experimental measurements (Ref 33,34). On the other hand, a TC with splat interfaces leads to early crack initiation and also to a larger crack length, which is expected because of the weak interface. From Table 3, it can be concluded that the presence of pores results in an improved resistance of the TBC against fracture whereas the presence of splat interfaces is detrimental to the TBC. To understand the effect of splat interfaces on TBC failure behavior, parametric simulations are carried out in the next section for different geometries and fracture properties of the splat interfaces.

Parametric Simulation
The parametric studies are carried out for three different parameters such as splat interface waviness, spacing and porosity, and the variations are listed in Table 4. The simulation results are discussed in terms of crack initiation temperature and total crack length. Similar to the study on effect of TBC microstructural features, five different realizations are utilized for each parametric case to obtain the  statistical variation due to the random distribution of pores and vertical splat interfaces.

Effect of Splat Interface Waviness
To investigate the effect of splat interface waviness, four different waviness values 0, 0.4, 0.8 and 1 are analyzed. The splat interface spacing and porosity are fixed with a value of 12.5 lm and 10 %, respectively. The waviness value of 0 results in flat interfaces in the computational domain, whereas the waviness of 1 corresponds to the cosine curve with specified amplitude of 10 lm for all the splat interfaces (see Eq 1). For other waviness values, the amplitude of the splats (cosine curve) is scaled linearly with the increase in number of splat interfaces. The results of the simulations are summarized in Fig. 4 to show the influence of splat waviness on initiation temperature and total crack length. From Fig. 4(a), it can be observed that the crack initiation temperature shows a weak dependency on splat interface waviness. This is due to the fact that the crack initiation is governed mostly by the free edge stress which is maximum at a region away from the TC/BC interface. In this region, the change in splat interface amplitude for waviness values less than 1 is minimal because, with increase in number of splat interfaces, there is a decrease in amplitude. For a waviness value of 1, the amplitude of the splat interfaces remain unaffected which results in a trend toward delay in crack initiation as shown in Fig. 4(a).
In terms of total crack length, the splat interface waviness shows a significant influence especially for the waviness of 0 and 1 as shown in Fig. 4(b). The crack length for the TBC with flat splat interfaces (waviness = 0) is sensitive to the distribution of pores. Two out of five realizations show complete failure of the TC with major cracking close to the TC/BC interface. Hence, the statistical variation of total crack length for zero waviness is considered as inconclusive. This behavior is due to the presence of a flat interface which favors the linking and propagation of microcracks. In addition, the formation of microcracks is influenced by the distribution of different pore features. Complete failure of the TC occurs if these microcracks reach a critical total length. For the splat interfaces with waviness of 0.4 and 0.8, the total crack length remains almost the same and it decreases by a factor of 2 for the waviness of 1. Thus, in general, it can be concluded that an increase in splat interface waviness improves the resistance against thermal fracture of the TBC. These findings are in agreement with indentation tests on APS TBC samples (Ref 21) that show that the fracture resistance is indeed improved significantly for TBCs with rough splat interfaces. Similarly, thermal cycling studies (Ref 35, 36) on TBC also conclude that the bond coat roughness has a significant effect on TBC lifetime with higher roughness leading to an increase in TBC lifetime.

Effect of Splat Interface Spacing
Another parameter of interest is the splat interface spacing, which is a direct measure of the splat thickness. In the current TBC modeling setup, splat interface spacing affects both the interface waviness and the number of weak interfaces within the TC layer (refer Eq 1). Hence, the variations in splat interface spacing may have a significant effect on the TBC fracture characteristics. In order to study this effect, four different splat interface spacing given by 6.25, 12.5, 25 and 50 lm are considered. The values for splat interface waviness and volume fraction are kept constant at 0.8 and 10 %, respectively. The results obtained from the analysis are summarized in terms of crack initiation temperature and total crack length and are shown in Fig. 5. As observed in Fig. 5(a), the crack initiation temperature remains almost constant for the range of splat interface spacing considered. Nevertheless, for the reference configuration with only pores (no splats) given in Table 3, the crack initiation temperature is significantly lower than the initiation values reported in Fig. 5(a). In terms of total crack length, the results indicate that the total crack length decreases nonlinearly with the increase in splat interface spacing as shown in Fig. 5(b). For a spacing of 6 lm, the final crack length is highly sensitive to the distribution of pores with three realizations predicting complete TC failure. This inconclusive behavior is governed by the stress distribution close to the TC/BC interface and the local stress concentration generated by the pores. The decrease in total crack length with the splat interface spacing is attributed to the reduced number of weak interfaces and the increase in splat interface waviness.

Effect of Pore Volume Fraction
In this section, the volume fraction of the pores is varied with a constant splat interface waviness of 0.8 and a spacing value of 12.5 lm. Four different pore volume fractions given by V f = 0, 10, 15 and 20 % are considered. Two distinct microstructure configurations are considered in order to illustrate the effect of splat interfaces on the failure behavior of TBC along with porosity. One configuration is created with only pores, and one configuration is formulated with combined splat interfaces and pores. The crack initiation temperature and total crack length for both configurations are reported in Fig. 6. For all porosity values considered, the presence of splat interfaces significantly influences the failure characteristics of the TBC as observed in Fig. 6. The presence of splat interfaces leads to early crack initiation when compared to the TC layer with only pores as shown in Fig. 6(a). The crack initiation temperature remains approximately constant for the TC layer with only pores, whereas the crack initiation temperature decreases nonlinearly with an increase in porosity for the microstructure with combined splat interfaces and pores. For the TC layer with combined splat interfaces and pores, the crack initiation temperature of 10% porosity decreases from 800 to 700°C when compared with the zero porosity case. Upon further increase in porosity, the crack initiation temperature remains almost constant up to V f = 15% after which it decreases to 550°C for V f = 20 %. In terms of total crack length, it can be seen that for both cases (with and without splat interfaces) the trend remains almost the same. The presence of pores improves the fracture resistance of the TC compared to the dense TBC (see Fig. 6b). This is due to the decrease in strain energy contribution to the crack driving force with porosity as explained in the section on the effect of TBC microstructural features. The scatter in total crack length values between individual simulations for a nominally identical sets of input parameters is attributed to the complex interaction of different pore features and also its interaction with the splat interface. Also, it is important to note that complete failure of TBC occurs for a critical value of the porosity is V f = 30% for TC with pores and V f = 20% for TC with combined splat interface and pores. The optimal porosity value remains constant at V f = 15% for both microstructural configurations.

Effect of Fracture Properties
As discussed earlier, the constitutive behavior of the splat interface is governed by its fracture parameters. Hence, the variation of these parameters significantly influences the fracture behavior of the TBC system. To study this effect, a parametric analysis is carried out for distinct material properties, expressed as ratios of a reference benchmark material with fracture strength r TC 0 and fracture energy G TC 0 , as follows: where f is the fracture ratio. Note that both properties are varied using the same fracture ratio (f).
The parametric study is done with fixed splat interface waviness of 0.8, splat interface spacing of 12.5 lm and porosity of 10 %. Two distinct normal (f n ) and shear fracture (f s ) ratios are considered to analyze the influence of normal and shear fracture properties of the splat interfaces. For each fracture ratio, the simulations were carried out for five different values given by 0.5, 0.75, 1, 1.5 and 2. While studying the influence of one fracture ratio, the other fracture ratio is fixed at the value of 1. In other words, to study the effect of splat interface normal fracture properties, the normal fracture ratio (f n ) is varied with values given by 0.5, 0.75, 1, 1.5 and 2 for a fixed shear ratio (f s ) of 1 and vice versa.
Crack initiation temperature and total crack length for both normal and shear fracture ratios are shown in Fig. 7. From Fig. 7(a), it can be observed that the crack initiation temperature decreases almost linearly with increase in normal fracture ratio. For the variations in shear fracture ratio (f s ), the crack initiation temperature remains unaffected as shown in Fig. 7(c). This is because the normal fracture properties, being always less than the shear fracture properties for the splat interface, are the most critical values that control the crack initiation. In terms of crack length, the results indicate that the variation of the total crack length with the normal fracture ratio is not linear as shown in Fig. 7(b). In particular, a drastic increase in damage is observed in the TC layer when the normal fracture ratio is decreased to a value of 1 [i.e., total crack length increases abruptly by a factor of 2.5 when compared to the fracture ratio (f n ) of 1.5]. Similarly, the results shown in Fig. 7(d) are parameter value sensitive particularly, when the shear fracture ratio is decreased below a certain critical value. For instance, the complete failure of the TBC occurs when the shear fracture ratio is decreased below the value of 1. The critical value for the complete TBC failure is decided by the value of the fracture parameters for which the crack coalescence is favoured.

Effect of TGO Thickness
The above reported simulations are carried out for the assprayed TBC system under single thermal cycle. However, during thermal cyclic operation, the oxidation of the bond coat results in growth of an additional layer called the thermally grown oxide layer (TGO). The growth of this additional layer affects both the stress distribution and the crack evolution pattern due to its strong CTE mismatch with the adjacent layers. Complete failure of the TBC occurs when the TGO thickness is in the range of 8 to 12 lm ( Ref 37,38). In this section, the influence of TGO on TBC fracture behavior is analyzed. Due to computational constraints, this effect is studied parametrically for a single cycle. This is based on the assumption that the TBC will undergo cracking only after the TGO thickness is increased to a specified value. Simulations are carried out for four different TGO thicknesses values given by t TGO = 0, 3, 6 and 9 lm. The microstructural features of the TC layer such as porosity, splat interface waviness and spacing are kept constant with a value of 10%, 0.8 and 12.5 lm, respectively. The pores are distributed randomly in the TC layer with pore characteristics specified in Table 1.
The normal stress distribution in the vertical direction (r yy ) and the crack pattern in the TBC are shown in Fig. 8 for two different TGO thicknesses t TGO = 0 lm (left inset) and 9 lm (right inset). From Fig. 8, it can be observed that the presence of TGO significantly influences the stress distribution and cracking pattern in the TBC. The TGO thickness t TGO = 9 lm leads to complete delamination failure of the TBC with majority of cracking on the splat interface close to the TC/TGO boundary, whereas for the TBC without a TGO layer(t TGO = 0 lm) only limited cracking is observed. This is due to the altered TBC stress fields especially close to the TC/TGO interface. The stress concentration in the TBC is governed by two distinct parameters, one due to the presence of pores and the other due to the CTE mismatch between the layers for both cases (with and without TGO layer). However, for the TBC with TGO layer, the stress concentration close to the interface is amplified as shown in Fig. 8. This is attributed to the strong CTE mismatch between the TGO and the other layers (TC and BC) compared to the TBC without TGO layer. Furthermore, the presence of TGO also influences the stress concentration in the vicinity of pores close to the TC/TGO interface. For instance, the vertical normal stress (r yy ) at the proximity of pores close to the TGO layer is largely tensile for TBC with TGO layer, whereas the stress fields are mostly compressive for TBC without TGO layer. The stress concentration close to the TC/TGO interface is further amplified with the TGO thickness. This results in early initiation and propagation of microcracks at the horizontal splat interface close to the TGO layer. Complete failure of the TBC occurs when these microcracks coalesce to form a major delamination crack. It is also important to note that the edge crack in the TC region is minimally influenced by the TGO layer when compared to the cracking in the TC region close to the TC/TGO interface. This is expected as the stress distribution in the TC region sufficiently away from the TC/TGO interface is less affected by the presence of the TGO layer. The influence of the TGO thickness on the TBC fracture performance is given in Fig. 9 in terms of the total crack length. The total crack length is plotted as a function of four considered TGO thicknesses. From Fig. 9, it can be seen that the increase in TGO thickness decreases the mechanical integrity of the TBC system. A TGO thickness of 9 lm results in complete failure. This behavior  is due to the increase in stress concentration with TGO layer thickness as discussed above.

Conclusions
Thermal fracture behavior of the APS TBC with explicitly modeled splat interface and pores is analyzed using a cohesive elements based finite element method. The influence of splat interface geometric characteristics (such as the roughness of the planes in which the splat interfaces are located and their spacing), porosity and splat interface material parameter (fracture properties) on TBC fracture behavior is studied through parametric simulations. The results of the parametric studies are reported in terms of the crack initiation temperature and total microcrack length in the computational domain upon cooling down from the peak temperature of 1100°C to room temperature (30°C), whereby both the splat interface and pores are found to exhibit a significant influence on the fracture behavior. The following conclusions are drawn from the parametric finite element investigations.
1. Microstructural features of the TBCs play a significant role in determining its failure characteristics. With reference to the fully dense TBC system, the presence of pores improves the fracture resistance whereas the presence of splat interface leads to early and more extensive failure of the TBC. 2. The influence of splat interface waviness on the crack initiation temperature is limited. The crack length for the flat splat interface is sensitive to the distribution of pores. In general, the increase in splat interface waviness improves the TBC fracture resistance. 3. Splat interface spacing did not have a significant influence on the initiation temperature, whereas the total crack length decreases with an increase in splat interface spacing. 4. The presence of pores improves the resistance against TBC fracture for both microstructural configurations considered. However, porosity beyond a critical value leads to early failure. 5. On the effect of splat interface fracture properties, the crack initiation temperature scales more or less linearly with normal tensile strength of the splat interfaces. In terms of crack length, the results are sensitive when the fracture ratio values are reduced to a value below a critical level. Complete failure of TBC occurs when shear fracture ratio is decreased below a value of 1. 6. The increase in TGO thickness proportionally decreases the mechanical integrity of the TBC system.
The results and insights as reported here aim to clarify understanding the effect of TBC microstructure especially splat interfaces and pores on failure behavior, which in turn can aid in optimizing the processing and spraying parameters for improved performance/lifetime. Based on the current simulations, TBC systems with an optimal porosity value, having the splat interfaces arranged in non-planar configuration, having a large splat interface spacing and a high splat interface tensile strength perpendicular to the splat interface plane, should give the best thermal cracking resistance.