Unstable System. Error Trying to model the Red Bull Stratosphere Free Fall
I tried to model a problem months ago and still I can't solve, I don't want forget this problem so I would like to share what I tried and if you know what is happening with my graph.
This is the problem is related with the world record free fall:
<PERSON> (m = 75 kg) recently stepped out of a balloon at an elevation of 40,000 m and parachuted to earth, setting a height record for skydivers. Here we simulate his jump from the time he steps out of the balloon until the time he opens the parachute. Assume the drag force FD is of the form:
$ F_D = -k|Vy|Vy $
where $k$ is the drag coefficient and $ V_y $ is the vertical velocity. The effect of elevation on the gravitational constant is minor. 229
However, because air density depends on elevation and he jumped from a very high elevation, the effect of elevation on the drag coefficient k should be accounted for. A crude model for k that takes into account the variation of air density with elevation is
$ k = k_{0}e^{-0.00011y} $
where $k$ is the drag coefficient at elevation $y$, $k_0$ is the drag coefficient at sea level $(y = O)$ and where $y$ is in $m$, and $k$ and $k_0$ are in $kg/m$. Assume $k_0 = 0.16 kg/m$ before the parachute opens.
a. Using a coordinate system where $y$ is positive upward, write the differential equations that represent the velocity $v_y$ and elevation $y$ of <PERSON> during the time between stepping out of the balloon and opening his chute. Include the initial conditions.
b.Simulate and graph <PERSON>'s velocity and elevation as functions of time. At what elevation will he have reached his maximum velocity? What is his maximum velocity? For comparison, measurements during the jump indicated that his maximum velocity was $1358 km/h$ and his free fall from $40,000$ to $2000m$ (when his chute opened) lasted for $260s$.
I tried this: $ a $
Initial Conditions: $ V(0) = 0\frac{m}{s} $
$ FD = -k|Vy|Vy $
$ FD = -kV^2 $
Before the parachute opens:
$ F = ma $
$ F = ma $
$ F = m\frac{dv}{dt} $
$ W + FD = ma $
$ W = mg $
System model:
$m\frac{dv}{dt} = -kv^2 - mg $
$\frac{dv}{dt} = \frac{-kv^2}{m} - g $
$\frac{dv}{dt} = \frac{-0.16e^{-0.00011y} v^2}{m} - g $
$a = \ddot{y} , v = \dot{y} , y=y.$
2nd order differential equation:
$ \ddot{y}= \frac{-0.16e^{-0.00011(y-y0)} \dot{y}^2}{m} - g $
This is my block diagram:
This graphic show my error:
What I am forgetting?
Bibliography:
A First Course in Differential Equations, Modeling, and Simulation - <PERSON>, <PERSON> 771