scholarly journals Three-dimensional instabilities in compressible flow over open cavities

2008 ◽  
Vol 599 ◽  
pp. 309-339 ◽  
Author(s):  
GUILLAUME A. BRÈS ◽  
TIM COLONIUS

Direct numerical simulations are performed to investigate the three-dimensional stability of compressible flow over open cavities. A linear stability analysis is conducted to search for three-dimensional global instabilities of the two-dimensional mean flow for cavities that are homogeneous in the spanwise direction. The presence of such instabilities is reported for a range of flow conditions and cavity aspect ratios. For cavities of aspect ratio (length to depth) of 2 and 4, the three-dimensional mode has a spanwise wavelength of approximately one cavity depth and oscillates with a frequency about one order of magnitude lower than two-dimensional Rossiter (flow/acoustics) instabilities. A steady mode of smaller spanwise wavelength is also identified for square cavities. The linear results indicate that the instability is hydrodynamic (rather than acoustic) in nature and arises from a generic centrifugal instability mechanism associated with the mean recirculating vortical flow in the downstream part of the cavity. These three-dimensional instabilities are related to centrifugal instabilities previously reported in flows over backward-facing steps, lid-driven cavity flows and Couette flows. Results from three-dimensional simulations of the nonlinear compressible Navier–Stokes equations are also reported. The formation of oscillating (and, in some cases, steady) spanwise structures is observed inside the cavity. The spanwise wavelength and oscillation frequency of these structures agree with the linear analysis predictions. When present, the shear-layer (Rossiter) oscillations experience a low-frequency modulation that arises from nonlinear interactions with the three-dimensional mode. The results are consistent with observations of low-frequency modulations and spanwise structures in previous experimental and numerical studies on open cavity flows.

1987 ◽  
Vol 184 ◽  
pp. 207-243 ◽  
Author(s):  
Ralph W. Metcalfe ◽  
Steven A. Orszag ◽  
Marc E. Brachet ◽  
Suresh Menon ◽  
James J. Riley

The three-dimensional stability of two-dimensional vortical states of planar mixing layers is studied by direct numerical integration of the Navier-Stokes equations. Small-scale instabilities are shown to exist for spanwise scales at which classical linear modes are stable. These modes grow on convective timescales, extract their energy from the mean flow and exist at moderately low Reynolds numbers. Their growth rates are comparable with the most rapidly growing inviscid instability and with the growth rates of two-dimensional subharmonic (pairing) modes. At high amplitudes, they can evolve into pairs of counter-rotating, streamwise vortices, connecting the primary spanwise vortices, which are very similar to the structures observed in laboratory experiments. The three-dimensional modes do not appear to saturate in quasi-steady states as do the purely two-dimensional fundamental and subharmonic modes in the absence of pairing. The subsequent evolution of the flow depends on the relative amplitudes of the pairing modes. Persistent pairings can inhibit three-dimensional instability and, hence, keep the flow predominantly two-dimensional. Conversely, suppression of the pairing process can drive the three-dimensional modes to more chaotic, turbulent-like states. An analysis of high-resolution simulations of fully turbulent mixing layers confirms the existence of rib-like structures and that their coherence depends strongly on the presence of the two-dimensional pairing modes.


2002 ◽  
Vol 450 ◽  
pp. 1-33 ◽  
Author(s):  
HERMANN F. FASEL

Direct numerical simulations (DNS) of the Navier–Stokes equations are used to investigate the role of the Klebanoff-mode in laminar–turbulent transition in a flatplate boundary layer. To model the effects of free-stream turbulence, volume forces are used to generate low-frequency streamwise vortices outside the boundary layer. A suction/blowing slot at the wall is used to generate a two-dimensional Tollmien–Schlichting (TS) wave inside the boundary layer. The characteristics of the fluctuations inside the boundary layer agree very well with those measured in experiments. It is shown how the interaction of the Klebanoff-mode with the two-dimensional TS-wave leads to the formation of three-dimensional TS-wavepackets. When the disturbance amplitudes reach a critical level, a fundamental resonance-type secondary instability causes the breakdown of the TS-wavepackets into turbulent spots.


1990 ◽  
Vol 217 ◽  
pp. 367-390 ◽  
Author(s):  
Philip Hall ◽  
Sharon O. Seddougui

