We investigate theoretically and experimentally the stability of two interlaced helical vortices with respect to displacement perturbations having wavelengths that are large compared to the size of the vortex cores. First, existing theoretical results are recalled and applied to the present configuration. Various modes of unstable perturbations, involving different phase relationships between the two vortices, are identified and their growth rates are calculated. They lead to a local pairing of neighbouring helix loops, or to a global pairing with one helix expanding and the other one contracting. A relation is established between this instability and the three-dimensional pairing of arrays of straight parallel vortices, and a striking quantitative agreement concerning the growth rates and frequencies is found. This shows that the local pairing of vortices is the driving mechanism behind the instability of the helix system. Second, an experimental study designed to observe these instabilities in a real flow is presented. Two helical vortices are generated by a two-bladed rotor in a water channel and characterised through dye visualisations and particle image velocimetry measurements. Unstable displacement modes are triggered individually, either by varying the rotation frequency of the rotor, or by imposing a small rotor eccentricity. The observed unstable mode structure, and the corresponding growth rates obtained from advanced processing of visualisation sequences, are in good agreement with theoretical predictions. The nonlinear late stages of the instability are also documented experimentally. Whereas local pairing leads to strong deformations and subsequent breakup of the vortices, global pairing results in a leapfrogging phenomenon, which temporarily restores the initial double-helix geometry, in agreement with recent observations from numerical simulations.

The need to numerically represent a free vortex system arises frequently in fundamental and applied research. Many possible techniques for realizing this vortex system exist but most tend to prioritize accuracy either inside or outside of the vortex core, which therefore makes them unsuitable for a stability analysis considering the entire flow field. In this article, a simple method is presented that is shown to yield an accurate representation of the flow inside and outside of the vortex core. The method is readily implemented in any incompressible Navier–Stokes solver using primitive variables and Cartesian coordinates. It can potentially be used to model a wide range of vortices but is here applied to the case of two helices, which is of renewed interest due to its relevance for wind turbines and helicopters. Three-dimensional stability analysis is performed in both a rotating and a translating frame of reference, which yield eigenvalue spectra that feature both mutual inductance and elliptic instabilities. Comparison of these spectra with available theoretical predictions is used to validate the proposed baseflow model, and new insights into the elliptic instability of curved Batchelor vortices are presented. Furthermore, it is shown that the instabilities in the rotating and the translating reference frames have the same structure and growth rate, but different frequency. A relation between these frequencies is provided.

In this study the origins of premature transition due to oversuction in boundary layers are studied. An infinite row of circular suction pipes that are mounted at right angles to a flat plate subject to a Blasius boundary layer is considered. The interaction between the flow originating from neighbouring holes is weak and for the parameters investigated, the pipe is always found to be unsteady regardless of the state of the flow in the boundary layer. A stability analysis reveals that the appearance of boundary layer transition can be associated with a linear instability in the form of two unstable eigenmodes inside the pipe that have weak tails, which extend into the boundary layer. Through an energy budget and a structural sensitivity analysis, the origin of this flow instability is traced to the structures developing inside the pipe near the pipe junction. Although the amplitudes of the modes in the boundary layer are orders of magnitude smaller than the corresponding amplitudes inside the pipe, a Koopman analysis of the data gathered from a nonlinear direct numerical simulation confirms that it is precisely these disturbances that are responsible for transition to turbulence in the boundary layer due to oversuction.

The need to numerically represent a free vortex system arises frequently in fundamental and applied research. Many possible techniques for realizing this vortex system exist but most tend to prioritize accuracy either inside or outside of the vortex core, which therefore makes them unsuitable to for a stability analysis considering the entire flow field. In this article, a simple method is presented that is shown to yield an accurate representation of the flow inside and outside of the vortex core. The method is readily implemented in any incompressible Navier–Stokes solver using primitive variables and Cartesian coordinates. It can potentially be used to model a wide range of vortices but is here applied to reproduce a recent experiment by Quaranta et al. (2017) considering two helices. A three-dimensional stability analysis is performed and yields an eigenvalue spectrum that features both long- and short-wave instabilities.

For problems governed by a non-normal operator, the leading eigenvalue of the operator is of limited interest and a more relevant measure of the stability is obtained by considering the harmonic forcing causing the largest system response. Various methods for determining this so-called optimal forcing exist, but they all suffer from great computational expense and are hence not practical for large-scale problems. In the present paper a new method is presented, which is applicable to problems of arbitrary size. The method does not rely on timestepping, but on the solution of linear systems, in which the inverse Laplacian acts as a preconditioner. By formulating the search for the optimal forcing as an eigenvalue problem based on the resolvent operator, repeated system solves amount to power iterations, in which the dominant eigenvalue is seen to correspond to the energy amplification in a system for a given frequency, and the eigenfunction to the corresponding forcing function. Implementation of the method requires only minor modifications of an existing timestepping code, and is applicable to any partial differential equation containing the Laplacian, such as the Navier-Stokes equations. We discuss the method, first, in the context of the linear Ginzburg-Landau equation and then, the two-dimensional lid-driven cavity flow governed by the Navier-Stokes equations. Most importantly, we demonstrate that for the lid-driven cavity, the optimal forcing can be computed using a factor of up to 500 times fewer operator evaluations than the standard method based on exponential timestepping.

