Flow instability is formally a linear concept, applicable only for infinitesimal perturbations to a steady or periodic solution to the governing equations. Nevertheless, it can be applied to turbulent flows with surprising success, yielding useful information about how to control such flows.

There is a standard procedure for performing a stability analysis on a flow. (i) Choose the governing equations and boundary conditions (e.g. the Navier--Stokes equations for the flow around a cylinder in a channel with no slip boundaries). (ii) Find a steady solution to the governing equations (e.g. by iterating with a Newton solver). (iii) Linearize the governing equations around this steady solution. This is achieved by perturbing each variable (e.g. U = U̅ + εu and P = P̅ + εp) and grouping terms by powers of ε. The terms multiplied by ε^0 are already satisfied by the steady solution. The terms multiplied by ε^2 (and higher order in ε) are very small. Therefore the linearized governing equations are the terms multiplied by ε. (iv) For a modal (eigenvalue) analysis, perform a Laplace/Fourier decomposition in time and write the linearized equations as an eigenvalue problem Aq = sBq for the Laplace/Fourier variable s. (v) Usually, the problem is discretized and solved as a generalized matrix eigenvalue problem. This is described in the following tutorial:

Modal Stability Theory
M. P. Juniper, A. Hanifi, and V. Theofilis
Applied Mechanics Review 66 024804 (2014) doi:10.1115/1.4026604
doi: https://doi.org/10.1115/1.4026604
Section 2.8 Plug flow Matlab tutorial
Section 2.9 to 2.10 Planar Poiseuille Flow Matlab tutorial
Section 4.6 Parabolized Stability Equation Matlab tutorial
Section 5.3.1 2D Helmholtz Equation Matlab tutorial
Section 5.3.2 2D Eigenvalue problem Matlab tutorial
Section 6.1 Linearized Navier-Stokes Equation Matlab tutorial

This article contains a review of modal stability theory. It covers local stability analysis of parallel flows including temporal stability, spatial stability, phase velocity, group velocity, spatio-temporal stability, the linearized Navier-Stokes equations, the Orr-Sommerfeld equation, the Rayleigh equation, the Briggs-Bers criterion, Poiseuille flow, free shear flows, and secondary modal instability. It also covers the Parabolized Stability Equation (PSE), temporal and spatial biglobal theory, 2D eigenvalue problems, 3D eigenvalue problems, spectral collocation methods, and other numerical solution methods. Computer codes are provided for tutorials described in the article. These tutorials cover the main topics of the article and can be adapted to form the basis of research codes.

As the tutorial shows, the stability analysis can be local or global. The local analysis is cheap and approximate, while the global analysis is expensive and exact. These analyses give the direct modes. For the global analysis, the adjoint modes can be obtained with a further expensive calculation. For the local analysis, however, the adjoint modes can be obtained at little extra cost. This is explained in the following paper, which (in the final paragraph of section 3) also explains how the concept of the wavemaker in the local analysis is related to that in the global analysis:

The structural sensitivity of open shear flows calculated with a local stability analysis
M. P. Juniper, B. Pier
European Journal of Mechanics B 49 426--437 (2014) doi:10.1016/j.euromechflu.2014.05.011
Open Access
doi: https://doi.org/10.1016/j.euromechflu.2014.05.011

The structural sensitivity shows where an instability of a fluid flow is most sensitive to changes in internal feedback mechanisms. It is formed from the overlap of the flow's direct and adjoint global modes. These global modes are usually calculated with 2D or 3D global stability analyses, which can be very computationally expensive. For weakly non-parallel flows the direct global mode can also be calculated with a local stability analysis, which is orders of magnitude cheaper. In this theoretical paper we show that, if the direct global mode has been calculated with a local analysis, then the adjoint global mode follows at little extra cost. We also show that the maximum of the structural sensitivity is the location at which the local k+ and k- branches have the same imaginary value. Finally, we use the local analysis to derive the structural sensitivity of two flows: a confined co-flow wake at Re = 400, for which it works very well, and the flow behind a cylinder at Re = 50, for which it works reasonably well. As expected, we find that the local analysis becomes less accurate when the flow becomes less parallel.

