In this investigation, the transient electroosmotic flow of multi-layer immiscible viscoelastic fluids in a slit microchannel is studied. Through an appropriate combination of the momentum equation with the rheological model for Maxwell fluids, an hyperbolic partial differential equation is obtained and semi-analytically solved by using the Laplace transform method to describe the velocity field. In the solution process, different electrostatic conditions and electro-viscous stresses have to be considered in the liquid-liquid interfaces due to the transported fluids content buffer solutions based on symmetrical electrolytes. By adopting a dimensionless mathematical model for the governing and constitutive equations, certain dimensionless parameters that control the start-up of electroosmotic flow appear, as the viscosity ratios, dielectric permittivity ratios, the density ratios, the relaxation times, the electrokinetic parameters and the potential differences. In the results, it is shown that the velocity exhibits an oscillatory behavior in the transient regime as a consequence of the competition between the viscous and elastic forces; also, the flow field is affected by the electrostatic conditions at the liquid-liquid interfaces, producing steep velocity gradients, and finally, the time to reach the steady-state is strongly dependent on the relaxation times, viscosity ratios and the number of fluid layers.

This work investigates the transient multilayer electro-osmotic flow of viscoelastic fluids through an annular microchannel. The dimensionless mathematical model of multilayer flow is integrated by the linearized Poisson-Boltzmann equation, the Cauchy momentum equation, the rheological Maxwell model, initial conditions, and the electrostatic and hydrodynamic boundary conditions at liquid-liquid and solid-liquid interfaces. Although the main force that drives the movement of fluids is due to electrokinetic effects, a pressure gradient can also be added to the flow. The semi-analytical solution for the electric potential distribution and velocity profiles considers analytical techniques as the Laplace transform method, with numerical procedures using the inverse matrix method for linear algebraic equations and the concentrated matrix exponential method for the inversion of the Laplace transform. The results presented for velocity profiles and velocity tracking at the transient regime reveal an interesting oscillatory behavior that depends on elastic fluid properties via relaxation times. The time required for the flow to reach steady-state is highly dependent on the viscosity ratios and the dimensionless relaxation times. In addition, the influence of other dimensionless parameters on the flow as the electrokinetic parameters, zeta potentials at the walls, permittivity ratios, ratio of pressure forces to electro-osmotic forces, number of fluid layers, and annular thickness are investigated. The findings of this study have significant implications for the precise control of parallel fluid transport in microfluidic devices for flow-focusing applications.

The present investigation analyzes the transient multilayer electro-osmotic flow through an annular microchannel with hydrophobic walls. The fluids are considered immiscible and viscoelastic, following the Maxwell rheological model. In the problem examined, the linearized Poisson–Boltzmann and Cauchy momentum equations are used to determine the electric potential distribution and the flow field, respectively. Here, different interfacial phenomena are studied through the imposed boundary conditions, such as the hydrodynamic slip and specified zeta potentials at solid–liquid interfaces, the velocity continuity, the electroviscous stresses balance, the potential difference, and the continuity of electrical displacements at the interfaces between fluids. The semi-analytic solution uses the Laplace transform theory. In the results, the velocity profiles and velocity tracking show the oscillatory behavior of flow, which strongly depends on the dimensionless relaxation time. Furthermore, the hydrodynamic slip on the channel walls contributes to the release of energy stored in the fluids due to elastic effects at the start-up of the flow. Similarly, other dimensionless parameters are also investigated. This research aims to predict the parallel flow behavior in microfluidic devices under electro-osmotic effects.

The present work analyzes the combined viscoelectric and steric effects on electroosmotic flow in a soft channel with polyelectrolyte coating. The structured channel surface, which controls the electric potential, creates two different flow regions: the electrolyte flow within the permeable polyelectrolyte layer (PEL) and the bulk electrolyte. Thus, this study discusses the interaction of various electrostatic effects to predict the electroosmotic flow field. The nonlinear governing equations describing the fluid flow are the modified Poisson–Boltzmann equation for the electric potential distribution, the mass conservation equation, and the modified Navier–Stokes equations for the flow field, which are solved numerically using a one-dimensional (1D) scheme. The results indicate that the flow enhances when increasing the electric potential magnitude across the channel cross-section via the rise in different dimensionless parameters, such as the PEL thickness, the steric factor, and the ratio of the electrokinetic parameter of the PEL to that of the electrolyte layer. This research demonstrates that the PEL significantly enhances control over electroosmotic flow. However, it is crucial to consider that viscoelectric effects at high electric fields and the friction generated by the grafted polymer brushes of the PEL can reduce these benefits.

This paper presents an analytical study of the transient electroosmotic flow for Newtonian fluids through a parallel flat plate microchannel with heterogeneous zeta potentials. The dimensionless mathematical model is based on the Poisson-Boltzmann, mass and momentum conservation governing equations together with the lubrication theory. The distribution of the zeta potentials at the walls obeys to a sinusoidal function, which includes dimensional parameters as Δζ that controls the magnitude and polarity of the zeta potentials, being capable to produce slanted velocity profiles and inverse flows. On the other hand, the combination of the phase angle between the sinusoidal functions of the zeta potentials ω, the dimensionless parameter of their amplitude Δζ, and the parameter that controls the frequency of the sinusoidal functions m, induce additional perturbations on the flow, which is directly related to the dimensionless pressure distribution and to the transient flow field. The transient behavior characteristics of the electroosmotic flow are discussed in terms of the zeta potential variations. It is demonstrated that the results for the transient electroosmotic flow, predict the influence of the main dimensionless parameters above mentioned on the velocity profiles and the streamlines. This work about the perturbations on the electroosmotic flow by heterogeneous zeta potentials, contributes to a better understanding of the transport phenomena in microfluidic devices for future mixing applications.