Numerical Methods L1 Ok
Numerical Methods L1 Ok
LECTURE 1
INTRODUCTION
Mathematical Model
Exercises
MATHEMATICAL MODEL
A mathematical model can be broadly defined as a formulation or equation that expresses the
essential features of a physical system or process in mathematical terms.
In a very general sense, it can be represented as a functional relationship of the form
(1.1)
The dependent variable is a characteristic that usually reflects the behavior or state of the system;
The independent variables are usually dimensions, such as time and space, along which the system’s
behavior is being determined;
The parameters are reflective of the system’s properties or composition; and
The forcing functions are external influences acting upon the system.
MATHEMATICAL MODEL
The actual mathematical expression of Eq. (1.1) can range from a simple algebraic relationship to large
complicated sets of differential equations.
For example, on the basis of his observations, Newton formulated his second law of motion, which
states that the time rate of change of momentum of a body is equal to the resultant force acting on it.
The mathematical expression, or model, of the second law is the well-known equation
(1.2)
MATHEMATICAL MODEL
(1.3)
There is no independent variable in this case because we are not yet predicting how acceleration
varies in time or space.
MATHEMATICAL MODEL
(1.3)
if the force on an object and the mass of an object are known, Eq. (1.3) can be used to compute
acceleration.
Because of its simple algebraic form, the solution of Eq. (1.2) can be obtained easily.
However, other mathematical models of physical phenomena may be much more complex, and either
cannot be solved exactly or require more sophisticated mathematical techniques than simple algebra
for their solution.
MATHEMATICAL MODEL
To illustrate a more complex model of this kind, Newton’s second law can be used
to determine the terminal velocity of a free-falling body near the earth’s surface. Our
falling body will be a parachutist.
A model for this case can be derived by expressing the acceleration as the time rate
of change of the velocity and substituting it into Eq. (1.3) to yield:
(1.4)
Thus, the mass multiplied by the rate of change of the velocity is equal to the net force
acting on the body.
If the net force is positive, the object will accelerate. If it is negative, the object will
decelerate.
If the net force is zero, the object’s velocity will remain at a constant level.
MATHEMATICAL MODEL
Next, we will express the net force in terms of measurable variables and parameters.
c is the drag coefficient (kg/s). Thus, the greater the fall velocity, the greater the upward force due to air
resistance.
The parameter c accounts for properties of the falling object, such as shape or surface roughness, that
affect air resistance.
For the present case, c might be a function of the type of jumpsuit or the orientation used by the
parachutist during free-fall.
MATHEMATICAL MODEL
The net force is the difference between the downward and upward force
Equation (1.9) is a model that relates the acceleration of a falling object to the forces acting on it. It
is a differential equation.
However, in contrast to the solution of Newton’s second law in Eq. (1.3), the exact solution of Eq.
(1.9) for the velocity of the falling parachutist cannot be obtained using simple algebraic
manipulation.
Rather, more advanced techniques, such as those of calculus, must be applied to obtain an exact or
analytical solution.
MATHEMATICAL MODEL
For example, if the parachutist is initially at rest (𝑣 = 0 𝑎𝑡 𝑡 = 0), calculus
can be used to solve Eq. (1.9) for
MATHEMATICAL MODEL
Notice that the term in brackets is the right-hand side of the differential equation itself [Eq. (1.9)]. That is,
it provides a means to compute the rate of change or slope of v.
Thus, the differential equation has been transformed into an equation that can be used to determine the
velocity algebraically at ti+1 using the slope and previous values of v and t.
MATHEMATICAL MODEL
If you are given an initial value for velocity at some time ti, you can easily compute
velocity at a later time ti+1.
This new value of velocity at ti+1 can in turn be employed to extend the computation
to velocity at ti+2 and so on. Thus, at any time along the way,
For the next interval (from t = 2 to 4 s), the computation is repeated, with the result:
MATHEMATICAL MODEL
The calculation is continued in a similar
fashion to obtain additional values:
MATHEMATICAL MODEL
The numerical method captures the essential
features of the exact solution.
Because we have employed straight-line
segments to approximate a continuously
curving function, there is some discrepancy.
One way to minimize such discrepancies is to
use a smaller step size.
For example, using 1-s intervals results in a
smaller error, as the straight-line segments
track closer to the true solution.
MATHEMATICAL MODEL
Hand calculations: using smaller and smaller step sizes
is impractical.
With computer, large numbers of calculations can be
performed easily.
Thus, you can accurately model the velocity of the
falling parachutist without having to solve the
differential equation exactly.
A computational price must be paid for a more
accurate numerical result (trade-off between accuracy
& computational effort)
Each halving of the step size leads to a doubling of the
number of computations.
CONSERVATION LAWS & ENGINEERING
Aside from Newton’s second law, there are other major organizing principles in
engineering.
Among the most important of these are the conservation laws.
(1.13)
Although simple, Eq. (1.13) embodies one of the most fundamental ways in which
conservation laws are used in engineering—that is, to predict changes with respect to time.
We give Eq. (1.13) the special name: time-variable (or transient) computation.
CONSERVATION LAWS & ENGINEERING
Aside from predicting changes, another way in which conservation laws are applied is for
cases where change is nonexistent. If change is zero, Eq. (1.13) becomes
Thus, at steady state, the downward and upward forces are in balance, and Eq. (1.15) can
be solved for the terminal velocity
The engineering applications are designed to illustrate how numerical methods are
actually employed in the engineering problem-solving process.
Example. A storage tank contains a liquid at depth y, where y = 0 when the tank is half
full. Liquid is withdrawn at a constant flow rate 𝑄 to meet demands. The contents are
resupplied at a sinusoidal rate 3𝑄𝑠𝑖𝑛2 (𝑡).
Use Euler’s method to solve for the depth y from t = 0 to 10 d with a step size of 0.5 d.
The parameter values are A = 1250 m2 and Q = 450 m3/d. Assume that the initial
condition is y = 0.
EXERCISES
(1) Repeat Example 1.2. Compute the velocity to t = 8 s, with a step size of (a) 1
and (b) 0.5 s. Can you make any statement regarding the errors of the calculation
based on the results?
EXERCISES
(1)
EXERCISES
(2) Rather than the linear relationship of Eq. (1.7), you might choose to model the
upward force on the parachutist as a second order relationship,