Rayleigh–Taylor instability

Water suspended atop oil is an everyday example of Rayleigh–Taylor instability, and it may be modeled by two completely plane-parallel layers of immiscible fluid, the denser fluid on top of the less dense one and both subject to the Earth's gravity. The equilibrium here is unstable to any perturbations or disturbances of the interface: if a parcel of heavier fluid is displaced downward with an equal volume of lighter fluid displaced upwards, the potential energy of the configuration is lower than the initial state. Thus the disturbance will grow and lead to a further release of potential energy, as the denser material moves down under the (effective) gravitational field, and the less dense material is further displaced upwards. This was the set-up as studied by Lord Rayleigh. The important insight by G. I. Taylor was his realisation that this situation is equivalent to the situation when the fluids are accelerated, with the less dense fluid accelerating into the denser fluid. This occurs deep underwater on the surface of an expanding bubble and in a nuclear explosion.{{cite web | url=http://www.dtic.mil/dtic/tr/fulltext/u2/737271.pdf | archive-url=https://web.archive.org/web/20121018020911/http://www.dtic.mil/dtic/tr/fulltext/u2/737271.pdf | archive-date=October 18, 2012 | title=EVALUATION OF VARIOUS THEORETICAL MODELS FOR UNDERWATER EXPLOSION | publisher=U.S. Government | date=1971 | access-date=October 9, 2012 | author=John Pritchett | page=86}}

As the RT instability develops, the initial perturbations progress from a linear growth phase into a non-linear growth phase, eventually developing "plumes" flowing upwards (in the gravitational buoyancy sense) and "spikes" falling downwards. In the linear phase, the fluid movement can be closely approximated by linear equations, and the amplitude of perturbations is growing exponentially with time. In the non-linear phase, perturbation amplitude is too large for a linear approximation, and non-linear equations are required to describe fluid motions. In general, the density disparity between the fluids determines the structure of the subsequent non-linear RT instability flows (assuming other variables such as surface tension and viscosity are negligible here). The difference in the fluid densities divided by their sum is defined as the Atwood number, A. For A close to 0, RT instability flows take the form of symmetric "fingers" of fluid; for A close to 1, the much lighter fluid "below" the heavier fluid takes the form of larger bubble-like plumes.

This process is evident not only in many terrestrial examples, from salt domes to weather inversions, but also in astrophysics and electrohydrodynamics. For example, RT instability structure is evident in the Crab Nebula, in which the expanding pulsar wind nebula powered by the Crab pulsar is sweeping up ejected material from the supernova explosion 1000 years ago.{{Cite journal

| last = Hester | first = J. Jeff

| date = 2008

| title = The Crab Nebula: an Astrophysical Chimera

| journal = Annual Review of Astronomy and Astrophysics

| volume = 46

| pages = 127–155

| doi = 10.1146/annurev.astro.45.051806.110608

| bibcode=2008ARA&A..46..127H

}} The RT instability has also recently been discovered in the Sun's outer atmosphere, or solar corona, when a relatively dense solar prominence overlies a less dense plasma bubble.{{Cite journal

| title = Quiescent Prominence Dynamics Observed with the Hinode Solar Optical Telescope. I. Turbulent Upflow Plumes

| journal = The Astrophysical Journal

| volume = 716

| pages = 1288–1307

| doi = 10.1088/0004-637X/716/2/1288

| bibcode=2010ApJ...716.1288B

| date = 2010 | first1 = Thomas E. | display-authors = 4 | last1 = Berger | first2 = Gregory | last2 = Slater | first3 = Neal | last3 = Hurlburt | first4 = Richard | last4 = Shine | first5 = Theodore | last5 = Tarbell | first6 = Alan | last6 = Title | first7 = Bruce W. | last7 = Lites | first8 = Takenori J. | last8 = Okamoto | first9 = Kiyoshi | last9 = Ichimoto | first10 = Yukio | last10 = Katsukawa | first11 = Tetsuya | last11 = Magara | first12 = Yoshinori | last12 = Suematsu | first13 = Toshifumi | last13 = Shimizu | issue = 2

