Optimization Design of Lattice Structures in Internal Cooling Channel of Turbine Blade

: Recently, the inlet temperatures in gas turbine units have been drastically increased, which extremely affects the lifespan of gas turbine blades. Traditional cooling structures greatly improve the high temperature resistance of the blade; however, these structures scarcely concern both heat transfer and mechanical performances. Lattice structure (LS) can realize these requirements because of its characteristics of light weight, high strength, and porosity. Although the topology of LS is complex, it can be manufactured with the 3D metal printing technology. In this study, an integral optimization method of lattice cooling structure, used at the trailing edge of turbine blades, concerned with heat transfer and mechanical performance, was presented. Firstly, functions between the ﬁrst-order natural frequency ( freq 1), elasticity modulus ( E ), relative density ( ρ ), and Nusselt number ( Nu ), and the geometric variables of pyramid type LS (PLS) and X-type LS (XLS) were established, and the reliability of these functions was veriﬁed. Then, a mathematical optimization model was developed based on these functions which contained two selected optimization problems. Finally, relations among objectives were analyzed; inﬂuence law of geometric variables to objectives were discussed, and the accuracy of the optimal LS was proved by experiment and numerical simulation. The optimization results suggest that, compared to the initial LS, Nu increases by 24.1% and ρ decreases by 31% in the optimal LS of the ﬁrst selected problem, and the Nu increases by 28.8% while freq 1 and ρ are almost unchanged in the optimal LS of the second selected problem compared to the initial LS. This study may provide a guidance for functions integration design of lattice cooling structures used at turbine blades based on 3D printing.


Introduction
The increase of turbine inlet temperature can improve the efficiency of turbomachinery and cause additional damage to turbine blades. To reach higher inlet temperature and protect turbine blades, many cooling structures have been developed, such as lattice structure (LS). For decades, great attention has been paid to LS because of its unique characteristics of high porosity [1,2], ultralight weight [3,4], high specific stiffness and strength [5], and high heat transfer [3,5]. Although traditional cooling structures of turbine blades can greatly improve the cooling efficiency, Han et al. [6] found that the flow resistance coefficient increases rapidly with the improvement of heat transfer capacity in different degrees, and the heat transfer coefficient distribution is inhomogeneous, causing excessive structural thermal stress. In order to solve these issues, Kim et al. [7] studied the heat transfer performance of pyramid-type lattice structure (PLS) by numerical simulation and experiments, and the results showed that PLS had good homogeneity of heat transfer distribution. Yun et al. [8] carried out an experiment to study face-centered cubic LS and found that this LS showed a good thermal performance compared to non-LS channel. Due to the complex topology of LS, metal additive manufacturing technology can easily solve this limit [9,10].
Numerous studies have been conducted on the mechanical performances of LS, including relative density and specific strength [11,12], effective modulus and Poisson's ratio [13,14], yield strength, buckling strength, and plateau strength [15]. These studies greatly enriched the mechanical theory for LS, supporting the mechanical optimization of LS. Hoang et al. [16] investigated LS with direct multiscale topology optimization method which provided an inexpensive and efficient method to obtain optimal LS. Zhang et al. [17] proposed topology optimization of LS to obtain optimal density distribution of LS. The results showed that the compressive and bending stiffness and energy absorption of the optimal LS improved. Wu et al. [18] provided a novel strengthening method of LS and presented a pyramidal LS with some mechanical performance advantages. Other optimization approaches for LS included semi-automated optimization approach [19] and inverse identification approach [20], which greatly improved the mechanical performances of LS. These studies almost focus on mechanical optimization with single objective for LS. The optimization studies on heat transfer or specific mechanical performance (lightweight, vibration, and bearing load capacity), or both, performances of LS are as follows.
Liu et al. [21] proposed a mathematical optimization method with two objectives of LS, which improved the structural performance of each panel under multiple loading cases and minimized structural mass of LS simultaneously. Valdevit et al. [22,23] optimized the design of various weight, mechanical, and heat transfer performance of LS, and they found that hexagon structure had the optimal comprehensive performance. Chen et al. [24] found that sandwich structure filled with sponge had the strongest absorbing capacity and reached the targets of minimum dimensionless mass and maximum absorbing capacity of LS. Roper et al. [25] investigated and analyzed the PLS and found that compression modulus, compressive strength, and maximum heat flux increased with the increase of density. Su et al. [26] presented a multiparameter optimization approach for lightweight fiber-reinforced polymer (FRP) composite triangular LS under nonlinear structural response constraints and achieved a significant improvement in terms of weight saving. Moon et al. [27] made a comparison among Kagome type LS, PLS, and diamond type LS for high bearing load and lightweight characteristics, and came to the conclusion that Kagome LS had the best performance, which was applied to the wings of unmanned aerial vehicles (UAV). Meng et al. [28] optimized the strength, vibration, and bird impact load of LS and realized the lightweight design of LS, which was applied to blades. Fazilati et al. [29] optimized the hexagonal multilayer honeycomb lattice and reached the targets of maximizing energy absorption, and minimizing impact vibration and overall size. The results showed that the performance of the optimized structure was significantly improved compared to that of the single layer of nonoptimized structure. Xu et al. [30] optimized PLS using a genetic algorithm to reduce weight and improve sound insulation performance, and also studied the effects of materials and upper and lower panel thicknesses on performance. Bailong et al. [31] carried out an optimal design of body-centered cubic (BCC) LS, considering the relative density, initial stiffness, and plastic failure strength. The research showed that, compared to BCC structure, the elastic modulus and yield limit of the optimized structure increased by 121.5% and 77.3%, respectively. Smardzewski et al. [32] provided the design of optimized or near-optimized LS by experimentally estimating the mechanical strengths and failure mechanisms and found that with the increase of inclination angle of struts and relative density of LS, the mechanical properties of beams increased. Zhao et al. [33] proposed a concurrent optimization method to optimize natural frequencies of additive manufacturing-fabricated LS, which increased the first natural frequency (freq1) of the optimized structure, 38.5% higher than that of the original structure. Akihiro et al. [34] used a simple basic lattice shape composed of pillars, and only optimized its density distribution by setting the pillar diameter as design variable, and treating steady-state pressure and temperature reductions as multiobjective functions. Adil et al. [35] investigated the optimal designs of novel LS filled with thin tubes based on multiobjective crashworthiness optimization procedure. Lattice member diameter and tube thickness were set as design variables, and minimizing peak crash force and maximizing specific energy absorption were chosen as design objectives. The research revealed that the optimized BCC LS hybrid designs generally had a better crashworthiness performance than BCC LS with vertical strut counterparts with the same design objectives.
As outlined above, optimization methods of LS with two or more objectives of mechanical or heat transfer, or both, greatly expanded the application of LS in heat pipes and tubes. However, few studies on optimization method of LS concerned both heat transfer and mechanical performances used at region "1" of turbine blades, as shown in Figure 1. Although, in [36], Gao et al. developed an effective optimization method of LS concerning both heat transfer and mechanical performances, the application background was for composite structure instead of turbine blades. Therefore, the aim of this study was to construct an optimization model of PLS and X-type LS (XLS) to obtain the optimal LS used at the trailing edge of turbine blades based on functions establishment between geometric parameters (diameter and inclination angle of LS) and representative parameters (Nusselt number, relative density, the first-order natural frequency, and elastic modulus). The boundary conditions were simplified from turbine blades, such as thin-wall structure; applying a constant heat flux on the target surface; mass flow inlet; and pressure outlet. The results show that the heat transfer and mechanical (lightweight, vibration, and bearing load) performances of the optimal LS are better than that of the initial LS obtained from [36] i.e., Nu increases by 24.1% and ρ decreases by 31% in the optimal LS of the first selected problem; and Nu increases by 28.8% while freq1 and ρ are almost unchanged in the optimal LS of the second selected problem, compared to initial LS.
al. [35] investigated the optimal designs of novel LS filled with thin tubes based on multiobjective crashworthiness optimization procedure. Lattice member diameter and tube thickness were set as design variables, and minimizing peak crash force and maximizing specific energy absorption were chosen as design objectives. The research revealed that the optimized BCC LS hybrid designs generally had a better crashworthiness performance than BCC LS with vertical strut counterparts with the same design objectives.
As outlined above, optimization methods of LS with two or more objectives of mechanical or heat transfer, or both, greatly expanded the application of LS in heat pipes and tubes. However, few studies on optimization method of LS concerned both heat transfer and mechanical performances used at region "1" of turbine blades, as shown in Figure 1. Although, in [36], Gao et al. developed an effective optimization method of LS concerning both heat transfer and mechanical performances, the application background was for composite structure instead of turbine blades. Therefore, the aim of this study was to construct an optimization model of PLS and X-type LS (XLS) to obtain the optimal LS used at the trailing edge of turbine blades based on functions establishment between geometric parameters (diameter and inclination angle of LS) and representative parameters (Nusselt number, relative density, the first-order natural frequency, and elastic modulus). The boundary conditions were simplified from turbine blades, such as thin-wall structure; applying a constant heat flux on the target surface; mass flow inlet; and pressure outlet. The results show that the heat transfer and mechanical (lightweight, vibration, and bearing load) performances of the optimal LS are better than that of the initial LS obtained from [36] i.e., Nu increases by 24.1% and ρ decreases by 31% in the optimal LS of the first selected problem; and Nu increases by 28.8% while freq1 and ρ are almost unchanged in the optimal LS of the second selected problem, compared to initial LS.

