Rayleigh–Taylor instability

From Wikipedia, the free encyclopedia
Hydrodynamics simulation of a single "finger" of the Rayleigh–Taylor instability.[1] Note the formation of Kelvin–Helmholtz instabilities, in the second and later snapshots shown (starting initially around the level ), as well as the formation of a "mushroom cap" at a later stage in the third and fourth frame in the sequence.
RT instability fingers evident in the Crab Nebula

The Rayleigh–Taylor instability, or RT instability (after Lord Rayleigh and G. I. Taylor), is an instability of an interface between two fluids of different densities which occurs when the lighter fluid is pushing the heavier fluid.[2][3][4] Examples include the behavior of water suspended above oil in the gravity of Earth,[3] mushroom clouds like those from volcanic eruptions and atmospheric nuclear explosions,[5] supernova explosions in which expanding core gas is accelerated into denser shell gas,[6][7] instabilities in plasma fusion reactors and[8] inertial confinement fusion.[9]

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 more dense 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 more dense 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.[3] 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 more dense fluid.[3] This occurs deep underwater on the surface of an expanding bubble and in a nuclear explosion.[10]

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.[2]

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.[11] 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.[12] This latter case resembles magnetically modulated RT instabilities.[13][14] [15]

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.

Stages of development and eventual evolution into turbulent mixing[]

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.[2] 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.[16]

Linear stability analysis[]

Base state of the Rayleigh–Taylor instability. Gravity points downwards.

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.[17] This is the equilibrium state that exists before any perturbation is added to the system, and is described by the mean 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[3]

where is the temporal growth rate, is the spatial wavenumber and is the Atwood number.

Details of the linear stability analysis[17] A similar derivation appears in,[13] §92, pp. 433–435.

The time evolution of the free interface elevation initially at is given by:

which grows exponentially in time. Here B is the amplitude of the initial perturbation, and denotes the real part of the complex valued expression between brackets.

In general, the condition for linear instability is that the imaginary part of the "wave speed" c is positive. Finally, restoring the surface tension makes c2 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.

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.[18]

Vorticity explanation[]

Visualization of an unstable Rayleigh–Taylor instability configuration where baroclinic torque at the interface creates vorticity and induces a velocity field that increases the baroclinic torque. Here ω is vorticity, p is pressure, ρ is density, u is velocity and g is gravity. The thick circular arrows represent the velocity field created by the vortex.

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, , 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.[16][19]

Late-time behaviour[]

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.

See also[]

Notes[]

  1. ^ Li, Shengtai & Hui Li. "Parallel AMR Code for Compressible MHD or HD Equations". Los Alamos National Laboratory. Retrieved 2006-09-05.
  2. ^ Jump up to: a b c Sharp, D.H. (1984). "An Overview of Rayleigh-Taylor Instability". Physica D. 12 (1): 3–18. Bibcode:1984PhyD...12....3S. doi:10.1016/0167-2789(84)90510-4.
  3. ^ Jump up to: a b c d e Drazin (2002) pp. 50–51.
  4. ^ David Youngs (ed.). "Rayleigh–Taylor instability and mixing". Scholarpedia.
  5. ^ "Why Nuclear Bombs Create Mushroom Clouds".
  6. ^ Wang, C.-Y. & Chevalier R. A. (2000). "Instabilities and Clumping in Type Ia Supernova Remnants". The Astrophysical Journal. 549 (2): 1119–1134. arXiv:astro-ph/0005105v1. Bibcode:2001ApJ...549.1119W. doi:10.1086/319439. S2CID 15244583.
  7. ^ Hillebrandt, W.; Höflich, P. (1992). "Supernova 1987a in the Large Magellanic Cloud". In R. J. Tayler (ed.). Stellar Astrophysics. CRC Press. pp. 249–302. ISBN 978-0-7503-0200-5.. See page 274.
  8. ^ Chen, H. B.; Hilko, B.; Panarella, E. (1994). "The Rayleigh–Taylor instability in the spherical pinch". Journal of Fusion Energy. 13 (4): 275–280. Bibcode:1994JFuE...13..275C. doi:10.1007/BF02215847. S2CID 122223176.
  9. ^ Betti, R.; Goncharov, V.N.; McCrory, R.L.; Verdon, C.P. (1998). "Growth rates of the ablative Rayleigh–Taylor instability in inertial confinement fusion". Physics of Plasmas. 5 (5): 1446–1454. Bibcode:1998PhPl....5.1446B. doi:10.1063/1.872802.
  10. ^ John Pritchett (1971). "EVALUATION OF VARIOUS THEORETICAL MODELS FOR UNDERWATER EXPLOSION" (PDF). U.S. Government. p. 86. Retrieved October 9, 2012.
  11. ^ Hester, J. Jeff (2008). "The Crab Nebula: an Astrophysical Chimera". Annual Review of Astronomy and Astrophysics. 46: 127–155. Bibcode:2008ARA&A..46..127H. doi:10.1146/annurev.astro.45.051806.110608.
  12. ^ Berger, Thomas E.; Slater, Gregory; Hurlburt, Neal; Shine, Richard; et al. (2010). "Quiescent Prominence Dynamics Observed with the Hinode Solar Optical Telescope. I. Turbulent Upflow Plumes". The Astrophysical Journal. 716 (2): 1288–1307. Bibcode:2010ApJ...716.1288B. doi:10.1088/0004-637X/716/2/1288.
  13. ^ Jump up to: a b Chandrasekhar, S. (1981). Hydrodynamic and Hydromagnetic Stability. Dover. ISBN 978-0-486-64071-6.. See Chap. X.
  14. ^ Hillier, A.; Berger, Thomas; Isobe, Hiroaki; Shibata, Kazunari (2012). "Numerical Simulations of the Magnetic Rayleigh-Taylor Instability in the Kippenhahn-Schl{\"u}ter Prominence Model. I. Formation of Upflows". The Astrophysical Journal. 716 (2): 120–133. Bibcode:2012ApJ...746..120H. doi:10.1088/0004-637X/746/2/120.
  15. ^ Singh, Chamkor; Das, Arup K.; Das, Prasanta K. (2016), "Single-mode instability of a ferrofluid-mercury interface under a nonuniform magnetic field", Physical Review E, 94 (1): 012803, Bibcode:2016PhRvE..94a2803S, doi:10.1103/PhysRevE.94.012803, PMID 27575198
  16. ^ Jump up to: a b Roberts, M.S.; Jacobs, J.W. (2015). "The effects of forced small-wavelength, finite-bandwidth initial perturbations and miscibility on the turbulent Rayleigh Taylor instability". Journal of Fluid Mechanics. 787: 50–83. Bibcode:2016JFM...787...50R. doi:10.1017/jfm.2015.599. OSTI 1436483. S2CID 126143908.
  17. ^ Jump up to: a b Drazin (2002) pp. 48–52.
  18. ^ Qin, H.; et al. (2019). "Kelvin-Helmholtz instability is the result of parity-time symmetry breaking". Physics of Plasmas. 26 (3): 032102. arXiv:1810.11460. Bibcode:2019PhPl...26c2102Q. doi:10.1063/1.5088498. S2CID 53658729.
  19. ^ Roberts, M.S. (2012). "Experiments and Simulations on the Incompressible, Rayleigh-Taylor Instability with Small Wavelength Initial Perturbations". University of Arizona Dissertations. Bibcode:2012PhDT.......222R. Cite journal requires |journal= (help)

References[]

Original research papers[]

Other[]

External links[]

Retrieved from ""