Effect of blowing rate on the film cooling coverage on a multi-holed plate: application on combustor walls

Similar documents
LDA-Measurements of Jets in Crossflow for Effusion Cooling Applications

THE EFFECT OF SAMPLE SIZE, TURBULENCE INTENSITY AND THE VELOCITY FIELD ON THE EXPERIMENTAL ACCURACY OF ENSEMBLE AVERAGED PIV MEASUREMENTS

The Effect of Endplates on Rectangular Jets of Different Aspect Ratios

University of Maiduguri Faculty of Engineering Seminar Series Volume 6, december Seminar Series Volume 6, 2015 Page 58

Evolution of the pdf of a high Schmidt number passive scalar in a plane wake

UNIT II CONVECTION HEAT TRANSFER

Vortex shedding from slender surface mounted pyramids

Study on Impingement of Air Jet from Orifice on Convex Surface for Unconfined Flow

AEROSPACE ENGINEERING DEPARTMENT. Second Year - Second Term ( ) Fluid Mechanics & Gas Dynamics

Applied Fluid Mechanics

Comparison Of Square-hole And Round-hole Film Cooling: A Computational Study

86400 Batu Pahat, Johor, Malaysia. Iwate University, Japan

Applied Fluid Mechanics

Full-Coverage Film Cooling With Short Normal Injection Holes

PDA Measurements of Single Point Injection in Cross-flow

Intensely swirling turbulent pipe flow downstream of an orifice: the influence of an outlet contraction

Journal of Fluid Science and Technology

Chapter 7: External Forced Convection. Dr Ali Jawarneh Department of Mechanical Engineering Hashemite University

PIV study for the analysis of planar jets in cross-flow at low Reynolds number

Experimental Investigation of Adiabatic Film Cooling Effectiveness and Heat Transfer Coefficients over a Gas Turbine Blade Leading Edge Configuration

Unsteady Volumetric Entropy Generation Rate in Laminar Boundary Layers

Chapter 3 NATURAL CONVECTION

Experimental Investigations on the Local Distribution of wall static pressure coefficient Due To an Impinging Slot Air Jet on a Confined Rough Surface

Effects of the Leakage Flow Tangential Velocity in Shrouded Axial Compressor Cascades *

Module 3: Velocity Measurement Lecture 16: Validation of PIV with HWA. The Lecture Contains: Hotwire Anemometry. Uncertainity

Experimental Study of the Holes Density Effect on Film Cooling Performance

Numerical Validation of Flow Through an S-shaped Diffuser

ASSESSMENT OF ANISOTROPY IN THE NEAR FIELD OF A RECTANGULAR TURBULENT JET

ENGR Heat Transfer II

Experimental Verification of CFD Modeling of Turbulent Flow over Circular Cavities using FLUENT

Principles of Convection

Periodic planes v i+1 Top wall u i. Inlet. U m y. Jet hole. Figure 2. Schematic of computational domain.

Fluctuating Heat Transfer to an Impinging Air Jet in the Transitional Wall Jet Region

Heat and Mass Transfer Prof. S.P. Sukhatme Department of Mechanical Engineering Indian Institute of Technology, Bombay

6. Laser Doppler Anemometry. Introduction to principles and applications

ENERGY PERFORMANCE IMPROVEMENT, FLOW BEHAVIOR AND HEAT TRANSFER INVESTIGATION IN A CIRCULAR TUBE WITH V-DOWNSTREAM DISCRETE BAFFLES

Jet Aircraft Propulsion Prof. Bhaskar Roy Prof. A.M. Pradeep Department of Aerospace Engineering

NUMERICAL HEAT TRANSFER ENHANCEMENT IN SQUARE DUCT WITH INTERNAL RIB

PIV measurements and convective heat transfer of an impinging air jet

Maximum Heat Transfer Density From Finned Tubes Cooled By Natural Convection

Lecture 22. Mechanical Energy Balance

Effects of Free-Stream Vorticity on the Blasius Boundary Layer

PROOF COPY JTM. Oguz Uzol Cengiz Camci Turbomachinery Heat Transfer Laboratory, The Pennsylvania State University, University Park, PA 16802

EFFECT OF VELOCITY RATIO ON FLOW AND HEAT TRANSFER CHARACTERISTICS OF AN IMPINGING JET IN CROSSFLOW

Defense Technical Information Center Compilation Part Notice

CHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE

ME332 FLUID MECHANICS LABORATORY (PART I)

External Forced Convection. Copyright The McGraw-Hill Companies, Inc. Permission required for reproduction or display.

Standard Practices for Air Speed Calibration Testing

EFFECT OF VORTICES ON JET IMPINGEMENT HEAT TRANSFER. Abstract

CONVECTION HEAT TRANSFER

FLOW CONTROL USING DBD PLASMA ON BACKWARD-FACING STEP

Heat Transfer F12-ENG Lab #4 Forced convection School of Engineering, UC Merced.

GAS TURBINE ENGINE DURABILITY IMPACTS OF HIGH FUEL-AIR RATIO COMBUSTORS: NEAR WALL REACTION EFFECTS ON FILM-COOLED BACKWARD-FACING STEP HEAT TRANSFER

Numerical and Experimental Study on the Effect of Guide Vane Insertion on the Flow Characteristics in a 90º Rectangular Elbow

MODELLING PARTICLE DEPOSITION ON GAS TURBINE BLADE SURFACES

Lecture 30 Review of Fluid Flow and Heat Transfer

HIGH INTENSITY, LARGE LENGTH-SCALE FREESTREAM TURBULENCE GENERATION IN A TRANSONIC TURBINE CASCADE

Active Control of Turbulence and Fluid- Structure Interactions

White Paper FINAL REPORT AN EVALUATION OF THE HYDRODYNAMICS MECHANISMS WHICH DRIVE THE PERFORMANCE OF THE WESTFALL STATIC MIXER.

Journal of Fluid Science and Technology

TUBE BANKS TEST PROBLEMS

COMPUTATIONAL FLUID DYNAMICS ANALYSIS OF A V-RIB WITH GAP ROUGHENED SOLAR AIR HEATER

Advances in Fluid Mechanics and Heat & Mass Transfer

IEEE TRANSACTIONS ON COMPONENTS, PACKAGING, AND MANUFACTURING TECHNOLOGY PART A, VOL. 20, NO. 4, DECEMBER

Boundary Layer Transition on the Suction Side of a Turbine Blade

CFD Analysis for Thermal Behavior of Turbulent Channel Flow of Different Geometry of Bottom Plate

10.52 Mechanics of Fluids Spring 2006 Problem Set 3

Flow Structure Investigations in a "Tornado" Combustor

GT Proceedings of ASME Turbo Expo 2013: Turbine Technical Conference and Exposition GT2013 June 3-7, 2013, San Antonio, Texas, USA

An adiabatic homogeneous model for the flow around a multi-perforated plate

INTERACTION OF AN AIR-BUBBLE DISPERSED PHASE WITH AN INITIALLY ISOTROPIC TURBULENT FLOW FIELD

An alternative turbulent heat flux modelling for gas turbine cooling application

GENERALISATION OF THE TWO-SCALE MOMENTUM THEORY FOR COUPLED WIND TURBINE/FARM OPTIMISATION

NUMERICAL SIMULATION OF STATIC INFLOW DISTORTION ON AN AXIAL FLOW FAN

3D hot-wire measurements of a wind turbine wake

NUMERICAL STUDY OF EFFUSION COOLING FLOW AND HEAT TRANSFER

Optimization of the Gas Flow in a GEM Tracker with COMSOL and TENDIGEM Development. Presented at the 2011 COMSOL Conference

Numerical simulation of fluid flow in a monolithic exchanger related to high temperature and high pressure operating conditions

CONVECTION HEAT TRANSFER

FLOW MALDISTRIBUTION IN A SIMPLIFIED PLATE HEAT EXCHANGER MODEL - A Numerical Study

Effect of the Liquid Injection Angle on the Atomization of Liquid Jets in Subsonic Crossflows

An Experimental Investigation to Control the Flow Emerging From a Wide Angle Diffuser

Meysam ATASHAFROOZ, Seyyed Abdolreza GANDJALIKHAN NASSAB, and Amir Babak ANSARI

Heat Transfer Enhancement of Solar Air Heater By Using Artificial Roughness double inclined ribs

Development of Velocity Profile Generating Screens for Gas Turbine Components

Local blockage effect for wind turbines

Effect of roughness shape on heat transfer and flow friction characteristics of solar air heater with roughened absorber plate

GÖRTLER VORTICES AND THEIR EFFECT ON HEAT TRANSFER

