Phase Plane — Question 7

PDF ↗

Question 7

Consider the damped nonlinear system x′=y,y′=−x−x3−cy,c>0,x'=y,\qquad y'=-x-x^3-cy,\qquad c>0, with initial state (1,0)(1,0). Define E(x,y)=12y2+12x2+14x4E(x,y)=\tfrac 12y^2+\tfrac 12x^2+\tfrac 14x^4. A level curve of EE need not itself be a solution curve.

Tasks

  1. Find the equilibrium and derive the exact change of EE along a solution.

  2. Use the initial energy to prove bounds on both coordinates for every forward time. Explain why the smooth system then has a global forward solution.

  3. At the initial state, E′=0E\prime=0. Does this mean that the solution is stationary or remains on its original energy level? Use the original equations to decide.

  4. Prove that no nonconstant periodic orbit exists, by integrating the energy balance over a hypothetical period. Explain what the enclosing energy curve represents in the phase picture.

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

Question 7 – Solution

Strategy. A decreasing energy supplies a containing region and rules out a closed return without requiring an explicit trajectory.

Step 1: Compute the equilibrium and dissipation. Equilibrium requires y=0y=0 and x+x3=0x+x^3=0, hence only (0,0)(0,0). Along solutions, E′=y(−x−x3−cy)+(x+x3)y=−cy2≤0.E'=y(-x-x^3-cy)+(x+x^3)y=\boxed{-cy^2\le 0}. Thus energy decreases except at instants when y=0y=0.

Step 2: Bound the forward solution. Initially E=3/4E=3/4. Since all energy terms are nonnegative, y2/2≤3/4y^2/2\le 3/4 gives |y|≤3/2|y|\le\sqrt{3/2}. Also x4+2x2≤3x^4+2x^2\le 3, so (x2+1)2≤4(x^2+1)^2\le 4 and |x|≤1|x|\le 1. Therefore |x(t)|≤1,|y(t)|≤3/2(t≥0).\boxed{|x(t)|\le 1,\qquad |y(t)|\le\sqrt{3/2}\quad(t\ge 0).} The polynomial field is locally Lipschitz. These bounds keep the state in a compact region, preventing finite-time escape and allowing continuation for all forward times.

Step 3: Interpret zero instantaneous energy loss. At (1,0)(1,0), the velocity is (0,−2)(0,-2), not zero. Moreover y(t)=−2t+o(t)y(t)=-2t+o(t) for small positive tt, so yy becomes nonzero and E′=−cy2<0E'=-cy^2<0. The solution enters lower energy levels. A horizontal zero of E′E' at one instant does not make energy constant on an interval or make its level curve a trajectory.

Step 4: Exclude nonconstant periodic motion. If a solution had period T>0T>0, integrating gives 0=E(T)−E(0)=−c∫0Ty(t)2dt0=E(T)-E(0)=-c\int_0^T y(t)^2\,dt. Continuity and c>0c>0 force y≡0y\equiv 0, then x′=0x'=0 and x+x3=0x+x^3=0, so the solution is the zero equilibrium. Thus no nonconstant periodic orbit exists. The drawn E=3/4E=3/4 curve bounds a forward-invariant region; it is not a closed orbit. Interior arrows show velocities for the illustrative choice c=1c=1.

See the diagram in the original worksheet below.

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

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