University of California, Davis Physics Department
High School Apprenticeship Program in Computational Physics: Summer 2021

Sponsored by:

The United States Department of Energy
"The mission of the Energy Department is to ensure America's security and prosperity by addressing its energy, environmental and nuclear challenges through transformative science and technology solutions. "

Oak Ridge Institute of Science Education
"Supporting DOE and other federal agencies' missions to strengthen the nation's science education and research initiatives"
(Special thanks to Jennifer Tyrell of ORAU.)

Our FANTASTIC 2021 Summer Interns!

Left to Right: Ivan Jimenez, Bryce Lu, Erwin Jimenez, Aiden Rijov, Jalan Masongsong, Diego Jimenez,
Mario Arzola, Naveen Bonoan, Jenessa Vasquez, Cara Mann, Aracely Magana, Summer 2031 recruit! Anisha Vikash

Our students come from River City High School in West Sacramento.

Our instructors

Richard Scalettar (and friends) Valentin Taufour Dubarrie Fagout Andrea Lopez-Arguello Chunhan Feng Yuxi Zhang Maxwell Wright

What we do- Part 1:

We learn to code in python! This program computes the trajectory of our rockets. It is the final version which includes
the variation of mass as the rocket engine burns and also the effects of air resistance!!

# THIS IS NEEDED TO USE ARRAYS IN PYTHON:
from array import *
# THIS IS NEEDED TO MAKE PLOTS IN PYTHON:
import matplotlib.pyplot as plt
# INPUTS
dt=float(input('Enter dt: '))
m0 =float(input('Enter starting mass: '))
y =0.001
ymax=0.
v =float(input('Enter v0: '))
g = -9.8
THR =float(input('Enter thrust: '))
tTHR =float(input('Enter duration of thrust: '))
pm =float(input('Enter propellant mass : '))
b =float(input('Enter air resistance: '))
# START OFF TIME, POSITION, VELOCITY ARRAYS
time =array('f', [0.])
ysave=array('f',[y])
vsave=array('f',[v])
# DO THE MD AND APPEND EACH NEW TIME,
POSITION, AND VELOCITY TO ITS ARRAY
t=0.0
while y > 0:
t=t+dt
y=y+v*dt
if (y > ymax):
ymax=y
if (t < tTHR):
m=m0-pm*t/tTHR
F=THR+m*g-b*v
v=v+(F/m)*dt
else:
m=m0-pm
F=m*g-b*v
v=v+(F/m)*dt
if (abs(t-tTHR)<=dt/2):
print('height when engines turn off',y)
time.append(t)
ysave.append(y)
vsave.append(v)
print('time of flight',t)
print('apogee ',ymax)
plt.plot(time, ysave)
plt.xlabel('t')
plt.ylabel('y')
plt.show()

We build Estes rockets and we test our code against real data!!! The python code works really well!!!

Athena Rocket Data comparison
with two engine types AStrocam and Phantom Blue
Rocket Data comparison Patriot Rocket Data comparison

Athena Rocket Data comparison
Now including air resistance The Teams and
all their data.

What we do- Part 2:

We build and fly rockets. This involves learning constuction techniques, basic principles of rocket engines, etc.
We then compare to codes which compute trajectories.

Rocket Construction Rocket Construction TA's Yuxi and Andrea

Rocket construction Doing some trajectory calculations At Russell Field

Ignition is imminent Enjoying the shade Two of the teams