Computation of high Reynolds number flow around a circular cylinder with surface roughness

Incompressible high-Reynolds-number flows around a circular cylinder are analyzed by direct integration of the Navier-Stokes equations using finite-difference method. A generalized coordinate system is used so that a sufficient number of grid points are distributed in the boundary layer and the wake. A numerical scheme which suppresses non-linear instability for calculations of high-Reynolds-number flows is developed. The computation of an impulsively started flow at Re = 1200 is compared with corresponding experimental observations, and excellent agreements are obtained.A series of computations are carried out on the flow around a circular cylinder with surface roughness. The height of the roughness in these computations is 0.5% of the diameter. The range of Reynolds numbers is from 10³ to 10⁵; no turbulence model is employed. Sharp reduction of drag coefficient is observed near Re = 2 × 10⁴, which indicates that the critical Reynolds number is captured in the present computation.     

Numerical simulation of Reynolds number effects on velocity shear flow around a circular cylinder

Three-dimensional Direct Numerical Simulation (DNS) and Large Eddy Simulation (LES) are performed to investigate the shear effects on flow around a circular cylinder at Reynolds numbers of Re=60–1000. The shear parameter, β, which is based on the velocity gradient, cylinder diameter and upstream mean velocity at the center plane of the cylinder, varies from 0 to 0.30. Variations of Strouhal number, drag and lift coefficients, and unsteady wake structures with shear parameter are studied, along with their dependence on Reynolds number. The presented simulation provides detailed information for the flow field around a circular cylinder in shear flow. This study shows that the Strouhal number exhibits no significant variation with shear parameter. The stagnation point moves to the high-velocity side almost linearly with shear parameter, and this result mainly influences the aerodynamic forces acting on a circular cylinder in shear flow. Both the Reynolds number and shear parameter influence the movement of the stagnation point and separation point. Mode A wake instability is suppressed into parallel vortex shedding mode at a certain shear parameter. The lift force increases with increasing shear parameter and acts from the high-velocity side to the low-velocity side. In addition, a simple method to estimate the lift force coefficient in shear flow is provided.     

A turbulence model for the heat transfer near stagnation point of a circular cylinder

A one-equation low-Reynolds number turbulence model has been applied successfully to the flow and heat transfer over a circular cylinder in turbulent cross flow. The turbulence length-scale was found to be equal 3.7y up to a distance 0.05 and then constant equal to 0.185 up to the edge of the boundary layer (wherey is the distance from the surface and is the boundary layer thickness).The model predictions for heat transfer coefficient, skin friction factor, velocity and kinetic energy profiles were in good agreement with the data. The model was applied for Re 250,000 and Tu0.07.Nomenclature µ,C _D Constants in the turbulence kinetic energy equation - C ₁,C ₂ Constants in the turbulence length-scale equation - Skin friction coefficient atx - D Cylinder diameter - F Dimensionless flow streamwise velocityu/u _e - k Turbulence kinetic energy =1/2 the sum of the squared three fluctuating velocities - K Dimensionless turbulence kinetic energyk/u _e ^/2 - I Dimensionless temperature (T–T _w )/(T –T _w ) - l Turbulence length-scale - l _e Turbulence length-scale at outer region - Nu _D Nusselt number - p Pressure - Pr Prandtl number - Pr _t Turbulent Prandtl number - Pr _k Constant in the turbulence kinetic energy equation - R Cylinder radius - Re _D Reynolds number u D/µ - Re _x Reynolds number u x/µ - R _K Reynolds number of turbulence - T Mean temperature - T Mean temperature at ambient - T _s Mean temperature at surface - Tu Cross flow turbulence intensity, - u Mean flow streamwise velocity - u Fluctuating streamwise velocity - u _e Mean flow velocity at far field distance - u Mean flow velocity at ambient - u* Friction velocity - v Mean velocity normal to surface - V Dimensionless mean velocity normal to surface - x,x ₁ Distance along the surface - y Distance normal to surface - Dimensionless pressure gradient parameter - Boundary layer thickness atu=0.9995u _e - Transformed coordinate iny direction - Fluid molecular viscosity - _t Turbulent viscosity - _eff + _t - µ Fluid molecular viscosity at ambient - Kinematic viscosity/ - Density - Density at ambient - _w Wall shear stress - _w,0 Wall shear stress at zero free stream turbulence     

Heat transfer and flow around a circular cylinder with tripping-wires