Hydrodynamic Measurements of the Flow Structure Emanating From A Multi-Row Film Cooling Configuration

Colloquium FLUID DYNAMICS 2013 Institute of Thermomechanics AS CR, v.v.i., Prague, October 23-25, 2013 p.1

APPLIED FLUID DYNAMICS HANDBOOK

TURBULENCE AND PRESSURE DROP BEHAVIORS AROUND SEMICIRCULAR RIBS IN A RECTANGULAR CHANNEL

ON IMPROVING FILM COOLING EFFECTIVENESS OF COMBUSTOR LINER PLATES OF GAS TURBINES BY USING PLATES

Convection. forced convection when the flow is caused by external means, such as by a fan, a pump, or atmospheric winds.

Heat Exchangers for Condensation and Evaporation Applications Operating in a Low Pressure Atmosphere

Infrared measurements of heat transfer in jet impingement on concave wall applied to anti-icing

Fall 2014 Qualifying Exam Thermodynamics Closed Book

Detailed Outline, M E 320 Fluid Flow, Spring Semester 2015

Keywords - Gas Turbine, Exhaust Diffuser, Annular Diffuser, CFD, Numerical Simulations.

Transcription:

Effect of blowing rate on the film cooling coverage on a multi-holed plate: application on combustor walls P. Miron 1,2, C. Berat 1 & V. Sabelnikov 3 1 TURBOMECA-Bordes, France 2 LaTEP, Université de Pau et des Pays de l'adour - France 3 ONERA-Palaiseau, France Abstract Multi-holed film cooling systems are the most effective and widely used cooling techniques applied to gas turbine blades and combustor liners to obtain resistance in extremely high temperature operating conditions. This paper reports on measurements in flows generated by jets issuing from a 10:1 scaled multi-holed plate for a blowing ratio varying from 0.2 to 11. The plate geometry is similar to geometries encountered within current combustion chamber walls. The test section has 138 holes arranged on 12 staggered rows. The holes are spaced 6.74 diameters apart in the spanwise direction and 5.84 diameters apart in the downstream direction. The holes angle inclination is 30 degrees with respect to the plate surface. Mean velocity vectors and turbulent Reynolds shear stresses are measured using a two-component Laser Doppler Anemometry system. An investigation is performed to study the effects of adjacent jets flow interactions and jet-to-cross-flow combinations on the film cooling coverage, downstream from the end of the perforated zone. Results obtained are intended to help improve the aerodynamic and thermodynamic heat-transfer laws for multi-holed and near wall flows. The acquired experimental database is useful to validate different physical models and numerical codes in multi-holed flow cooling techniques. Keywords: LDA, perforated plate, combustion chambers, blowing rate, turbulence intensity, Reynolds shear stresses, mean flow velocity, film cooling. é Heat Transfer VIII, B. Sund

508 Heat Transfer VIII 1 Introduction Cooling of the hot section in gas turbine engines has become more difficult due to the rise of compressor exit temperature that has doubled over the last 30 years. In addition, gas turbine manufacturers are also striving to reach higher component lifetimes. These two factors lead to the utilisation of state of the art wall cooling methods. Film cooling techniques are used for various gas turbine components, including combustion chambers, disc cavities and turbine blades. Multi-holed plates are used in different forms, such as impinging jets, wall jets, or cross-flow jets. In the film cooling application, usually, arrays of discrete holes are used to provide the cold air flow. In fact, only a few percent of the cold air coming from compressor, passes through the perforated regions of the combustor walls. The rest of the cold air is conducted to the end of the combustion chamber, where is used for the dilution of the burned gases, and also for the post-combustion process. In the perforated regions, jets requires a distance to merge together and to form a uniform cooling film. Before the jet flows merge, hot spots exist between the jets, therefore even downstream of the jet-merging location, cooling effectiveness is not necessarily uniform. This nonuniformity may lead to thermal fatigue and thermal barrier coating spalling. Today experiments are important to validate predictive methods for turbulent flows and especially detailed measurements of the velocity distribution in the near wall flow field are much needed for the improvement of turbulence models and the development of physically correct models. There are many film cooling studies using a single hole (Fric and Roshko [1]) or one to three rows of holes (Goldstein [2], Brown and Saluja [3]). However, little research has been performed using staggered rows of holes with a realistic cold air supply method for the perforated plate as correctly simulated by our experimental setup. 2 Experimental facility and instrumentation A special wind tunnel working at ambient temperature was built to investigate aerodynamically the film cooling through a multi-holed plate as an application of combustor walls cooling techniques. Two multi-holed plates, 60 and 30 degrees holes angle inclination with respect to plate surface, are available and can be mounted into the test section. In this study, only the results of the second perforated plate are shown. The low-speed wind tunnel consists of four parts : the main flow blowers, the development region, the test section and the exhaust system as described next in detail. 2.1 Main flow blowers, development flow region and exhaust system The experiments are conducted in an open-circuit, blowing-type wind tunnel as shown in Figure 1. Two identical fan blowers provide the primary and the secondary flows through the two rectangular ducts. Before entering these ducts, the two flows pass through a 50-cm contraction segment with a ratio of 2:1 that contains a 2 cm long honeycomb section and a screen to provide uniform flows.

