Next Article in Journal / Special Issue
Influence of Tip Clearance on Flow Characteristics of Axial Compressor
Previous Article in Journal
Heat Transfer Improvement in MHD Natural Convection Flow of Graphite Oxide/Carbon Nanotubes-Methanol Based Casson Nanofluids Past a Horizontal Circular Cylinder
Previous Article in Special Issue
Experimental and Comparative RANS/URANS Investigations on the Effect of Radius of Volute Tongue on the Aerodynamics and Aeroacoustics of a Sirocco Fan
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Large Eddy Simulation of Film Cooling with Triple Holes: Injectant Behavior and Adiabatic Film-Cooling Effectiveness

School of Mechanical Engineering, Kookmin University, 77 Jeongneung-ro, Seongbuk-gu, Seoul 02707, Korea
*
Author to whom correspondence should be addressed.
Processes 2020, 8(11), 1443; https://doi.org/10.3390/pr8111443
Submission received: 19 October 2020 / Revised: 7 November 2020 / Accepted: 9 November 2020 / Published: 11 November 2020

Abstract

:
This study investigated the effect of adding two sister holes placed downstream the main hole on film cooling by employing large eddy simulation. Here, film-cooling flow fields from a triple-hole system inclined by 35° to a flat plate at blowing ratios of M = 0.5 and unity were simulated. Each sister hole supplies a cooling fluid at a flow rate that is a quarter of that for the main hole. The simulations were conducted using the Smagorinsky–Lilly model as the subgrid-scale model, and the results were compared with those for a single-hole system for the same amount of total cooling air and same cross-sectional area of the holes. Relative to the single-hole system, the spanwise-averaged film-cooling effectiveness in the triple-hole configuration at M = 1.0 increased by as much as 345%. The subsequent proper orthogonal decomposition analysis showed that the kinetic energy of a counter-rotating vortex pair in the triple-hole system dropped by 30–40% relative to that of the single-hole system. This indicates that the additional sister holes promoted the suppression of the mixing of the coolant jet with the mainstream flow, thus keeping the temperature of the latter low. Cross-sectional views of the root-mean-square temperature contours were also analyzed; with the results confirming that the effect of the sister holes on the jet trajectory greatly contributes in promoting film-cooling effectiveness as compared to the effect of the reduced mixing.

1. Introduction

The principle of idealized Brayton cycle states that gas turbine efficiency can be improved by increasing the inlet temperature of the turbine [1]. Here, to avoid the excessive thermal stress build-up on the turbine blades, the surface temperature of the blades must remain below a tolerable limit. Various cooling methods can be employed, among which the film-cooling technique has been widely used. In this technique, the cooling air which is bled from the compressor is injected through holes on the turbine blade surface—this reduces the surface temperature and protects the surface from the main flow. For the holes, a cylindrical shape is often appropriate because of its simplicity and ease of manufacture.
However, substantial investigations of the physics of film cooling in such a cylindrical hole system have revealed that a single cylindrical hole can be relatively vulnerable to the generation of the jet liftoff [2] and that a specific hole arrangement could produce film cooling that provides better thermal protection than single cylindrical holes do. Practically, specific hole arrangements make the intended vortex interactions between film cooling jets to increase film-cooling effectiveness.
Film cooling in turbine engines can be effectively studied at a reasonable cost using computational fluid dynamics (CFD). High turbulence is generated in the film-cooling flow field, after which a model is constructed for the turbulence. Another option is the well-known large eddy simulation (LES), which can predict the mixing between the main flow and coolant jet better than Reynolds-averaged Navier–Stokes Simulation (RANS) [3,4,5], although such requires much higher computational costs. Comparatively speaking, RANS cannot accurately predict the complex flow structures induced in film-cooling flow fields as all of the turbulent fluctuations are ensemble-averaged, whereas LES resolves large-scale eddies directly in the turbulent flow, resulting in more accurate predictions of the complex flow [6,7].
Heidmann et al. [8] showed in a numerical investigation that relative to a one-hole system, three holes (consisting of the main hole and two side holes) can increase film-cooling effectiveness, as well as decrease the heat transfer coefficients at a blowing ratio of unity, using RANS. They also reported that the fabrication of cylindrical triple holes is economically more viable than that of other shapes. Meanwhile, Javadi et al. [9] constructed square sister holes to control vortex interactions between the jet and crossflow at M = 0.5 by employing RANS. They found that two sister holes with a square cross-section reduce the mixing between the coolant and crossflow, thus improving film-cooling effectiveness. Ely et al. [10] also used RANS to investigate an increase in film-cooling effectiveness with two cylindrical sister holes. Using blowing ratios of 0.2, 0.5, 1.0, and 1.5, they demonstrated that definite improvements in film-cooling effectiveness occur at high blowing ratios.
Furthermore, Wu et al. [11] conducted both experimental and numerical studies of film cooling given three sister-hole configurations and confirmed an improvement in the performance as above. Here, they compared the experimental and numerical data only in terms of spatially averaged film-cooling effectiveness, which showed a difference of more than 30% under blowing ratios of unity and 1.5.
Farhadi-Azar et al. [12] introduced LES to study the effects of the triple jets on film cooling. On the basis of their demonstration, triple jets from the main hole along with two small sister holes increased film-cooling effectiveness because of weaker counter-rotating vortex pairs in the jet from the main hole, which led to less mixing between the coolant jets and main flow. However, several limitations could be found in their study. For instance, they employed an injection angle (α) of 90° on their hole, which is very uncommon for a gas turbine design because of the strong jet liftoff. Aside from this, the plenum was not included in the simulation, and a flow with a uniform velocity was applied to the hole inlet, which is unrealistic in gas turbine scenarios.
Considering the above, this paper aims to investigate the effect of a triple-hole system on film-cooling performance on a flat plate at blowing ratios of 0.5 and 1.0. The simulation was conducted using LES under a more realistic flow condition and at an injection angle (α) of 35° for the hole. The plenum is included, and the velocity in the hole is nonuniform with both the high-momentum region generated by the jetting effect and low-momentum region.
The LES results are used to understand the effect of the flow generated in the triple hole on the trajectory of the film-cooling jet and mixing with the mainstream flow. Vortex structures that affect the mixing are identified through time-averaged velocity fields and instantaneous flow fields, whereas the contribution of these vortices is determined through a proper orthogonal decomposition (POD) analysis. The film-cooling performance, degree of mixing with the main flow, and trajectory of the film-cooling jet are also observed through the time-averaged temperature fields. On the basis of the analysis results, the effect and mechanism of the triple hole on the film-cooling performance are identified.

