FLOW STRUCTURE AROUND A HORIZONTAL CYLINDER AT DIFFERENT ELEVATIONS IN SHALLOW WATER

Similar documents
Pulsating flow around a stationary cylinder: An experimental study

PASSIVE CONTROL OF FLOW STRUCTURE INTERACTION BETWEEN A SPHERE AND FREE-SURFACE

Figure 1. Schematic representation of the flow past a finite-height square prism mounted normal to a ground plane and partially immersed in a flat-pla

9.9 Coherent Structure Detection in a Backward-Facing Step Flow

THE FLUCTUATING VELOCITY FIELD ABOVE THE FREE END OF A SURFACE- MOUNTED FINITE-HEIGHT SQUARE PRISM

Time-resolved PIV measurements with CAVILUX HF diode laser

Hydrodynamic Instability and Particle Image Velocimetry

FLOW CHARACTERISTICS AROUND A CIRCULAR CYLINDER NEAR A PLANE BOUNDARY

COMPARISON OF FLOW STRUCTURE BEHAVIORS AROUND A SPHERE AND ITS PASSIVE CONTROL IN A BOUNDARY LAYER FLOW

Experimental investigation of Shedding Mode II in circular cylinder wake

PIV and LDV measurements behind a backward facing step

Assessment of Turbulence Models for Flow around a Surface-Mounted Cube

Interaction between a tethered sphere and a free surface flow

FLOW VISUALISATION AROUND A SOLID SPHERE ON A ROUGH BED UNDER REGULAR WAVES

ANALYSIS OF THE NON-STATIONARY FLOW AROUND A RECTANGULAR CYLINDER

INVESTIGATION OF FLOW BEHAVIOR PASSING OVER A CURVETURE STEP WITH AID OF PIV SYSTEM

A Passive Method to Control the Wake Flow behind a Circular Cylinder

Particle Velocimetry Data from COMSOL Model of Micro-channels

MODELLING THE FLOW AROUND AN ISLAND AND A HEADLAND: APPLICATION OF A TWO MIXING LENGTH MODEL WITH TELEMAC3D. Nicolas Chini 1 and Peter K.

Laser speckle based background oriented schlieren measurements in a fire backlayering front

ACTIVE SEPARATION CONTROL WITH LONGITUDINAL VORTICES GENERATED BY THREE TYPES OF JET ORIFICE SHAPE

Flow Field of Truncated Spherical Turrets

The effects of free surface and end cell on flow around a finite circular cylinder with low aspect ratio

Turbulencja w mikrokanale i jej wpływ na proces emulsyfikacji

CFD Analysis of 2-D Unsteady Flow Past a Square Cylinder at an Angle of Incidence

Module 3: Velocity Measurement Lecture 14: Analysis of PIV data. The Lecture Contains: Flow Visualization. Test Cell Flow Quality

Experimental and Numerical Investigation of Flow Structures around Cylindrical Bluff Bodies

Flow structure and air entrainment mechanism in a turbulent stationary bore

Effect of Leading Edge Porosity on the Flow Field of an Air Launched Grenade

7/13/05 Rough Draft. Copyright L. Pratt and J. Whitehead. Not for distribution.

Dispersion of rod-like particles in a turbulent free jet

Particle Image Velocimetry for Fluid Dynamics Measurements

MOMENTUM AND HEAT TRANSPORT INSIDE AND AROUND

GLASGOW 2003 INTEGRATING CFD AND EXPERIMENT

Direct numerical simulations of flow and heat transfer over a circular cylinder at Re = 2000

ScienceDirect. Physical Modeling of Flow around the Underwater Tidal Power

Measurement Techniques. Digital Particle Image Velocimetry

The Influence of End Conditions on Vortex Shedding from a Circular Cylinder in Sub-Critical Flow

THE INTERACTION OF A COLD ATOMISED SPRAY WITH A CIRCULAR CYLINDER

