Perturbations at a two-fluid interface

figure19243
Figure 10.7: An incompressible fluid of density tex2html_wrap_inline12965 overlays a incompressible second fluid of density tex2html_wrap_inline12967 . The interface is distorted so that a fluid element of tex2html_wrap_inline12969 is raised a distance tex2html_wrap_inline12971 . If tex2html_wrap_inline12973 it is energetically favorable for the perturbation to grow.

We have all noticed how the surface of a lake gets ruffled by a gentle breeze. To understand such phenomena, we have to study perturbations at an interface between two fluids. Let us consider a plane horizontal interface separating two fluids, with gravity acting downward. If this interface is disturbed, we expect that, under certain circumstances, the interface may oscillate in the form of waves and, under other circumstances, the disturbance may grow to give rise to an instability.

We shall first develop a general mathematical analysis to show how a perturbation at such an interface evolves. Then several special cases of waves and instabilities will follow from the general dispersion relation we derive. To simplify the analysis, we assume that the fluids on both sides of the interface are incompressible and ideal. Hence, if there is no vorticity inside one of these fluids at the beginning, Kelvin's vorticity theorem asserts that the velocities induced inside the fluid as a result of the perturbation remain irrotational. So velocity potentials can be introduced. We neglect the surface tension at the interface between the two fluids.

We have previously derived Bernoulli's principle for steady flows. We now obtain a similar result for a flow which is incompressible and irrotational, but varying in time. Writing tex2html_wrap_inline19148 and substituting in the Euler equation, we get

equation6555

where tex2html_wrap_inline13075 is the the potential for the body force (gravity in the present case). This equation gives rise to the integral

equation6562

where F(t) is constant in space, but can be a function of time. It may be noted that in deriving Bernoulli's theorem we did not restrict ourselves to irrotational flows. Hence, in order to obtain Bernoulli's theorem, we had to integrate along streamlines to make tex2html_wrap_inline19154 zero. In the present case, tex2html_wrap_inline13746 by the assumption so that it is not necessary to integrate along a streamline. Hence Eq (10.60) should ho1d between points not lying on a streamline.

figure19244
Figure 10.8: Clouds over Denver, Co. The pattern is similar to that seen in the Kelvin-Helmholtz instability. This instability that accompanies wind shear occurs frequently in the atmosphere, however those currents are not normally visible. This is called clear air turbulence, and is a significant hazard to the flight of airplanes.

Figure 10.7 shows the configuration we want to study. The horizontal plane z = 0 separates the two fluids with densities tex2html_wrap_inline12967 (fluid below) and tex2html_wrap_inline12965 (fluid above). Let us also assume that the two fluids have uniform velocities U and U' in the x-direction. It is not difficult to show that this configuration satisfies steady state equations (i.e. hydrodynamic equations with tex2html_wrap_inline18784 ). Let the interface between the two fluids be perturbed from its undisturbed position z = 0 and let us denote the position of the interface by tex2html_wrap_inline19174 . The subscript 1 is to follow the convention that all perturbed quantities have subscript 1. Our aim is to find out whether this perturbation tex2html_wrap_inline19176 grows with time or decays or oscillates. The velocity potential in the fluid below can be written as

equation6578

where the unperturbed part -Ux would give the uniform velocity U in the x direction and tex2html_wrap_inline19184 is the perturbed part satisfying

equation6581

as a result of the incompressibility condition tex2html_wrap_inline13641 . Similarly the velocity potential in the fluid above can be written as

equation6584

where tex2html_wrap_inline19188 is the perturbed part satisfying

equation6587

Since the velocity perturbations are caused by the displacements of the interface, we have to connect the perturbed parts of the velocity potential with the displacements of the interface. To do this, consider a fluid element of the lower fluid lying infinitesimally close to the interface and find out its vertical velocity. In terms of the velocity potential of the lower fluid, the vertical velocity is given by tex2html_wrap_inline19190 . On the other hand, the vertical velocity of the fluid element is also given by the Lagrangian derivative of the displacement tex2html_wrap_inline19192 . Equating these two, we have