It is clear that linear stability analysis is valid when the steady solution is at the threshold of instability, for example at Reynolds number = 45 for the flow around a cylinder. It predicts the frequency and wavelength of these oscillations remarkably accurately.

As the Reynolds number increases, however, this periodic solution becomes unstable to a span-wise motion, which is the start of a cascade of instabilities that eventually lead to fully turbulent flow as the Reynolds number increases further.

Linear stability analysis may seem hopeless at higher Reynolds numbers. Experiments show, however, that the low dimensional behaviour returns in a fully turbulent flow:

Global instability in the flow behind a cylinder
This movie shows a the Von-Karman vortex shedding behind a cylinder.
Credit: Physics Graphics
Vortex shedding around Heard Island (bottom left of image) in the Southern Indian Ocean
Although the flow is strongly turbulent, the low-dimensional dynamics of vortex shedding is easy to see in the flow behind Heard Island in the Southern Indian Ocean. These oscillations occur at a Strouhal number around 0.2. This frequency and the mode shape can, surprisingly, be predicted from a linear stability analysis around the turbulent mean flow.
Credit: from Wikepedia

It turns out that a linear stability analysis performed about a tubulent mean flow predicts the frequency and wavelength of oscillations remarkably accurately. This is curious and, to my mind, is one of the most important questions in flow instability at the moment. My first instinct was to reject the concept. My second was to adopt the strategy I take when watching a film involving time-travel: not to think about it too much. More recently, researchers have made significant progress in understanding why this should be the case. Section 2 of the paper below contains our current thoughts and key references on this subject:

Coherent structures in a swirl injector at Re = 4800 by nonlinear simulations and linear global modes
O. Tammisola, Juniper, M. P.
Journal of Fluid Mechanics 792 620--657 (2016) doi:10.1017/jfm.2016.86
Open Access
doi: https://doi.org/10.1017/jfm.2016.86

The large-scale coherent motions in a realistic swirl fuel injector geometry are analysed by direct numerical simulations (DNS), proper orthogonal decomposition (POD), and linear global modes. The aim is to identify the origin of instability in this turbulent flow in a complex internal geometry.

The flow field in the nonlinear simulation is highly turbulent, but with a distinguishable coherent structure: the precessing vortex core (a spiraling mode). The most energetic POD mode pair is identified as the precessing vortex core. By analysing the FFT of the time coefficients of the POD modes, we conclude that the first four POD modes contain the coherent fluctuations. The remaining POD modes (incoherent fluctuations) are used to form a turbulent viscosity field, using the Newtonian eddy model.

The turbulence sets in from convective shear layer instabilities even before the nonlinear flow reaches the other end of the domain, indicating that equilibrium solutions of the Navier?Stokes are never observed. Linear global modes are computed around the mean flow from DNS, applying the turbulent viscosity extracted from POD modes. A slightly stable discrete m = 1 eigenmode is found, well separated from the continuous spectrum, in very good agreement with the POD mode shape and frequency. The structural sensitivity of the precessing vortex core is located upstream of the central recirculation zone, identifying it as a spiral vortex breakdown instability in the nozzle. Furthermore, the structural sensitivity indicates that the dominant instability mechanism is the Kelvin-Helmholtz instability at the inflection point forming near vortex breakdown. Adjoint modes are strong in the shear layer along the whole extent of the nozzle, showing that the optimal initial condition for the global mode is localized in the shear layer.

We analyse the qualitative influence of turbulent dissipation in the stability problem (eddy viscosity) on the eigenmodes by comparing them to eigenmodes computed without eddy viscosity. The results show that the eddy viscosity improves the complex frequency and shape of global modes around the fuel injector mean flow, while a qualitative wavemaker position can be obtained with or without turbulent dissipation, in agreement with previous studies.

