Numerical investigation of turbulent supersonic axisymmetric wakes

2012 ◽  
Vol 702 ◽  
pp. 488-520 ◽  
Author(s):  
Richard D. Sandberg

AbstractNumerical experiments are conducted of turbulent supersonic axisymmetric wakes at Mach number $M= 2. 46$ and Reynolds number, based on free-stream velocity and base diameter, ${\mathit{Re}}_{D} = 1\ensuremath{\times} 1{0}^{5} $. Direct numerical simulations (DNS) are used to study the effect of approach flow conditions, and of specific azimuthal modes, on the near-wake behaviour. To that end, DNS are performed with laminar and turbulent approach boundary layers, and additional turbulent approach flow DNS with reduced circumferential size are conducted to deliberately eliminate certain azimuthal/helical modes. DNS with turbulent approach flow show an increased turning angle and increased growth of the separating shear layer, leading to a shorter recirculation region, a stronger recompression shock system, and ultimately good agreement with experimental data at considerably higher Reynolds number. A similar wake structure is found for laminar and turbulent inflow conditions, giving further evidence of the wake structure being a consequence of the global near-wake instabilities and not a result of upstream conditions. Stability analyses of two-dimensional basic states are carried out by computing the temporal pulse response using forced Navier–Stokes simulations to investigate which azimuthal modes are dominant for fully turbulent wakes and how the stability behaviour is influenced by the choice of basic state. Using the time- and azimuthally averaged data from three-dimensional DNS with turbulent inflow as basic state, an absolute instability of the axisymmetric mode was found and helical modes $m= 4, 5, 6$ were found to be linearly most unstable, in contrast to results obtained earlier using an axisymmetric flow solution as the basic state. The addition of a turbulence viscosity in the forced DNS retains most of the stability characteristics but reduces the wavenumber of the linearly most-amplified modes.

2005 ◽  
Vol 127 (6) ◽  
pp. 1085-1094 ◽  
Author(s):  
Alan L. Kastengren ◽  
J. Craig Dutton

The near wake of a blunt-base cylinder at 10° angle-of-attack to a Mach 2.46 free-stream flow is visualized at several locations to study unsteady aspects of its structure. In both side-view and end-view images, the shear layer flapping grows monotonically as the shear layer develops, similar to the trends seen in a corresponding axisymmetric supersonic base flow. The interface convolution, a measure of the tortuousness of the shear layer, peaks for side-view and end-view images during recompression. The high convolution for a septum of fluid seen in the middle of the wake indicates that the septum actively entrains fluid from the recirculation region, which helps to explain the low base pressure for this wake compared to that for a corresponding axisymmetric wake.


1991 ◽  
Vol 113 (3) ◽  
pp. 384-398 ◽  
Author(s):  
M. P. Arnal ◽  
D. J. Goering ◽  
J. A. C. Humphrey

The characteristics of the flow around a bluff body of square cross-section in contact with a solid-wall boundary are investigated numerically using a finite difference procedure. Previous studies (Taneda, 1965; Kamemoto et al., 1984) have shown qualitatively the strong influence of solid-wall boundaries on the vortex-shedding process and the formation of the vortex street downstream. In the present study three cases are investigated which correspond to flow past a square rib in a freestream, flow past a rib on a fixed wall and flow past a rib on a sliding wall. Values of the Reynolds number studied ranged from 100 to 2000, where the Reynolds number is based on the rib height, H, and bulk stream velocity, Ub. Comparisons between the sliding-wall and fixed-wall cases show that the sliding wall has a significant destabilizing effect on the recirculation region behind the rib. Results show the onset of unsteadiness at a lower Reynolds number for the sliding-wall case (50 ≤ Recrit ≤100) than for the fixed-wall case (Recrit≥100). A careful examination of the vortex-shedding process reveals similarities between the sliding-wall case and both the freestream and fixed-wall cases. At moderate Reynolds numbers (Re≥250) the sliding-wall results show that the rib periodically sheds vortices of alternating circulation in much the same manner as the rib in a freestream; as in, for example, Davis and Moore [1982]. The vortices are distributed asymmetrically downstream of the rib and are not of equal strength as in the freestream case. However, the sliding-wall case shows no tendency to develop cycle-to-cycle variations at higher Reynolds numbers, as observed in the freestream and fixed-wall cases. Thus, while the moving wall causes the flow past the rib to become unsteady at a lower Reynolds number than in the fixed-wall case, it also acts to stabilize or “lock-in” the vortex-shedding frequency. This is attributed to the additional source of positive vorticity immediately downstream of the rib on the sliding wall.


2009 ◽  
Vol 643 ◽  
pp. 349-362 ◽  
Author(s):  
DAVID LO JACONO ◽  
JUSTIN S. LEONTINI ◽  
MARK C. THOMPSON ◽  
JOHN SHERIDAN

