0% found this document useful (0 votes)
29 views17 pages

CMP 2

This document discusses numerical methods for solving differential equations that cannot be solved analytically. It provides an example of using Euler's method to model the velocity of a falling parachutist. The method approximates the derivative as the slope between discrete time steps. Using a time step of 2 seconds, the numerical solution tracks the essential features of the exact solution, though with some discrepancy due to the linear approximations between time steps. Using a smaller time step can minimize this discrepancy. Students are assigned related problems applying these concepts to other systems modeled by differential equations.
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PPTX, PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
29 views17 pages

CMP 2

This document discusses numerical methods for solving differential equations that cannot be solved analytically. It provides an example of using Euler's method to model the velocity of a falling parachutist. The method approximates the derivative as the slope between discrete time steps. Using a time step of 2 seconds, the numerical solution tracks the essential features of the exact solution, though with some discrepancy due to the linear approximations between time steps. Using a smaller time step can minimize this discrepancy. Students are assigned related problems applying these concepts to other systems modeled by differential equations.
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PPTX, PDF, TXT or read online on Scribd
You are on page 1/ 17

COEP Technological University

Department of Mechanical Engineering

Computational Methods &


Programming

Welcome!!
Session 2
Summary of last session..
A simple mathematical model
F=ma (1)
a=F/m (2)
dv / dt = F / m (3)
F = FD + FU (4)
FD = mg (5)
FU = cv (6)
but F = mg – cv
dv / dt = ( mg – cv ) / m (7)
dv / dt = g – (c/ m) * v (8)

After a sufficiently long time, a constant velocity,


called the terminal velocity, of 53.39 m/s is reached.
Limitations of exact approach..
• As we studied that, solving the mathematical models using
the exact approach involves the calculus.

• To find out rate of change of velocity w.r.t. time, we applied


the integration to the mathematical model which was
merely consist of fewer factors/parameters.

• But in real-time engineering practice, unfortunately, there


are many mathematical models that cannot be solved
exactly.

• The computation of calculus becomes tedious and time


consuming.
Introduction to
numerical/approximation
approach
A numerical/approximation approach:
• As mentioned previously, numerical methods are those in
which the mathematical problem is reformulated so it can
be solved by arithmetic operations.
Approximated
Velocity response
Exact / True
response
B
v(ti + 1)
dv/dt

∆𝒗
A
v(ti)

ti ∆𝒕 ti + 1 time
Continue…
• So you might come with the doubt that, why is
approximated response is considered as straight line???

• Remember that slope depends on time interval

• Hence by considering the slope and previous values of v


and t, value of v(ti + 1) can be easily determined

• If you reduce the step size of time, automatically the


accuracy of the response of velocity will increase

• Hence lesser the value of step size, more accurate will be


the velocity at v(ti + 1)
Continue…
• Also remember from calculus that

• Equation (10) represents the reverse process.

(10)

• Equation (10) is called a finite divided difference approximation


of the derivative at time ti .
Continue:

dv / dt = g – (c/ m) * v (8)
• Notice that the term in brackets is the right-hand side of the
differential equation itself [Eq. 8].

• That is, it provides a means to compute the rate of change or


slope of v.

• Thus, to transform the differential equation into simple


algebraic equation Eq. 10 can be substituted into Eq. 8, we
get
Continue:

• This equation can then be rearranged as we aim to calculate


the velocity at time instance v (ti+1).

• Hence term v (ti+1) can only be kept on left hand side and by
rearranging other terms we get,

(11)

• 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.
Continue:

• 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,

• New value = old value + slope × step size

• Note that this approach is formally called Euler’s method.


Example 2: Numerical Solution to the Falling Parachutist Problem
• Problem Statement. Perform the same computation as in
Example 1 but use Eq. (11) to compute the velocity. Employ a
step size of 2 s for the calculation. Solution. At the start of the
computation (ti = 0), the velocity of the parachutist is zero.

• Using the same values from Example 1, Eq. (11) can be used
to compute velocity at ti+1 = 2 s:
Example 2: Numerical Solution to the Falling Parachutist Problem
Example 2: Numerical Solution to the Falling Parachutist Problem
• It can be seen that the numerical method captures the essential features
of the exact solution.

• However, because we have employed straight-line segments to


approximate a continuously curving function, there is some discrepancy
between the two results.

• One way to minimize such discrepancies is to use a smaller step size. For
example, applying Eq. (11) at l-s intervals results in a smaller error, as the
straight-line segments track closer to the true solution.

• However, with the aid of the 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.
Assignment for the 1st week
1. (A & B batch) Use calculus to solve Eq. (8) for the case where the initial
velocity, v(0) is non zero.

2. (C & D batch) Rather than the linear relationship of Eq. (6), you might
choose to model the upward force on the parachutist as a second order
relationship,
FU = − c′ *
where c′ = a second-order drag coefficient (kg/m).
a) Using calculus, obtain the closed-form solution for the case where the
jumper is initially at rest (v = 0 at t = 0).
b) Repeat the numerical calculation in Example 3 with the same initial
condition and parameter values. Use a value of 0.225 kg/m for c′.

3. (G & H batch) Compute the velocity of a free-falling parachutist using


Euler’s method for the case where m = 80 kg and c = 10 kg/s. Perform the
calculation from t = 0 to 20 s with a step size of 1 s. Use an initial condition
that the parachutist has an upward velocity of 20 m/s at t = 0. At t = 10 s,
assume that the chute is instantaneously deployed so that the drag
coefficient jumps to 50 kg/s.
Assignment for the 1st week
5. (E & F batch) 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 Q to meet
demands. The contents are resupplied at a sinusoidal rate 3Q*(t).
The above equation can be written for this system as

or, since the surface area A is constant,

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 = 1200 m2 and Q = 500 m3/d. Assume
that the initial condition is y = 0.
Submission mode for submission:
• You can submit either handwritten or typed scanned copy (PDF
compulsory) of your solution.

• The name of PDF should be your MIS number (Do not write name or
division).

• One student should make the power point presentation of a solution.


(He/she shall be excluded from submission of handwritten or typed
scanned copy)

• In addition to that, another one student should deliver the same


presentation and capture the video. (He/she shall be excluded from
submission of handwritten or typed scanned copy)

• The link for submission shall be shared to your CR.


Thank you!

You might also like