Convective instability and transient growth in flow over a backward-facing step

Transient energy growths of two- and three-dimensional optimal linear perturbations to two-dimensional flow in a rectangular backward-facing-step geometry with expansion ratio two are presented. Reynolds numbers based on the step height and peak inflow speed are considered in the range 0–500, which is below the value for the onset of three-dimensional asymptotic instability. As is well known, the flow has a strong local convective instability, and the maximum linear transient energy growth values computed here are of order 80×103 at Re = 500. The critical Reynolds number below which there is no growth over any time interval is determined to be Re = 57.7 in the two-dimensional case. The centroidal location of the energy distribution for maximum transient growth is typically downstream of all the stagnation/reattachment points of the steady base flow. Sub-optimal transient modes are also computed and discussed. A direct study of weakly nonlinear effects demonstrates that nonlinearity is stablizing at Re = 500. The optimal three-dimensional disturbances have spanwise wavelength of order ten step heights. Though they have slightly larger growths than two-dimensional cases, they are broadly similar in character. When the inflow of the full nonlinear system is perturbed with white noise, narrowband random velocity perturbations are observed in the downstream channel at locations corresponding to maximum linear transient growth. The centre frequency of this response matches that computed from the streamwise wavelength and mean advection speed of the predicted optimal disturbance. Linkage between the response of the driven flow and the optimal disturbance is further demonstrated by a partition of response energy into velocity components

Convective instability and transient growth in steady and pulsatile stenotic flows

We show that suitable initial disturbances to steady or long-period pulsatile flows in a straight tube with an axisymmetric 75%-occlusion stenosis can produce very large transient energy growths. The global optimal disturbances to an initially axisymmetric state found by linear analyses are three-dimensional wave packets that produce localized sinuous convective instability in extended shear layers. In pulsatile flow, initial conditions that trigger the largest disturbances are either initiated at, or advect to, the separating shear layer at the stenosis in phase with peak systolic flow. Movies are available with the online version of the paper

Encapsulated formulation of the Selective Frequency Damping method

We present an alternative "encapsulated" formulation of the Selective Frequency Damping method for finding unstable equilibria of dynamical systems, which is particularly useful when analysing the stability of fluid flows. The formulation makes use of splitting methods, which means that it can be wrapped around an existing time-stepping code as a "black box". The method is first applied to a scalar problem in order to analyse its stability and highlight the roles of the control coefficient $\chi$ and the filter width $\Delta$ in the convergence (or not) towards the steady-state. Then the steady-state of the incompressible flow past a two-dimensional cylinder at $Re=100$, obtained with a code which implements the spectral/hp element method, is presented

High-order methods for the numerical solution of the BiGlobal linear stability eigenvalue problem in complex geometries.

A high-order computational tool based on spectral and spectral/hp elements (J. Fluid. Mech. 2009; to appear) discretizations is employed for the analysis of BiGlobal fluid instability problems. Unlike other implementations of this type, which use a time-stepping-based formulation (J. Comput. Phys. 1994; 110(1):82–102; J. Fluid Mech. 1996; 322:215–241), a formulation is considered here in which the discretized matrix is constructed and stored prior to applying an iterative shift-and-invert Arnoldi algorithm for the solution of the generalized eigenvalue problem. In contrast to the time-stepping-based formulations, the matrix-based approach permits searching anywhere in the eigenspace using shifting. Hybrid and fully unstructured meshes are used in conjunction with the spatial discretization. This permits analysis of flow instability on arbitrarily complex 2-D geometries, homogeneous in the third spatial direction and allows both mesh (h)-refinement as well as polynomial (p)-refinement. A series of validation cases has been defined, using well-known stability results in confined geometries. In addition new results are presented for ducts of curvilinear cross-sections with rounded corners

Spectral/hp element methods: recent developments, applications, and perspectives

The spectral/hp element method combines the geometric flexibility of the classical h-type finite element technique with the desirable numerical properties of spectral methods, employing high-degree piecewise polynomial basis functions on coarse finite element-type meshes. The spatial approximation is based upon orthogonal polynomials, such as Legendre or Chebychev polynomials, modified to accommodate C0-continuous expansions. Computationally and theoretically, by increasing the polynomial order p, high-precision solutions and fast convergence can be obtained and, in particular, under certain regularity assumptions an exponential reduction in approximation error between numerical and exact solutions can be achieved. This method has now been applied in many simulation studies of both fundamental and practical engineering flows. This paper briefly describes the formulation of the spectral/hp element method and provides an overview of its application to computational fluid dynamics. In particular, it focuses on the use the spectral/hp element method in transitional flows and ocean engineering. Finally, some of the major challenges to be overcome in order to use the spectral/hp element method in more complex science and engineering applications are discussed