Direct numerical simulation of three-dimensional flow past a yawed circular cylinder of infinite length

Direct numerical simulation of flow past a stationary circular cylinder at yaw angles (α) in the range of 0–60° was conducted at Reynolds number of 1000. The three-dimensional (3-D) Navier–Stokes equations were solved using the Petrov–Galerkin finite element method. The transition of the flow from 2-D to 3-D was studied. The phenomena that were observed in flow visualization, such as the streamwise vortices, the vortex dislocation and the instability of the shear layer, were reproduced numerically. The effects of the yaw angle on wake structures, vortex shedding frequency and hydrodynamic forces of the cylinder were investigated. It was found that the Strouhal number at different yaw angles (α) follows the independence principle. The mean drag coefficient agrees well with the independence principle. It slightly increases with the increase of α and reaches a maximum value at α=60°, which is about 10% larger than that when α=0°. The root-mean-square (r.m.s.) values of the lift coefficient are noticeably dependent on α.     

A note on the numerical simulations of flow past a wavy square-section cylinder

The flow past a square-section cylinder with a geometric disturbance is investigated by numerical simulations. The extra terms, due to the introduction of mapping transformation simulating the effect of disturbance into the transformed Navier-Stokes equations, are correctly derived, and the incorrect ones in the previous literature are pointed out and analyzed. Furthermore, the relationship between the vorticity, especially on the cylinder surface, and the disturbance is derived and explained theoretically. The computations are performed at two Reynolds numbers of 100 and 180 and three amplitudes of waviness of 0.006, 0.025 and 0.167 with another aim to explore the effects of different Reynolds numbers and disturbance on the vortex dynamics in the wake and forces on the body. Numerical results have shown that, at the mild waviness of 0.025, the Kairmain vortex shedding is suppressed completely for Re = 100, while the forced vortex dislocation is appeared in the near wake at the Reynolds number of 180. The drag reduction is up to 21.6% at Re = 100 and 25.7% at Re = 180 for the high waviness of 0.167 compared with the non-wavy cylinder. The lift and the Strouhal number varied with different Reynolds numbers and the wave steepness are also obtained.     

Supercritical regimes of hypersonic flow past a planar yawed delta wing

Hypersonic viscous perfect gas flow past a planar delta wing in the viscous-inviscid interaction regime is considered. The effect of the yaw angle on the parameters of the laminar boundary layer on the cold wing and the formation of subcritical and supercritical flow regions is studied.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 4, pp. 151–158, July–August, 1995.     

Transverse subsonic flow past a cylinder

At around the critical Reynolds number Re = (1.5–4.0)·10⁵ there is an abrupt change in the pattern of transverse subsonic flow past a circular cylinder, and the drag coefficient C_x decreases sharply [1]. A large body of both experimental and computational investigations has now been made into subsonic flow past a cylinder [1–4]. A significant contribution to a deeper understanding of the phenomenon was made by [4], which gives a physical interpretation of a number of theoretical and experimental results obtained in a wide range of Re. Nevertheless, the complicated nonstationary nature of flow past a cylinder with separation and the occurrence of three-dimensional flows when two-dimensional flow is simulated in wind tunnels do not permit one to regard the problem as fully studied. The aim of the present work was to make additional experimental investigations into transverse subsonic flow past a cylinder and, in particular, to study the possible asymmetric stable flow regimes near the critical Reynolds number.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 154–157, March–April, 1980.     

Drag force on a free surface-piercing yawed circular cylinder in steady flow

The force distribution on a surface-piercing yawed cylinder surface differs significantly from that on a surface-piercing vertical cylinder. The established numerical model for flow past the surface-piercing yawed cylinder with yaw angles from −45° to 45° was solved by the standard large-eddy simulation (LES) methodology. Six cases at intervals of ±15° relative to the vertical were studied at the Reynolds number of 27 000 and the Froude number of 0.8 based on the cylinder diameter and free-stream velocity, among which the drag forces on four cylinders with yaw angles from −15° to 30° were tested for the validation of the LES approach. The results revealed that the time-averaged total drag coefficient for all cases increases with the increase of yaw angle compared to that of the surface-piercing vertical cylinder, even over 2.5 for the ±45°-yawed cylinders. The sectional drag coefficients for the negatively yawed cylinders are much greater than that for the vertical cylinder, and much less for the positively yawed cylinders. The unbalanced hydrostatic pressures on the inclined section are mainly responsible for those increment and decrement. Once the hydrostatic pressure was removed, the sectional drag coefficient on the mid-span of the positively yawed cylinder increases from the top section to the bottom, and decreases for the negatively yawed cylinder. The corresponding integrated total drag coefficient decreases with the increase of the yaw angle to ±15°, then increases with the further increase of the magnitude of yaw angle.     

