Instabilities in Fluids - Semester 2, 2018, ANU
This repository contains notes and material for the class "Instabilities in Fluids" (EMSC3050/EMSC4050/EMSC8014) taught in 2018 at the Australian National University.
The main subject of the class is hydrodynamic instabilities. It is aimed to provide with the students both analytical and numerical skills to attack problems in hydrodynamic stability and also understand the underlying physical mechanisms responsible for those instabilities.
The first 2-3 lectures consist a gentle stroll through the Kelvin-Helmholtz instability, a manifestation of which is seen in the clouds over Mount Duval, NSW below:
Then each student will undertake a 2-week project with a theme from a list (to be finalized). The students will present their findings in the rest of the class. Throughout the course each student will do three projects.
Some project examples
A few example of proposed projects include:
Stability analysis of classical textbook hydrodynamic instabilities with basic states that consist of a flow u(x, y, z) = (U(z), 0, 0) and density ρ(x, y, z) = ρ(z), like:
[Courtesy of the Institute of Theoretical Physics, WWU Münster]
[Courtesy of NOAA / GOES / NASA]
- Intro to Stability Theory
- Hydrodynamic Stability: Parrallel stratified shear flows
- Notes on exponential solution in linear systems
- Notes on material lines and vorticity lines in an inviscid fluid
- Kelvin-Helmholtz instability: analytical study
- Kelvin-Helmholtz instability: numerical study
- Intro to Generalized Stability Theory
- The Eady model of baroclinic instability: modal and non-modal energy growth
In addition, the
lectures folder also contains a reactive notebook named
nonnormalgrowth101.jl that was created with Pluto.jl and was used during an introductory discussion regarding non-normality at the Climate and Fluid Physics Group at ANU (August 2020).
Below is a list for the first projects (final)
|1||Study of stratified Kelvin-Helmholtz profile: Analytical, numerical. Mechanistic interpretation through edge-wave interactions. (Nick - project report)||Drazin & Reid (2004) ch. 1, pp. 14-22; Vallis (2017), ch. 9.1-9.2, pp. 335-345 (or in the 2006 edition, ch. 6.1-6.2)|
|2||Study of Rayleigh's profile: Analytical, numerical. Mechanistic interpretation through edge-wave interactions. (Ruth - project report)||Drazin & Reid, pp. 144-147; Vallis (2017), ch. 9.1-9.2, pp. 335-345 (or in the 2006 edition, ch. 6.1-6.2)|
|3||Study of Taylor's profile: Analytical, numerical. Mechanistic interpretation through edge-wave interactions. (Fergus - project report)||Drazin & Reid, pp. 144-147; Vallis (2017), ch. 9.1-9.2, pp. 335-345 (or in the 2006 edition, ch. 6.1-6.2)|
|4||Study of Holmboe's profile: Analytical, numerical. Mechanistic interpretation through edge-wave interactions.|
|5||General instability criteria for parallel shear flows (stratified and unstratified): Squire's theorem. The Rayleigh and Fjortoft's criteria for unstratified instabilities and the Miles & Howard criterion for stratified. (Kieran - project report)||Drazin & Reid (2004) ch. 4, pp. 124-143; Charru (2011), ch. 4, pp. 104-116; Howard, JFM (1961)|
|Geophysical fluid dynamics: baroclinic instability|
|6||Eady's model of baroclinic instability (Xihan - project report)||Vallis 2017, ch. 9, pp. 349-369 (or in the 2006 edition, ch. 6.4-6.7)|
|7||Phillips' model of baroclinic instability (Shangyu - project report)||Vallis 2017, ch. 9, pp. 349-369 (or in the 2006 edition, ch. 6.4-6.7)|
|8||Rayleigh–Bernard convection||Drazin & Reid book, Gayen lecture notes, Doering WHOI lectures 2008, Rayleigh 1916 paper ()?|
Projects 1-4 are closely related to the unstratified Kelvin-Helmholtz instability we discussed in class. However, I'd like you to further understand the instabilities that rise through interaction of edge-waves (i.e., the waves that "live" on the discontinuities of the velocity shear profiles or the discontinuity of the density profiles). Project 4 might be slightly harder than Projects 1-3 in terms of algebra but not in terms of concepts.
Project 5 is for the more mathematically inclined. It involves more pen-and-paper work rather than numerics. You are asked to study some theorems that help us decide a priori whether a general velocity and density basic profile is amenable to modal instability. I will suggest some numerical investigation of some basic states to demonstrate the validity of the theorems.
Projects 6-7 are a first step in studying cyclogenesis in the mid-latidute atmospheres. These projects require to work with the, so-called, quasigeostrophic equations on a beta-plane. Beta-plane is a simplified approximation of a flow on a rotating sphere; it keeps the geometry simple (cartesian) but incorporates the basic physics due to rotation. I can help you get a hold of the quasigeostrophic equation if you haven't been familiarized with them before.
Project 8 involves the classical Rayleigh-Benard convection instability (2D) for two horizontal plates (hot plate below; cold plate on top).
Problem set for midterm break
Next follows a (non-finalized) list for projects to come:
|8||Study of the Couette profile: The illusive theorem by Squire. Transient perturbation growth of 2D perturbations and the Orr mechanism. (Shangyu - project report)||Boyd 1983, Farrell & Ioannou 1993|
|9||Transient perturbation growth in Rayleigh's profile.|
|10||The problems arising from eigenanalysis of non-normal operators & Singular Value Decomposition (Fergus - project report)|
|11||Hydrodynamic stability without eigenvalues. Spectra & pseudospectra and subcritical transition to turbulence. (Kieran - project report)||Trefethen & Embree 2005, Trefethen et al. 1993|
|12||Transient perturbation growth and optimal excitation of Eady's model. (Nick: transient growth of 3D perturbations - project report; Xihan: Contribution of the Continuous versus Discrete spectrum on transient non-normal growth - project report; Ruth: Modal and non-normal growth in the Eady problem with a topographic slope - project report)||lecture notes|
|13||From Rayleigh–Bernard to chaos: Lorenz’s butterfly effect||Lorenz 1963|
Project 8 deals with the Couette profile. As we discussed briefly in class the Couette profile is modally stable but still it does show finite-time perturbation growth. This project will involves first the discussion of Squire's theorem (see Lecture 3). Then finding the analytic solutions for the initial value problem for 2D and 3D perturbations of the Couette profile. Lastly, some discussion of the relations of the energy evolution of 2D and 3D perturbations and the connection with Squire's theorem.
Projects 9 studies the transient energy growth and optimal excitation (adjoints versus eigenmodes) of the classic piecewise profiles in projects 2. Similarly, projects 12 and 13 for the baroclinic instability problems studies in projects 6 and 7.
Projects 10 presents the problems arrising from eigenanalysis of non-normal operators (i.e. when the eigenvectors are not orthogonal to each other) and shows how singular value decomposition provides for a robust, orthonormal basis related with linear operators, which also has physical meaning..
Project 13 studies the classic paper by Lorenz in which he discovered chaos in a low-order truncated system of Rayleigh-Benard convection.