Rayleigh–Taylor instability
![](http://upload.wikimedia.org/wikipedia/commons/thumb/0/00/Crab_Nebula.jpg/200px-Crab_Nebula.jpg)
The Rayleigh–Taylor instability, or RT instability (after Lord Rayleigh and G. I. Taylor), occurs any time a dense, heavy fluid is being accelerated into a light fluid. This is the case with a cloud and shock system, or when a fluid of a certain density floats above a fluid of lesser density, such as dense oil floating above water.
Two completely plane-parallel layers of immiscible fluid are in equilibrium, but the slightest perturbation leads to release of potential energy, as the heavier material moves down under the (effective) gravitational field, and the lighter material is displaced upwards. As the instability develops, downward-moving irregularities ('dimples') are quickly magnified into sets of inter-penetrating RT fingers. The upward-moving, lighter material behaves like 'Spherical Cap Bubbles'.
This process is evident not only in many terrestrial examples, from salt domes to weather inversions, but also in astrophysics and electrohydrodynamics. RT fingers are especially obvious 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.
Note that the RT instability is not to be confused with the Rayleigh instability (or Plateau-Rayleigh instability) of a liquid jet. This latter 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 volume but lower surface area.
Linear stability analysis
The inviscid two-dimensional Rayleigh--Taylor (RT) instability provides an excellent springboard into the mathematical study of stability because of the exceptionally simple nature of the base state. This is the equilibrium state that exists before any perturbation is added to the system, and is described by the velocity field , where the gravitational
field is . An interface at separates the fluids of densities in the upper region, and in the lower region. In this section it is shown that when the heavy fluid sits on top, the growth of a small perturbation at the interface is exponential, and takes place at the rate[1]
- .
The perturbation introduced to the system is described by a velocity field of infinitesimally small amplitude, . Because the fluid is assumed incompressible, this velocity field has the streamfunction representation
- ,
where the subscripts indicate partial derivatives. Moreover, in an initially stationary incompressible fluid, there are no vortices, and the fluid is irrotational, hence . In the streamfunction representation, . Next, because of the translational invariance of the system in the x-direction, it is possible to make the ansatz
- ,
where is a spatial wavenumber. Thus, the problem reduces to solving the equation
The domain of the problem is the following: the fluid with label `L' lives in the region , while the fluid with the label `G' lives in the upper half-plane . To specify the solution fully, it is necessary to fix conditions at the boundaries and interface. This determines the wave speed c, which in turn determiens the stabilty properties of the system.
The first of these conditions is provided by details at the boundary. The perturbation velocities should satisfy a no-flux condition, so that fluid does not leak out at the boundaries . Thus, on , and on . In terms of the streamfunction, this is
The other three conditions are provided by details at the interface .
Continuity of vertical velocity: At , the vertical velocities match, . Using the streamfunction representation, this gives
Expanding about gives
where H.O.T. means `higher-order terms'. This equation is the required interfacial condition.
The free-surface condition: At the free surface , the kinematic condition holds:
Linearizing, this is simply
where the velocity is linearized on to the surface . Using the normal-mode and streamfunction representations, this condition is , the second interfacial condition.
Pressure relation across the interface: To derive this condition, take the full Navier--Stokes equation
and dot it with the normal vector to obtain the relation
where is the force acting normal to the interface (, by Newton's second law), and is the derivative in a direction normal to the interface. Furthermore, is the surface tension coefficient, is the curvature of the interface, which in a linear approximation is
and is the Dirac delta function, which localizes the surface tension force to the interface.
Next, integrate this equation across the interface, from to (now the interface is at in the n-coordinate system).
For the force acting across the interface to be finite, the left-hand side vanishes as , and similarly for the gravitational term on the right-hand side. Thus,
which in statics is called the Young-Laplace condition. However, this condition refers to the total pressure (base+perturbed), thus
(As usual, The perturbed quantities can be linearized onto the surface z=0.) Using hydrostatic balance, in the form
this becomes
The perturbed pressures are evaluated in terms of streamfunctions to yield
Putting this last equation and the jump condition together,
Substituting the second interfacial condition and using the normal-mode representation, this relation becomes
where there is no need to label (only its derivatives) because at .
Solution
Now that the model of stratified flow has b een set up, the solution is at hand. The streamfunction equation , with the boundary conditions has the solution
The first interfacial condition states that at , which forces . The third interfacial condition states that
Plugging the solution into this equation gives the relation
The A cancels from both sides and we are left with
To understand the implications of this result in full, it is helpful to consider the case of zero surface tension. Then,
![](http://upload.wikimedia.org/wikipedia/commons/thumb/d/d6/HD-Rayleigh-Taylor.gif/400px-HD-Rayleigh-Taylor.gif)
and clearly
- If , and c is real. This happens when the
lighter fluid sits on top;
- If , and c is purely imaginary. This happens
when the heavier fluid sits on top.
Now, when the heavier fluid sits on top, , and
where is the unsigned Atwood number. By taking the positive solution, we see that the solution has the form
which grows exponentially in time. In general, the condition for linear instability is that the imaginary part of the wave speed be positive. Finally, restoring the surface tension makes less negative and is therefore stabilizing. Indeed, there is a range of short waves for which the surface tension stabilizes the system and prevents the instability forming.
Late-time behaviour
The analysis of the previous section breaks down when the amplitude of the perturbation is large. Then, as in the figure, numerical simulation of the full problem is required to describe the system.
Footnotes
- ^ A similar derivation appears in the book by S. Chandrasekhar, 'Hydrodynamic and Hydromagnetic Stability', Dover Publications (1981).
- ^ Li, Shengtai and Hui Li. "Parallel AMR Code for Compressible MHD or HD Equations". Los Alamos National Laboratory. Retrieved 2006-09-05.
See also
- Richtmyer-Meshkov instability
- Kelvin–Helmholtz instability
- Mushroom cloud
- Plateau-Rayleigh instability
- Salt fingering
- Kármán vortex street
External links
- Java demonstration of the RT instability in fluids
- Actual images and videos of RT fingers
- Experiments on Rayleigh-Taylor experiments at the University of Arizona
References
Original research papers
- Rayleigh, Lord (John William Strutt) (1883). "Investigation of the character of the equilibrium of an incompressible heavy fluid of variable density". Proceedings of the London Mathematical Society. 14: 170–177. doi:10.1112/plms/s1-14.1.170. (Original paper is available at: https://www.irphe.univ-mrs.fr/~clanet/otherpaperfile/articles/Rayleigh/rayleigh1883.pdf .)
- Taylor, Sir Geoffrey Ingram (1950). "The instability of liquid surfaces when accelerated in a direction perpendicular to their planes". Proceedings of the Royal Society of London. Series A, Mathematical and Physical Sciences. 201 (1065): 192–196. doi:10.1098/rspa.1950.0052.