THREE DIMENSIONALITY IN THE NEAR FIELD OF A ROUND JET

Leading edge. Leading edge vortex, LE-V. Closed-type Boundary layer separation. Unsteady shedding by axial velocity

EXPERIMENTAL INVESTIGATION OF THE FLOW PATTERN BEHIND CYLINDER. Ing. Rut Vitkovičová, Ing. Vladislav Skála, Ing. Jan Čížek Ph.D.

Visualization Study on Flow Structures past a Rotating Circular Cylinder in a Uniform Stream

Separation in three-dimensional steady flows. Part 2 : DETACHMENT AND ATTACHMENT SEPARATION LINES DETACHMENT AND ATTACHMENT SEPARATION SURFACES

McNair Scholars Research Journal

APPLICATION OF WAVELET VECTOR MULTI-RESOLUTION TECHNIQUE TO PIV MEASUREMENTS

Cloud Cavitating Flow around an Axisymmetric Projectile in the shallow water

An Experimental and Computational Investigation of a 3D, l/h=5 Transonic Cavity Flow. Prof Kevin Knowles Dr Simon Ritchie Dr Nick Lawson

USE OF PROPER ORTHOGONAL DECOMPOSITION TO INVESTIGATE THE TURBULENT WAKE OF A SURFACE-MOUNTED FINITE SQUARE PRISM

Investigation of cross flow over a circular cylinder at low Re using the Immersed Boundary Method (IBM)

PARTICLE IMAGE VELOCIMETRY (PIV) AND VOLUMETRIC VELOCIMETRY (V3V) SYSTEMS

chinaxiv: v1

Vortex Generator Induced Flow in a High Re Boundary Layer

Volumetric 3-Component Velocity Measurements of Vortex Rings from Inclined Exits

EXPERIMENTELLE UNTERSUCHUNG SELBSTINDUZIERTER SCHWINGUNGEN FLEXIBLER STRUKTUREN IN LAMINARER STRÖMUNG

Particle Image Velocimetry Part - 3

Human beings are extremely interested in the observation of nature, as this was and still is of utmost importance for their survival.

The viscous forces on the cylinder are proportional to the gradient of the velocity field at the

Passive jet control of flow around a circular cylinder

IMAGE PROCESSING IN PIV

Surface flow visualisation over forward facing steps with varying yaw angle

Viscous/Potential Flow Coupling Study for Podded Propulsors

Measurements of Three-Dimensional Velocity Fields Under Breaking Waves

Visualization of three-dimensional vortex dynamics and fluid transport in translating plates, using defocusing DPIV

WONG HSI, J. J. MIAU,

Characterisation of an idealised offshore wind farm foundation, under waves and the combination of waves and currents

SIMULATION OF FLOW FIELD AROUND AND INSIDE SCOUR PROTECTION WITH PHYSICAL AND REALISTIC PARTICLE CONFIGURATIONS

Numerical Modeling of Flow Around Groynes with Different Shapes Using TELEMAC-3D Software

4. ANKARA INTERNATIONAL AEROSPACE CONFERENCE AIAC September, METU, Ankara

Vector Visualization. CSC 7443: Scientific Information Visualization

LAGRANGIAN VISUALIZATIONS and. EULERIAN DIAGNOSTICS of VORTEX BREAKDOWN

Learn more at

Effect of Position of Wall Mounted Surface Protrusion in Drag Characteristics At Low Reynolds Number

TURBULENT SHALLOW FLOW THROUGH VEGETATION

3D Dynamics of Vortex Structures in a Quasi Two-dimensional Jet

Backward facing step Homework. Department of Fluid Mechanics. For Personal Use. Budapest University of Technology and Economics. Budapest, 2010 autumn

Lecture # 11: Particle image velocimetry

On the flow and noise of a two-dimensional step element in a turbulent boundary layer

Slip effects in vortical structure behind cavitating propeller wake

Three-dimensional direct numerical simulation of wake transitions of a circular cylinder

