We analyze the pressure-driven flow of a viscoelastic fluid in arbitrarily shaped, narrow channels and present a theoretical framework for calculating the relationship between the flow rate $q$ and pressure drop $Delta p$. We utilize the Oldroyd-B mo
del and first identify the characteristic scales and dimensionless parameters governing the flow in the lubrication limit. Employing a perturbation expansion in powers of the Deborah number ($De$), we provide analytical expressions for the velocity, stress, and the $q-Delta p$ relation in the weakly viscoelastic limit up to $O(De^2)$. Furthermore, we exploit the reciprocal theorem derived by Boyko $&$ Stone (Phys. Rev. Fluids, vol. 6, 2021, pp. L081301) to obtain the $q-Delta p$ relation at the next order, $O(De^3)$, using only the velocity and stress fields at the previous orders. We validate our analytical results with two-dimensional numerical simulations of the Oldroyd-B fluid in a hyperbolic, symmetric contracting channel and find excellent agreement. For the flow-rate-controlled situation, both our theory and simulations reveal weak dependence of the velocity field on the Deborah number, so that the velocity can be approximated as Newtonian. In contrast to the velocity, the pressure drop strongly depends on the viscoelastic effects and decreases with $De$. Elucidating the relative importance of different terms in the momentum equation contributing to the pressure drop, we identify that a pressure drop reduction for narrow contracting geometries is primarily due to gradients in the viscoelastic shear stresses, while viscoelastic axial stresses have a minor effect on the pressure drop along the symmetry line.
We present an analytical study, validated by numerical simulations, of electroosmotic flow in a Hele-Shaw cell with non-uniform surface charge patterning. Applying the lubrication approximation and assuming thin electric double layer, we obtain a pai
r of uncoupled Poisson equations which relate the pressure and the stream function, respectively, to gradients in the zeta potential distribution parallel and perpendicular to the applied electric field. We solve the governing equations for the fundamental case of a disk with uniform zeta potential and show that the flow-field in the outer region takes the form of a pure dipole. We illustrate the ability to generate complex flow-fields around smooth convex regions by superposition of such disks with uniform zeta potential and a uniform pressure driven flow. This method may be useful for future on-chip devices, allowing flow control without the need for mechanical components.
