Roughness induced forced convective laminar-transitional micropipe flow: energy and exergy analysis

Variable fluid property continuity, Navier–Stokes and energy equations are solved for roughness induced forced convective laminar-transitional flow in a micropipe. Influences of Reynolds number, heat flux and surface roughness, on the momentum-energy transport mechanisms and second-law of thermodynamics, are investigated for the ranges of Re = 1–2,000, Q = 5–100 W/m² and ε = 1–50 μm. Numerical investigations put forward that surface roughness accelerates transition with flatter velocity profiles and increased intermittency values (γ); such that a high roughness of ε = 50 μm resulted in transitional character at Re _tra = 450 with γ = 0.136. Normalized friction coefficient (C _f^*) values showed augmentation with Re, as the evaluated C _f^* are 1.006, 1.028 and 1.088 for Re = 100, 500 and 1,500, respectively, at ε = 1 μm, the corresponding values rise to C _f^* = 1.021, 1.116 and 1.350 at ε = 50 μm. Heat transfer rates are also recorded to rise with Re and ε; moreover the growing influence of ε on Nusselt number with Re is determined by the Nu _ε=50 μm/Nu _ε=1 μm ratios of 1.086, 1.168 and 1.259 at Re = 500, 1,000 and 1,500. Thermal volumetric entropy generation values decrease with Re and ε in heating; however the contrary is recorded for frictional volumetric entropy generation data, where the augmentations in are more considerable when compared with the decrease rates of      

Experiments on the lift of a spinning sphere in a range of intermediate Reynolds numbers

 The lift force experienced by a spinning sphere moving in a viscous fluid, with constant linear and angular velocities, is measured by means of a trajectographic technique. Measurements are performed in the range of dimensionless angular velocities γ=aω/V lying between 1 and 6, and in the range of Reynolds numbers Re=2aV/ν lying between 10 and 140 (a sphere radius, ω angular velocity, V relative velocity of the sphere centre, ν fluid kinematic viscosity). A notable departure from the theoretical relationship at low Reynolds number, C _L =2γ, is obtained, the ratio C _L /γ being found to significantly decrease with increasing γ and increasing Re. The following correlation is finally proposed to estimate the lift coefficient in the range 10<Re<140: C _L ≅0.45+(2γ−0.45) exp (−0.075γ^0.4 Re ^0.7) Received: 20 May 1996/Accepted: 9 November 1997     

Some Reynolds number effects on two- and three-dimensional turbulent boundary layers

Experimental data for a two-dimensional (2-D) turbulent boundary layer (TBL) flow and a three-dimensional (3-D) pressure-driven TBL flow outside of a wing/body junction were obtained for an approach Reynolds number based on momentum thickness of Re _θ =23,200. The wing shape had a 3:2 elliptical nose, NACA 0020 profiled tail, and was mounted on a flat wall. Some Reynolds number effects are examined using fine spatial resolution (Δy ⁺=1.8) three-velocity-component laser-Doppler velocimeter measurements of mean velocities and Reynolds stresses at nine stations for Re _θ =23,200 and previously reported data for a much thinner boundary layer at Re _θ =5,940 for the same wing shape. In the 3-D boundary layers, while the stress profiles vary considerably along the flow due to deceleration, acceleration, and skewing, profiles of the parameter correlate well and over available Reynolds numbers. The measured static pressure variations on the flat wall are similar for the two Reynolds numbers, so the vorticity flux and the measured mean velocities scaled on wall variables agree closely near the wall. The stresses vary similarly for both cases, but with higher values in the outer region of the higher Re _θ case. The outer layer turbulence in the thicker high Reynolds number case behaves similarly to a rapid distortion of the flow, since stream-wise vortical effects from the wall have not diffused completely through the boundary layer at all measurement stations. Received: 9 June 2000/Accepted: 26 January 2001     

An experimental study of a two-dimensional plane turbulent wall jet

 Laser-Doppler measurements were conducted in a plane turbulent wall jet at a Reynolds number based on inlet velocity, Re ₀, of 9600. The initial development as well as the fully developed flow was studied. Special attention was given to the near-wall region, including the use of small measuring volumes and the application of specific near-wall data corrections, so that wall shear stresses were determined directly from the mean velocity gradient at the wall using only data below y ⁺=4. It was possible to resolve the inner peak in the streamwise turbulence intensity as well as the inner (negative) peak in the shear stress. Limiting values of (u′)⁺ and uv ⁺ were determined. Turbulence data from the outer region of the flow were compared to earlier hot wire measurements and large differences in the normal turbulence intensity and the shear stress were found. These differences can be attributed to high turbulence intensity effects on the hot-wires. Received: 17 October 1996 / Accepted: 8 December 1997     