Array of cylinders for heat transfer augmentation
Internal turbulence cooling and impingement cooling Showerhead

Selected Model
The schematic of cooling arrangement in a turbine blade is shown in Figure 1. In regions "1" and "2" are set some traditional internal cooling structures that could improve cooling efficiency and protect the blade from high temperature; however, these structures scarcely concern both heat transfer and mechanical performances. Herein, considering better performances of heat transfer and mechanics, LS was set at thin walls of the blade.

Selected Model
The schematic of cooling arrangement in a turbine blade is shown in Figure 1. In regions "1" and "2" are set some traditional internal cooling structures that could improve cooling efficiency and protect the blade from high temperature; however, these structures scarcely concern both heat transfer and mechanical performances. Herein, considering better performances of heat transfer and mechanics, LS was set at thin walls of the blade. Some new internal cooling channels would be established at the thin walls of the blade and LS was placed upon these channels as a component of cooling structure. PLS and XLS used herein are shown in Figure 2. Boundary conditions of LS were given as follows: In both PLS and XLS, the height of the cooling channel (H) was considered as a constant in order to simplify optimization calculation. The left and right surfaces were velocity inlet and pressure outlet, respectively. The bottom wall was applied with a constant heat flux which corresponded to the outer heat source surface of the thin wall of the blade, and was called target surface. Other surfaces of PLS and XLS were adiabatic. Fundamental material of PLS and XLS was nickel-base superalloy IN718 with the density of 8240 kg/m 3 , the Young's modulus of 199.9 GPa, and the Poisson's ratio of 0.3 (Table 1).
Some new internal cooling channels would be established at the thin walls of the blade and LS was placed upon these channels as a component of cooling structure. PLS and XLS used herein are shown in Figure 2. Boundary conditions of LS were given as follows: In both PLS and XLS, the height of the cooling channel (H) was considered as a constant in order to simplify optimization calculation. The left and right surfaces were velocity inlet and pressure outlet, respectively. The bottom wall was applied with a constant heat flux which corresponded to the outer heat source surface of the thin wall of the blade, and was called target surface. Other surfaces of PLS and XLS were adiabatic. Fundamental material of PLS and XLS was nickel-base superalloy IN718 with the density of 8240 kg/m 3 , the Young's modulus of 199.9 GPa, and the Poisson's ratio of 0.3 (Table 1). The size of the cooling channel was 80 × 80 × 40 mm 3 . Diameter (D) and inclination angle (ω) of optimal LS were geometric variables. According to the study shown in reference [36], when ω of PLS was selected as 45°, the heat transfer performance and mechanical properties of cooling structure were relatively optimum. Therefore, parameters of nonoptimal PLS and XLS (initial LS) were selected as D = 2 mm and ω = 45°, as shown in Table 1.

The Function of Heat Transfer Performance
In the analysis of cooling performance, boundary condition of a constant heat flux was employed [37], which was consistent with the boundary condition of the target surface of the LS selected. In order to characterize the heat transfer performance of LS, the Nu formula was derived according to the fin method [38].
Nu is used to evaluate heat transfer performance. Functional relationship between geometric parameters of LS and Nu is established: The size of the cooling channel was 80 × 80 × 40 mm 3 . Diameter (D) and inclination angle (ω) of optimal LS were geometric variables. According to the study shown in reference [36], when ω of PLS was selected as 45 • , the heat transfer performance and mechanical properties of cooling structure were relatively optimum. Therefore, parameters of nonoptimal PLS and XLS (initial LS) were selected as D = 2 mm and ω = 45 • , as shown in Table 1.

The Function of Heat Transfer Performance
In the analysis of cooling performance, boundary condition of a constant heat flux was employed [37], which was consistent with the boundary condition of the target surface of the LS selected. In order to characterize the heat transfer performance of LS, the Nu formula was derived according to the fin method [38].
Nu is used to evaluate heat transfer performance. Functional relationship between geometric parameters of LS and Nu is established:

of 23
where h is the overall convective heat transfer coefficient of structure, D h is hydraulic diameter, k f is the thermal conductivity of air, q is the average heat flux of the target surface, T w is the average temperature of the target surface, and T in is the average inlet temperature of air. Heat transfer of each member in LS channels was considered as the same, so a single member in LS was analyzed. The temperature control equation is as follows: where η is the coordinates along member direction; n is the number of LS cores along the flow direction in LS channel; T f (η,n) is the temperature of the η coordinate in member of the nth row of LS core; T f (n) is the average temperature of air in the nth row of LS core; h ∞ is the heat transfer coefficient between member surface, and according to the research in literature [39], h ∞ = 1.19(Re) 0.6 ; P is the section perimeter of the member; S is the cross-sectional area of the member; k s is the thermal conductivity of the member. Initial heat flux (q 0 ) flowed into member from the target surface is calculated: where q 0 is the heat flux coming from the member in the place of η = 0. By combining Equations (4) and (5): where l is the length of a single member; l = H/ sin ω, and H is the height of the LS channel; m = 2 h ∞ /(k s D), with the dimension of L −1 . When η = 0, the temperature of the member is equal to that of the wall, i.e., T s (η,n) = T f (0,n) = T w (n), according to Equation (6): From the law of conservation of energy, the following function is obtained: where ρ f is air density, u m is the average velocity of air, u m = Reµ/ρ f D h ; Re is Reynolds number, µ is dynamic viscosity of air, and C P is the specific heat capacity of air. In order to get q 0 , it is necessary to obtain heat flux q 1 which enters internal fluid directly from the target surface: where h 1 is the heat transfer coefficient between air and target surface. In order to simplify the calculation, thermal resistance between member and target surface is ignored, so h 1 and h ∞ are considered to be the same. Therefore, the heat flux q entering the nth row of LS core is: Appl. Sci. 2021, 11, 5838 6 of 23 By combining Equations (9) and (10), the following equation can be obtained: In conclusion: By combining Equations (7), (11), and (12), the final equations can be obtained. The formula of Nu for a single core of PLS is: The formula of Nu for a single core of XLS is: The function of average Nu for a whole LS channel is: where N was 5, which meant that the whole LS channel had five LS cores. f 1 (x) was the first objective. This function was validated as follows. Material used in the simulation was IN 718, as shown in Table 1. It can be observed from Figure 3 that with the increase of Re, the trend for the growth of Nu was consistent between the experiment value of Kim [38] and the calculated value of PLS channel with the fin method. Meanwhile, the maximum error between the calculated value with the fin method and the experiment value of Kim [38] was approximately 17%. Therefore, the given function of Nu was reliable.
Furthermore, compared with 5 LS cores, when 10 LS cores were set, the average Nu decreased about 23.1%; when 15 LS cores were set, the Nu decreased about 47.9%. Therefore, five LS cells were set in flow direction of LS channel to reduce the overall size and improve heat transfer performance of the cooling structure, as shown in Figure 2. Generally, the pressure drop of the cooling channel gradually increases with the increase of Re. However, when Re reached up to 10,000, the pressure drop was only around 300 Pa. In addition, the changing D and ω of LS had little effect on the changing of pressure drop of LS channel. Thus, pressure drop was not further investigated.

Re
Experimental value of Kim [38] Calculated value with the fin method  [38] and calculated value with the fin method.

The Function of Vibration Performance
Low-frequency vibration and vibration causing structural resonance of turbine blades should be avoided. The natural frequency of turbine blades should be increased. To achieve these objectives, LS cores were set at internal cooling channel of turbine blades, and the finite element analysis method (FEA) was used to set up the relevance between the first-order natural frequency (freq1) and geometric parameters of LS channel. This relevance was determined by an implicit function box by numerical calculation using ABAQUS software.
Therefore, the function of freq1 can be expressed as: This value was given by ABAQUS directly by numerical calculation. The boundary condition of LS channel was fixed support on one side and free on the other side. Herein, vibration of PLS channel and XLS channel were mainly investigated. However, the XLS had symmetrical characteristic in topology. Therefore, the value of Equation (16) was based on the XLS channel as well as this value validation. f2(x) was the second objective, and this function was validated as follows.
The XLS channel is shown in Figure 4a. Material used for the printed XLS channel was PLA, as shown in Table 1. The modified quadric tetrahedron element (C3D10M) was used in the process of mesh generation. After the verification of grid independence, the first six-order natural frequencies (freq1-6) tended to be stable when the grid number was approximately 100,000. Therefore, the computational grid number of the XLS channel was approximately 100,000. Relevant experiments were carried out to further verify the reliability of the finite element results. In the experiment, a SIEMENS SCM205 dynamic signal analyzer was used; the model of force hammer was PCB086C03 with a sensitivity of 2.25 mV/N; a PCB356A66 piezoelectric accelerometer sensor was used with a sensitivity of 1.034 mV/(m·s 2 ), and a vice was used to fix the test sample so as to reduce its impact on the accuracy of the experimental results. Figure 4a shows the sample of XLS channel and the layout of the relevant measuring points. The sample was divided into five parts in the length direction and two parts in the width direction with equal proportion. The nodes were numbered in sequence, and the geometric model, established in LMS Test.Lab, was consistent with those settings. Acceleration sensors were set at points 5, 11, and 14 of the sample to pick up vibration signals, and a hammer was used to knock at point 17. The  [38] and calculated value with the fin method.

