SUMMARYA cell-vertex hybrid ÿnite volume=element method is investigated that is implemented on triangles and applied to the numerical solution of Oldroyd model uids in contraction ows. Particular attention is paid to establishing high-order accuracy, whilst retaining favourable stability properties. Elevated levels of elasticity are sought. The main impact of this study reveals that switching from quadratic to linear ÿnite volume stress representation with discontinuous stress gradients, and incorporating local reduced quadrature at the re-entrant corner, provide enhance stability properties. Solution smoothness is achieved by adopting the non-conservative ux form with area integration, by appealing to quadratic recovered velocity-gradients, and through consistency considerations in the treatment of the time term in the constitutive equation. In this manner, high-order accuracy is maintained, stability is ensured, and the ÿner features of the ow are conÿrmed via mesh reÿnement. Lip vortices are observed for We¿1, and a trailing-edge vortex is also apparent. Loss of evolution and solution asymptotic behaviour towards the re-entrant corner are also discussed.