Phase Plane — Question 6

PDF ↗

Question 6

For r2=x2+y2r^2=x^2+y^2, consider x′=x(1−r2)−y,y′=y(1−r2)+x.x'=x(1-r^2)-y,\qquad y'=y(1-r^2)+x. For r>0r>0, use polar coordinates with a continuously followed angle along each trajectory. A periodic orbit is a nonconstant closed trajectory.

Tasks

  1. Find every equilibrium. Derive differential equations for rr and the angular coordinate.

  2. Identify the invariant unit circle, its orientation and least positive period. Determine whether nearby trajectories move toward it from each side.

  3. For any initial radius r0>0r_0>0, solve the equation for u=r2u=r^2 and prove convergence to the unit circle in forward time. State separately what happens at r0=0r_0=0.

  4. Does a nonzero solution converge to one particular point on the unit circle? Prove your answer and distinguish approaching an orbit from approaching an equilibrium.

Original worksheet page 1: question and worked solution for 5-6-006
Show solutionHide solution

Question 6 – Solution

Strategy. Separate the radial attraction from the angular motion; approaching a circle does not require stopping on it.

Step 1: Find the radial and angular equations. For a nonzero state, the equilibrium coefficient matrix has determinant (1−r2)2+1>0(1-r^2)^2+1>0, so it cannot annihilate that state. The only equilibrium is the origin. Direct differentiation gives r′=r(1−r2),θ′=1(r>0).\boxed{r'=r(1-r^2),\qquad\theta'=1\quad(r>0).} The formulas follow from rr′=xx′+yy′rr'=xx'+yy' and r2θ′=xy′−yx′r^2\theta'=xy'-yx'.

Step 2: Identify the periodic orbit and radial directions. At r=1r=1, the radius is constant and the angle increases at unit rate. This circle is a counterclockwise periodic orbit of least period 2π2\pi. For 0<r<10<r<1, r′>0r'>0; for r>1r>1, r′<0r'<0. Both sides move radially toward the circle rather than away from it.

Step 3: Prove the forward radial limit. The squared radius satisfies u′=2u(1−u)u'=2u(1-u), with solution u(t)=11+(r0−2−1)e−2t(t≥0).\boxed{u(t)=\frac 1{1+(r_0^{-2}-1)e^{-2t}}\quad(t\ge 0).} The denominator is positive for every t≥0t\ge 0, and direct differentiation and evaluation at zero verify the formula. It tends to 11, hence r→1r\to 1. With θ=θ0+t\theta=\theta_0+t, this also supplies a forward-global solution. At r0=0r_0=0 the solution remains at the equilibrium; the reciprocal formula is not used there.

Step 4: Rule out convergence to a single point. Along times t=2πnt=2\pi n the states tend to (cos⁡θ0,sin⁡θ0)(\cos\theta_0,\sin\theta_0). Along t=2πn+πt=2\pi n+\pi they tend to its negative. These are distinct, so no nonzero solution has a single limiting point. Its distance to the unit circle, |r−1||r-1|, tends to zero while it continues rotating. The figure shows the unit orbit and sample inward/outward spirals; all arrows represent increasing time.

See the diagram in the original worksheet below.

Original worksheet page 2: question and worked solution for 5-6-006

Original worksheet layout. Use Enlarge or open the PDF for a closer view.