The Function of Vibration Performance
Low-frequency vibration and vibration causing structural resonance of turbine blades should be avoided. The natural frequency of turbine blades should be increased. To achieve these objectives, LS cores were set at internal cooling channel of turbine blades, and the finite element analysis method (FEA) was used to set up the relevance between the firstorder natural frequency (freq1) and geometric parameters of LS channel. This relevance was determined by an implicit function box by numerical calculation using ABAQUS software.
Therefore, the function of freq1 can be expressed as: This value was given by ABAQUS directly by numerical calculation. The boundary condition of LS channel was fixed support on one side and free on the other side. Herein, vibration of PLS channel and XLS channel were mainly investigated. However, the XLS had symmetrical characteristic in topology. Therefore, the value of Equation (16) was based on the XLS channel as well as this value validation. f 2 (x) was the second objective, and this function was validated as follows.
The XLS channel is shown in Figure 4a. Material used for the printed XLS channel was PLA, as shown in Table 1. The modified quadric tetrahedron element (C3D10M) was used in the process of mesh generation. After the verification of grid independence, the first six-order natural frequencies (freq1-6) tended to be stable when the grid number was approximately 100,000. Therefore, the computational grid number of the XLS channel was approximately 100,000. Relevant experiments were carried out to further verify the reliability of the finite element results. In the experiment, a SIEMENS SCM205 dynamic signal analyzer was used; the model of force hammer was PCB086C03 with a sensitivity of 2.25 mV/N; a PCB356A66 piezoelectric accelerometer sensor was used with a sensitivity of 1.034 mV/(m·s 2 ), and a vice was used to fix the test sample so as to reduce its impact on the accuracy of the experimental results. Figure 4a shows the sample of XLS channel and the layout of the relevant measuring points. The sample was divided into five parts in the length direction and two parts in the width direction with equal proportion. The nodes were numbered in sequence, and the geometric model, established in LMS Test.Lab, was consistent with those settings. Acceleration sensors were set at points 5, 11, and 14 of the sample to pick up vibration signals, and a hammer was used to knock at point 17. The average value was obtained after striking the point five times. The signals were captured by acceleration sensor and force hammer sensor and were transmitted to multi-channel data acquisition apparatus for processing. Then, the data were transmitted to PC end for relevant data processing, and the relevant parameters were obtained. average value was obtained after striking the point five times. The signals were captured by acceleration sensor and force hammer sensor and were transmitted to multi-channel data acquisition apparatus for processing. Then, the data were transmitted to PC end for relevant data processing, and the relevant parameters were obtained.  Figure 4b shows that in freq1-3, the experimental value and calculated value showed good consistency, of which the maximum error was within 10.5%. However, when it came to freq4-6, the error between them was relatively large, and the maximum error reached 24.9%. The reason for the analysis error lay in two aspects. On the one hand, due to the limitation of printing accuracy of the printer, the printed sample had certain defects. Particularly, the printing mode was adopted without internal support at the top of the channel, which made the resin material on the top of the sample not fit tightly enough. On the other hand, because the sample was a thin-walled structure, it was inevitable that some deformation would occur in the sample fixing with the vice. However, in the simulation process using FEA, there were no defects in the three-dimensional model, and its deformation was not considered when the constraints were applied. Meanwhile, the maximum errors between the modal experimental values and the simulated values caused by some defects of the test piece were 17.7% and 15.21%, respectively [40,41]. Therefore, it can be considered that the error between the experimental value and the simulation value will not greatly affect the subsequent simulation results. Additionally, Figure 4b shows that the freq1-6 of IN 718 was basically linear with that of the resin material. Therefore, the experimental results can provide guidance for research of LS channels made by IN 718. Figure 4c presents the mode shapes of modal experiment and simulation results. It showed a good consistency between modal experiment result and simulation result. Therefore, the calculation model of vibration function established by FEA is reliable.  shows that in freq1-3, the experimental value and calculated value showed good consistency, of which the maximum error was within 10.5%. However, when it came to freq4-6, the error between them was relatively large, and the maximum error reached 24.9%. The reason for the analysis error lay in two aspects. On the one hand, due to the limitation of printing accuracy of the printer, the printed sample had certain defects. Particularly, the printing mode was adopted without internal support at the top of the channel, which made the resin material on the top of the sample not fit tightly enough. On the other hand, because the sample was a thin-walled structure, it was inevitable that some deformation would occur in the sample fixing with the vice. However, in the simulation process using FEA, there were no defects in the three-dimensional model, and its deformation was not considered when the constraints were applied. Meanwhile, the maximum errors between the modal experimental values and the simulated values caused by some defects of the test piece were 17.7% and 15.21%, respectively [40,41]. Therefore, it can be considered that the error between the experimental value and the simulation value will not greatly affect the subsequent simulation results. Additionally, Figure 4b shows that the freq1-6 of IN 718 was basically linear with that of the resin material. Therefore, the experimental results can provide guidance for research of LS channels made by IN 718. Figure 4c presents the mode shapes of modal experiment and simulation results. It showed a good consistency between modal experiment result and simulation result. Therefore, the calculation model of vibration function established by FEA is reliable.

The Function of Bearing Load Performance
Bearing load performance was improved in turbine blades using LS cores. Elastic modulus (E) of the LS channel was used to evaluate its bearing load capacity. The relevance Appl. Sci. 2021, 11, 5838 9 of 23 between E and geometric parameters of LS core was derived using PLS core, and one member of PLS core was analyzed, as shown in Figure 5. An outer load F was added at the upper end of the member, and the member had a displacement δ, which caused axial force F A and shear force F S . As all forces in the member were balanced, F can be obtained along the vertical direction: According to the compression and bending function of the member, F A and F S can be obtained: where E c is the elastic modulus of the material of PLS channel. The equivalent plane compressive stress and strain can be defined as: where l is the length of the member, as shown in Figure 5, l = H/sin ω. A is the area of bottom surface of the PLS core, A = √ 2l cos ω + 2t 2 ; D is the diameter of the member; t is a small spacing between two members at the top wall of the PLS core and also is a small spacing between one of the side walls of the PLS core and the member at the bottom wall of the core.
Combining Equations (20) and (21), the plane elastic modulus of PLS can be obtained: where ρ P is the relative density of the PLS core, and ρ P = πD 2 / sin √ 2l cos ω + 2t 2 . As shown in Figure 5b, D/(2l) is generally a small value; therefore, for Equation (22), the second part can be ignored.
Thus, the equivalent elastic modulus can be expressed as: The function of E can be obtained as: where f 3 (x) is the third objective of the present study. Equation (24) is also used for the E value of XLS channel. The Euler buckling critical load of the member is: where k is the end restraint factor for buckling. For LS, due to the large slenderness ratio of the member and the thin panel, the constraint at the end of the member was weak. Therefore, k was taken as 2, and f 3 (x) was not used as an objective of the present study, while it was used as a constraint in order to simplify the optimization calculation.
For LS, due to the large slenderness ratio of the member and the thin panel, the constraint at the end of the member was weak. Therefore, k was taken as 2, and f3(x) was not used as an objective of the present study, while it was used as a constraint in order to simplify the optimization calculation.