2. Geometry and Boundary Conditions

The geometry of the single-hole system in the test case was taken from Seo et al. [13], whose experimental apparatus consisted of a row of five cooling holes. Herein, a single-hole configuration was adopted to simplify the computational domain, and a periodic boundary condition was applied on the mesh sides to reduce the computational cost. The computational domains of the single- and triple-hole systems are delineated by the orange dashed lines in Figure 1, whereas a schematic of the cooling holes in each case is shown in Figure 2. In the single-hole geometry (Figure 2a), the hole diameter D was 25 mm. The triple holes comprised two sister holes which are placed slightly downstream of the primary hole (Figure 2b) as described in Ely et al. [10]. In the triple-hole geometry, the hole diameters were calculated to match their total cross-sectional areas to the cross-sectional area of the single-hole. Thus, the primary hole diameter D′ was 20.4124 mm, and the diameter of each sister hole was D′/2 = 10.2062 mm. The hole length-to-diameter ratio (L/D), injection angle (α), and pitch-to-hole diameter (P/D) were set to 1.6, 35°, and 3.0, respectively. In addition, the holes were cylindrical, with no compound angles.
The boundary conditions in each computational domain are specified in Table 1. The turbulence intensity of the mainstream flow was 0.2% at the main inlet (as reported experimentally [13]), and the main flow velocity was 10 m/s. The temperatures of the main flow and coolant at the inlets were 313 and 293 K, respectively. Moreover, the boundary layer thickness at the cooling hole center is about D.
As is well-known, injecting the film-cooling jet induces a counter-rotating vortex pair (CRVP) in the system (see Figure 3a). This vortex pair is sourced from the vortices generated in the boundary layer of the coolant tube wall and around the leading edge of the hole [12,14,15]. Interactions between the jet and main flow cause the injected coolant to split into two directions, whereas the vortices generated in the coolant are stretched and turned. As the CRVP vortices draw closer together, their mutual induction increases, promoting coolant jet liftoff. The consequent strong entrainment of the hot main flow under the jet increases the adiabatic temperature on the test plate and decreases film-cooling effectiveness. Moreover, when two sister holes are placed slightly downstream of the primary hole (Figure 2b), their jets generate an anti-CRVP, as illustrated in Figure 3b. Therefore, film-cooling effectiveness is increased by establishing vortex interactions between the jets and by reducing the intensity of the main CRVP [9,10,11,12].

3. Validation of Numerical Methods

