Project2.ipynb

Kernel: Python 2 (SageMath)

** Instructions: Write up your answers for parts a,b, and c and turn in a single copy for the entire group. For part c, include a description of what code you added. I'll collect your code electronically. Include all group members names on the sheet you turn in and let me know which group member has the correct code in their SageMathCloud account.**

Project 2: Part A: Escape from ice station McMurdo

---

You are a member of a scientific research team at McMurdo ice station which is funded by the Carver Media Group in Antarctica. Two members of your research team have recently returned from investigating an incident at a Norwegian research facility. They brought with them a burnt humanoid body with two faces. Since the disturbing discovery several inhabitants of the ice station have disappeared. Frightened, a member of your team decided to flee the station on a fan powered hovercraft but you receive a distress call not long after their escape that their steering and acceleration controls have been jammed and they need your help.

Time	Your Team Member's Position	Your Position
0 s	2536.40 m	10.47 m
10 s	3072.80 m	41.88 m
20 s	3609.20 m	94.22 m

You decide to attempt a rescue in another hovercraft. You must decide how many members of your team help in the rescue operation. The hovercrafts do not have a velocity or acceleration gauge but they do have GPS locators and you possess your trusty stop watch. The GPS locator tells you the exact position of both your craft and other team members craft relative to the ice station. You are following their path. You collect the following data for the first 20 seconds of your journey.

You need to tell the runaway researcher the exact time from your starting time to jump onto your hovercraft as you may only have one shot at this rescue.

• Example on page 72 of text is helpful here. Ask the professor for help if you don't have your text with you.

Project 2: Part B: Escape from ice station McMurdo

---

Just as you are about to radio the time to jump to the runaway researcher, you realize the steering and acceleration controls have become frozen on your hovercraft and so it continues to accelerate and you cannot change direction. 200m ahead of the point at which you were going to tell the researcher to jump is an ice ravine. At the bottom of the ice ravine, 400m below, is an unfrozen salt water pool surrounded by stalagmites. From the ravine's edge to the pool is 490m and the pool stretches for 900m. You are moving too quickly to survive jumping off the hovercraft, but might survive the fall into the pool by staying on the hovercraft; it has seat belts. You now have a choice to make, to stay on your hovercraft or jump to the runaway researcher's hovercraft. One or both may make it to the pool. Your choice may be the difference between life and death.

Project 2: Part C: Escape from ice station McMurdo

---

Surprisingly enough hovercrafts are an expensive piece of kit. Your employer, the Carver Media Group, is concerned by the happenings at the McMurdo ice station and would like you to produce an accident report detailing the events after you lost control of your hovercraft. The accident report should include a detailed computational model that provides the projected motion of the runaway hovercraft.

from __future__ import division
from vpython import *
from physutil import *
 
#Window setup
scene.width = 1024
scene.height = 768
scene.center = vector(600,0,0)
 
#Objects
cliff = box(pos=vector(-100,0,0), size=vector(200,800,0), color=color.white)
ravine = box(pos=vector(245,-200, 0), size=vector(490,400,0), color=color.white)
lake = box(pos=vector(940, -200, 0), size=vector(900,400,0), color=color.blue)
runawaycraft = sphere(pos=vector(-200,400,0), radius=10, color=color.red)

#Parameters and Initial Conditions
g = vector(0,-9.81,0)
b = 0  #Drag coefficient

runawaycraftm = 1500
runawaycraftv = vector(10,0,0)
runawaycraftp = runawaycraftm*runawaycraftv

#Time and time step
t=0
tf=5
dt = 0.01

#MotionMap/Graph
runawaycraftMotionMap = MotionMap(runawaycraft, tf, 5, markerScale=1, labelMarkerOrder=False, markerColor=color.orange)

#Calculation Loop
while runawaycraft.pos.x < 0:
    rate(500)

    Fgrav = runawaycraftm*g
    Fground = -Fgrav
    Fnet = Fgrav + Fground
 
    runawaycraftp = runawaycraftp + Fnet*dt
    runawaycraft.pos = runawaycraft.pos + (runawaycraftp/runawaycraftm)*dt

    runawaycraftMotionMap.update(t, runawaycraftp/runawaycraftm)

    t = t + dt

MIME type unknown not supported

MIME type unknown not supported

MIME type unknown not supported

MIME type unknown not supported

MIME type unknown not supported

MIME type unknown not supported

MIME type unknown not supported

MIME type unknown not supported