STEADY FLUID FLOW INVESTIGATION USING PIV IN A MULTI-PASS COOLANT CHANNEL

Transcription

1 11 th International Symposium on Applications of Laser Techniques to Fluid Mechanics, 8 11 July 2, Lisbon, Portugal, Paper No. 17- (COMPLEX WALL FLOWS) STEADY FLUID FLOW INVESTIGATION USING PIV IN A MULTI-PASS COOLANT CHANNEL Dr.-Ing. Marc P. JARIUS German Aerospace Center Institute of Propulsion Technology Linder Höhe, D-1147 Köln, Germany Fax Number: address: marc.jarius@dlr.de Dipl.-Ing. Martin ELFERT German Aerospace Center Institute of Propulsion Technology Linder Höhe, D-1147 Köln, Germany Fax Number: address: martin.elfert@dlr.de KEY WORDS: MULTI-PASS COOLANT CHANNEL, PARTICLE IMAGE VELOCIMETRY, LASER-2- FOCUS VELOCIMETRY. ABSTRACT In a stationary two-pass coolant channel system the fluid flow is investigated experimentally. One of the main objectives of this paper is to describe the capability of the measuring technique Particle-Image Velocimetry (PIV). Results of heat exchange and pressure losses are not the objectives of this paper. The investigated system has an enginenear lay-out with an 18 turn and has smooth walls in the beginning of this project, later on ribbed walls will be investigated. As a first step, the system is analyzed in non-rotating mode. During future work it will rotate about an axis orthogonal to the center-line of the straight passes. The results shown in this paper demonstrate the effect of the 18 bend with isothermal flow condition excluding any buoyancy. Turbulent channel flow with a REYNOLDS number of. and., derived with the hydraulic diameter of the first pass, was investigated. At the test rig at DLR several test models were mounted to investigate pressure drop behavior and fluid motion separately. Numerical investigations regarding flow structure and pressure drop have been carried out in order to validate a CFD code against the experimental data. The numerical solution is compared with streamline pictures, obtained from flow visualization of all walls using oil flow technique, and near wall PIV results. Mid-span flow distribution is received from laser light sheet visualization technique with oil fog as indicator and also from PIV. Figure 1 shows a comparison of the measured mean flow velocity distribution at mid-span (cut 1 & camera position 1) with and without turning vane. The results presented in this paper clarify the complex flow situation given by the two pass system with inherent turn. Especially in the bend region appear separation regions and vortices with high local turbulence. These very demanding measuring task represents a benchmark test case for the used measuring techniques z inmm z inmm Figure 1. Comparison of the mean flow velocity distribution at mid-span (cut 1 & camera position 1) with and without turning vane 1

2 1 INTRODUCTION The power generation industry and the aero-engine industry operate in a highly competitive market. Major objectives are high requirements in terms of high technology development as well as cost and development time reduction. On the other hand, environmental and safety constraints are an increasingly stringent necessity, which enforces the demand for new technologies. Currently the industry relies on expensive and time consuming rig test programs, whereas the existing numerical design tools have a number of deficiencies in accurately describing the complex multipass coolant channel flow. Under these conditions the Institute of Propulsion Technology is involved in national and European research programs aimed to providing the industry with high quality experimental data from the flow field for CFD validation. In past projects the flow behavior in rotating passages was analyzed at DLR. Laser-2-Focus velocimetry (L2F) was used to obtain flow velocity components and fluctuations. Although time-consuming, this non-intrusive, singlepoint measurement technique worked very well within straight and smooth duct flows, which generally have a moderate degree of turbulence. However, state of the art, serpentine shaped, multi-pass systems are equipped with ribbed walls, in order to improve heat exchange which is typical in realistic configurations. In this case L2F velocimetry was not able to measure accurate flow properties due to the increased turbulence intensities in the vicinity of the ribs as in the bend region and further downstream. The dividing wall separating the two passages forces the flow into a sharp turn which generally results in a flow separation. The flow within the separation bubble itself is very unsteady. An important work package within these research programs is to provide detailed information of the flow field inside the multi-pass coolant channel using planar techniques such as Particle Image Velocimetry (PIV) [ADRIAN, 1991, WILLERT, GHARIB, 1991, RAFFEL et al, 1998] or Planar Doppler Velocimetry (PDV) [RÖHLE, ]. Modern planar measurement techniques such as PIV are capable of obtaining complete maps of flows even at high turbulence. As a first step toward applying this technique, a multi-pass cooling system with and without turning vane is investigated in stationary (e.g. non-rotating) mode using two-component PIV. The results are compared with results from L2F (multipass cooling system with turning vane), flow visualization and CFD. The high quality of the obtained new results encourage the application of two-component PIV to the rotating system. As a logical consequence, the application of three-component PIV will be necessary to obtain the complete flow field information within the complex flow passages. This paper s intention is to report on the status of applicability of PIV in a multi-pass coolant channel. 2 NOMENCLATURE d h m Hydraulic diameter f mm Lens focal length f# - Lens f-number ν m²/s Kinematic viscosity n rpm Speed Re - REYNOLDS number, Re=ρv d/η τ µs Pulse delay v m/s Volumetric absolute velocity v m/s Intake velocity x, y, z mm CARTESIAN co-ordinates Abbreviations p.s. s.s. t.e. t.v. pressure side suction side trailing edge turning vane 2