| doi-access = free }} This latter case resembles magnetically modulated RT instabilities.{{sfn|Chandrasekhar|1981|loc=Chap. X}}{{Cite journal

| last1 = Hillier | first1 = A.

| title = Numerical Simulations of the Magnetic Rayleigh–Taylor Instability in the Kippenhahn-Schlüter Prominence Model. I. Formation of Upflows

| journal = The Astrophysical Journal

| volume = 716

| issue = 2

| pages = 120–133

| doi = 10.1088/0004-637X/746/2/120

| bibcode=2012ApJ...746..120H

|

last2=Berger|

first2=Thomas|

last3=Isobe|

first3=Hiroaki|

last4=Shibata|

first4=Kazunari| year = 2012

| doi-access =free}}

{{citation|title=Single-mode instability of a ferrofluid-mercury interface under a nonuniform magnetic field|journal=Physical Review E|volume=94|year=2016|author=Singh, Chamkor|author2=Das, Arup K.|author3=Das, Prasanta K.|issue=1|page=012803|doi=10.1103/PhysRevE.94.012803|pmid=27575198|bibcode=2016PhRvE..94a2803S}}

Note that the RT instability is not to be confused with the Plateau–Rayleigh instability (also known as Rayleigh instability) of a liquid jet. This instability, sometimes called the hosepipe (or firehose) instability, occurs due to surface tension, which acts to break a cylindrical jet into a stream of droplets having the same total volume but higher surface area.

Many people have witnessed the RT instability by looking at a lava lamp, although some might claim this is more accurately described as an example of Rayleigh–Bénard convection due to the active heating of the fluid layer at the bottom of the lamp.

[[File:Evolution_of_the_Rayleigh_Taylor_Instability.xcf|thumb|left|400px|This figure represents the evolution of the Rayleigh–Taylor instability from small wavelength perturbations at the interface (a) which grow into the ubiquitous mushroom shaped spikes (fluid

structures of heavy into light fluid) and bubbles (fluid structures of

light into heavy fluid) (b) and these fluid structures interact due to

bubble merging and competition (c) eventually developing into a mixing region (d). Here ρ2 represents the heavy fluid and ρ1 represents

the light fluid. Gravity is acting downward and the system is RT

unstable.]]

The evolution of the RTI follows four main stages. In the first stage, the perturbation amplitudes are small when compared to their wavelengths, the equations of motion can be linearized, resulting in exponential instability growth. In the early portion of this stage, a sinusoidal initial perturbation retains its sinusoidal shape. However, after the end of this first stage, when non-linear effects begin to appear, one observes the beginnings of the formation of the ubiquitous mushroom-shaped spikes (fluid structures of heavy fluid growing into light fluid) and bubbles (fluid structures of light fluid growing into heavy fluid). The growth of the mushroom structures continues in the second stage and can be modeled using buoyancy drag models, resulting in a growth rate that is approximately constant in time. At this point, nonlinear terms in the equations of motion can no longer be ignored. The spikes and bubbles then begin to interact with one another in the third stage. Bubble merging takes place, where the nonlinear interaction of mode coupling acts to combine smaller spikes and bubbles to produce larger ones. Also, bubble competition takes places, where spikes and bubbles of smaller wavelength that have become saturated are enveloped by larger ones that have not yet saturated. This eventually develops into a region of turbulent mixing, which is the fourth and final stage in the evolution. It is generally assumed that the mixing region that finally develops is self-similar and turbulent, provided that the Reynolds number is sufficiently large.{{Cite journal | first1=M.S. | last1=Roberts | first2=J.W. | last2=Jacobs | title=The effects of forced small-wavelength, finite-bandwidth initial perturbations and miscibility on the turbulent Rayleigh Taylor instability | journal= Journal of Fluid Mechanics | volume=787 | date=2015 | pages=50–83 | doi=10.1017/jfm.2015.599| bibcode=2016JFM...787...50R | osti=1436483 | s2cid=126143908 }}

{{clear}}

File:Rti base.svg