We investigate theoretically and experimentally the stability of two interlaced helical vortices with respect to displacement perturbations having wavelengths that are large compared to the size of the vortex cores. First, existing theoretical results are recalled and applied to the present configuration. Various modes of unstable perturbations, involving different phase relationships between the two vortices, are identified and their growth rates are calculated. They lead to a local pairing of neighbouring helix loops, or to a uniform pairing with one helix expanding and the other one contracting. A relation is established between this instability and the three-dimensional pairing of arrays of straight parallel vortices, and a striking quantitative agreement concerning the growth rates is found, showing that the local pairing of vortices is the driving mechanism behind the instability of the helix system. Second, an experimental study designed to observe these instabilities in a real flow is presented. Two helical vortices are generated by a two-bladed rotor in a water channel and characterised through dye visualisations and PIV measurements. Unstable displacement modes are triggered individually, either by varying the rotation frequency of the rotor, or by imposing a small rotor eccentricity. The observed unstable mode structure, and the corresponding growth rates obtained from advanced processing of visualisation sequences, are in good agreement with theoretical predictions. The non-linear late stages of the instability are also documented experimentally. Whereas local pairing leads to strong deformations and subsequent break-up of the vortices, uniform pairing results in a leapfrogging phenomenon, which intermittently restores the initial double-helix geometry, in agreement with recent observations from numerical simulations.

In this work the problem of premature transition in boundary layers due to localized suction is revisited. A thorough study involving nonlinear direct numerical simulations, a three-dimensional linear stability analysis, a sensitivity study and a Koopman analysis is presented. The ensemble of these different techniques enables the origins of oversuction to be studied in great detail and provides new insight into the transition process of the flow. The configuration considered consists of an infinite row of widely separated suction pipes that are mounted to the plate at right angles. For the parameter range investigated, the flow inside the pipe is seen to bifurcate at a lower suction ratio than the boundary layer and thus act as an oscillator that forces the external flow over the plate. At low levels of suction, this forcing is not enough to cause transition in the boundary layer, but as the suction level is increased beyond criticality, modes originating from the pipe and extending into the boundary layer are seen to destabilize as well. These modes enable the perturbations forced in the pipe to also amplify in the boundary layer, which leads to a rapid breakdown to turbulence in the wake of the suction hole.

With the motivation of determining the critical roughness size, a global stability and sensitivity analysis of a three-dimensional Falkner-Skan-Cooke (FSC) boundary layer with a cylindrical surface roughness is performed. The roughness size is chosen such that breakdown to turbulence is initiated by a global version of traditional secondary instabilities of the cross-flow (CF) vortices instead of an immediate flow tripping at the roughness. The resulting global eigenvalue spectra of the systems are found to be very sensitive to numerical parameters and domain size. This sensitivity to numerical parameters is quantified using the epsilon-pseudospectrum, and the dependency on the domain is analysed through an impulse response, structural sensitivity analysis and an energy budget. It is shown that while the frequencies remain relatively unchanged, the growth rates increase with domain size, which originates from the inclusion of stronger CF vortices in the baseflow. This is reflected in a change in the rate of advective energy transport by the baseflow. It is concluded that the onset of global instability in a FSC boundary layer as the roughness height is increased does not correspond to an immediate flow tripping behind the roughness, but occurs for lower roughness heights if sufficiently long domains are considered. However, the great sensitivity results in an inability to accurately pinpoint the exact parameter values for the bifurcation, and the large spatial growth of the disturbances in the long domains eventually becomes larger than can be resolved using finite-precision arithmetic.

We present a comparison between experimental observations and theoretical predictions concerning long-wave displacement instabilities of the helical vortices in the wake of a two-bladed rotor. Experiments are performed with a small-scale rotor in a water channel, using a set-up that allows the individual triggering of various instability modes at different azimuthal wave numbers, leading to local or global pairing of successive vortex loops. The initial development of the instability and the measured growth rates are in good agreement with the predictions from linear stability theory, based on an approach where the helical vortex system is represented by filaments. At later times, local pairing develops into large-scale distortions of the vortices, whereas for global pairing the non-linear evolution returns the system almost to its initial geometry.

With the motivation of determining the critical roughness size, a global stability and sensitivity analysis of a three-dimensional Falkner–Skan–Cooke (FSC) boundary layer with a cylindrical surface roughness is performed. The roughness size is chosen such that breakdown to turbulence is initiated by a global version of traditional secondary instabilities of the crossflow (CF) vortices, instead of an immediate flow tripping at the roughness. The resulting global eigenvalue spectra of the systems are found to be very sensitive to numerical parameters and domain size. This sensitivity to numerical parameters is quantified using the "-pseudospectrum, and the dependency on the domain is analysed through an impulse response and an energy budget. It is shown that the growth rates increase with domain size, which originates from the inclusion of stronger CF vortices in the baseflow. This is reflected in a change in the rate of advective energy transport by the baseflow. It is concluded that the onset of global instability in a FSC boundary layer as the roughness height is increased does not correspond to an immediate flow tripping behind the roughness, but occurs for lower roughness heights if su ciently long domains are considered. However, the great sensitivity results in an inability to accurately pinpoint the exact parameter values for the bifurcation, and the large spatial growth of the disturbances in the long domains eventually becomes larger than what can be resolved using finite precision arithmetics.