Viscoelastic fluids are highly challenging from the rheological standpoint, and their discretization demands robust, efficient numerical solvers. Simulating viscoelastic flows requires combining Navier–Stokes system with a dynamic tensorial equation, increasing mathematical computational demands. Hence, fractional-step methods decoupling calculation of flow quantities an attractive option. In c...