This preview shows pages 1–2. Sign up to view the full content.
This preview has intentionally blurred sections. Sign up to view the full version.
View Full Document
Unformatted text preview: Ph507. Homework 9 (due: Wednesday, April 20).
PROBLEM 81 (15 pts.) Perform numerical study of the pendulum subjected to timedependent torque: H= p2 + 2
2 0 cos + cos t a) Investigate the transition from periodic to chaotic behavior by constructing a Poincare map of the system for several values of . (one way of doing this is to run the simulation several times with random initial conditions). b) Study the sensitivity of the dynamics to the initial conditions. In particular, consider simultaneous motion of two or more identical systems which are originally very close. Follow the evolution of their relative "distance" (in phase space). Compare results for chaotic and periodic orbits. In the chaotic case, estimate the Lyapunov exponents, as functions of . c) After running the simulation for certain time T , run the system "backward in time" over the same time interval T . Similarly to part (b) follow the evolution of the "distance" between two points in phase space. Perform this study both in periodic and chaotic regimes. In the latter case, investigate the dependence of the behavior on the choice of T . Interpret your results. Solution a) Typical Poincare map is below. The chaotic sea emerges near the separatrix, and expands with .
2.5 2 1.5 1 0.5 0 0.5 1 1.5 2 2.5 3 2 1 0 1 2 3 b) Coevolution:
10
0 =.1 dashed line: = /5
0 =5 10 2 distance 10 4 10 6 =0.01 10 8 10 10 10 12 periodic orbit
14 10 0 50 100 t
0 150 200 250 1 Typically, Lyapunov exponent increases for small and then saturates. At the saturation, 1 1/T were T is the characteristic period of the oscillations. 2 = 1 , and 3 = 0. c) Timereversal:
10
2 10 0 backward
10
2 distance 10 4 10 6 10 8 forward
10
10 10 12 0 10 20 30 40 50 60 70 80 90 100 t 0 While the periodic orbits are fully reversible, the chaotic ones become effectively irreversible after certain time. That is because the numerical error gets exponentially enhanced, as exp (1 t). 2 ...
View
Full
Document
 Spring '04
 ANON
 mechanics, Work

Click to edit the document details