This study shows how sensitivity analysis can identify which parts of the flow in a complex geometry need to be altered in order to change its hydrodynamic stability characteristics.

Many applications are far from the threshold of instability, but in a regime that nonetheless appears to be amenable to linear stability analysis around the mean flow. My first published paper on this subject was the application of local linear stability analysis to the flow in gas turbine fuel injectors:

Absolute and convective instability in gas turbine fuel injectors
M. P. Juniper
ASME Turbo Expo, Copenhagen, Denmark, 11-15 June 2012, GT2012-68253, (2012)

Hydrodynamic instabilities in gas turbine fuel injectors help to mix the fuel and air but can sometimes lock into acoustic oscillations and contribute to thermoacoustic instability. This paper describes a linear stability analysis that predicts the frequencies and strengths of hydrodynamic instabilities and identifies the regions of the flow that cause them. It distinguishes between convective instabilities, which grow in time but are convected away by the flow, and absolute instabilities, which grow in time without being convected away. Convectively unstable flows amplify external perturbations, while absolutely unstable flows also oscillate at intrinsic frequencies. As an input, this analysis requires velocity and density fields, either from a steady but unstable solution to the Navier--Stokes equations, or from time-averaged numerical simulations. In the former case, the analysis is a predictive tool. In the latter case, it is a diagnostic tool. This technique is applied to three flows: a swirling wake at Re = 400, a single stream swirling fuel injector at Re ~ 10^6, and a lean premixed gas turbine injector with five swirling streams at Re ~ 10^6.

Its application to the swirling wake demonstrates that this technique can correctly predict the frequency, growth rate and dominant wavemaker region of the flow. It also shows that the zone of absolute instability found from the spatio-temporal analysis is a good approximation to the wavemaker region, which is found by overlapping the direct and adjoint global modes. This approximation is used in the other two flows because it is difficult to calculate their adjoint global modes.

Its application to the single stream fuel injector demonstrates that it can identify the regions of the flow that are responsible for generating the hydrodynamic oscillations seen in LES and experimental data. The frequencies predicted by this technique are within a few percent of the measured frequencies. The technique also explains why these oscillations become weaker when a central jet is injected along the centreline. This is because the absolutely unstable region that causes the oscillations becomes convectively unstable.

Its application to the lean premixed gas turbine injector reveals that several regions of the flow are hydrodynamically unstable, each with a different frequency and a different strength. For example, it reveals that the central region of confined swirling flow is strongly absolutely unstable and sets up a precessing vortex core, which is likely to aid mixing throughout the injector. It also reveals that the region between the second and third streams is slightly absolutely unstable at a frequency that is likely to coincide with acoustic modes within the combustion chamber. This technique, coupled with knowledge of the acoustic modes in a combustion chamber, is likely to be a useful design tool for the passive control of mixing and combustion instability.

This summarized work that had been done in 2008 and that won The Engineer magazine 2009 award for Environmental Technology. In particular, this revealed the region of the gas turbine fuel injector that causes the vortex core to precess at a well-defined frequency. This is important for mixing of the fuel and air and may also be important for lock-in with thermoacoustic oscillations.

Next, we performed a more thorough study of the hydrodynamic stability of the flow from a gas turbine fuel injector. We performed Direct Numerical Simulations (DNS) and extracted the mean flow field and the turbulence statistics. We then performed a global linear stability analysis using only the mean fields and compared our results with those from POD analysis of the original DNS data. The results were almost identical:

Comparison of global modes with POD modes for the flow from a gas turbine fuel injector
We performed Direct Numerical Simulations (DNS) of the flow from a gas turbine fuel injector at Re = 4800. From these simulations we extracted the turbulent mean flow, the turbulence statistics, the local eddy viscosity, and the Proper Orthogonal Decomposition (POD) modes. We then performed a linear stability analysis on the turbulent mean flow, using the local turbulence viscosity. The most unstable global mode is shown in the top line (axial, radial, and azimuthal velocity). Its shape and frequency is very similar to the first POD mode, which is shown in the bottom row. This shows that the linear analysis, performed around the turbulent mean flow, gives useful results in this application.
Credit: Outi Tammisola
Jump to publication (will be at top of next screen)

