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. Linear stability analysis 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. 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 where is the temporal growth rate, is the spatial wavenumber and is the Atwood number. 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 is no vorticity, and the fluid stays 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 determines the stability 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: