Applications of Differential Equations in Engineering
Applications of Differential Equations in Engineering
OF
NUMERICAL
ANALYSIS
(MTH 204)
Submitted To:
Ms. Ritu Sharma Submitted By:
Ishan Arora
B.Tech-M.B.A (IT)
Section: D2702
Roll No. : 21
CONTENTS
Nomenclature
The theory of differential equations is quite developed and the methods used to study
them vary significantly with the type of the equation.
• An ordinary differential equation (ODE) is a differential equation in which the
unknown function (also known as the dependent variable) is a function of a
single independent variable. In the simplest form, the unknown function is a real
or complex valued function, but more generally, it may be vector-valued or
matrix-valued: this corresponds to considering a system of ordinary differential
equations for a single function. Ordinary differential equations are further
classified according to the order of the highest derivative with respect to the
dependent variable appearing in the equation. The most important cases for
applications are first order and second order differential equations. In the classical
literature also distinction is made between differential equations explicitly solved
with respect to the highest derivative and differential equations in an implicit
form.
Both ordinary and partial differential equations are broadly classified as linear and
nonlinear. A differential equation is linear if the unknown function and its derivatives
appear to the power 1 (products are not allowed) and nonlinear otherwise. The
characteristic property of linear equations is that their solutions form an affine subspace
of an appropriate function space, which results in much more developed theory of linear
differential equations. Homogeneous linear differential equations are a further subclass
for which the space of solutions is a linear subspace i.e. the sum of any set of solutions or
multiples of solutions is also a solution. The coefficients of the unknown function and its
derivatives in a linear differential equation are allowed to be (known) functions of the
independent variable or variables; if these coefficients are constants then one speaks of a
constant coefficient linear differential equation.
There are very few methods of explicitly solving nonlinear differential equations; those
that are known typically depend on the equation having particular symmetries. Nonlinear
differential equations can exhibit very complicated behavior over extended time intervals,
characteristic of chaos. Even the fundamental questions of existence, uniqueness, and
extendability of solutions for nonlinear differential equations, and well-posedness of
initial and boundary value problems for nonlinear PDEs are hard problems and their
resolution in special cases is considered to be a significant advance in the mathematical
theory (cf Navier–Stokes existence and smoothness).
Examples
In the first group of examples, let u be an unknown function of x, and c and ω are known
constants.
In the next group of examples, the unknown function u depends on two variables x and t
or x and y.
We present examples where differential equations are widely applied to model natural
phenomena, engineering systems and many other situations.
Radioactive Decay
Many radioactive materials disintegrate at a rate proportional to the amount present. For
example, if X is the radioactive material and Q(t) is the amount present at time t, then the
rate of change of Q(t) with respect to time t is given by
where r is a positive constant (r>0). Let us call the initial quantity of the
material X, then we have
Clearly, in order to determine Q(t) we need to find the constant r. This can be done using
what is called the half-life T of the material X. The half-life is the time span needed to
Example: A radioactive isotope has a half-life of 16 days. You wish to have 30 g at the
end of 30 days. How much radioisotope should you start with?
Solution: Since the half-life is given in days we will measure time in days. Let Q(t) be
the amount present at time t and the amount we are looking for (the initial amount).
We know that
Hence, since
we get
Hence,
which implies
This equation makes it possible to find k if the interval of time is known and vice-
versa.
Example: Time of Death Suppose that a corpse was discovered in a motel room at
midnight and its temperature was . The temperature of the room is kept constant at
. Two hours later the temperature of the corpse dropped to . Find the time of
death.
Solution: First we use the observed temperatures of the corpse to find the constant k. We
have
In order to find the time of death we need to remember that the temperature of a corpse at
time of death is (assuming the dead person was not sick!). Then we have
where m and C are parameters. Clearly, we may change the names of the variables and
still have the same geometric curves. For example, the above families define the same
geometric object as
Note that the first family describes all the lines passing by the origin (0,0) while the
second the family describes all the circles centered at the origin (including the limit case
when the radius 0 which reduces to the single point (0,0)) (see the pictures below).
and
In this page, we will only use the variables x and y. Any family of curves will be written
as
One may ask whether any family of curves may be generated from a differential
equation? In general, the answer is no. Let us see how to proceed if the answer were to be
yes. First differentiate with respect to x, and get a new equation involving in general x, y,
, and C. Using the original equation, we may able to eliminate the parameter C from
the new equation.
Since we have
then we get
You may want to do some algebra to make the new equation easy to read. The next step is
to rewrite this equation in the explicit form
Example. Find the differential equation (in the explicit form) satisfied by the family
Answer. We have already found the differential equation in the implicit form
Definition. Consider two families of curves and . We say that and are
orthogonal whenever any curve from intersects any curve from , the two curves are
orthogonal at the point of intersection.
For example, we have seen that the families y = m x and are orthogonal.
One may then ask the following natural question:
Consider the family of curves . We assume that an associated differential equation may
be found, say
We know that for any curve from the family passing by the point (x,y), the slope of the
tangent at this point is f(x,y). Hence the slope of the line perpendicular (or orthogonal) to
this tangent is which happens to be the slope of the tangent line to the
orthogonal curve passing by the point (x,y). In other words, the family of orthogonal
curves are solutions to the differential equation
From this we see what we have to do. Indeed consider a family of curves . In order to
find the orthogonal family, we use the following practical steps
Step 3. Write down the differential equation associated to the orthogonal family
Step 4. Solve the new equation. The solutions are exactly the family of
orthogonal curves.
Step 5. You may be asked to give a geometric view of the two families. Also you
may be asked to find a specific curve from the orthogonal family (something like
an IVP).
which gives
We recognize the family of lines and we confirm our earlier observation (that the two
families are indeed orthogonal).
This example is somehow easy and was given here to illustrate the technique.
We recognize an homogeneous equation. Let us use the technique developed to solve this
kind of equations. Consider the new variable (or equivalently y = x z). Then we
have
and
Hence we have
which gives z=0. The non-constant solutions are found once we separate the variables
Before we perform the integration for the left-hand side, we need to use partial
decomposition technique. We have
We will leave the details to you to show that A = 1, B=-2, and C=0. Hence we have
Hence
which is equivalent to
which is equivalent to
We recognize a family of circles centered on the y-axis and the line y=0 (the x-axis which
was easy to guess, isn't it?)
If we put both families together, we appreciate better the orthogonality of the curves (see
the picture below).
Population Dynamics
Here are some natural questions related to population problems:
More can be said about the problem but, in this little review we will not discuss them in
detail. In order to illustrate the use of differential equations with regard to this problem
we consider the easiest mathematical model offered to govern the population dynamics of
a certain species. It is commonly called the exponential model, that is, the rate of change
of the population is proportional to the existing population. In other words, if P(t)
measures the population, we have
,
where the rate k is constant. It is fairly easy to see that if k > 0, we have growth, and if k
<0, we have decay. This is a linear equation which solves into
• if k>0, then the population grows and continues to expand to infinity, that is,
• if k<0, then the population will shrink and tend to 0. In other words we are facing
extinction.
Clearly, the first case, k>0, is not adequate and the model can be dropped. The main
argument for this has to do with environmental limitations. The complication is that
population growth is eventually limited by some factor, usually one from among many
essential resources. When a population is far from its limits of growth it can grow
exponentially. However, when nearing its limits the population size can fluctuate, even
chaotically. Another model was proposed to remedy this flaw in the exponential model. It
is called the logistic model (also called Verhulst-Pearl model). The differential equation
for this model is
where M is a limiting size for the population (also called the carrying capacity). Clearly,
when P is small compared to M, the equation reduces to the exponential one. In order to
solve this equation we recognize a nonlinear equation which is separable. The constant
solutions are P=0 and P=M. The non-constant solutions may obtained by separating the
variables
and integration
The partial fraction techniques gives
which gives
which, once substituted into the expression for P(t) and simplified, we find
However, this is still not satisfactory because this model does not tell us when a
population is facing extinction since it never implies that. Even starting with a small
population it will always tend to the carrying capacity M.
Some other Applications to Engineering and
Sciences
Historically, it has been the needs of the physical sciences which have driven
the development of many parts of mathematics, particularly analysis. The
applications are sometimes difficult to classify mathematically, since tools
from several areas of mathematics may be applied. We focus on these
applications not by discussing the nature of their discipline but rather their
interaction with mathematics.