Heat Transfer VIII 509 Next, two rectangular ducts (40 cm x 12 cm) of 2.5 m in length are connected to the test section. Flow rates through the ducts are controlled by changing the rotational speed of the blowers and a calibration with the LDA measurements velocity profiles is done. At the end of the test section, a pressure screen is used to create a pressure difference to facilitate the passage through holes into the primary flow, and therefore, to create a film cooling downstream the end of the perforated zone. Depending on the secondary flow rate and the flows pressure difference, a small part of the secondary flow entering the test section passes through inclined holes, and the rest exits, through the pressure screen, to the exhaust system. The primary flow passes through the test section at the atmospheric pressure and then directly into the exhaust system. A smoke generator, using white F5-Smoggy fluid is seeding both flows at the inlet blowers to ensure homogenous particle concentration in the test section. Figure 1: Experimental setup. 2.2 Test section The test section was designed to duplicate existing manufacturing and operating conditions of perforated plates in today s combustion chambers. Two groups of parameters, geometrical and aerodynamic, were investigated to obtain the optimum configuration. The selection process is explained in the following sections. 2.2.1 Geometrical parameters Many factors affect how the mixing process of the jets occurs in the film cooling region of the perforated plate. Some of the geometric parameters, (see Figure 2) like the shape of the jet issuing holes, the diameter of the jet (d), the inclination angle (α), the spacing between two adjacent jets (s), the inter-row distance (r), and the plate thickness (e), form the design parameters with respect to the jets. These variables require a large number of combinations to be investigated. But, since our perforated plate is a 10:1 test scaled model with respect to a real perforated plate of combustor walls, the holes diameter is fixed, and the plate thickness is two-times the hole diameter. The perforated zone has 138 holes arranged on 12 staggered rows starting 20 diameters downstream the beginning of the test section. The step in the spanwise direction is s = 6.74 d and r = 5.84 d in the downstream direction. The inclination angle of the holes is α = 30 degrees with respect to the plate surface.

510 Heat Transfer VIII Figure 2: Perforated plate geometric parameters. 2.2.2 Aerodynamic flow parameters The flow parameters were identified as the blowing ratio M = ρ jet V jet / ρ 1 V 1, the maximum pressure difference between secondary flow and primary flow P = (P 2 P 1 ) [Pa], and the Reynolds number based on the duct centerline velocity and the half height of the rectangular duct Re = (V 1 H/2) / ν. In order to simulate real flows inside combustors, Reynolds similarity criteria was respected. The maximum mean velocities for the flows inside the two rectangular ducts are varying from 1 to 11 m/s (Re H/2 = 8000 to 42000) and the static pressure difference between the secondary flow and the primary flow varies between 20 and 130 Pa. 2.3 Assumptions and simplifications The named secondary flow is supposed to represent the cold air coming from the compressor of a gas turbine engine. Both the primary flow, which represents the burned gases flow inside the combustion chamber, and the secondary flow, are at the same ambient temperature, in our case approximately 27 C. Therefore, the blowing ratio is reduced to velocity ratio, and we have: ρ 1 = ρ 2 = ρ jet = ρ air.

