Title: Nonlinear Mechanics and Chaos
1Nonlinear Mechanics and Chaos
12
- Dr. Jen-Hao Yeh
- Prof. Anlage
2This is a brief introduction to the ideas and
concepts of nonlinear mechanics, and a
discussion of various quantitative methods for
analyzing such problems We will focus on the
driven damped pendulum (DDP)
3Linear Nonlinear
easy hard
special general
analytical numerical
Superposition principle
Chaos
4Nonlinear
Chaotic
Linear
5Driven, Damped Pendulum (DDP)
We expect something interesting to happen as g ?
1, i.e. the driving force becomes comparable to
the weight
6A Route to Chaos
7NDSolve in Mathematica
8Driven, Damped Pendulum (DDP)
For all following plots
Period 2p/w 1
9Small Oscillations of the Driven, Damped Pendulum
g ltlt 1 will give small oscillations
Linear
g 0.2
f(t)
t
After the initial transient dies out, the
solution looks like
Periodic attractor
10Small Oscillations of the Driven, Damped Pendulum
g ltlt 1 will give small oscillations
- The motion approaches a unique periodic attractor
- independent of initial conditions
- The motion is sinusoidal with the same frequency
as the drive
11Moderate Oscillations of the Driven, Damped
Pendulum
g lt 1 and the nonlinearity becomes significant
Try
This solution gives from the f3 term
Since there is no cos(3wt) on the RHS, it must be
that all develop a cos(3wt) time dependence.
Hence we expect
We expect to see a third harmonic as the driving
force grows
12Harmonics
13Moderate Driving The Nonlinearity Distorts the
cos(wt)
cos(wt)
f(t)
f(t)
g 0.9
t
t
-cos(3wt)
The motion is periodic, but
The third harmonic distorts the simple
14Even Stronger Driving Complicated Transients
then Periodic!
After a wild initial transient, the motion
becomes periodic
f(t)
g 1.06
t
After careful analysis of the long-term motion,
it is found to be periodic with the same period
as the driving force
15Slightly Stronger Driving Period Doubling
After a wilder initial transient, the motion
becomes periodic, but period 2!
f(t)
g 1.073
t
The long-term motion is TWICE the period of the
driving force!
A SUB-Harmonic has appeared
16Harmonics and Subharmonics
17Slightly Stronger Driving Period 3
The period-2 behavior still has a strong period-1
component Increase the driving force slightly and
we have a very strong period-3 component
g 1.077
f(t)
Period 3
t
18Multiple Attractors
The linear oscillator has a single attractor for
a given set of initial conditions
For the drive damped pendulum Different
initial conditions result in different long-term
behavior (attractors)
g 1.077
f(t)
Period 3
t
Period 2
19Period Doubling Cascade
f(t)
f(t)
Early-time motion
g 1.06
t
Close-up of steady-state motion
Period 1
t
g 1.078
Period 2
g 1.081
Period 4
g 1.0826
Period 8
20Period Doubling Cascade
g 1.06
g 1.078
g 1.081
g 1.0826
21Bifurcation Points in the Period Doubling
Cascade
Driven Damped Pendulum
n period gn interval
(gn1-gn)
1 1 ? 2 1.0663 0.0130 2
2 ? 4 1.0793 0.0028 3
4 ? 8 1.0821 0.0006 4
8 ? 16 1.0827
The spacing between consecutive bifurcation
points grows smaller at a steady rate
? as n ? 8
d 4.6692016 is called the Feigenbaum number
The limiting value as n ? 8 is gc 1.0829.
Beyond that is chaos!
22Period Doubling Cascade
Period doubling continues in a sequence of
ever-closer values of g
Such period-doubling cascades are seen in many
nonlinear systems Their form is essentially the
same in all systems it is universal
23Period infinity
24Chaos!
g 1.105
f(t)
t
The pendulum is trying to oscillate at the
driving frequency, but the motion remains erratic
for all time
25Chaos
- Nonperiodic
- Sensitivity to initial conditions
26Sensitivity of the Motion to Initial Conditions
Start the motion of two identical pendulums with
slightly different initial conditions Does their
motion converge to the same attractor? Does it
diverge quickly?
27Convergence of Trajectories in Linear Motion
Take the logarithm of Df(t) to magnify small
differences.
Plotting log10Df(t) vs. t should be a
straight line of slope b, plus some wiggles from
the lncos(w1t d1) term
Note that log10x log10e lnx
28Convergence of Trajectories in Linear Motion
Log10Df(t)
t
g 0.1
Df(0) 0.1 Radians
The trajectories converge quickly for the small
driving force ( linear) case
This shows that the linear oscillator is
essentially insensitive to its initial conditions!
29Convergence of Trajectories in Period-2 Motion
Log10Df(t)
t
g 1.07
Df(0) 0.1 Radians
The trajectories converge more slowly, but still
converge
30Divergence of Trajectories in Chaotic Motion
g 1.105
Log10Df(t)
Df(0) 0.0001 Radians
t
If the motion remains bounded, as it does in this
case, then Df can never exceed 2p. Hence this
plot will saturate
The trajectories diverge, even when very close
initially Df(16) p, so there is essentially
complete loss of correlation between the pendulums
Extreme Sensitivity to Initial Conditions Practica
lly impossible to predict the motion
31The Lyapunov Exponent
l Lyapunov exponent
l lt 0 periodic motion in the long term
l gt 0 chaotic motion
32Linear Nonlinear Nonlinear
Linear Chaos
Drive period Harmonics, Subharmonics, Period-doubling Nonperiodic, Extreme sensitivity
l lt 0 l lt 0 l gt 0
33What Happens if we Increase the Driving Force
Further? Does the chaos become more intense?
Df(0) 0.001 Radians
g 1.13
f(t)
Log10Df(t)
t
t
Period 3 motion re-appears!
With increasing g the motion alternates between
chaotic and periodic
34What Happens if we Increase the Driving Force
Further? Does the chaos re-appear?
Df(0) 0.001 Radians
g 1.503
f(t)
Log10Df(t)
t
t
Chaotic motion re-appears!
This is a kind of rolling chaotic motion
35Divergence of Two Nearby Initial Conditions for
Rolling Chaotic Motion
g 1.503
f(t)
t
Df(0) 0.001 Radians
Chaotic motion is always associated with extreme
sensitivity to initial conditions
Periodic and chaotic motion occur in narrow
intervals of g
36Bifurcation Diagram
37Period Doubling Cascade
Period doubling continues in a sequence of
ever-closer values of g
Such period-doubling cascades are seen in many
nonlinear systems Their form is essentially the
same in all systems it is universal
Sub-harmonic frequency spectrum
Driven Diode experiment
F0 cos(wt)
w/2
38Period Doubling Cascade
Period doubling continues in a sequence of
ever-closer values of g
Such period-doubling cascades are seen in many
nonlinear systems Their form is essentially the
same in all systems it is universal
A period doubling cascade in convection of
mercury in a small convection cell. The plots
show the temperature at one fixed point in the
cell as a function of time, for four successively
larger temperature gradients as given by the
parameter R/Rc
Fig. 12.9, Taylor
The Brain-behaviour Continuum The Subtle
Transition Between Sanity and Insanity By Jose
Luis. Perez Velazquez
39Bifurcation Diagram
Used to visualize the behavior as a function of
driving amplitude g
- Choose a value of g
- Solve for f(t), and plot a periodic sampling of
the function - t0 chosen at a time after the attractor behavior
has been achieved - Move on to the next value of g
1.0793
1.0663
Fig. 12.17
40Construction of the Bifurcation Diagram
f(t)
g 1.06
Period 1
t
g 1.078
Period 2
g 1.081
Period 4
Period 6 window
g 1.0826
Period 8
41The Rolling Motion Renders the Bifurcation
Diagram Useless
g 1.503
f(t)
t
Df(0) 0.001 Radians
As an alternative, plot
42Bifurcation Diagram Over a Broad Range of g
Period-1 followed by period doubling bifurcation
Previous diagram range
Mostly chaos
Mostly chaos
Mostly chaos
Period-3
Rolling Motion (next slide)
43Period-1 Rolling Motion at g 1.4
g 1.4
t
f(t)
t
Even though the pendulum is rolling, is
periodic
44An Alternative View State Space Trajectory
Plot vs. with time as a parameter
f(t)
t
g 0.6
periodic attractor
start
First 20 cycles
Cycles 5 -20
Fig. 12.20, 12.21
45An Alternative View State Space Trajectory
Plot vs. with time as a parameter
g 0.6
start
periodic attractor
Cycles 5 -20
First 20 cycles
The periodic attractor , is an
ellipse
The state space point moves clockwise on the orbit
Fig. 12.22
46State Space Trajectory for Period Doubling Cascade
g 1.078
g 1.081
Period-2
Period-4
Plotting cycles 20 to 60
Fig. 12.23
47State Space Trajectory for Chaos
g 1.105
Cycles 14 - 21
Cycles 14 - 94
The orbit has not repeated itself
48State Space Trajectory for Chaos
Chaotic rolling motion Mapped into the interval
p lt f lt p
Cycles 10 200
This plot is still quite messy. Theres got to
be a better way to visualize the motion
49The Poincaré Section
Similar to the bifurcation diagram, look at a
sub-set of the data
- Solve for f(t), and construct the state-space
orbit - Plot a periodic sampling of the orbit
- with t0 chosen after the attractor behavior has
been achieved
Samples 10 60,000
Enlarged on the next slide
50The Poincaré Section is a Fractal
The Poincaré section is a much more elegant way
to represent chaotic motion
51The Superconducting Josephson Junction as a
Driven Damped Pendulum
2
1
I
(Tunnel barrier)
The Josephson Equations
f f2-f1 phase difference of SC wave-function
across the junction
52Radio Frequency (RF) Superconducting Quantum
Interference Devices (SQUIDs)
Flux Quantization in the loop
Ic
F
I(t)
53Single rf-SQUID
D. Zhang, et al., Phys. Rev. X (in press),
arXiv1504.08301
54D. Zhang, et al., Phys. Rev. X (in press),
arXiv1504.08301
55THz Emission from the Intrinsic Josephson
Effect A classic problem in nonlinear physics
F0 h/2e 2.07 x 10-15 Tm2
DC voltage on junction creates an oscillating
f(t), which in turn creates an AC current
that radiates
L. Ozyuzer, et al., Science 318, 1291 (2007)
Best emission is seen when the crystal is
partially heated above Tc! Results are extremely
sensitive to details (number of layers, edge
properties, type of material, width of mesa,
etc.) Many competing states do not show
emission Emission enhanced near cavity mode
resonances ? requires non-uniform current
injection, assisted by inhom. heating, p-phase
kinks, crystal defects
56Chaos in Newtonian Billiards
Imagine a point-particle trapped in a 2D
enclosure and making elastic collisions with the
walls
Describe the successive wall-collisions with a
mapping function
- The Chaos arises due to the shape of the
boundaries enclosing the system.
Computer animation of extreme sensitivity to
initial conditions for the stadium billiard