Combined Experimental and CFD Investigation of Flat Plate Film Cooling through Fan Shaped Holes

The present paper reports the results of an experimental and computational investigation of flat plate film cooling jets discharged from three fan-shaped holes. Measurements have been carried out at near unity density ratio in a low-speed wind tunnel, at low inlet turbulence intensity, with blowing ratios (BR) of 1 and 2. Aerodynamic results have shown that the jet stays attached to the flat plate. Thermal measurements have revealed that film cooling effectiveness decreases downstream of the holes, and BR equal to 1 provides the best trade-off between cooling air consumption and thermal protection. Consequently, BR = 1 was selected for assessing the performance of different turbulence models, implemented in STAR-CCM+, according with the steady Reynolds-averaged Navier–Stokes (RANS) approach. Predictions from realizable k-ε (RKE), shear stress transport k-ω (SST KW) and Reynolds stress model (RSM) were compared against measurements of laterally averaged and centerline adiabatic effectiveness, as well as off-the-wall velocity maps and profiles of stress components. RSM provided the most accurate predictions.


Introduction
There is no doubt that film cooling is an integral part of cooling techniques in the gas turbine industry to protect gas turbine components from high turbine inlet temperature.Hole geometry, of all aspects that influence the film cooling performance, has been extensively studied in the past decades.These studies confirmed the better performance of fan-shaped holes with respect to the cylindrical ones, due to diffusing expansion at the exit of the hole as mentioned by Goldstein et al. [1], and Gritsch et al. [2].However, the majority of these studies investigated film-cooling performance of fan-shaped holes in terms of film cooling effectiveness and heat transfer coefficient, (see for example Gritsch et al. [3,4]), whereas little attention has been paid to the flow field in the near hole-region which determines the thermal coverage.There are few experimental studies that measure the flow field in the mixing region between jet and mainstream for fan-shaped holes.Thole et al. [5] found that velocity gradients were reduced for the fan-shaped holes compared to the round holes, thus resulting in a decrease in turbulence levels downstream of the hole.Porter et al. [6] reported the velocity and turbulence data in the near field of round and laterally expanded fan-shaped holes.Their results are consistent with those discussed in [5].However, neither Thole et al. [5] nor Porter et al. [6] combined thermal and aerodynamic performance of film cooling.
From the numerical point of view, investigations of fan shaped holes are less common than those on cylindrical holes and, typically, make use of two-equation turbulence models to provide closure for steady RANS simulations.One of the most thorough numerical studies was carried out by Hyams and Leylek [7].They simulated four different shaped hole configurations, together with round holes as the reference case, using the k-ε model with generalized wall functions.They found that the laterally expanded fan-shaped holes ensure the highest film cooling effectiveness.Kohli and Thole [8] investigated the effect of inlet crossflow on cooling performance of a round and fan-shaped hole inclined at 35 • .This study revealed that, in case of shaped holes, there is no production of counter-rotating vortices at the exit so that turbulence peak occurs at the downstream edge of the hole.Porter et al. [6] also validated predictions of fan-shaped film cooling holes from two different Computational Fluid Dynamic (CFD) codes, ANSYS CFX 10.0 and HYDRA.Good agreement between numerical and experimental data was obtained in terms of centerline velocity profiles for a steady state simulation.Silieti et al. [9] performed a steady CFD modelling of a fan-shaped hole using three different turbulence models: The realizable k-ε (RKE) model, Menter SST k-ω (SST KW) as well as the v2-f turbulence model.Results were compared to experimental data in terms of centerline film cooling effectiveness and RKE resulted to be the most accurate turbulence model.On the other hand, some studies on cylindrical holes point to the fact that the Reynolds stress model is potentially the best model to predict anisotropic turbulence flow fields.Gustafsson and Johansson [10] investigate the flow field around a slanted jet in cross flow both experimentally and numerically.The Reynolds stress model (RSM) was found to be the best option for predicting the mean velocity flow field compared to RKE and SST KW.Moreover, it captures the features of the turbulent stresses reasonably well.Similarly, results of Javadi et al. [11] confirmed the capability of RSM to predict the flow field velocity and stress components.
This paper presents a detailed experimental characterization of film cooling jets discharged from a row of three fan-shaped holes, considering both thermal and aerodynamic aspects.Then, CFD simulations were performed and compared with experimental data.Three turbulence models were tested using STAR-CCM+, i.e., RKE, SST KW and RSM, to identify the most appropriate option for the current test case, through a severe process of validation against measurements.

