• Use vectors to mathematically determine:
  • Relative velocity within inner tube-river system
  • Launch angle
• Use code to model relative motion

• Relative velocity
• Vector addition
• Vector decomposition
  • $\sin$ and $\cos$ functions
  • $\arctan$ function

• Interpret code
• Modify variables within existing code

Part 1
It is a hot summer day, and you are meeting your friends for a picnic lunch down by the river. When you get to the river, you discover you are on the wrong side, 20 meters from the picnic site. As an experienced inner tuber, you know that you can paddle your tube at a speed of 2.0 m/s. Additionally, you know the river moves westward at a speed of 1.5 m/s. Your friends have already started the grill and your meal will be ready in 15 seconds.

You have an app that lets you enter your velocity and the velocity of the river to determine your net velocity. For the following questions, calculate the answer and check your results with the code.

1. If you were to get in the river and begin paddling directly across, what velocity (magnitude and direction) do your friends see you moving relative to the shore?
2. How long does it take to cross the river? If the river had a faster current, how would this affect the time it takes you to get across?

Part 2
You would like to paddle across the river to get to your friends at the picnic site, but you don't want to hike back upriver once you paddle across.

1. At what angle with the shore must you paddle in order to get directly across the river to the picnic area?
2. How long will this take? Why is this different from your answer in Part 1?

Part 3
After lunch, you decide to cool off by floating in the river. You want to chat with your friends on the shore while you float, so you need to paddle in such a way that you remain stationary.

1. How might you modify your code so that you model the scenario of floating in the center of the river, in front of the picnic site?

Part 4
Coding questions:

1. What is a while loop, and what quantities/values must go inside a while loop?
2. What does “dt” stand for (line 37)? What happens when you make “dt” larger?
3. Why is the width of the river divided in half in “tube = ring(pos=vector(0,-w/2,0)…” (line 16)?
4. What is the purpose of “w/2” in “while tube.pos.y<w/2:” (line 46)?
5. If you modified you code correctly in Part 3, the text “time to cross” and “tube velocity” will not appear. Why?

Link

GlowScript 2.7 VPython
get_library('https://rawgit.com/perlatmsu/physutil/master/js/physutil.js')
 
#Window setup
scene.range = 20
scene.width = 900
scene.height = 500
scene.background=vector(0,0.7,0.1)
 
 
#Objects
w=20
river = box(pos=vector(0,0,-1), length=200, height=w, width=1, color=color.blue, texture=textures.rough)
#southshore = box(pos=vector(0,-100,0), length=600, height=50, width=1, color=vector(0.4,0.6, 0.2))
#northshore = box(pos=vector(0,100,0), length=600, height=50, width=1, color=vector(0.4,0.6, 0.2))
tube = ring(pos=vector(0,-w/2,0), axis=vector(0,0,180), radius=2, thickness=0.8, color=color.black, make_trail=False,)
#target= text(text='X', align='center', color=color.red)
label(pos=vec(tube.pos.x,w/2,0), text='picnic site', color=color.black )
label(pos=vec(-10,-10,0), text='click to run', color=color.black )
 
 
#Parameters and Initial Conditions
tubes=0
riverv = vector(0.5,0,0)
degrees=0
theta=degrees*pi/180
tubesx=tubes*sin(theta)
tubesy=tubes*cos(theta)
tubev=vec(tubesx,tubesy,0)
scale=1.0
 
 
 
#Time and time step
t = 0
tf = 3
dt = 0.01
 
tubeMotionMapx = MotionMap(tube, tf, 2, markerScale=scale, markerColor=color.blue, labelMarkerOrder=False)
tubeMotionMapy= MotionMap(tube, tf, 2, markerScale=scale, markerColer=color.red, labelMarkerOrder=False, )
tubeMotionMapr= MotionMap(tube, tf, 2, markerScale=scale, markerColor=color.green, labelMarkerOrder=False)
 
ev = scene.waitfor('click')
 
#Calculation Loop
while tube.pos.y<w/2: 
    rate(500)
 
    tube.pos = tube.pos + tubev*dt
    tube.pos = tube.pos + riverv*dt
 
#Set up motion map
    relv = riverv + tubev
    tubeMotionMapr.update(t, relv)
    tubeMotionMapx.update(t, tubev)
    tubeMotionMapy.update(t, riverv)
 
 
    t = t + dt
 
#determine angle
angle=atan(mag(riverv)/mag(tubev))*(180/pi)
 
#label
print("time to cross =", t , "s")
print("tube velocity=",mag(relv),"m/s at", angle, "degrees")

Part 1

1. You can determine your velocity with respect to the shore through vector addition:
   1. $\vec{v}_{tube/shore}=\vec{v}_{river/shore}+\vec{v}_{tube/river}$
   2. $\vec{v}_{tube/shore}= < -1.5, 0, 0>+<0, 2, 0> = < -1.5, 2, 0>$ m/s
      1. $\text{Magnitude}=\sqrt{(-1.5)^2+(2^2)}=\sqrt{6.25}=2.5$ m/s
      2. $\text{Angle}=\arctan(\dfrac{1.5}{2})=36.9°$
2. The movement of the tube in the y-direction is completely independent of the river's movement in the x-direction. Because of this, we just have to look at the y-coordinate of the tube's velocity in order to calculate how long it takes to cross
   1. $\text{Time}=\dfrac{\text{Distance}}{\text{Velocity}}$
   2. $t=\dfrac{20\text{m}}{2.0\text{m/s}}=10\text{s}$

• Ball Launch
• Block on a Ramp