Three-dimensional instability and vorticity patterns in the wake of a flat plate

International audienceWe investigated experimentally the dynamics of the three-dimensional secondary instability developing in the wake of a thin flat plate at moderate Reynolds numbers. The wake is formed as the two laminar boundary layers developing on each side merge at the trailing edge of the flat plate. Both the spatial and temporal evolution of the two- and three-dimensional instabilities are analysed by means of laser-induced visualizations of the deformation of the interface separating the two streams. It was found that although the wake may exhibit two distinct three-dimensional modes with different symmetry characteristics, Modes 1 and 2 (Lasheras & Meiburg 1990), the latter appears to be amplified first, thereafter dominating the evolution of the near wake. By varying the forcing frequency of the primary two-dimensional instability, we found that the wavelength of the three-dimensional mode is selected by the wavelength of the two-dimensional Karman vortices, with a ratio (lambda(3D)/lambda(2D)) of order one. In the far-wake region, both modes appear to grow and co-exist. Furthermore, by analysing the response of the wake to spanwise-periodic and impulsive perturbations applied at the trailing edge of the plate, we demonstrate that the nature of the secondary instability of the wake behind a thin flat plate is convective. In addition, both modes are shown to have comparable wavelengths and to be the result of the same instability mechanism

Development and stability of gyrotactic plumes in bioconvection

Using the continuum model of Pedley, Hill and Kessler (1988) for bioconvection in a suspension of swimming, gyrotactic micro-organisms, we investigate the existence and stability of a two-dimensional plume in tall, narrow chambers with stress-free sidewalls. The system is governed by the Navier–Stokes equations for an incompressible fluid coupled with a micro-organism conservation equation. These equations are solved numerically using a conservative finite-difference scheme. In sufficiently deep chambers, the plume is always unstable to both varicose and meandering modes. A linear stability analysis for an infinitely long plume predicts the growth rates of these instabilities, explains the mechanisms, and is in good agreement with the numerical results

Global stability analysis and direct numerical simulation of boundary layers with an isolated roughness element

Global stability analysis and direct numerical simulation (DNS) are performed to study boundary layer flows with an isolated roughness element. Wall-attached cuboids with aspect ratios $\eta=1$ and $\eta=0.5$ are investigated for fixed ratio of roughness height to displacement boundary layer thickness $h/\delta^*=2.86$. Global stability analysis is able to capture the frequency of the primary vortical structures. For $\eta=1$, only varicose instability is seen. For the thinner roughness element ($\eta=0.5$), the varicose instability dominates the sinuous instability, and the sinuous instability becomes more pronounced as $Re_h$ increases, due to increased spanwise shear in the near-wake region. The unstable modes mainly extract energy from the central streak, although the lateral streaks also contribute. The DNS results show that different instability features lead to different behavior and development of vortical structures in the nonlinear transition process. For $\eta=1$, the varicose mode is associated with the shedding of hairpin vortices. As $Re_h$ increases, the breakdown of hairpin vortices occurs closer to the roughness and sinuous breakdown behavior promoting transition to turbulence is seen in the farther wake. A fully-developed turbulent flow is established in both the inner and outer layers farther downstream when $Re_h$ is sufficiently high. For $\eta=0.5$, the sinuous wiggling of hairpin vortices is prominent at higher $Re_h$, leading to stronger interactions in the near wake, as a result of combined varicose and sinuous instabilities. A sinuous mode captured by dynamic mode decomposition (DMD) analysis, and associated with the `wiggling' of streaks persists far downstream

Nonlinear input/output analysis: application to boundary layer transition

We extend linear input/output (resolvent) analysis to take into account nonlinear triadic interactions by considering a finite number of harmonics in the frequency domain using the harmonic balance method. Forcing mechanisms that maximise the drag are calculated using a gradient-based ascent algorithm. By including nonlinearity in the analysis, the proposed frequency-domain framework identifies the worst-case disturbances for laminar-turbulent transition. We demonstrate the framework on a flat-plate boundary layer by considering three-dimensional spanwise-periodic perturbations triggered by a few optimal forcing modes of finite amplitude. Two types of volumetric forcing are considered, one corresponding to a single frequency/spanwise wavenumber pair, and a multi-harmonic where a harmonic frequency and wavenumber are also added. Depending on the forcing strategy, we recover a range of transition scenarios associated with K-type and H-type mechanisms, including oblique and planar Tollmien–Schlichting waves, streaks and their breakdown. We show that nonlinearity plays a critical role in optimising growth by combining and redistributing energy between the linear mechanisms and the higher perturbation harmonics. With a very limited range of frequencies and wavenumbers, the calculations appear to reach the early stages of the turbulent regime through the generation and breakdown of hairpin and quasi-streamwise staggered vortices

Rayleigh-Benard Convection with a Radial Ramp in Plate Separation

Pattern formation in Rayleigh-Benard convection in a large-aspect-ratio cylinder with a radial ramp in the plate separation is studied analytically and numerically by performing numerical simulations of the Boussinesq equations. A horizontal mean flow and a vertical large scale counterflow are quantified and used to understand the pattern wavenumber. Our results suggest that the mean flow, generated by amplitude gradients, plays an important role in the roll compression observed as the control parameter is increased. Near threshold the mean flow has a quadrupole dependence with a single vortex in each quadrant while away from threshold the mean flow exhibits an octupole dependence with a counter-rotating pair of vortices in each quadrant. This is confirmed analytically using the amplitude equation and Cross-Newell mean flow equation. By performing numerical experiments the large scale counterflow is also found to aid in the roll compression away from threshold but to a much lesser degree. Our results yield an understanding of the pattern wavenumbers observed in experiment away from threshold and suggest that near threshold the mean flow and large scale counterflow are not responsible for the observed shift to smaller than critical wavenumbers.Comment: 10 pages, 13 figure