Experimental Set-Up
The low-speed flat plate wind tunnel available at the Energy Systems and Turbomachinery Laboratory of the University of Bergamo was used for this study.Tests have been carried out at low speed, about 15 m/s, and low inlet turbulence intensity level below 1.0%, for blowing ratio (BR) values of 1 and 2. A detailed description of both wind tunnel and test conditions is provided in [12].The film cooling scheme consisted of a row of three fan-shaped holes with 5 mm diameter (Figure 1).Holes were sharp edge without using any fillets and inclined at 30 • to the surface; pitch-to-diameter ratio was P/D = 6 and length to diameter ratio was L/D = 8.The lateral expansion angle of the fan-shaped hole was 14 • , resulting in a hole width of 13.98 mm at the hole exit with t/P of 0.47 and a cross-sectional area ratio of 2.74.The chosen geometry of the fan-shaped hole was based on Gritsch et al. [2].In addition, Table 1 reports the cross-sectional planes that have been selected for the experimental measurements and their positions according to the reference coordinate system shown in Figures 1 and 2. It should be noted that origin of the reference coordinate system was located at the downstream edge of the middle hole (Z/D = 0).laterally expanded fan-shaped holes ensure the highest film cooling effectiveness.Kohli and Thole [8] investigated the effect of inlet crossflow on cooling performance of a round and fan-shaped hole inclined at 35°.This study revealed that, in case of shaped holes, there is no production of counterrotating vortices at the exit so that turbulence peak occurs at the downstream edge of the hole.Porter et al. [6] also validated predictions of fan-shaped film cooling holes from two different Computational Fluid Dynamic (CFD) codes, ANSYS CFX 10.0 and HYDRA.Good agreement between numerical and experimental data was obtained in terms of centerline velocity profiles for a steady state simulation.Silieti et al. [9] performed a steady CFD modelling of a fan-shaped hole using three different turbulence models: The realizable k-ε (RKE) model, Menter SST k-ω (SST KW) as well as the v2-f turbulence model.Results were compared to experimental data in terms of centerline film cooling effectiveness and RKE resulted to be the most accurate turbulence model.On the other hand, some studies on cylindrical holes point to the fact that the Reynolds stress model is potentially the best model to predict anisotropic turbulence flow fields.Gustafsson and Johansson [10] investigate the flow field around a slanted jet in cross flow both experimentally and numerically.The Reynolds stress model (RSM) was found to be the best option for predicting the mean velocity flow field compared to RKE and SST KW.Moreover, it captures the features of the turbulent stresses reasonably well.Similarly, results of Javadi et al. [11] confirmed the capability of RSM to predict the flow field velocity and stress components.This paper presents a detailed experimental characterization of film cooling jets discharged from a row of three fan-shaped holes, considering both thermal and aerodynamic aspects.Then, CFD simulations were performed and compared with experimental data.Three turbulence models were tested using STAR-CCM+, i.e., RKE, SST KW and RSM, to identify the most appropriate option for the current test case, through a severe process of validation against measurements.

