Linear instability in planar viscoelastic Taylor-Couette flow with and without explicit polymer diffusion

Elastic turbulence has been found in computations of planar viscoelastic Taylor-Couette flow using the Oldroyd-B model, apparently generated by a linear instability (van Buel et al. Europhys. Lett., 124, 14001, 2018). We demonstrate that no such linear instability exists in the governing equations used unless some diffusion is added to the polymer conformation tensor equation, as might occur through a diffusive numerical scheme. With this addition, the polymer diffusive instability (PDI) (Beneitez et al. (Phys. Rev. Fluids, 8, L101901, 2023)) exists and leads to chaotic flows resembling those found by van Buel et al. (2018). We show how finite volume or finite-difference discretisations of the governing equations can naturally introduce diffusive errors near boundaries which are sufficient to trigger PDI. This suggests that PDI could well be important in numerical solutions of wall-bounded viscoelastic flows modelled using Oldroyd-B and FENE-P even with no polymer stress diffusion explicitly included.