3 3 TEST FACILITY AND INSTRUMENTATION 3.1. Test Facility A large-scale coolant multi-pass system rotating spanwise was investigated during the past experiments. Air is supplied in the rig through a rotary sealing assembly, a second one is used for venting after the air has passed through the test section. Both are mounted to the end of the double hollow shaft. A schematic of the test rig is given in Figure 2. More detailed information about the test facility are presented in [RATHJEN et al, 1999]. Reynolds number and rotation number can be adjusted appropriately to the models inserted. For the PIV investigations a non-rotating test rig is installed at a laboratory and shown in Figure Test Model Geometry Figure 4 shows a schematic draw of the coolant multi-pass-model with turning vane. The test section consists of a leading edge duct (first pass) with a trapezoidal cross-section extending radially outward, a 18 bend with a 9 turning vane and a second pass with a different trapezoidal cross-section extending radially inward. The ratio between the hydraulic diameter of the second pass and the hydraulic diameter of the first pass is 1.. The length of each pass is 9.97 d h, where d h denotes the hydraulic diameter of the first passage. In the bend region the flow is not only directed through the 18 u-bend but also perpendicular to that following the turbine airfoil shape. The test section is 12 d h long. The distance between the tip plate and the divider wall is about 1.87 d h. The divider wall has a thickness of.3 d h. The model has a small plenum chamber and an inlet and outlet length, both are 3 d h long. A schematic draw of a second cooling multi-pass-model without turning vane and with a modified second pass is shown in Figure. The length of the wall 1* is reduced to 1.3 d h, the angle is reduced to 8 and the wall 1* is now parallel to 3*. The length of each pass is increased to.6 d h and consequently the distance between the tip plate and the divider wall is reduced to 1.2 d h. A CAD-picture of the model geometry without turning vane is presented in Figure Measuring Techniques PIV Particle Image Velocimetry (Flow) The flow field is measured first in steady mode by means of PIV. This modern technique is a scientific tool for qualitative and quantitative investigations of flow fields. The principle of the PIV with laser, light sheet optic, laser light sheet and test section is shown in Figure 7 and the existing test set-up in the laboratory is shown in Figure 2. PIV is based on the principle to capture the image of the flow field two times (pulse delay in the range of µs). Tracer particles (aerosol) with an average size of.1 µm are added to the flow. These particles reflect in the laser light sheet, which is formed with a light sheet optic consisting of a system of one spherical and two cylindrical lenses. It generates a light sheet with a thickness of 1 mm and a divergent angle of 24. The laser light sheet is adjustable and visualizes the first pass as well as the second pass, due to a mirror. Illumination was provided by a standard, frequency-doubled, double-cavity Nd:YAG laser (NewWave, Gemini PIV) with a pulse energy of up to 1 mj per pulse at 32 nm. On the recording side, a thermo-electrically cooled, interline transfer CCD camera (PCO, 128 x 24 pixel resolution) with a f = mm, f # 2.8 lens (Nikon) was used. A bandpass filter with center frequency 32 nm and width nm (FWHM) placed in front of the lens rejected most of the unwanted radiation. The pulse delay between the laser pulses varied between τ = 6 µs (investigation of mean flow) and τ = 12 µs (investigation of secondary flow). The high resolution CCD 1 camera takes two pictures, depending on the pulsing laser. Each picture has a size of 2. MB. The CCD camera, perpendicularly positioned to the light sheet, is a so-called cross correlation camera having double-frame single-exposure evaluation. Now it is possible to calculate the value, direction and orientation of the absolute velocities with the Cross Correlation Function (CCF). pictures at each camera position are taken to calculate the mean value of the absolute velocity. For the secondary flow investigations (cut 1...6) into the first pass from the multi-pass coolant channel one camera positions is adjusted (Figure 8a) and for the investigations of the mean flow (cut 1 & 2) into the second pass of the multi-pass coolant channel five camera positions are necessary (Figure 8b). 1 Charge Coupled Device 3

