1. Introduction
The flow around a circular cylinder is encountered in many engineering applications as well as in nature and in fundamental research. Some examples are the air flow around a cooling tower, as well as missiles and aircrafts in the transonic regime. In transonic flows and fluid dynamical situations in general, a variety of typical phenomena, such as shock waves, vortices, boundary layers, flow separation and shear layers, arise. Studying the interaction of those phenomena is of large interest in fluid–structure interactions, because resonance frequencies can be excited by flow instabilities. The question arises as to what frequency is triggered by the flow around a cylinder, which could provoke an oscillation. In dimensionless terms, this means that the Strouhal number becomes a function not only of the Reynolds number , but also of the Mach number . In practical situations embedded in a larger project context, where questions like the mentioned one arise, one often does not have the time nor the resources to develop one’s own numerical tool for simulation and calculation. Rather, one has to refer to methods and tools available on the market; that is why the investigation described here is based on a commercially available CFD tool. It is clear, though, that such a tool has to be checked appropriately before the corresponding findings can be used.
There are various studies describing the flow phenomena occurring in incompressible flow. For small Reynolds numbers in the range of
, a Kármán vortex street is formed in the wake of the cylinder described by Schlichting and Gersten [
1]. Although the flow around a circular cylinder has received a lot of attention in recent decades, very few experimental and numerical studies are available about the compressible and in particular transonic flow around a circular cylinder. One of the reasons might be the difficulties in performing experiments of simple flow situations such as the planar flow around a circular cylinder. However, planar aerodynamical investigations are also of interest, as they allow one, for example, to estimate in a simplified manner the force on a body in a surrounding transonic flow. The most important experimental studies and their findings are summarised below.
Macha [
2] performed wind tunnel tests in order to determine the drag coefficient
for a Reynolds number range of
and a Mach number range of
. One of the most important findings is the reduction in
from
to
, caused by the formation of shock waves. In addition, Murthy and Rose [
3] performed a series of wind tunnel tests with
and
. The increase in
, as
reaches sonic conditions and agrees with the findings of Macha [
2]. Furthermore, it was found that the detectable vortex shedding ceases at
. The ranges
and
were investigated by Rodriguez [
4] using a wind tunnel. It was found that the coupling between the near wake and the vortex street increases with increasing
. As soon as local regions of the flow reach sonic conditions and
shocks occur, the coupling between the vortex street and the near wake is cut off. The upstream flow field is now independent of the vortex street. In addition, the Strouhal number
is approximately 0.2, except for a rise when the quasi-steady regime is reached. In addition, the drag and lift coefficients,
and
, were calculated from the pressure measurements. Ackerman et al. [
5] experimentally investigated time-resolved pressure distributions at
and
. From these measurements, the surface pressure fluctuations,
,
and the occurring flow regimes were evaluated. For
, local regions of flow around the cylinder reach sonic conditions, but only on one side of the cylinder at a time. The flow enters the intermittent shock wave regime. As
increases beyond 0.4,
increases. The region downstream of the cylinder, in which the vortices are formed, shortens. Beyond around
, the flow enters the permanent shock wave regime and
decreases. Once the flow enters the wake shock wave regime below
, the vortex formation region becomes elongated. A normal shock grows at the point of vortex roll up and
increases. Nagata et al. [
6] used a low-density wind tunnel with time-resolved Schlieren visualisations, pressure and force measurements, in order to characterise the flow for
and
. The trend of the
effect on the flow field,
and the maximum width of the recirculation change at approximately
.
increases as
increases and the increment becomes larger as
increases. For
,
is independent of
. For
and
at
,
decreases and increases, respectively. Furthermore, it is observed that
increases as
or
increase. Gowen and Perkins [
7] measured the pressure distribution around a circular cylinder in subsonic and supersonic flows and calculated
for
and
. It is shown that
is not influenced by
under the supersonic conditions investigated.
In addition to the experimental investigations, numerical simulations of the compressible flow around a circular cylinder have been increasingly carried out over the last few decades. Some examples of numerical investigations are described below. Botta [
8] integrated the Euler equations numerically to investigate the inviscid flow for
, that means for a Reynolds number
. The time-dependent
and
are evaluated in order to determine
. Furthermore, the distributions of the vorticity, the entropy deviation, pressure coefficient, as well as the velocity fields are provided. Two transitions over the investigated
range were observed, the transition to a chaotic turbulent regime and from this to a quasi-steady flow. In the range
, the solution shows a periodic behaviour. Bobenrieth Miserda and Leal [
9] performed numerical Detached Eddy Simulations of the unsteady transonic flow at
and
500,000, where several complex viscous shock interactions were observed. The frequency of
corresponds to the vortex-shedding frequency, whereas the frequency of the
characterises the viscous shock interaction. In addition, Xu et al. [
10] performed Detached Eddy Simulations for
and various Mach numbers
. Two flow states are found, an unsteady one for
and a quasi-steady flow state for
. The unsteady flow state is characterised by the interaction of moving shock waves, the turbulent boundary layer on the cylinder wall and the vortex shedding in the cylinders near the wake. In the quasi-steady flow state, strong oblique shock waves are formed and the vortex shedding is suppressed. Furthermore, the local supersonic zone, the separation angle and
are evaluated and analysed. Hong et al. [
11] studied
and
using constrained Large Eddy Simulations. The effects of
on the flow patterns and state variables such as the pressure, the skin friction,
and the cylinder surface temperature are studied. Non-monotonic behaviour of the pressure and skin friction distributions are observed with increasing
. The minimum mean separation angle occurs at
. Canuto and Taira [
12] performed Direct Numerical Simulations of
and
. The wake is characterised using different lengths and
, and
and some examples of the pressure distribution are provided. Furthermore, a stability analysis is performed. It is shown that
increases and
decreases with increasing
for constant
. Xia et al. [
13] performed constrained Large Eddy Simulations for
and
and various Mach numbers of
. The separation angle,
, the pressure distribution and the skin friction coefficient were evaluated and analysed. Furthermore, the density gradient
was used to identify four different flow regimes. Shirani [
14] simulated
and
solving the two-dimensional time averaged Navier–Stokes equations numerically. The behaviour of the time averages of
and
and their fluctuation frequencies were evaluated. Matar et al. [
15] investigated the real gas flow around a circular cylinder at high Reynolds numbers
and Mach numbers
between 0.7 and 0.9 using wall-resolved implicit Large Eddy Simulations (iLESs). For experimental validation at
, Background Oriented Schlieren (BOS) visualisations are used. The flow phenomena of a Kármán vortex street, acoustic waves and compression waves are observed. In addition, the Strouhal number
, the wall pressure and the mean pressure drag coefficients are provided and are compared with literature and URANS simulation results. In addition, Linn and Awruch [
16] performed Large Eddy Simulations (LESs) of the two- and three-dimensional flow around a circular cylinder at
500,000 and
using various tetrahedron-adapted meshes. The density gradient
, the streamlines and Q-criterion isosurfaces are used to visualise the flow behaviour. Moreover, the drag and the lift coefficients
and
, the Strouhal number
, the mean surface pressure coefficient and the angle of the boundary layer separation point are provided.
The present study describes the transonic planar flow around a circular cylinder at conditions
and
where only a few investigations have been carried out (
Figure 1). The investigation of the planar situation enables one to analyse the phenomena uncoupled from the influence of potential three-dimensional effects. This investigation gives an overview of the flow phenomena occurring in a wide range of Mach numbers of
to 2 and Reynolds numbers of
to 80,000. Within the first few chapters, the fundamentals and the numerical implementation are briefly introduced. The simulations of the compressible flow around a circular cylinder are verified by applying the code used to the flow in a Laval nozzle. For validation, the results obtained for the flow around a cylinder are compared to the results from other authors. Different regimes with phenomena such as shock waves, sound waves, flow separation, vortex shedding, shear layers and tangential discontinuities are identified and some of them are analysed in more detail. To capture the different flow phenomena, different regions of interest are used for the numerical procedure. This investigation focuses on the behaviour in the wake of the cylinder. In addition, the critical Mach number
is evaluated and the averaged drag coefficient
and the Strouhal number
are provided. Finally, polar diagrams are presented, that means the time-resolved drag and lift coefficients,
and
, are plotted as a function of time in a
-
-diagram, to analyse their phase shift and their frequency ratio.