The CFD calculations were performed in ANSYS Fluent v.19.1 [16], and the meshes were generated in Pointwise v.18.1 [17]. The numerical simulations were executed in the LES using the Smagorinsky–Lilly model as the subgrid-scale model. The time step was set to 6.25 × 10−6 s. Assuming the mainstream convected the hole diameter after 400 time steps, the simulations were carried out with the computational time step Δ t = 0.0025 D / U [18,19]. Once the system had reached a statistical steady state, the statistics of the solution was collected in multiples of the period. In each time step, approximately 10 subiterations were executed to obtain well-resolved data [20]. The LES runtime using 18 cores of an Intel Xeon Gold 6148 processor was 1–2 months in each case.
The fluid was assumed as Newtonian and incompressible and had temperature-dependent variable properties. The compressibility effect was negligible because the main flow velocity was 10 m/s (Mach 0.029) and the jet injection velocities were 5 and 10 m/s at M = 0.5 and M = 1.0, respectively; that is, both velocities were far below Mach 0.3 [21]. The governing equations of the CFD simulations are the continuity and momentum equations. The filtered Navier–Stokes equations in the LES approach are given by [22]
ρ t + x i ( ρ u ¯ i ) = 0
and
( ρ u ¯ ) t + x j ( ρ u ¯ i u ¯ j ) = x j [ μ ( u i ¯ x j + u j ¯ x i )   2 3 μ u l ¯ x l δ i j ] d p ¯ d x + τ i j x j
where τij, the subgrid-scale turbulent stress, needs modeling using the Boussinesq hypothesis like RANS models as
τ i j 1 3 τ k k δ i j = μ t ( u i ¯ x j + u j ¯ x i )
where μt is the turbulent viscosity on the subgrid scale. In the Smagorinsky–Lilly model, the turbulent viscosity is modeled as [18]
μ t = ρ L s 2 1 2 ( u i ¯ x j + u j ¯ x i ) ( u i ¯ x j + u j ¯ x i )
Figure 4 and Figure 5 show the meshes of the single- and triple-hole systems in the xy and yz planes, respectively. The computational domains of the systems were composed of 2.04 (single-hole) and 2.69 (triple-hole) million hexahedron cells. Figure 6 shows a close-up of the mesh near the hole, where P2 and P1 are the static pressures at the inlet and outlet of the hole, respectively, as mentioned earlier.
Figure 7 and Figure 8 show the results of the grid sensitivity tests in each system at a blowing ratio of 0.5, with the grid descriptions given in Table 2 and Table 3. The results were obtained using the Smagorinsky–Lilly model in LES, and the tests were performed on five different grids. In the single-hole system, the effectiveness value along the centerline (z = 0, y = 0) on the third grid almost matched those of the finer fourth and fifth grids.
Therefore, the adiabatic film-cooling effectiveness in the single-hole system was evaluated on the 2.04 million-cell grid. In the triple-hole system, the third grid with 2.69 million cells again showed almost the same centerline effectiveness as the fourth and fifth grids. Therefore, the 2.69 million-cell grid was selected for modeling the film-cooling effectiveness in the triple-hole system.

4. Results and Discussion

4.1. Centerline and Spanwise-Averaged Film Cooling Effectiveness

Figure 9 shows the centerline and spanwise-averaged effectiveness values in the single- and triple-hole systems at M = 0.5 and 1.0. The centerline effectiveness decreased as the main flow traveled downstream. This degradation is attributable to the turbulence generation and rapid mixing between the main flow and jet, which increases the test plate temperature. For the same total cross-sectional hole area and same amount of total cooling air, the triple holes improved the centerline and spanwise-averaged effectiveness values, respectively, by approximately 36% and 45% over the single-hole configuration at M = 0.5 and approximately 250% and 345% at M = 1.0. In Figure 9a–d, the LES-predicted centerline film-cooling effectiveness in the single-hole system showed a good match to the experimental data by Seo et al. [13], whereas the LES model under-predicted the triple-hole, spanwise-averaged effectiveness experimental data by Cao et al. [22] because the data were measured for L/D = 10. Thus, compared to lower L/D cases, higher L/D values of the hole increase the adiabatic film-cooling effectiveness under the same operating condition as the coolant velocity at the hole exit becomes more uniform [13].
The ratios of the spanwise-averaged effectiveness values in the triple-hole system to those in the single-hole system as functions of x/D for M = 0.5 and 1.0 are plotted in Figure 10. Note that increasing the blowing ratio improves the ratio of the spanwise-averaged effectiveness. This indicates that under high blowing ratios, configuring the main hole and two sister holes is more effective for film cooling than using an ordinary single hole mainly because the jet trajectory in the single-hole system is mostly governed by the jet liftoff at a high blowing ratio of M = 1.0. The result is very low adiabatic film-cooling effectiveness, as indicated in Figure 9c,d, and very high ratios of the spanwise-averaged effectiveness values, as shown in Figure 10. The maximum value of the ratio at M = 1.0 is observed around at x/D = 4 since the minimum value of the spanwise-averaged effectiveness in the single-hole system is shown around at x/D = 4 due to the jet lift off, as shown in Figure 9d.

4.2. Contours of Film Cooling Effectiveness and Dimensionless Temperature