This also enabled us to identify the wavemaker region of the flow, which sits just upstream of the injector exit and is the region that must be altered in order to alter the behaviour of this hydrodynamic mode.

We also performed a local stability analysis on experimental data from a reacting flow:

Local stability analysis and eigenvalue sensitivity of reacting bluff body wakes
B. Emerson, T. C. Lieuwen, M. P. Juniper
Journal of Fluid Mechanics 788 549-575 (2016) doi:10.1017/jfm.2015.724
Open Access
doi: https://doi.org/10.1017/jfm.2015.724

This paper presents an experimental and theoretical investigation of high Reynolds number, low density reacting wakes near a hydrodynamic Hopf bifurcation. This configuration is applicable to the wake flows that are commonly used to stabilize flames in high velocity flows. First, an experimental study is conducted to measure the limit cycle oscillation of this reacting bluff body wake. The experiment is repeated while independently varying the bluff body lip velocity and the density ratio across the flame. In all cases, the wake exhibits a sinuous oscillation. Linear stability analysis is performed on the measured time-averaged velocity and density fields. In the first stage of this analysis, a local, spatio-temporal stability analysis is performed on the measured, time averaged velocity and density fields. The stability analysis results are compared to the experimental measurement, and demonstrate that the local stability analysis correctly captures the influence of the lip velocity and density ratio parameters on the sinuous mode. In the second stage of this analysis, the linear direct and adjoint global modes are estimated by combining the local results. The sensitivity of the eigenvalue to changes in intrinsic feedback mechanisms is found by combining the direct and adjoint global modes. This is referred to as the eigenvalue sensitivity throughout the paper for reasons of brevity. The predicted global mode frequency is consistently within 10 % of the measured value, and the linear global mode shape closely resembles the measured nonlinear oscillations. The adjoint global mode reveals that the oscillation is strongly sensitive to open loop forcing in the shear layers. The eigenvalue sensitivity identifies a wavemaker in the recirculation zone of the wake. A parametric study shows that these regions change little when the density ratio and lip velocity change. In the third stage of the analysis, the stability analysis is repeated for the varicose hydrodynamic mode. Although not physically observed in this unforced flow, the varicose mode can lock into longitudinal acoustic waves and cause thermoacoustic oscillations to occur. This paper shows that the local stability analysis successfully predicts the global hydrodynamic stability characteristics of this flow and shows that experimental data can be post-processed with this method in order to identify the wavemaker regions and the regions that are most sensitive to external forcing, for example from acoustic waves.

We extracted the direct and adjoint global modes from a local analysis, using the technique in (Juniper and Pier, Eur. J. Mech. B 49, 426--437) and found that the direct modes were very similar to the oscillations observed experimentally. With the adjoint modes, we could also identify the wavemaker regions, which show how to control these flows.

We have since applied local and global stability analysis to cyclone injectors, where the linear analysis around the mean flow gives useful information about the source of noise:

Local linear stability analysis of cyclone separators
T. Grimble, A. Agarwal, M. P. Juniper
Journal of Fluid Mechanics 816 507--538 (2017) doi:10.1017/jfm.2017.89
Open Access
doi: https://doi.org/10.1017/jfm.2017.89

Local linear stability analysis is applied to the flow inside a cyclone separator to investigate the unsteady precession of the vortex core. The results of the stability analysis are compared with experimental measurements of the vortex oscillations using high speed photography with particle seeding, and hot wire anemometry. The experiments reveal distinct spatial variation in the oscillation behaviour within the cyclones. The unsteady motion is focused at each end of the device, at both the narrow cone tip and just below the exhaust duct at the top of the cone, which is known as a vortex finder. The local stability analysis shows that an absolute instability is present throughout the flow for some non-zero azimuthal wavenumbers. The unsteady flow is observed to be driven by coupling between the shear layer and inertial waves confined within the vortex core. Comparing the stability analysis with experiments shows the same frequency and mode shape behaviour and suggests that the local analysis accurately predicts the unstable modes of the system. The precessing vortex core is responsible for a narrow-band acoustic noise. Comparisons are also drawn with acoustic measurements made on cyclones in which the system is defined by key non-dimensional parameters, such as the swirl number and outlet diameter ratio. The results in this study demonstrate the applicability of local stability analysis to a complex swirling system and yield credible details about the underlying mechanisms of the unstable flow inside the cyclone.