Heat Transfer VIII 511 Two-dimensionality is a concession to theory - it is rare in nature and is very difficult to obtain in an experiment. However, at the entrance of the test section we have measured with the LDA system, a two-dimensional uniform flow, which spreads over 65 % of each rectangular duct width. The development region upstream the test section was designed to be long enough to obtain a fully developed and turbulent duct flow characteristic, and thus the boundary layer thickness for both rectangular ducts is equal to the half height of the duct, δ = H/2. The turbulence intensity on the duct centreline has been measured to be around 5 %. To reduce the edge effects, the LDA measurements were conducted for three adjacent jets in the spanwise direction centred at the middle of the test section width. 2.4 Laser Doppler Anemometry system A 5 W Coherent Argon Ion two-component laser system is used for the LDA measurements. The system operates with two blue beam lights of wavelength 488 nm and two green beam lights of 514.5 nm wavelength and allows direct measurement of the two velocity components: u, in the x direction, and v on the y direction. The transmitter and the receiver are situated in a single probe head containing all the necessary optics to form the control probe volume and also to measure the back-scattered signal. This signal is detected by photo-multipliers by means of multimode optic fiber and is then directed as an electrical signal to the data acquisition system. These measurements are treated to correct for fast particles bias. As fast particles will pass the measurement control volume more rapidly than slow particles, there is a trend to measure more frequently fast particles, and the simple mean velocity will be higher than the true mean. The correction treatment accounts for this bias by weighting the measured velocity with the gate time of each particle passing through the probe volume. Table 1 shows the nominal optical characteristics of the LDA system used in our laboratory. Table 1: Nominal optical characteristics of the LDA system. Transmitter lens focal length (mm) 250 Probe volume length (mm) 1.5 Beam waist diameter (µm) green / blue 117 / 111 Fringe spacing (µm) green / blue 3.22 / 3.06 3 Velocity measurements and uncertainties The LDA system is equipped with 40 MHz frequency shifting process in order to distinguish between forward and reverse flow and to avoid fringe bias errors in measurements. The LDA probe volume is precisely positioned by a 3D computer controlled traverse system, and the positional errors are estimated to less than 1/60 mm. The measurement time at one point varied from about 10 s to 120 s, and the mean values of u and v velocity components were obtained from

512 Heat Transfer VIII minimum 2 10 valid particles. Uncertainties in the mean velocity can be estimated for up to 1% for u component and up to 2% for v, at the ducts centreline regions. Uncertainties in the turbulence can be estimated up to 5 percent for both u rms and v rms. Velocity profiles were measured with a very high density of measuring points in the near-wall region at a minimum wall distance of y + < 3. 4 Results and discussions As shown in Figure 3, LDA measurements were conducted in three principal zones: upstream the perforated zone, downstream the perforated zone and in the middle of the perforated zone. Figure 3: Test Section flows, dimensions and LDA Measurement Zones. Instantaneous velocities u and v, as well as turbulent Reynolds shear stresses are measured simultaneously in vertical planes, in the longitudinal and then normal directions to the flow. Since the primary flow is fully developed and has perfect longitudinal centreline symmetry, the experimental grid domain was reduced to the superior half of the primary rectangular duct. At the entrance of the test section, LDA measurements were conducted ten diameters, upstream the perforated zone, in order to verify the fully developed duct flows hypothesis. Figure 4 shows the mean velocity profiles and the turbulent Reynolds stresses normalised by the friction velocity against y +.

Heat Transfer VIII 513 Figure 4: Mean velocity profiles and normalised turbulent Reynolds stresses. These results are totally in agreement with the universal von Karman log-law velocity profile for an incompressible turbulent boundary layer near the plate surface (see, e.g. theoretical fundamentals in texts of Pope [4] and Rohsenow et al [5]). In the perforated zone, for the mean velocity profiles there is a high similarity starting from the seventh row of inclined holes. All velocity profiles measured at three hole diameters downstream hole centre axes, and for each even row, are characterised by two velocity picks, like in Figure 5. Figure 5: Mean velocity contours a), and velocity profiles for even rows b). The first pick represents the jet core and is obtained as a result of the interaction between the jet issued from the hole and the primary free stream flow. The normal coordinate of the maximum velocity represents the jet penetration. The second pick represents the film cooling core, which is obtained from the previous interactions between the upstream jets and the primary free stream flow. The film cooling thickness is defined as the normal coordinate where the local mean velocity starts to be greater than the unperturbed upstream mean velocity.

