Developed by Deva O'Neil  Published July 23, 2017
Subject Area  Mechanics 

Level  Beyond the First Year 
Available Implementations  Glowscript, IPython/Jupyter Notebook, and Mathematica 
Learning Objectives 
* Use numerical method (EulerCromer) to model oscillating systems and produce plots of the motion. (**Exercise 1**)
* Predict how motion responds to changes in initial conditions (**Exercise 1**)
* Apply concept of energy conservation to undamped and damped oscillators (**Exercise 2**)
* Make connection between the velocitydependence of damping force and the loss of total mechanical energy in the damped oscillator(**Exercise 2**)
* Interpret phase space diagrams (**Exercise 3**)

Time to Complete  100 min 
These exercises are not tied to a specific programming language. Example implementations are provided under the Code tab, but the Exercises can be implemented in whatever platform you wish to use (e.g., Excel, Python, MATLAB, etc.).
Test your predictions. If your results do not agree with your predictions, determine the source of the discrepancy before continuing. * Finally, what do you predict would happen if $b$ is set to zero? Record you prediction below, then test it.
* Explain the role that $b$ plays in this system:
#Exercise 2 The goal of this exercise is to understand how the energy of a damped harmonic oscillator changes with time. We will assume that our block starts out displaced to the right of its equilibrium position, with $v(0) = 0$. We will use the classical equation for kinetic energy, $K = \frac{1}{2}m~v^2$. The potential energy for a spring stretched a distance $\vec{r}$ from (0,0,0) is $U = \frac{1}{2} k~ r^2$. * Taking into account the initial conditions, sketch your predictions for $K(t)$, $U(t)$ for the case of $b = 0$ (no damping): ![](images/Damping/K_and_U_axes.png "") * With no damping, how do you expect the total mechanical energy to evolve with time?
Open the template for Exercise 2 and test your predictions. There are blanks for you to fill in the kinetic and potential energies. (Because this template is also used in Exercise 3, you'll see blocks of code commented out  leave those alone for now.) In testing your predictions to the graphs generated by the program, compare both the starting point of each graph and the overall behavior. If your results do not agree with your predictions, determine the source of the discrepancy before continuing. Test your prediction about total energy by changing the graph of U to a graph of total mechanical energy (K+U). * In Exercise 1, you should have found that giving the system an initial velocity increases the amplitude of the motion. Explain this result using energy reasoning.
**Damped Oscillator** We will now put a nonzero damping constant back in. Before running the code, record your predictions: * If there was damping, how would you expect the total mechanical energy to evolve with time? Make a sketch to support your answer: ![](images/Damping/total_mechanical_axes.png "") * How would you expect kinetic energy to evolve with time? Make a sketch to support your answer: ![](images/Damping/Kinetic_axes.png "") Test your predictions  your code should exhibit graphs of $K$ and $K + U$. * Explain why the decay of the total mechanical energy is not a smooth curve. (This is not a numerical artifact; it has a physical explanation.) *Hint:* At what times in the block's motion is the loss of mechanical energy the greatest? Explain your answer by revisiting the damping force equation.
#Exercise 3 To fully describe the state of a system at a point in time requires knowing both position and velocity (or, equivalently, momentum). One can describe the system at a point in the "space" of these two variables. For onedimensional motion, a graph of the motion in this "phase space" is a graph of $p$ vs. $x$. Both depend on time, so as time evolves, a path in the phase space is traced out. In the template, set $b = 0$. We'll start with the undamped oscillator. Sketch a prediction for the phase space graph for a damped oscillator ($ b = 0.1$) ![](images/Damping/phase.JPG "") In your program, remove the comments that disable the phase space graph. Don't forget to also uncomment the line in the while loop that adds to the phase space data. Compare the results to your prediction. Experiment with higher values of $b$. * Explain how to tell from the phase space plot what level of damping the block is subjected to.
* Explain how to tell from a phase space plot whether the system returns to its original state.
Download Options
Share a Variation
Credits and Licensing
Deva O'Neil, "Energy and Phase Space of a Damped Oscillator," Published in the PICUP Collection, July 2017, https://doi.org/10.1119/PICUP.Exercise.Damping.
DOI: 10.1119/PICUP.Exercise.Damping
The instructor materials are ©2017 Deva O'Neil.
The exercises are released under a Creative Commons AttributionNonCommercialShareAlike 4.0 license