This work is on-going and its publication is usually delayed for commercial reasons.

The software suite developed by the group for local stability analysis is called Instaflow. Broad details can be found in two theses by the group:

Hydrodynamic instability of confined jets and wakes and implications for gas turbine fuel injectors
S. J. Rees
University of Cambridge, (2009), examined by B. Pier and N. Peake

This dissertation investigates the stability of injector flows. This is carried out both theoretically and numerically.

In injector flows three main features are identified which affect the stability of the flow. These are: shear, geometry and density and are given in the relative order of im- portance for the consideration of this dissertation.

Shear is the primary instability mechanism within an injector flow. In order to capture this physical mechanism the simplest flow with shear is considered: the inviscid single vortex sheet. This is unstable due to the Kelvin--Helmholtz instability and forms the building block with which to construct various models of injector flows. Variants of this construct include the inclusion of surface tension at the interface and a finite thickness shear layer. Injector flows are most simply modelled by considering two shear layers interacting. Depending upon the relative velocity of the different streams the flow can describe a jet or a wake.

The second feature, geometry, is introduced into the model by placing confining walls either side of the two shear layers. It is shown that the configuration of these confining walls has a profound effect on the instability of the flow and can in some case make the flow much more unstable. Further realism is added by introducing curvature by considering a round geometry. Many of the results in the planar case are carried over into the round case.

The third feature, density, is explored briefly in this dissertation and is found to also have a profound effect on the stability. In particular low density jets and high density wake configurations are found to be strongly unstable. Density does not receive nearly as much attention as does shear and geometry since in practical terms it is largely fixed with little scope for wide-scale variation. The other two parameters by comparison can be chosen over a wide range of values in a practical setup.

Even these simple models are still capable of producing very complex stability characteristics. These models, however, represent the limit of the theoretical studies. In order to progress any further and add more realism to the model, either in the form of viscosity or smooth velocity profiles it was necessary to adopt a numerical approach. This has led to the develop of FLOWTOOL, a piece of software capable of calculating a spatio-temporal analysis of a given velocity profile and determining the local stability properties. The code is successfully demonstrated on a real injector flow. Excellent agreement is found between the predicted frequencies and those obtained from global methods, namely a Large Eddy Simulation.

Linear amplification analysis for extraction of coherent structures in wall-bounded turbulent flows
V. Gupta
University of Cambridge, (2014), examined by A. Sharma and R. S. Cant

Coherent structures in turbulent flows provide a means of understanding turbulence in terms of large organised motions. Understanding the mechanism of formation of coherent structures can be helpful in suppressing or enhancing the turbulence in a flow by means of active or passive control devices. Knowledge of the Reynolds number scaling of the size and energy content of coherent structures can extend the knowledge to high Reynolds number flows, which are out of reach of the present computational and experimental facilities.

In this thesis, linear amplification and eigenvalue stability analyses are performed by linearising the Navier?-Stokes or Reynolds-averaged Navier-?Stokes (RANS) equations over the mean flow profiles in several wall-bounded turbulent shear flows. It is investigated whether the linear optimal modes or the leading eigenmodes approximate the coherent structures in fully nonlinear turbulent flows. This is done by comparing various kinematic properties of the optimal modes, such as the shape and energy spectra, with those of the observed coherent structures in turbulent channel and pipe flows in the first half of the thesis. The use of the linearised Navier-?Stokes equations in the regions of high mean shear in the flows is justified based on rapid distortion theory. In the linearised RANS equations-based analysis, turbulence models are used to account for the effect of wave-induced perturbations in the Reynolds stress on the behaviour of small external wave motions. The turbulence models used in this thesis are the eddy viscosity model (EVM) and the explicit algebraic Reynolds stress model (EARSM). The focus of this thesis is to investigate whether this effect of wave-induced perturbations in the Reynolds stress needs to be included in stability analysis of wall-bounded turbulent flows.