AN EXPERIMENTAL STUDY AND SIMULATION WITH FLUENT

Multi-Scale Stereoscopic PIV measurements of a Jet in a Cross-Flow

VISUALIZATION OF MULTI-SCALE TURBULENT STRUCTURE IN LOBED MIXING JET USING WAVELETS

CFD MODELING FOR PNEUMATIC CONVEYING

Strengths of horseshoe vortices around a circular cylinder with an upstream cavity

CFD PREDICTION OF WIND PRESSURES ON CONICAL TANK

Measurements in Fluid Mechanics

Mathematical Modeling of Drag Coefficient Reduction in Circular Cylinder Using Two Passive Controls at Re = 1000

Vortex Generator Induced Flow in a High Re Boundary Layer

Numerical Study on the Ellipsoid

THE NUMERICAL INVESTIGATION ON THE EFFECTS OF ASPECT RATIO AND CROSS-SECTIONAL SHAPRE ON THE WAKE STRUCTURE BEHIND A CANTILEVER

Quantifying Three-Dimensional Deformations of Migrating Fibroblasts

Particle Image Velocimetry for Measuring Water Flow Velocity King Kuok Kuok, Po Chan Chiu

SPC 307 Aerodynamics. Lecture 1. February 10, 2018

COMPUTATIONAL FLUID DYNAMICS ANALYSIS OF ORIFICE PLATE METERING SITUATIONS UNDER ABNORMAL CONFIGURATIONS

Measurement of liquid film thickness in a micro parallel channel with interferometer and laser focus displacement meter

ITTC Recommended Procedures and Guidelines

Prerequisites: This tutorial assumes that you are familiar with the menu structure in FLUENT, and that you have solved Tutorial 1.

Recent Progress of NPLS Technique and Its Applications. in Measuring Supersonic Flows

3D Investigation of Seabed Stress Around Subsea Pipelines

Transcription:

FLOW STRUCTURE AROUND A HORIZONTAL CYLINDER AT DIFFERENT ELEVATIONS IN SHALLOW WATER 1 N.Filiz OZDIL, * 2 Huseyin AKILLI 1 Department of Mechanical Engineering, Adana Science and Technology University, Adana, Turkey 2 Department of Mechanical Engineering, Cukurova University, Adana, Turkey Abstract In this study, the flow characteristics around a horizontal circular cylinder were investigated using PIV technique in shallow water. The diameter of circular cylinder, height of shallow water and free stream velocity were kept constant during the experimental study as D=30 mm, h w =60 mm and U=167 mm/sn, respectively. Particle Image Velocimetry (PIV) was used to measure the instantaneous velocity vector field in the wake region at Reynolds number Re D =5000 based on the cylinder diameter. In order to investigate the effect of submergency, the cylinder was positioned at different elevations between bottom and free-surface (from 7.5 to 60 mm). The mean velocity vector field, corresponding streamline topology and Reynolds Stress correlation were obtained using 500 instantaneous images. As the blockage rate increases, the magnitude of jet-like flow velocity increases for h D /D=0.25, 0.5, 0.75 and 1 cases. Key words: horizontal cylinder, flow around the pipeline, particle image velocimetry 1. Introduction Vortex shedding behind a circular cylinder is an important fundamental question in various areas of engineering and science because of the practical importance. These areas include mechanical and aerospace engineering, power and process industries (turbine blades, heat exchanger tubes, cooling systems for nucleer power plants, power transmission lines), civil engineering (chimney stacks, bridges, buildings, radio telescopes, ice dams, offshore structures, power lins), and undersea technology (offshore drilling rigs, underwater pipelines, marine cables). The wake flow behavior behind bluff bodies in shallow water flow helps us to understand the transport characteristics of the flow. Akıllı and Rockwell [1] investigated the near wake of a vertical, circular cylinder in shallow water using a combination of visulation marker and particle image velocimetry tecnique (PIV). At the bed, the time- averaged streamline topology downstream of the base of the cylinder takes on a form known as an owl face of the first kind, which was originally defined for a completely different exterior flow. Immediately adjacent to *Corresponding author: Address: Adana Bilim ve Teknoloji Üniversitesi, Makine Mühendisliği Bölümü Adana 34355 Turkey. E-mail address: hakilli@cu.edu.tr

