next up previous
Next: Taylor-Couette Instability Up: No Title Previous: Introduction

Rayleigh's Inviscid Instability Argument

To study the stability of inviscid circular flow between concentric cylinders, Lord Rayleigh developed an analogy with the stability of a fluid of variable density under the force of gravity. Taylor summarizes this as: ``the varying centrifugal force of the different layers of fluid plays the part of gravity and the resulting condition for stability is that the square of the circulation must increase continuously in passing from the inner to the outer cylinder, just as the density of a fluid under gravity must decrease continuously with height in order that it may be in stable equilibrium.''gif

To examine this argument, consider a ring of fluid in a flow field with a velocity distribution tex2html_wrap_inline823 . The angular momentum per unit volume measured about the central axis is given by tex2html_wrap_inline825 . Expressing the velocity tex2html_wrap_inline823 in terms of the angular velocity, tex2html_wrap_inline829 , gives tex2html_wrap_inline831 , and therefore

equation70

One can also consider the circulation of the fluid along a circular path,

equation74

and notice that the circulation has the dimensions of angular momentum, tex2html_wrap_inline833 .

In an inviscid fluid, the angular momentum of a small fluid element is conserved. This can be seen from the Euler equations for inviscid flow,

equation78

We assume the motion is axisymmetric, so that the tex2html_wrap_inline835 -component of the Euler equations becomes

equation84

Multiplication by r leads to

eqnarray94

This is simply the material derivative of tex2html_wrap_inline839 , so that

equation107

and we see that tex2html_wrap_inline839 or tex2html_wrap_inline833 is constant and angular momentum is conserved.

Next we consider the perturbation resulting when two rings of fluid at different radii are interchanged while conserving angular momentum. The change in angular momentum is

equation114

from which it follows that the corresponding change in angular velocity for a displaced fluid ring is

  equation116

We now consider the stability of such a perturbation. The rotating fluid is stable if the displaced fluid ring experiences a force which tends to return it to its original position. Thus, we can imagine displacing a small fluid element from r to tex2html_wrap_inline847 such that angular momentum is conserved, and then asking whether the pressure in the surrounding fluid is able to maintain equilibrium. We seek a comparison between the pressure in the equilibrium state, and the pressure in the disturbed state. If the pressure in the undisturbed case is greater than in the disturbed state, it will tend to push the disturbed ring of fluid back to its equilibrium position and maintain stability in the flow. If, on the other hand, the pressure of the disturbed ring of fluid is less than the pressure in the undisturbed fluid surrounding it, instability will set in.

The r-component of the Euler equation gives the centrifugal pressure gradient as

equation121

so that tex2html_wrap_inline851 is the pressure to within a constant. For fluid elements in the equilibrium state, the angular velocity at r is tex2html_wrap_inline855 , and at tex2html_wrap_inline847 it is tex2html_wrap_inline859 , or tex2html_wrap_inline861 to first order. The pressure in the undisturbed flow at tex2html_wrap_inline847 is then

  equation133

In the disturbed flow case, the angular velocity at r is tex2html_wrap_inline855 and at tex2html_wrap_inline847 it is tex2html_wrap_inline871 . Since angular momentum is conserved, tex2html_wrap_inline873 (from Equation 11), so we then have

equation144

The pressure in the disturbed flow at tex2html_wrap_inline847 is therefore given by

  equation148

Taking the difference between Eqs. 13 and 15 we obtain

equation158

so that

equation168

Stability requires that tex2html_wrap_inline877 , which implies the following conditions:

eqnarray182

If we think of tex2html_wrap_inline833 as the circulation of the flow, we come to Rayleigh's conclusion that ``A circulation always increasing outwards ... ensures stability.''

Returning to the case of Couette flow, we have tex2html_wrap_inline881 so that tex2html_wrap_inline883 . Substituting the value given for A in Equation 2, we require that for stability

equation194

Since tex2html_wrap_inline887 by definition, we must have tex2html_wrap_inline889 for stable flow, or

equation200

In tex2html_wrap_inline891 -space, the stable region is therefore bounded by the positive tex2html_wrap_inline817 -axis and the line tex2html_wrap_inline895 . From this we see that if the outer cylinder is fixed ( tex2html_wrap_inline897 ), the flow is always unstable.

When tex2html_wrap_inline899 and the cylinders are rotating in opposite directions, instability is predicted for all rotation speeds in the absence of viscosity. As explained by Chandrasekhargif, however, the instability is only partial. It occurs between the interior cylinder and the nodal surface, defined by tex2html_wrap_inline901 , which separates the flow into two parts. The radius of this surface is given by

equation207

so there is stability for tex2html_wrap_inline903 and instability for tex2html_wrap_inline905 . Thus, according to Taylor, ``In the inner region the square of the circulation decreases outwards, so that centrifugal force tends to make the flow unstable. In the outer region the square of the circulation increases so that centrigual force tends to make the flow stable''gif (p. 69). It is interesting to note that as tex2html_wrap_inline907 , tex2html_wrap_inline909 , and the area where the fluid is unstable in the inviscid case becomes infinitesimally small.

When viscous effects are included, there is an added area of stability in tex2html_wrap_inline891 -space due to viscous stabilization. Stable flow with tex2html_wrap_inline897 , for instance, is possible for tex2html_wrap_inline815 less than some critical value. In general, flow is unstable only if

equation214

Taylor was the first to investigate this case theoretically and experimentally, and determine the new criterion for the onset of instability. His experiments confirmed his findings, and showed that the marginal state was stationary and showed a steady axisymmetric cellular patterngif. He found that when the two cylinders rotate in the same direction, Rayleigh's criterion for inviscid flow still holds, with a small added region of stability. When the cylinders rotate in opposite directions, however, Taylor found a large region of stable flow possible, unlike Rayleigh's prediction that all rotation speeds lead to instability.


next up previous
Next: Taylor-Couette Instability Up: No Title Previous: Introduction

brenner@math.mit.edu