Numerical study on the chaotic stirring of viscous flow in an alternately driven cavity has been performed. Even under the Stokes-flow assumption, the inherent singularity at the corners made the problem not so easily accessible. With some special treatments to the region near the corners, the biharmonic equation was solved numerically by using the fully implicit method. The velocity field was then used in obtaining the trajectories of passive particles for studying the stirring effect. The three tools developed in the field of the nonlinear dynamics and chaos, that are the Poincare sections, the unstable manifolds, and the Lyapunov exponents, were used in analysing the stirring effect. It was shown that the unstable manifolds obtained in this study well fit the experimental results given by the previous investigators. It is predicted that the best stirring can be obtained when the aspect ratio a is near 0.8 and the dimensionless period T is in the range 4.3 - 4.7.