N.F. OZDIL et al./ ISITES2014 Karabuk - TURKEY 2003 1 the base of the cylinder, an additional system of saddle points is located at either end of a nodal line. At locations above bed, one of the two principle saddle points of the owl face of the first kind disappears and the principle foci are transformed from a stable to an unstable state. Vortex formation from a horizontal cylinder coincident with a free surface of shallow water has been examined by Kahraman et al. [2] using PIV technique. They found that the variation of reattachment location of the separated flow to the free surface is a strong function of the cylinder diameter and the Froude number. On the other hand, the flow characteristic of the horizontal cylinder placed on the plane boundary has been investigated by Akoz et al. [3] using particle image velocimetry technique. They found that intersection of the bed surface and cylinder enhanced the burial mechanisms hydrodynamically even in wake flow region. And the wake flow region was shortened in size in longitudinal direction as a function of Reynolds number. In the present study, the flow characteristics around a horizontal cylinder have been represented for different elevations. Aim of this experimental study first of all is to demonstrate the flow structure away from both the free and bottom surfaces in the near wake region. 2. Experimental set-up The experiments was carried out in a closed-loop water channel. The overall dimensions of the water channel were 8000 mm in length, 1000 mm in width and 750 mm in depth. In this study, the investigation of flow characteristics around a horizontal circular cylinder was carried out changing immersion level (h D /D) of the cylinder from the surface ranged from 0.25 to 2 with 0,25 increments. All experiments were carried out in shallow water. The water level of shallow water (h w ) was 60 mm. The diameter of horizontal cylinder was 30 mm. The schematic of experimental set up was shown in figure 1. The free stream velocity was U=167 mm/sec, which represented a value of Reynolds number based on cylinder diameter of between Re D = 5000. Using the PIV technique, instantaneous velocity vectors was measured in a region illuminated by a two-dimensional laser sheet. Velocity vector measurements was carried out using Dantec PIV system. The flow field illumination was provided by two Nd: Yag pulsed laser sources of a wavelength of 532 nm, each with a maximum energy output of 120 mj. The image capturing was performed by an 8-bit cross-correlation charge-coupled device (CCD) camera having a resolution of 1,600 pixels x 1,200 pixels, equipped with a Nikon AF Micro 60 f/2.8d

N.F. OZDIL et al./ ISITES2014 Karabuk - TURKEY 2004 lens. In the image processing, 32 x32 pixels rectangular effective interrogation windows was used. The total 3,168 (99 x 32) velocity vectors were obtained for an instantaneous velocity field at a rate of 15 frames/s. Figure 1. Experimental Set up The time interval between pulses was 1.750, and the thickness of the laser sheet illuminating the measuring plane was nearly 2 mm throughout the experiments. Erroneous vectors will be removed (less than 2%) and replaced by using interpolation between surrounding vectors in the post-processing step. Streamlines and circulation will be obtained by post processing of the velocity data. The overall field of view was 180 x 180 mm 2. In each experiment, 500 instantaneous images were captured and recorded. 3. Results and Discussion Fig. 2 and 3 show the flow characteristics which includes the time-averaged velocity vectors distribution (<V>) and Reynolds Stress correlation (<uv/u 2 >). In order to investigate the effect of the different value of dimensionless height, the cylinder was immersed in shallow water ranging from h D /D=0.25 to h D /D=2 with 0.25 increments in each case. As seen in Fig. 2 and 3, developing focus, F, and a primary circulating bubble are observed for h D /D=2 case, as mentioned by Kahraman et al. [2]. The primary circulating bubble spreads in a very large area in the wake region. In Fig. 3, positive values are characterized by solid lines, while negative values

