Because of the form of the linear part of the vector field it seems advisable to explore the dynamic of the Euclidean radius. For simplicity of computation, use $E=\frac12r^2=\frac12(x^2+y^2)$ to get $$ \frac{d}{dt}E=x\dot x+y\dot y=-2E+2E(x^2+2y^2) $$ so that $$ 2E(2E-1)\le\dot E\le 2E(4E-1) $$ This means that $\dot E$ is negative for $0< E< \frac14$ and positive for $E>\frac12$. Looking in the direction of negative times, the time reversed ODE, this means that the annulus $\frac14
![diagram showing limit cycle]( some numerical solutions, the circles have radii $0.6$ to $1.1$ by $0.1$.