equation6590

where the RHS is just the expansion of the Lagrangian derivative, in which we have kept only the terms linear in the perturbed quantities. Similarly, consider a fluid element infinitesimally above the interface and obtain

equation6600

Since we are linearizing the perturbation equations, we can again write any arbitrary perturbation as a superposition of Fourier components. Because of the symmetry in the x-direction, a Fourier component of the displacement tex2html_wrap_inline19196 can be written as

equation6610

The corresponding Fourier components of tex2html_wrap_inline19184 and tex2html_wrap_inline19188 should have similar x and t dependences. The z dependences should be such that the Laplace equations (10.62) and (10.64) are satisfied. Hence, they have to be of the form

equation6615

equation6618

where the signs before kz have been chosen so that the perturbations vanish as we go far away from the interface. Substituting Eq (10.67-10.69) into Eq (10.65-10.66), we get

equation6625

equation6628

These are two equations relating the amplitudes of perturbed quantities A, C and C'. We need a third equation to solve the problem, since we are dealing with three quantities. This is provided by the condition that the pressure has to be continuous across the interface. From Eq (10.60), the pressure inside the lower fluid at a point infinitesimally close to the interface can be written as

equation6632

where we have written tex2html_wrap_inline19216 for the gravitational potential. Writing a similar expression for the pressure infinitesimally above the interface and equating them,

equation6638

Although tex2html_wrap_inline19218 can in principle be a function of time, here it has to be a constant due to the boundary condition that the perturbations vanish far away from the interface at all times. We can find K by considering the unperturbed configuration for which u = U, u'= U' and tex2html_wrap_inline19184tex2html_wrap_inline19188tex2html_wrap_inline19196 are all zero. This gives

equation6648

We note that

equation6653

if we keep only the linear terms in the perturbed quantities. Substituting for tex2html_wrap_inline19232 from Eq (10.75) and making a similar substitution for tex2html_wrap_inline19234 , we find from Eq (10.73) that

equation6660

after making use of Eq (10.74). Substituting from Eq (10.67-10.69) into this equation, we finally get

equation6675

which is the third equation connecting the amplitudes of perturbed quantities. Combining Eq (10.70), (10.71) and (10.77), we obtain

equation6681

An equation connecting the frequency and wavenumber of a Fourier component of the perturbation is called a dispersion relation. For the two-fluid interface problem, Eq (10.78) then is the dispersion relation. We note that Eq (10.78) is a quadratic equation for tex2html_wrap_inline13820 if k is regarded as given. The solution of this quadratic equation is

equation6686
 
 

Kelvin-Helmholtz instability

We now consider what happens if U and U' are not zero. Let us assume tex2html_wrap_inline19248 so that the system is Rayleigh-Taylor stable. If the expression within the square root in Eq (10.79) is negative, then we note that tex2html_wrap_inline13820 has an imaginary part (we are again taking k to be real).

We have already seen that a positive imaginary part in tex2html_wrap_inline13820 leads to instability. Hence the instability occurs when the expression within the square root in Eq (10.79) is negative, i.e. when

equation6733

This instability is known as the Kelvin-Helmholtz instability (Helmholtz 1868; Kelvin 1871) and causes an interface between two fluids to wrinkle if the fluids on the two sides are moving with different speeds. The most common example of the Kelvin-Helmholtz instability is provided by the observation that a wind blowing over a water surface causes the water surface to undulate.

Even when |U - U'| is arbitrarily small, we notice from Eq (10.82) that the Fourier components satisfying

equation6738

become unstable. Since large k means small wavelengths, we conclude that only perturbations with small wavelengths are unstable if | U - U'| is sufficiently small. This conclusion, however, gets completely modified when we include surface tension, because surface tension tends to suppress very small wavelengths. On inclusion of surface tension, it is found that the Kelvin-Helmholtz instability takes place only when |U - U'| is larger than a critical value.

Back to m2.gif - 1967 Bytes or back to 
 

james r. graham 649

Tuesday June 13 09:09:24 PDT 2000