Chapter1
Chapter1
ESOF 3558
Numerical Methods & Modeling
Part 1
Chapter 1: Mathematical Modeling, Numerical
Methods, and Problem Solving
All images copyright © The McGraw-Hill Companies, Inc. Permission required for reproduction or display.
1
2025-01-07
Chapter Objectives
• Define numerical methods and their problem-solving role.
• Develop models of physical systems based on conservation laws and
scientific principles.
• Learn how numerical methods generate solutions.
• Examine the difference between steady-state and dynamic solutions.
• Provide an overview various numerical methods in the book.
2
2025-01-07
Model Function
• Dependent variable
• Reflects the behavior or state of the system
• Independent variables
• Dimensions, like time and space, determining the system’s behavior
• Parameters
• Constants reflective of the system’s properties
• Forcing functions
• External influences acting upon the system
3
2025-01-07
Analytical Solution
• If the jumper starts at rest, calculus
provides a closed-form solution:
𝑔𝑚 𝑔𝑐
𝑣 𝑡 𝐭𝐚𝐧𝐡 𝑡
𝑐 𝑚
• A computer can be used to generate a
graph of the system.
• Example: graph representing 𝑣 𝑡 with
𝑣 0 0, 𝑚 68.1 kg, 𝑐 0.25 kg/m,
𝑔 = 9.81 m/s2
Numerical Methods
• To solve the problem numerically, approximate the rate of change of
velocity as its average change over a small-time interval:
𝑑𝑣 𝑐
𝑔 𝑣
𝑑𝑡 𝑚
𝑑𝑣 Δ𝑣 𝑣 𝑡 𝑣 𝑡
lim lim
𝑑𝑡 ∆ → Δ𝑡 → 𝑡 𝑡
𝑑𝑣 Δ𝑣 𝑣 𝑡 𝑣 𝑡
𝑑𝑡 Δ𝑡 𝑡 𝑡
4
2025-01-07
Euler's Method
• Substituting the finite difference equation into the differential
equation gives:
| 𝑔 𝑣 𝑡
• After rearranging:
𝑣 𝑡 𝑣 𝑡 𝑔 𝑣 𝑡 𝑡 𝑡 “Euler’s method”
New value = Old value + Slope Step
t0 t1 t2 t3 … ti ti+1 t
9
Numerical Results
• Apply Euler's method with ∆𝑡 2s intervals for numerical solution
• Numerical solution captures key features of the exact solution
• Improve numerical solution by using smaller steps
Approximation 𝑐
Error 𝑣 𝑡 𝑣 𝑡 𝑔 𝑣 𝑡 𝑡 𝑡
𝑚
Initial condition: 𝑣 0 0
𝑚 68.1 kg
𝑐 0.25 kg/m
𝑔 = 9.81 m/s2
10
5
2025-01-07
• Predicting Changes:
• Conservation laws predict changes over time Transient
11 11
12
6
2025-01-07
13