An exact solution of the lubrication equations for the Oldroyd-B model in a hyperbolic channel

An exact similarity solution of the lubrication equations for the steady flow of a viscoelastic Oldroyd-B fluid in a contracting and symmetric hyperbolic channel is derived. The solution is valid for small values of the Deborah number, De (the ratio of the polymer's longest relaxation time to a characteristic residence time of the fluid in the channel), all values of the polymer viscosity ratio, η (the ratio of polymer viscosity to the total viscosity of the fluid), as well as to typical values of the contraction ratio, Λ (the ratio of the channel height at the inlet to the channel height at the outlet). The solution, the significance of which for the hyperbolic geometry is analogous to the classic Poiseuille solution in a straight channel, precisely satisfies the exact analytical solution of the Oldroyd-B model along the walls, but does not satisfy any initial conditions for the polymer extra-stresses at the inlet of the hyperbolic section of the channel. The exact solution is given in terms of the streamfunction which is used in the momentum balance to derive a non-linear ordinary differential equation with an unknown function which corresponds to a modified fluid velocity. The final equation is solved numerically using a fully spectral method with a Galerkin-type approach to calculate the unknown function and the pressure gradient. The exact solution for the polymer extra-stresses allows for deriving a variety of expressions for the average pressure-drop in the channel. The most accurate expression is that resulting from the mechanical energy of the flow. Its range of validity is determined in terms of De, Λ and η, and confirmed by numerical pseudospectral simulations of the lubrication equations. In all cases, a decrease in the pressure drop compared to the Newtonian value with increasing De and/or η is predicted.