  The ability to generate complete, or almost complete, chaotic mixing is of
great interest in numerous applications, particularly for microfluidics. For
this purpose, we propose a strategy that allows us to quickly target the
parameter values at which complete mixing occurs. The technique is applied to a
time periodic, two-dimensional electro-osmotic flow with spatially and
temporally varying Helmoltz-Smoluchowski slip boundary conditions. The strategy
consists of following the linear stability of some key periodic pathlines in
parameter space (i.e., amplitude and frequency of the forcing), particularly
through the bifurcation points at which such pathlines become unstable.
