Numerical characterisation of folding flow microchannel mixers

Micromixers have been considered in numerous recent studies with the aim of mixing different liquid streams for the common circumstance of non-inertial flow, i.e., in the Stokes flow regime. Under such conditions, the diffusion of momentum is dominant but the diffusion of species remains weak because the Schmidt number of liquids is large. Most mixers that have potential for application in the Stokes regime make use of a folding flow pattern that approximates the baker's transformation. In the work presented here, the general scaling of mixers of this type is developed from the exact equation for species transport and computations are made for a specimen mixer geometry to test the effectiveness of the resulting scaling. The scaling relation developed is found to give an excellent representation of the actual mixing characteristics of the specimen mixer over the entire range of Péclet number of practical interest. Finite volume computations are employed to solve the governing equations up to around Pe=10³. At higher Péclet numbers, where finite volume numerical solution becomes inaccurate with affordable mesh sizes, the species equation is solved using a Monte Carlo method instead. Finally, the scaling relation is used to develop the design relations needed to determine the number of mixing elements, the pressure drop incurred and the Péclet number of operation to achieve a given mixture uniformity within a specified mixing time.