Investigation of the vortex induced unsteadiness of a separation bubble via time-resolved and scanning PIV measurements

A transitional separation bubble on the suction side of an SD7003 airfoil is considered. The transition process that forces the separated shear layer to reattach seems to be governed by Kelvin–Helmholtz instabilities. Large scale vortices are formed due to this mechanism at the downstream end of the bubble. These vortices possess a three-dimensional structure and detach from the recirculation region, while other vortices are formed within the bubble. This separation of the vortex is a highly unsteady process, which leads to a bubble flapping. The structure of these vortices and the flapping of the separation bubble due to these vortices are temporally and spatially analyzed at angles of attack from 4° to 8° and chord-length based Reynolds numbers Re _c = 20,000–60,000 using time-resolved PIV measurements in a 2D and a 3D set-up, i.e., stereo-scanning PIV measurements are done in the latter case. These measurements complete former studies at a Reynolds number of Re _c = 20,000. The results of the time-resolved PIV measurements in a single light-sheet show the influence of the angle of attack and the Reynolds number. The characteristic parameters of the separation bubble are analyzed focusing on the unsteadiness of the separation bubble, e.g., the varying size of the main recirculation region, which characterizes the bubble flapping, and the corresponding Strouhal number are investigated. Furthermore, the impact of the freestream turbulence is investigated by juxtaposing the current and former results. The stereo-scanning PIV measurements at Reynolds numbers up to 60,000 elucidate the three-dimensional character of the vortical structures, which evolve at the downstream end of the separation bubble. It is shown that the same typical structures are formed, e.g., the c-shape vortex and the screwdriver vortex at each Reynolds number and angle of attack investigated and the occurrence of these patterns in relation to Λ-structures is discussed. To evidence the impact of the freestream turbulence, these results are compared with findings of former measurements.     

Approach to local axisymmetry in a turbulent cylinder wake

The objective of this experimental study is to characterise the small-scale turbulence in the intermediate wake of a circular cylinder using measured mean-squared velocity gradients. Seven of the twelve terms which feature in ε, the mean dissipation rate of the turbulent kinetic energy, were measured throughout the intermediate wake at a Reynolds number of Re _d  ≈ 3000 based on the cylinder diameter (d). Earlier measurements of the nine major terms of ε by Browne et al. (J Fluid Mech 179: 307–326 1987) at a downstream distance (x) of x = 420d and Re _d  ≈ 1170 are also used. Whilst departures from local isotropy are significant at all locations in the wake, local axisymmetry of the small-scale turbulence with respect to the mean flow direction is first satisfied approximately at x = 40d. The approach towards local axisymmetry is discussed in some detail in the context of the relative values of the mean-squared velocity gradients. The data also indicate that axisymmetry is approximately satisfied by the large scales at x/d ≥ 40, suggesting that the characteristics of the small scales reflect to a major extent those of the large scales. Nevertheless, the far-wake data of Browne et al. (1987) show a discernible departure from axisymmetry for both small and large scales.     

Conjugate heat transfer study of incompressible turbulent offset jet flows

In the present case, the conjugate heat transfer involving a turbulent plane offset jet is considered. The bottom wall of the solid block is maintained at an isothermal temperature higher than the jet inlet temperature. The parameters considered are the offset ratio (OR), the conductivity ratio (K), the solid slab thickness (S) and the Prandtl number (Pr). The Reynolds number considered is 15,000 because the flow becomes fully turbulent and then it becomes independent of the Reynolds number. The ranges of parameters considered are: OR = 3, 7 and 11, K = 1–1,000, S = 1–10 and Pr = 0.01–100. High Reynolds number two-equation model (k–ε) has been used for turbulence modeling. Results for the solid–fluid interface temperature, local Nusselt number, local heat flux, average Nusselt number and average heat transfer have been presented and discussed.     

Testing of Model Equations for the Mean Dissipation using Kolmogorov Flows