Experimental Set-Up
The low-speed flat plate wind tunnel available at the Energy Systems and Turbomachinery Laboratory of the University of Bergamo was used for this study.Tests have been carried out at low speed, about 15 m/s, and low inlet turbulence intensity level below 1.0%, for blowing ratio (BR) values of 1 and 2. A detailed description of both wind tunnel and test conditions is provided in [12].The film cooling scheme consisted of a row of three fan-shaped holes with 5 mm diameter (Figure 1).Holes were sharp edge without using any fillets and inclined at 30° to the surface; pitch-to-diameter ratio was P/D = 6 and length to diameter ratio was L/D = 8.The lateral expansion angle of the fanshaped hole was 14°, resulting in a hole width of 13.98 mm at the hole exit with t/P of 0.47 and a crosssectional area ratio of 2.74.The chosen geometry of the fan-shaped hole was based on Gritsch et al. [2].In addition, Table 1 reports the cross-sectional planes that have been selected for the experimental measurements and their positions according to the reference coordinate system shown in Figures 1  and 2. It should be noted that origin of the reference coordinate system was located at the downstream edge of the middle hole (Z/D = 0).models, such as RKE, SST KW and RSM, were studied within the steady RANS approach.The working fluid was assumed to be incompressible air with temperature-dependent properties.Profiles of velocity and turbulence intensity, measured by the LDV technique, were applied at the inlet of the wind tunnel.The turbulent length scale for the mainstream was set to λ = 2.4 mm whereas, for the coolant, turbulence intensity (Tuc = 0.16•Rec −1/8 ) and length scale (λc = 0.07 Din) were determined from correlations of fully-developed pipe flow commonly used in CFD applications (see the user's guide by ANSYS Fluent [13]).Atmospheric pressure was set at the outlet.Smooth wall, no-slip boundary conditions were applied to the tunnel walls.A mass flow boundary condition was applied at the plenum inlet: A coolant mass flow rate of 0.000957 kg/s corresponds to BR = 1.Temperature for the mainstream and coolant flow was set at 288 K and 323 K, respectively.The computational domain, as shown in Figure 2, includes the main test section, the cooling hole and the plenum.In the spanwise direction, the domain extends from one mid-pitch plane to another, containing a full injection hole.All the residuals were kept below 10 −5 .1500 iterations were required to get convergence.

Computational Grid
The polyhedral mesh model was chosen in this study.The base mesh size was initially set to 0.03 m and refinement through the domain was controlled by growth parameters.Eight layers of prismatic cells were added near the wall, whose total height is 0.44D from the plate surface.Moreover, a cylindrical volumetric control was used to refine the mesh near the desired hole (Figure 3).Three grids were investigated for the sensitivity analysis, whose features are summarized in Table 2. Figure 4 compares the results from medium and fine grids, in terms of velocity profiles 5D downstream of injection holes, at BR = 1.Though the number of cells is quite different between these two grids, very good agreement is shown among the numerical velocity predictions.Accordingly,

Numerical Method and Boundary Conditions
RANS computations of a single jet injected through a fan-shaped hole into the mainstream were carried out at BR = 1 using the commercial CFD code, STAR-CCM+.Three different turbulence models, such as RKE, SST KW and RSM, were studied within the steady RANS approach.The working fluid was assumed to be incompressible air with temperature-dependent properties.Profiles of velocity and turbulence intensity, measured by the LDV technique, were applied at the inlet of the wind tunnel.The turbulent length scale for the mainstream was set to λ = 2.4 mm whereas, for the coolant, turbulence intensity (Tu c = 0.16•Re c −1/8 ) and length scale (λ c = 0.07 D in ) were determined from correlations of fully-developed pipe flow commonly used in CFD applications (see the user's guide by ANSYS Fluent [13]).Atmospheric pressure was set at the outlet.Smooth wall, no-slip boundary conditions were applied to the tunnel walls.A mass flow boundary condition was applied at the plenum inlet: A coolant mass flow rate of 0.000957 kg/s corresponds to BR = 1.Temperature for the mainstream and coolant flow was set at 288 K and 323 K, respectively.The computational domain, as shown in Figure 2, includes the main test section, the cooling hole and the plenum.In the span-wise direction, the domain extends from one mid-pitch plane to another, containing a full injection hole.All the residuals were kept below 10 −5 .1500 iterations were required to get convergence.

Computational Grid
The polyhedral mesh model was chosen in this study.The base mesh size was initially set to 0.03 m and refinement through the domain was controlled by growth parameters.Eight layers of prismatic cells were added near the wall, whose total height is 0.44D from the plate surface.Moreover, a cylindrical volumetric control was used to refine the mesh near the desired hole (Figure 3).Three grids were investigated for the sensitivity analysis, whose features are summarized in Table 2. Figure 4 compares the results from medium and fine grids, in terms of velocity profiles 5D downstream of injection holes, at BR = 1.Though the number of cells is quite different between these two grids, very good agreement is shown among the numerical velocity predictions.Accordingly, the medium grid was employed for all the RANS simulations.Since all turbulence models have been used with a two-layer zonal model for near-wall treatment, values of y + were on the order of unity.
Int. J. Turbomach.Propuls.Power 2019, 4, x FOR PEER REVIEW 4 of 13 the medium grid was employed for all the RANS simulations.Since all turbulence models have been used with a two-layer zonal model for near-wall treatment, values of y + were on the order of unity.

Experimental Results
Experimental results provided information regarding the aerodynamic and thermal behavior of the film cooling jet exiting three fan-shaped holes at blowing ratios of BR = 1 and 2. Aerodynamic investigations have been performed by means of laser doppler velocimetry (LDV), particle image velocimetry (PIV) and hot wire anemometer (HW) techniques.Concerning thermal measurements, the binary pressure sensitive paint (PSP) technique has been employed to measure adiabatic film cooling effectiveness.

Aerodynamic Aspects of Film Cooling
The aerodynamic results start with the characterization of the injection condition through measuring the discharge coefficient.Figure 5 shows the Cd distribution versus the ratio of coolant pressure pt,c to the free stream static pressure, p.As pressure ratio increases, the flow entering the  the medium grid was employed for all the RANS simulations.Since all turbulence models have been used with a two-layer zonal model for near-wall treatment, values of y + were on the order of unity.

Experimental Results
Experimental results provided information regarding the aerodynamic and thermal behavior of the film cooling jet exiting three fan-shaped holes at blowing ratios of BR = 1 and 2. Aerodynamic investigations have been performed by means of laser doppler velocimetry (LDV), particle image velocimetry (PIV) and hot wire anemometer (HW) techniques.Concerning thermal measurements, the binary pressure sensitive paint (PSP) technique has been employed to measure adiabatic film cooling effectiveness.

Aerodynamic Aspects of Film Cooling
The aerodynamic results start with the characterization of the injection condition through measuring the discharge coefficient.Figure 5 shows the Cd distribution versus the ratio of coolant pressure pt,c to the free stream static pressure, p.As pressure ratio increases, the flow entering the

Experimental Results
Experimental results provided information regarding the aerodynamic and thermal behavior of the film cooling jet exiting three fan-shaped holes at blowing ratios of BR = 1 and 2. Aerodynamic investigations have been performed by means of laser doppler velocimetry (LDV), particle image velocimetry (PIV) and hot wire anemometer (HW) techniques.Concerning thermal measurements, the binary pressure sensitive paint (PSP) technique has been employed to measure adiabatic film cooling effectiveness.

Aerodynamic Aspects of Film Cooling
The aerodynamic results start with the characterization of the injection condition through measuring the discharge coefficient.Figure 5 shows the C d distribution versus the ratio of coolant pressure p t,c to the free stream static pressure, p.As pressure ratio increases, the flow entering the diffuser separates more and more from the diffuser wall, thus leading to a reduction in the pressure recovery as well as the discharge coefficient.diffuser separates more and more from the diffuser wall, thus leading to a reduction in the pressure recovery as well as the discharge coefficient.The next step consisted in measuring the boundary layer downstream of the mid-hole at hole centerline by means of the LDV technique.Details on LDV measurements are given in [12].Results from LDV measurements at hole centerline are reported in Figures 6 and 7 in terms of profiles of time-averaged streamwise velocity component as well as streamwise and wall normal rms velocities at X/D = 1 and = 5, respectively (Figure 1).Data were normalized with respect to the mainstream velocity Ue.There was no indication of a strong shear layer at the hole exit, for both BR.Also, the velocity profile was relatively flat downstream of the hole.PIV technique (for more details refer again to Rouina et al. [12]) was used to characterize the mean and turbulent flow distribution over a stream-wise plane located at mid hole centerline (Z/D = 0).Normalized mean velocity flow fields in the near-hole region are shown in Figure 8, for BR of 1 The next step consisted in measuring the boundary layer downstream of the mid-hole at hole centerline by means of the LDV technique.Details on LDV measurements are given in [12].Results from LDV measurements at hole centerline are reported in Figures 6 and 7 in terms of profiles of time-averaged streamwise velocity component as well as streamwise and wall normal rms velocities at X/D = 1 and = 5, respectively (Figure 1).Data were normalized with respect to the mainstream velocity U e .There was no indication of a strong shear layer at the hole exit, for both BR.Also, the velocity profile was relatively flat downstream of the hole.diffuser separates more and more from the diffuser wall, thus leading to a reduction in the pressure recovery as well as the discharge coefficient.The next step consisted in measuring the boundary layer downstream of the mid-hole at hole centerline by means of the LDV technique.Details on LDV measurements are given in [12].Results from LDV measurements at hole centerline are reported in Figures 6 and 7 in terms of profiles of time-averaged streamwise velocity component as well as streamwise and wall normal rms velocities at X/D = 1 and = 5, respectively (Figure 1).Data were normalized with respect to the mainstream velocity Ue.There was no indication of a strong shear layer at the hole exit, for both BR.Also, the velocity profile was relatively flat downstream of the hole.PIV technique (for more details refer again to Rouina et al. [12]) was used to characterize the mean and turbulent flow distribution over a stream-wise plane located at mid hole centerline (Z/D = 0).Normalized mean velocity flow fields in the near-hole region are shown in Figure 8, for BR of 1 diffuser separates more and more from the diffuser wall, thus leading to a reduction in the pressure recovery as well as the discharge coefficient.The next step consisted in measuring the boundary layer downstream of the mid-hole at hole centerline by means of the LDV technique.Details on LDV measurements are given in [12].Results from LDV measurements at hole centerline are reported in Figures 6 and 7 in terms of profiles of time-averaged streamwise velocity component as well as streamwise and wall normal rms velocities at X/D = 1 and = 5, respectively (Figure 1).Data were normalized with respect to the mainstream velocity Ue.There was no indication of a strong shear layer at the hole exit, for both BR.Also, the velocity profile was relatively flat downstream of the hole.PIV technique (for more details refer again to Rouina et al. [12]) was used to characterize the mean and turbulent flow distribution over a stream-wise plane located at mid hole centerline (Z/D = 0).Normalized mean velocity flow fields in the near-hole region are shown in Figure 8, for BR of 1 PIV technique (for more details refer again to Rouina et al. [12]) was used to characterize the mean and turbulent flow distribution over a stream-wise plane located at mid hole centerline (Z/D = 0).Normalized mean velocity flow fields in the near-hole region are shown in Figure 8, for BR of 1 and 2. It is evident that the coolant jet remained attached to the surface, in agreement with LDV results, even at the higher BR.It is worth mentioning that, increasing the BR to 2 did not lead to jet separation from the flat plate surface.
Int. J. Turbomach.Propuls.Power 2019, 4, x FOR PEER REVIEW 6 of 13 and 2. It is evident that the coolant jet remained attached to the surface, in agreement with LDV results, even at the higher BR.It is worth mentioning that, increasing the BR to 2 did not lead to jet separation from the flat plate surface.Flow visualizations (see again [12]) were also performed on the same plane to investigate the mixing process between the coolant and the mainstream (Figure 9).The unsteady nature of the jet interaction with the mainstream is documented by the presence of large coherent structures at the jet boundaries.

Thermal Aspects of Film Cooling
PSP technique was used to measure adiabatic film cooling effectiveness (η) on the surface of the flat plate (see Abdeh et al. [14]).Contours of η obtained by injecting nitrogen as the coolant stream, are depicted in Figure 10.For both BR, the jet is attached to the flat plate surface and spreads out laterally to ensure a more effective thermal protection compared to the cylindrical hole [12].The centerline effectiveness, as a function of the distance from the hole, indicates that the fan-shaped hole provides values of ηcl > 0.7 within X/D < 5, for both BR (Figure 11).Conversely, the BR influences the ηcl decay at X/D > 7: The higher BR is more effective than the lower one at downstream locations.Figure 12 shows the laterally averaged film-cooling effectiveness (ηave) plotted versus the streamwise distance from the hole exit.Effectiveness decreases monotonously with the streamwise distance from the hole, for both BR, but with different slope.The highest ηave was found at BR = 1, in the near vicinity of the hole.Doubling the BR caused a slight reduction in ηave within X/D < 7, despite some minor improvements in ηave far from the hole exit.Flow visualizations (see again [12]) were also performed on the same plane to investigate the mixing process between the coolant and the mainstream (Figure 9).The unsteady nature of the jet interaction with the mainstream is documented by the presence of large coherent structures at the jet boundaries.and 2. It is evident that the coolant jet remained attached to the surface, in agreement with LDV results, even at the higher BR.It is worth mentioning that, increasing the BR to 2 did not lead to jet separation from the flat plate surface.Flow visualizations (see again [12]) were also performed on the same plane to investigate the mixing process between the coolant and the mainstream (Figure 9).The unsteady nature of the jet interaction with the mainstream is documented by the presence of large coherent structures at the jet boundaries.

Thermal Aspects of Film Cooling
PSP technique was used to measure adiabatic film cooling effectiveness (η) on the surface of the flat plate (see Abdeh et al. [14]).Contours of η obtained by injecting nitrogen as the coolant stream, are depicted in Figure 10.For both BR, the jet is attached to the flat plate surface and spreads out laterally to ensure a more effective thermal protection compared to the cylindrical hole [12].The centerline effectiveness, as a function of the distance from the hole, indicates that the fan-shaped hole provides values of ηcl > 0.7 within X/D < 5, for both BR (Figure 11).Conversely, the BR influences the ηcl decay at X/D > 7: The higher BR is more effective than the lower one at downstream locations.Figure 12 shows the laterally averaged film-cooling effectiveness (ηave) plotted versus the streamwise distance from the hole exit.Effectiveness decreases monotonously with the streamwise distance from the hole, for both BR, but with different slope.The highest ηave was found at BR = 1, in the near vicinity of the hole.Doubling the BR caused a slight reduction in ηave within X/D < 7, despite some minor improvements in ηave far from the hole exit.

Thermal Aspects of Film Cooling
PSP technique was used to measure adiabatic film cooling effectiveness (η) on the surface of the flat plate (see Abdeh et al. [14]).Contours of η obtained by injecting nitrogen as the coolant stream, are depicted in Figure 10.For both BR, the jet is attached to the flat plate surface and spreads out laterally to ensure a more effective thermal protection compared to the cylindrical hole [12].The centerline effectiveness, as a function of the distance from the hole, indicates that the fan-shaped hole provides values of η cl > 0.7 within X/D < 5, for both BR (Figure 11).Conversely, the BR influences the η cl decay at X/D > 7: The higher BR is more effective than the lower one at downstream locations.Figure 12 shows the laterally averaged film-cooling effectiveness (η ave ) plotted versus the streamwise distance from the hole exit.Effectiveness decreases monotonously with the streamwise distance from the hole, for both BR, but with different slope.The highest η ave was found at BR = 1, in the near vicinity of the hole.Doubling the BR caused a slight reduction in η ave within X/D < 7, despite some minor improvements in η ave far from the hole exit.

Numerical Results
A preliminary RANS-based CFD study was conducted to validate and compare the above mentioned three turbulence models in order to determine the best option for predicting the film cooling jet behavior at BR = 1, i.e., the condition providing the highest effectiveness while minimizing the cooling air consumption.

Turbulence Model Assessment
Laterally averaged and centerline film cooling effectiveness predicted by RKE, SST KW and RSM turbulence models were plotted and compared to PSP measurements in Figures 13 and 14

Numerical Results
A preliminary RANS-based CFD study was conducted to validate and compare the above mentioned three turbulence models in order to determine the best option for predicting the film cooling jet behavior at BR = 1, i.e., the condition providing the highest effectiveness while minimizing the cooling air consumption.

Turbulence Model Assessment
Laterally averaged and centerline film cooling effectiveness predicted by RKE, SST KW and RSM turbulence models were plotted and compared to PSP measurements in Figures 13 and 14

Numerical Results
A preliminary RANS-based CFD study was conducted to validate and compare the above mentioned three turbulence models in order to determine the best option for predicting the film cooling jet behavior at BR = 1, i.e., the condition providing the highest effectiveness while minimizing the cooling air consumption.

Turbulence Model Assessment
Laterally averaged and centerline film cooling effectiveness predicted by RKE, SST KW and RSM turbulence models were plotted and compared to PSP measurements in Figures 13 and 14

Numerical Results
A preliminary RANS-based CFD study was conducted to validate and compare the above mentioned three turbulence models in order to determine the best option for predicting the film cooling jet behavior at BR = 1, i.e., the condition providing the highest effectiveness while minimizing the cooling air consumption.

Turbulence Model Assessment
Laterally averaged and centerline film cooling effectiveness predicted by RKE, SST KW and RSM turbulence models were plotted and compared to PSP measurements in Figures 13 and 14, respectively.The decreasing trend of film cooling effectiveness was best predicted by RSM and RKE.Surprisingly, results from the SST KW model were completely different from the experimental data as well as from the predictions by the other two models.This was due to SST KW tendency to overestimate separation of the cooling jet, thus leading to a narrower trace on the surface and, consequently, lower η. Figure 13 highlights that RSM predictions of η ave are the closest to PSP data: They are within 4% of measurements up to X/D = 4 whereas the deviation from the experimental data raises to about 15% at X/D = 18 (measured η ave = 0.195 against a predicted η ave = 0.225).Figure 14 shows that RSM slightly over-predicts η cl while RKE under-predicts it, within X/D < 9.At X/D > 12, RKE and RSM profiles of η cl collapse onto a single curve, thereby over-predicting the measured data.
respectively.The decreasing trend of film cooling effectiveness was best predicted by RSM and RKE.Surprisingly, results from the SST KW model were completely different from the experimental data as well as from the predictions by the other two models.This was due to SST KW tendency to overestimate separation of the cooling jet, thus leading to a narrower trace on the surface and, consequently, lower η. Figure 13 highlights that RSM predictions of ηave are the closest to PSP data: They are within 4% of measurements up to X/D = 4 whereas the deviation from the experimental data raises to about 15% at X/D = 18 (measured ηave = 0.195 against a predicted ηave = 0.225).Figure 14 shows that RSM slightly over-predicts ηcl while RKE under-predicts it, within X/D < 9.At X/D > 12, RKE and RSM profiles of ηcl collapse onto a single curve, thereby over-predicting the measured data.
Similarly, Figure 15 contains numerical profiles of the normalized streamwise velocity at 1D and 5D downstream of the holes for comparison against PIV measurements.RSM provides the most accurate prediction of U/Ue at 1D downstream of hole.Moving to X = 5D, although all three turbulence models predicted a steeper velocity profile compared to measurements, RSM and RKE results were anyhow reliable.Please note that flow separation occurred on the leeward side of the hole, according with RSM prediction.Considering both the reported thermal and aerodynamic aspects, one could conclude that the RSM is the most appropriate model to choose for fan-shaped holes in this set of experiments, thus minimizing the error in the prediction of ηave and U/Ue, at least in the near vicinity of the hole exit (X/D = 1).respectively.The decreasing trend of film cooling effectiveness was best predicted by RSM and RKE.Surprisingly, results from the SST KW model were completely different from the experimental data as well as from the predictions by the other two models.This was due to SST KW tendency to overestimate separation of the cooling jet, thus leading to a narrower trace on the surface and, consequently, lower η. Figure 13 highlights that RSM predictions of ηave are the closest to PSP data: They are within 4% of measurements up to X/D = 4 whereas the deviation from the experimental data raises to about 15% at X/D = 18 (measured ηave = 0.195 against a predicted ηave = 0.225).Figure 14 shows that RSM slightly over-predicts ηcl while RKE under-predicts it, within X/D < 9.At X/D > 12, RKE and RSM profiles of ηcl collapse onto a single curve, thereby over-predicting the measured data.
Similarly, Figure 15 contains numerical profiles of the normalized streamwise velocity at 1D and 5D downstream of the holes for comparison against PIV measurements.RSM provides the most accurate prediction of U/Ue at 1D downstream of hole.Moving to X = 5D, although all three turbulence models predicted a steeper velocity profile compared to measurements, RSM and RKE results were anyhow reliable.Please note that flow separation occurred on the leeward side of the hole, according with RSM prediction.Considering both the reported thermal and aerodynamic aspects, one could conclude that the RSM is the most appropriate model to choose for fan-shaped holes in this set of experiments, thus minimizing the error in the prediction of ηave and U/Ue, at least in the near vicinity of the hole exit (X/D = 1).Similarly, Figure 15 contains numerical profiles of the normalized streamwise velocity at 1D and 5D downstream of the holes for comparison against PIV measurements.RSM provides the most accurate prediction of U/U e at 1D downstream of hole.Moving to X = 5D, although all three turbulence models predicted a steeper velocity profile compared to measurements, RSM and RKE results were anyhow reliable.Please note that flow separation occurred on the leeward side of the hole, according with RSM prediction.Considering both the reported thermal and aerodynamic aspects, one could conclude that the RSM is the most appropriate model to choose for fan-shaped holes in this set of experiments, thus minimizing the error in the prediction of η ave and U/U e , at least in the near vicinity of the hole exit (X/D = 1).

RSM Validation Against Experimental Data
Figure 16 presents the η contour predicted by RSM for validation against those shown in Figure 10a.Simulation results were in good agreement with the experimental one.The jet stays close to the surface according with the RSM model leading to high value of η in the near hole region, consistently with the PSP data.Considering the lateral spreading of coolant, the experimental pattern of η shows a slight narrowing of the coolant footprint on the surface, at increasing distance from the hole exit, whereas the RSM model provided a constant width coolant trace travelling in the stream wise direction.This indicates that lateral mixing was underestimated by RSM simulation, consistently with the above-mentioned profiles of ηave.
Special care was taken to assess the potential of RSM when dealing with off-the-wall features of the flow in the mixing region.2D planes of normalized streamwise velocity flow field obtained from PIV measurements and RSM simulation are reported in Figure 17.Again, the measured velocity field matched well with RSM prediction.The same conclusion can be drawn looking at Figure 18 where the streamwise velocity in a vertical-lateral plane at 5D downstream of the hole is reported from both RSM and HW.However, it seems that RSM slightly overpredicts the vertical spreading of the jet when comparing to HW results, notwithstanding the impossibility of collecting measurements very close to the wall.

RSM Validation Against Experimental Data
Figure 16 presents the η contour predicted by RSM for validation against those shown in Figure 10a.Simulation results were in good agreement with the experimental one.The jet stays close to the surface according with the RSM model leading to high value of η in the near hole region, consistently with the PSP data.Considering the lateral spreading of coolant, the experimental pattern of η shows a slight narrowing of the coolant footprint on the surface, at increasing distance from the hole exit, whereas the RSM model provided a constant width coolant trace travelling in the stream wise direction.This indicates that lateral mixing was underestimated by RSM simulation, consistently with the above-mentioned profiles of η ave .

RSM Validation Against Experimental Data
Figure 16 presents the η contour predicted by RSM for validation against those shown in Figure 10a.Simulation results were in good agreement with the experimental one.The jet stays close to the surface according with the RSM model leading to high value of η in the near hole region, consistently with the PSP data.Considering the lateral spreading of coolant, the experimental pattern of η shows a slight narrowing of the coolant footprint on the surface, at increasing distance from the hole exit, whereas the RSM model provided a constant width coolant trace travelling in the stream wise direction.This indicates that lateral mixing was underestimated by RSM simulation, consistently with the above-mentioned profiles of ηave.
Special care was taken to assess the potential of RSM when dealing with off-the-wall features of the flow in the mixing region.2D planes of normalized streamwise velocity flow field obtained from PIV measurements and RSM simulation are reported in Figure 17.Again, the measured velocity field matched well with RSM prediction.The same conclusion can be drawn looking at Figure 18 where the streamwise velocity in a vertical-lateral plane at 5D downstream of the hole is reported from both RSM and HW.However, it seems that RSM slightly overpredicts the vertical spreading of the jet when comparing to HW results, notwithstanding the impossibility of collecting measurements very close to the wall.Special care was taken to assess the potential of RSM when dealing with off-the-wall features of the flow in the mixing region.2D planes of normalized streamwise velocity flow field obtained from PIV measurements and RSM simulation are reported in Figure 17.Again, the measured velocity field matched well with RSM prediction.The same conclusion can be drawn looking at Figure 18 where the streamwise velocity in a vertical-lateral plane at 5D downstream of the hole is reported from both RSM and HW.However, it seems that RSM slightly overpredicts the vertical spreading of the jet when comparing to HW results, notwithstanding the impossibility of collecting measurements very close to the wall.Computed and measured normalized stress profiles along the jet center plane (Z/D = 0) are shown in Figure 19 at X/D = 1 and X/D = 5, for BR = 1.On the one hand, RSM predicts the experimental trend of stress components accurately.On the other hand, it provides lower levels of streamwise stress components close to the wall than expected.Actually, the level of the normal stresses in the RSM model is controlled by their respective production terms, pressure-strain correlation terms and dissipation rate correlation.As documented by Hoda et al. [15], the failure of RSM to capture the experimental trends specifically in the wake of the jet is largely due to under predicted production levels.
Please also note that Figure 19 shows the experimental stress profiles at BR = 2, with the aim of assessing the influence of increasing BR on rms velocities.As expected, higher turbulent activity can be found inside the shear layer between the coolant and the mainstream when doubling the BR, with streamwise velocity fluctuations as large as 50% of the freestream velocity.A noteworthy degree of anisotropy is also evident, thus emphasizing the importance of using RSM.Computed and measured normalized stress profiles along the jet center plane (Z/D = 0) are shown in Figure 19 at X/D = 1 and X/D = 5, for BR = 1.On the one hand, RSM predicts the experimental trend of stress components accurately.On the other hand, it provides lower levels of streamwise stress components close to the wall than expected.Actually, the level of the normal stresses in the RSM model is controlled by their respective production terms, pressure-strain correlation terms and dissipation rate correlation.As documented by Hoda et al. [15], the failure of RSM to capture the experimental trends specifically in the wake of the jet is largely due to under predicted production levels.
Please also note that Figure 19 shows the experimental stress profiles at BR = 2, with the aim of assessing the influence of increasing BR on rms velocities.As expected, higher turbulent activity can be found inside the shear layer between the coolant and the mainstream when doubling the BR, with streamwise velocity fluctuations as large as 50% of the freestream velocity.A noteworthy degree of anisotropy is also evident, thus emphasizing the importance of using RSM.Computed and measured normalized stress profiles along the jet center plane (Z/D = 0) are shown in Figure 19 at X/D = 1 and X/D = 5, for BR = 1.On the one hand, RSM predicts the experimental trend of stress components accurately.On the other hand, it provides lower levels of streamwise stress components close to the wall than expected.Actually, the level of the normal stresses in the RSM model is controlled by their respective production terms, pressure-strain correlation terms and dissipation rate correlation.As documented by Hoda et al. [15], the failure of RSM to capture the experimental trends specifically in the wake of the jet is largely due to under predicted production levels.
Please also note that Figure 19 shows the experimental stress profiles at BR = 2, with the aim of assessing the influence of increasing BR on rms velocities.As expected, higher turbulent activity can be found inside the shear layer between the coolant and the mainstream when doubling the BR, with streamwise velocity fluctuations as large as 50% of the freestream velocity.A noteworthy degree of anisotropy is also evident, thus emphasizing the importance of using RSM.

Conclusions
The interaction between the mainstream and the film cooling jets exiting three fan-shaped holes has been investigated experimentally and numerically in terms of aerodynamic and thermal aspects, by using a flat plate wind tunnel.The operating conditions included low-speed and low inlet turbulence intensity, with BR equal to 1 and 2. Experimental investigations confirmed previous data

Conclusions
The interaction between the mainstream and the film cooling jets exiting three fan-shaped holes has been investigated experimentally and numerically in terms of aerodynamic and thermal aspects, by using a flat plate wind tunnel.The operating conditions included low-speed and low inlet turbulence intensity, with BR equal to 1 and 2. Experimental investigations confirmed previous data available in the open literature regarding fan-shaped holes, i.e., the coolant jet stayed attached to the flat plate surface even at the higher BR of 2. Concerning numerical investigations, predictions from RKE, SST KW and RSM were compared against measurements of laterally and centerline adiabatic effectiveness, as well as off-the-wall velocity maps and profiles of stress components.RSM provided the most accurate predictions, considering both aerodynamic and thermal aspects of the coolant interaction with the mainstream.Solving individual components of the Reynolds stress tensor was rewarding to simulate the current set of experiments, as compared to the isotropic eddy viscosity assumption, while saving computational effort compared to scale resolving simulations.

Figure 3 .
Figure 3. Prism layers and near wall mesh refinement through cylindrical volumetric control.

Figure 3 .
Figure 3. Prism layers and near wall mesh refinement through cylindrical volumetric control.

Figure 3 .
Figure 3. Prism layers and near wall mesh refinement through cylindrical volumetric control.
GridBase Size (m) Volumetric Control Grid Size (m) y + Cells Count
GridBase Size (m) Volumetric Control Grid Size (m) y + Cells Count