The supercritical Strouhal and critical Reynolds numbers associated with the drag crisis of circular cylinders in cross-flow are known to be dependent on surface roughness and turbulence levels. To better understand the role of free-stream turbulence, we consider the traditional model of oscillating lift forcing together with a randomly fluctuating disturbance associated with fluid turbulence. To the extent that these two components are additive, statistically independent, and Gaussian, we are able to derive a closed-form expression describing the relationship between the effective shedding rate and free-stream turbulence. A semi-empirical model of the critical Reynolds number accounting for fluid turbulence is also presented.
I. INTRODUCTION
Flow tests with subresonant cylinders, more commonly known as flexible cylinders with shedding rates well below resonance conditions,1–5 have suggested that the reported supercritical Strouhal numbers associated with the drag crisis decrease in a near-linear fashion with increasing turbulence intensity, while another series of tests6,7 have indicated a nonlinear increase.
In this paper, a simplified model of vortex shedding in the presence of turbulent velocity fluctuations is presented, and this is followed by an investigation of the integral scale length of the fluctuations to quantify its spectral content. Stochastic analysis is then applied to develop an expression for the effective shedding rate and its relationship with turbulence intensity and its spectral content.
In collecting the data associated with the drag crisis, we find that two categories of cylinders are involved. First, we have those with measurably smooth surfaces in pressurized wind tunnels and pipe flow tests,1–5 and, second, those in atmospheric wind tunnel flow tests with less than ideal surfaces.6,7 In the first case, with integral scale lengths of free-stream turbulence greater than the cylinder diameter, the supercritical Strouhal number is elevated over a significant range of Reynolds numbers at low turbulence, but decreases with increased turbulence intensity.1–5 In the second case, flows with integral scales less than the cylinder diameter, the supercritical Strouhal number is increased in the presence of moderate turbulence level.6,7 The critical Reynolds number data in both groups of tests are found to decrease with increasing turbulence intensity.
Occasionally, the supercritical shedding spectrum lacks an identifiable peak in reduced frequency coordinates.4,5 When this occurs, we take the upper frequency half-power point of the oscillating lift spectrum as the effective Strouhal number5 that would be associated with resonance of more flexible or longer cylinders.
II. FLUID FORCING
Unlike earlier models involving vortex shedding in the presence of free-stream turbulence, such as that as in Ref. 8, we consider shedding rates well below the onset of flow-induced resonance. Avoidance of resonance is mandated in flow tests used to characterize Reynolds-number-dependent vortex shedding parameters1–7 and the role played by free-stream turbulence.
The Reynolds-number-dependent forcing parameters are (cL, cR), with low turbulence shedding rate ω0, uniformly distributed phase disturbances ϕ(t), and zero mean velocity fluctuations v′(t). This model of vortex shedding with additive turbulent forcing is analogous to those suggested previously, although in spectral form.9,10 The phase disturbances are considered to be independent of the velocity fluctuations. Were ϕ(t) deterministic, its rate of change would act to modulate the effective shedding rate directly.
The root mean square (RMS) of the fluctuations, normalized by the mean fluid velocity, defines the turbulence intensity as , while the spectral bandwidth of the fluctuations is inferred from reported integral scale lengths using autocorrelation methods.5,9–11 The magnitude of the turbulence forcing parameter may be estimated by examining the changes in the oscillating lift parameter that result from increased turbulence intensity, as in Refs. 6 and 7. As most flow tests lack such data, the forcing parameter ratio cR/cL, is considered an adjustable parameter.
III. TURBULENCE SPECTRA
This parameter proves useful in interpreting Strouhal number behavior in a given Reynolds number range. It also addresses the question raised by Zdravkovich15 regarding the lack of adequate spectral details of free-stream turbulence,6 where the Strouhal numbers are found to increase with increasing turbulence. Unlike models and tests involving flow-induced resonances,8,16,17 Eq. (7) is defined in terms of the turbulence-free shedding rate of a subresonant cylinder to facilitate identification of potential resonance conditions that must be avoided at higher velocities.
In Sec. IV, we address the effective shedding rate in the presence of free-stream turbulence. This is accomplished by using random noise theory to quantify the zero-crossing rate of the shedding forces, and hence the effective Strouhal number. While the theory can be applied in other Reynolds regimes, the present focus is on Strouhal numbers associated with the drag crisis.
IV. THE ZERO-CROSSING PROBLEM
Traditional forcing models17 are used routinely for engineering risk assessments of intrusive pipe fittings. However, consideration of elevated, supercritical Strouhal numbers with isolated cylinders has remained controversial, in spite of documented failures, and subsequent flow tests involving the drag crisis.5,18 To better understand the role of free-stream turbulence in such cases, consider the following.
When the leading term in Eq. (2) is characterized by stationary random fluctuations in magnitude with mean frequency and uniformly distributed phase disturbances {0, …, 2π}, it is formally described by Gaussian statistics.19–21 Under steady flow conditions, the central limit theorem commonly applies, and oscillating lift forcing may be considered as Gaussian, if only as an approximation.
The partial derivative of Eq. (13), ∂St/∂TI, is proportional to (α2 − 3), revealing a bifurcation or branching in the shedding rate dependence on TI. As a result, we have a decreasing St(TI) for α2 < 3, and increasing St(TI) for α2 > 3, with St(TI) independent of TI for α2 = 3. When higher-order integrable spectra are considered, both the bandwidth parameter and the branching condition of Strouhal number behavior are altered with slightly improved data fits in Fig. 1.
It is readily shown that the steady-state energy transfer between the fluid forcing and the cylinder is proportional to the shedding rate. As a consequence, both experimentally and within the framework of the model presented here, we find that turbulence causes a decrease in the shedding rate for α2 < 3 and an increase for α2 > 3.
V. FLOW TEST DATA AND MODEL PREDICTIONS
Figure 1 compares the supercritical Strouhal numbers1–7 with the results of Eq. (13) using the Strouhal number data associated with the drag crisis in each test together with the respective turbulence intensities and reported integral scale lengths. The data in Figs. 1–3 are summarized in Table I.
Supercritical Strouhal number vs TI for two datasets. The dashed lines represent Eq. (13) for the parameters shown and the solid lines highlight the datasets involved.
Supercritical Strouhal number vs TI for two datasets. The dashed lines represent Eq. (13) for the parameters shown and the solid lines highlight the datasets involved.
Turbulence bandwidth parameters vs TI using reported data to calculate the bandwidth parameter according to Eq. (7).
Turbulence bandwidth parameters vs TI using reported data to calculate the bandwidth parameter according to Eq. (7).
Reference . | TI (%) . | Λ/d . | α . | St(0) . | St(TI) . | ⟨Rdcrit⟩ . |
---|---|---|---|---|---|---|
Bearman2 | 0.20 | ⋯ | ⋯ | 0.46 | 0.46 | 346 410 |
Schewe1 | 0.40 | ⋯ | ⋯ | 0.48 | 0.48 | 324 037 |
Achenbach and Heinecke3 | 0.46 | ⋯ | ⋯ | 0.51 | 0.51 | 273 861 |
So and Savakar4 | 0.50 | ⋯ | 1.18 | 0.50 | 0.49 | 256 905 |
Kawamura et al.5 | 1 | 5.7 | 0.11 | 0.48 | 0.48 | 223 607 |
Kawamura et al.5 | 5 | 0.53 | 1.18 | 0.48 | 0.42 | 146 969 |
So and Savakar4 | 9.5 | ⋯ | 1.18 | 0.23 | 0.34 | 60 000 |
Kawamura et al.5 | 13 | 0.089 | 7.93 | 0.48 | 0.281 | 63 246 |
Blackburn and Melbourne7 | 18 | 0.53 | 2.43 | 0.23 | 0.23 | ⋯ |
Cheung and Melbourne6 | 0.4 | 0.521 | 2.47 | 0.23 | 0.23 | 229 169 |
Cheung and Melbourne6 | 1.6 | 0.521 | 2.47 | 0.23 | 0.25 | 141 421 |
Cheung and Melbourne6 | 4.4 | 0.521 | 2.47 | 0.23 | 0.30 | 109 545 |
Cheung and Melbourne6 | 6.8 | 0.521 | 2.47 | 0.23 | 0.32 | 69 282 |
Cheung and Melbourne6 | 9.1 | 0.521 | 2.47 | 0.23 | 0.34 | 77 460 |
Reference . | TI (%) . | Λ/d . | α . | St(0) . | St(TI) . | ⟨Rdcrit⟩ . |
---|---|---|---|---|---|---|
Bearman2 | 0.20 | ⋯ | ⋯ | 0.46 | 0.46 | 346 410 |
Schewe1 | 0.40 | ⋯ | ⋯ | 0.48 | 0.48 | 324 037 |
Achenbach and Heinecke3 | 0.46 | ⋯ | ⋯ | 0.51 | 0.51 | 273 861 |
So and Savakar4 | 0.50 | ⋯ | 1.18 | 0.50 | 0.49 | 256 905 |
Kawamura et al.5 | 1 | 5.7 | 0.11 | 0.48 | 0.48 | 223 607 |
Kawamura et al.5 | 5 | 0.53 | 1.18 | 0.48 | 0.42 | 146 969 |
So and Savakar4 | 9.5 | ⋯ | 1.18 | 0.23 | 0.34 | 60 000 |
Kawamura et al.5 | 13 | 0.089 | 7.93 | 0.48 | 0.281 | 63 246 |
Blackburn and Melbourne7 | 18 | 0.53 | 2.43 | 0.23 | 0.23 | ⋯ |
Cheung and Melbourne6 | 0.4 | 0.521 | 2.47 | 0.23 | 0.23 | 229 169 |
Cheung and Melbourne6 | 1.6 | 0.521 | 2.47 | 0.23 | 0.25 | 141 421 |
Cheung and Melbourne6 | 4.4 | 0.521 | 2.47 | 0.23 | 0.30 | 109 545 |
Cheung and Melbourne6 | 6.8 | 0.521 | 2.47 | 0.23 | 0.32 | 69 282 |
Cheung and Melbourne6 | 9.1 | 0.521 | 2.47 | 0.23 | 0.34 | 77 460 |
The dimensionless turbulence forcing parameter cR/cL for the two groups of flow tests is adjusted for nominal fits with the data as shown. The Strouhal number behavior for α2≅ 3, commonly associated with extreme surface roughness3 or with localized turbulence,22 is also seen as an asymptotic limit of the Strouhal number data at extreme turbulence,5,7 as in cross-flow heat exchangers.
In the following, we consider three classes of turbulence intensity: low turbulence (TI < 1%), moderate turbulence (TI < 5%–8%), commonly present in straight pipe flow, and extreme turbulence (TI > 10%), found in highly disturbed flows adjacent to piping elbows, heat exchangers, turbines, and combustion processes. The role played by surface roughness will become more apparent when we examine the critical Reynolds number data.
Representative supercritical Strouhal numbers vs TI are shown in Fig. 1 for the two groups of flow tests and are clearly distinguished by their bandwidth parameters at low to moderate TI. At high levels of turbulence, the two datasets appear to merge, followed by decreasing Strouhal numbers with increased turbulence, regardless of the bandwidth parameter. This apparent merger marks the onset of turbulence-dominated forcing. For example, arbitrarily considering a doubling of the relative RMS forcing as marking the onset of turbulence-dominated forcing or , turbulence domination is expected to occur at TI ≅ 10% and TI ≅ 6%, respectively, for the two datasets in Fig. 1. As only subresonant cylinders are tested, the displacement is linearly related to the applied force, and so the forcing and displacement ratios are identical. The RMS error of the data fit for the two datasets is 3%.
The representative supercritical Strouhal numbers shown in Fig. 1 suggest the following. Narrowband turbulence, presently defined as α2 < 3, tends to interfere with the shedding process, with a reduction in the shedding rate with increasing turbulence intensity.5 Broadband turbulence, however, with α2 > 3, tends to energize the shedding process6 at low to moderate turbulence, with the shedding rate increasing with increasing turbulence intensity. At extreme levels of TI, the shedding rate in both cases decreases, approaching that of subcritical Reynolds numbers or highly roughened surfaces.
Figure 2, and also Table I, show the bandwidth parameters using the reported St0 and longitudinal integral scales Λ/d at velocities well into the drag crisis.
At low to moderate turbulence levels, the datasets are separated by the Strouhal bifurcation, but for TI greater than 8%–10%, we find a combination of elevated bandwidth parameters with a turbulence-dominated reduction in Strouhal numbers5,7 for both liquid and atmospheric wind tunnel tests.
Finally, the mean critical Reynolds numbers determined from reported drag force measurements, in the presence of free-stream turbulence, are shown in Fig. 3, with error bars marking the transition associated with the drag crisis.
VI. CONCLUSION
Application of Bendat’s analysis19 to vortex shedding was motivated by a study of thermowell failure in pipe flow,25 where consideration of a range of Reynolds-number-dependent supercritical Strouhal numbers is required for proper risk assessment. While this is hardly a solution of the Navier–Stokes equations, it helps resolve various interpretations of the supercritical Strouhal numbers associated with wind-induced vortex shedding rates6,7,24 of architectural and engineering structures, whether due to turbulence parameters, surface roughness, or combinations thereof.
The use of thermowell calculations based on traditional sinusoidal forcing models is still possible. First, the Reynolds numbers should be checked under maximum flow conditions. If the fitting is exposed to supercritical Reynolds numbers, then two calculations should be performed over the expected flow range, one with and one without elevated Reynolds-number-dependent Strouhal numbers, as in Refs. 17, 25, and 26. Should potential resonance conditions be identified, a reduction in the length of the intrusive fitting will often eliminate the risk of resonance. Such changes have little effect on measurement error or speed of response, as can be calculated analytically or by simulation.
Subject to additional flow tests and simulations that account for fluid turbulence, preliminary evaluations suggest that the proposed model with suitable Reynolds-number-dependent forcing parameters may be used in characterizing the Strouhal number for subcritical Reynolds numbers.27 The theory may also provide additional insight into the application of straked thermowells,28 which exhibit contrasting behavior to that observed in tests with roughened cylinders,3,29 since weakened resonances are observed, with elevated Strouhal numbers marking the postcritical Reynolds transition in high-pressure low-turbulence gas flows.30 The analysis presented here also explains the behavior of vortex flow meters31 and in so doing suggests that errors associated with free-stream turbulence may be corrected with suitable flow conditioners.
AUTHOR DECLARATIONS
Conflict of Interest
The authors have no conflicts to disclose.
Author Contributions
Dave Bartran: Conceptualization (equal); Data curation (equal); Formal analysis (equal); Funding acquisition (equal); Investigation (equal); Methodology (equal); Project administration (equal); Resources (equal); Software (equal); Supervision (equal); Validation (equal); Visualization (equal); Writing – original draft (equal); Writing – review & editing (equal).
DATA AVAILABILITY
The data that supports the findings of this study are available within the article.
NOMENCLATURE
- cL
mean RMS oscillating lift force parameter
- cR
mean RMS turbulence forcing parameter
- d
cylinder diameter
dimensionless forcing and its rate of change
- g(ω)
single-sided bandlimited white noise spectra
- h(ω)
single-sided higher-order monotonic spectra
level-crossing rate (twice the effective frequency)
- Rd
Reynolds number, Vd/ν
- St
Strouhal number, (ω/2π) d/V
- St0
Strouhal number at minimal turbulence
- TI
turbulence intensity,
- V, v′(t)
mean fluid velocity and its fluctuations in time
Greek
- α ≡ ωBW/ω0
turbulence bandwidth parameter
level-crossing conditions for forcing and its rate of change
- Λ
integral scale length of turbulent fluctuations
- ρf
fluid density
variance of lift forcing and of its rate of change
- τΛ
integral time scale of turbulent eddies
- ϕ(t)
uniformly distributed phase fluctuations of sinusoidal forcing
- ω
frequency (rad/s)
- ω0
vortex shedding rate at minimal turbulence
- ωBW
spectral bandwidth of velocity fluctuations
NOMENCLATURE
- cL
mean RMS oscillating lift force parameter
- cR
mean RMS turbulence forcing parameter
- d
cylinder diameter
dimensionless forcing and its rate of change
- g(ω)
single-sided bandlimited white noise spectra
- h(ω)
single-sided higher-order monotonic spectra
level-crossing rate (twice the effective frequency)
- Rd
Reynolds number, Vd/ν
- St
Strouhal number, (ω/2π) d/V
- St0
Strouhal number at minimal turbulence
- TI
turbulence intensity,
- V, v′(t)
mean fluid velocity and its fluctuations in time
Greek
- α ≡ ωBW/ω0
turbulence bandwidth parameter
level-crossing conditions for forcing and its rate of change
- Λ
integral scale length of turbulent fluctuations
- ρf
fluid density
variance of lift forcing and of its rate of change
- τΛ
integral time scale of turbulent eddies
- ϕ(t)
uniformly distributed phase fluctuations of sinusoidal forcing
- ω
frequency (rad/s)
- ω0
vortex shedding rate at minimal turbulence
- ωBW
spectral bandwidth of velocity fluctuations