Figure 11 presents top views of the time-averaged film-cooling effectiveness contours on the test plate in each hole system. The additional sister holes largely increased the amount of coolant in contact with the test plate, which boosts film-cooling effectiveness. Moreover, the spread of the coolant in the spanwise direction in the triple-hole system was much larger than that in the single-hole configuration as the main CRVP weakens. Furthermore, lower film-cooling effectiveness was predicted at M of unity than at M = 0.5 because of strong jet liftoff.
The cross-sectional views of the dimensionless mean temperature contours at x/D = 2, 5, 10, 15, and 20 in the two systems are illustrated in Figure 12. Initially, the dimensionless temperature near the test plate was higher in the triple-hole system than in the single-hole configuration, which increased film-cooling effectiveness on the test plate. As the main flow went downstream, the dimensionless temperatures near the test plate in the triple-hole system became similar to those in the single-hole system, leading to almost identical film-cooling effectiveness in the range of x/D > 15 (Figure 8 and Figure 9) because of intensive mixing between the mainstream flow and coolant jet.
Figure 13 shows the central cross-sectional contours of the dimensionless mean temperatures in the z = 0 plane for both hole systems. The entrainment of the main flow under the cooling jet in the triple-hole system was noticeably weaker than that in the single-hole system. Many researchers have argued that the jet lifted off because of the high momentum of the coolant and reattached downstream to the plate at a high blowing ratio [2,22,23]. However, as illustrated in Figure 13b, the dispersion of the coolant jet increases without reattachment with the test plate as the flow goes downstream after the jet liftoff, followed by the entrainment of the main flow under the jet. In addition, the dimensionless temperature near the test plate in the near-hole region in the triple-hole system was much higher than that in the single-hole configuration as vortex interactions between the jets were established, thereby reducing the intensity of the main CRVP at M = 0.5 and 1.0 and increasing film-cooling effectiveness on the test plate.
Figure 14 illustrates the streamlines of cooling jets to a cross-sectional plane at x/D = 5 and film-cooling effectiveness on the test plate in the single- and triple-hole systems at M = 0.5 and unity. Unlike in a single hole, more streamlines of the coolant head to the region near the test plate in the x/D = 5 plane appeared in the triple-hole system, with these streamlines more dispersed, which resulted in higher adiabatic film-cooling effectiveness on the plate given that most streamlines from the sister holes have converged to the centerline. At M = 1.0, the strong jet liftoff caused the streamlines to cover smaller regions laterally in the x/D = 5 plane even in the triple-hole system.

4.3. Contours of x-Vorticity and Mean Velocity Vectors

Figure 15 shows the cross-sectional contours of the mean x-vorticity at the x/D = 2 and x/D = 5 planes in the single- and triple-hole systems at M = 0.5 and unity. The mean velocity vectors are superimposed on the contours. As shown in Figure 15b,d of the triple-hole system, distinct anti-CRVPs were found from the sister holes around the x/D = ±0.5 and y/D = 0.2 regions in the x/D = 2 plane. The anti-CRVPs weakened the main CRVP and reduced the main jet liftoff, thereby reducing the entrainment of the hot mainstream gas under the coolant and increasing the dimensionless mean temperature and film-cooling effectiveness, as illustrated in Figure 11 and Figure 13. However, the intensity of the anti-CRVPs became very small in the x/D = 5 plane because the anti-CRVP from the sister holes was not as strong as the main CRVP. In the triple-hole system, the upward velocity vectors around the line of z/D = 0 were less concentrated than those in the single-hole system. Thus, coolant liftoff was weakened by the generated anti-CRVPs. At M = 1.0, the positions of the core for the CRVPs were higher than those at M = 0.5 because of the strong jet liftoff generated.
Figure 16 shows the cross-sectional velocity vectors of POD Mode 1 at x/D = 2 in the two-hole systems at M = 0.5 and unity. The POD was introduced by Dr. John Lumley in 1967 in order to decompose a random velocity vector field of turbulence into a set of portions of the total kinetic energy [24,25]. On the basis of the figure, CRVP was most distinct in the velocity vector field at x/D = 2. Figure 17 illustrates the portion of the total kinetic energy at x/D = 2 for Modes 1, 2, 3, and 4 in the single- and triple-hole systems at M = 0.5 and M = 1.0, which showed Mode 1 as the most dominant. At M = 0.5, POD Mode 1 contributed to about 60% and 36% of the total kinetic energy in the single- and triple-hole systems, respectively, which means that CRVP was dominant for both, although relative to the former, the portion of the kinetic energy in the latter system was largely decreased by about 40% because of a weakened main CRVP by anti-CRVPs from the two sister holes. Similarly, at M = 1.0, POD Mode 1 contributed to about 71% and 49% of the total kinetic energy in the single- and triple-hole systems, respectively, indicating that the portion of the kinetic energy was decreased by nearly 30%. Indeed, it could be inferred that compared to that in the single-hole system, the kinetic energy of the CRVP in the triple-hole system was reduced by about 30–40%.

4.4. Root-Mean-Square Temperature Contours

