Processing math: 100%

Numerical Methods

Outline

Introduction

Linear
Equations

Interpolation

Numerical
Solutions

Computer
Measurement

      

Rocket launch

A model rocket with a mass m is launched with a motor that provides a upward force of Fthrust for 3 seconds. This force can be expressed mathematically as F=FthrustH(3t)ˆz where H(x) is the Heaviside step function. As the rocket uses fuel, its mass decreases. In the example below, the rocket looses half of its mass in the three seconds that it accelerates and then the mass remains constant, m=0.1(2H(3t)t/3H(t3)). Other forces that act on the rocket are gravity mgˆz and a drag force that can be described by,

Ffric=a(vvwind)b(vvwind)|(vvwind)|,

where a and b are constants and vwind is the velocity of the wind which can depend on position and time.

m=  kg Fthrust=  N a=  N s/m b=  N s²/m²

The three components of the wind vectors can be functions of space and time.
vwind,x=  m/s vwind,y=  m/s vwind,z=  m/s

The initial conditions at time t=0 are,
x=  m  y=  m  z=  m  vx=  m/s  vy=  m/s  vz=  m/s

 Numerical 6th order differential equation solver 

dxdt=

vx

dvxdt=

dydt=

vy

dvydt=

dzdt=

vz

dvzdt=

Initial conditions:

t0=

Δt=

x(t0)=

Nsteps

vx(t0)=

Plot:

vs.

y(t0)=

vy(t0)=

z(t0)=

vz(t0)=

z
 
2
4
6
8
0
5
10
15
20
25
30
35
-5

t

t   x   vx   y   vy   z   vz