4 OFV Oil Flow Visualization (Wall Flow) To analyze the flow phenomena, several kinds of flow visualization techniques have been applied. Using the oil flow visualization method, wall streamline patterns could be achieved within the duct. It is possible to show and localize the boundary layer separation, the reverse flow and vortex regions next to the surfaces. To some extent it is possible to draw conclusions from the two-dimensional method for application to the strongly three-dimensional character of the flow field. The duct wall has been previously painted with a mixture consisting of oil of appropriate viscosity, oil acid and Titanium dioxide (TiO 2 ) for contrast purpose. Then it has been exposed to the flow for a time period (approx. minutes) long enough to dry the oil paint. The duct is dividable in two halves and is opened after the test run for viewing. A detailed description of the range of application is given in [MALTBY, KEATING, 1962] L2F Laser-Two-Focus Velocimetry (Flow) For the non-intrusive measurement of flow velocities, angles and fluctuations inside the rotating duct the Laser-2- Focus (L2F) technique [BEVERSDORFF, HEIN, SCHODL, 1992] has been applied using a newly developed optical unit to direct the laser beams into the rotating frame of reference. The optical device mounted in front of the rotor enables the system to stationary non-triggered signal processing. Following the beam path from the laser to the probe volume, the laser light is emitted by an Ar+ ion laser of 4 Watt power output which is connected to the launching unit mounted to the test rig by an optical fiber for vibration decoupling purpose. The L2F method is based on time of flight measurements of co-flowing particles. Small oil droplets with a size less than. µm were added to the flow before entering the shaft. The flow-following behavior of these particles even at high accelerations has been successfully demonstrated and thereby a correlation between particle and flow velocity is assured [SCHODL, 1977] LLSV Laser-Light-Sheet Visualization (Flow) To obtain more information about the flow around the turning vane and the separation zone, a camera system has been installed to the Perspex duct observing inside the second pass with the mid plane lighted with a laser light sheet. In order to visualize the flow structures, seeding of the flow with oil steam is provided by an oil fog generator. With frame grabbing on a PC the video pictures were captured. The digitized pictures have been improved with regard to luminance and brightness. 4 RESULTS In the following, some exemplary experimental results are presented. All results are presented in the form of isolines and vectors. For a better understanding of the isolines, red areas represent high velocities, blue areas represent areas with lower velocities. The shown vectors all have a uniform size to represent only the orientation of the flow field. With regard to PIV processing, standard cross correlation algorithms were used with an interrogation window size of 32 x 32 pixel and a % overlap. The velocity of the intake flow into the first pass amounts to v = 1. m/s, and consequently the REYNOLDS number is about. for the PIV measurements. The comparison of the axial velocity component at cut 1 and 2 between L2F and PIV measurement in the second pass at the position z = 193 mm is presented in Figure 9. The axial velocity component at the z-position is extracted off the measured PIV results. For the case with t.v. the results shows relatively good agreement with the L2F results. The position of the jet and the wake is slightly different, which depends on deviations in model manufacturing. Measurements in the region of the separation bubble were not possible with the L2F system due to too high turbulence. For the case without t.v. the flow separation on the dividing wall is much larger then with t.v. and the value of the axial velocity component inside the bubble is even higher. By means of the thickness of the separation bubble, the open cross section is contracted, due to that fact the value of the axial velocity component is increased. At mid-span the axial velocity component is smaller and the separation region is nearly not remarkable. The measured flow field of the secondary flow within the first pass in six cuts (z = mm, step size mm) for a REYNOLDS number of. are shown in Figure. At a distance of 19 mm from the duct entry (Figure a, cut 1) a relatively big and strong vortex occurs (x = 16. mm, y = 13. mm). A smaller vortex exists at the channel too (x = 8. mm, y =. mm). An interpretation of that phenomenon is, that the flow is disturbed due to the inlet distortions, which are still visible at this position. The same result was found during CFD simulation where the complete inlet geometry with settling chamber was modeled (Figure 1) at cut 1. The mean value of the secondary flow velocity is in the range of about 1 m/s. At the position of cut 2 (Figure b, z = mm) the big vortex 4