Figure 18 illustrates the cross-sectional views of the root-mean-square (RMS) temperature contours at x/D = 2, 5, 10, 15, and 20 for the mixing region between the mainstream flow and coolant in the single- and triple-hole systems. The mixing area was more expanded in the triple-hole system in the spanwise direction, especially in the near-hole region as x/D = 2, because the weakened main CRVP allowed wider spreading and closer contact of the coolant jet with the test plate. Although no significant difference in the mixing intensity was observed in both hole systems, the Trms in the x/D = 2 plane reached maximum values of 6.25, 6.24, 6.08, and 5.19, as shown in Figure 18a–d, respectively, which indicates a similar but slightly decreased mixing intensity in the triple-hole system relative to the single-hole system. This reduced mixing intensity is beneficial to film cooling because increased mixing between the coolant and mainstream flow consequently decreases film-cooling effectiveness [22]. Nevertheless, it could be stated that in the triple-hole system, the effect of the sister holes on the jet trajectory contributes greatly to improve film-cooling effectiveness as compared to the effect of less mixing between the main flow and coolant.

4.5. Q-Criterion Iso-Surfaces

In the Q-criterion, Q is defined as the difference between the magnitudes of the rotation rate and strain rate, representing the dominance of the former over the latter [26,27,28]. The instantaneous contours of a Q-criterion with iso-surfaces of level 80,000 at M = 0.5 and unity are displayed in Figure 19. Q-criterion contours in the single-hole system clearly showed various vortex structures, such as the CRVP, hairpin and jet shear layer vortices, and horseshoe vortex, whereas those in the triple-hole configuration had smaller hairpin and jet shear layer vortices and were not very distinguishable. These smaller and indistinct hairpin vortices represent the weakened CRVP as the hairpin vortex is closely related to the generation of the CRVP [22]. A weakened CRVP positively affects the cooling performance on the plate because it increases film-cooling effectiveness on the test plate.

4.6. Cross-Sectional Contours of the Mean Velocity Magnitude in the Hole

Figure 20 shows the cross-sectional contours of the mean velocity magnitude in the hole measured around the hole exit (y/D = −0.2). As mentioned earlier, the ratio of the delivery tube length to the hole diameter in the single-hole system was 1.6, whereas it was 1.8 and 3.6 for the main hole and sister holes in the triple-hole system, respectively, because the triple holes had the same total cross-sectional hole area as the single hole. Therefore, the mean velocity magnitude in the main hole is not uniform but has high- and low-momentum regions, whereas that in the sister holes shows more developed coolant flow. The cooling jet having a uniform velocity profile is advantageous to film-cooling performance because of less jet liftoff. Thus, it could be stated that relative to the single-hole system, the sister holes improve film-cooling effectiveness by the generation of the anti-CRVP, as well as higher L/D ratios.

5. Conclusions

Film-cooling simulations of the single- and triple-hole systems on a flat plate were carried out in this paper at blowing ratios of M = 0.5 and unity by employing LES. On the basis of these simulations, when two sister holes which were placed slightly downstream the main hole, the coolant jets from the former generated an anti-CRVP, which increased film-cooling effectiveness. Thus, for the same total cross-sectional hole area and same amount of total cooling air, the triple-hole system increased the centerline and spanwise-averaged effectiveness values by approximately 36% and 45% at M = 0.5 and by approximately 250% and 345% at M = 1.0 as compared to that done by the single-hole configuration. Moreover, the ratios of the spanwise-averaged effectiveness values in the triple-hole system indicate that it is the more effective system for film cooling under a high blowing ratio.
Meanwhile, the time-averaged film-cooling effectiveness contours on the test plate and the cross-sectional views of the dimensionless mean temperature contours in the two hole systems showed that additional sister holes can increase the amount of coolant contacting with the test plate and boost film-cooling effectiveness. Specifically, the mean dimensionless temperature contours in the z/D = 0 plane showed that in the single-hole system at M of unity, as the cooling jet goes downstream, the jet gets dispersed without reattachment with the test plate, and the entrainment of the main flow under the jet is generated.
Furthermore, the cross-sectional mean x-vorticity contours in the triple-hole system showed that the generated anti-CRVPs weaken the main CRVP. The results of the POD analysis confirmed that the kinetic energy of the CRVP in the triple-hole system was reduced by about 30–40% relative to that in the single-hole system. Streamlines of the cooling jets to the cross-sectional plane at x/D = 5 also showed that most streamlines from the sister holes converge to the centerline. As with the Q-criterion contours in the triple-hole system, the hairpin and jet shear layer vortices were smaller and indistinguishable, unlike in the single-hole system.
Lastly, on the basis of the cross-sectional views of the RMS temperature contours, the effect of the sister holes on the cooling jet trajectory contributes significantly to the increase in film-cooling effectiveness, as compared to the effect of less mixing. Accordingly, the sister holes improve film-cooling effectiveness using higher L/D ratios.

Author Contributions

S.I.B. carried out the simulations, analyzed the CFD data and wrote the paper. J.A. supervised the research, analyzed the CFD data and wrote the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declared no conflict of interest.

Nomenclature

