


\documentclass[11pt,letterpaper]{article}

%\usepackage{natbib}
%\usepackage{cite}
\usepackage{biblatex}
\usepackage{graphicx}
\usepackage[margin=1.in,centering]{geometry}
\usepackage{hyperref}
\usepackage{caption}
\usepackage[export]{adjustbox}
\usepackage{float}


\bibliography{/home/caes/wmu/phy-4660/adv_lab.bib}


\begin{document}


\title{Lab 2: Chaos}
\author{Otho Ulrich, Eugene Kopf}
\maketitle

\begin{abstract}
The chaotic behaviour of a driven pendulum is explored. Phase space behaviours of a random response, non-chaotic periodic response, and chaotic response are generated from computational models. The functional form of each is compared, demonstrating the strange attractor as a distinguishing feature of a chaotic response. The motion of a physical pendulum is observed under damped and damped-driven conditions. The motion is characterized according to its phase space output, and determined to be chaotic.
\end{abstract}


\section{Chaos}
\label{sec:chaos}
Chaos is observed in many non-linear physical systems. It is the condition that a system's outcome is strongly sensitive to initial conditions. The changing conditions as the system evolves affect the outcome such that predicting the future state becomes impossible. The motion of a driven oscillator, such as a driven pendulum, for example, becomes unpredictable as the driving frequency and natural frequency of the pendulum interact. Damping can constrain the motion, and we find that while the motion is unpredictable, it still displays certain characteristics that can be analyzed.

For a finite period of time, chaotic behaviour isn't completely discernable from periodic behaviour, because the possibility exists that the function may repeat itself at some future time. To identify chaos, one makes a judgment after enough time has elapsed to assume for practical purposes the function will not repeat. Distinctions are seen between a random variable, periodic variable, and chaotic variable in phase space and poincar\'{e} sections. \cite{TANGLEDTALE}


    \subsection{Model of a Driven Pendulum}
    \label{sec:modeling}
    The driven pendulum exemplifies chaotic motion. In the angular coordinate $\theta$, the equation of motion of a driven pendulum is 

    \begin{center}
        $\frac{d^2}{dt^2} = \frac{{\omega_0}^2}{I} sin(\theta) - \frac{{\alpha}}{I} \frac{d\theta}{dt} + \frac{f}{I} cos(\omega t + \phi)$.
    \end{center}

    Here, $\omega_0$ represents the natural frequency of the pedulum, also its resonant frequency. The system will respond most strongly to the driver at this frequency. $\alpha$ is a damping term -- this can take a variety of forms, and in the experiment of section \ref{sec:experiment} is produced by a neodymium magnet interacting with the metal wheel of the pendulum. $f$ is the forcing amplitude where $\omega$ is the forcing frequency, offset from the angular coordinate by a phase $\phi$. $I$ is the moment of inertia of the pendulum.


\section{Phase Space}
\label{sec:phasespace}
A 2-dimensional phase space is a useful environment in which to identify the chaos in the motion in a coordinate. For an oscillator, the convenient coordinate is its angular position $\theta$, with the angular velocity response $\theta^\prime$. In figure \ref{fig:model_periodic}, observe how a periodic variable can be identified in phase space.

