We outline a comprehensive numerical procedure for modeling of species transport and surface reaction kinetics in electrokinetically actuated microfluidic devices of rectangular cross section. Our results confirm the findings of previous simplified approaches that a concentration wave is created for sufficiently long microreactors. An analytical solution, developed for the wave propagation speed, shows that, when normalizing with the fluid mean velocity, it becomes a function of three parameters comprising the channel aspect ratio, the relative adsorption capacity, and the kinetic equilibrium constant. Our studies also reveal that the reactor geometry idealized as a slit, instead of a rectangular shape, gives rise to the underestimation of the saturation time. The extent of this underestimation increases by increasing the Damkohler number or decreasing the dimensionless Debye–Hückel parameter. Moreover, increasing the values of the Damkohler number, the dimensionless Debye–Hückel parameter, the relative adsorption capacity, and the velocity scale ratio results in lower saturation times.