5 moved closer to the side wall (x = 17. mm, y = 11. mm) and the smaller vortex is already disappearing. The mean value of the secondary flow velocity is in the same range like before. If the position of the light sheet is varied to the positions of cut 3, 4, and 6 (Figure c, d, e and f) all vortices are disappeared. The mean value of the secondary flow velocity is increasing more and more from cut to cut and is finally reaching a value of 13 m/s at cut 6 (Figure f) and is becoming now the main stream velocity component due to the flow turning process. At cut 6 at the leading edge the beginning of the recirculation from the upper passage corner can be observed. Within the second pass with turning vane the measured flow field of the main flow in two cuts 2 mm (right, wall section) and mm (left, mid-span) for a REYNOLDS number of. are shown in Figure 11. At the mid-span section (left) a relative stationary separation bubble occurs at the upper surface (s.s.) of the t.v.. Behind the separation bubble exists a relative high flow deviation with high velocities about 24 m/s close to the wall. An other separation bubble with reverse flow is visible on the dividing wall. Due to that flow recirculation a jet flow with velocities about 24 m/s exists on the pressure side at the t.e. of the t.v.. The thickness of the separation bubble is decreasing in downstream direction. The form and extension of the separation bubble with a length of approximately d h is similar to those found by visualization (Figure 16a). The areas of high velocity are coinciding further downstream (z = 19 mm) where mixing of both jets leads to an uniform velocity profile. The mean value of the main flow velocity further downstream is in the range of about m/s. At the wall section (right) the mean value of the velocity is increased to m/s on the upper side of the t.v.. Turning of the flow leads to an impingement effect to the leading wall of the second leg due to the inclination of the center-lines of both cross-sections. The separation bubble on the dividing wall is much larger in longitudinal extension than before. At the s.s. of the t.v. appears a boundary layer separation (x = mm, z = mm) due to side wall effects. The jet flow at the pressure side of the t.v. is still existing. But the mean value of the main flow velocity distribution is decelerated further downstream to a range of about 8... m/s. The measured flow field of the main flow within the second pass without turning vane in two cuts 2 mm (right, wall section) and mm (left, mid-span) for a REYNOLDS number of. are shown in Figure 12. At the mid-span section (left) a separation bubble with reverse flow is visible on the dividing wall. The thickness of the separation bubble is decreasing more and more in downstream direction. The form and extension of the separation bubble with a length of approximately 2 d h is similar to those found by visualization (Figure 16b). At the top of the duct exists a relative high flow deviation with the highest velocities with a value of over 24 m/s. This area of high velocity (> 24 m/s) reaches up to z = 1 mm downstream. At the wall section (right) the area with high velocity over 24 m/s is increased and the extension of the separation zone at the dividing wall is larger as at the other cases. The streamline patterns on the pressure side (Figure 13) and the suction side (Figure 14) obtained from oil flow visualization experiments for both cases, with and w/o turning vane are compared with a numerical solution inclusive t.v. close to the wall. The streamlines on both walls show the combined effect of the 18 turn, the influence of the turning vane and clearly the extension and thickness of the separation bubble at the top of the dividing wall in the second channel. The main phenomena in this engine-near configuration are: (i) The separation bubble at the dividing wall is larger at the pressure side than at the opposite wall. (ii) The deflection of the flow within the bend due to an inclination of the orientation axis of the two passes forms a strong secondary motion after the bend running from the concave wall (p.s.) around the rear wall to the s.s. of the 2 nd pass. (iii) The bend induces a counter-rotating pair of secondary vortices, These are affected by the secondary motion described above. The vortex along the p.s. has an adverse orientation and is strongly weakened whereas the other along the s.s. is enlarged recognizable at the slope of streamlines against the pass axis. The saddle lines between both vortices were expected normally at the center of the side walls in a rectangular U-turn. In this case of an inclined-oriented two-pass system the saddle lines were shifted to the p.s. and to the corner of the rear side, respectively. After the separation bubble at the top of the divider wall, almost cross-oriented streamlines are visible over a curtain area at the pressure side confirming the dominance of the coflowing secondary vortices of (ii) and (iii). Both vortex generation mechanisms depend on inertia effects in viscous flows. A comparison of a calculated result (Figure 1a) and measured results with PIV with t.v. (Figure 1b) and without t.v. (Figure 1c) is presented. The relatively big vortex close to the dividing wall is visible in all results, but the position of this vortex is a little different for the case without t.v.. The core of the vortex is centred more into the middle of the channel and the orientation is changed.