CsSmagorinsky constant
Dhole diameter for a single hole [mm]
D’diameter of primary hole for triple holes [mm]
E i j rate of deformation
Ldelivery tube length [mm]
Lsmixing length of sub-grid scales = min ( κ d ,   C s Δ )
Mblowing ratio = ( ρ C U C ) / ( ρ G U G )
Ppitch between holes [mm]
Ttemperature [K]
ttime [s]
Uflow velocity [m/s]
Vmimain flow velocity at the main inlet [m/s]
Vpicoolant velocity at the plenum inlet [m/s]
xstreamwise coordinate
ywall-normal coordinate
zspanwise coordinate
Greek symbols
κvon Karman’s universal constant = 0.41
η adiabatic film cooling effectiveness =   ( T G T a w ) T G   T C
η C centerline film cooling effectiveness
η m spanwise-averaged film cooling effectiveness
ρ density [kg/m3]
τijsub-grid scale turbulent stress =   ρ u i u j ¯   ρ u i ¯ u j ¯
μtsub-grid scale turbulent viscosity [kg/(m·s)]
νlocal kinematic viscosity [m2/s]
local grid scale
Θdimensionless temperature =   ( T G T ) T G   T C
Subscripts
awadiabatic wall
ccenterline
Ccoolant
Gmainstream gas
mspanwise averaged
μ local dynamic viscosity [m2/s]

References

  1. Moran, M.; Shapiro, H.; Boettner, D.; Bailey, M. Fundamentals of Engineering Thermodynamics, 8th ed.; Wiley: Hoboken, NJ, USA, 2014; p. 532. [Google Scholar]
  2. Bogard, D. Airfoil Film Cooling. In The Gas Turbine Handbook; National Energy Technology Laboratory: Albany, NY, USA, 2006. [Google Scholar]
  3. Leedom, D.H.; Acharya, S. Large Eddy Simulations of Film Cooling Flow Fields From Cylindrical And Shaped Holes. In Proceedings of the ASME Turbo Expo 2008, Berlin, Germany, 9–13 June 2008; GT2008-51009. pp. 865–877. [Google Scholar]
  4. Acharya, S.; Leedom, D.H. Large Eddy Simulations of Discrete Hole Film Cooling With Plenum Inflow Orientation Effects. J. Heat Transf. 2013, 135, 011010. [Google Scholar] [CrossRef]
  5. Tyagi, M.; Acharya, S. Large Eddy Simulation of Film Cooling Flow From an Inclined Cylindrical Jet. J. Turbomach. 2003, 125, 734–742. [Google Scholar] [CrossRef]
  6. Fujimoto, S. Large Eddy Simulation of Film Cooling Flows Using Octree Hexahedral Meshes. In Proceedings of the ASME Turbo Expo 2012, Copenhagen, Denmark, 11–15 June 2012. GT-2012-7009. [Google Scholar]
  7. Rozati, A.; Tafti, D. Large-eddy simulations of leading edge film cooling: Analysis of flow structures, effectiveness, and heat transfer coefficient. Int. J. Heat Fluid Flow 2008, 29, 1–17. [Google Scholar] [CrossRef]
  8. Heidmann, J.; Ekkad, S. A Novel Anti-Vortex Turbine Film Cooling Hole Concept. In Proceedings of the ASME Turbo Expo 2007, GT2007-27528, Montreal, QC, Canada, 14–17 May 2007. [Google Scholar] [CrossRef]
  9. Javadi, K.; Taeibi-Rahni, M.; Darbandi, M.; Javadi, K. Jet-into-Crossflow Boundary Layer Control: Innovation in Gas Turbine Blade Cooling. AIAA J. 2007, 45, 2910–2925. [Google Scholar] [CrossRef]
  10. Ely, M.J.; Jubran, B.A. A Numerical Evaluation on the Effect of Sister Holes on Film Cooling Effectiveness and the Surrounding Flow Field. Heat Mass Transf. 2009, 45, 1435–1446. [Google Scholar] [CrossRef]
  11. Wu, H.; Cheng, H.; Li, Y.; Rong, C.; Ding, S. Effects of Side Hole Position and Blowing Ratio on Sister Hole Film Cooling Performance in a Flat Plate. Appl. Therm. Eng. 2016, 93, 718–730. [Google Scholar] [CrossRef]
  12. Farhadi-Azar, R.; Ramezanizadeh, M.; Taeibi-Rahni, M.; Salimi, M. Compound Triple Jets Film Cooling Improvements via Velocity and Density Ratios: Large Eddy Simulation. J. Fluids Eng. 2011, 133, 031202. [Google Scholar] [CrossRef]
  13. Seo, H.; Leea, J.; Ligrani, P. The Effect of Injection Hole Length on Film Cooling with Bulk Flow Pulsations. Int. J. Heat Mass Transf. 1998, 41, 3515–3528. [Google Scholar] [CrossRef]
  14. Fric, T.F.; Roshko, A. Vortical structure in the wake of a transverse jet. J. Fluid Mech. 1994, 279, 1–47. [Google Scholar] [CrossRef] [Green Version]
  15. Haven, B.A.; Kurosaka, M. Kidney and anti-kidney vortices in crossflow jets. J. Fluid Mech. 1997, 352, 27–64. [Google Scholar] [CrossRef]
  16. ANSYS Fluent Theory Guide version 19. Available online: https://www.ansys.com/products/fluids/ansys-fluent (accessed on 7 October 2020).
  17. Pointwise version 18. Available online: http://www.pointwise.com/ (accessed on 7 October 2020).
  18. Renze, P.; Schröder, W.; Meinke, M. Large-eddy Simulation of Film Cooling Flows with Variable Density Jets. Flow Turbul. Combust. 2008, 80, 119–132. [Google Scholar] [CrossRef]
  19. Peet, Y.; Lele, S. Near field of film cooling jet issued into a flat plate boundary layer: LES STUDY. In Proceedings of the ASME Turbo Expo 2008, GT2008-50420, Berlin, Germany, 9–13 June 2008. [Google Scholar]
  20. Yao, Y.; Zhang, J.; Yang, Y. Numerical study on film cooling mechanism and characteristics of cylindrical holes with branched jet injections. Propuls. Power Res. 2013, 2, 30–37. [Google Scholar] [CrossRef] [Green Version]
  21. White, F. Fluid Mechanics, 4th ed.; McGraw-Hill: New York, NY, USA, 1999. [Google Scholar]
  22. Cao, N.; Li, X.; Wu, Z.; Luo, X. Effect of film hole geometry and blowing ratio on film cooling performance. Appl. Therm. Eng. 2020, 165, 114578. [Google Scholar] [CrossRef]
  23. Bacci, T.; Picchi, A.; Facchini, B. Flat plate and turbine vane film cooling performance with laid-back fan shaped holes. Int. J. Turbomach. Propuls. Power 2019, 4, 14. [Google Scholar] [CrossRef] [Green Version]
  24. Lumley, J. The structure of inhomogeneous turbulent flows. In Atmospheric Turbulence and Radio Wave Propagation, Proceedings of the International Colloquium, Moscow, Russia, 15–22 June 1965; Publishing House Nauka: Moscow, Russia, 1967; Available online: https://www.amazon.com/Atmospheric-Turbulence-Radio-Wave-Propagation/dp/B007V5IG6W (accessed on 7 October 2020).
  25. Weiss, J. A tutorial on the proper orthogonal decomposition. In Proceedings of the AIAA Aviation Forum, Dallas, TX, USA, 17–21 June 2019. [Google Scholar]
  26. Kolar, V. Vortex identification: New requirements and limitations. Int. J. Heat Fluid Flow 2007, 28, 638–652. [Google Scholar] [CrossRef]
  27. Schroder, A.; Willert, C. Particle Image Velocimetry: New Developments and Recent Applications; Springer: New York, NY, USA, 2008; p. 382. [Google Scholar]
  28. Haydt, S.; Lynch, S. Flowfield of a shaped film cooling hole over a range of compound angles. In Proceedings of the ASME Turbo Expo proceeding, GT2018-75728, Lillestrom, Norway, 11–15 June 2018. [Google Scholar]