The inviscid two-dimensional Rayleigh–Taylor (RT) instability provides an excellent springboard into the mathematical study of stability because of the simple nature of the base state.Drazin (2002) pp. 48–52. Consider a base state in which there is an interface, located at $z=0$ that separates fluid media with different densities, $\backslash rho\_1$ for and $\backslash rho\_2$ for $z>0$. The gravitational acceleration is described by the vector $\backslash mathbf\; g\; =\; -g\backslash ,\; \backslash mathbf\{e\}\_z$. The velocity field and pressure field in this equilibrium state, denoted with an overbar, are given by

:

where the reference location for the pressure is taken to be at $z=0$. Let this interface be slightly perturbed, so that it assumes the position $z=f(x,t)$. Correspondingly, the base state is also slightly perturbed. In the linear theory, we can write

:$\backslash mathbf\{v\}\; =\; \backslash overline\{\backslash mathbf\; v\}\; +\; \backslash hat\{\backslash mathbf\{v\}\}(z)\; e^\{ikx+\backslash sigma\; t\},\; \backslash quad\; p\; =\; \backslash overline\{p\}\; +\; \backslash hat\; p(z)\; e^\{ikx+\backslash sigma\; t\},\; \backslash quad\; f\; =\backslash hat\{f\}\; e^\{ikx+\backslash sigma\; t\}$

where $k$ is the real wavenumber in the $x$-direction and $\backslash sigma$ is the growth rate of the perturbation. Then the linear stability analysis based on the inviscid governing equations shows that

:$\backslash sigma^2\; =\; \backslash frac\{\backslash rho\_2-\backslash rho\_1\}\{\backslash rho\_2+\backslash rho\_1\}\; gk.$

Thus, if , the base state is stable and while if $\backslash rho\_2>\backslash rho\_1$, it is unstable for all wavenumbers. If the interface has a surface tension $\backslash gamma$, then the dispersion relation becomes

:$\backslash sigma^2\; =\; \backslash frac\{\backslash rho\_2-\backslash rho\_1\}\{\backslash rho\_2+\backslash rho\_1\}gk\; -\; \backslash frac\{\backslash gamma\; k^3\}\{\backslash rho\_2+\backslash rho\_1\},$

which indicates that the instability occurs only for a range of wavenumbers

:$\backslash sigma\_m^2\; =\; \backslash frac\{2\backslash gamma\}\{\backslash rho\_2-\backslash rho\_1\}\backslash left[\backslash frac\{(\backslash rho\_2-\backslash rho\_1)g\}\{3\backslash gamma\}\backslash right]^\{3/2\}.$

{{hidden begin

|toggle = left

|bodystyle = font-size: 100%

|title = Details of the linear stability analysis. A similar derivation appears in {{harvnb|Chandrasekhar|1981|at= §92, pp. 433–435}}.

}}

The perturbation introduced to the system is described by a velocity field of infinitesimally small amplitude, $(u\text{'}(x,z,t),w\text{'}(x,z,t)).\backslash ,$ Because the fluid is assumed incompressible, this velocity field has the streamfunction representation