The Function of Lightweight Performance
Lightweight performance was considered to be improved in turbine blades using LS cores. Relative density ( ρ ) of the LS channel was used to evaluate the structural characteristic of being lightweight. The relevance between ρ and geometric parameters of PLS channel can be expressed as: where Vmembers is the volume of the member, Vcore is the volume of the single PLS core, and the function of relative density for the XLS core was , and Sy is the length of the XLS core.

The Mathematical Optimization Model and Two Selected Optimization Problems
Previous studies of LS were mainly on the application of LS's mechanical performances, and few studies used LS in turbine blades to improve both thermal and mechanical performances of blades. We simplified the application of LS on turbine blades as an LS channel with variable aspect ratio and developed a corresponding optimization method. As shown in Figure 1, to set more LS on thin walls of regions "1" and "2", bearing load capacity, vibration, and heat transfer performances of turbine blade would be better. However, it would enlarge weight of turbine blades and pressure drop of the channel. Therefore, the optimization method aimed to seek a reasonable topology of LS which can improve both heat transfer and mechanical performances. Herein, based on the given functions which represented relevance between heat transfer and mechanical performances and geometric parameters of the LS channel, a mathematical optimization model was established; two selected optimization problems were proposed; NSGA-II algorithm was used to solve these problems, and optimal structural parameters were obtained.

The Mathematical Optimization Model
An LS channel may improve both heat transfer and mechanical performances of turbine blades. Relevance between heat transfer and mechanical performances and geometric

The Function of Lightweight Performance
Lightweight performance was considered to be improved in turbine blades using LS cores. Relative density (ρ) of the LS channel was used to evaluate the structural characteristic of being lightweight. The relevance between ρ and geometric parameters of PLS channel can be expressed as: where V members is the volume of the member, V core is the volume of the single PLS core, and the function of relative density for the XLS core was ρ = 3πD 2 / 4S 2 y sin ω , where sin ω = H/(2l), and S y is the length of the XLS core.

The Mathematical Optimization Model and Two Selected Optimization Problems
Previous studies of LS were mainly on the application of LS's mechanical performances, and few studies used LS in turbine blades to improve both thermal and mechanical performances of blades. We simplified the application of LS on turbine blades as an LS channel with variable aspect ratio and developed a corresponding optimization method. As shown in Figure 1, to set more LS on thin walls of regions "1" and "2", bearing load capacity, vibration, and heat transfer performances of turbine blade would be better. However, it would enlarge weight of turbine blades and pressure drop of the channel. Therefore, the optimization method aimed to seek a reasonable topology of LS which can improve both heat transfer and mechanical performances. Herein, based on the given functions which represented relevance between heat transfer and mechanical performances and geometric parameters of the LS channel, a mathematical optimization model was established; two selected optimization problems were proposed; NSGA-II algorithm was used to solve these problems, and optimal structural parameters were obtained.

The Mathematical Optimization Model
An LS channel may improve both heat transfer and mechanical performances of turbine blades. Relevance between heat transfer and mechanical performances and geometric parameters of a LS channel was established, and their functions were obtained. Based on these functions, functional integration design and the multidisciplinary performance requirements were taken into account, and a mathematical optimization design model was established at the start of structural design.
Diameter (D) and inclination angle (ω) of the LS channel were taken as geometric variables; Nu, freq1, E, and ρ were selected as objectives to characterize the performance of heat transfer, vibration, load bearing, and light weight, respectively, and the corresponding mathematical optimization model was developed as follows.
Variables of the model are two geometric variables of the LS channel, which are: Three functions of heat transfer and mechanical performance of mathematical optimization model are given as: Constraints of the variables of the model are: Additionally, its equivalent elastic modulus represented bearing load capacity of LS channel, and the value of equivalent elastic modulus was taken greater than the value of its crushing strength, which was also a constraint for LS. The function is expressed as:

Two Selected Optimization Problems
According to the different emphases of specific engineering applications on turbine blades, two kinds of combinatorial optimization problems were proposed. The optimization problem I (Op-I) was a problem based on Nu and ρ, focusing on the optimization of structural performance of heat transfer and light weight of LS channel. The corresponding mathematical model is as follows: Op-I: where "Max" and "Min" mean that the structural performances of heat transfer and light weight of the LS channel should reach relatively maximum and minimum, respectively, and "s.t." refers to the constraint of the optimization problem. The optimization problem II (Op-II) was a problem based on Nu, ρ, and freq1, focusing on the optimization of the structural performances of heat transfer, light weight, and vibration. The corresponding mathematical model is as follows: Op-II: where "Max" and "Min" mean that the structural performances of heat transfer, light weight, and vibration of LS channel should reach relatively maximum and minimum, respectively.

Establishment and Procedure of Solution for Mathematical Optimization Model
The ISIGHT2018 optimization software was used to build the mathematical optimization model for PLS and XLS channel. As shown in Figure 6, models of PLS and XLS core were established by SOLIDWORKS; value of freq1 was calculated by ABAQUS; and Nu, relative density (ρ), and E were directly calculated by the given functions. Herein, we mainly investigated the heat transfer and mechanical performances of the LS channel. Moreover, in future study, pressure loss, porosity, and friction coefficient of the LS channel can be added at this mathematical optimization model.

Establishment of Approximate Model
In ISIGHT, model objectives (freq1, Nu and ρ) were the input values. Before the optimization work, the relevance between model variables (D and ω) and these model objectives was established by the third-order response surface model (RSM), as shown in Figure 6a, and RSM was used to establish the approximate model to improve the calculation efficiency. At the same time, in order to make all the test points evenly distributed in the design space, the DOE optimal Latin hypercube (opt. LHD) sampling mode was selected, with the sample number of 50, and the sample number for error analysis was 25.
The RSM cannot guarantee to pass all the sample points; therefore, there was a certain error compared to the actual model. In order to ensure the calculation accuracy of the approximate model, the fitting accuracy R 2 was used to verify the accuracy of the approximate model ( Table 2).
The calculation formula of fitting accuracy R 2 is: where SS R represents the sum of regression squares, SS R = ∑ k l=1 (ŷ l − y) 2 ; SS T represents the sum of total squares, SS T = ∑ k i=1 (ŷ i − y) 2 ; y is the average value of response;ŷ l is the predicted value at the design point; y i is the real value of response, and k is the number of sample points. It can be observed from Table 2 that the fitting accuracy values for the four optimization objectives were all above 0.97. Therefore, the model can well approximate the actual model, and this model can be used for the following optimization work of the LS channel.

