 Open Access
 Total Downloads : 944
 Authors : Nidhul K, Sunil A S, Benphil C M
 Paper ID : IJERTV3IS081045
 Volume & Issue : Volume 03, Issue 08 (August 2014)
 Published (First Online): 01092014
 ISSN (Online) : 22780181
 Publisher Name : IJERT
 License: This work is licensed under a Creative Commons Attribution 4.0 International License
Effect of Domain Size and Grid Spacing on Flow Past a Circular Cylinder at Low Reynolds Number
Nidhul K 
Sunil A S 
Benphil C M 

Assistant Professor 
Associate Professor 
Assistant Professor 

Dept. of Mechanical Engg. 
Dept. of Mechanical Engg. 
Dept. of Mechanical Engg. 

Royal college of Engineering and Technology Thrissur, India 
Govt. 
Engineering College Thrissur Thrissur, India 
College of Engineering Adoor Pathanamthitta, Thrissur 
Abstract In this paper, the effect of size of rectangular flow domain and grid spacing on flow past a circular cylinder at low Reynolds number is numerically investigated using ANSYS Fluent. Two dimensional steady incompressible flow past a circular cylinder is solved numerically using finite volume method with second order accuracy in space. Pressurevelocity coupling scheme employed is SIMPLEC. A domain independent study was conducted by varying the size of flow domain to ensure that the simulation results were independent of domain size. Then on the flow domain chosen, grid independence study is conducted so as to optimize computational hardware resource and time while ensuring that the results predicted are accurate. The flow pattern and flow characteristics past a circular cylinder at low Reynolds number were studied by means of vorticity magnitude contour, flow streamlines, coefficient of drag and variation of xcomponent of velocity.
Keywords grid; circular cylinder; Reynolds number; incompressible flow; coefficient of drag, coefficient of pressure

INTRODUCTION
Vortex shedding behind bluff bodies is of concern for many engineering applications. Bluff bodies are structures with shapes that significantly disturb the flow around them, as opposed to flow around a streamlined body. A bluff body is one in which the length in the flow direction is close to or equal to the length perpendicular to the flow direction. Examples of bluff bodies include circular cylinders, square cylinders and rectangular cylinders. Flow around a circular cylinder is a fundamental fluid mechanics problem of practical importance such as in submarines, off shore structures, bridge piers, pipelines etc. The laminar steady viscous flow behind a circular cylinder has been the subject of numerous experimental and numerical studies.
Zakir Faruquee et.al [1] conducted domain independent study for flow past circular cylinder using circular domain whose farfield boundary radius was varied from 10D60D. In all cases, the same cell distribution was used for common subdomains. The influence of the farfield boundary on the drag coefficient was studied. Behr et al. [2] studied the influence of the location of the lateral boundaries on the unsteady flow past a circular cylinder at Re = 100. They have suggested that the lateral boundaries should be away from the cylinder by a distance of at least 8 cylinder diameters. Inoue et.al [3] constructed a nonuniform mesh but divided the computational domain into three regions, each with a different
grid ratio. Mohamed Sukri Mat Ali et.al [4] numerically investigated the sensitivity of the computed flow field to flow parameters for a flow with Reynolds number 150.They constructed computational meshes based on reasonable estimates of cell size and grid stretching ratios. Sumer et.al [5] investigated and reported that the flow field over the circular cylinder is symmetric at low values of Reynolds number. As the Reynolds number increases, flow begins to separate behind the cylinder causing vortex shedding which is an unsteady phenomenon. Kovasznay [6] and Roshko [7] determined experimentally that for Re less than 40, the flow remains laminar, twodimensional and steady, but can become unstable under the influence of external excitation. Braza et.al
[8] showed that when the Reynolds number is increased beyond 40, not every destabilization of the flow, random or not, can be damped out. 
NUMERICAL SIMULATION PROCEDURE