Direct Numerical Simulations (DNS) of Kolmogorov flows are performed at three different Reynolds numbers Re _λ between 110 and 190 by imposing a mean velocity profile in y-direction of the form U(y) = F sin(y) in a periodic box of volume (2π)³. After a few integral times the turbulent flow turns out to be statistically steady. Profiles of mean quantities are then obtained by averaging over planes at constant y. Based on these profiles two different model equations for the mean dissipation ε in the context of two-equation RANS (Reynolds Averaged Navier–Stokes) modelling of turbulence are compared to each other. The high Reynolds number version of the k-ε-model (Jones and Launder, Int J Heat Mass Transfer 15:301–314, 1972), to be called the standard model and a new model by Menter et al. (2006), to be called the Menter–Egorov model, are tested against the DNS results. Both models are solved numerically and it is found that the standard model does not provide a steady solution for the present case, while the Menter–Egorov model does. In addition a fairly good quantitative agreement of the model solution and the DNS data is found for the averaged profiles of the kinetic energy k and the dissipation ε. Furthermore, an analysis based on flow-inherent geometries, called dissipation elements (Wang and Peters, J Fluid Mech 608:113–138, 2008), is used to examine the Menter–Egorov ε model equation. An expression for the evolution of ε is derived by taking appropriate moments of the equation for the evolution of the probability density function (pdf) of the length of dissipation elements. A term-by-term comparison with the model equation allows a prediction of the constants, which with increasing Reynolds number approach the empirical values.     

Tailoring turbulence with an active grid

Using an active grid in a wind tunnel, we generate homogeneous shear turbulence and initiate turbulent boundary layers with adjustable properties. Homogeneous shear turbulence is characterized by a constant gradient of the mean velocity and a constant turbulence intensity. It is the simplest anisotropic turbulent flow thinkable, and it is generated traditionally by equipping a wind tunnel with screens which have a varying transparency and flow straighteners. This is not done easily, and the reachable turbulence levels are modest. We describe a new technique for generating homogeneous shear turbulence using an active grid only. Our active grid consists of a grid of rods with attached vanes which can be rotated by servo motors. We control the grid by prescribing the time-dependent angle of each axis. We tune the vertical transparency profile of the grid by setting appropriate angles of each rod such as to generate a uniform velocity gradient, and set the rods in flapping motion around these angles to tailor the turbulence intensity. The Taylor Reynolds number reached was R _λ = 870, the shear rate S = ∂U/∂y = 9.2 s⁻¹, the nondimensional shear parameter S ^*≡ Sq ²/ε = 12 and u = 1.4 ms⁻¹. As a further application of this idea we demonstrate the generation of a simulated atmospheric boundary layer in a wind tunnel which has tunable properties. This method offers a great advantage over the traditional one, in which vortex-generating structures need to be placed in the wind tunnel to initiate a fat boundary layer.     

On a Generalized Nonlinear K–ε Model and the Use of Extended Thermodynamics in Turbulence

A resent extension of the nonlinear K–ε model is critically discussed from a basic theoretical standpoint. While it was said in the paper that this model was formulated to incorporate relaxation effects, it will be shown that the model is incapable of describing one of the most basic such turbulent flows as is obvious but is described for clarity. It will be shown in detail that this generalized nonlinear K–ε model yields erroneous results for the Reynolds stress tensor when the mean strains are set to zero in a turbulent flow – the return-to-isotropy problem which is one of the most elementary relaxational turbulent flows. It is clear that K–ε type models cannot describe relaxation effects. While their general formalism can describe relaxation effects, the nonlinear K–ε model – which the paper is centered on – cannot. The deviatoric part of the Reynolds stress tensor is predicted to be zero when it actually only gradually relaxes to zero. Since this model was formulated by using the extended thermodynamics, it too will be critically assessed. It will be argued that there is an unsubstantial physical basis for the use of extended thermodynamics in turbulence. The role of Material Frame-Indifference and the implications for future research in turbulence modeling are also discussed. Received 19 February 1998 and accepted 23 October 1998     

