← Classical Mechanics
🌀

Nonlinear Dynamics & Chaos

Most real mechanical systems are nonlinear, and nonlinearity can produce chaos — deterministic yet unpredictable motion with exponential sensitivity to initial conditions. Understanding chaos requires new tools: Lyapunov exponents, Poincaré sections, and bifurcation diagrams.

Key Concepts

Sensitive Dependence
Two nearby trajectories in a chaotic system diverge exponentially: δx(t)δx0eλt|\delta x(t)|\approx|\delta x_0|e^{\lambda t}, where λ>0\lambda>0 is the Lyapunov exponent. This makes long-term prediction impossible even with excellent initial knowledge.
Fixed Points & Stability
A fixed point satisfies x˙=0\dot{\vec x}=0. It is stable (attractor) if nearby trajectories converge, unstable (repeller) if they diverge. Stability is determined by the eigenvalues of the Jacobian f/x\partial\vec f/\partial\vec x at the fixed point.
Period Doubling
Many systems transition to chaos through a cascade of period-doubling bifurcations. The ratio of successive bifurcation parameters converges to the Feigenbaum constant δ4.669\delta\approx4.669.
Poincaré Section
A 2D slice through the phase space taken stroboscopically (once per forcing period for driven systems). Periodic orbits appear as points; chaotic orbits fill a fractal strange attractor.

Key Equations

Lyapunov exponent
λ=limt1tlnδx(t)δx0\lambda = \lim_{t\to\infty}\frac{1}{t}\ln\frac{|\delta x(t)|}{|\delta x_0|}
Positive λ indicates chaos; largest Lyapunov exponent characterizes the system.
Logistic map
xn+1=rxn(1xn)x_{n+1} = rx_n(1-x_n)
Simplest chaotic map; chaos appears for r > 3.57... Feigenbaum constant δ ≈ 4.669.
Pendulum (nonlinear)
θ¨+glsinθ=βθ˙+Acos(ωdt)\ddot{\theta} + \frac{g}{l}\sin\theta = -\beta\dot{\theta} + A\cos(\omega_d t)
Driven damped pendulum: exhibits chaos for large driving amplitude A.
Worked Example

Stability of a Fixed Point

Problem

The Lorenz system has a fixed point at the origin. The Jacobian there has eigenvalues λ1=22.8\lambda_1=-22.8, λ2=11.8\lambda_2=11.8, λ3=2.67\lambda_3=-2.67 (for standard parameters σ=10\sigma=10, ρ=28\rho=28, β=8/3\beta=8/3). Is the origin stable?

Solution

A fixed point is stable only if all eigenvalues have negative real parts.

λ2=11.8>0\lambda_2=11.8>0: at least one eigenvalue is positive.

Therefore the origin is an unstable saddle point — trajectories near it will be repelled in the direction of the positive eigenvalue.

Answer Unstable — the positive eigenvalue λ2=11.8\lambda_2=11.8 causes nearby trajectories to diverge from the origin.
Practice

Exercises

7 problems
1 of 7

Two chaotic trajectories start δx0=1010|\delta x_0|=10^{-10} m apart with Lyapunov exponent λ=2.0\lambda=2.0 /s. After t=10t=10 s, find δx|\delta x| (in m). Round to one decimal.

m
2 of 7

In the logistic map xn+1=rxn(1xn)x_{n+1}=rx_n(1-x_n), find the non-zero fixed point xx^* for r=3.0r=3.0. (x=11/rx^*=1-1/r)

(dimensionless)
3 of 7

The first period doubling of the logistic map occurs at r1=3.0r_1=3.0, the second at r2=3.449r_2=3.449, the third at r3=3.544r_3=3.544. Estimate the Feigenbaum constant δ=(r2r1)/(r3r2)\delta=(r_2-r_1)/(r_3-r_2).

Unlock Exercise 3

Subscribe to PhysWeb Pro to access all exercises and track your progress.

Upgrade to Pro →
4 of 7

A pendulum has l=1.0l=1.0 m, g=9.8g=9.8 m/s². Find the small-angle frequency ω0\omega_0 (in rad/s).

Unlock Exercise 4

Subscribe to PhysWeb Pro to access all exercises and track your progress.

Upgrade to Pro →
5 of 7

For a 1D map with Jacobian eigenvalue μ=0.5|\mu|=0.5 at a fixed point, how many iterations nn does it take for a perturbation to shrink by a factor of 100? Use μn=0.01|\mu|^n=0.01.

Unlock Exercise 5

Subscribe to PhysWeb Pro to access all exercises and track your progress.

Upgrade to Pro →
6 of 7

The Hénon map has a strange attractor. Its Lyapunov exponent is λ0.42\lambda\approx0.42 bits/iteration. How many iterations for an uncertainty to grow by a factor of e10e^{10}?

Unlock Exercise 6

Subscribe to PhysWeb Pro to access all exercises and track your progress.

Upgrade to Pro →
7 of 7

A system has N=3N=3 positive Lyapunov exponents with values 2.1, 0.8, and 0.3 /s. What is the Kaplan-Yorke dimension contribution from these? (Sum of positive exponents / largest.) This estimates the information dimension growth rate.

Unlock Exercise 7

Subscribe to PhysWeb Pro to access all exercises and track your progress.

Upgrade to Pro →

Key Takeaways

  • Positive Lyapunov exponents signal chaos — exponential divergence of nearby trajectories.
  • Period doubling is the most common route to chaos; the Feigenbaum constant δ4.669\delta\approx4.669 is universal.
  • Poincaré sections reveal the underlying structure: periodic orbits → points, chaos → strange attractors.
  • Chaos is deterministic but practically unpredictable beyond the Lyapunov time tL=1/λt_L=1/\lambda.