Numerical study of electroosmotic slip flow of fractional Oldroyd-B fluids at high zeta potentials

In this paper, an investigation of the electroosmotic flow of fractional Oldroyd-B fluids in a narrow circular tube with high zeta potential is presented. The Navier linear slip law at the walls is considered. The potential field is applied along the walls described by the nonlinear Poisson–Boltzmann equation. It's worth noting here that the linear Debye–Hückel approximation can't be used at the condition of high zeta potential and the exact solution of potential in cylindrical coordinates can't be obtained. Therefore, the Matlab bvp4c solver method and the finite difference method are employed to numerically solve the nonlinear Poisson–Boltzmann equation and the governing equations of the velocity distribution, respectively. To verify the validity of our numerical approach, a comparison has been made with the previous work in the case of low zeta potential and the excellent agreement between the solutions is clear. Then, in view of the obtained numerical solution for the velocity distribution, the numerical solutions of the flow rate and the shear stress are derived. Furthermore, based on numerical analysis, the influence of pertinent parameters on the potential distribution and the generation of flow is presented graphically.