Laminarisation and re-transition of a turbulent boundary layer subjected to favourable pressure gradient

 Experimental results are reported for the response of an initially turbulent boundary layer (Re _θ≈1700) to a favourable pressure gradient with a peak value of K≡(−υ/ρU ³ _E ) dp/dx equal to 4.4×10^-6. In the near-wall region of the boundary layer (y/δ<0.1) the turbulence intensity u′ scales roughly with the free-stream velocity U _E until close to the location where K is a maximum whereas in the outer region u′ remains essentially frozen. Once the pressure gradient is relaxed, the turbulence level increases throughout the boundary layer until K falls to zero when the near wall u′ levels show a significant decrease. The intermittency γ is the clearest indicator of a fundamental change in the turbulence structure: once K exceeds 3×10^-6, the value of γ in the immediate vicinity of the wall γ _s falls rapidly from unity, reaches zero at the location where K again falls below 3×10^-6 and then rises back to unity. Although γ is practically zero throughout the boundary layer in the vicinity of γ _s =0, the turbulence level remains high. The explanation for what appears to be a contradiction is that the turbulent frequencies are too low to induce turbulent mixing. The mean velocity profile changes shape abruptly where K exceeds 3×10^-6. Values for the skin friction coefficient, based upon hot-film measurements, peak at the same location as K and fall to a minimum close to the location where K drops back to zero. Received: 28 January 1998/Accepted: 8 April 1998     

A rotating elliptic airfoil in fluid at rest and in a parallel freestream

This paper reports results of DPIV measurements on a two-dimensional elliptic airfoil rotating about its own axis of symmetry in a fluid at rest and in a parallel freestream. In the former case, we examined three rotating speeds (Re _c,Ω = 400, 1,000 and 2,000), and in the later case, four rotating speeds (Ro _c,Ω = 2.4, 1.2, 0.6 and 0.4), together with two freestream velocities (Re _c,u  = 200 and 1,000) and two starting configurations of the airfoil (i.e., chord parallel to (α ₀ = 0°) or normal (α ₀ = 90°) to the freestream). Results show that a rotating airfoil in a stationary fluid produces two distinct types of vortex structures depending on the Reynolds number. The first type occurs at the lowest Reynolds number (Re _c,Ω = 400), where vortices shed from the two edges or tips of the airfoil dissipated quickly, resulting in the airfoil rotating in a layer of diffused vorticity. The second type occurs at higher Reynolds numbers (i.e., Re _c,Ω = 1,000 and 2,000), where the corresponding vortices rotated together with the airfoil. Due to the vortex suction effect, the torque characteristics are likely to be heavily damped for the first type because of the rapidly subsiding vortex shedding, and more oscillatory for the second type due to persistent presence of tip vortices. In a parallel freestream, increasing the tip-speed ratio (V/U) of the airfoil (i.e., decreasing the Rossby number, Ro _c,Ω) transformed the flow topology from periodic vortex shedding at Ro _c,Ω = 2.4 to the generation of a “hovering vortex” at Ro _c,Ω = 0.6 and 0.4. The presence of the hovering vortex, which has not been reported in literature before, is likely to enhance the lift characteristics of the airfoil. Freestream Reynolds number is found to have minimal effect on the vortex formation and shedding process, although it enhances shear layer instability and produces more small-scale flow structures that affect the dynamics of the hovering vortex. Likewise, initial starting configuration of the airfoil, while affecting the flow transient during the initial phase of rotation, has insignificant effect on the overall flow topology. Unfortunately, technical constraint of our apparatus prevented us from carrying out complimentary force measurements; nevertheless, the results presented herein, which are more extensive than those computed by Lugt and Ohring (1977), will provide useful benchmark data, from which more advanced numerical calculations can be carried out to ascertain the corresponding force characteristics, particularly for those conditions with the presence of hovering vortex.     

Reynolds Number Effects on the Coherent Dynamics of the Turbulent Horseshoe Vortex System

The adverse pressure gradient induced by a surface-mounted obstacle in a turbulent boundary layer causes the approaching flow to separate and form a dynamically rich horseshoe vortex system (HSV) in the junction of the obstacle with the wall. The Reynolds number of the flow (Re) is one of the important parameters that control the rich coherent dynamics of the vortex, which are known to give rise to low-frequency, bimodal fluctuations of the velocity field (Devenport and Simpson, J Fluid Mech 210:23–55, 1990; Paik et al., Phys Fluids 19:045107, 2007). We carry out detached eddy simulations (DES) of the flow past a circular cylinder mounted on a rectangular channel for Re = 2.0 × 10⁴ and 3.9 × 10⁴ (Dargahi, Exp Fluids 8:1–12, 1989) in order to systematically investigate the effect of the Reynolds number on the HSV dynamics. The computed results are compared with each other and with previous experimental and computational results for a related junction flow at a much higher Reynolds number (Re = 1.15 × 10⁵) (Devenport and Simpson, J Fluid Mech 210:23–55, 1990; Paik et al., Phys Fluids 19:045107, 2007). The computed results reveal significant variations with Re in terms of the mean-flow quantities, turbulence statistics, and the coherent dynamics of the turbulent HSV. For Re = 2.0 × 10⁴ the HSV system consists of a large number of necklace-type vortices that are shed periodically at higher frequencies than those observed in the Re = 3.9 × 10⁴ case. For this latter case the number of large-scale vortical structures that comprise the instantaneous HSV system is reduced significantly and the flow dynamics becomes quasi-periodic. For both cases, we show that the instantaneous flowfields are dominated by eruptions of wall-generated vorticity associated with the growth of hairpin vortices that wrap around and disorganize the primary HSV system. The intensity and frequency of these eruptions, however, appears to diminish rapidly with decreasing Re. In the high Re case the HSV system consists of a single, highly energetic, large-scale necklace vortex that is aperiodically disorganized by the growth of the hairpin mode. Regardless of the Re, we find pockets in the junction region within which the histograms of velocity fluctuations are bimodal as has also been observed in several previous experimental studies.     

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.     

PIV–LES analysis of channel flow rotating about the streamwise axis

The flow field of a channel rotating about the streamwise axis is analyzed experimentally and numerically. The current investigations were carried out at a bulk velocity based Reynolds number of Re_m = 2850 and a friction velocity based Reynolds number of Re_τ = 180, respectively. Particle-image velocimetry (PIV) measurements are compared with large-eddy simulation data to show earlier direct numerical simulation findings to generate too large a reverse flow region in the center region of the spanwise flow. The development of the mean spanwise velocity distribution and the influence of the rotation on the turbulent properties, i.e., the Reynolds stresses and the two-point correlations of the flow, are confirmed in both investigations. The rotation primarily influences those components of the Reynolds shear stresses, which contain the spanwise velocity component. The size of the correlation areas and thus the length scales of the flow generally grow in all three coordinate directions leading to longer structures. Furthermore, experimental results of the same channel flow at a significantly lower bulk Reynolds number of Re_{m, l} = 665, i.e., a laminar flow in a non-rotating channel, are introduced. The experiments show the low Reynolds number flow to become turbulent under rotation and to develop the same characteristics as the high Reynolds number flow.     

Secondary flow patterns and mixing in laminar pulsating flow through a curved pipe

Mixing by secondary flow is studied by particle image velocimetry (PIV) in a developing laminar pulsating flow through a circular curved pipe. The pipe curvature ratio is η = r ₀/r _c  = 0.09, and the curvature angle is 90°. Different secondary flow patterns are formed during an oscillation period due to competition among the centrifugal, inertial, and viscous forces. These different secondary-flow structures lead to different transverse-mixing schemes in the flow. Here, transverse mixing enhancement is investigated by imposing different pulsating conditions (Dean number, velocity ratio, and frequency parameter); favorable pulsating conditions for mixing are introduced. To obviate light-refraction effects during PIV measurements, a T-shaped structure is installed downstream of the curved pipe. Experiments are carried out for the Reynolds numbers range 420 ≤ Re_st ≤ 1,000 (Dean numbers 126.6 ≤ Dn ≤ 301.5) corresponding to non-oscillating flow, velocity component ratios 1 ≤ (β = U _max,osc/U _m,st) ≤ 4 (the ratio of velocity amplitude of oscillations to the mean velocity without oscillations), and frequency parameters 8.37 < (α = r ₀(ω/ν)^0.5) < 24.5, where α² is the ratio of viscous diffusion time over the pipe radius to the characteristic oscillation time. The variations in cross-sectional average values of absolute axial vorticity (|ζ|) and transverse strain rate (|ε|) are analyzed in order to quantify mixing. The effects of each parameter (Re_st, β, and α) on transverse mixing are discussed by comparing the dimensionless vorticities (|ζ _P |/|ζ _S |) and dimensionless transverse strain rates (|ε _P |/|ε _S |) during a complete oscillation period.     

