Quantification of chaotic strength and mixing in a micro fluidic system

Abstract

Comparative studies of five different techniques commonly employed to identify the chaotic strength and mixing efficiency in micro fluidic systems are presented to demonstrate the competitive advantages and shortcomings of each method. The ‘chaotic electroosmotic stirrer’ of Qian and Bau (2002 Anal. Chem. 74 3616–25) is utilized as the benchmark case due to its well-defined flow kinematics. Lagrangian particle tracking methods are utilized to study particle dispersion in the conceptual device using spectral element and fourth-order Runge–Kutta discretizations in space and time, respectively. Stirring efficiency is predicted using the stirring index based on the box counting method, and Poincaré sections are utilized to identify the chaotic and regular regions under various actuation conditions. Finite time Lyapunov exponents are calculated to quantify the chaotic strength, while the probability density function of the stretching field is utilized as an alternative method to demonstrate the statistical analysis of chaotic and partially chaotic cases. Mixing index inverse, based on the standard deviation of scalar species distribution, is utilized as a metric to quantify the mixing efficiency. Series of numerical simulations are performed by varying the Peclet number (Pe) at fixed kinematic conditions. The mixing time (tm) is characterized as a function of the Pe number, and tm ∝ ln(Pe) scaling is demonstrated for fully chaotic cases, while tm ∝ Peα scaling with α ≈ 0.33 and α = 0.5 are observed for partially chaotic and regular cases, respectively. Employing the aforementioned techniques, optimum kinematic conditions and the actuation frequency of the stirrer that result in the highest mixing/stirring efficiency are identified.