Order-of-magnitude speedup for steady states and traveling waves via Stokes preconditioning in Channelflow and Openpipeflow

Steady states and traveling waves play a fundamental role in understanding hydrodynamic problems. Even when unstable, these states provide the bifurcation-theoretic explanation for the origin of the observed states. In turbulent wall-bounded shear flows, these states have been hypothesized to be saddle points organizing the trajectories within a chaotic attractor. These states must be computed with Newton's method or one of its generalizations, since time-integration cannot converge to unstable equilibria. The bottleneck is the solution of linear systems involving the Jacobian of the Navier-Stokes or Boussinesq equations. Originally such computations were carried out by constructing and directly inverting the Jacobian, but this is unfeasible for the matrices arising from three-dimensional hydrodynamic configurations in large domains. A popular method is to seek states that are invariant under numerical time integration. Surprisingly, equilibria may also be found by seeking flows that are invariant under a single very large Backwards-Euler Forwards-Euler timestep. We show that this method, called Stokes preconditioning, is 10 to 50 times faster at computing steady states in plane Couette flow and traveling waves in pipe flow. Moreover, it can be carried out using Channelflow (by Gibson) and Openpipeflow (by Willis) without any changes to these popular spectral codes. We explain the convergence rate as a function of the integration period and Reynolds number by computing the full spectra of the operators corresponding to the Jacobians of both methods.Comment: in Computational Modelling of Bifurcations and Instabilities in Fluid Dynamics, ed. Alexander Gelfgat (Springer, 2018

Effect of distributed surface roughness on the stability of Falkner-Skan-Cooke boundary layers

Direct numerical simulations (DNS) and large eddy simulations (LES) have been carried out to investigate the impact of surface roughness with various heights on the the most unstable/least stable eigenmode developing in a Falkner-Skan-Cooke boundary layer. Crossflow vortices are excited using a spanwise array of cylindrical roughness elements that are introduced by application of the virtual boundary method (VBM). The height of the roughness elements are in the range of medium to large, meaning that they are completely submerged into the boundary layer but close to the critical roughness height needed to trigger transition to turbulence (roughness Reynolds number, Rek for the investigated roughnesses are about 146 and 284). An impulse response is used to effectively trigger the eigenmodes of the flow, and the most unstable/least stable mode is obtained by using the outflow solution that has been filtered through the fringe region as input perturbation, and iterating the response in time. The method of forcing crossflow vortices using VBM proves to work very well with the present pseudo-spectral code and the results are promising. It turns out that when the height of the roughness is close to the critical height for transition, only a slight increase in roughness height will change the stability characteristics drastically, exciting a high-frequency secondary instability as opposed to an otherwise excited low-frequency mode. Differences in the excited modes are traced back to the shear layers developing in the near-field region of the roughness wake. It is concluded that a critical level of shear stress must be surpassed as the roughness size is increased and that this is responsible for triggering the secondary instability.Validerat; 20121120 (global_studentproject_submitter

Global stability analysis of three-dimensional boundary layer flows

This thesis considers the stability and transition of incompressible boundary layers. In particular, the Falkner–Skan–Cooke boundary layer subject to a cylindrical surface roughness, and the Blasius boundary layer with applied localized suction are investigated. These flows are of great importance within the aviation industry, feature complex transition scenarios, and are strongly three-dimensional in nature. Consequently, no assumptions regarding homogeneity in any of the spatial directions are possible, and the stability of the flow is governed by an extensive three-dimensional eigenvalue problem. The stability of these flows is addressed by high-order direct numerical simulations using the spectral element method, in combination with a Krylov subspace projection method. Such techniques target the long-term behavior of the flow and can provide lower limits beyond which transition is unavoidable. The origin of the instabilities, as well as the mechanisms leading to transition in the aforementioned cases are studied and the findings are reported. Additionally, a novel method for computing the optimal forcing of a dynamical system is developed. This type of analysis provides valuable information about the frequencies and structures that cause the largest energy amplification in the system. The method is based on the inverse power method, and is discussed in the context of the one-dimensional Ginzburg–Landau equation and a two-dimensional flow case governed by the Navier–Stokes equations.QC 20151015</p

Studies on instability and optimal forcing of incompressible flows

This thesis considers the hydrodynamic instability and optimal forcing of a number of incompressible flow cases. In the first part, the instabilities of three problems that are of great interest in energy and aerospace applications are studied, namely a Blasius boundary layer subject to localized wall-suction, a Falkner–Skan–Cooke boundary layer with a localized surface roughness, and a pair of helical vortices. The two boundary layer flows are studied through spectral element simulations and eigenvalue computations, which enable their long-term behavior as well as the mechanisms causing transition to be determined. The emergence of transition in these cases is found to originate from a linear flow instability, but whereas the onset of this instability in the Blasius flow can be associated with a localized region in the vicinity of the suction orifice, the instability in the Falkner–Skan–Cooke flow involves the entire flow field. Due to this difference, the results of the eigenvalue analysis in the former case are found to be robust with respect to numerical parameters and domain size, whereas the results in the latter case exhibit an extreme sensitivity that prevents domain independent critical parameters from being determined. The instability of the two helices is primarily addressed through experiments and analytic theory. It is shown that the well known pairing instability of neighboring vortex filaments is responsible for transition, and careful measurements enable growth rates of the instabilities to be obtained that are in close agreement with theoretical predictions. Using the experimental baseflow data, a successful attempt is subsequently also made to reproduce this experiment numerically. In the second part of the thesis, a novel method for computing the optimal forcing of a dynamical system is developed. The method is based on an application of the inverse power method preconditioned by the Laplace preconditioner to the direct and adjoint resolvent operators. The method is analyzed for the Ginzburg–Landau equation and afterwards the Navier–Stokes equations, where it is implemented in the spectral element method and validated on the two-dimensional lid-driven cavity flow and the flow around a cylinder.QC 20171124</p

Long-wave instabilities of two interlaced helical vortices

International audienc

Long-wave instabilities of two interlaced helical vortices

International audienc

Local and global pairing instabilities of two interlaced helical vortices

International audienceWe 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 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, global pairing results in a leapfrogging phenomenon, which temporarily restores the initial double-helix geometry, in agreement with recent observations from numerical simulations

Long-wave instabilities of two interlaced helical vortices

International audienc

Stability and sensitivity of a crossflow-dominated Falkner–Skan–Cooke boundary layer with discrete surface roughness

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. QC 20161125</p

Long-wave instabilities in the wake of a two-bladed rotor

International audienc