6 Two captures of the flow around the bend with and without t.v. obtained by laser-light-sheet visualization is shown in Figure 16a & b. This technique confirms the already recognized and already described flow phenomena. SUMMARY This paper describes experimental investigations of a typical two-pass cooling channel system of a turbine blade with respect to fluid flow phenomena affected by geometry effects and flow turning. A complex flow situation is present showing several kind of separation and vortices: (i) At the upper surface of the t.v. due to a to high angle of attack versus the incoming flow, (ii) further downstream at the s.s. of the t.v. due to too strong curvature, (iii) at the dividing wall due to the sharp edge where the flow is incapable to follow. The jet region under the t.v. and the wake from the t.v. are leading to high velocity gradients and corresponding high shear stresses at the contact layer producing very high turbulence levels there. Higher turbulence levels than 3 % leads to uncertainties in L2F results. With supposed levels of % and more the PIV technique has no difficulties. But in order to give an accurate value of the turbulence level by PIV, you will need about pictures or more. This relatively big amount of data demand high storage capacities. The performed flow visualization and simulation are very helpful in that case to understand the interaction of all apparent effects. The presented experiments were mainly intended to assess the feasibility of applying PIV in a multipass cooling channel used for applied cooling research. The results have shown that it is possible to apply PIV in a multi-pass cooling channel. The application of PIV is approved and will be adopted to the rotating system in the near future. As a logical consequence, the application of three-component PIV will be necessary to obtain the complete flow field information within the complex flow passages. 6 REFERENCES ADRIAN, R. J., 1991, Particle-Imaging Techniques for Experimental Fluid Mechanics, Annual Review of Fluid Mechanics 23. BEVERSDORFF, M., HEIN, O., SCHODL, R., 1992, An L2F-Measurement Device with Image Rotator Prism for Flow Velocity Analysis in Rotating Coolant Channels. 8 th Symp. On Heat Transfer and Cooling in Gas Turbines, AGARD PEP, Antalya, Turkey. ELFERT, M., HOEVEL, H., TOWFIGHI, K., 1996, The Influence of Rotation and Buoyancy on Radially Inward and Outward Directed Flow in a Rotating Circular Coolant Channel, Proc. th ICAS Conference, Sorrento, Italy, 249- MALTBY, R. L., KEATING, R. F. A., 1962, The Surface Oil Flow Technique for Use in Low Speed Wind Tunnels, AGARDograph 7. RATHJEN, L., HENNECKE, D.K., ELFERT, M., BOCK, S., HENRICH, E., 1999, Investigation of Fluid Flow, Heat Transfer and Pressure loss in a Rotating Multi-Pass Coolant Channel with an Engine Near Geometry, ISABE- Conference, Florence, ISABE- IS-216 RAFFEL, M., WILLERT, C., KOMPENHANS, J., 1998, Particle Image Velocimetry, Springer Verlag Berlin Heidelberg. RÖHLE I., Doppler global velocimetry, in RTO Lecture Series 217, Planar Optical Measurement Methods for Gas Turbine Components, Neuilly-Sur-Seine Cedex, France. SCHODL, R., 1977, Entwicklung des Laser-Zwei-Focus-Verfahrens für die berührungslose Messung der Strömungsvektoren, Dissertation, Tech. Univ. Aachen, Germany. WILLERT, C., GHARIB, M., 1991, Digital Particle Image Velocimetry, Experiments in Fluids, No., Springer Verlag. 6