$$\backslash textbf\{u\}\text{'}=(u\text{'}(x,z,t),w\text{'}(x,z,t))=(\backslash psi\_z,-\backslash psi\_x),\backslash ,$$

where the subscripts indicate partial derivatives. Moreover, in an initially stationary incompressible fluid, there is no vorticity, and the fluid stays irrotational, hence $\backslash nabla\backslash times\backslash textbf\{u\}\text{'}=0\backslash ,$. In the streamfunction representation, $\backslash nabla^2\backslash psi=0.\backslash ,$ Next, because of the translational invariance of the system in the x-direction, it is possible to make the ansatz

$$\backslash psi\backslash left(x,z,t\backslash right)=e^\{i\backslash alpha\backslash left(x-ct\backslash right)\}\backslash Psi\backslash left(z\backslash right),\backslash ,$$

where $\backslash alpha\backslash ,$ is a spatial wavenumber. Thus, the problem reduces to solving the equation

$$\backslash left(D^2-\backslash alpha^2\backslash right)\backslash Psi\_j=0,\backslash ,\backslash ,\backslash ,\backslash \; D=\backslash frac\{d\}\{dz\},\backslash ,\backslash ,\backslash ,\backslash \; j=L,G.\backslash ,$$

The domain of the problem is the following: the fluid with label 'L' lives in the region c, which in turn determines the stability properties of the system.

The first of these conditions is provided by details at the boundary. The perturbation velocities $w\text{'}\_i\backslash ,$ should satisfy a no-flux condition, so that fluid does not leak out at the boundaries $z=\backslash pm\backslash infty.\backslash ,$ Thus, $w\_L\text{'}=0\backslash ,$ on $z=-\backslash infty\backslash ,$, and $w\_G\text{'}=0\backslash ,$ on $z=\backslash infty\backslash ,$. In terms of the streamfunction, this is

$$\backslash Psi\_L\backslash left(-\backslash infty\backslash right)=0,\backslash qquad\; \backslash Psi\_G\backslash left(\backslash infty\backslash right)=0.\backslash ,$$

The other three conditions are provided by details at the interface $z=\backslash eta\backslash left(x,t\backslash right)\backslash ,$.

Continuity of vertical velocity: At $z=\backslash eta$, the vertical velocities match, $w\text{'}\_L=w\text{'}\_G\backslash ,$. Using the stream function representation, this gives

$$\backslash Psi\_L\backslash left(\backslash eta\backslash right)=\backslash Psi\_G\backslash left(\backslash eta\backslash right).\backslash ,$$

Expanding about $z=0\backslash ,$ gives

$$\backslash Psi\_L\backslash left(0\backslash right)=\backslash Psi\_G\backslash left(0\backslash right)+\backslash text\{H.O.T.\},\backslash ,$$

where H.O.T. means 'higher-order terms'. This equation is the required interfacial condition.

The free-surface condition: At the free surface $z=\backslash eta\backslash left(x,t\backslash right)\backslash ,$, the kinematic condition holds:

$$\backslash frac\{\backslash partial\backslash eta\}\{\backslash partial\; t\}+u\text{'}\backslash frac\{\backslash partial\backslash eta\}\{\backslash partial\; x\}=w\text{'}\backslash left(\backslash eta\backslash right).\backslash ,$$

Linearizing, this is simply

$$\backslash frac\{\backslash partial\backslash eta\}\{\backslash partial\; t\}=w\text{'}\backslash left(0\backslash right),\backslash ,$$

where the velocity $w\text{'}\backslash left(\backslash eta\backslash right)\backslash ,$ is linearized on to the surface $z=0\backslash ,$. Using the normal-mode and streamfunction representations, this condition is $c\; \backslash eta=\backslash Psi\backslash ,$, the second interfacial condition.

Pressure relation across the interface: For the case with surface tension, the pressure difference over the interface at $z=\backslash eta$ is given by the Young–Laplace equation:

$$p\_G\backslash left(z=\backslash eta\backslash right)\; -\; p\_L\backslash left(z=\backslash eta\backslash right)\; =\; \backslash sigma\backslash kappa,\backslash ,$$

where σ is the surface tension and κ is the curvature of the interface, which in a linear approximation is

$$\backslash kappa=\backslash nabla^2\backslash eta=\backslash eta\_\{xx\}.\backslash ,$$

Thus,

$$p\_G\backslash left(z=\backslash eta\backslash right)-p\_L\backslash left(z=\backslash eta\backslash right)=\backslash sigma\backslash eta\_\{xx\}.\backslash ,$$

However, this condition refers to the total pressure (base+perturbed), thus

$$\backslash left[P\_G\backslash left(\backslash eta\backslash right)+p\text{'}\_G\backslash left(0\backslash right)\backslash right]-\backslash left[P\_L\backslash left(\backslash eta\backslash right)+p\text{'}\_L\backslash left(0\backslash right)\backslash right]=\backslash sigma\backslash eta\_\{xx\}.\backslash ,$$

(As usual, The perturbed quantities can be linearized onto the surface z=0.) Using hydrostatic balance, in the form

$$P\_L=-\backslash rho\_L\; g\; z+p\_0,\backslash qquad\; P\_G=-\backslash rho\_G\; gz\; +p\_0,\backslash ,$$

this becomes

$$p\text{'}\_G-p\text{'}\_L=g\backslash eta\backslash left(\backslash rho\_G-\backslash rho\_L\backslash right)+\backslash sigma\backslash eta\_\{xx\},\backslash qquad\backslash text\{on\; \}z\; =\; 0.\backslash ,$$

The perturbed pressures are evaluated in terms of streamfunctions, using the horizontal momentum equation of the linearised Euler equations for the perturbations,

$$\backslash frac\{\backslash partial\; u\_i\text{'}\}\{\backslash partial\; t\}\; =\; -\; \backslash frac\{1\}\{\backslash rho\_i\}\backslash frac\{\backslash partial\; p\_i\text{'}\}\{\backslash partial\; x\}\backslash ,$$

with $i=L,G,\backslash ,$ to yield

$$p\_i\text{'}=\backslash rho\_i\; c\; D\backslash Psi\_i,\backslash qquad\; i=L,G.\backslash ,$$

Putting this last equation and the jump condition on $p\text{'}\_G-p\text{'}\_L$ together,

$$c\backslash left(\backslash rho\_G\; D\backslash Psi\_G-\backslash rho\_L\; D\backslash Psi\_L\backslash right)=g\backslash eta\backslash left(\backslash rho\_G-\backslash rho\_L\backslash right)+\backslash sigma\backslash eta\_\{xx\}.\backslash ,$$

Substituting the second interfacial condition $c\backslash eta=\backslash Psi\backslash ,$ and using the normal-mode representation, this relation becomes

$$c^2\backslash left(\backslash rho\_G\; D\backslash Psi\_G-\backslash rho\_L\; D\backslash Psi\_L\backslash right)=g\backslash Psi\backslash left(\backslash rho\_G-\backslash rho\_L\backslash right)-\backslash sigma\backslash alpha^2\backslash Psi,\backslash ,$$

where there is no need to label $\backslash Psi\backslash ,$ (only its derivatives) because $\backslash Psi\_L=\backslash Psi\_G\backslash ,$

at $z=0.\backslash ,$

; Solution

Now that the model of stratified flow has been set up, the solution is at hand. The streamfunction equation $\backslash left(D^2-\backslash alpha^2\backslash right)\backslash Psi\_i=0,\backslash ,$ with the boundary conditions $\backslash Psi\backslash left(\backslash pm\backslash infty\backslash right)\backslash ,$ has the solution

$$\backslash Psi\_L=A\_L\; e^\{\backslash alpha\; z\},\backslash qquad\; \backslash Psi\_G\; =\; A\_G\; e^\{-\backslash alpha\; z\}.\backslash ,$$

The first interfacial condition states that $\backslash Psi\_L=\backslash Psi\_G\backslash ,$ at $z=0\backslash ,$, which forces $A\_L=A\_G=A.\backslash ,$ The third interfacial condition states that

$$c^2\backslash left(\backslash rho\_G\; D\backslash Psi\_G-\backslash rho\_L\; D\backslash Psi\_L\backslash right)=g\backslash Psi\backslash left(\backslash rho\_G-\backslash rho\_L\backslash right)-\backslash sigma\backslash alpha^2\backslash Psi.\backslash ,$$

Plugging the solution into this equation gives the relation

$$Ac^2\backslash alpha\backslash left(-\backslash rho\_G-\backslash rho\_L\backslash right)=Ag\backslash left(\backslash rho\_G-\backslash rho\_L\backslash right)-\backslash sigma\backslash alpha^2A.\backslash ,$$

The A cancels from both sides and we are left with

$$c^2=\backslash frac\{g\}\{\backslash alpha\}\backslash frac\{\backslash rho\_L-\backslash rho\_G\}\{\backslash rho\_L+\backslash rho\_G\}+\backslash frac\{\backslash sigma\backslash alpha\}\{\backslash rho\_L+\backslash rho\_G\}.\backslash ,$$

To understand the implications of this result in full, it is helpful to consider the case of zero surface tension. Then,

$$c^2=\backslash frac\{g\}\{\backslash alpha\}\backslash frac\{\backslash rho\_L-\backslash rho\_G\}\{\backslash rho\_L+\backslash rho\_G\},\backslash qquad\; \backslash sigma=0,\backslash ,$$

and clearly

• If , $c^2>0\backslash ,$ and c is real. This happens when the lighter fluid sits on top;
• If $\backslash rho\_G>\backslash rho\_L\backslash ,$, and c is purely imaginary. This happens when the heavier fluid sits on top.

Now, when the heavier fluid sits on top, , and

$$c=\backslash pm\; i\; \backslash sqrt\{\backslash frac\{g\backslash mathcal\{A\}\}\{\backslash alpha\}\},\backslash qquad\; \backslash mathcal\{A\}=\backslash frac\{\backslash rho\_G-\backslash rho\_L\}\{\backslash rho\_G+\backslash rho\_L\},\backslash ,$$

where $\backslash mathcal\{A\}\backslash ,$ is the Atwood number. By taking the positive solution, we see that the solution has the form

$$\backslash Psi\backslash left(x,z,t\backslash right)\; =\; A\; e^\{-\backslash alpha|z|\}\; \backslash exp\; \backslash left[i\backslash alpha\backslash left(x-ct\backslash right)\backslash right]\; =\; A\; \backslash exp\backslash left(\backslash alpha\backslash sqrt\{\backslash frac\{g\backslash tilde\{\backslash mathcal\{A\}\}\}\{\backslash alpha\}\}t\backslash right)\backslash exp\backslash left(i\backslash alpha\; x\; -\; \backslash alpha|z|\backslash right)\backslash ,$$

and this is associated to the interface position η by: $c\backslash eta=\backslash Psi.\backslash ,$ Now define $B=A/c.\backslash ,$

{{hidden end}}

When the two layers of the fluid are allowed to have a relative velocity, the instability is generalized to the Kelvin–Helmholtz–Rayleigh–Taylor instability, which includes both the Kelvin–Helmholtz instability and the Rayleigh–Taylor instability as special cases. It was recently discovered that the fluid equations governing the linear dynamics of the system admit a parity-time symmetry, and the Kelvin–Helmholtz–Rayleigh–Taylor instability occurs when and only when the parity-time symmetry breaks spontaneously.{{cite journal|last=Qin|first=H.|display-authors=et al.|journal=Physics of Plasmas| date=2019| volume=26| issue=3| title=Kelvin–Helmholtz instability is the result of parity-time symmetry breaking |page=032102|doi=10.1063/1.5088498|arxiv= 1810.11460 | bibcode=2019PhPl...26c2102Q|s2cid=53658729}}

File:Image_for_explanation_of_Rayleigh_Taylor_Instability_by_using_vorticity.xcf

The RT instability can be seen as the result of baroclinic torque created by the misalignment of the pressure and density gradients at the perturbed interface, as described by the two-dimensional inviscid vorticity equation, $\backslash frac\{D\backslash omega\}\{Dt\}\; =\; \backslash frac\{1\}\{\backslash rho\; ^2\}\backslash nabla\; \backslash rho\; \backslash times\; \backslash nabla\; p$, where ω is vorticity, ρ density and p is the pressure. In this case the dominant pressure gradient is hydrostatic, resulting from the acceleration.

When in the unstable configuration, for a particular harmonic component of the initial perturbation, the torque on the interface creates vorticity that will tend to increase the misalignment of the gradient vectors. This in turn creates additional vorticity, leading to further misalignment. This concept is depicted in the figure, where it is observed that the two counter-rotating vortices have velocity fields that sum at the peak and trough of the perturbed interface. In the stable configuration, the vorticity, and thus the induced velocity field, will be in a direction that decreases the misalignment and therefore stabilizes the system.{{cite thesis |degree=PhD |url=| title=Experiments and Simulations on the Incompressible, Rayleigh–Taylor Instability with Small Wavelength Initial Perturbations | publisher= University of Arizona Dissertations | date=2012 | last= Roberts | first=M.S. | bibcode=2012PhDT.......222R |hdl=10150/265355 |hdl-access=free}}

A much simpler explanation of the basic physics of the Rayleigh–Taylor instability was published in 2006.{{Cite journal |last1=Piriz |first1=A. R. |last2=Cortázar |first2=O. D. |last3=López Cela |first3=J. J. |last4=Tahir |first4=N. A. |date=2006-12-01 |title=The Rayleigh-Taylor instability |url=https://pubs.aip.org/aapt/ajp/article-abstract/74/12/1095/1042158/The-Rayleigh-Taylor-instability?redirectedFrom=fulltext |journal=American Journal of Physics |volume=74 |issue=12 |pages=1095–1098 |doi=10.1119/1.2358158 |issn=0002-9505|url-access=subscription }}

The analysis in the previous section breaks down when the amplitude of the perturbation is large. The growth then becomes non-linear as the spikes and bubbles of the instability tangle and roll up into vortices. Then, as in the figure, numerical simulation of the full problem is required to describe the system.

{{clear}}

• Saffman–Taylor instability
• Darrieus–Landau instability
• Richtmyer–Meshkov instability
• Kelvin–Helmholtz instability
• Mushroom cloud
• Plateau–Rayleigh instability
• Salt fingering
• Hydrodynamic stability
• Kármán vortex street
• Fluid thread breakup
• Rayleigh–Bénard convection

{{Reflist}}

• {{cite journal| author=Rayleigh, Lord (John William Strutt) | author-link=John Strutt, 3rd Baron Rayleigh | title=Investigation of the character of the equilibrium of an incompressible heavy fluid of variable density | journal=Proceedings of the London Mathematical Society | volume=14 | pages=170–177 | date=1883 |doi=10.1112/plms/s1-14.1.170 | url=https://zenodo.org/record/1447724 }} (Original paper is available at: https://www.irphe.fr/~clanet/otherpaperfile/articles/Rayleigh/rayleigh1883.pdf .)
• {{cite journal| author=Taylor, Sir Geoffrey Ingram | author-link=Geoffrey Ingram Taylor | title=The instability of liquid surfaces when accelerated in a direction perpendicular to their planes | journal=Proceedings of the Royal Society of London. Series A, Mathematical and Physical Sciences | volume=201 | issue=1065 | pages=192–196 | date=1950 | doi=10.1098/rspa.1950.0052 |bibcode = 1950RSPSA.201..192T | s2cid=98831861 }}

• {{cite book| last=Chandrasekhar | first = Subrahmanyan | author-link=Subrahmanyan Chandrasekhar | title=Hydrodynamic and Hydromagnetic Stability | publisher=Dover Publications | date=1981 | isbn=978-0-486-64071-6 }}
• {{cite book| title=Introduction to hydrodynamic stability | first=P. G. | last=Drazin | author-link=Philip Drazin | publisher=Cambridge University Press | date=2002 | isbn=978-0-521-00965-2 }} xvii+238 pages.
• {{cite book | author= Drazin, P. G. |author2=Reid, W. H. | title= Hydrodynamic stability | url= https://archive.org/details/hydrodynamicstab0000draz_02ed | url-access= registration | date=2004 | publisher=Cambridge University Press | location=Cambridge | isbn=978-0-521-52541-1 |edition=2nd }} 626 pages.

{{Commons category}}

• [https://web.archive.org/web/20050404031045/http://acg.media.mit.edu/people/fry/mixing/ Java demonstration of the RT instability in fluids]
• [http://www.enseeiht.fr/hmf/travaux/CD0001/travaux/optmfn/hi/01pa/hyb72/rt/rt.htm Actual images and videos of RT fingers]
• [http://fluidlab.arizona.edu/ Experiments on Rayleigh–Taylor instability at the University of Arizona]
• [http://media.caltech.edu/press_releases/13496 plasma Rayleigh–Taylor instability experiment at California Institute of Technology]

{{Authority control}}

{{DEFAULTSORT:Rayleigh-Taylor instability}}

Category:Fluid dynamics

Category:Fluid dynamic instabilities

Category:Plasma instabilities