6
$\begingroup$

Help! I finished my numerical simulation of falcon 9 launch in google sheets, where I calculate all the value changes with time. But, when I looked at real launch to compare the values, I got really big, and non-understandable for me difference in values

On 5 second of flight, my rocket had velocity of (75,6) kmph, when on the video it only reached 26 kmph On 10 second values were 156 to 95 kmph. Here is the link to the launch video I used to compare values :

And my numerical simulations, with formulas inside : https://docs.google.com/spreadsheets/d/1oBAyEQGJUU6Kpze6J2j8qlJ9HxCKDXyiTmY2s5b2q_E/edit?usp=sharing

I just cant stand where is the problem... All the equations seem to be correct, and values too...

 Time          Predicted          Launch video       
              H (m)     v (m/s)    H         v
  0 s         0         0          0         0 
2.5 s        15        10.7        -         2
  5 s        58        21.5        -         7.8
7.5 s        124.7     32.5        -        16.5 
 10 s        221.5     44          -        26.3

*On the video the height is too low to be displayed accurately

$\endgroup$
  • 1
    $\begingroup$ I'm not able to view that drive link. Could you include the relevant values and formulas of your simulation in the post? $\endgroup$ – SE - stop firing the good guys Oct 4 at 10:41
  • 1
    $\begingroup$ That was the wrong link, my bad. Edited and attached right one $\endgroup$ – mad.redhead Oct 4 at 10:50
  • 3
    $\begingroup$ You are doing so many things wrong I barely where to start. As a starting point, (1) You are ignoring gravity. (2) This is a three dimensional problem. You only have one. (3) Your use of units is inconsistent. (4) Your time step is far too large. (5) Your atmosphere model is very wrong. $\endgroup$ – David Hammen Oct 4 at 10:59
  • $\begingroup$ Im not ignoring gravity at this point. I know that I need to change (g) and thrust with altitude, but I mean the start altitudes are low, I ignore their change. What do you mean under "three dimensional problem"? In my case I wanted only to simulate 2d launch. Rocket goes up or falls down. $\endgroup$ – mad.redhead Oct 4 at 11:01
  • 1
    $\begingroup$ A rocket that launches things into orbit has to achieve two things. One is to take the payload above the bulk of the Earth's atmosphere. (This was point #5, "your atmosphere model is very wrong".) The other is to increase the horizontal component of velocity from less than half of a kilometer per second (the Earth's rotational velocity at the equator) to about 7.8 kilometers per second. (This was point #2, "this is a three dimensional problem".) Rockets that launch things into space do not go straight up, $\endgroup$ – David Hammen Oct 4 at 11:12
8
$\begingroup$

File:Au pfeffernuesse 01.jpg

Source (cropped)

Why do plastic rulers make you measure between two lines, instead of starting at the edge?

tl;dr: Never use the edge of a ruler or the "zero" of any analog measuring device when you are making careful measurements if you can help it. Instead try measuring between two equivalent demarkations.

from this answer; doesn't start moving until T+00:00:07!

Better yet, measure several times, plot it, and draw a line through your measurements and determine the slope.

Newton's $F = ma$ means that $a = F/m$. But $F$ has to be the total force, and there are two big ones in opposite directions, and that's what lifting off the Earth is all about!

acceleration pointing up

This is the thrust of the rocket and Wikipedia's Falcon 9 says that for F9 v1.2 or Full Thrust it is +7,607,000 Newtons where the plus sign means "up".

The loaded mass is 549,054 kg, so

$a_{thrust} = F_{thrust}/m = \text{13.85 m/s}^2$

acceleration pointing down

At liftoff and while it's not going really really fast yet, so gravity is the only important one. Earth's standard gravityis 9.80665 m/s2.

So the rocket will initially accelerate at roughly 13.85 - 9.81 = 4.04 m/s2.

Using $v = at$ and assuming the mass doesn't change too much in the first ten seconds (it does a bit of course) then I get 72.7 and 145.4 kph after 5 and 10 seconds, just like you.

what could be the difference?

Never start a measurement at zero. Never use the zero of a ruler; if something is about 10 cm long, measure from 1 cm to 11 cm, that way if the manufacturer screwed up slightly cutting it it won't make a difference.

Play the video slowly, like 1/4 or 1/8 speed, and stop it each time the timer changes from one seconds digit to the next, and write down the velocity.

Try this:

Plot v vs t on a graph, and you'll see it is probably a fairly straight line but it hits zero at around 1 or 1.5 seconds, not zero.

Then read all the answers and all the comments on What "actually" happens at T-minus-0


update:

You're doing just fine!

I see you are active in Stack Overflow so I'll add some Python script. I took your numbers and did two things:

  1. you have the h and v switched in the 7.5 second line so I switched them back
  2. included a 3.5 second offset between whatever T=0 is in the video and time since the rocket starts moving.
  3. plotted your simulation as solid black line
  4. plotted your video measurements as blue dots
  5. Added a red dashed line for a simple 4.04 m/s2 acceleration

The resulting fit is quite nice; all three agree fairly well already! You are doing just fine and you can extend this further in time now, and check more effects.

quick calc

import numpy as np
import matplotlib.pyplot as plt

time = np.array([0, 2.5, 5, 7.5, 10])
time_smooth = np.arange(0, 16, 0.1)

model_h = np.array([0, 15, 58, 124.7, 221.5])
model_v = np.array([0, 10.7, 21.5, 32.5, 44])

video_v = np.array([0, 2, 7.8, 16.5, 26.3])

T_start = 3.5

a_guess = 4.04

plt.figure()
plt.subplot(2, 1, 1)
plt.plot(time +T_start, model_v, '-k')
plt.plot(time, video_v, 'ob')
plt.plot([-1, 14], [0, 0], '-k', linewidth=0.5)
plt.plot(time_smooth, a_guess * (time_smooth - T_start).clip(0, None), '--r')
plt.xlim(-1, 14)
plt.ylim(-2, 42)
plt.title('vertical speed (m/s)')
plt.subplot(2, 1, 2)
plt.plot(time + T_start, model_h, '-k')
plt.plot(time_smooth, 0.5 * a_guess * (time_smooth - T_start).clip(0, None)**2, '--r')
plt.plot([-1, 14], [0, 0], '-k', linewidth=0.5)
plt.xlim(-1, 14)
plt.ylim(-10, 210)
plt.title('vertical height (m)')
plt.xlabel('time (sec)')
plt.show()
| improve this answer | |
$\endgroup$
  • 1
    $\begingroup$ Edited post with tables for 0, 2.5, 5, 7.5 and 10 seconds $\endgroup$ – mad.redhead Oct 4 at 12:43
  • 1
    $\begingroup$ @OrganicMarble it's simply an extreme example of what I've explained. The answer is rocket-agnostic; check the sentence directly above the video. (basically T=0 doesn't mean anything, the rocket starts moving when it decides to start moving) $\endgroup$ – uhoh Oct 4 at 15:46
  • 1
    $\begingroup$ Wow... I really surprised you put THAT much effort in explaining some simple (at first look) things in rocket science to a newbie, like me... The world is not without good people :). Im going to look thru this, and write you, if I`ll have more issues. But with answer like this, I guess i wont have any problems :). Thanks also for including code, im familliar with that library for plots. **Im still amazed, and very grateful!!! $\endgroup$ – mad.redhead Oct 4 at 17:58
  • 2
    $\begingroup$ "Play the video slowly, like 1/4 or 1/8 speed, and stop it each time the timer changes from one seconds digit to the next, and write down the velocity.". Youtube lets you even step through videos frame-by-frame with . and , ;) (forward and reverse) $\endgroup$ – Polygnome Oct 4 at 19:11
  • 1
    $\begingroup$ @mad.redhead I don't know if I'm a good person but I'm basically just addicted to Stack Exchange ;-) You're welcome to follow up here in comments with more clarifications on this question, but if you have a new issue I recommend that you ask about it by posting a new question and linking back to your old question for background. That way it has maximum visibility to everyone which in turn maximizes the quantity/quality/diversity of answers. $\endgroup$ – uhoh Oct 4 at 22:54

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.