Simulation of Taylor-Couette flow. Part 1. Numerical methods and comparison with experiment

We present a numerical method that allows us to solve the Navier-Stokes equation with boundary conditions for the viscous flow between two concentrically rotating cylinders as an initial-value problem. We use a pseudospectral code in which all of the time-splitting errors are removed by using a set of Green functions (capacitance matrix) that allows us to satisfy the inviscid boundary conditions exactly. For this geometry we find that a small time-splitting error can produce large errors in the computed velocity field. We test the code by comparing our numerically determined growth rates and wave speeds with linear theory and by comparing our computed torques with experimentally measured values and with the values that appear in other published numerical simulations. We find good agreement in all of our tests of the numerical calculation of wavy vortex flows. A test that is more sensitive than the comparison of torques is the comparison of the numerically computed wave speed with the experimentally observed wave speed. The agreements between the simulated and measured wave speeds are within the experimental uncertainties; the best-measured speeds have fractional uncertainties of less than 0.2%.