The three-dimensional boundary layer on a swept wing can support different types of hydrodynamic instability. Here attention is focused on the so-called ‘spanwise instability’ problem which occurs when the attachment-line boundary layer on the leading edge becomes unstable to Tollmien–Schlichting waves. In order to gain insight into the interactions that are important in that problem a simplified basic state is considered. This simplified flow corresponds to the swept attachment-line boundary layer on an infinite flat plate. The basic flow here is an exact solution of the Navier–Stokes equations and its stability to two-dimensional waves propagating along the attachment line can be considered exactly at finite Reynolds number. This has been done in the linear and weakly nonlinear regimes by Hall, Malik & Poll (1984) and Hall & Malik (1986). Here the corresponding problem is studied for oblique waves and their interaction with two-dimensional waves is investigated. In fact oblique modes cannot be described exactly at finite Reynolds number so it is necessary to make a high-Reynolds-number approximation and use triple-deck theory. It is shown that there are two types of oblique wave which, if excited, cause the destabilization of the two-dimensional mode and the breakdown of the disturbed flow at a finite distance from the leading edge. First a low-frequency mode closely related to the viscous stationary crossflow mode discussed by Hall (1986) and MacKerrell (1987) is a possible cause of breakdown. Secondly a class of oblique wave with frequency comparable with that of the two-dimensional mode is another cause of breakdown. It is shown that the relative importance of the modes depends on the distance from the attachment line.


2021 ◽  
Vol 0 (0) ◽  
Author(s):  
Niklas Ericsson

Abstract We develop a framework for solving the stationary, incompressible Stokes equations in an axisymmetric domain. By means of Fourier expansion with respect to the angular variable, the three-dimensional Stokes problem is reduced to an equivalent, countable family of decoupled two-dimensional problems. By using decomposition of three-dimensional Sobolev norms, we derive natural variational spaces for the two-dimensional problems, and show that the variational formulations are well-posed. We analyze the error due to Fourier truncation and conclude that, for data that are sufficiently regular, it suffices to solve a small number of two-dimensional problems.


1995 ◽  
Vol 291 ◽  
pp. 369-392 ◽  
Author(s):  
Ronald D. Joslin

The spatial evolution of three-dimensional disturbances in an attachment-line boundary layer is computed by direct numerical simulation of the unsteady, incompressible Navier–Stokes equations. Disturbances are introduced into the boundary layer by harmonic sources that involve unsteady suction and blowing through the wall. Various harmonic-source generators are implemented on or near the attachment line, and the disturbance evolutions are compared. Previous two-dimensional simulation results and nonparallel theory are compared with the present results. The three-dimensional simulation results for disturbances with quasi-two-dimensional features indicate growth rates of only a few percent larger than pure two-dimensional results; however, the results are close enough to enable the use of the more computationally efficient, two-dimensional approach. However, true three-dimensional disturbances are more likely in practice and are more stable than two-dimensional disturbances. Disturbances generated off (but near) the attachment line spread both away from and toward the attachment line as they evolve. The evolution pattern is comparable to wave packets in flat-plate boundary-layer flows. Suction stabilizes the quasi-two-dimensional attachment-line instabilities, and blowing destabilizes these instabilities; these results qualitatively agree with the theory. Furthermore, suction stabilizes the disturbances that develop off the attachment line. Clearly, disturbances that are generated near the attachment line can supply energy to attachment-line instabilities, but suction can be used to stabilize these instabilities.


Author(s):  
Susanne Horn ◽  
Peter J. Schmid ◽  
Jonathan M. Aurnou

Abstract The large-scale circulation (LSC) is the most fundamental turbulent coherent flow structure in Rayleigh-B\'enard convection. Further, LSCs provide the foundation upon which superstructures, the largest observable features in convective systems, are formed. In confined cylindrical geometries with diameter-to-height aspect ratios of Γ ≅ 1, LSC dynamics are known to be governed by a quasi-two-dimensional, coupled horizontal sloshing and torsional (ST) oscillatory mode. In contrast, in Γ ≥ √2 cylinders, a three-dimensional jump rope vortex (JRV) motion dominates the LSC dynamics. Here, we use dynamic mode decomposition (DMD) on direct numerical simulation data of liquid metal to show that both types of modes co-exist in Γ = 1 and Γ = 2 cylinders but with opposite dynamical importance. Furthermore, with this analysis, we demonstrate that ST oscillations originate from a tilted elliptical mean flow superposed with a symmetric higher order mode, which is connected to the four rolls in the plane perpendicular to the LSC in Γ = 1 tanks.


2018 ◽  
Vol 2018 ◽  
pp. 1-16 ◽  
Author(s):  
Fabio Nardecchia ◽  
Annalisa Di Bernardino ◽  
Francesca Pagliaro ◽  
Paolo Monti ◽  
Giovanni Leuzzi ◽  
...  

Computational fluid dynamics (CFD) is currently used in the environmental field to simulate flow and dispersion of pollutants around buildings. However, the closure assumptions of the turbulence usually employed in CFD codes are not always physically based and adequate for all the flow regimes relating to practical applications. The starting point of this work is the performance assessment of the V2F (i.e., v2¯ − f) model implemented in Ansys Fluent for simulating the flow field in an idealized array of two-dimensional canyons. The V2F model has been used in the past to predict low-speed and wall-bounded flows, but it has never been used to simulate airflows in urban street canyons. The numerical results are validated against experimental data collected in the water channel and compared with other turbulence models incorporated in Ansys Fluent (i.e., variations of both k-ε and k-ω models and the Reynolds stress model). The results show that the V2F model provides the best prediction of the flow field for two flow regimes commonly found in urban canopies. The V2F model is also employed to quantify the air-exchange rate (ACH) for a series of two-dimensional building arrangements, such as step-up and step-down configurations, having different aspect ratios and relative heights of the buildings. The results show a clear dependence of the ACH on the latter two parameters and highlight the role played by the turbulence in the exchange of air mass, particularly important for the step-down configurations, when the ventilation associated with the mean flow is generally poor.


1989 ◽  
Vol 199 ◽  
pp. 403-440 ◽  
Author(s):  
E. Laurien ◽  
L. Kleiser

The laminar-turbulent transition process in a parallel boundary-layer with Blasius profile is simulated by numerical integration of the three-dimensional incompressible Navier-Stokes equations using a spectral method. The model of spatially periodic disturbances developing in time is used. Both the classical Klebanoff-type and the subharmonic type of transition are simulated. Maps of the three-dimensional velocity and vorticity fields and visualizations by integrated fluid markers are obtained. The numerical results are compared with experimental measurements and flow visualizations by other authors. Good qualitative and quantitative agreement is found at corresponding stages of development up to the one-spike stage. After the appearance of two-dimensional Tollmien-Schlichting waves of sufficiently large amplitude an increasing three-dimensionality is observed. In particular, a peak-valley structure of the velocity fluctuations, mean longitudinal vortices and sharp spike-like instantaneous velocity signals are formed. The flow field is dominated by a three-dimensional horseshoe vortex system connected with free high-shear layers. Visualizations by time-lines show the formation of A-structures. Our numerical results connect various observations obtained with different experimental techniques. The initial three-dimensional steps of the transition process are consistent with the linear theory of secondary instability. In the later stages nonlinear interactions of the disturbance modes and the production of higher harmonics are essential.We also study the control of transition by local two-dimensional suction and blowing at the wall. It is shown that transition can be delayed or accelerated by superposing disturbances which are out of phase or in phase with oncoming Tollmien-Schlichting instability waves, respectively. Control is only effective if applied at an early, two-dimensional stage of transition. Mean longitudinal vortices remain even after successful control of the fluctuations.


2014 ◽  
Vol 16 (5) ◽  
pp. 901-918 ◽  

<div> <p>Three-dimensional calculations were performed to simulate the flow around a cylindrical vegetation element using the Scale Adaptive Simulation (SAS) model; commonly, this is the first step of the modeling of the flow through multiple vegetation elements. SAS solves the Reynolds Averaged Navier-Stokes equations in stable flow regions, while in regions with unstable flow it goes unsteady producing a resolved turbulent spectrum after reducing eddy viscosity according to the locally resolved vortex size represented by the von Karman length scale. A finite volume numerical code was used for the spatial discretisation of the rectangular computational domain with stream-wise, cross-flow and vertical dimensions equal to 30D, 11D and 1D, respectively, which was resolved with unstructured grids. Calculations were compared with experiments and Large Eddy Simulations (LES). Predicted overall flow parameters and mean flow velocities exhibited a very satisfactory agreement with experiments and LES, while the agreement of predicted turbulent stresses was satisfactory. Calculations showed that SAS is an efficient and relatively fast turbulence modeling approach, especially in relevant practical problems, in which the very high accuracy that can be achieved by LES at the expense of large computational times is not required.</p> </div> <p>&nbsp;</p>


2008 ◽  
Vol 614 ◽  
pp. 315-327 ◽  
Author(s):  
UWE EHRENSTEIN ◽  
FRANÇOIS GALLAIRE

A separated boundary-layer flow at the rear of a bump is considered. Two-dimensional equilibrium stationary states of the Navier–Stokes equations are determined using a nonlinear continuation procedure varying the bump height as well as the Reynolds number. A global instability analysis of the steady states is performed by computing two-dimensional temporal modes. The onset of instability is shown to be characterized by a family of modes with localized structures around the reattachment point becoming almost simultaneously unstable. The optimal perturbation analysis, by projecting the initial disturbance on the set of temporal eigenmodes, reveals that the non-normal modes are able to describe localized initial perturbations associated with the large transient energy growth. At larger time a global low-frequency oscillation is found, accompanied by a periodic regeneration of the flow perturbation inside the bubble, as the consequence of non-normal cancellation of modes. The initial condition provided by the optimal perturbation analysis is applied to Navier–Stokes time integration and is shown to trigger the nonlinear ‘flapping’ typical of separation bubbles. It is possible to follow the stationary equilibrium state on increasing the Reynolds number far beyond instability, ruling out for the present flow case the hypothesis of some authors that topological flow changes are responsible for the ‘flapping’.


Sign in / Sign up

Export Citation Format

Share Document