Optimization Procedure
As shown in Figure 6a, the NSGA-II algorithm, with good exploration performance, was selected for optimization. The specific procedure is shown in Figure 6b and described as follows:

1.
At the beginning of the solution, an accurate approximate model was established with RSM.

2.
The population was initialized. Value of the objective function was solved.

3.
Pareto sorting was carried out to obtain its elite individuals. 4.
Step 3 was iterated repeatedly until the optimal solution set was found.

5.
At the end of the solution, the optimum results were obtained.

NSGA-II Algorithm Setting and Stability Verification
As mentioned above, the NSGA-II algorithm was used to solve the mathematical optimization model of the LS channel, and the PLS channel was used to describe this algorithm as follows: The specific parameters of this algorithm were set as sizes of population, 80; algebras of heredity, 100; crossover of probability, 0.9; and index for distribution of mutation, 20. Figure 7 reveals that when the sizes of population increased from 12 to 24, the Pareto front tended to be consistent. When the sizes of population were 20, the Pareto fronts of PLS and XLS were basically stable, as shown in Figure 7a. With the increase of algebras of heredity from 20 to 40, the Pareto frontier tended to be consistent. When the algebras of heredity were 30, the Pareto fronts of PLS and XLS were basically stable, as shown in Figure 7b. Therefore, to avoid only obtaining the local optimal solution, and to improve the accuracy of the optimization results, the sizes of population and the algebras of heredity were selected as 20 and 30, respectively.
The specific parameters of this algorithm were set as sizes of population, 80; algebras of heredity, 100; crossover of probability, 0.9; and index for distribution of mutation, 20. Figure 7 reveals that when the sizes of population increased from 12 to 24, the Pareto front tended to be consistent. When the sizes of population were 20, the Pareto fronts of PLS and XLS were basically stable, as shown in Figure 7a. With the increase of algebras of heredity from 20 to 40, the Pareto frontier tended to be consistent. When the algebras of heredity were 30, the Pareto fronts of PLS and XLS were basically stable, as shown in Figure 7b. Therefore, to avoid only obtaining the local optimal solution, and to improve the accuracy of the optimization results, the sizes of population and the algebras of heredity were selected as 20 and 30, respectively.

Results and Analysis
In this study, functions representing relevance between heat transfer and mechanical performances and geometric parameters were given. Based on these functions, a mathematical optimization model was built. The given functions were input values of the mathematical optimization model; two selected optimization problems were solved by NSGA-II algorithm, and results are shown as follows.

Sensitivity Analysis
The two selected optimization problems had two and three objectives, respectively, and the values of these objectives were expected to be relatively maximum or minimum.

Results and Analysis
In this study, functions representing relevance between heat transfer and mechanical performances and geometric parameters were given. Based on these functions, a mathematical optimization model was built. The given functions were input values of the mathematical optimization model; two selected optimization problems were solved by NSGA-II algorithm, and results are shown as follows.

Sensitivity Analysis
The two selected optimization problems had two and three objectives, respectively, and the values of these objectives were expected to be relatively maximum or minimum. Therefore, the objective of this study was to reveal the law of influence between the objectives and geometric parameters and to obtain a relatively optimal LS structure. The sensitivity between them was analyzed. The overall sensitivity reveals the relations among objectives, and a higher overall sensitivity means that the solution set of the optimization problem is more concentrative. Single sensitivity reveals the law of influence of geometric variables to optimization objectives, and a higher single sensitivity means that change of variables has a larger influence on objectives. Figure 7a shows that both for PLS or XLS, the Pareto front was an exponential distribution, approximately. Nu of XLS was always bigger than that of PLS, under the same ρ, as XLS was regarded as two PLS composed together. Thus, the heat conduction of XLS was better than PLS. Tables 3 and 4 show the distribution for the Pareto solutions of Nu and ρ for Op-I, respectively. It can be found that the Pareto front was mainly concentrated in the region of higher Nu both for PLS and XLS. For PLS, the optimal Pareto solutions of Nu between 130 and 148 were 73.45% of the total set of solutions of Nu. For XLS, the optimal Pareto solutions of Nu between 150 and 172 were 79.1% of the total set of solutions of Nu. It can also be found that the set of Pareto solutions of ρ were relatively uniform in both PLS and XLS. Therefore, Op-I is more sensitive to Nu than to ρ. As discussed above, under the same ρ, Nu of XLS was greater than Nu of PLS, so only XLS was selected for the analysis of single sensitivity.