7 7 Figures Light sheet optic Laser Acrylic model CCD-camera Flow meter Figure 2. Schematics of the test rig at DLR-Cologne Figure 3. Test set-up for stationary PIV-measurements 12dh 1.87dh 9.97dh.1dh 1.1dh.12dh x/dh = dh Leading wall 4 1. dh 12 dh.6 dh 1.2 dh.1 dh dh dh dh 1 1.dh axis of rotation 4 3 axis of rotation dh = hydraulic diameter of the first pass dh = hydraulic diameter of the first pass 3dh x/dh = Inlet.3dh Outlet x/dh = 22.9 Trailing wall 3 dh x/dh = Inlet.3dh Outlet x/dh = 22.9 Figure 4. Schematic of test model geometry with turning vane Figure. Schematic of test model geometry without turning vane 7

8 Light sheet optic Laser Mirror Laser light sheet Tracer particle into light sheet Particle image from first light pulse at t Imaging optic Image plane Flow direction Particle image from second light pulse at t 1 Figure 6. Test model geometry without turning vane Figure 7. Principle of Particle-Image Velocimetry 8

9 a b Figure 8. a Camera and cut positions. a. Secondary flow investigations into first pass (cut 1...6); b. Mean flow investigation into second pass (cut 1 & 2) b 3 3 Wallsection, PIV without tv Wall section, PIV with tv Wall section, L2F Wall section, L2F, corr. L2F vinm/s 1 L2F PIV vinm/s 1 PIV Region of Backflow 1 3 Mid-span, PIV without tv Mid-span, PIV with tv Mid-span, L2F Mid-span, L2F, corr. 1 3 Figure 9. Comparison of the axial velocity component between L2F and PIV measurement in the second pass at the position z = 193 mm. a cut 1; b cut 2; (according Figure 8b) 9

10 a `ìí=ni=ò=z=nvm=ãã 3 Cut 1, z = 19 mm b `ìí=oi=ò=z=omm=ãã 3 Cut 2, z = mm y 1 y 1 c 1 3 `ìí=pi=ò=z=onm=ãã Cut 3, z = 2 mm 1 d 3 `ìí=qi=ò=z=oom=ãã Cut 4, z = 2 mm y 1 y 1 e 1 3 `ìí=ri=ò=z=opm=ãã Cut, z = 23 mm f 1 3 `ìí=si=ò=z=oqm=ãã Cut 6, z = 24 mm y 1 y Figure. Secondary flow velocity distribution in six cross cuts of the first pass measured with PIV (Re =.)

11 Mid-span y = mm 24 vin m/s: Pos 1 24 vin m/s: Wall section y = 2 mm z in mm 23 2 z in mm Pos Pos Pos Pos Figure 11. Mean flow velocity distribution in two longitudinal cuts of the second pass with turning vane measured with PIV (Re =.) 11

12 Mid-span y = mm v in m/s: Pos 1 v in m/s: Wall section y = 2 mm Pos 2 Pos 3 Pos 4 Pos Figure 12. Mean flow velocity distribution in two longitudinal cuts of the second pass without turning vane measured with PIV (Re =.) 12

13 Suction side Figure 13. Comparison between calculated and visualized wall velocities and streamlines on suction side, nonrotating (Re =.) 13

14 Pressure side Figure 14. Comparison between calculated and visualized wall velocities and streamlines on pressure side, nonrotating (Re =.) 14

15 a b c Figure 1. Comparison of calculated and measured results. A Numerical simulation at cut 1 in first pass; b Measured at cut 1 in first pass with turning vane; c Measured at cut 1 in first pass without turning vane a b Figure 16. Pattern of laser light sheet visualization at the bend region. a With turning vane; b Without turning vane 1