This week's sets of classical pen and paper and computational exercises are a continuation of the topics from the previous homework set and lectures from the last two weeks. We keep dealing with simple motion problems and conservation laws; energy, momentum and angular momentum. These conservation laws are central in Physics and understanding them properly lays the foundation for understanding and analyzing more complicated physics problems.
The relevant reading background is
In both textbooks there are many nice worked out examples. Malthe-Sørenssen's text contains also several coding examples you may find useful.
The numerical homework focuses on another motion problem where you can use the code you developed in homework sets 3 and 4, almost entirely. Please take a look at the posted solutions (jupyter-notebook) for homework sets 3 and 4. You need only to change the forces at play, the one code to rule all exercises. The problem at hand is a classic, the gravitational force acting between the Sun and the Earth. Here you will notice also that the standard Euler-integration algorithm is not the best choice and we will introduce the so-called Euler-Cromer method and the Velocity-Verlet method. These methods will give much more stable numerical results with only few additions to your code.
The code you develop here will also be reused when we analyze energy conservation in homework set 6. And for those of you doing the honors project, it serves as a starting point for the solar system variant.
This exercise was partly discussed during the lectures of the weeks of January 30-February 3 and February 6-10.
We will study a classical electron which moves in the \( x \)-direction along a surface. The force from the surface is
$$ \boldsymbol{F}(x)=-F_0\sin{(\frac{2\pi x}{b})}\boldsymbol{e}_x. $$The constant \( b \) represents the distance between atoms at the surface of the material, \( F_0 \) is a constant and \( x \) is the position of the electron.
Taylor exercise 3.11. This exercise is meant to illustrate momentum conservation or not. This exercise was partly discussed during the lectures, see the notes from the week of February 6-10 on Energy and Momentum etc, see the part on Momentum conservation. Taylor's chapter 3.2 covers also this example.
Taylor exercises 3.13 (5pt) and 3.14 (5pt). This is a continuation of the previous exercise and most of the relevant background material can be found in Taylor chapter 3.2.
Taylor exercise 3.20. Here Taylor's chapter 3.3 can be of use. This relation will turn out to be very useful when we discuss systems of many classical particles.
The aim of this exercise (as well as the next) is to study the motion of objects under the influence of the gravitational force. We will limit ourselves to the Earth-Sun system. Here we will scale the equations and sketch our first algorithm for solving the equations, namely using Euler's method again, as we did in homework sets 3 and 4. This part together with the numerical part forms also the entry point for the solar system honors project. Furthermore, we will reuse parts of these results when analyzing energy conservation in homework 6.
We will limit ourselves (in order to test the algorithm) to a hypothetical solar system with the Earth only orbiting around the sun. The only force in the problem is gravity. Newton's law of gravitation is given by a force \( F_G \)
$$ F_G=\frac{GM_{\odot}M_{\mathrm{Earth}}}{r^2}, $$where \( M_{\odot} \) is the mass of the Sun and \( M_{\mathrm{Earth}} \) is the mass of the Earth. The gravitational constant is \( G \) and \( r \) is the distance between the Earth and the Sun. The Sun has a mass which is much larger than that of the Earth. We can therefore safely neglect the motion of the Sun in this problem.
We assume that the orbit of the Earth around the Sun is co-planar, and we take this to be the \( xy \)-plane. Using Newton's second law of motion we get the following equations
$$ \frac{d^2x}{dt^2}=\frac{F_{G,x}}{M_{\mathrm{Earth}}}, $$and
$$ \frac{d^2y}{dt^2}=\frac{F_{G,y}}{M_{\mathrm{Earth}}}, $$where \( F_{G,x} \) and \( F_{G,y} \) are the \( x \) and \( y \) components of the gravitational force.
We will use so-called astronomical units when rewriting our equations. Using astronomical units (AU as abbreviation)it means that one astronomical unit of length, known as 1 AU, is the average distance between the Sun and Earth, that is \( 1 \) AU = \( 1.5\times 10^{11} \) m. It can also be convenient to use years instead of seconds since years match better the time evolution of the solar system. The mass of the Sun is \( M_{\mathrm{sun}}=M_{\odot}=2\times 10^{30} \) kg. The mass of Earth is \( M_{\mathrm{Earth}}=6\times 10^{24} \) kg with a distance 1AU.
In setting up the equations we limit ourselves to a co-planar motion and use only the \( x \) and \( y \) coordinates. But you should feel free to extend your equations to three dimensions, it is not very difficult and the data from NASA are all in three dimensions.
NASA has an excellent site at http://ssd.jpl.nasa.gov/horizons.cgi#top. From there you can extract initial conditions in order to start your differential equation solver. At the above website you need to change from OBSERVER to VECTOR and then write in the planet you are interested in. The generated data contain the \( x \), \( y \) and \( z \) values as well as their corresponding velocities. The velocities are in units of AU per day. Alternatively they can be obtained in terms of km and km/s.
For the system below involving only the Earth and the Sun, you could just initialize the position with say \( x=1 \) AU and \( y=0 \) AU.
We assume that mass units can be obtained by using the fact that Earth's orbit is almost circular around the Sun.
For circular motion we know that the force must obey the following relation
$$ F_G= \frac{M_{\mathrm{Earth}}v^2}{r}=\frac{GM_{\odot}M_{\mathrm{Earth}}}{r^2}, $$where \( v \) is the velocity of Earth. The latter equation can be used to show that
$$ v^2r=GM_{\odot}=4\pi^2\mathrm{AU}^3/\mathrm{yr}^2. $$This exercise should be handed in as a jupyter-notebook at D2L. Remember to write your name(s).
In homework sets 3 and 4 we:
This week we will reuse our code from homework 4 (exercises 6) and replace the air resistance and force from the ground with the gravitational force. Then we will study the stability of the system as function of initial conditions and the time step length.
We start by importing some standard packages
# let's start by importing useful packages we are familiar with
import numpy as np
from math import *
import matplotlib.pyplot as plt
%matplotlib inline
Here we add a code-example which may aid in the above studies using Euler's forward method.
DeltaT = 0.001
#set up arrays
tfinal = 10 # in years
n = ceil(tfinal/DeltaT)
# set up arrays for time t, velocity v, and position r
t = np.zeros(n)
v = np.zeros((n,2))
r = np.zeros((n,2))
# Initial conditions as compact 2-dimensional arrays
r0 = np.array([1.0,0.0])
v0 = np.array([0.0,2*pi])
r[0] = r0
v[0] = v0
Fourpi2 = 4*pi*pi
# Start integrating using Euler's method
for i in range(n-1):
# Set up the acceleration
# Here you could have defined your own function for this
rabs = sqrt(sum(r[i]*r[i]))
a = -Fourpi2*r[i]/(rabs**3)
# update velocity, time and position using Euler's forward method
v[i+1] = v[i] + DeltaT*a
r[i+1] = r[i] + DeltaT*v[i]
t[i+1] = t[i] + DeltaT
The following gives you an opportunity to earn five extra credit points on each of the remaining homeworks and ten extra credit points on the midterms and finals. This assignment also covers an aspect of the scientific process that is not taught in most undergraduate programs: scientific writing. Writing scientific reports is how scientist communicate their results to the rest of the field. Knowing how to assemble a well written scientific report will greatly benefit you in you upper level classes, in graduate school, and in the work place.
The full information on extra credits is found at https://github.com/mhjensen/Physics321/blob/master/doc/Homeworks/ExtraCredits/. There you will also find examples on how to write a scientific article. Below you can also find a description on how to gain extra credits by attending scientific talks.
This assignment allows you to gain extra credit points by practicing your scientific writing. For each of the remaining homeworks you can submit the specified section of a scientific report (written about the numerical aspect of the homework) for five extra credit points on the assignment. For the two midterms and the final, submitting a full scientific report covering the numerical analysis problem will be worth ten extra points. For credit the grader must be able to tell that you put effort into the assignment (i.e. well written, well formatted, etc.). If you are unfamiliar with writing scientific reports, see the information here
The following table explains what aspect of a scientific report is due with which homework. You can submit the assignment in any format you like, in the same document as your homework, or in a different one. Remember to cite any external references you use and include a reference list. There are no length requirements, but make sure what you turn in is complete and through. If you have any questions, please contact Julie Butler at butler@frib.msu.edu.
| HW/Project | Due Date | Extra Credit Assignment |
|---|---|---|
| HW 3 | 2-8 | Abstract |
| HW 4 | 2-15 | Introduction |
| HW 5 | 2-22 | Methods |
| HW 6 | 3-1 | Results and Discussion |
| Midterm 1 | 3-12 | Full Written Report |
| HW 7 | 3-22 | Abstract |
| HW 8 | 3-29 | Introduction |
| HW 9 | 4-5 | Results and Discussion |
| Midterm 2 | _4-16_ | Full Written Report |
| HW 10 | 4-26 | Abstract |
| Final | 4-30 | Full Written Report |
You can also gain extra credits if you attend scientific talks. This is described here.
This opportunity will allow you to earn up to 5 extra credit points on a Homework per week. These points can push you above 100% or help make up for missed exercises. In order to earn all points you must:
Approved talks: Talks given by researchers through the following clubs:
If you have any questions please consult Jeremy Rebenstock, rebensto@msu.edu.
All the material on extra credits is at https://github.com/mhjensen/Physics321/blob/master/doc/Homeworks/ExtraCredits/.