N.F. OZDIL et al./ ISITES2014 Karabuk - TURKEY 2005 are described by dash lines. For h D /D=1.75 and 1.25 configurations, the Reynolds Stress (<uv/u 2 >) contours, are almost symmetrical and get closer to each other. A saddle point, S, and two foci, F, occur in close region of the cylinder. On the other hand, due to the effects of freesurface and bottom surface at h D /D=1.5 case, the minumum wake width is observed at h D /D=1.5 case because of the increase in momentum transfer. Velocity vector map (<V>), presented in the first column, shows that low velocity region which represents the wake region just downstream of the cylinder for this case. The velocity vectors (<V>) exterior to the bubble type wake region are much greater than those of in the wake region. The near wake is relatively narrow and symmetrical. Reynolds Stress concentrations (<uv/u 2 >), pronounced concentrations of Reynolds Stress (<uv/u 2 >) occur adjacent to the base of the cylinder. As the flow goes downstream, the Reynolds stress (<uv/u 2 >) concentration of the wake downstream of the cylinder increases initially, reaches a peak value of 0.09 and then decreases. Besides the well-defined Reynolds stress (<uv/u 2 >) pattern due to fluctuations in the shear layers, a weaker Reynolds stress (<uv/u 2 >) region very close to the base of the cylinder occurs as a result of the entrainment of free-stream flow into this wake region. However, for h D /D=1, 0.75, 0.5 cases, a well-defined swirl takes place downstream of the cylinder. And also, because of the extent of this primary circulating bubble, second plan measurements are taken. As seen in Fig. 2, the primary circulating bubble elongates in the second plan for h D /D=1, 0.75 cases. The other important observation is that there is a separation in the forward face of the horizantal cylinder for h D /D=1 and 0.75 cases as reported by Akoz et al. [3] Because the stagnation point is moved to upper side of the cylinder and this causes a perturbation at free-surface of the flow. Therefore, a small-scale focus, F, arises; lives only two immersion level, disappears immediately after h D /D=1, 0.75 cases. Besides, the time-averaged velocity vector (<V>), presented in Fig. 2, shows also that there is an upstream wake region on the shoulder of the cylinder and a high magnitude of velocity vectors (<V>) occurs along the shear layers. The Reynolds Stress (<uv/u 2 >), clearly indicates that positive contours dominate and negative contours disappear due to the movement of flow in the upward direction for h D /D=1, 0.75 cases. Namely, the shear layer can not form at the upper left quarter of the cylinder because the cylinder is positioned coincident with free-surface. This means that the momentum of flow decreases and Karman Street vanishes. For h D /D=0.5 and 0.25 configurations, well defined clusters of contour structures are observed in Reynolds Stress (<uv/u 2 >) correlations. However, the most important point is that the wake region expands in the streamwise direction

N.F. OZDIL et al./ ISITES2014 Karabuk - TURKEY 2006 as the dimensionless value of immersion level ratio (h D /D) increases. An effected area of flow region decrements because of the falling of the dimensionless immersion level ratio (h D /D). Positive contours are observed due to the presence of free-surface effect, whereas negative contours are not observed because of the lack of the viscous effect. At the h D /D=0.25 case, there is a smaller reverse flow than at the upper h D /D level. As the h D /D level rises, the center of the focus, F, in the wake region is moved further downward. The attached boundary layer on the bottom side of the cylinder elongates further downstream. This gives rise to a smaller wake at h D /D=0.25 case compared with upper h D /D level. Besides, there is an increase in flow velocity due to the reduction of the blockage rate for h D /D=0.25 case. Small wake doesn t subscribe to the clear transport of fluid in the upstream direction for this case. Spectra of streamwise velocity fluctuations evaluated at various locations in the flow field for h D /D=1.25 and 1.5, 1.75 case and bare cylinder case demonstrate different behavior of the flow. The locations taken for spectra are designated by numbers as indicated in Fig. 4. The spectra taken from the downstream locations of the cylinder (points 1-2-3) does not show a clear and distinct peak for h D /D=1.25 due to the effect of free surface. However, the spectra taken from the downstream locations of the cylinder (points 1-2-3) show a clear and distinct peak at the frequencies of f=0.909 Hz, 0.909 Hz and 0.909 Hz respectively, corresponding to a Strouhal number (St) of 0.163 for h D /D=1.5. Moreover, the spectra taken from the downstream locations of the cylinder (points 1-2-3) show a clear and distinct peak at the frequencies of f=0.928 Hz, 0.928 Hz and 0.341 Hz respectively, corresponding to a Strouhal number (St) of 0.166.