Figure 1. Computational domains of the single- and triple-hole configurations (orange dashed lines).
Figure 1. Computational domains of the single- and triple-hole configurations (orange dashed lines).
Processes 08 01443 g001
Figure 2. Schematics of the single hole and triple holes cooling hole configurations.
Figure 2. Schematics of the single hole and triple holes cooling hole configurations.
Processes 08 01443 g002
Figure 3. Cross-sectional schematic of a counter rotating vortex pair generated in the single- and triple-hole configurations.
Figure 3. Cross-sectional schematic of a counter rotating vortex pair generated in the single- and triple-hole configurations.
Processes 08 01443 g003
Figure 4. CFD meshes of the single- and triple-hole systems in the xy plane.
Figure 4. CFD meshes of the single- and triple-hole systems in the xy plane.
Processes 08 01443 g004
Figure 5. CFD meshes in the yz plane.
Figure 5. CFD meshes in the yz plane.
Processes 08 01443 g005
Figure 6. Close-up of the mesh near the hole.
Figure 6. Close-up of the mesh near the hole.
Processes 08 01443 g006
Figure 7. Grid sensitivity test of the large eddy simulation (LES) calculation in the single-hole system.
Figure 7. Grid sensitivity test of the large eddy simulation (LES) calculation in the single-hole system.
Processes 08 01443 g007
Figure 8. Grid sensitivity test of the LES calculation in the triple-hole system.
Figure 8. Grid sensitivity test of the LES calculation in the triple-hole system.
Processes 08 01443 g008
Figure 9. Film cooling effectiveness values of the single- and triple-hole systems.
Figure 9. Film cooling effectiveness values of the single- and triple-hole systems.
Processes 08 01443 g009
Figure 10. Ratios of spanwise-averaged effectiveness in the triple hole system to those in the single-hole system.
Figure 10. Ratios of spanwise-averaged effectiveness in the triple hole system to those in the single-hole system.
Processes 08 01443 g010
Figure 11. Time averaged film cooling effectiveness on the test plate.
Figure 11. Time averaged film cooling effectiveness on the test plate.
Processes 08 01443 g011
Figure 12. Cross-sectional contours of the dimensionless mean temperatures.
Figure 12. Cross-sectional contours of the dimensionless mean temperatures.
Processes 08 01443 g012aProcesses 08 01443 g012b
Figure 13. Central cross-sectional contours of the dimensionless mean temperatures in the z = 0 plane.
Figure 13. Central cross-sectional contours of the dimensionless mean temperatures in the z = 0 plane.
Processes 08 01443 g013
Figure 14. The streamlines of cooling jets to cross-sectional plane at x/D = 5.
Figure 14. The streamlines of cooling jets to cross-sectional plane at x/D = 5.
Processes 08 01443 g014
Figure 15. Cross-sectional contours of the x-vorticity and velocity vectors.
Figure 15. Cross-sectional contours of the x-vorticity and velocity vectors.
Processes 08 01443 g015aProcesses 08 01443 g015b
Figure 16. Cross-sectional velocity vectors at x/D = 2 for POD Mode 1.
Figure 16. Cross-sectional velocity vectors at x/D = 2 for POD Mode 1.
Processes 08 01443 g016
Figure 17. Portions of the total kinetic energy for Mode 1, 2, 3, and 4 in the plane of x/D = 2.
Figure 17. Portions of the total kinetic energy for Mode 1, 2, 3, and 4 in the plane of x/D = 2.
Processes 08 01443 g017
Figure 18. Cross-sectional contours of the Trms.
Figure 18. Cross-sectional contours of the Trms.
Processes 08 01443 g018
Figure 19. Q-criterion contours with iso-surfaces of level 80,000.
Figure 19. Q-criterion contours with iso-surfaces of level 80,000.
Processes 08 01443 g019
Figure 20. Cross-sectional contours of the mean velocity magnitude in the hole at y/D = −0.2.
Figure 20. Cross-sectional contours of the mean velocity magnitude in the hole at y/D = −0.2.
Processes 08 01443 g020
Table 1. Boundary conditions of the computational fluid dynamics (CFD) domain.
Table 1. Boundary conditions of the computational fluid dynamics (CFD) domain.
SurfaceBoundary Condition
Main inletVelocity inlet
Plenum inletVelocity inlet
TopSymmetry
Test plateAdiabatic wall
OutflowPressure outlet
Main sidesPeriodic
Sides of plenumWall
Tube wallWall
Table 2. Specifications of mesh arrangements in the grid sensitivity test of the single-hole system.
Table 2. Specifications of mesh arrangements in the grid sensitivity test of the single-hole system.
GridNumber of Cells in the x DirectionNumber of Cells in the y DirectionNumber of Cells in the z DirectionNumber of Cells in Cross Flow Block (million)Total Cells Number (million)
First32050320.521.14
Second33460480.981.60
Third35280501.422.04
Fourth36494561.942.56
Fifth390110642.783.4
Table 3. Specifications of mesh arrangements in the grid sensitivity test of the triple-hole system.
Table 3. Specifications of mesh arrangements in the grid sensitivity test of the triple-hole system.
GridNumber of Cells in the x DirectionNumber of Cells in the y DirectionNumber of Cells in the z DirectionNumber of Cells in Crossflow Block (million)Total Cells Number (million)
First30046651.011.78
Second33060751.622.39
Third35070801.922.69
Fourth37080852.413.19
Fifth39090902.943.71
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Baek, S.I.; Ahn, J. Large Eddy Simulation of Film Cooling with Triple Holes: Injectant Behavior and Adiabatic Film-Cooling Effectiveness. Processes 2020, 8, 1443. https://doi.org/10.3390/pr8111443

AMA Style

Baek SI, Ahn J. Large Eddy Simulation of Film Cooling with Triple Holes: Injectant Behavior and Adiabatic Film-Cooling Effectiveness. Processes. 2020; 8(11):1443. https://doi.org/10.3390/pr8111443

Chicago/Turabian Style

Baek, Seung Il, and Joon Ahn. 2020. "Large Eddy Simulation of Film Cooling with Triple Holes: Injectant Behavior and Adiabatic Film-Cooling Effectiveness" Processes 8, no. 11: 1443. https://doi.org/10.3390/pr8111443

APA Style

Baek, S. I., & Ahn, J. (2020). Large Eddy Simulation of Film Cooling with Triple Holes: Injectant Behavior and Adiabatic Film-Cooling Effectiveness. Processes, 8(11), 1443. https://doi.org/10.3390/pr8111443

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop