L0-IntroductionNMsnew
L0-IntroductionNMsnew
NUMERICAL METHODS
Introduction
➢ Observations and experiments are used to model the physical world and the
processes of an engineering system can be determined empirically.
➢ Over years and years of observation and experiment engineers and scientists
have noticed that certain aspects of their empirical studies occur repeatedly.
Such general behaviour can then be expressed as fundamental laws that
essentially embody the cumulative wisdom of past experience.
FD = mg
• Analytical Solution:
FU = −cv
If the parachutist is initially at rest (v=0 at t=0),
dv mg − cv
= using calculus dv/dt can be solved to give the result:
dt m
Independent variable
Dependent variable
dv c
dt
=g− v
m v(t ) =
gm
c
(
1− e −( c / m ) t
)
14
Forcing function Parameters
Analytical Solution
v(t ) =
gm
c
(
1− e −( c / m ) t
) If v(t) could not be solved analytically, then
we need to use a numerical method to solve it
t (sec.) V (m/s)
0 0
2 16.40
4 27.77
8 41.10
10 44.87
12 47.49
∞ 53.39
What if we can’t find an analytical solution?
• How do you get a computer to solve the differential
equation?
• Use numerical methods
v(ti +1 ) − v(ti ) c
= g − v(ti )
ti +1 − ti m
This equation can be rearranged to yield
c
v(ti +1 ) = v(ti ) + [ g − v(ti )](ti +1 − ti )
m
t (sec.) V (m/s)
0 0
2 19.60 ∆t = 2 sec
4 32.00
8 44.82
10 47.97
12 49.96 To minimize the error, use a smaller step size, ∆t
∞ 53.39 No problem, if you use a computer!
Numerical Solution
Assume Δt=2s
t (sec.) V (m/s)
0 0
2 19.60
4 32.00
8 44.82
10 47.97
12 49.96
∞ 53.39
18
Numerical methods