Shiyuan Zhang, Beilin Zhengand Pengfei He

(1. Department of Civil Engineering, Fujian University of Technology, Fuzhou 350108, Fujian, China;2. School of Aerospace Engineering and Applied Mechanics, Tongji University, Shanghai 200092, China)

Abstract: Temperature filed, thermal stress, especially tensile stress and J-integral are important for thermal barrier coatings (TBCs) under thermal shock. At the micro- and nano-scale, the energy transport mechanisms are significantly different from those at the macro-scale. The temperature fields, which are obtained by combining the Equation of Phonon Radiative Transport (EPRT) (for the nano-scale ceramic TBCs)and the Fourier law (for the substrate), are used as the thermal loading in the thermal stress and J-integral of an edge in the TBCs analysis by the finite element method. The temperature field and thermal stresses as well as J-integral are compared with those which are calculated by applying the Fourier law to both the TBCs and the substrate. The influence of the physical heat properties of the TBCs on the temperature field and thermal stress and J-integral have been analyzed in this paper. It is concluded that the temperature, thermal stress, including the tensile and compressive components, and J-integral which are calculated with the EPRT, are lower than that calculated with the Fourier law in the TBCs. Moreover, thermal stress in the TBCs increase with increasing phonon speed and relaxation time, but J-integral at the crack tip is in the opposite

Keywords: nano-ceramic TBCs; EPRT; the Fourier law; thermal stress

0 Introduction

Ceramic thermal barrier coatings (TBCs) have been applied to metal substrates to make them resistant to thermal shock in the aerospace and electronic devices, due to its many well-known advantageous properties, such as insulation and corrosion resistance. The accurate prediction of heat propagation in these devices may be important in ensuring their efficient work and design. Its heat transfer and thermal stress properties have been analyzed by the Fourier law, which is traditionally used to analyze the behavior of macroscopic heat transport. Recently, as nanotechnology and ultrafast processes continue to develop, TBCs with thicknesses between 10 nm and 100 nm are being used extensively in devices (Zeng and Chen[1]). The average free path of the energy carrier has been reached (microscale effect in size), resulting in the discontinuity of heat conduction because of the insufficient number of collisions of the energy-required heat conduction. In addition to microscale effect in time, the high-speed thermal shock reaches femtosecond or picosecond order[2], which is shortened into the range of the heat carrier mean free time, resulting in the premise of Fourier heat conduction theorem, in which the heat conduction velocity is infinite, thus not established. Therefore, there are obvious differences between the heat conduction theory at the microscopic scale and those at the macroscopic scale. New research indicates that the heat transfer mechanics at the micro-and nano-scales may deviate significantly from those at the macro-scale. For example, by the experiment of short-pulse laser heating on metals, Qiu and Tien[3]have verified that the actual heat flux is not as fast as the Fourier law prediction.

Many scholars have studied and published non-Fourier heat transport models that are based on different physical theories. Maxwell[4]first reported that the heat propagation speed is limited and modified the Fourier law by adding a relaxation time. Landau[5]proved heat waves in liquid helium Ⅱ propagate with a finite speed. The CV model by modifying the Fourier law to interpret the thermal wave phenomenon was recently applied by some researchers (Cattaneo[6]and Vernotte[7]). Guyer and Krumhansl[8]brought forward the phonon scattering model by directly solving the linearized Boltzmann equation for heat transport by phonon collision and scattering. This model was improved later by Joseph and Preziosi[9], who emphasized the interrelation with a Jeffrey’s type of heat flux equation. Tzou[10]proposed the dual-phase-lag model with an emphasis on the lagging behavior at the micro-scale. Joshi and Majumdar[11-12]solved the linearized Boltzmann transport equation to derive the EPRT and analyzed heat transfer in dielectric thin films. The EPRT is valid for microscale in both time and size. Qiu et al.[13-14]proposed a two-step model with a quantum mechanical and statistical basis for a short-pulsed laser heating a metal film. Choudhuri[15]established the three-phase delay hyperbolic heat conduction equation (TPL) by joining the phase delay of heat flow, temperature gradient and thermal displacement gradient.

At present, lots of studies researched the objects in the fields of biology and semiconductor. Li[16]compared the analytical and numerical solutions of non-Fourier temperature field of laminated plate and laminated cylinder, and proved that the numerical solution obtained by backward difference separation is absolutely stable. Atfi and Talaee[17-18]studied the non-Fourier thermal response of finite hollow cylinder subjected to periodic thermal shock based on C-V model. Cao[19]researched the dynamic response of piezoelectric pole under thermal shock of heat source based on fractional order theory. Yu[20]solved the temperature field and stress of the laminated material. On this basis, he studied the stress distribution at the crack tip under thermal shock. Li[21]analyzed the fracture behavior of solid plates under thermal shock by the non-Fourier. Zhang and Li[22-23]analyzed transient thermal stress and intensity factors for a circumferential crack in a hollow cylinder by generalized fractional heat conduction. Guo[24]analyzed the fracture mechanics of plate, sandwich plate, coating, cylinder and other structures with the two-phase delayed model. Zhang and Chen[25]studied the temperature field, displacement field and stress intensity factor of cracked hollow cylinder under thermal shock by using single-phase delayed non-Fourier heat conduction mode. Mao et al.[26]studied silicon crystal with point defect scatter by molecular dynamics. Those researches are mainly aimed at only on the TBCs rather than the TBCs structure (thermal barrier coating and substrate), and kingdom is mainly in heat conduction.