Governing equations
A single circular cylinder subjected to a uniform upstream air flow in an open atmosphere is the physical model of the problem. Flow over the cylinder is governed by the partial differential equations derived from the laws of conservation of mass (continuity equation) and momentum (NavierStokes equations), which can be written as follows:
(1)
(2)
(3)
Where is the fluid density, Âµ is the fluid viscosity, V is the velocity vector of the flow field, p is the pressure, and u and v are the velocity components in the x and ydirections, respectively. Bx and By are also the body forces per unit volume, which are negligible in the present study. The fluid is assumed to be incompressible, and its properties has been taken as =1.225 kg/m3 and Âµ=1.7894Ã—105 kg/m s.

Computational domain
The governing equations (1)(3) have been applied in the computational domain (Fig.1) of the physical model. A rectangular domain is employed in this study. To minimize the effect of the farfield boundary conditions, the outer boundary of the computational domain must be placed sufficiently far away from the cylinder.
Symmetry
Pressure
Velocity outlet
The application of appropriate boundary conditions at the cylinder surface and at the outer boundary is crucial for accurate simulation. A uniform approach velocity was imposed on the upstream of the outer boundary of the computational domain, while symmetry condition was imposed on top and bottom freestreams and a pressure outlet condition was imposed on the downstream. Noslip condition was applied along the cylinder wall. In this study, the QUICK method of interpolation was used for solving velocity and momentum. The SIMPLEC method was employed for the pressure velocity coupling for all simulations. The second
inlet
re 1: Schematic of computational dom
Figu
Noslip wall
Symmetry
boundary conditions.
ain and
order interpolation scheme was applied for pressure. Convergence criterion employed is 106.


GRID INFLUENCE
Within the scope of predicting the flow field around a
A computational domain extending up to a length of 4D (i.e., 4 times cylinder diameter) in the upstream, to a length of 4D in transverse directions and 9D in the downstream direction from the center of the cylinder was initially employed. Then, to determine an appropriate location for the outer boundaries, length of boundaries was varied upto 16D in upstream and transverse direction and 21D in downstream direction. In all cases, the same cell distribution was used for common subdomains. The influence of the farfield boundaries on the drag coefficient is given in Table 1 at Re = 40.
TABLE 1
Effect of farfield boundaries on the drag coefficient for a circular cylinder at Re = 40
Farfield boundaries
No. of cells
CD
4D, 4D and 9D
33120
1.89
5D, 5D and 10D
40392
1.79
6D, 6D and 11D
49704
1.72
7D, 7D and 12D
60480
1.68
8D, 8D and 13D
74200
1.65
9D, 9D and 14D
91104
1.63
10D, 10D and 15D
110424
1.61
11D, 11D and 16D
134140
1.60
12D, 12D and 17D
162960
1.59
13D, 13D and 18D
197274
1.58
14D, 14D and 19D
240640
1.57
15D, 15D and 20D
291588
1.56
16D, 16D and 21D
353918
1.56
As seen from this table, moving the outer boundaries from 15D, 15D and 20D, does not have any noticeable effect on the flow parameter. Therefore, computational domain having an upstream length of 15D, transverse length of 15D and downstream length of 20D was chosen to simulate the unbounded flow past circular cylinder.
square cylinder using numerical analysis, many similar investigations have been made, but the results always show small discrepancies even though the overall global trends are similar. One of the reasons for these discrepancies is the difference in the construction of the mesh. Generally, the accuracy of a numerical solution increases as the number of cells increases. But employing a larger number of cells requires more computer hardware resource and computing time. All simulations performed in this study employed quadrilateral cells. Grid spacing was gradually increased from the cylinder surface towards the outer boundaries in order to avoid sudden distortion and skewness, and to provide a sufficiently clustered mesh near the cylinder wall where the flow gradients are large and to capture flow regime within the boundary layer.
Grid independency was achieved in this study according to the procedure outlined by Yunus A. cengel and Michael A. Boles [10].Several meshes, with increasing refinement, were tested to ensure that the solution was independent of the mesh. These meshes and the drag coefficient (CD) predictions are reported in Table 2 for the case of a circular cylinder at Re = 40. It is observed that initially with 27300 cells, the value of coefficient of drag is overpreicted. Whereas the value does not vary for mesh finer than 81952 cells.
TABLE 2
Effects of cell density on drag coefficient for a circular cylinder at Re = 40
No. of cells
CD
27300
1.56
0
39312
1.56
2
56560
1.56
3
81952
1.56
4
116392
1.56
4

