We construct self-similar solutions to the 2D Navier–Stokes equations evolving from arbitrarily large $-1$–homogeneous initial data and present numerical evidence for their non-uniqueness.