In this paper, the equations combining the EPRT (for the TBCs) and the Fourier law (for the substrate) are numerically solved by the temperature field under thermal shock. The thermal stresses in the ceramic TBCs and J-integral of crack at the ceramic TBCs surface are studied by the finite element method, and the results are compared with those obtained by the Fourier law on both the TBCs and the substrate. The influence of the physical heat properties of the TBCs (such as the relaxation time and phonon speed) on the temperature field and the thermal stress as well as J-integral at an edge crack tip is analyzed.

1 Calculated Temperature Field and Stress and J-Integral

1.1 Heat Transportation Theoretical Model

Thermal shock in the direction of thex-axis has been studied (Fig. 1).

Because the thickness size of the ceramic TBCs is comparable with the mean free path(MFP) of typical phonons, which are considered to be the major heat carriers in the dielectric material, the phonon intensity for the ceramic TBCs ((0≤x≤b)) is numerically predicted by the EPRT[11]:

(1)

Fig.1 Structure of the nano-scale ceramic TBCs and the substrate

Once, the phonon intensity is numerically solved by Eqs. (1), the temperature can be calculated by the Bose-Einstein distribution function at an equilibrium state[11]:

(2)

Because the thickness of the substrate (b≤x≤b+h) is much larger than the MFP, the temperature field of the substrate is predicted by the Fourier law[25]:

(3)

whereαis the thermal diffusivity.

1.2 Initial Conditions and Boundary Conditions

Att=0, the ceramic TBCs and the substrate are assumed to be at a uniform temperatureT0and remain adiathermal.

Att>0, the temperature atx=0 is instantaneously raised toT1.For EPRT, the intensity function for incoming phonon remains constant, and for Fourier’s law, the boundary condition is constant temperatureT1.An idealized thermal shock on the surface of the TBCs is assumed, i.e., the heat transfer coefficient is infinitely large. Under this condition, the thermal stress is the greatest.

Because the substrate material is a metal or an alloy, the phonons are absolutely reflected at the bottom of the TBCs.

One assumes that there is no heat loss, and the two materials are in perfect thermal contact at the interface of the ceramic TBCs and substrate, and at the bottom of the substrate the temperature is a constantT0.

The above conditions of the phonon and heat transport are expressed as follows:

(4)

1.3 Temperature Field Numerical Calculation

Eq. (1) are two wave equations in the integrodifferential form. Eq. (1) are numerically solved by an explicit upstream differencing method, i.e., by backward differencing when the phonon intensity wave is moving in the positivex-direction (the first equation in Eq.(1)) and by forward differencing when it is moving in the negativex-direction (the second equation in Eq. (1)). To guarantee the stability of the numerical computation, the forward difference in time is applied while ensuring that Δt≤Δx/v|μ| in Eq. (1) and Eq. (3). The range of the cosine of the angleμ(-1≤μ≤1) is divided intoKdivisions. The integral in Eq. (1) can be approximated by a Gaussian quadrature[12]as

(5)

whereμishows discrete directions andwiis the weighting factor. The integrodifferential equations can be converted into the differential equation of the unknownIi. In Eq. (3), a complete implicit difference method is used to calculate temperature field. Eq. (1) and Eq. (3) can be represented in the finite difference expressions as follows:

(6)

where the superscriptnrepresentsnth time step, subscriptiis the index of the unit direction vector, subscriptjrepresents thejth node,A=vΔt/Δx,B=Δt/τ,Dis the last node of TBCs, andLis the last node of the substrate. Material parameters[27]in the TBCs and substrate are shown in Table 1.

1.4 Thermal Stress and Fracture Calculation by Finite Element

In this work, the thermal stress and J-integral at the crack tip in the TBCs are analyzed using the finite element method after numerical calculation of the temperature field. The numerical calculations make use of the element software MSC.MARC, which can be applied in the thermal stress and J-integral[28]. The stress and temperature analyses are uncoupled. The TBCs and the substrate are modeled by a four-node plane strain and quadratic elements in the stress analysis, because the model is a good elastic body without plastic deformation. In fracture mechanics analysis, it is assumed that the fracture is an edge crack that opens to the surface and its length is 2×10-8m. Quarter point 8-node isoparametric degenerate singular elements are employed in modeling the region near the crack tip, as shown in Fig.2.

Table 1 Parameters and material properties relevant to heat transport

Fig.2 Crack at the ceramic TBCs

The assumed boundary conditions put constraints on the displacement field in they-axis. The displacements only take place along thex-axis, but at the bottom of substrate, thex-displacement is assumed to be zero. The initial temperature should beT0at all nodes. The thermal loads are derived from temperature fields. The parameters and the material properties used in the heat transport, thermal stress, and J-integral numerical calculations are listed in Tables 1 and 2[29-32]. The parameters and properties remain constant unless otherwise specified in the numerical calculations and subsequent analysis.

