03_01_oscillations.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Oscillatory Motion\n",
    "==================\n",
    "\n",
    "The harmonic oscillator\n",
    "-----------------------\n",
    "\n",
    "We know that simple oscillatory motion is described by a general\n",
    "quadratic potential of the form $$V=\\frac{1}{2}kx^2,$$ where $k$ is a\n",
    "“force constant”. In a spring, this is a measure of the stiffness of the\n",
    "spring. In a pendulum, it is given by $$k=\\frac{mg}{L}$$ where $m$ is\n",
    "the mass of the particle and $L$ its length. From this potential we can\n",
    "extract the force $$F=-\\frac{dV}{dx}=-kx$$ and the equation of motion\n",
    "can be written as \n",
    "$$\\frac{d^2x}{dt^2}=-\\omega _0 ^2x,\n",
    "$$ \n",
    "where the quantity $\\omega _0$ is defined by\n",
    "$$\\omega _0^2=\\frac{k}{m}$$ The dynamical motion described by these\n",
    "equations is called “simple harmonic motion” and can be solved\n",
    "analytically: $$x(t)=A\\cos{(\\omega _0t+\\delta)},\n",
    "$$ \n",
    "where $a$ is the amplitude of th oscillations and\n",
    "$\\delta$ is a phase that depends on the initial conditions.\n",
    "\n",
    "We can now define a period $T$ such as $$x(t+T)=x(t).$$ Since\n",
    "$\\omega _0T$ corresponds to one cycle, we obtain\n",
    "$$T=\\frac{2\\pi}{\\omega _0}=\\frac{2\\pi}{\\sqrt{k/m}}.$$ The frequency\n",
    "$\\nu$ is the number of cycles per unit of time and is given by\n",
    "$\\nu=1/T$. Note that the period of the motion depends on the ratio\n",
    "$k/m$ and is completely independent of the initial conditions.\n",
    "\n",
    "### Exercise 3.1: Energy conservation \n",
    "\n",
    "1.  In Exercise 1.6 we wrote program to simulate simple harmonic motion.\n",
    "    Now use the program to calculate the relative change in the total\n",
    "    energy during one cycle $\\Delta _n=(E_n-E_0)/E_0$. Is the function\n",
    "    $\\Delta$ uniformly small during the cycle? Choose $x_0=1$, $v_0=0$,\n",
    "    and $\\omega_0^2=9$ Compare using Euler’s method, 2nd.\n",
    "    order Runge-Kutta, and Verlet's algorithm.\n",
    "\n",
    "2.  Plot position and velocity as a function of the time $t$.\n",
    "\n",
    "3.  Compute the amplitude $A$ for the initial conditions $x_0=4$,\n",
    "    $v_0=0$, and $x_0=0$, $v_0=4$; choose $\\omega _0^2=4$ in both cases.\n",
    "    What quantity determines the value of $A$?.\n",
    "\n",
    "4.  Compute the average value of the kinetic energy and the potential\n",
    "    energy during a complete cycle. Is there a relation between the two\n",
    "    averages?\n",
    "\n",
    "5.  Plot the path of the oscillator in phase space $(x,v)$. Set\n",
    "    $\\omega _0^2=9$ and use different initial conditions\n",
    "    $(x_0,v_0)=(1,0);(0,1);(4,0)$. Do you find different paths for each\n",
    "    of them? What physical quantity distinguishes or characterizes each\n",
    "    path? What is the shape of the phase paths? Is the motion of a\n",
    "    representative point $(x,v)$ always in the clockwise or\n",
    "    counterclockwise direction?\n",
    "\n",
    "### Challenge 3:1\n",
    "\n",
    "Modify the code for Exercise 1.6 to solve Exercise 3.1\n",
    "\n",
    "### Exercise 3.2: Large oscillations \n",
    "\n",
    "The pendulum responds to the equation of motion for simple harmonic motion only in\n",
    "the limit of small angles. In the case of large oscillations, the\n",
    "equation has to be modified becoming\n",
    "$$\\frac{d^2\\theta}{dt^2}=-\\frac{g}{L}\\sin{\\theta}.$$\n",
    "\n",
    "The energy of the pendulum is then given by\n",
    "$$E=\\frac{1}{2}mL^2(\\frac{d\\theta}{dt})^2+mgL(1-\\cos{\\theta})$$\n",
    "\n",
    "1.  Modify your program to simulate large amplitude oscillations in\n",
    "    a pendulum. Set $g/L=9$ and choose $\\Delta t$ so that the numerical\n",
    "    solution is stable, <span>*i.e.*</span> it does no diverge with time\n",
    "    from the “true” solution. Check the stability by calculating the\n",
    "    total energy and ensuring that it does not drift from its\n",
    "    initial value.\n",
    "\n",
    "2.  Set $d\\theta /dt|_{t=0}=0$ and make plots of $\\theta (t)$ and\n",
    "    $d\\theta /dt(t)$ for the initial conditions $\\theta(t=0)=0.1$, 0.2,\n",
    "    0.4, 0.8, 1.0. Describe the qualitative behavior of $\\theta$ and\n",
    "    $d\\theta /dt$. What is the period $T$ and the maximum amplitude\n",
    "    $\\theta _{max}$ in each case? Plot $T$ versus $\\theta _{max}$ and\n",
    "    discuss the qualitative dependence of the period on the amplitude.\n",
    "    How do the results compare in the linear and non-linear cases,\n",
    "    <span>*e.g.*</span> which period is larger? Explain the relative\n",
    "    values of $T$ in physical terms.\n",
    "    \n",
    "### Challenge 3.2:\n",
    "\n",
    "Modify the code and solve Exercise 3.2\n",
    "\n",
    "### Exercise 3.3: Damped oscillations \n",
    "\n",
    "If a drag force is included in the problem, the equation of motion\n",
    "becomes: $$\\frac{d^2x}{dt^2}=-\\omega _0^2x-\\gamma \\frac{dx}{dt},$$ where\n",
    "the “damping coefficient” gamma is a measure of the friction. Note that\n",
    "the drag force opposes the motion.\n",
    "\n",
    "1.  Incorporate the effects of damping in your program and plot the time\n",
    "    dependece of position and velocity. Make runs for $\\omega _0^2=9$,\n",
    "    $x_0=1$, $v_0=0$, and $\\gamma=0.5$.\n",
    "\n",
    "2.  Compare the period and angular frequency to the undamped case. Is the\n",
    "    period longer or shorter? Make additional runs for $\\gamma = \n",
    "    1,2,3$. Does the frequency increase or decrease with greater\n",
    "    damping?\n",
    "\n",
    "3.  Define the amplitude as the maximum value of $x$ in one cycle.\n",
    "    Compute the “relaxation time” $\\tau$, the time it takes for the\n",
    "    amplitude to change from its maximum to $1/e\\approx 0.37$ of its\n",
    "    maximum value. To do this, plot the maximum amplitud of each cycle,\n",
    "    and fit it with an exponential of the form $A_0 \\exp{-t/\\tau}$.\n",
    "    Compute $\\tau$ for each of the values of $\\gamma$ used in the previous\n",
    "    item and discuss the qualitative dependence of $\\tau$ with $\\gamma$.\n",
    "\n",
    "4.  Plot the total energy as a function of time for the values of\n",
    "    $\\gamma$ considered previously. If the decrease of the energy is not\n",
    "    monotonic, explain the cause of the time-dependence.\n",
    "\n",
    "5.  Compute the average value of the kinetic energy, potential energy,\n",
    "    and total energy over a complete cycle. Plot these averages as a\n",
    "    function of the number of cycles. Due to the presence of damping,\n",
    "    these averages decrease with time. Is the decrease uniform?\n",
    "\n",
    "6.  Compute the time-dependence of $x(t)$ and $v(t)$ for\n",
    "    $\\gamma=4,5,6,7,8$. Is the motion oscillatory for all $\\gamma$?\n",
    "    Consider a condition for equilibrium $x<0.0001$; how quickly does\n",
    "    $x(t)$ decay to equilibrium? For fixed $\\omega _0$, the oscillator\n",
    "    is said to be “critically damped” at the smallest values of $\\gamma$\n",
    "    for which the decay to equilibrium is monotonic. For what value of\n",
    "    $\\gamma$ does critical damping occur for $\\omega_0^2=9$ and\n",
    "    $\\omega^2 _0=4$? For each value of $\\omega _0$ compute the value of\n",
    "    $\\gamma$ for which the system approaches equilibrium more quickly.\n",
    "\n",
    "7.  Construct the phase space diagram for cases $\\omega_0 ^2=9$ and\n",
    "    $\\gamma=0.5$, 2, 4, 6, 8. Area the qualitative features of the paths\n",
    "    independent of $\\gamma$? If not, discuss the\n",
    "    qualitative differences.\n",
    "\n",
    "### Challenge 3.3:\n",
    "\n",
    "Modify the code and solve Exercise 3.3\n",
    "\n",
    "### Exercise 3.3: Linear response to external forces \n",
    "\n",
    "How can we determine the period of a pendulum that is not already in\n",
    "motion? The obvious way is to disturb the system, for instance, applying a\n",
    "small displacement and observe the resulting motion. We will finr that\n",
    "the “response” of a system is actually an intrinsic property of the\n",
    "system and can tell us about its nature in the absence of perturbations.\n",
    "\n",
    "Consider a damped linear oscillator with an external force $F(t)$\n",
    "$$\\frac{d^2x}{dt^2}=-\\omega _0^2x-\\gamma\\frac{dx}{dt}+\\frac{1}{m}F(t).$$\n",
    "It is customary to interpret the reponse of the system in terms of the\n",
    "displacement $x$. The time dependence in $F(t)$ is arbitrary. A\n",
    "particular case is when the force is harmonic:\n",
    "$$\\frac{1}{m}F(t)=A_0 \\cos{\\omega t},\n",
    "$$ where $\\omega$ is the angular frequency of the\n",
    "driving force.\n",
    "\n",
    "1.  Modify your program so that an external force of the form\n",
    "    (\\[driving force\\]) is included. Set $\\omega _0^2=9$, $\\gamma=0.5$,\n",
    "    $A_0=1$ and $\\omega = 2$ (we’ll use these values for the rest of\n",
    "    the exercise). These values correspond to a lightly\n",
    "    damped oscillator. Plot $x(t)$ versus $t$ for the initial conditions\n",
    "    $(x_0=1,v_0=0)$. How does the qualitative behavior differ from the\n",
    "    unperturbed case? What is the period and angular frequency of $x(t)$\n",
    "    after several oscillations? Obtain a similar plot with\n",
    "    $(x_0=0,v_0=1)$. What is the period and angular frequency after\n",
    "    several oscillations? Does $x(t)$ approach a limiting behavior that\n",
    "    is independent of the initial conditions? Identify a “transient”\n",
    "    part of $x(t)$ which depends on the initial conditions and decays in\n",
    "    time, and a “steady state” part which dominates at longer times and\n",
    "    which is independent of the initial conditions.\n",
    "\n",
    "2.  Compute $x(t)$ for $\\omega = 1$ and $\\omega = 4$. What is the period\n",
    "    and angular frequency of the steady state in each case?\n",
    "\n",
    "3.  Compute $x(t)$ for $\\omega _0=4$. What is the angular frequency of\n",
    "    the steady state motion? On the basis of these results, explain which\n",
    "    parameters determine the frequency of the steady state behavior.\n",
    "\n",
    "4.  Verify that the steady state behavior is given by\n",
    "    $$x(t)=A(\\omega )\\cos{\\omega t+\\delta},$$ where $\\delta$ is the\n",
    "    phase difference between the applied force and the steady\n",
    "    state motion. Compute $\\delta$ for\n",
    "    $\\omega_0^2=9,\\gamma=0.5$,$\\omega=0$, $1.0$, $2.0$, $2.2$, $2.4$,\n",
    "    $2.6$, $2.8$, $3.0$, $3.2$, $3.4$. Repeat the computation for\n",
    "    $\\gamma=1.5$ and plot $\\delta$ versus $\\omega$ for the two values of\n",
    "    $\\gamma$. Discuss the qualitative dependence of $\\delta(\\omega )$ in\n",
    "    the two cases.\n",
    "\n",
    "### Exercise 3.5: Resonance \n",
    "\n",
    "The long term behavior of the driven harmonic oscillator depends on the\n",
    "frequency of the driving force. One measure of this behavior is the\n",
    "maximum of the steady state displacement $A(\\omega )$.\n",
    "\n",
    "1.  Adopt the initial condition $(x_0=0,v_0=0)$. Compute $A(\\omega )$\n",
    "    for $\\omega =0$, $1.0$, $2.0$, $2.2$, $2.4$, $2.6$, $2.8$, $3.0$,\n",
    "    $3.2$, $3.4$ with $\\omega _0=3$ and $\\gamma =0.5$. Plot $A(\\omega )$\n",
    "    versus $\\omega$ and describe its qualitative behavior. If\n",
    "    $A(\\omega )$ has a maximum, determine the “resonance angular\n",
    "    frequency” $\\omega _{max}$, which is the frequency at the maximum of\n",
    "    $A$. Is the value of $\\omega _{max}$ close to the natural angular\n",
    "    frequency $\\omega _0$?\n",
    "\n",
    "2.  Compute $A_{max}$, the value of the amplitude at $\\omega _{max}$,\n",
    "    and the ratio $\\Delta \\omega /\\omega _{max}$, where $\\Delta \\omega$\n",
    "    is the “width” of the resonance. Define $\\Delta \\omega$ as the\n",
    "    frequency interval between points on the amplitude curve which are\n",
    "    $1/\\sqrt{2}A_{max}$. Set $\\omega _0=3$ and consider $\\gamma=0.1$,\n",
    "    0.5, 1.0, 2.0. Describe the qualitative dependence of $A_{max}$ and\n",
    "    $\\Delta \\omega/\\omega _{max}$ on $\\gamma$. The quantity\n",
    "    $\\Delta \\omega /\\omega_{max}$ is proportional to $1/Q$, where $Q$ is\n",
    "    the “quality factor” of the oscillator.\n",
    "\n",
    "3.  Decribe the qualitative behavior of the steady state amplitude\n",
    "    $A(\\omega )$ near $\\omega=0$ and $\\omega\\gg \\omega _0$. Why is\n",
    "    $A(\\omega =0) > A(\\omega )$ for small $\\omega$? Why does\n",
    "    $A(\\omega ) \n",
    "    \\rightarrow 0$ for $\\omega \\gg \\omega _0$?\n",
    "\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 2",
   "language": "python",
   "name": "python2"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 2
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython2",
   "version": "2.7.10"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}