An experimental study was conducted on the heat transfer under the condition of constant heat flux and the flow around a circular cylinder with tripping-wires, which were affixed at ± 65° from the forward stagnation point on the cylinder surface. The testing fluid was air and the Reynolds number Red, based on the cylinder diameter, ranged from 1.2 × 10⁴ to 5.2×10⁴. Especially investigated are the interactions between the heat transfer and the flow in the critical flow state, in relation to the static pressure distribution along the cylinder surface and the mean and turbulent fluctuating velocities in the wake. It is found that the heat transfer from the cylinder to the cross flow is in very close connection with the width of near wake.     

A comparative assessment of fluid flow and heat transfer parameters in laminar flow around a circular cylinder and a sphere

Theoretical studies have been carried out for a comparative assessment of hydrodynamic boundary layer thickness, displacement thickness and shear stress at the wall for laminar flow around a circular cylinder and a sphere with the help of the approximate method due to Karman and Pohlhausen for two dimensional flow and the method as applied to bodies of revolution based on the work of F. W. Scholkemeier, respectively. Thermal boundary layer thickness and Nusselt number have been evaluated around the surface of the solids. Comparison is made with available solutions. The graphical presentation of the results depicts a concise and relative assessment of fluid flow and heat-transfer parameters for flow around cylinder and sphere.     

Analysis of buoyancy-aided convection heat transfer from a horizontal cylinder in a vertical duct at low Reynolds number

This paper reports a numerical study on buoyancy-aided steady convection heat transfer from a horizontal cylinder situated in a vertical adiabatic duct. Numerical results have been generated forH ₁/D=2.5, 4, 8,H/D=8, 16, 24,S/D=2, 4, 6, 20Re60, andRi up to 4. The placing of a horizontal cylinder in a vertical duct of smaller width results in significantly enhanced pure forced convection due to the blockage effect, but degrades appreciably the extent of buoyancy-aided enhancement in the heat transfer rate. Nevertheless, the presence of a vertical duct leads to an overall enhancement of mixed convection heat transfer coefficient relative to that without the confining duct. Moreover, the average Nusselt number is rather insensitive to the variation of either the position of the cylinder in the duct or the duct height in the investigated ranges of these geometric parameters.Dieser Artikel beschreibt eine numerische Studie über auftriebsunterstützte konvektive Wärmeübertragung von einem horizontalen Zylinder der in einem vertikalen adiabaten Kanal positioniert ist. Die numerischen Ergebnisse sind fürH ₁/D=2, 5, 4, 8,H/D=8, 16, 24,S/D=2, 4, 6 sowie 20Re60 undRi bis 4 berechnet worden. Die Anordnung des horizontalen Zylinders in einem schmaleren vertikalen Kanal führt auf Grund des Blockierungs-effektes zu einem deutlichen Anstieg der reinen Zwangskonvektion. Aber sie verschlechtert deutlich den Betrag der auftriebsbedingten Steigerung in der Wärmeübergangsrate. Trotzdem führt die Anwesenheit des vertikalen Kanals insgesamt zu einer Steigerung des Wärmeübergangskoeffizienten bei Mischkonvektion im Vergleich zur Abwesenheit des begrenzenden Kanals. Des weiteren ist die durchschnittliche Nusseltzahl von der Variation der Zylinderposition im Kanal oder der Kanalhöhe abhängig.     

A three-dimensional simulation of a steady approach flow past a circular cylinder at low Reynolds number

The three-dimensional (3D) unsteady viscous wake of a circular cylinder exposed to a steady approach flow is calculated using a fractional-step finite-difference/spectral-element method. The calculated flow fields at Reynolds numbers of 100 (2D) and 200 (3D) are examined in detail. The flow field at Re = 100 is 2D as expected, while the flow field at Re = 200 has distinct 3D features, with spanwise wavelengths of about 3.75 cylinder diameters. The calculated results produce drag and lift coefficients and Strouhal numbers that agree extremely well with the experimental values. These 3D values at Re = 200 are in better agreement with experimental values than the results of a 2D calculation at Re = 200, which is expected. © 1998 John Wiley & Sons, Ltd.     

Heat transfer in a second-order fluid for flow around a circular cylinder