\begin{figure}
    \includegraphics[width=6.5in]{model_periodic.png}
    \caption{Computer-generated model of a periodic function in phase space. A pendulum's angular coordinate would correspond to $\theta=x$. The motion is predictable and orbits a single point.}
    \label{fig:model_periodic}
\end{figure}

    \subsection{Chaotic Attractor}
    Figure \ref{fig:model_damped_phase} demonstrates a damped oscillator, which exhibits stable critical points where the velocity goes to zero. Under forced conditions, the velocity does not converge to zero, but the motion produces orbits about these critical points in phase space, and we call these points chaotic attractors. Poincar\'{e} first postulated that chaos would be exemplified by complicated paths that roughly follow one of these orbits (the apex) about these attractors. \cite{CHAOSAT50} Attractors are a primary identifying characteristic of chaos, and should be observable in the chaotic motion of a forced pendulum. \cite{CHAOSDYNAMICS} A driven oscillator's path orbits around these critical points but can be seen to jump between them in an unpredictable way along the position coordinate; observe figure \ref{fig:model_driven_phase}.

    Damping still plays an important role in this chaotic motion; figure \ref{fig:model_driven_nodamp_phase} shows a driven pendulum with no damping, where the motion freely exhibits the effects of both the natural and forcing frequency, and becomes extremely complicated. The natural (un-driven) response of the pendulum is seen in the tall cirular strokes, which represent the pendulum's weight attempting to bring the pendulum to equilibrium. The driving frequency produces small-amplitude variations when the pendulum has a high angular speed and when the pendulum has a low angular speed, it can easily reverse the motion of the pendulum. The oscillator jumps around to many possible states, making it difficult to discern attractors. Chaotic attractors are much more clear when forcing interacts with damping in a system; figure \ref{fig:model_driven_phase} shows this very well.

    Computer models were generated using the Chaos for Java program written by Brian Davies. \cite{CHAOSFORJAVA} The path is computed from the circle to the triangle. 

    \begin{figure}
        \includegraphics[width=6.5in]{model_damped_phase.png}
        \caption{Model of a damped oscillator in phase space. The critical points are stable, since the velocity approaches zero from all points within the associated region. \cite{CHAOSDYNAMICS}}
        \label{fig:model_damped_phase}
    \end{figure}

    \begin{figure}
        \includegraphics[width=6.5in]{model_driven_phase.png}
        \caption{Computer-generated model of a damped and driven pendulum's angular motion. The motion jumps between critical points (called attractors in this context), and exhibits a semi-stable low-energy state around one of attractors. Arguments: $\theta=x, f=1, I=1, \omega=2/3, \alpha=0.7, k=\frac{\omega_0^2}{I}=2, \phi=0, \theta(0)=1, \theta^\prime(0)=1$,  150 time steps.}
        \label{fig:model_driven_phase}
    \end{figure}

    \begin{figure}
        \includegraphics[width=6.5in]{no_drag.png}
        \caption{Computer-generated model of a driven pendulum with no damping. The path through phase space is unstable and could easily take off toward either extreme of the $\theta$ coordinate. This sensitivity is characteristic of chaotic motion. Arguments: $\theta=x, f=1, I=1, \omega=0.5, \alpha=0, k=1.5, \phi=0$, 180 time steps.}
        \label{fig:model_driven_nodamp_phase}
    \end{figure}

    \subsection{Poincar\'{e} Map}
    A Poincar\'{e} section is a 2D phase space cross-section; in the case of the driven pendulum, the cross-section at a phase $\phi$ from the forcing term. One draws a map between each successive point to create a Poincar\'e map, which is a useful representation of a system's behaviour in phase space. A poincar\'e section of the driven pendulum model is shown in figure \ref{fig:model_driven_poincare}. An undriven and undamped oscillator will always return to the same point after one period, so a Poincar\'e map sampled using the period corresponding to the oscillator's natural frequency will consist of a single dot. The existence of multiple points at the same phase indicates chaotic motion.

    \begin{figure}
        \includegraphics[width=6.5in]{model_driven_poincare.png}
        \caption{Poincar\'e section of the damped and driven pendulum from figure \ref{fig:model_driven_phase} at $\phi=0$. Several dots are observed at this sampling phase, so the motion can be deemed chaotic.}
        \label{fig:model_driven_poincare}
    \end{figure}


\section{Experimental Driven Pendulum}
\label{sec:experiment}
An experimental driven pendulum is built in order to determine whether its motion is chaotic. The apparatus is shown in figure \ref{fig:apparatus}. The springs are connected so that they both hold about equal tension while the pendulum is at its top-most position, and so that they can stretch over at least one full turn of the pendulum in either direction. The tension was not measured precisely, nor the weight of the pendulum or the damping of the pendulum. The amplitude of the driving arm was $1.45" \pm 0.1"$. A photogate is used to determine the driving period when necessary.

\begin{figure}
    \includegraphics[width=6.5in]{appartus.png}
    \caption{Experimental driven pendulum. This angular motion of the round metal pendulum at the top is expected to exhibit chaotic motion. \cite{ADVLABCHAOS}}
    \label{fig:apparatus}
\end{figure}

    \subsection{Resonant Frequency}
    To determine the natural frequency of the pendulum, it was released from the top position under damped conditions and its frequency of oscillation was measured. The waveform is shown in figure \ref{fig:exp_damped_time}. Fourier transformations of the waveform in figure \ref{fig:exp_damped_fourier} reveal the natural frequency as the frequency where the power spectral density peaks, discounting the low-frequency peak from the flattened waveform due to damping. The natural frequency is estimated to be $0.93\pm0.6 Hz$. 
    
    \begin{figure}
        \includegraphics[width=6.5in]{exp_damped_time.png}
        \caption{The damped oscillations of the experimental pendulum. The fourier transforms of this function reveal the natural frequency of the pendulum, in figure \ref{fig:exp_damped_fourier}.}
        \label{fig:exp_damped_time}
    \end{figure}


    \begin{figure}
        \includegraphics[width=6.5in]{exp_damped_fourier.png}
        \caption{The fourier frequency-space representation of the damped oscillating. The peak near 0.8 Hz is the natural/resonant frequency of this pendulum.}
        \label{fig:exp_damped_fourier}
    \end{figure}




    \subsection{Periodic Motion}
    To produce periodic motion, the driving arm is run with a driving frequency $1.45 \pm 0.02s$. The observed motions are plotted in phase space in figure \ref{fig:exp_periodic_phase}. The motion still suggests some chaotic behaviour by the variation in the path taken about the critical point, but the poincar\'e section in figure \ref{fig:exp_periodic_poincare} shows that the orbit is converging.

    \begin{figure}
        \includegraphics[width=6.5in]{exp_periodic_phase.png}
        \caption{The periodic motion at driving period $1.45 s$is easily recognizable in phase space. It follows a stable orbit around the attractor. This plot still suggests chaotic behaviour, as the }
        \label{fig:exp_periodic_phase}
    \end{figure}

    \begin{figure}
        \includegraphics[width=6.5in]{exp_periodic_poincare.png}
        \caption{The periodic motion at driving period $1.45 s$is easily recognizable in phase space. It follows a stable orbit around the attractor. This plot still suggests chaotic behaviour, as the }
        \label{fig:exp_periodic_poincare}
    \end{figure}


\begin{figure}
    \includegraphics[width=6.5in]{chaotic_b_time.png}
    \caption{Driven oscillating with chaotic motion.}
    \label{fig:chaotic_b_time}
\end{figure}



\section{Discussion}
\label{sec:discussion}



\printbibliography


\end{document}
