This paper reports the results of a theoretical modeling of the fully developed electroosmotic flow in a rectangular microchannel with a high-density polyelectrolyte layer (PEL) attached to the walls. At these conditions, the ions are partitioned between the PEL and the fluid outside the PEL owing to the difference between the permittivities of the two media. It is taken into account that the dynamic viscosity is higher within the PEL because of hydration effects. Solutions are obtained for the electric potential and velocity distributions as well as the mean velocity by making use of a variational approach, applied to the linearized form of the governing equations, that treats the whole area under consideration as a single domain with variable physical properties. The resulting equations are solved using a spectral method. Closed-form analytical expressions are obtained for a slit geometry, representing the case of high aspect ratios. The solutions obtained are validated by comparing to finite-element simulations of the full non-linear equations. It is shown that the electric potential drop inside the channel increases due to the depletion of the counterions within the PEL caused by the ion partitioning effect. This effect, surprisingly, magnifies the electroosmotic flow rate because of the increase of the space charge outside the PEL. As expected, the hydration effects reduce the flow rate, especially for thick PELs.