Overall and Single Sensitivity Analysis for the Results of Op-I
As shown in Figure 7b, the optimal solutions of D concentrated between 2.3 mm and 2.9 mm, and the total distribution was relatively uniform. Figure 8c shows that the optimal solution of ω was between 52 • and 66 • , which was mainly concentrated between 58 • and 66 • . The diffusion level of ω was bigger than D. Therefore, D had a greater impact on Op-I than ω. In addition, the correlation coefficient r was used to represent the single sensitivity of each parameter, and its calculation formula is as follows: At the same time, the larger these values, the stronger is the correlation between two variables, and the stronger is the single sensitivity, and vice versa. Figure 8d shows that the set of Pareto optimization solutions of ω was the most sensitive for Nu, while the set of Pareto optimization solutions of D was the most sensitive for ρ.
In conclusion, the overall performance of Op-I was that the objective function of Nu had the highest sensitivity for Op-I, ω had the highest correlation or sensitivity for Nu, and D had the highest correlation or sensitivity for ρ.
Among them, the positive sign indicates a positive correlation, and the negative sign indicates a negative correlation. At the same time, the larger these values, the stronger is the correlation between two variables, and the stronger is the single sensitivity, and vice versa. Figure 8d shows that the set of Pareto optimization solutions of ω was the most sensitive for Nu, while the set of Pareto optimization solutions of D was the most sensitive for ρ .
In conclusion, the overall performance of Op-I was that the objective function of Nu had the highest sensitivity for Op-I, ω had the highest correlation or sensitivity for Nu, and D had the highest correlation or sensitivity for ρ . Figure 9a shows the Pareto fronts of Op-II. The Pareto fronts distributed at a curved surface, approximately. The values of freq1 had a great impact on the other two parameters. At the same time, the optimal solutions of the freq1 were all concentrated between 269 Hz and 271 Hz. Therefore, freq1 was the highest for the overall sensitivity of Op-II. Tables 5 and 6 show that the set of Pareto solution of XLS was mainly concentrated in the region with higher Nu. The optimal solutions with Nu between 140 and 172 were 95.3% of the total set of solutions of Nu. The optimal solutions with ρ between 3 and 9  Figure 9a shows the Pareto fronts of Op-II. The Pareto fronts distributed at a curved surface, approximately. The values of freq1 had a great impact on the other two parameters. At the same time, the optimal solutions of the freq1 were all concentrated between 269 Hz and 271 Hz. Therefore, freq1 was the highest for the overall sensitivity of Op-II. Tables 5 and 6 show that the set of Pareto solution of XLS was mainly concentrated in the region with higher Nu. The optimal solutions with Nu between 140 and 172 were 95.3% of the total set of solutions of Nu. The optimal solutions with ρ between 3 and 9 were 56.9% of the total set of solutions of ρ. Therefore, the overall sensitivity of Nu was higher than ρ. Thus, Op-II was more sensitive to freq1 than to Nu and least sensitive to ρ.    It can be found from Figure 9b,c that the optimal solutions of D concentrated between 2.3 mm and 2.9 mm, were generally distributed at a curved surface, and were relatively uniform. As shown in Figure 9d,e, the optimal solution of ω was between 52 • and 66 • . Figure 9f reveals that D and ω had negative correlation with freq1. Nu was more sensitive to ω than to D, while ρ was more sensitive to D than to ω. However, freq1 was not sensitive to D and ω, which means that the change in D and ω has little effect on freq1. freq1 may be most sensitive to H (height of cooling channel), while H herein remains unchanged.

The Representative Function for Optimization Results of Op-I and Op-II
As the set of Pareto optimization solutions obtained above were a series of discrete points, and the value in practical application was not included in these discrete points, it was necessary to provide a continuous set of optimization solution as a reference to the actual situation. Nu = a · e bρ + c (33) Figure 10a,b show that the distribution form of the Pareto set of solutions for Nu and ρ in the two optimization problems presented strong convexity. Therefore, the first-order exponential function was used to fit them. The function equation is shown in Equation (33), and the values of each coefficient in the formula are shown in Table 7.
As the set of Pareto optimization solutions obtained above were a series of discre points, and the value in practical application was not included in these discrete points, was necessary to provide a continuous set of optimization solution as a reference to t actual situation. ρ in the two optimization problems presented strong convexity. Therefore, the first-o der exponential function was used to fit them. The function equation is shown in Equatio (33), and the values of each coefficient in the formula are shown in Table 7.   Table 8 shows the comparison between the initial model and the optimized mod In order to reveal the performance improvement of the optimized structure compar with the original structure, a representative Pareto solution set of Op-I or Op-II was s lected. As can be found from Table 8, for Op-I, when the Nu reached 150, it was abo 24.1% higher than that of the initial model. At the same time, the relative density d creased by 31%. For Op-II, when the Nu reached 160, it was 28.8% higher than that of t initial model, and the relative density was basically the same as the initial model. Ther fore, the optimization results were very effective.   Table 8 shows the comparison between the initial model and the optimized model. In order to reveal the performance improvement of the optimized structure compared with the original structure, a representative Pareto solution set of Op-I or Op-II was selected. As can be found from Table 8, for Op-I, when the Nu reached 150, it was about 24.1% higher than that of the initial model. At the same time, the relative density decreased by 31%. For Op-II, when the Nu reached 160, it was 28.8% higher than that of the initial model, and the relative density was basically the same as the initial model. Therefore, the optimization results were very effective.

Validation of Mechanical and Heat Transfer Performances
As shown in Table 8, a group of optimal structural parameters of the XLS channel to develop a relevant XLS model were obtained, which had relatively excellent mechanical and heat transfer performances by the mathematical optimization model. Two validation methods were employed to verify the accuracy of this mathematical optimization model.

Validation of Bearing Load Performance
Although function of bearing load performance was used as a constraint, its value can be obtained by corresponding experiment to verify the accuracy of the mathematical optimization model. Equivalent elastic modulus (E) was used to evaluate bearing load capacity of XLS cores. A relevant experiment was developed to obtain E.
Firstly, both the initial model and the optimal model were built. Because XLS has the property of periodicity, only a single XLS core was made, as shown in Figure 11a, and two side walls of LS core were removed. Secondly, according to the GB/T 1041-2008 test standard, the unidirectional compression test of the structure was carried out by a ctm2000 universal material testing machine. The pre-cyclic loading method was adopted in the experiment, and the loading rate was 0.5 mm/min. The force displacement curve of the XLS core was obtained (Figure 11b). It can be found that the black curve was the force displacement curve of the initial model of the XLS core, and the maximum force was 8543.267 N, which was converted into the elastic modulus of 680.196 MPa. The theoretical calculation value of equivalent elastic modulus was 830.265 MPa. The error between the maximum valve of the XLS core and the theoretical calculation value was 10.3%. Therefore, E can accurately evaluate the bearing load capacity of the XLS core. Moreover, it can be observed from Figure 11b that the maximum bearing capacities of the optimization result of Op-I and Op-II were 11,089.41 N and 12,404.22 N, respectively. Compared with the initial model, the maximum bearing load capacity of the optimized structure was significantly improved. From the results of Figure 10c, the maximum bearing load capacity for the optimal structure of Op-I was increased by 29.8% compared to that of the initial model. The maximum capacity of bearing load for the optimal structure of Op-II was increased by 45.2% compared to that of the initial model. Therefore, the mathematical optimization model was accurate for bearing load performance.

Validation of Heat Transfer Performance
Nusselt number (Nu) was used to evaluate heat transfer performance of the XLS channel. Relevant numerical simulation was developed to obtain values of Nu for both the initial model and the optimization models to verify accuracy of the mathematical optimization model as follows.
As shown in Figure 2b, the XLS model was established for numerical simulation. The