The temperature distribution within the thermal boundary layer region due to the flow of an incompressible second-order fluid around a heated circular cylinder, maintained at a constant temperature higher than that of the fluid at infinity, has been obtained near the forward stagnation point by series expansion. The graph of the Nusselt number Nu for the Prandtl number P = 25 and the Eckert number E = 0.1 indicates that the non-Newtonian effect is to increase the heat flux from the cylinder to the liquid in the region 0° ? θ ? 15.7° and to decrease it in the region 15.7° < θ ? 0₀ where θ is the angular distance on the cylinder measured from the forward stagnation point. The critical point θ₀ at which Nu = 0, that is, where the effect of the frictional heating balances the effect of the temperature difference and there is no heat flux either from the cylinder or from the liquid, shifts towards the forward stagnation point with the increase of non-Newtonian effects.     

Coherent and turbulent process analysis in the flow past a circular cylinder at high Reynolds number

This article presents an experimental study of the turbulent flow past a circular cylinder at high Reynolds number by means of advanced optical measurements techniques. Following previous studies using standard PIV and stereoscopic PIV (3C PIV), TRPIV and 3C-TRPIV have been employed in low subsonic wind tunnel environment. The database consisting of statistical and time-dependent fields aims at providing a physical analysis of the coherent and turbulent part, as well as a proper basis for validation and improvement of recent turbulence modelling approaches for strongly detached flow at high Reynolds number. As the nonlinear interaction between the coherent and turbulent dynamics have to be taken into account in a model, particular attention is paid to a decomposition of the flow into a coherent and a turbulent part, and to the analysis of their dynamics. This is achieved both using phase averaging and Proper Orthogonal Decomposition. For phase averaging, the two first POD coefficients are used for the evaluation of the vortex shedding phase angle. Furthermore, selected results of a Detached Eddy Simulation which had been validated by means of the experiment, are also presented to contribute to the physical analysis. The present study's experimental data resolved in space and time allow the confirmation of the conditional averaging for the turbulent stresses evaluation, by alleviating their overestimation due to phase jitter that occurs between the trigger signal and the velocity, when the phase angle is determined from a wall pressure signal. A more accurate physical analysis of the flow is achieved, particularly regarding the occurrence of irregular vortex shedding.     

Fluid flow and heat transfer around circular cylinder – flat and curved plates combinations

Experimental studies were carried out to investigate the fluid flow and heat transfer around a heated circular cylinder which was placed at various distances of a wall boundary with different geometries (flat or curved plate) with subcritical Reynolds number ranging from 3.5×10³ to 10⁴. The effects of plate geometry (aspect ratio: W|H=1.0,1.5 and 2.0, and rim angle, φ=0°,60°,90°, and 120°) and gap ratio, (G|D=0.0,0.86,2.0,7.0,10.0) on the fluid flow and heat transfer characteristics (static pressure around cylinder surface, wake width, base pressure, pressure drag coefficients, velocity distribution, and both local and mean Nusselt numbers) were presented. Also flow visualization was carried out to illustrate the flow patterns around the cylinder at various gap ratios (G|D). It was found that the heat transfer and fluid flow characteristics are dependent on the plate geometry at all tested gap ratios, except for G|D=7.0 and 10.0, they are independent of the plate geometry.     

Numerical investigation on flow around circular cylinders in tandem arrangement at a subcritical Reynolds number

Three-dimensional fluid computations have been performed to investigate the flows around two circular cylinders in tandem arrangements at a subcritical Reynolds number, Re=2.2×10⁴. The center-to-center space between the cylinders was varied from twice the cylinder diameter to five times that, and the flows and fluid-dynamic forces obtained from the simulations are compared with the experimental results reported in the literature. Special attention is paid to the characteristics of the vortices shed from the upstream cylinder such as the convection, the impingement onto the downstream cylinder and the interaction with the vortices from the downstream cylinder. The effects of the vortices from the upstream cylinder on the fluid-dynamic forces acting on the downstream cylinder are discussed.     

Vortex shedding and aerodynamic forces on a circular cylinder in linear shear flow at subcritical Reynolds number

Vortex shedding and aerodynamic forces on a circular cylinder in a linear shear flow with its axis normal to the plane of the velocity shear profile at subcritical Reynolds number are investigated experimentally. The shear parameter β, which is based on the velocity gradient, cylinder diameter and upstream mean velocity at the center plane of the cylinder, varies from 0 to 0.27. The Strouhal number has no significant variation with the shear parameter. The time-mean base pressure increases and the fluctuating component of the base pressure decreases significantly with increasing shear parameter. Vortex shedding is suppressed by the velocity shear. Dislocation of the stagnation point takes place and this influences the pressure distribution around the cylinder together with the velocity shear. A mean lift force arises in the shear flow due to asymmetry of the pressure distribution, and it acts from the high velocity side to the low velocity side. In addition, the lift coefficient increases and the drag coefficient decreases with increasing shear parameter.     

Large-eddy simulation of circular cylinder flow at subcritical Reynolds number: Turbulent wake and sound radiation

The flows past a circular cylinder at Reynolds number 3900 are simulated using large-eddy simulation(LES) and the far-field sound is calculated from the LES results. A low dissipation energy-conserving finite volume scheme is used to discretize the incompressible Navier–Stokes equations. The dynamic global coefficient version of the Vreman's subgrid scale(SGS) model is used to compute the sub-grid stresses. Curle's integral of Lighthill's acoustic analogy is used to extract the sound radiated from the cylinder. The profiles of mean velocity and turbulent fluctuations obtained are consistent with the previous experimental and computational results. The sound radiation at far field exhibits the characteristic of a dipole and directivity. The sound spectra display the-5/3 power law. It is shown that Vreman's SGS model in company with dynamic procedure is suitable for LES of turbulence generated noise.     

Momentum and heat transfer from an asymmetrically confined circular cylinder in a plane channel

Unsteady momentum and heat transfer from an asymmetrically confined circular cylinder in a plane channel is numerically investigated using FLUENT for the ranges of Reynolds numbers as 10≤Re≤500, of the blockage ratio as 0.1≤β≤0.4, and of the gap ratio as 0.125≤γ≤1 for a constant value of the Prandtl number of 0.744. The transition of the flow from steady to unsteady (characterized by critical Re) is determined as a function of γ and β. The effect of γ on the mean drag and lift coefficients, Strouhal number (St), and Nusselt number (Nu _w ) is studied. Critical Re was found to increase with decreasing γ for all values of β. and St were found to increase with decreasing values of γ for fixed β and Re. The effect of decrease in γ on was found to be negligible for all blockage ratios investigated.     

Reynolds number dependence of the inertial range scaling of energy dissipation rate and enstrophy in a cylinder wake

The full energy dissipation rate and enstrophy are measured simultaneously using a probe consisting of four X-wires in the intermediate region of a cylinder wake for Taylor microscale Reynolds number in the range of 120–320. Longitudinal and transverse velocity increments are also obtained temporally using Taylor’s hypothesis. The inertial range scaling exponents indicate that the full enstrophy field has a stronger intermittency than does the full dissipation field for all the Reynolds numbers considered. The approximations of the energy dissipation rate and enstrophy based on isotropy are more intermittent than their corresponding true values. While the scaling exponents of the full energy dissipation rate remain approximately constant for different Reynolds numbers, those of the enstrophy decrease slightly and consistently with the increase of Reynolds number. It is conjectured that the scaling of the energy dissipation rate and the enstrophy may be the same when Reynolds number is extremely high, a trend that is consistent with that suggested by Nelkin (Phys Fluids 11:2202–2204, 1999; Am J Phys 68:310–318, 2000).     

Two-degree-of-freedom vortex induced vibration of low-mass horizontal circular cylinder near a free surface at low Reynolds number

Two-degree-of-freedom vortex induced vibration (VIV) of a low-mass zero-damping circular cylinder horizontally placed near a free surface at Re = 100 was numerically studied with an adaptive Cartesian cut-cell/level-set method. Two Froude numbers and various normalized submergence depths were considered. The results reveal that the Froude number affects the critical normalized submergence depth and possible physical mechanisms are proposed. The in-line vibration amplitude cannot be neglected. Proximity to a free surface strengthens and suppresses the VIV for low and high Froude numbers, respectively; increases the occurrence of amplitude modulation; and in general enhances the magnitude of the time-averaged lift coefficient, which is always negative. The phase lag of the transverse displacement behind the lift coefficient jumps at some reduced velocity, which strongly depends on the Froude number and normalized submergence depth. Regular trajectories exist only in cases with a small vibration amplitude or a large normalized submergence depth. The vortex structures in any case with large transverse amplitude basically originate from the alternative vortex shedding with the negative vortex weaker than the positive one. For the higher Froude number, an extra free surface positive vortex interacts with the vortices from the cylinder surface. The vibration frequency deviates from the natural structure frequency in fluids in the large-amplitude regime.