N.F. OZDIL et al./ ISITES2014 Karabuk - TURKEY 2007 Figure 2. Velocity vector field for horizontal circular cylinder

N.F. OZDIL et al./ ISITES2014 Karabuk - TURKEY 2008 Figure 3. Reynolds stress correlation

N.F. OZDIL et al./ ISITES2014 Karabuk - TURKEY 2009 h D /D=1.75 h D /D=1.5 h D /D=1.25 Figure 4. Spectra of streamwise velocity fluctuations evaluated at various locations in the flow field for h D /D=1.25 and 1.5, 1.75

N.F. OZDIL et al./ ISITES2014 Karabuk - TURKEY 2010 4. Conclusion The flow structure around the horizontal circular cylinder has been represented for different elevations. This study has helped understanding the flow structure away from the freesurface and wall and near free-surface and wall at the same time. Remarkable jet-like flow was observed for h D /D=0.25, 0.5, 0.75 and 1 cases. Only one swirl was remarkably observed for h D /D= 0.5, 0.75, 1 and 2 immersion level ratios. However, two foci and a saddle point was obtained for h D /D= 1.25, 1.5 and 1.75 cases. The overall form of the time-averaged velocity vector field are almost similar for these cases. Regarding the patterns of Reynolds shear stress, positive or negative contours remain horizontal and extend in the streamwise direction for h D /D= 0.5, 0.75 and 1 while positive and negative contours are symmetrical for h D /D= 1.25, 1.5 and 1.75 cases. Reynolds stress contours elongate along the bed due to the bed friction for h D /D=2. On the other hand, Reynolds shear stress contours for h D /D= 1 elongate along the free-stream direction due to the effect of the free surface. References [1] Akilli, H. and Rockwell, D., 2002. Vortex formation from a cylinder in shallow water. Physics of Fluids, Vol. 14, No. 9, pp. 2957-2967. [2] Kahraman, A., Ozgoren, M., Sahin, B., 2012. Flow Structure from a Horizantal Cylinder Coincident with a Free Surface in Shallow Water Flow. Thermal Science, 16(1): 93-107. [3] Akoz, MS., Sahin, B., Akilli, H., 2010. Flow characteristic of the horizontal cylinder placed on the plane boundary. Flow Measurement and Instrumentation, 21: 476-487. [4] Wolanski, E., Imberger, I. and Heron, M. L., 1984. Island Wakes in shallow coastal waters. Journal of Geophysical Research, 89(C6): 10553-10569. [5] Ozkan, GM., Oruc, V., Akilli, H., 2012. Flow around a cylinder surrounded by a permeable cylinder in shallow water. Exp Fluids, 53: 1751-1763 [6] Akilli, H., Sahin, B., Tumen, NF., 2005 Suppression of vortex shedding of circular cylinder in shallow water by a splitter plate. Flow Meas Instrum 16:211 219 [7] Akilli, H., Karakus, C., Akar, A., Sahin, B., Tumen, NF., 2008 Control of vortex shedding of circular cylinder in shallow water flow using an attached splitter plate. J Fluid Eng-T ASME 130:041401