Accurate solution of the Orr–Sommerfeld stability equation

The Orr-Sommerfeld equation is solved numerically using expansions in Chebyshev polynomials and the QR matrix eigenvalue algorithm. It is shown that results of great accuracy are obtained very economically. The method is applied to the stability of plane Poiseuille flow; it is found that the critical Reynolds number is 5772·22. It is explained why expansions in Chebyshev polynomials are better suited to the solution of hydrodynamic stability problems than expansions in other, seemingly more relevant, sets of orthogonal functions.