Influence of wavy surfaces on coherent structures in a turbulent flow

We describe how outer flow turbulence phenomena depend on the interaction with the wall. We investigate coherent structures in turbulent flows over different wavy surfaces and specify the influence of the different surface geometries on the coherent structures. The most important contribution to the turbulent momentum transport is attributed to these structures, therefore this flow configuration is of large engineering interest. In order to achieve a homogeneous and inhomogeneous reference flow situation two different types of surface geometries are considered: (1) three sinusoidal bottom wall profiles with different amplitude-to-wavelength ratios of α = 2a/Λ = 0.2 (Λ = 30 mm), α = 0.2 (Λ = 15 mm), and α = 0.1 (Λ = 30 mm); and (2) a profile consisting of two superimposed sinusoidal waves with α = 0.1 (Λ = 30 mm). Measurements are carried out in a wide water channel facility (aspect ratio 12:1). Digital particle image velocimetry (PIV) is performed to examine the spatial variation of the streamwise, spanwise and wall-normal velocity components in three measurement planes. Measurements are performed at a Reynolds number of 11,200, defined with the half channel height h and the bulk velocity U _B. We apply the method of snapshots and perform a proper orthogonal decomposition (POD) of the streamwise, spanwise, and wall-normal velocity components to extract the most dominant flow structures. The structure of the most dominant eigenmode is related to counter-rotating, streamwise-oriented vortices. A qualitative comparison of the eigenfunctions for different sinusoidal wall profiles shows similar structures and comparable characteristic spanwise scales Λ _z = 1.5 H in the spanwise direction for each mode. The scale is observed to be slightly smaller for α = 0.2 (Λ = 15 mm) and slightly larger for α = 0.2 (Λ = 30 mm). This scaling for the flow over the basic wave geometries indicates that the size of the largest structures is neither directly linked to the solid wave amplitude, nor to the wavelength. The characteristic spanwise scale of the dominant eigenmode for the developed flow over the surface consisting of two superimposed waves reduces to 0.85 H. However, a scale in the order of 1.3 H is identified for the second mode. The eigenvalue spectra for the superimposed waves is much broader, more modes contribute to the energy-containing range. The turbulent flow with increased complexity of the bottom surface is characterized by an increased number of dominant large-scale structures with different spanwise scales.     

Capillary cavity flow past a circular cylinder

Cavity flow past a circular cylinder is considered accounting for the surface tension on the cavity boundary. The fluid is assumed to be inviscid and incompressible, and the flow is assumed to be irrotational. The solution is based on two derived governing expressions, which are the complex velocity and the derivative of the complex potential defined in an auxiliary parameter region. An integral equation in the velocity magnitude along the free surface is derived from the dynamic boundary condition. The Brillouin–Villat criterion is employed to determine the location of the point of flow separation. The cases of zero surface tension and zero cavitation number are obtained as limiting cases of the solution. Numerical results concerning the effects of surface tension and cavitation development on the cavity detachment, the drag force and the geometry of the free boundaries are presented over a wide range of the Weber and the cavitation numbers.     

Investigation of scale efect for the computation of turbulent flow around a circular cylinder

In order to investigate the scale efect of turbulent flow around a circular cylinder, two similarity numbers(criteria) based on turbulent kinetic and dissipation rates associated with the fluctuation characteristics of turbulence wake are deduced by analyzing the Reynolds averaged Navier–Stokes equations(RANS). The RNG k–ε models and finite volume method are used to solve the governing equations and the second-order implicit time and upwind space discretization algorithms are used to discrete the governing equations. A numerical computation of flow parameters around a two-dimensional circular cylinder with Reynolds numbers ranging from 102to 107is accomplished and the result indicates that the fluctuation of turbulence flow along the center line in the wake of circular cylinder can never be changed with increasing Reynolds numbers when Re3 × 106. This conclusion is useful for controlling the scale of numerical calculations and for applying model test data to engineering practice.     

Investigation of scale effect for the computation of turbulent flow around a circular cylinder

In order to investigate the scale effect of turbulent flow around a circular cylinder, two similarity numbers （criteria） based on turbulent kinetic and dissipation rates associ- ated with the fluctuation characteristics of turbulence wake are deduced by analyzing the Reynolds averaged NavierStokes equations （RANS）. The RNG k-s models and finite volume method are used to solve the governing equations and the second-order implicit time and upwind space discretization algorithms are used to discrete the governing equations. A numerical computation of flow parameters around a two-dimensional circular cylinder with Reynolds numbers ranging from 102 to l07 is accomplished and the result indicates that the fluctuation of turbulence flow along the center line in the wake of circular cylinder can never be changed with increasing Reynolds numbers when Re ≥ 3 × 10^6. This conclusion is useful for controlling the scale of numerical calculations and for applying model test data to engineering practice.     

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.     

The turbulent flow field around a circular cylinder

The flow field around a circular cylinder mounted vertically on a flat bottom has been investigated experimentally. This type of flow occurs in several technical applications, e.g. local scouring around bridge piers. Hydrogen bubble flow visualization was carried out for Reynolds numbers ranging from 6,600 to 65,000. The main flow characteristic upstream of the cylinder is a system of horse-shoe vortices which are shed quasi-periodically. The number of vortices depends on Reynolds number. The vortex system was found to be independent of the vortices that are shed in the wake of the cylinder. The topology of the separated flow contains several separation and attachment lines which are Reynolds number dependent. In the wake region different flow patterns exist for each constant Reynolds number.     

PIV measurements of flow past a confined cylinder

This study investigates the flow past a confined circular cylinder built into a narrow rectangular duct with a Reynolds number range of 1,500 ≤ Re _d ≤ 6,150, by employing the particle image velocimetry technique. In order to better explain the 3-D flow behaviour in the juncture regions of the lower and upper plates and the cylinder, respectively, as well as the dynamics of the horseshoe vortex system, both time-averaged and instantaneous flow data are presented for regions upstream and downstream of the cylinder. The size, intensity and interaction of the vortex systems vary substantially with the Reynolds number. Although the narrow rectangular duct with a single built-in cylinder is a geometrically symmetrical arrrangement, instantaneous flow data have revealed that the flow structures in both the lower and upper plate–cylinder junction regions are not symmetrical with respect to the centreline of the flow passage. The vortical flow structures obtained in side-view planes become dominant sometimes in the lower juncture region and sometimes in the upper juncture region in unsteady mode.     

The boundary layer on a slightly yawed cylinder

The velocity in a turbulent boundary layer on a long cylinder at a slight yaw to the free stream was measured using hot wire anemometry for yaw angles of ?0.55°?α?0.55°. The mean velocity profile retains a log region regardless of yaw with a slope that is slightly dependent upon the yaw angle. The boundary layer thickness increases nonlinearly with yaw angle, but the dimensionless distance from the wall of the maximum turbulence intensity is independent of yaw angle.     

Fluctuating flow past a cylinder at high frequencies

The flow in a laminar boundary layer for an arbitrary periodic main stream is considered at high frequencies when the fluid is incompressible. The analysis, incorporating length scales appropriate to the thin “Stokes” layer immediately adjacent to the surface and to the outer “Prandtl” boundary layer, involves expressing the dependent variables as mean parts plus superimposed periodic parts and expanding these in inverse powers of the frequency parameter in the two layers. Thus earlier approaches based on physical arguments are placed in the context of a systematic mathematical expansion scheme which is itself formulated for more general main stream velocities than hitherto. Expressions for skin friction and heat transfer are obtained and briefly discussed.     

Supersonic flow past a cylinder with a fluid flare

The aerodynamic parameters and the pressure distribution over the surface of a cylinder in a steady axisymmetric supersonic flow is studied within the framework of the inviscid perfect gas model in the absence and the presence of combined intense air injection fromthe flat face and the lateral surface into the shock layer. The purpose of the study is to investigate the effect of gas blowing from different regions of the cylindrical surface on the supersonic axisymmetric flow past the body.     

Laminar boundary layer on a yawed infinite cylinder

Using mixed momentum and energy integral equations, a simple quadrature method is developed to compute incompressible laminar boundary layer on a yawed infinite cylinder. As an illustration, the results — including various boundary layer thicknesses, form parameters and potential and surface streamlines — are obtained for a circular cylinder and compared with a known solution.     

Laminar flow of power-law fluids past a rotating cylinder

In this work, the continuity and momentum equations have been solved numerically to investigate the flow of power-law fluids over a rotating cylinder. In particular, consideration has been given to the prediction of drag and lift coefficients as functions of the pertinent governing dimensionless parameters, namely, power-law index (1 ≥ n ≥ 0.2), dimensionless rotational velocity (0 ≤ α ≤ 6) and the Reynolds number (0.1 ≤ Re ≤ 40). Over the range of Reynolds number, the flow is known to be steady. Detailed streamline and vorticity contours adjacent to the rotating cylinder and surface pressure profiles provide further insights into the nature of flow. Finally, the paper is concluded by comparing the present numerical results with the scant experimental data on velocity profiles in the vicinity of a rotating cylinder available in the literature. The correspondence is seen to be excellent for Newtonian and inelastic fluids.