6th International Advanced Technologies Symposium (IATS’11), 16-18 May 2011, Elazığ, Turkey
Experimental Investigation of the Flow Structure around a Sphere and Its Control with Jet Flow via PIV M. Ozgoren1, A. Okbaz2, A. Kahraman3, R. Hassanzadeh4, B. Sahin5, H. Akıllı6 and S. Dogan7 1,2,3,7
Selcuk University, Konya, Turkey,
[email protected],
[email protected],
[email protected],
[email protected] 4,5,6 Cukurova University, Adana, Turkey,
[email protected],
[email protected],
[email protected] structure around a smooth and vented sphere using PIV technique.
Abstract—The flow characteristics of a smooth sphere and a vent sphere located in a uniform flow for the Reynolds number range 2500≤Re≤10,000 are studied quantitatively using the particle image velocimetry technique and qualitatively with dye visualization. The flow phenomena in the downstream regions of the sphere increase the instability of the vortical flow structure significantly. The present investigation concentrates on the characterization of the flow structures using the time-averaged and instantaneous flow data in a plan view at laser sheet location passing through the sphere equator. Distributions of velocity vectors, patterns of streamlines and vorticity contours characterize the flow structure, in detail. It is observed that the modified flow structure of the near wake of the vented sphere may be characterized by a pair of counter-rotating ring vortices, which have the effect of aerodynamically streamlining the sphere. The obtained results can be helpful for developing and validating numerical predictions as well as designing. Keywords—Flow control, , Instability, PIV, Sphere, Vorticity, Turbulence, Vortex shedding, Wake flow
I. INTRODUCTION
F
low structures in the downstream of a single sphere are encountered so frequently in various industrial applications such as various multiphase flow system, air pollution, nuclear and thermal power plants, pneumatic and hydraulic conveying, chemical and food processing, combustion systems, sport balls and bomb. Large amount of researches have been conducted flow around a sphere for different Reynolds number and massive amounts of data for vortex shedding frequency, pressure coefficient, drag and lift force coefficients have been accumulated as cited in the references and cited therein. However, at moderate Reynolds numbers such as 2,500≤Re≤10,000 due to lack of data under various conditions, more work is required to fully understand the effects of the smooth and vented sphere. Some of studies on flow characteristics and control around sphere and other bluff bodies were given and cited therein 125. None of these studies have not been compared the flow
II. EXPERIMENTAL SETUP Experiments were performed in a large-scale open water channel with a test section length of 8000 mm and a width of 1000 mm at the Department of Mechanical Engineering at Cukurova University, Turkey. To perform the present experimental study, the test section made from 15 mm thick transparent Plexiglas sheet, which had a total height of 750 mm, was filled with water to a level of only 450 mm. Before reaching the test chamber, the water was pumped into a settling chamber and passed through a honeycomb section and a two-to-one channel contraction. An overview of experimental system of the sphere is shown in Fig. 1. Free stream turbulence intensity of the flow is less than 0.5% in the range of the present Reynolds numbers, Re U D v , based on the sphere diameter (D). Here, and D are kinematics viscosity and diameter of the sphere (D), respectively. U∞ is free stream velocity taken as 117.6 mm/s. The sphere with a diameter of 42.5 mm was made of Plexiglas so that the laser light propagates easily from them. In addition, water cell segment of the sphere equator with a diameter of 38.5mm and a wall thickness of 2.0 mm wascreated. It was filled with distilled water and had a total height of 8 mm in order to reduce largely the laser light deflection on the sphere. The sphere surface was highly polished to avoid effects of surface roughness. The laser sheet was located at 225 mm above the bottom surface of the channel while the water height hw was 450 mm in all cases. To support the sphere in the water channel, a circular bar with a 5 mm diameter was connected to the sphere from the top. Disturbance effect of the support bar on the laser sheet location of the measurement plane that was observed by dye injection was negligible in the consideration of support diameter with respect to the sphere diameter. The solid blockage ratio of the sphere including support was 1.3 %. Sphere models are presented in Figure 2b. The vented sphere made of solid plexiglas and did not permit to pass the laser light. The ratio of the smooth sphere diameter
272
Experimental Investigation of the Flow Structure around a Sphere and Its Control with Jet Flow via PIV (D=42.5 mm) to vented hole diameter (d hole =6.4 mm) is dhole /D=0.15, which provides a vent area of 2.26%.
(a)
An interrogation window of 32x32 pixels in the image was selected and converted to grid size approximately 1.44x1.44 mm2 for the single sphere (0.034Dx0.034D). The overall fields of physical view were for both spheres, yielding to 7,227 (99x73) velocity vectors for whole taken images. During the interrogation process, an overlap of 50% was employed in order to satisfy the Nyquist criterion. Patterns of instantaneous particle images with a total of 350 images for a continuous series were taken at the rate of 15 Hz, thereby spanning 23.27 sec. Averaged patterns of the flow structure were calculated from all of the instantaneous images. The laser sheet was generated from a dual pulsed Nd:YAG system, having the maximum output of 120 mJ per pulse, which had time delays t =1.0-1.7 ms for the present experiments. Inappropriate displacement vectors caused by shadows, reflections, or laser sheet distortions in the flow field replaced by using bilinear interpolation between surrounding vectors in the postprocessing step. This algorithm included magnification factor and image captured rate to calculate velocities from the valid vectors. The field was then smoothed by a Gaussian weighted averaging technique. To minimize distortion of the velocity field, a smoothing parameter of 1.3 was chosen. After having been vector field, the vorticity patterns of the wake flow were determined from the velocity field using a finite difference scheme with an in-house software. The details of the effects of these factors can be found in the studies of Hart 26, Fouras and 27, Westerweel 28 and Adrian 29. The factors contributing to uncertainty in the velocity measurement using the PIV technique were critically assessed by Westerweel 28 concluding that uncertainty estimation in the velocity measurement was less than ±2%. III. RESULTS AND DISCUSSION
(b)
Figure 1. Schematic view of the experimental set up and sphere models
Nd:YAG laser was used to generate a laser sheet that was perpendicular to the axis for the sphere and the symmetry axis (i.e. equator of the sphere) was passed through them. A CCD camera having a resolution of 1,600 x 1,186 pixels was used to record the images. The seeding particles with a diameter of 10 µm in the flow were silver metallic coated hollow plastic spheres. The densities of the particles and water are close enough so that the distribution of particles in suspension remains uniform for several hours. The high-image-density criterion was satisfied by ensuring that a minimum of approximately 20-30 particles was contained within the interrogation area. The illuminating laser sheet thickness in the flow field was approximately 1.5 mm. As shown in Fig. 1, the camera was mounted in a fixed position beneath the water tank. Dantec Flow Grabber digital PIV software employing the cross-correlation algorithm was used to compute the raw displacement vector field from the particle image data.
Figures 2a and 2b show typical visualization images of instantaneous flow fields around a smooth sphere and a vented sphere for 5,000 in which dye ports with 0.7mm diameter are located on equator of the sphere at angle values with respect to the flow direction as 0o, 70o,90o, 110o, 180o, 290o, 270o and 250o. Dye visualization representative images are is presented with two different images to show evolution and progress of the small scale vortices designated by A to H. The separated and recirculating flow in the near-wake region with the help of visualizing are clearly seen for the sphere with laser illumination using the Rhodamine dye injection technique in the near wake region. Small scale vortices around the wake region are formed around larger vortices with a wavy appearance due to Kelvin Helmholtz instability. Formation of the spiral vortices begins to occur in the very close region of the sphere. As the flow travels in the downstream direction, the dimensions of the vortices increase around the bluff body. Then these vortices are shed from the periphery of the sphere directly to the inward wake region. The large eddies are formed at a regular frequency, and they produce pressure disturbances in the flow. The flow patterns in Figs. 2a and b show that the laminar boundary layer separates at around =85o±5o for the present Reynolds number range, where is measured from the front stagnation point. Shedding shear layer
273
M.Ozgoren, A.Okbaz, A.Kahraman, R.Hassanzadeh, B.Sahin, H.Akıllı, S.Dogan becomes unstable due to the Kelvin-Helmholtz instability caused by the large velocity difference at the interface between the free-stream and sphere wake flow regions. Thereafter, the laminar shear layer turns into a powerful turbulent flow structure. Several vortex-ring shaped protrusions appear as an indication of the shear-layer instability along the borders between the wake and free stream regions, as also observed by Jang and Lee 5.
(a)
(b) Figure 2. Comparison of flow visualization of flow structure with laser illumination of Rhodamine dye injection technique around the smooth sphere (a) and vented sphere (b) at 5,000.
It is seen that the near wake circulating region is large and the wake generates a progressive wave motion. Weaker and small-scale vortices occur around the primary vortices. In the case of the vented sphere, the energy level of the flow decreases and then free shear layer emanating from the separated flow region around the sphere has a less tendency to travel in the near wake region of the sphere than smooth sphere case. This movement generates vortex induced vibration and provides a loading effect on the sphere body. On occasion, the shedding vortices have taken place at both side of the central axis of sphere symmetrically and simultaneously. These small size vortices rotate about their axis and move on
in a wavy form in the wake region for both spheres as seen in Fig 2. Shear layers emanating from the both side of the sphere merge at a location approximately two sphere diameter (2D) from the central point of the sphere and instability rises after this point as seen in Fig 2. The vortices generated by the shearlayer instability travel downstream and eventually compose a large-scale waviness of vortical structures in the wake. These instabilities and disorderness in the flow retain further downstream in the free stream flow direction having an “S” form like von Karman vortex streets approximately after a distance of 2.5D from the central point of the sphere. Figures 3 and 4 show comparison of instantaneous velocity fields V and * around the smooth sphere (left column) and vented sphere (right column) for 2500≤Re≤10,000. Instantaneous vorticity contours of the wake structure is normalized as * (i.e., *=D/U∞). All figure dimensions are normalized with the sphere diameter designated as x/D and y/D. Even though small-scale Kelvin Helmholtz’s vortices are clearly shown for dye visualization results, PIV results display more rounded vortices. However, numerous eddy vortices occur due to the three dimensional and complex flow structures. Jet flow through the ventilation hole make the flow structure more symmetric and increase the fluctuations in the base region of the sphere. The vortices produced from the flow separation around the periphery of the sphere have a tendency to move inwards because of the lower pressures prevailing within the wake. This situation is counter-balanced by the growing wake size, which shifted the vortex centerline outwards. Regarding the onset and development of small-scale vortical structures in the separating shear layer, regions of lowlevel vorticity concentration are discernible in the pattern of instantaneous vorticity for all patterns. The wake region accommodates velocity vectors with very small magnitude in the downstream region of the smooth and vented sphere, which are the source of small-scale secondary vortices, as seen in Fig. 3. The streamwise separation of successive vorticity peaks in the near wake region for the smooth sphere is larger than that for the vented one. The flow is three-dimensional, and shedding vortices convey fresh fluid into the wake flow region, magnifying the entrainment processes, which occurs with higher order of magnitude for the vented sphere. However, it is clear that there is only a small change on the flow patterns in the presence of ventilation. On the other hand, Suryanarayana and Prabhu (2000) stated that the observed drag reduction of wake unsteadiness, presumably caused by the stabilization/weakening of the randomly rotating 3-D vortical structure of the smooth sphere, which was also observed in the present study. As seen in Fig 4, secondary vortices forms in the wake region due to jet flow through the vent. It may be noted that any reduction in the wake unsteadiness may be expected to result in drag reduction. The modified flow structure of the near wake of the vented sphere may be characterized by a pair of counter-rotating ring vortices, which have the effect of aerodynamically streamlining the sphere as expressed by Suryanarayana and Prabhu (2000).
274
Experimental Investigation of the Flow Structure around a Sphere and Its Control with Jet Flow via PIV be concluded that the flow around spheres is three dimensional. The time-averaged streamline pattern shapes extended upstream to near base region of the body are elliptical and have a limited spiral cycle. Comparison of the time-averaged patterns shows that flow structures of the wake are almost equally symmetrical with respect to the centerline of the smooth and vented sphere. As the Reynolds number increases, saddle point S1 becomes closer to the base of the sphere due to the increased momentum transfer. For smooth sphere case, two foci F1 and F2 occur while flow structure of the vented sphere case exerts four foci as F1, F2, F3 and F4 due to the jet flow through the ventilation hole.
Figure 3. Comparison of instantaneous velocity fields V around the smooth sphere (left column) and vented sphere (right column) for 2,500≤Re≤10,000.
In the present study, the Reynolds number range is limited due to restriction of free stream flow velocity in the water channel and system capability. Therefore, all experiments were performed at the subcritical Re. As pointed by Suryanarayana and Prabhu (2000), the effect of natural ventilation can be categorized in two broad regimes as weak and strong interaction regimes. At subcritical Reynolds number (Re4x105), the weak interaction regime occurs and the broad features of the basic unvented sphere are largely unaltered despite the larger addition of mass in the near wake. Strong interaction is promoted by the closer proximity of the inner and outer shear layers at the supercritical Re, which results in a modified/weakened and steady near-wake flow characterized by reduced unsteadiness. Comparison of time-averaged streamline topology around the smooth sphere (left column) and vented sphere (right column) for 2500≤Re≤10,000 is displayed in Fig 5. The junction point of the streamlines separated from the surface of the sphere is denoted with “S”. This junction point is called as Saddle point S which shows the merging point of the shear layers emanating from both sides of sphere. Locations of foci designated with “F1 and F2” that are rolling in the clockwise and in counter clockwise direction on with respect to the sphere central axis exhibit well-defined critical points. It can
Figure 4. Comparison of instantaneous vorticity contours * around the smooth sphere (left column) and vented sphere (right column) for 2,500≤Re≤10,000.
The effective length of the jet flow decreases with increasing the Reynolds number. Two saddle points S1 and S2 for the vented sphere case occur owing to the presence of jet flow. The saddle point on the time-averaged streamline topology generally forms later for the vented sphere for all Reynolds numbers. For example, at Re=5,000, the nondimensional wake lengths (L/D) measured from the base of the sphere to the saddle point S1 in the averaged streamline topology are approximately 1.15D and 1.22D for the smooth sphere and vented sphere, respectively.
275
M.Ozgoren, A.Okbaz, A.Kahraman, R.Hassanzadeh, B.Sahin, H.Akıllı, S.Dogan
Figure 5. Comparison of time-averaged streamline topology around the smooth sphere (left column) and vented sphere (right column) for 2,500≤Re≤10,000.
Figure 6. Comparison of Reynolds stress contour correlations u' v' U 2 around the smooth sphere (left column) and vented sphere (right column) for 2,500≤Re≤10,000.
Variation of Reynolds stress correlations u' v' U 2 for the smooth sphere and vented sphere is displayed in Fig 6. Welldefined Reynolds stress patterns due to fluctuations along the shear layers produce a weaker Reynolds stress region very close to the base of the vented sphere case, which occurs as a result of the entrainment of free-stream flow into this wake region. This behavior of the flow provides the momentum transfer from the free-stream flow into the wake region due to pressure difference between the wake and free-stream flow. The upstream extension of each large-scale cluster occurs further upstream, and the tip of the large-scale clusters is located closer to the base of the sphere for all Reynolds numbers as shown in Fig. 6. For the vented sphere case, peak values of Reynolds stress correlations are not varied considerably due to less effect of the jet flow in the subcritical Reynolds number range. Therefore, surface modifications as in examined in the present study may not always results effectively in all cases because there could be a need to retain the structural geometry due to design considerations as well as application of Reynolds number range. The vented sphere model could significantly suppress vortex formation by diffusing the concentrated vorticity in the shear layers behind the body for supercritical Reynolds numbers.
IV. CONCLUSIONS The flow structure in the downstream region of a smooth sphere and vented sphere for 2500≤Re≤10,000 was investigated using a PIV and dye techniques in a circulating open water channel. The following results are derived: Small and large-scale waviness of vortical structures occur in the wake region of the sphere in the streamwise direction. The formation of these small and large-scale waviness of vortical structures is assessed to be closely associated with the temporal evolution of vortices generated by the shear-layer instability. It is observed that the shedding location of the large-scale vortices rotates slowly and irregularly, and they rotate at random about the streamwise axis when they travel downstream for both smooth and vented spheres. It is observed that the wake is unstable with vortex shedding. The mean velocity field showed two peculiar large-scale recirculation vortices downstream of the sphere model, and the length of the circulation region was approximately one diameter in the wake region for the spheres. The PIV results for the instantaneous velocity fields of the smooth and vented spheres show an unsteady vortex formation mechanism. The PIV results of time-averaged turbulent
276
Experimental Investigation of the Flow Structure around a Sphere and Its Control with Jet Flow via PIV structures are consistent with the visualized flow field showing the onset of shear layer instability. Time-averaged and instantaneous flow data reveal that the size of the structure of the turbulent flow, the size of wake flow region, the location of singular and double points, the peak values of turbulence quantities as a function of Reynolds numbers. As observed by dye visualizations, flow structures from the vented sphere exhibit more parallel vortex streets that keeps their form without combination in downstream direction. The shear layer surrounding the recirculation bubble behind the sphere has a region of intensified velocity fluctuations with higher values of the Reynolds stress correlations than the vented sphere because of the three-dimensional vortex interactions. However, the concentration of small scale vortices (eddies) is more dominant in the wake of the vented sphere than that of the single sphere. The instantaneous vorticity fields reveal an unsteady wavy structure of the sphere wake. The vortical flow structure in the near-wake region of the sphere show an unsteady behavior and the onset of shear-layer instability. The appearance of the unsteadiness in the sphere wake could be the result of instability of the two streamwise counter-rotating vortices. Because of the jet flow through the ventilation hole at the sphere center, the wake region of the sphere becomes more symmetrical. Mass and momentum transfer due to the jet flow provides fresh flow in the wake region of the sphere. Because of variation of flow structure with time, the flow structure downstream of the sphere has an unsteady motion. The modified flow structure of the near wake of the vented sphere can be characterized by a pair of counter rotating ring vortices, which have the effect of aerodynamically streamlining the sphere. Experimental data obtained in this study will enhance our understanding of the smooth and vented sphere wake characteristics to develop a strategy for the validation of numerical models and design in the future.
[7]
[8]
[9] [10]
[11] [12]
[13] [14]
[15]
[16]
[17]
[18]
[19]
[20] [21] [22] [23]
ACKNOWLEDGMENT The authors would like to acknowledge the funding of the Scientific and Technological Research Council of Turkey (TÜBİTAK) under contract no:109R028, Coordinatorship of Selcuk University’s Scientific Research Project No.2004/131 and Scientific Research Projects Office of Cukurova University.
[24]
[25]
[26] [27]
REFERENCES [1] [2] [3]
[4]
[5] [6]
E. Achebach, “Experiments on the flow past spheres at very high Reynolds numbers,” J. Fluid Mech. 54, 565–575, 1972. E. Achebach, “Vortex shedding from spheres,” J. Fluid Mech. 62(2), 209–221, 1974. G. S. Constantinescu, K. D. Squires, “LES and DES investigation of turbulent flow over a sphere at Re=10,000,” Flow Turbul. Combust., vol.70, 267–298, 2003. I. Hadzic, V. Bakic, M. Peric, V. Sajn, F. Kosel, “Experimental and numerical studies of flow around sphere at sub-critical Reynolds number, ” Eng .Turbul. Model Exp. 5, 667–676, 2002. Y. II, J. Jang, S. J. Lee, “PIV analysis of near-wake behind a sphere at a subcritical Reynolds number,” Exp. Fluids. 44, Issue 6, 905-914, 2008. T. A. Johnson, V. C. Patel, “Flow past a sphere up to a Reynolds number of 300,” J. Fluid Mech. 378, 19–70, 1990.
[28] [29]
277
H. J. Kim, P. A. Durbin, “Observations of the frequencies in a sphere wake and of drag increase by acoustic excitation,” Phys. Fluids 31, issue.11, 3260–3265, 1988. M. Kiya, H. Ishikawa, H. Sakamoto, “Near-wake instabilities and vortex structures of three-dimensional bluff bodies, ” A review. J. Wind Eng. Ind. Aerod. 89, 1219–1232, 2001. A. Leder, D. Geropp, “The unsteady flow structure in the wake of the sphere,” SPIE 2052, 119–126, 1993. T. Leweke, M. Provansal, D. Ormie`res, R. Lebescond, “Vortex dynamics in the wake of a sphere,” Phys. Fluid, vol.11, Issue9, 12, 1999. I. Nakamura, “Steady wake behind a sphere,” Phys. Fluids 19, Issue.1, 5–8, 1976. M. Ozgoren, B. Sahin, H. Akilli, C. Karakus, S. Yayla, A. Ozturk, A. Akar, “Investigation of flow structure around a sphere via PIV,” I. Ankara National Aerospace Conference, UHUK 06-028, METU Ankara, 21-23 Sept., Turkey, 2006. H. Sakamoto, H. Haniu, “A study on vortex shedding from spheres in a uniform flow,” J. Fluids Eng. 112, 386–392, 1990. G. K. Suryanarayana, A. Prabhu, “Effect of natural ventilation on the boundary layer separation and near-wake vortex shedding characteristics of a sphere, ” Exp. Fluids. 29, Issue.7, 582–591, 2000. S. Taneda, “Experimental investigation of the wake behind a sphere at low Reynolds number,” J. Phys. Soc. Japan. 11, Issue.10, p.1104–1108, 1956. J. S. Wu, G. M. Faeth, “Sphere wakes in still surroundings at intermediate Reynolds numbers,” AIAA J., vol.3, Issue.8, 1448–1455, 1993. Y. Tsuji, N. Kato and T. Tanaka, “Experiments on the unsteady Drag and wake of a sphere at high Reynolds numbers,” Int. J. Multiphase Flow, vol.17, No.3, pp.343-354, 1991. Kim D. Yung, H. Choi, “Vortical structures behind a sphere at subcritical Reynolds numbers,” Phys Fluids 18(1):015102, doi:10.1063/1.2166454, 2006. D. S. Dandy, H. A. Dwyer, “A sphere in shear flow at finite Reynolds number: effect of shear on particle lift, drag, and heat transfer,” J. Fluid Mech 216, 381-410, 1990. T. Igarashi, “Flow characteristics of a circular cylinder with a vent,” Bulletin of the JSME 21, 154, 656-664, 1978. G. K. Suryanarayana, GEA, Meier “Effect of ventilation on the flow field around a sphere,” Exp Fluids 19, 78-88, 1995. G. K. Suryanarayana, H. Pauer, GEA Meier, “Bluff-body drag reduction by passive ventilation,” Exp Fluids 19, 73-81, 1993. M. M. Zdravkovich, “Review and classification of various aerodynamic and hydrodynamic means for suppressing vortex shedding,” J. Wind Eng Ind. Aerodyn, 7, 145-189, 1981. M. Ozgoren, B. Şahin, E. Pınar, H. Akıllı, “Experimental Investigation of Flow structure around a Sphere and Cylinder via Flow Visualization” 5. Ankara International Aerospace Conference, AIAC-2009-012, Ankara, 17-19 August, 2009. M. Ozgoren, B. Şahin, E. Pınar, H. Akıllı, “Küre ve Silindir Etrafındaki Akış Karakteristiklerinin Karşılaştırılması,” 17. Ulusal Isı Bilimi ve Tekniği Kongresi, s.73-79,24-27, Sivas, Haziran 2009. D. P. Hart, PIV error correction. Exp. Fluids. 29, 13–22, 2000. A. Fouras, J. Soria, “Accuracy of out-of-plane vorticity measurements derived from in-plane velocity field data,” Exp. Fluids, vol.25, 409–30, 1998. J. Westerweel, “Digital particle image velocimetry, theory and application,” Delft University Press, 1993. R. J. Adrian, “Particle-imaging techniques for experimental fluid mechanics,” Annu. Rev. Fluid Mech., vol.23, 261–304, 1991.