Lecture Slides AMM Week 1 - Introduction
Lecture Slides AMM Week 1 - Introduction
Dr Anna Huber
Applied Mathematical Methods
Dr Anna Huber
Stephenson Building Room IC 0.08
A.Huber@tees.ac.uk
timetable on MyTU
https://mytu.tees.ac.uk/
Tutorials (Seminars): 1h
Alex, PT tutorials
Behnam, mechanical
Kin, chemical
Teaching Team
Anna, lectures
Roshan, civil
Lectures: Tuesday, 9-10:30
Tutorials (Seminars): 1h
⚫ Important announcements
⚫ Learning materials: Slides, recordings, …
⚫ Tutorial sheet
• Reading list
• Teaching schedule
Academic Support Hub
Dr Anna Huber
• Tuesday 12:00 - 13:00, Room IC0.08 or IC1.43
Dr Anna Huber
Teaching
5 Numerical Integration
6 EXAM (ICA, 40%, online)
Schedule 7
8
Roots of Equations
Optimisation
9 Independent Study Week
10 Probability
11 Statistics
12 Revision
13/14 EXAM (ECA, 60%, on Campus)
Mathematical
Modelling
The Modelling Cycle
𝐹 = 𝐹𝐷 + 𝐹𝑈
𝐹𝐷 = 𝑚𝑔
𝐹𝑈 = −𝑐𝑣
𝐹=0
𝑚𝑔 − 𝑐𝑣 = 0
21
Mathematical Model
• Model is represented as a functional relationship of the form:
𝐹 = 𝐹𝐷 + 𝐹𝑈
𝐹𝐷 = 𝑚𝑔
𝐹𝑈 = −𝑐𝑣
𝐹=0
𝑚𝑔 − 𝑐𝑣 = 0
23
Velocity
Some complex models may not be solved
exactly or require more sophisticated
mathematical techniques than simple • This is a differential equation and is
algebra written in terms of the differential
dv F rate of change dv/dt of the variable
=
dt m that we are interested in predicting.
F = FD + FU
FD = mg
FU = −cv dv c
dv mg − cv =g− v
= dt m
dt m
24
What do we want to solve in a differential
equation?
• What do we know ?
Differential Equations
An equation which contains a dependent
variable and its derivatives is called a
Differential Equation.
• Examples:
Differential Equations—DE’s
An equation which contains a dependent
variable and its derivatives is called a
Differential Equation.
F = FD + FU
• If the parachutist is initially at rest (𝑣 = 0 at
FD = mg
𝑡 = 0), using calculus 𝑑𝑣/𝑑𝑡 can be solved to give
FU = −cv
the result:
dv mg − cv
= Dependent Independent
dt m variable variable
dv
dt
c
=g− v
m
v(t ) =
gm
c
1− e (
−( c / m ) t
)
Parameters
Forcing function
31
Analytical Solution
( )
gm This formula can be obtained analytically.
−( c / m ) t
v(t ) = 1− e But it may not always be possible.
Then we would need to use a numerical
c 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
32
Numerical Solution
dv v v(ti +1 ) − v(ti ) dv v
= ........ = lim
dt t ti +1 − ti dt t →0 t
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 = 2 sec
t (sec.) V (m/s)
0 0
2 16.40
4 27.77
8 41.10
10 44.87
12 47.49
∞ 53.39
34
Numerical Solution
dv v v(ti +1 ) − v(ti ) dv v
= ........ = lim
dt t ti +1 − ti dt t →0 t
v(ti +1 ) − v(ti ) c
= g − v(ti )
ti +1 − ti m
This equation can be rearranged to
c
v(ti +1 ) = v(ti ) + [ g − v(ti )](ti +1 − ti )
m
t (sec.) V (m/s) ∆t = 2 sec
0 0
2 19.6
g = 9.8 m/s2
4 32
c =12.5 kg/s
8 39.85 m = 68.1 kg
10 44.82
12 49.96 To minimize the error, use a smaller step size, ∆t
No problem, if you use a computer!
∞ 53.39 35
Summary
• Many mathematical models cannot be solved analytically
• Numerical solutions can be found that approximate the
exact solution
Numerical Methods
Question Example Analytical Method Numerical
Method
𝑓 ??? Velocity for Solve Differential Differential
parachutist Equation Equations
Weekly
Feedback
Opportunity