2 Results and Discussion

2.1 Temperature Distribution and Comparison

Figs. 3 and 4 show temperature field from EPRT and the Fourier law, respectively. It is clear that the temperature jumps on the TBCs surface are calculated by the EPRT, which are not shown by the Fourier law. In the TBCs and substrate, the amplitude of temperature attenuation obtained from the EPRT is lower than the attenuation from the Fourier law. Because the heat transfer velocity is finite according to the EPRT, where the heat transport velocity is represented by income and out-coming of phonon, it is infinite for the Fourier law.

Table 2 Material properties relevant to the analysis

Fig. 3 Temperature field obtained from EPRT

Fig. 4 Temperature field obtained from Fourier heat conduction model

2.2 Results of Thermal Stress and J-integral

The comparison between the thermal stress in they-direction obtained from the EPRT and that obtained from the Fourier law is presented in Fig.5. Both of the thermal stress profiles along thex-direction are similar, i.e., the thermal stresses gradually change from tensile stress to compressive stress along thex-direction in the TBCs and are only compressive stresses along thex-direction within the substrate. But the thermal stress, special tensile stress in the TBCs calculated with the EPRT, is lower than that calculated from the Fourier law, because temperature contrast from the surface of the TBCs to the bottom of the substrate, which is obtained from the EPRT, is lower than that obtained from the Fourier law. At the interface, a thermal stress discontinuity occurs for the difference in the characters between the two materials.

Fig. 5 Comparison of the thermal stress in the y-axis obtained from the EPRT with that obtained from the Fourier law along the x-axis

For reasons given above, the amplitude of temperature attenuation from the EPRT is lower than the attenuation from the Fourier law, and the J-integral from the EPRT is lower than the results from the Fourier law in the TBCs and changes slowlier with time, and J-integral from the Fourier law decreases with time in Fig.6, because temperature difference between the top of TCBs and bottom of substrate from EPRT changes slowlier with time than that from the Fourier law.

Fig. 6 Comparison of the J-integral obtained from the EPRT with that obtained from the Fourier

3 Influence of the TBCs Heat Physical Properties on Temperature, Thermal Stress and J-Integral

3.1 Influence of Phonon Speed

The temperature field, thermal stress profiles and J-Integral at the crack tip predicted by the EPRT with different phonon speeds (v1=8×102m/s,v2=1.3804×103m/s,v3=5×103m/s) are shown in Figs. 3, 7, 8, 9 and 10.

Fig.7 Change of temperature fields with the phonon speed v1

Fig. 8 Change of temperature fields with the phonon speed v3

At the surface of the TBCs, the amplitude of temperature attenuation decreases with increasing phonon speed. However, in the TBCs, the amplitude of temperature attenuation increases with increasing phonon speed, because the intensity function of incoming phonon remains constant at the surface of the TBCs for the EPRT, and the intensity of out-coming phonon increases with increasing phonon speed, so the temperature jump increases, whereas in the TBCs, the phonon transport increases with increasing phonon speed, so the heat transport increases.

Fig. 9 Thermal stress profiles predicted by the EPRT for different phonon speeds

Fig.10 J-integral predicted by the EPRT for different phonon speeds

At the surface of the TBCs, the magnitude of the thermal stress decreases as the phonon speed increases. But in the TBCs, the thermal stress rises with the rise of phonon speed. The J-integral at the crack tip decreases with increasing phonon speed.

3.2 Influence of Relaxation Time

The temperature field, thermal stress profile and J-Integral at the crack tip predicted by the EPRT with different relaxation times (τ1=3×10-12s,τ2=6×10-12s,τ3=1.2×10-11s) are shown in Figs.3, 11, 12, 13 and 14.

At the surface of the TBCs, the amplitude of temperature attenuation descends with the increase of relaxation time. But the amplitude of temperature attenuation increases with increasing relaxation time in the TBCs.

On the surface of the TBCs, the thermal stress decreases with the increase of relaxation time. Nevertheless, in the TBCs, the thermal stress and compressive stress increase for the same reason as stated above. The J-integral at the crack tip decreases with increasing relaxation time.

Fig. 11 Change of temperature fields with relaxation time τ1

Fig.12 Chang of temperature fields with relaxation time τ3

Fig.13 Thermal stress profiles predicted by the EPRT for different relaxation times

Fig.14 J-integral predicted by the EPRT for different relaxation times

4 Conclusions

This study numerically analyzed temperature field, thermal stress and J-integral obtained from the EPRT and from the Fourier law. The following conclusions can be drawn:

1)For nano-ceramic TBCs under thermal shock, the temperature field, thermal stress and J-integral derived from the EPRT are lower than those derived from the Fourier law.

2)The temperature attenuation and thermal stress decrease on the surface of the TBCs with increasing phonon speed and relaxation time, but the amplitude of temperature attenuation and thermal stress in the TBCs increase with increasing phonon speed and relaxation time.

3)J-integral at the crack tip decreases with increasing phonon speed and relaxation time.