Strojniški vestnik - Journal of Mechanical Engineering 66(2020)10, 581-590 © 2020 Journal of Mechanical Engineering. All rights reserved. D0l:10.5545/sv-jme.2020.6883 Original Scientific Paper Received for review: 2020-07-25 Received revised form: 2020-09-14 Accepted for publication: 2020-09-21 Effect of Thermal Barrier Coating on the Thermal Stress of Gas Microturbine Blades and Nozzles Oscar Tenango-Pirin1 - Elva Reynoso-Jardón1 - Juan Carlos García2* -Yahir Mariaca1 - Yuri Sara Hernández3 - Raúl Ñeco1 - Omar Dávalos1 1 Departamento de Ingeniería Industrial y Manufactura, Universidad Autónoma de Ciudad Juárez, México 2 Centro de Investigación en Ingeniería y Ciencias Aplicadas, Universidad Autónoma del Estado de Morelos, México 3 Tecnológico Nacional de México/Campus Pachuca, Pachuca de Soto, Hidalgo, México Thermal barrier coatings play a key roie in the operational life of microturbines because they reduce thermal stress in the turbine components. In this work, numerical computations were carried out to assess new materials developed to be used as a thermal barrier coating for gas turbine blades. The performance of the microturbine components protection is also evaluated. The new materials were 8YSZ, Mg2SiO4, Y3Ce7Ta2O235, and Yb3Ce7Ta2O235. For testing the materials, a 3D gas microturbine model is developed, in which the fluid-structure interaction is solved using CFD and FEM. Temperature fields and stress magnitudes are calculated on the nozzle and blade, and then these are compared with a case in which no thermal barrier is used. Based on these results, the non-uniform temperature distributions are used to compute the stress levels in nozzles and blades. Higher temperature gradients are observed on the nozzle; the maximum temperature magnitudes are observed in the blades. However, it is found that Mg2SiO4 and Y3Ce7Ta2O235 provided better thermal insulation for the turbine components compared with the other evaluated materials. Mg2SiO4 and Y3Ce7Ta2O235 presented the best performance regarding stress and thermal insulation for the microturbine components. Keywords: thermal barrier coating, gas microturbine, turbine blade, thermal stress Highlights • An investigation of the effectiveness of novel ceramics for TBC applications is carried out for their use in gas microturbine blades. • The Mg2SiO4 and Y3Ce7Ta2O23 provided the best performance on thermal insulation under operational environments of the turbine. • The Mg2SiO4 and Y3Ce7Ta2O23 showed similar thermal and mechanical stress magnitudes on the blades, which were the lowest among the materials studied. • The use of those ceramics led diminishing the temperature and stress developed on the blades, which in turn, enables an increase in the operating life of the turbine. 0 INTRODUCTION Gas microturbines (GMT) are small turbomachines that work using gases at high temperatures, with power capacities ranging from 15 kW to 300 kW and offer variable speeds from 30,000 rpm up to 120,000 rpm [1]. They operate with the same operation principle of all the conventional gas turbines; therefore, the efficiency of these devices depends on the gas temperature, which can become higher than 1000 K [2] and [3]. Blades and nozzles of the turbines are subjected to different loads like high temperature, corrosion, centrifugal forces, etc., which could lead to failures [4] and [5]. Regarding high temperature, one of the main drawbacks of GMTs is their small size, which augments heat transfer among their components, leading to failures by burning out or highly stressed zones. Highly stressed zones are often located near the root of the blades as a consequence of non-uniform temperature fields because of sudden changes in geometry and restrictions at the root [6] to [9]. In conventional turbines, in order to cool the blades, internal cooling passages are manufactured; however, this method can not be implemented for microturbines, given the size of such machines [10]. Therefore, the thermal barrier coating (TBC) is used to protect turbine blades and to resist high temperatures environments. A TBC often has a cover that is composed of three layers: the first is a ceramic topcoat (TC) layer, which has direct contact with hot gases; the second is a bond coat (BC) layer, which offers corrosion resistance; the third is a thermally grown oxide (TGO) layer, which is frequently formed between the TC and BC. Ni-based superalloys are often used as a substrate in gas turbines where TBC provides them with thermal insulation. In this way, TBC allows reducing substrate temperature, prolonging the operation life of the turbine and improving turbine efficiency by increasing its operating temperature [11]. *Corr. Author's Address: Universidad Autonoma del Estado de Morelos, Av. Universidad 1001, Cuernavaca, Mor., México, jcgarcia@uaem.mx 581 Strojniški vestnik - Journal of Mechanical Engineering 66(2020)10, 581-590 The design of proper TBC plays a key role in the operating life of turbine blades; as a consequence, several studies are focused on the study of TBC characteristics. Li et al. [12] studied the thickness of TBC for a gas turbine blade. The materials used were ZrO2-8 wt% Y2O3 (8YSZ) as a TC layer, a-Al2O3 as TGO layer and NiCrAlY as a BC layer. TC thickness was varied from 100 ^m to 1000 ^m. It was observed that when increasing the TC layer thickness, the thermal insulation capability and stress levels within the coatings are enhanced. Radwan & Elusta [13] performed one-dimensional (1D) calculations to assessing a four-TBC layer made of Zirconia. They concluded that the TBC layer allowed reducing the blade temperature, thus enhancing the blade durability. Also, the material influenced the temperature distribution. Thickness optimization of a barrier coating of partially stabilized zirconia (PSZ) of a turbine blade was executed by Sankar et al. [6]. TBC thickness was varied from 100 ^m to 800 ^m, and they found a critical thickness that occurred when thickness reached 550 ^m, where heat transfer rate was the lowest. In other research [14], the temperature distribution and thermal stress field of TBC were obtained by employing a 2D decoupled method. TBC included the BC, TGO and TC, where the TC was made of 8YSZ. Non-uniform temperature fields were obtained, and zones with high stress were detected at the suction and the leading edge of the blade. In another study, it was shown that the impact of foreign object damage (FOD) could cause erosion on blade samples with TBC, and was most dangerous as it becomes perpendicular to the surface [5]. Several materials have been developed to be used as TBC and studies have been performed to compare their effectiveness of substrate protection. In a review of the main materials used as TC and TGO presented by Sahith et al. [15], it was concluded that the most common material used as substrate material and bond coat was nickel-based superalloys. Meanwhile, yttria-stabilized zirconia (YSZ) (7 %) was most commonly used as the topcoat. YSZ appears as one of the best options to thermally protect the substrate blade, given its thermo-physical and mechanical properties, such as low thermal conductivity. Regarding the gas microturbines, YSZ was also the preferred material to be used as TBC [16]. In recent studies, some researchers have proposed new materials to be used as TBC for gas turbine applications. Chen et al. [17] proposed the synthesized forsterite-type Mg2SiO4 material as an alternative to zirconia. They showed a comparison with zirconia (8YSZ) in terms of mechanical properties. The new material proved to have a lower thermal conductivity at 1073 K (1.5 W/(mK) and 2.2 W/(mK) for the new material and 8YSZ, respectively) and better thermal-shock resistance than those made of 8YSZ. Other mechanical properties, such as hardness, fracture toughness and Young's modulus were similar to those of zirconia. Shi-min et al. [18] also introduced two novel ceramics for thermal barrier coatings. The proposed synthesized materials, Y3Ce7Ta2O23.5 and Yb3Ce7Ta2O23.5, were new rare-earth tantalite oxides with thermal conductivities lower than that of 8YSZ at 1000 K. Beyond this temperature, the new materials showed good stability which makes them appropriate for high-temperature applications. In another study, Yang et al. [19] synthesized high-purity Dy0 02Gd002 5Yb0.025Y0.05Zr0.88O1.94 (DZX Ti0.02Dy0.02Gd0.025Yb 0.025Y0.05Zr0.86O194 (TZ), and the YSZ powder and coating. According to their results, the TZ TBCs could effectively protect the superalloy substrate at 1873 K. Also, the thermal conductivity of the TZ coating was lower than both DZ and YSZ, showing its potential to be used as TBC. However, most of those new materials have been tested under controlled conditions in a laboratory, and their effectiveness under realistic turbine operating conditions need to be considered. Numerical methods are the preferred tools to predict the thermal and structural fields on blades coated with TBC. Abubakar et al. [20] performed a general review of some methods for predicting residual stress in thermal spray coatings, concluding that some finite element method (FEM) schemes provide results of those stresses reasonably close to experiments; thus, they can be used to predict them in coatings. In contrast, Zhu et al. [7], evaluated the effectiveness of three versions of the k-e turbulence model to predict temperature fields by means of computational fluid dynamics (CFD). The k-e realizable offered the most accurate results when modelling blades with one-layer TBC. Li et al. [12] employed the FEM to design the TBC for a gas turbine blade. They found that thermal insulation was enhanced with the increase of the TC. Also, as mentioned before, Tang et al. [14] carried out a fluid-structure interaction (FSI) coupling method to predict stress fields in a turbine blade. However, their model was restricted to a 2D model. In other works [12] and [21], the uniform temperature on blade surfaces and static blade boundary conditions were used. However, these simplified conditions could drive to imprecise results since temperature fields on blade surfaces, induced by high-temperature combustion gases and the blade, are highly three dimensional and non-uniform. 582 Tenango-Pirin, O. - Reynoso-Jardón, E. - García, J.C. - Mariaca, Y. - Hernández, Y.S. - Ñeco, R. - Dávalos, O. Strojniški vestnik - Journal of Mechanical Engineering 66(2020)10, 581-590 In this work, a 3D FSI decoupling method of a gas microturbine is carried out to assess novel TBC materials proposed in the technical literature for applications in turbine blades. The coating materials studied are taken from literature: 8YSZ, Mg2SiO4, Y3Ce7Ta2O23.5 and Yb3Ce7Ta2O23 5. To accurately predict the non-uniform temperature fields, stationary guide vanes are taken into account for the numerical domain, giving rise to a single axial passage that includes a static guide vane (nozzle), a rotating blade and their corresponding fluid domains. Results in terms of temperature and stress are discussed to identify the most effective thermal barrier coating for the microturbine. 1 METHODS 1.1 Microturbine Characteristics The main components of a gas microturbine engine are the expansion microturbine, the combustion chamber and the compressor. The expansion turbine section is referred to here as the microturbine. Computational modelling of an axial gas microturbine is performed using Ansys (Fluent and Mechanical structural). The microturbine geometry is taken from the technical literature [22], which is designed to supply output power of about 29 kW at the rated speed of 76,000 rpm. The microturbine geometry consists of one stage with 17 guide vanes of the stator (nozzle) and 32 rotating blades. Given its axisymmetric geometry, one 3D single passage that includes one nozzle and one blade is constructed for CFD and FEM computations. Nozzle Interface Fig. 1. Turbine geometry and boundary conditions In Fig. 1, the numeric domain is depicted; the shroud is hidden to provide better visualization of the components. The whole computational domain is integrated by both fluid and solid domains corresponding to air and substrate (and TBC) domains. The dimensions of the turbine are 6.82 mm in nozzle height, 7 mm in blade height, 12 mm in nozzle chord, 8 mm in blade chord, and 77.17 mm in maximum turbine diameter. 1.2 Methodology for the CFD and Mesh Characteristics The numerical domain is discretized to generate the mesh needed for calculations. A mesh dependence analysis is carried out to obtain the optimal mesh for the passage. Given that the nozzle and the rotor domains have meshed separately, an interface boundary condition is used to join both domains. It is noteworthy than both solid and fluid domains been have meshed to predict temperature fields and heat transfer flux. Element sizes of 0.1 mm and 0.2 mm and hexahedral element types are used. A mesh refinement to model the boundary layer was used in the near-wall region of nozzles and blades. Refinement is defined using 5 to 20 layers with a growth rate of 1.2. As a result, nine different meshes with densities ranging from 1,376,014 to 7,627,342 elements are constructed. After the analysis, a mesh with 6,233,945 (Fig. 2) elements is selected to carry out all computations because it has a variation of about 1 % of the computed substrate temperature with respect to a finer mesh. It should be mentioned that the mesh region corresponding to substrates (nozzles and blades) is used for the FEM solution to solve the stress generation. Fig. 2. Computational mesh for flow field and substrate domains 1.3 Boundary Conditions of the CFD and Cases of Study As a consequence of modelling one single passage of the turbine (Fig. 1), boundary conditions are defined as follows: periodic boundary conditions are used at lateral sides of passage, and an interface is created between the nozzle and blade domains to ensure fluid continuity through the passage. The hub and shroud Effect of Thermal Barrier Coating on the Thermal Stress of Gas Microturbine Blades and Nozzles 583 Strojniški vestnik - Journal of Mechanical Engineering 66(2020)10, 581-590 are defined as walls with a no-slip condition. When modelling TBC, thickness and type of material are assigned at coupling walls to model the heat transfer through the coating material. Mass flow rate and pressure outlet conditions are assigned at the inlet and outlet of the passage, respectively. The total mass flow rate incoming to the turbine is equal to 0.23 kg/s at rated speed. Since one single passage is used, this flow rate is divided by the 17 guide vanes resulting in 0.0135 kg/s per passage. The fluid is considered as high-temperature compressible air, and a temperature of 1035 K is also specified at the inlet. The solid domain is defined according to the substrate material (nozzles or blades) and their corresponding TBC. The investigated TBCs are Mg2SiO4, Y3Ce7Ta2O235 and Yb3Ce7Ta2O235, which are compared to 8YSZ. Coating material properties are taken from [17] to [19], and their properties are defined based on a temperature range of 973 K to 1073 K to match with the operating temperature of the microturbine. The TBC properties used to predict heat transfer by CFD are indicated in Table 1, which also contains substrate (nozzles and blades) properties [23]. Coating thickness is considered uniform and constant for the whole cases with a value of 550 ^m based on some investigations [6] and [12], and it is added as a virtual layer into the CFD software. This virtual technique uses the TBC properties (thermal conductivity, specific heat, density and thickness) to compute heat transfer to the blades. Table 1. Top coating and substrate properties Material Thermal conductivity [W/(mK)] Density [kg/m3] Specific heat [J/(kgK)] Nimonic 105 22.23 8010 628 (substrate) [23] 8YSZ [19] 2.2 3610 505 Mg2SiO4 [17] 1.76 3210 177000 Y3Ce7Ta2Ü23.5 [18] 1.78 7245 472.1 Yb3Ce7Ta2O23.5 [18] 1.4 6321 431.1 1.4 Governing Equations and Turbulence Model A numerical FSI analysis is conducted using Ansys Workbench commercial code, which is used to solve thermal (Ansys Fluent) and structural (Ansys-Structural) analyses. First, the steady-state Reynolds-averaged Navier Stokes equations (RANS) approach is used to solve governing equations using CFD. Then the FEM is employed to calculate thermal stresses. CFD is used to compute the flow and temperature fields at steady-state conditions, taking into account the turbine operating conditions, like mass flow, pressure, temperature, wall thermal condition and rotating speed. Governing equations of continuity and momentum, given by Eqs. (1) and (2), are solved using Ansys Fluent. To compute the conjugate heat transfer (CHT), the energy equation, depicted by Eq. (3), is also solved. f),- », dv- P-+ Pv¡ ivi 8s -Pi -PFi = 0 p — + ps¡v + pvt i - T^Vj i + qt i - pr = 0, (1) (2) (3) where t is the time, p is the density per unit mass, vi is the i component of the velocity vector, e is the internal energy per unit mass, Fi is the i component of the body force vector, p is the pressure, qi is the heat flux, r is the heat supply per unit mass and tij is the viscous stress tensor, respectively, and the comma means is partially derived with respect to the independent variables [24]. As stated before, since the air is considered as compressible, continuity and momentum equations are coupled to the energy equation, which provided the temperature distribution in the flow field. Therefore, in order to solve turbulence and accurately predict the CHT, the realizable k - e turbulence model is selected based on the study of Zhu et al. [7]. This model is based on the k - e model [25] and differs in the formulation for the turbulent viscosity and for the transport equation for the dissipation rate as follows: f (pk )+8r( pku 8xj 'gb-ps~ ym 8x- h .h. a 8k k y pC,Ss - pC2 8x- .ËL a„ 8x- 8s (4) 8x- k + -Jvs where C1 = max 0.43, + C, - C^sGb + S s n = S , 5 = (5) -5 k is the turbulent kinetic energy and e is its rate of dissipation. Gk represents the generation of turbulence kinetic energy due to the mean velocity gradients, Gb is the generation of turbulence kinetic energy due to buoyancy, YM represents the contribution of the fluctuating compressible turbulence to the overall 2 584 Tenango-Pirin, O. - Reynoso-Jardón, E. - García, J.C. - Mariaca, Y. - Hernández, Y.S. - Ñeco, R. - Dávalos, O. Strojniški vestnik - Journal of Mechanical Engineering 66(2020)10, 581-590 dissipation rate, C2, C1e and C3e are constants, ak and ae are the turbulent Prandtl numbers for k and e, respectively. Finally, Sk and Se are user-defined source terms. In contrast, to model the rotation of blades, the multiple reference frame (MRF) model is used. This model is suggested when a steady-state approach is used, and it allows to observe the instantaneous flow field of the moving part as to freezing it. Both stationary and moving domains are governed by continuity and momentum equations (Eqs. (1) and (2)). To calculate the heat flux transferred from the flow to solid, a conductive and convective problem is modelled by solving the CHT. Newton's law of cooling and Fourier's law are solved to obtain the energy exchange by convection and conduction mechanisms, respectively. Fourier's law to compute the three-dimensional heat conduction is described by Eq. (6): d2T d 2T d2T q 1 dT —2 +—2 +—T +--> (6) dx dy dz k a dr where q is the heat transfer rate, k is thermal conductivity, a is the thermal diffusivity of the material, T is the temperature, and x, y and z, are coordinate directions. The convection is solved using Newton law which is given by Eq. (7) as follows: q = hA (Tw - T,), (7) where h is the convection heat-transfer coefficient, A is the surface area, and Tw and Tx are the surface and the free stream fluid temperature, respectively. In this case, h is highly influenced by the fluid motion driving to a forced convection case. In contrast, the thickness of the coating is specified on the walls to obtain the thermal insulation and temperature decrement on turbine components. Temperature fields are needed to calculate stress fields in the blades. Mechanisms for corrosion and erosion, phase transformation of the ceramic coating or its fracture behaviour are neglected to simplify the computations of the effectiveness of thermal insulation of TBCs and the thermal stress induced on nozzles and blades. 1.5 Finite Element Procedure Once temperature fields are solved using CFD, FEM is employed to compute stress on the components using Ansys Mechanical software. A decoupling method is used, in which the temperature in the components computed by CFD is sent to FEM software. The surface and internal temperature fields of substrates are used as a boundary condition in conjunction with restrictions to elongation of turbine components. Constraints are specified as follows: the nozzle is fixed at the top and bottom surfaces, and the blade is fixed only at its root. The properties of nozzle and blade are assigned based on Nimonic alloy 105 substrate material, as indicated in Table 1. This material is assumed to be homogeneous, continuous, and isotropic. Other properties used are the Young modulus of 139 GPa, the Poisson's Ratio of 0.3 and the thermal expansion coefficient of 15.3 * 10-6/K [23]. The same solid meshes depicted in Fig. 2 are used to predict three-dimensional stress fields accurately. 1.6 Governing Equations Thermal stress distribution was calculated following Eqs. (8) to (10); non-uniform and three-dimensional stress fields were obtained. Those equations describe three thermal stress components in the tangential, radial, and longitudinal directions, respectively [26]. The equations allow taking into account partial mechanical constraints and internal constraints due to differences in thermal expansion of elements due to different temperatures. Also, when solving them, it was assumed that thermal equilibrium is reached at the rated speed of the turbine. Ea 1 ( 2 r