The hydrodynamic instability of a viscous incompressible flow with a f
ree surface is studied both numerically and experimentally. While the
free-surface flow is basically two-dimensional at low Reynolds numbers
, a three-dimensional secondary flow pattern similar to the Taylor vor
ticies between two concentric cylinders appears at higher rotational s
peeds. The secondary flow has periodic velocity components in the axia
l direction and is characterized by a distinct spatially periodic vari
ation in surface height similar to a standing wave. A numerical method
, using boundary-fitted coordinates and multigrid methods to solve the
Navier Stokes equations in primitive variables, is developed to treat
two-dimensional free-surface flows. A similar numerical technique is
applied to the linearized three-dimensional perturbation equations to
treat the onset of secondary flows. Experimental measurements have bee
n obtained using light sheet techniques to visualize the secondary flo
w near the free surface. Photographs of streak lines were taken and co
mpared to the numerical calculations. It has been shown that the solut
ion of the linearized equations contains most of the important feature
s of the nonlinear secondary flows at Reynolds number higher than the
critical value. The experimental results also show that the numerical
method predicts well the onset of instability in terms of the critical
wavenumber and Reynolds number.