A study of the flow past an oscillatory rotating cylinder has been conducted, where the frequency of oscillation has been matched to the natural frequency of the vortex street generated in the wake of a stationary cylinder, at Reynolds number 300. The focus is on the wake transition to three-dimensional flow and, in particular, the changes induced in this transition by the addition of the oscillatory rotation. Using Floquet stability analysis, it is found that the fine-scale three-dimensional mode that typically dominates the wake at a Reynolds number beyond that at the second transition to three-dimensional flow (referred to as mode B) is suppressed for amplitudes of rotation beyond a critical amplitude, in agreement with past studies. However, the rotation does not suppress the development of three-dimensionality completely, as other modes are discovered that would lead to three-dimensional flow. In particular, the longer-wavelength mode that leads the three-dimensional transition in the wake of a stationary cylinder (referred to as mode A) is left essentially unaffected at low amplitudes of rotation. At higher amplitudes of oscillation, mode A is also suppressed as the two-dimensional near wake changes in character from a single- to a double-row wake; however, another mode is predicted to render the flow three-dimensional, dubbed mode D (for double row). This mode has the same spatio-temporal symmetries as mode A.


2011 ◽  
Vol 681 ◽  
pp. 411-433 ◽  
Author(s):  
HEMANT K. CHAURASIA ◽  
MARK C. THOMPSON

A detailed numerical study of the separating and reattaching flow over a square leading-edge plate is presented, examining the instability modes governing transition from two- to three-dimensional flow. Under the influence of background noise, experiments show that the transition scenario typically is incompletely described by either global stability analysis or the transient growth of dominant optimal perturbation modes. Instead two-dimensional transition effectively can be triggered by the convective Kelvin–Helmholtz (KH) shear-layer instability; although it may be possible that this could be described alternatively in terms of higher-order optimal perturbation modes. At least in some experiments, observed transition occurs by either: (i) KH vortices shedding downstream directly and then almost immediately undergoing three-dimensional transition or (ii) at higher Reynolds numbers, larger vortical structures are shed that are also three-dimensionally unstable. These two paths lead to distinctly different three-dimensional arrangements of vortical flow structures. This paper focuses on the mechanisms underlying these three-dimensional transitions. Floquet analysis of weakly periodically forced flow, mimicking the observed two-dimensional quasi-periodic base flow, indicates that the two-dimensional vortex rollers shed from the recirculation region become globally three-dimensionally unstable at a Reynolds number of approximately 380. This transition Reynolds number and the predicted wavelength and flow symmetries match well with those of the experiments. The instability appears to be elliptical in nature with the perturbation field mainly restricted to the cores of the shed rollers and showing the spatial vorticity distribution expected for that instability type. Indeed an estimate of the theoretical predicted wavelength is also a good match to the prediction from Floquet analysis and theoretical estimates indicate the growth rate is positive. Fully three-dimensional simulations are also undertaken to explore the nonlinear development of the three-dimensional instability. These show the development of the characteristic upright hairpins observed in the experimental dye visualisations. The three-dimensional instability that manifests at lower Reynolds numbers is shown to be consistent with an elliptic instability of the KH shear-layer vortices in both symmetry and spanwise wavelength.


2017 ◽  
Vol 822 ◽  
pp. 813-847 ◽  
Author(s):  
Azan M. Sapardi ◽  
Wisam K. Hussam ◽  
Alban Pothérat ◽  
Gregory J. Sheard

This study seeks to characterise the breakdown of the steady two-dimensional solution in the flow around a 180-degree sharp bend to infinitesimal three-dimensional disturbances using a linear stability analysis. The stability analysis predicts that three-dimensional transition is via a synchronous instability of the steady flows. A highly accurate global linear stability analysis of the flow was conducted with Reynolds number $\mathit{Re}<1150$ and bend opening ratio (ratio of bend width to inlet height) $0.2\leqslant \unicode[STIX]{x1D6FD}\leqslant 5$. This range of $\mathit{Re}$ and $\unicode[STIX]{x1D6FD}$ captures both steady-state two-dimensional flow solutions and the inception of unsteady two-dimensional flow. For $0.2\leqslant \unicode[STIX]{x1D6FD}\leqslant 1$, the two-dimensional base flow transitions from steady to unsteady at higher Reynolds number as $\unicode[STIX]{x1D6FD}$ increases. The stability analysis shows that at the onset of instability, the base flow becomes three-dimensionally unstable in two different modes, namely a spanwise oscillating mode for $\unicode[STIX]{x1D6FD}=0.2$ and a spanwise synchronous mode for $\unicode[STIX]{x1D6FD}\geqslant 0.3$. The critical Reynolds number and the spanwise wavelength of perturbations increase as $\unicode[STIX]{x1D6FD}$ increases. For $1<\unicode[STIX]{x1D6FD}\leqslant 2$ both the critical Reynolds number for onset of unsteadiness and the spanwise wavelength decrease as $\unicode[STIX]{x1D6FD}$ increases. Finally, for $2<\unicode[STIX]{x1D6FD}\leqslant 5$, the critical Reynolds number and spanwise wavelength remain almost constant. The linear stability analysis also shows that the base flow becomes unstable to different three-dimensional modes depending on the opening ratio. The modes are found to be localised near the reattachment point of the first recirculation bubble.