RESULTS AND DISCUSSIONS
From the domain independence study, it was observed that rectangular domain with dimensions 15D in upstream direction, 15D in transverse direction and 20D in the downstream direction predicts the results accurately for flow past a circular cylinder at Re=40. Then grid independence study was conducted on this domain with grid refinement starting from 27300 cells. It was observed that grid refinement beyond 81952 cells did not make any significant changes on flow parameters. Variation in the value of drag coefficient was taken for grid refinement study.
2D steady and incompressible flow over a circular cylinder was simulated at Re=40 using the above mentioned computational domain and grid. Flow at Re=40 is laminar and steady. As flow approaches the leading point of a circular cylinder, the fluid is brought to rest, causing a pressure rise and a boundary layer develops along the surface. The leading point is called the forward stagnation point. For very low Reynolds number (Re < 5), the flow splits into two symmetrically equal parts and passes around the cylinder to the rear stagnation point. When Reynolds number is increased, small size eddies develop in the immediate vicinity of the rear stagnation point of the cylinder, where the velocity is extremely small and pressure gradient is unfavorable. The recirculation zone remains steady, two dimensional and symmetrical about the streamwise center line of the flow. The contours of vorticity magnitude at Re=5 and Re=40 are shown in fig.2.
(a)
(b)
Figure 2: Contours of vorticity magnitude

At Re=5 and (b) At Re=40
The flow streamlines for flow past circular cylinder at Re=5 and Re=40 are shown in fig.3. At Re=40, only recirculation eddies are formed immediately downstream of the cylinder, near the rear stagnation point as indicated by streamlines.
(a)

Figure 3: Flow streamlines

At Re=5 and (b) At Re=40


Variation of xcomponent along the wake centerline is shown by nondimensionalised parameter of xcomponent velocity. Fig.4 shows the variation of xcomponent velocity along wake centerline at Re=5 and Re=40.
Figure 3: U/Uref along the wake centerline
REFERENCES

Zakir Faruquee, David SK. Ting, Amir Fartaj, Ronald M. Barron and Rupp Carriveau, The effects of axis ratio on laminar fluid flow around an elliptical cylinder, International Journal of Heat and Fluid Flow 28 (2007) 11781189.

Behr, M., Hastreiter, D., Mittal, S., Tezduyar, T.E., flow past a circular cylinder: dependence of the computed flow field on the location of the lateral boundaries, Comput. Methods Appl. Mech. Eng.(1995) 123, 309316.

Inoue, Iwakam and Hatakeyama,Aeolian tones radiated from flow past two square cylinders in a sidebyside arrangement. Physics of Fluids,(2006) 18(4), 046104.

Ali, M.S.M., Doolan, C.J., Wheatley, V., Grid convergence study for a two dimensional simulation of flow around a square cylinder at a low Reynolds number, In: Seventh International Conference on CFD in The Minerals and Process Industries.(2009) CSIRO, Melbourne, Australia.

D. Sumner, J. L. Heseltine and O. J. P. Dansereau, Wake structure of a finite circular cylinder of small aspect ratio, Experiments in Fluids, November 2004, Volume 37, Issue 5, pp 720730.

Kovasznay, L.S.G., Hotwire investigation of the wake behind cylinders at low Reynolds number, Proc. Roy. Soc. A (1949) 198, 174.190.

Roshko, A., On the development of turbulent wakes from vortex streets. National Advisory Committee for Aeronautics,(1954) Report 1191.

Braza, M., Chassaing, P. and Ha Minh, H., Numerical study and physical analysis of the pressure and velocity fields in the near wake of a circular cylinder. J. Fluid Mech.(1986) 165, 79130.

Okajima, A., Numerical analysis of the flow around an oscillating cylinder, In P.W. Bearman (Ed.), Proc. 6th Int. Conf. FlowInduced Vibration, London, UK,(1995) 1012 April, pp. 17. Balkema, Rotterdam.

Yunus A. Cengel & Michael A. Boles, Fluid Mechanics:Fundamentals and Applications, Published by McGrawHill.