1. Introduction
Due to the merits of ceramic oxides, such as high thermal, wear and electrical resistances, thermal barrier coating (TBC) composed of ceramic oxides is a common technique to keep the substrate from mechanical wear and thermal shock [
1,
2,
3]. TBC is a multi-layer composite material including topcoat (TC), bond coat (BC), and substrate. For TBC working on cycling thermal shock, it is almost unavoidable to generate thermally grown oxide (TGO) between TC layer and BC layer due to the oxidation of the metallic element in the BC layer [
4,
5]. Furthermore, the thickness of TGO will increase with its working time. The generated TGO will introduce extra stress into TBC. It has been proven that the stress induced by TGO is an important factor of TBC failure. In order to predict the service life of TBC without destroying its structure, scholars use various techniques to research TBC and TGO.
In 2004, Sugasawa et al. established a model to analyze the behavior of bulk ultrasonic waves in sprayed coating layers. Using this model, the sound speed can be determined if the coating thickness is known, and vice versa. According to their reports, the spatial resolution of this method is about tens of micros. Thus, this method should be applicable to measure the thickness of TBC, but not TGO. This is because the typical thickness of TGO is about several micros [
6]. In 2017, Shrestha et al. used thermal wave imaging to evaluate non-uniform coating thickness, with a measurement accuracy of about several microns [
7,
8]. However, they did not demonstrate the ability of this method on the thickness measurement of TGO. In 2019, Wang et al. developed a method to non-destructively evaluate the thickness of TGO by using impedance spectroscopy [
9].
Terahertz time-domain spectroscopy (THz-TDS) has recently emerged as a new nondestructive testing method. Due to excellent penetrability of terahertz (THz) radiation for ceramics, it is a powerful tool for non-destructive analysis of the coatings in TBC [
10,
11]. In 2009, White et al. adopted THz-TDS imaging to map the thickness of yttria stabilized zirconia (YSZ) in TBC, and thermal degradation was observed [
12]. In 2010, Chen et al. observed the relevance between the phase delay of THz pulse and the SEM-determined TGO thickness [
13]. In 2014, Fukuchi et al. measured the refractive index and thickness of the TC layer in TBC using reflective THz radiation [
14].
Although this excellent work has demonstrated the measurement capability of THz radiation on the structure of TBC, the accurate thickness of the TGO layer has not been steadily measured by the THz method until now. Based on the previous work, to obtain more accurate information regarding the TGO layer, this paper analyses the experimental signal and the simulated signal, combines with SEM results, and speculates on the thickness and refractive index of the TGO layer. The rest of this paper is organized as follows.
Section 2 introduces the measurement principle and experimental process. The experimental results and simulated results are provided in
Section 3.
Section 4 presents the conclusions.
2. Methods
In this work, a typical all-fiber reflective THz-TDS system was used. The working parameters ((Menlo Systems, Munich, Germany)manufacturer, city, state (only for USA and Canada), country) of the femtosecond laser adopted in the THz-TDS system were as follows: repetition rate 100 MHz (±1 MHz), wavelength 780 nm (±30 nm), power more than 330 mw, pulse width less than 100 fs. Two photoconductive antennas were used to emit and receive THz radiation. The spot diameter of the THz radiation was about 5 mm, and its reliable frequency was in the range 0.3 THz~1.5 THz. The experimental temperature was controlled at 20 °C. A double-sided polished high resistance silicon was used as a beam splitter.
Figure 1 shows the schematic diagram of the experimental setup.
Electron Beam Physical Vapor Deposition (EBPVD) technology was used to prepare the TBC samples. The material of the TC layer in this work was 8YSZ (ZrO2-8 wt.%Y2O3), and nickel-based alloy which contains nickel, chromium, aluminum, yttrium, and hafnium composed the BC layer. To prepare the samples containing TGO layer, the initial TBC samples were subjected to 500, 1000, 2000, 5000 thermal cycles. As the time of thermal cycle increases, the TGO layer continues growing. This newly generated TGO layer changes TBC from the initial three-layer structure to a four-layer structure.
The schematic diagram of the THz propagation path in a sample is shown in
Figure 2 [
15]. Both of the distances between
R0 and
R1,
R1 and
R2 mean the transmission process of THz radiation in TC layer, so that the
between them are equivalent. The formula between thickness and the refractive index of TC layer can be expressed as
where
represents the refractive index,
denotes the thickness,
denotes the speed of the light in vacuum and
denotes the time delay of THz radiation in the TC layer.
Figure 3 shows the initial and after thermal cycling multi-layer structures of the TBC samples, and the different propagation path of THz radiation in the samples. The samples used in this paper are prepared in the same batch. TC layers are considered to have the same thickness. Therefore, to get the information regarding the TGO layer with high reliability, the thickness of the TC layer is firstly determined in advance using the as-prepared samples.
As shown in
Figure 3a, in the initial TBC, the THz radiation emitted from the emitting antenna is represented as
E0. The signal obtained after reflection on the surface of the TBC is represented as
R0.
R1 and
R2 are the first and second reflected signals from the interface between TC layer and BC layer, respectively.
Similarly, in
Figure 3b,
R0 represents the surface reflection of TBC after thermal cycle.
E1 is the reflected signal from the interface between TC and TGO layer.
R1 represents the interface reflection of TGO and BC layer, and
R2 represents the second reflection which travels through the ceramic layers twice.
According to Fresnel formula, we have
and
where
and
(
,
= 0, 1, 2, 3) are the Fresnel transmissive and reflective coefficients, respectively, of each interface and 0 represents air, 1 the TC layer, 2 the BC layer (in
Figure 3a) and the TGO layer (in
Figure 3b), 3 the BC layer (in
Figure 3b). Let
(
= 1, 2) represent the transmission coefficient of each layer, where 1 represents the TC layer, and 2 the TGO layer (in
Figure 3b). Further,
, where
denotes the absorption coefficient per unit length in the sample,
denotes the thickness.
Using the transmissive and reflective coefficients, the signal from the as-prepared sample can be expressed as follows
The signals such as the signal reflected from the interface between TC and TGO is too small, so that it is difficult to distinguish the interface signal of TGO/BC from the noise (a conclusion derived from experimental results). Meanwhile, because of the relatively thin thickness of TGO layer, the time delay between
E1 and
R1 is too close. Therefore, the primary and secondary reflection of the interface are mainly studied in this paper. The expressions of the subsequent signals are omitted. Finally, the expressions of the samples after the thermal cycles are as follows
The optimization algorithm adopted in this paper to obtain the minimum error between the experimental signal and the simulation signal is defined as follows
where
represents the sum of the squares of the difference between experimental results and simulation results,
represents the experimental results, and
represents the simulation results. By the optimization search on the error function, the thickness and refractive indexes of TC and TGO can be determined. Here, we have obtained the simulated signals of the as-prepared sample and the sample after thermal cycle, respectively, and then collected the corresponding experimental signals.
By selecting the suitable range of refractive index and thickness, the mentioned previously can be calculated by the algorithm. After comparing the different , the minimum can be obtained and the corresponding indexes are considered to be the final result.
3. Results and Discussion
In this work, to determine the thickness and refractive index of the TC and TGO layers, the experimental measurement is carried out in two steps. The as-prepared sample is measured at first to obtain the thickness and refractive index of the TC layer. Then, indexes obtained by last step are used as the known quantities for the samples after thermal cycle. Repeating the calculating process in samples after thermal cycle, the relevant indexes of TGO layer can be finally acquired.
In order to simulate the reflective signals from the samples, an air signal of the THz-TDS system is needed as the reference signal, as shown in
Figure 4a. Then, the samples are fixed on the sample holder. By adjusting the position of holder, the focus of the THz radiation is just on the surface of the sample. A captured THz signal of the as-prepared sample and the simulated result are shown in the
Figure 4e.
Figure 4f shows the difference between simulation and measurement.
Figure 4b–d shows the simulated signals
R0,
R1 and
R2, whose combination obtains the closest fitting to the experimentally captured reflected signal from the as-prepared sample. The closest fitting is obtained when the refractive index takes 4.45, and the thickness of the TC layer takes 95 μm [
15]. As indicated in
Figure 4e,f, the simulated signal is well fitted with the experimental signal. In addition, from
Figure 4e, we know that TC layer is thick enough so that the peaks from the different reflected signals are distinguishable.
To ensure the accuracy of the thickness acquired by the algorithm, all the samples used in this paper are cut for the SEM observation. HITACHI S-4800 Field Emission Scanning Electron Microscope (Hitachi Limited, Tokyo, Japan) is used for microscopic observation. The microscope image of the as-prepared sample is shown in
Figure 5 that shows the thickness of the TC layer is about 94.1 μm. The SEM result matches the results obtained by the proposed THz method well. Each sample was tested five times, and the obtained results were analyzed with the SEM measurement results. Meanwhile, to minimize the error, the experimental signals are averaged after 100 scans. Finally, the measurement standard deviation of the as-prepared sample is 0.6893.
Figure 6 shows the experimental and simulated signals from the TBC samples after 1000 thermal cycles. As shown in
Figure 6c, the signal of
E1 represents the reflection from the TC/TGO interface. Because the refractive indices of the TC and TGO layers are relatively close to each other, the Fresnel reflective coefficient between them is low. This is the reason why the amplitude of
E1 is small. Therefore, the second reflected signal from the TC/TGO interface is not taken in account in the theoretical model of this work. The closest model fitting result is obtained when the refractive index and the thickness of TGO layers takes 4.54 and 4.71 μm, respectively.
Other thermal cycling samples are also analyzed using the above experimental process, and the optimization algorithm is also used in the samples to obtain the refractive index and thickness. Meanwhile, to verify the refractive index, the SEM results are put into the optimization algorithm. The results are listed in
Table 1. As shown in
Table 1, the maximum relative error of the refractive index is not more than 5%, which shows a good consistency. The error may be caused by the experimental environment, the experimental process and the algorithm itself. THz radiation is easily affected by the water vapor in the air which may lead to the occurrence of scattering. A slight shift in the position of the experimental platform during sample replacement will also bring in some measurement errors.
Meanwhile, the comparison between THz measurement and SEM observation regarding the thickness of the TGO layer is shown in
Figure 7. Both results indicate that the increasing rate of the TGO thickness decreases gradually with the increase of the thermal cycling times. The reason is that in the initial stage of the thermal cycle, the contact area between BC layer and external environment is still relatively large. When the alumina grows to a certain extent, the contact area between the BC layer and air shrinks which causes the reaction speed to slow down [
16]. The refractive indices of the generated TGO layer are almost constant, as shown in
Table 1.
Figure 8 presents the SEM images of those samples.
Figure 9 shows the distribution of the
in Equation (6). As shown in
Figure 9, it can be seen that the distribution is continuous, smooth, and has only one local minimum. It is proved that the
has a unique convergence value. Therefore, the global minimum of the
which represents the best fitting result can be found. The proposed algorithm to determine the thickness and the refractive index should be robust.
4. Conclusions
In this work, a reflective THz-TDS system is used to obtain the reference signal and experimental signal by measuring air and the TBC samples with different times of thermal cycle, respectively. Then, a model which shows the transmission of THz radiation in the TBCs is established. After simulating the transmission, the simulation signals of the samples before and after thermal cycle have been obtained. By comparing the simulated signals and experimental signals, an optimization algorithm which can get two kinds of signal interpolation squares is developed to evaluate the quality of the simulated signal. The minimum of the optimization algorithm means the simulation results attain the best effect. At this time, the refractive index and thickness in the simulation results are considered to be the real refractive index and thickness of TGO. To verify the validity of the algorithm, all samples are detected by the SEM. Before the SEM experiment, to get a better revivification of the working environment of the TBC, a wire-electrode cutting is used to cut the samples to expose the internal structure. By comparing the SEM results and experimental results, the accuracy of the algorithm is confirmed. Finally, the thickness of the TGO layer can be obtained by the optimization algorithm. Because of its nondestructive and nonionizing features, THz radiation is expected to be used in further research regarding TBC, such as predicting the service life and evaluating failure.