1992 ◽  
Vol 238 ◽  
pp. 1-30 ◽  
Author(s):  
George Em Karniadakis ◽  
George S. Triantafyllou

The wakes of bluff objects and in particular of circular cylinders are known to undergo a ‘fast’ transition, from a laminar two-dimensional state at Reynolds number 200 to a turbulent state at Reynolds number 400. The process has been documented in several experimental investigations, but the underlying physical mechanisms have remained largely unknown so far. In this paper, the transition process is investigated numerically, through direct simulation of the Navier—Stokes equations at representative Reynolds numbers, up to 500. A high-order time-accurate, mixed spectral/spectral element technique is used. It is shown that the wake first becomes three-dimensional, as a result of a secondary instability of the two-dimensional vortex street. This secondary instability appears at a Reynolds number close to 200. For slightly supercritical Reynolds numbers, a harmonic state develops, in which the flow oscillates at its fundamental frequency (Strouhal number) around a spanwise modulated time-average flow. In the near wake the modulation wavelength of the time-average flow is half of the spanwise wavelength of the perturbation flow, consistently with linear instability theory. The vortex filaments have a spanwise wavy shape in the near wake, and form rib-like structures further downstream. At higher Reynolds numbers the three-dimensional flow oscillation undergoes a period-doubling bifurcation, in which the flow alternates between two different states. Phase-space analysis of the flow shows that the basic limit cycle has branched into two connected limit cycles. In physical space the period doubling appears as the shedding of two distinct types of vortex filaments.Further increases of the Reynolds number result in a cascade of period-doubling bifurcations, which create a chaotic state in the flow at a Reynolds number of about 500. The flow is characterized by broadband power spectra, and the appearance of intermittent phenomena. It is concluded that the wake undergoes transition to turbulence following the period-doubling route.


2001 ◽  
Vol 438 ◽  
pp. 363-377 ◽  
Author(s):  
A. Yu. GELFGAT ◽  
P. Z. BAR-YOSEPH ◽  
A. SOLAN

The axisymmetry-breaking three-dimensional instability of the axisymmetric flow between a rotating lid and a stationary cylinder is analysed. The flow is governed by two parameters – the Reynolds number Re and the aspect ratio γ (=height/radius). Published experimental results indicate that in different ranges of γ axisymmetric or non-axisymmetric instabilities can be observed. Previous analyses considered only axisymmetric instability. The present analysis is devoted to the linear stability of the basic axisymmetric flow with respect to the non-axisymmetric perturbations. After the linearization the stability problem separates into a family of quasi-axisymmetric subproblems for discrete values of the azimuthal wavenumber k. The computations are done using the global Galerkin method. The stability analysis is carried out at various densely distributed values of γ in the range 1 < γ < 3.5. It is shown that the axisymmetric perturbations are dominant in the range 1.63 < γ < 2.76. Outside this range, for γ < 1.63 and for γ > 2.76, the instability is three-dimensional and sets in with k = 2 and k = 3 or 4, respectively. The azimuthal periodicity, patterns, characteristic frequencies and phase velocities of the dominant perturbations are discussed.


Author(s):  
Ryoji Kojima ◽  
Donghi Lee ◽  
Tomoaki Tatsukawa ◽  
Taku Nonomura ◽  
Akira Oyama ◽  
...  

The effects of aspect ratio and Reynolds number on aerodynamic characteristics of three-dimensional rectangular wing at low Reynolds number of 103 to 105, are investigated with Reynolds-averaged Navier-Stokes solver with the Baldwin-Lomax model. Present results show that lift coefficient decreases drastically at lower aspect ratio than 4. Besides, the much larger viscous drag coefficient is obtained at the lower Reynolds number, especially lower than 104. In order to focus on designing practical wings, the particular cases under the condition of fixed wing-surface area and fixed main stream velocity are conducted. The results show that there is trade-off between the decrease in viscous drag coefficient with increasing Reynolds number and the increase in lift coefficient with increasing aspect ratio. At the lower Reynolds number condition, as the former effect is stronger than the latter one, maximum lift-to-drag ratio is obtained at lower aspect ratio.


Author(s):  
A. Inasawa ◽  
K. Toda ◽  
M. Asai

Disturbance growth in the wake of a circular cylinder moving at a constant acceleration is examined experimentally. The cylinder is installed on a carriage moving in the still air. The results show that the critical Reynolds number for the onset of the global instability leading to a self-sustained wake oscillation increases with the magnitude of acceleration, while the Strouhal number of the growing disturbance at the critical Reynolds number is not strongly dependent on the magnitude of acceleration. It is also found that with increasing the acceleration, the Ka´rma´n vortex street remains two-dimensional even at the Reynolds numbers around 200 where the three-dimensional instability occurs to lead to the vortex dislocation in the case of cylinder moving at constant velocity or in the case of cylinder wake in the steady oncoming flow.


Sign in / Sign up

Export Citation Format

Share Document