Dispersion coefficient in an electro-osmotic flow of a viscoelastic fluid through a microchannel with a slowly varying wall zeta potential

Abstract

The dispersion coefficient of a passive solute in a steady-state pure electro-osmotic flow (EOF) of a viscoelastic liquid, whose rheological behaviour follows the simplified Phan-Thien–Tanner (sPTT) model, along a parallel flat plate microchannel, is studied. The walls of the microchannel are assumed to have modulated and low $\unicode[STIX]{x1D701}$ potentials, which vary slowly in the axial direction in a sinusoidal manner. The flow field required to obtain the dispersion coefficient was solved using the lubrication approximation theory (LAT). The solution of the electric potential is based on the Debye–Hückel approximation for a symmetric $(z:z)$ electrolyte. The viscoelasticity of the fluid is observed to notably amplify the axial distribution of the effective dispersion coefficients due to the variation in the $\unicode[STIX]{x1D701}$ potentials of the walls. The problem was formulated for two cases: when the Debye layer thickness (EDL) was on the order of unity (thick EDL) and in the limit where the thickness of the EDL was very small compared with the height of the microchannel (thin EDL limit). Due to the coupling between the nonlinear governing equations and the sPTT fluid model, they were replaced by their approximate linearized forms and solved in the limit of $\unicode[STIX]{x1D700}\ll 1$ using the regular perturbation technique. Here $\unicode[STIX]{x1D700}$ is the amplitude of the sinusoidal function of the $\unicode[STIX]{x1D701}$ potentials. Additionally, the numerical solution of the simplified governing equations was also obtained for $\unicode[STIX]{x1D700}=O(1)$ and compared with the approximate solution, showing excellent agreement for $0\leqslant \unicode[STIX]{x1D700}\leqslant 0.3$. Note that the dispersion coefficient primarily depends on the Deborah number, on the ratio of the half-height of the microchannel to the Debye length, and on the assumed variation in the $\unicode[STIX]{x1D701}$ potentials of the walls.