Transient Growth in Exactly Counter-Rotating Couette–Taylor Flow

Transient growth due to non-normality is investigated for the Couette-Taylor problem with counter-rotating cylinders as a function of aspect ratio η and Reynolds number Re. For all Re≤500, transient growth is enhanced by curvature, i.e. is greater for η<1 than for η=1, the plane Couette limit. For fixed Re>130, it is found that the greatest transient growth is achieved for η on the linear stability boundary. Transient growth is approximately 20% higher near the Couette-Taylor linear stability boundary at Re=310, η=0.986 than at Re=310, η=1, near the threshold observed for transition in plane Couette flow. For 106<Re<130, the greatest transient growth occurs for a value of η between the linear stability boundary and one. For Re<106, the flow is linearly stable and the greatest transient growth occurs for a value of η less than one. The energy in the optimal inputs is primarily meridional; that in the optimal outputs is primarily azimuthal. Pseudospectra are calculated for two contrasting cases. Received 5 November 2001 and accepted 29 March 2002 Published online: 2 October 2002 Communicated by H.J.S. Fernando     

Simulations of Turbulent Flow in a Plane Asymmetric Diffuser

Large-eddy simulations (LES) of a planar, asymmetric diffuser flow have been performed. The diverging angle of the inclined wall of the diffuser is chosen as 8.5°, a case for which recent experimental data are available. Reasonable agreement between the LES and the experiments is obtained. The numerical method is further validated for diffuser flow with the diffuser wall inclined at a diverging angle of 10°, which has served as a test case for a number of experimental as well as numerical studies in the literature (LES, RANS). For the present results, the subgrid-scale stresses have been closed using the dynamic Smagorinsky model. A resolution study has been performed, highlighting the disparity of the relevant temporal and spatial scales and thus the sensitivity of the simulation results to the specific numerical grids used. The effect of different Reynolds numbers of the inflowing, fully turbulent channel flow has been studied, in particular, Re _b  = 4,500, Re _b  = 9,000 and Re _b  = 20,000 with Re _b being the Reynolds number based on the bulk velocity and channel half width. The results consistently show that by increasing the Reynolds number a clear trend towards a larger separated region is evident; at least for the studied, comparably low Reynolds-number regime. It is further shown that the small separated region occurring at the diffuser throat shows the opposite behaviour as the main separation region, i.e. the flow is separating less with higher Re _b . Moreover, the influence of the Reynolds number on the internal layer occurring at the non-inclined wall described in a recent study has also been assessed. It can be concluded that this region close to the upper, straight wall, is more distinct for larger Re _b . Additionally, the influence of temporal correlations arising from the commonly used periodic turbulent channel flow as inflow condition (similar to a precursor simulation) for the diffuser is assessed.     

Thermal characteristics of the wake shear layers from a slightly heated circular cylinder

In this paper, we investigate the thermal characteristics of wake shear layers generated by a slightly heated circular cylinder. Measurements of the fluctuating temperature were made in the region x/d = 0.6 to x/d = 3 (where x is the downstream distance from the cylinder axis and d is the cylinder diameter) using a single cold-wire probe. The Reynolds number Re was varied in the range 2,600–8,600. For Re = 5,500, simultaneous measurements were made with a rake of 16 cold wires, aligned in the direction of the mean shear, at x/d = 2 and 3. The results indicate that the passive temperature can be an effective marker of various instabilities of the wake shear layers, including the Kelvin–Helmholtz (KH) instability. The temperature data have confirmed the approximate Re ^m dependence of the KH instability frequency (f _KH) with different values of m over different ranges of Re, as reported previously in the literature. However, it is found that this power-law dependence is not exact, and a third-order polynomial dependence appears to fit the data well over the full range of Re. Importantly, it is found that the wake shear-layer instabilities can be grouped into three categories: (1) one with frequencies much smaller than the Bénard–Kármán-vortex shedding frequency, (2) one associated with the vortex shedding and (3) one related to the KH instability. The low-frequency shear-layer instabilities from both sides of the cylinder are in-phase, in contrast to the anti-phase high-frequency KH instabilities. Finally, the observed streamwise decrease in the mean KH frequency provides strong support for the occurrence of vortex pairing in wake shear layers from a circular cylinder, thus implying that both the wake shear layer and a mixing layer develop in similar fashion.