The Schmidt number \(Sc\) , appearing solely in the polymer diffusion term and defined as the ratio between the solvent kinematic viscosity and polymer diffusivity [1]},
is typically of order \(O (10^6)\)  in physical applications.
In this work, enhanced diffusion (i.e. lower \(Sc\) ) had to be employed to regularize the hyperbolic equation (), as is customarily done in other works involving viscoelastic direct numerical simulations [2]}, [1]}.
This is also necessitated by the spectral methods embedded in our branch continuation scheme.
which behave slightly worse than finite difference methods in this respect,
pushing the maximum admissible Schmidt number down to \(Sc = 250\)  from typically 1000
[2]}.