Validation of Heat Transfer Performance
Nusselt number (Nu) was used to evaluate heat transfer performance of the XLS channel. Relevant numerical simulation was developed to obtain values of Nu for both the initial model and the optimization models to verify accuracy of the mathematical optimization model as follows.
As shown in Figure 2, the XLS model was established for numerical simulation. The grid of XLS was finished by ICEM CFD. The commercial software CFX18.0 was used for simulation, and a standard k-ε model was used for turbulent model. The pressure based on simple algorithm was adopted, and the second-order upwind scheme was selected.
The boundary conditions were as follows. The inlet was a mass flow inlet of 7.67 × 10 −3 Kg/s, the inlet Re was 10,000, the inlet air temperature was set at 300 K, the turbulent degree was 5%, and the heat flux q of the heated wall was constant at 3000 W/m 2 . The data processing methods involved were as follows: The definition of Re: The definition of average Nu: where d h is the hydraulic diameter of the lattice channel, d h = 50.33 mm; q is target heat flux; T W is target surface temperature; T air is inlet airflow temperature; H is height of channel; λ is air thermal conductivity; µ is aerodynamic viscosity, and . m is mass flow rate. In order to verify the reliability of the model, the experimental data of XLS from Kim et al. [38] was selected. The Nu values under different turbulent models are shown in Figure 12. It can be seen that the values of the standard k-ε turbulent model were very consistent with the experimental values, with an error of 9.6%. Therefore, the standard k-ε was used in the subsequent numerical simulation.
where dh is the hydraulic diameter of the lattice channel, dh = 50.33 mm; q is target heat flux; TW is target surface temperature; Tair is inlet airflow temperature; H is height of channel; λ is air thermal conductivity; μ is aerodynamic viscosity, and m  is mass flow rate.
In order to verify the reliability of the model, the experimental data of XLS from Kim et al. [38] was selected. The Nu values under different turbulent models are shown in Figure 12. It can be seen that the values of the standard k-ε turbulent model were very consistent with the experimental values, with an error of 9.6%. Therefore, the standard k-ε was used in the subsequent numerical simulation. Simulation results for fluid flow and heat transfer of the XLS channel were obtained. Figure 13a-c reveal that the XLS had a great disturbing effect on the air flow in the channel, and the flow velocity in the rear along the flow direction was significantly reduced. For the initial model, because the diameter of the member was smaller than that of the optimal structure, the low velocity area formed behind the member was small. For the Op-I model and the Op-II model, due to the increase of member diameter, airflow disturbance in the channel was all strengthened. Figure 13d-f show that the distribution of Nu of the optimal model was more uniform than that of the initial model, and Figure 13g shows that the Nu of the initial model obtained by the numerical calculation was close to the calculated value obtained by the optimized fin method. Moreover, pressure drops under the condition of Re = 10,000 were 302.179 Pa, 294.323 Pa, and 302.558 Pa for the initial  Figure 13a-c reveal that the XLS had a great disturbing effect on the air flow in the channel, and the flow velocity in the rear along the flow direction was significantly reduced. For the initial model, because the diameter of the member was smaller than that of the optimal structure, the low velocity area formed behind the member was small. For the Op-I model and the Op-II model, due to the increase of member diameter, airflow disturbance in the channel was all strengthened. Figure 13d-f show that the distribution of Nu of the optimal model was more uniform than that of the initial model, and Figure 13g shows that the Nu of the initial model obtained by the numerical calculation was close to the calculated value obtained by the optimized fin method. Moreover, pressure drops under the condition of Re = 10,000 were 302.179 Pa, 294.323 Pa, and 302.558 Pa for the initial model, the model of Op-I, and the model of Op-II, respectively. The pressure drops of these three models were similar, and acceptable for the LS channel which had only five simple LS cores. Therefore, pressure drop was not considered as a main parameter in this study.
Therefore, the mathematical optimization model was accurate and reasonable for the heat transfer performance.

Conclusions
1. The present study established a mathematical optimization model based on LS used in thin walls of turbine blades, which was simplified as an LS channel using an integral optimization design method. Two selected optimization problems (Op-I and Op-II) were proposed to investigate optimum structures of LS, and the optimal results calculated by NSGA-II algorithm were proved to be accurate and reasonable.
2. For Op-I (Nu and ρ were the objectives), the overall sensitivity of D and ω for Nu was higher than ρ . The single sensitivity of D for ρ was higher than Nu; the single sensitivity of ω for Nu was higher than ρ . Compared to the initial LS, Nu increased by 24.1% and ρ decreased by 31%, respectively. The optimal structural parameters were 2.3 mm ≤ D ≤ 2.9 mm, and 52° ≤ ω ≤ 66°. 3. For Op-II (Nu, freq1, and ρ were the objectives), the overall sensitivity of D and ω for the objectives was ranked as freq1, Nu, and ρ . The single sensitivity of D for the objectives was ranked as freq1, ρ , and Nu; the single sensitivity of ω for the objectives was ranked as Nu, ρ , and freq1. Compared to the initial LS, Nu increased by 28.8%, while freq1 and ρ were almost unchanged. The optimum structural param-

1.
The present study established a mathematical optimization model based on LS used in thin walls of turbine blades, which was simplified as an LS channel using an integral optimization design method. Two selected optimization problems (Op-I and Op-II) were proposed to investigate optimum structures of LS, and the optimal results calculated by NSGA-II algorithm were proved to be accurate and reasonable.

2.
For Op-I ( Nu and ρ were the objectives), the overall sensitivity of D and ω for Nu was higher than ρ. The single sensitivity of D for ρ was higher than Nu; the single sensitivity of ω for Nu was higher than ρ. Compared to the initial LS, Nu increased by 24.1% and ρ decreased by 31%, respectively. The optimal structural parameters were 2.3 mm ≤ D ≤ 2.9 mm, and 52 • ≤ ω ≤ 66 • .

3.
For Op-II (Nu, freq1, and ρ were the objectives), the overall sensitivity of D and ω for the objectives was ranked as freq1, Nu, and ρ. The single sensitivity of D for the objectives was ranked as freq1, ρ, and Nu; the single sensitivity of ω for the objectives was ranked as Nu, ρ and freq1. Compared to the initial LS, Nu increased by 28.8%, while freq1 and ρ were almost unchanged. The optimum structural parameters were 2.0 mm ≤ D ≤ 3.5 mm, and 56 • ≤ ω ≤ 70 • . 4.
The fitting function of Pareto fronts of Op-I and Op-II was obtained, which may provide guidance for structural parameters design for LS channels used in turbine blades.