The linear amplification analysis based on the Navier?Stokes equations finds three main types of structures in turbulent channel flows. The first type are the small streamwise wavelength (lambda x + = 200 ? 800) structures, which are found to scale in inner units and have preferred spanwise wavelength equal to around one hundred wall-units. These properties match well with those of observed near-wall structures. The second type are the intermediate streamwise wavelengths (from lambda x + > 800 to lambda x < 3) structures which correspond to hairpin vortical and large-scale streaky like structures. The peak in energy amplification in this wavelength range found from the analysis matches well with that from DNS. Various kinematic properties, such as the inclination angle of streaks with the wall, also match with those of large-scale-motions (LSMs) observed in experiments. The third type are the large streamwise wavelength (lambda x >= 6) structures. The preferred spanwise wavelength of these structures (lambda z peak ~ 2), their scaling in outer units, and the fact that they extend to the wall match with the observed features of very-large-scale-motions (VLSMs). All these results show that the most optimal modes obtained from the linearised Navier?Stokes equations, without any turbulence model or eddy viscosity, share many important features with those of observed coherent structures in turbulent channel flows.

In comparison, the results from the EVM- and EARSM-based linear amplification analyses find only two types of coherent structures. One type are of the small wavelengths, which correspond to the near-wall structures, and the other type are of the large wavelengths, which correspond to the VLSMs. These analyses, however, find minima in energy spectra in the intermediate wavelength region, where DNS and the Navier?Stokes equations-based analysis find maxima in energy spectra.

In axially rotating turbulent pipe flows, it is found from the linearised Navier?-Stokes equations-based analysis that rotation causes the widening of streaks and prevents the formation of quasi-streamwise vortices. These results match well with observations from DNS, which further shows the usefulness of the linearised Navier?Stokes equations.

In the second part of the thesis, stability analyses based on the linearised Navier-?Stokes and RANS equations are applied in more complex flows. Based on the results from the stability analyses for flows in gas-turbine systems, it is found that for such flows the inclusion of turbulence models in stability analysis has no significant qualitative effect on the results. This is because these instabilities are driven by regions of high mean shear for which analysis based on the linearised Navier?-Stokes equations is sufficient. It is also found from stability analysis that an expansion at the nozzle exit and swirl in the flow are destabilising, and therefore increase hydrodynamic instability.

Based on the preliminary comparisons of stability results and observations from DNS in Taylor-Couette flows, it is again concluded that the linearised Navier?-Stokes equations-based analysis is better at capturing intermittent coherent structures as compared to the linearised RANS equations-based analysis.

It is concluded in this thesis that the linearised Navier?-Stokes equations-based analysis, which does not require any turbulence model, can be used to find information about coherent structures in high mean shear flows, such as the flows in gas-turbine fuel injectors or wall-bounded turbulent flows.

Finally, we were able to use linear stability analysis to show why the strakes around a chimney stabilize the vortex shedding mode.

Strakes around thin-walled chimneys
Strakes are often placed around thin-walled chimneys in order to disrupt vortex shedding. Experiments show that, for best results, the strakes must wrap round the chimney roughly once every 6 chimney diameters. This stabilization can be explained as a second-order effect in a linear stability analysis.
Credit: Matthew Juniper
Jump to publication (will be at top of next screen)
Second-order perturbation of global modes and implications for spanwise wavy actuation
O. Tammisola, F. Giannetti, V. Citro and M. P Juniper
Journal of Fluid Mechanics 755 314--335 (2014) doi:10.1017/jfm.2014.415
Open Access
doi: https://doi.org/10.1017/jfm.2014.415

