Orr–Sommerfeld equation: Difference between revisions
No edit summary |
|||
Line 65: | Line 65: | ||
<math>\left|\begin{array}{cccc}\chi_1\left(0\right)&\chi_2\left(0\right)&\chi_3\left(0\right)&\chi_4\left(1\right)\\ |
<math>\left|\begin{array}{cccc}\chi_1\left(0\right)&\chi_2\left(0\right)&\chi_3\left(0\right)&\chi_4\left(1\right)\\ |
||
\chi_1'\left(0\right)&\chi_2'\left(0\right)&\chi_3'\left(0\right)&\chi_4'\left(0\right)\\ |
\chi_1'\left(0\right)&\chi_2'\left(0\right)&\chi_3'\left(0\right)&\chi_4'\left(0\right)\\ |
||
\Omega_1\left(1\right)&\Omega_2\left(1\right)&\Omega_3\left(1\right)&\Omega_4\left( |
\Omega_1\left(1\right)&\Omega_2\left(1\right)&\Omega_3\left(1\right)&\Omega_4\left(1\right)\\ |
||
\chi_1''\left(1\right)+\alpha^2\chi_1\left(1\right)&\chi_2''\left(1\right)+\alpha^2\chi_2\left(1\right)&\chi_3''\left(1\right)+\alpha^2\chi_3\left(1\right)&\chi_4''\left(1\right)+\alpha^2\chi_4\left(1\right)\end{array}\right|=0 |
\chi_1''\left(1\right)+\alpha^2\chi_1\left(1\right)&\chi_2''\left(1\right)+\alpha^2\chi_2\left(1\right)&\chi_3''\left(1\right)+\alpha^2\chi_3\left(1\right)&\chi_4''\left(1\right)+\alpha^2\chi_4\left(1\right)\end{array}\right|=0 |
||
</math> |
</math> |
Revision as of 09:38, 28 June 2008
The Orr-Sommerfeld equation, in fluid dynamics, is an eigenvalue equation describing the linear two-dimensional modes of disturbance to a viscous parallel flow. A parallel, laminar flow that solves the Navier-Stokes equations can become unstable if certain conditions on the flow are satisfied, and the Orr-Sommerfeld equation determines precisely what these conditions are.
Formulation
![](http://upload.wikimedia.org/wikipedia/commons/thumb/c/c3/Os_schematic.png/300px-Os_schematic.png)
The equation is derived by solving a linearized version of the Navier-Stokes equation for the perturbation velocity field
- ,
where is the unperturbed or basic flow. The perturbation velocity has the wave-like solution (real part understood). Using this knowledge, and the streamfunction representation for the flow, the following dimensional form of the Orr-Sommerfeld equation is obtained:
- ,
where is the dynamic viscosity of the fluid, is its density, and is the streamfunction. The equation can be written in non-dimensional form by measuring velocities according to a scale set by some characteristic velocity , and by measuring lengths according to channel depth . Then the equation takes the form
- ,
where
is the Reynolds number of the basic flow. The relevant boundary conditions are the no-slip boundary conditions at the channel top and bottom and ,
- at and ,
The eigenvalue parameter of the problem is and the eigenvector is . If the imaginary part of the wave speed is positive, then the base flow is unstable, and the small perturbation introduced to the system is amplified in time.
Solutions
For all but the simplest of velocity profiles , numerical or asymptotic methods are required to calculate solutions. Some typical flow profiles are discussed below.
![](http://upload.wikimedia.org/wikipedia/commons/thumb/5/57/Spectrum_OS.jpg/300px-Spectrum_OS.jpg)
![](http://upload.wikimedia.org/wikipedia/commons/thumb/9/9c/Poiseuille_dispersion_curve_OS.jpg/300px-Poiseuille_dispersion_curve_OS.jpg)
For Poiseuille flow, it has been shown that the flow is unstable (i.e. one or more eigenvalues has a positive imaginary part) for some when and the neutrally stable mode at having , . [1] To see the stability properties of the system, it is customary to plot a dispersion curve, that is, a plot of the growth rate as a function of the wavenumber .
The first figure shows the spectrum of the Orr--Sommerfeld equation at the critical values listed above. This is a plot of the eigenvalues (in the form ) in the complex plane. The rightmost eigenvalue is the most unstable one. At the critical values of Reynolds number and wavenumber, the rightmost eigenvalue is exactly zero. For higher (lower) values of Reynolds number, the rightmost eigenvalue shifts into the positive (negative) half of the complex plane. Then, a fuller picture of the stability properties is given by a plot exhibiting the functional dependence of this eigenvalue; this is shown in the second figure.
On the other hand, the spectrum of eigenvalues for Couette flow indicates stability, at all Reynolds numbers [2]. However, in experiments, Couette flow is found to be unstable to small perturbations. This is explained through the non-normality of the eigenvalue problem associated with Couette (and indeed, Poiseuille) flow [3]. That is, the eigenfunctions of the operator are complete but non-orthogonal. Thus, if a number of modes are excited, it is possible for the system energy to be transferred to the most dangerous of these modes. This cascade of energy destabilizes the system. Nevertheless, in this case, the Orr--Sommerfeld analysis still demonstrates that uni-model disturbances are stable.
Mathematical methods for free-surface flows
For Couette flow, it is possible to make mathematical progress in the solution of the Orr--Sommerfeld equation. In this section, a demonstration of this method is given for the case of free-surface flow, that is, when the upper lid of the channel is replaced by a free surface. Note first of all that it is necessary to modify upper boundary conditions to take account of the free surface. In non-dimensional form, these conditions now read
at ,
, at .
The first free-surface condition is the statement of continuity of tangential stress, while the second condition relates the normal stress to the surface tension. Here
are the surface tension and gravity parameters respectively (or the Froude and inverse Weber numbers.)
For Couette flow , the four linearly independent solutions to the non-dimensional Orr--Sommerfeld equation are [4],
- ,
where is the Airy function of the first kind. Substitution of the superposition solution into the four boundary conditions gives four equations in the four unknown constants . For the equations to have a non-trivial solution, the determinant condition
must be satisfied. This is a single equation in the unknown c, which can be solved numerically or by asymptotic methods. It can be shwon that for a range of wavenumbers and for sufficiently large Reynolds numbers, the growth rate is positive.
References
- ^ Orszag S. A. (1971) 'Accurate solution of the Orr-Sommerfeld stability equation' J. Fluid. Mech. 50, 689-703
- ^ P. G. Drazin and W. H. Reid (1981) 'Hydrodynamic Stability' Cambridge University Press
- ^ N. L. Trefethen, A. E. Trefethen, S. C. Teddy and T. A. Driscoll (1993) 'Hydrodynamic stability without eigenvalues' Science, 261, 578-584
- ^ R. Miesen and B. J. Boersma (1995) 'Hydrodynamic stability of a sheared liquid film' Journal of Fluid Mechanics, 301, 175-202