514 Heat Transfer VIII Downstream the perforated zone, after the twelve staggered rows of holes, the jets are totally merged and a film thickness of about six times hole diameters is obtained and maintained for more than 24 diameters in the downstream direction. Figure 6: Mean velocity contours in: a) normal and b) longitudinal to flow plan. Figure 6-a shows the well-known kidney shape of jet in cross flow for a normal plan (y, z) at 3d downstream the number eight row and Figure 6-b shows a mean velocity contour plot of the film, downstream the number twelve row. To compare different aerodynamic flow configurations and to study the influence of the blowing ratio on the film cooling characteristics, a normalisation by the duct centerline velocity of the primary flow was used. Figure 7-b, shows normalised velocity profiles for some flow configurations with the blowing ratio varying from M = 1.69 to 5.33. Figure 7: Film thickness evolution a), and downstream velocity profiles b). Different studies have related that the aerodynamic jet flows are very sensitive to the hole angle inclination and the hole exit shape as well as to the hole diameter and compound angle injection of the cold air, see Lee et al [6]. As

Heat Transfer VIII 515 one can see, in Figure 7-a, increasing blowing rate has little influence on jet penetration in the primary flow. Meanwhile, in the downstream zone of the perforated region, the film cooling aerodynamic characteristics parameters, such as the maximum velocity or the film thickness, are directly related to the blowing ratio. Also, flow configurations with high blowing ratios are dangerous for gas turbine engines because the cold air, which penetrates into the flame tube, could locally stop the combustion reactions and create unburned gas components. With the same cold air consumption, at small blowing ratios, better results for the film cooling efficiency can be obtained using different geometrical plate configurations. So, there is no need in using too much cold air coming from the compressor for the cooling process, but it is suitable rather to ameliorate the geometrical parameters of the multi-holed plates cooling configurations. These results are totally in agreement with experimental investigations of Emidio et al [7] conducted on multiholed flat plates used in today s combustor, at a 1:1 scale and in realistic aerodynamic and thermal environment conditions. 5 Conclusions Two components LDA measurements were conducted on a 10:1 perforated plate model in order to investigate the influence of aerodynamic flow parameters, like blowing ratio, flows pressure difference and Reynolds number, as an application to combustor walls cooling techniques. Except particular cases, the blowing ratio has little influence on the film cooling thickness downstream the perforated zone. That implies to choose for the multiholed plate flows, the minimum cold air fraction coming from compressor since a major part of the cold air is used for the combustion process, the gas dilution and the post-combustion process of the gas turbine engine. The present work reports on partial results of LDA measurements and the achieved database will be used to validate different physical models and numerical codes for multi-holed flow cooling techniques, and will serve also for future turbulent aerodynamic investigations in combustors. Nomenclature d circular diameter of the jet holes (mm) e perforated plate thickness (mm) H rectangular duct height (mm) M blowing ratio = (ρv) jet / (ρv) 1 P static pressure (Pa) r inter-row distance (mm) Re Reynolds number s spanwise holes distance (mm) u, v, w velocity components (m/s) V mean total velocity (m/s) x, y, z Cartesian coordinates x local row coordinate

516 Heat Transfer VIII Subscripts 1 primary flow 2 secondary flow air air flow atm atmosphere jet jet flow rms root mean square Symbols α hole angle inclination (degree) δ boundary layer thickness = H/2 (mm) pressure difference (Pa) ρ density (kg/m 3 ) ν cinematic viscosity (m 2 /s) References [1] Fric T.F., Roshko, A., Vortical structures in the wake of a transverse jet, J. Fluid Mech. 279, pp. 1-47, 1994. [2] Goldstein R. J., Effect of mainstream variables on jets issuing from a row of inclined round holes, J. Eng. Power 101, pp. 289-304, 1979 [3] Brown A., Saluja C.L., Film cooling from three rows of holes on adiabatic constant heat flux and isothermal surfaces in the presence of variable freestream velocity gradient and turbulence intensity, ASME 79, GT 24, 1979. [4] Pope B. Stephen, Turbulent Flows, Cornell Cambridge University Press [5] Rohsenow W. M., Hartnett J. P., Ganic E. N., Handbook of Heat Transfer Fundamentals, Second Edition, McGraw-Hill Book Company, 1985 [6] Lee H.W., Jung J. P., Lee J. S., Flow visualization and film cooling effectiveness measurements around shaped holes with compound angle orientations, Int. J. of Heat and Mass Transfer 45, pp. 145-156, 2002 [7] Emidio J. M., Leger B., Miron P., Geometric and aero-thermal influences on multiholed plate temperature: application on combustor wall, Int. J. of Heat and Mass Transfer 46, pp. 1215-1222, 2002