Sensitivity analysis has successfully located the most efficient regions in which to apply passive control in many globally unstable flows. As is shown here and in previous studies, the standard sensitivity analysis, which is linear (first order) with respect to the actuation amplitude, predicts that steady spanwise wavy alternating actuation/modification has no effect on the stability of planar flows, because the eigenvalue change integrates to zero in the spanwise direction. In experiments, however, spanwise wavy modification has been shown to stabilize the flow behind a cylinder quite efficiently. In this paper, we generalize sensitivity analysis by examining the eigenvalue drift (including stabilization/destabilization) up to second order in the perturbation, and show how the second-order eigenvalue changes can be computed numerically by overlapping the adjoint eigenfunction with the first-order global eigenmode correction, shown here for the first time. We confirm the prediction against a direct computation, showing that the eigenvalue drift due to a spanwise wavy base flow modification is of second order. Further analysis reveals that the second-order change in the eigenvalue arises through a resonance of the original (2-D) eigenmode with other unperturbed eigenmodes that have the same spanwise wavelength as the base flow modification. The eigenvalue drift due to each mode interaction is inversely proportional to the distance between the eigenvalues of the modes (which is similar to resonance), but also depends on mutual overlap of direct and adjoint eigenfunctions (which is similar to pseudoresonance). By this argument, and by calculating the most sensitive regions identified by our analysis, we explain why an in-phase actuation/modification is better than an out-of-phase actuation for control of wake flows by spanwise wavy suction and blowing. We also explain why wavelengths several times longer than the wake thickness are more efficient than short wavelengths.

This turns out to be a second order effect. This effect is also instrumental at lowering the critical transition Reynolds number in a nearly-axisymmetric stenotic flow:

Diagram of a stenosis with slight excentricity (breaking of axisymmetry)
A stenosis is an almost axisymmstric constriction in a round pipe flow. Experiments show transition to turbulence at a rather low Reynolds number. Numerical simulations, however, show transition to turbulence at a rather high Reynolds number. Through a second order stability analysis, we were able to show that even a slight breaking of the axisymmetry of the problem leads to a dramatic reduction in the Reynolds number at which the flow becomes unstable.
Credit: John Samuelsson
Jump to publication (will be at top of next screen)
Breaking axi-symmetry in stenotic flow lowers the critical transition Reynolds number
J. Samuelson, O. Tammisola, M. P. Juniper
Physics of Fluids 27 104103 (2015) doi:10.1063/1.4934530
Open Access
doi: https://doi.org/10.1063/1.4934530

Flow through a sinuous stenosis with varying degrees of non-axisymmetric shape variations and at Reynolds number ranging from 250 to 750 is investigated using direct numerical simulation (DNS) and global linear stability analysis. At low Reynolds numbers (Re < 390), the flow is always steady and symmetric for an axisymmetric geometry. Two steady state solutions are obtained when the Reynolds number is increased: a symmetric steady state and an eccentric, non-axisymmetric steady state. Either one can be obtained in the DNS depending on the initial condition. A linear global stability analysis around the symmetric and non-axisymmetric steady state reveals that both flows are linearly stable for the same Reynolds number, showing that the first bifurcation from symmetry to antisymmetry is subcritical. When the Reynolds number is increased further, the symmetric state becomes linearly unstable to an eigenmode, which drives the flow towards the non-axisymmetric state. The symmetric state remains steady up to Re = 713, while the non-axisymmetric state displays regimes of periodic oscillations for Re >= 417 and intermittency for Re >~ 525. Further, an offset of the stenosis throat is introduced through the eccentricity parameter E. When eccentricity is increased from zero to only 0.3% of the pipe diameter, the bifurcation Reynolds number decreases by more than 50%, showing that it is highly sensitive to non-axisymmetric shape variations. Based on the resulting bifurcation map and its dependency on E, we resolve the discrepancies between previous experimental and computational studies. We also present excellent agreement between our numerical results and previous experimental results.