Open source simulation tool for electrophoretic stacking, focusing, and separation

We present the development, formulation, and performance of a new simulation tool for electrophoretic preconcentration and separation processes such as capillary electrophoresis, isotachophoresis, and field amplified sample stacking. The code solves the one-dimensional transient advection-diffusion equations for multiple multivalent weak electrolytes (including ampholytes) and includes a model for pressure-driven flow and Taylor–Aris dispersion. The code uses a new approach for the discretization of the equations, consisting of a high resolution compact scheme which is combined with an adaptive grid algorithm. We show that this combination allows for accurate resolution of sharp concentration gradients at high electric fields, while at the same time significantly reducing the computational time. We demonstrate smooth, stable, and accurate solutions at current densities as high as 5000 A/m2 using only 300 grid points, and a 75-fold reduction in computational time compared with equivalent uniform grid techniques. The code is available as an open source for free at http://microfluidics.stanford.edu.