As every engineer or engineering student knows, almost everything in the world can be translated into numbers, this is especially helpful with physical problems, because using several tools that math gives us, we can find solutions that can be reflected in the real life later, and thus the problems are solved. These tools come from knowing in extend some branches of math like, algebra, geometry, trigonometry, calculus, etc.

In this post I´ll talk about a very useful aspect of calculus and it´s applications, the Ordinary Differential Equations (ODE) of first order, which are expressions where the unknown variable is a function, and include only the first derivative of said function (Thus of first order). Of course, there are ODEs of second or even higher orders (where the second, third or so on derivatives can be present in the equation) and Partial differential equations (PDE) that depend on several variables, but those are a topics a bit more complex.

https://steemitimages.com/DQmeqn2wvL6EA8ZcthiwxoESLiXKR2pbvwkahfNwPoin8C3/435px-GodfreyKneller-IsaacNewton-1689.jpg Isaac Newton, Father of the Differential calculus by Sir Godfrey Kneller

First we have to understand the concept of Modeling in engineering, this is when we take a problem and we give it a numerical representation, this is known as mathematical model, once we have the model ready, we proceed to find a mathematical solution, and later we give the solution a physical interpretation. Most programs and calculators can give you a help in finding a numerical answer, but the trick is in creating the mathematical model, and bringing the answer to real life.

To fully understand the equations and principles that I’ll discuss here, you´ll need some basic knowledge of Integral calculus, more specifically an understanding of derivatives (their real world, or physical applications) and some methods of integration.

A First-order ODE can be written in an explicit form like this

http://quicklatex.com/cache3/1b/ql_6f101cf92529418a0969f3501ec20f1b_l3.png

A general solution for this is any function i.e. g(x) that on some open interval a < x < b is defined and differentiable in that interval, and is such that replacing g and g’ in y and y’ the equations becomes an identity.

There are some ODEs that can be reduced to with common algebra to this form

http://quicklatex.com/cache3/9f/ql_f5a0e999bc981bb04328b9c6cea6fd9f_l3.png

Where we can further develop and integrate each side, because in calculus y’ = dy/dx

http://quicklatex.com/cache3/55/ql_0a6de4340131e8d948f6cb14b05c3355_l3.png

We should note that each integral throws a different constant “c”, one in the right and another in the left side of the equation, but these constants can be considered together as one. This method is called separating variables, because the variables are in different sides, “y” only appear at the left side, and “x” in the right.

Let´s do an example of this method with the following ODE, we´ll find a general solution

http://quicklatex.com/cache3/de/ql_45d5634b7eb1776606cbae76454a5ede_l3.png

First we transform y’ into dy/dx, and then we separate both variables.

http://quicklatex.com/cache3/4f/ql_a50fd405cd63b776211965101661ae4f_l3.png

Now, we integrate both sides

http://quicklatex.com/cache3/ec/ql_6f0d565dffc12658e522ceb8d47c3bec_l3.png

Rearranging the formula, we can see that it belongs to a large array of ellipses, which vary with the constant “c”

https://steemitimages.com/DQmdEnVPDBpytDMWcJzP34L4uQ3ohvKx6WpxbqxdyvX7km9/Elipses.jpeg

http://quicklatex.com/cache3/5c/ql_b18ad4a76f6593afca8e58007ce3105c_l3.png

Given this powerful method for finding solutions to first-order ODE we can continue into the real deal, the applications of these equations.

Radioactive Decay

Experiments show that at each instant a radioactive substance decomposes—and is thus decaying in time—proportional to the amount of substance present.

If we chose to represent y(t) as the amount of substance present at any time “t” of a radioactive element. We know by physical law that the time rate change is represented by y’=dy/dt.

Giving the stated (by experiments) proportionality, we define the next first order ODE

http://quicklatex.com/cache3/48/ql_51626b29c15edd16a77687f500c02b48_l3.png

This is because in math, proportionality between two magnitudes is represented by a constant “k” giving it a minus in the formula, we´ll get the decaying implied. The value of "k" is found with experiments of the radioactive substance of interest.

Applying the method of separating variables we get easily the general solution of the ODE

http://quicklatex.com/cache3/d0/ql_3a73f7946a6c3fd8c73f9d8d7d27b7d0_l3.png

By properties of the natural logarithm and “e” constant

http://quicklatex.com/cache3/b4/ql_6ab3031aa5f9b967418eb4ac31e336b4_l3.png

Where “c” is the initial amount of the radioactive substance of interest, say “y₀” the other constant “k” should have a compatible unit with “t”

http://quicklatex.com/cache3/d2/ql_ff891d2c54860c4e82bd5728a88956d2_l3.png

Given 0,9 g of ²³²U, find the amount present at any later time, given that it´s constant "k" is 1x10^-2 years^-1

We start with the formula of Radioactive decay, and then we substitute our initial values to find the answer with relative ease

https://steemitimages.com/DQmPTE3UP4aZVueyousQjbZEj6bhgTj3zeHsYqKfpXdbp66/Descomposici%C3%B3n.jpeg

http://quicklatex.com/cache3/fa/ql_f4c3af23fddf4295d2271d307dba56fa_l3.png

This formula will give us any amount of ²³²U in grams when introducing a time in years given the initial amount of 0,9 g. For example, in a hundred years there will be

http://quicklatex.com/cache3/9f/ql_3dd79f2180f33cdab2418c8fe2253a9f_l3.png

Radiocarbon Dating

This is a specific case of study that derives from Radioactive decay, there is radioactive carbon ¹⁴C in the atmosphere, and it gets into the blood stream of any living being by breathing and eating, but when the organism dies, this absorption halts, effectively stopping the gathering of new ¹⁴C molecules in the body. The isotope ¹⁴C decays into ¹²C, so by measuring the rate of ¹⁴C/¹²C in a relic, we can determine the age of said remain.

For example, we´ll suppose there is a relic with a rate of ¹⁴C/¹²C = 58 %, and we need to know how old it is, for this we are given the half-life of ¹⁴C, which is 5715 years

Taking the equation that we found in Radioactive Decaying we start with

http://quicklatex.com/cache3/d2/ql_ff891d2c54860c4e82bd5728a88956d2_l3.png

Half-life is a measure of time, which says the amount of time needed to decompose a certain quantity of a radioactive element to the half, having the half-life (H = 5715) for ¹⁴C we can found “k” by simply taking “H” as “t” and “y” as y₀/2

http://quicklatex.com/cache3/19/ql_b7e6488e16e5f253df1ac9ad27ca9319_l3.png

By properties of the natural logarithm and e constant

http://quicklatex.com/cache3/86/ql_d6ce714e6a21696e863b07a29951f886_l3.png

Now, we use the ratio of carbon ¹⁴C - ¹²C = y/ y₀ = 58 %

http://quicklatex.com/cache3/8f/ql_bb629106397c49d3d56f66a3cd54298f_l3.png

Some more logarithm properties and we get the result in years

http://quicklatex.com/cache3/92/ql_672adb623dac96698c81f8ae1370cd92_l3.png

Mixing Problem

We can apply ODE of first grade to solve a mathematical model of a mixing tank. This is something very common in the chemical industry.

https://steemitimages.com/DQmTDuetgcSh97S6chDStdoJSS8ZRwBCPFpTxLeNyVCYPA8/image.png

For this example, we have a tank that contains 3800 litres of water, were 45 kg of salt were deluded initially (we won´t take into consideration the volume of salt). Brine gets into the tank through a pipe at an estimated rate of 38 litres/min, and each litre has 0,5 Kg of salt in it. Brine also runs out from another pipe at 38 litres/min. Find the amount of salt in the tank at any time "t"

Assuming y(t) as the amount of salt in the tank at a time “t”. The rate of change will be

http://quicklatex.com/cache3/06/ql_17739328ed24a04dc394f6f70ebb4906_l3.png

With 0,5 Kg of salt times 38 litres we have 19 Kg of salt entering the tank, while the outflow is 38 litres of brine, this represents 38/3800 = 0.01 (1 %) of the total volume of the tank, that is 0,01 of the salt in the tank at any time y(t), so the outflow can be written as 0,01y, and the ODE for this problem is

http://quicklatex.com/cache3/33/ql_69d73cb16c28a8398b17673120b7d833_l3.png

Let´s use the separating variable method now

http://quicklatex.com/cache3/12/ql_cf2ebaeed9b1b97bac976b8352987312_l3.png

As always, using properties of the natural logarithm and e constant

http://quicklatex.com/cache3/37/ql_ca9fea62ec50405c6a1321708812fc37_l3.png

Now we take our initial values, with t₀ there are 45 kg of salt in the thank, so

https://steemitimages.com/DQmaZrnj6THLkUVbMJ76AiYgA1v2zxK7Ddyh8ohAGsNDXuQ/Mezcla%20de%20tanque.jpeg

http://quicklatex.com/cache3/13/ql_b06a4a639c23113964906fc839a64c13_l3.png

http://quicklatex.com/cache3/8a/ql_07f95fccec079fc546bc3879f28cf48a_l3.png

This equation will grant us the amount of salt in the tank (in kg) in any given moment, by just using it with the time in minutes.

Newton’s law of cooling

Parting with the assumption that the time rate of change of the temperature of a body (T for the temperature) that conducts heat well, like copper, is proportional to the difference of the temperature of the body and the surrounding medium, we can go on with the next application of ODEs, referred as Newton´s law of cooling.

In our next example we have an office room with an air conditioner set to 16 °C, the office (and thus the AC also) works from 6 a.m. to 8 p.m. In summer the temperature around the office reached 30 °C (We will assume it stayed like this the whole day, because it´s simpler to model) and in a given day at 1 a.m. the temperature of the building was found to be 23°C, What was the temperature of the office at opening time (6 p.m.)?

Let T(t) be the temperature inside the building and 30 the outside temperature “T_A” (assumed to be constant in Newton’s law). Then by Newton’s law

http://quicklatex.com/cache3/0a/ql_d5674a4a1a835c6e190bdc2b69015c0a_l3.png

Applying the method of separating variables and logarithm properties

http://quicklatex.com/cache3/5f/ql_71d962904fdcd16cc96b1b03a0f4785f_l3.png

http://quicklatex.com/cache3/b3/ql_2af3ca8be8046cd1958b47556f1e0bb3_l3.png

Taking “t₀” as 6 p.m. we assume the temperature in the building at the end of day is the one the AC was set to function, so T(t₀) = 16 °C

http://quicklatex.com/cache3/ba/ql_3217467962821e04dcf869467a1b29ba_l3.png

With the other data “at 1 a.m. the temperature of the building was found to be 23°C” we can assume that T(t₇) = 23 °C

http://quicklatex.com/cache3/4a/ql_c9cfe6ee37a2430f7b5c79f5a44fd54a_l3.png

By properties of the natural logarithm and e constant

http://quicklatex.com/cache3/83/ql_808ed910888a91eb1d8d1e3c4185d683_l3.png

Finally with our fully defined solution we can find the value of T(t₁₀) (the temperature at 6 a.m.)

http://quicklatex.com/cache3/4e/ql_1fca820bc13325c8bb24f794a7e30e4e_l3.png

Leaking Tank (Torricelli’s Law)

Maybe we want to find the height of water in a tank at any time (t), Knowing that the tank has a 2 m diameter, a hole with a diameter of 1 cm, and the initial height of the water when the hole was drilled was 2,25 m. Of course this problem can be solved with our current understanding of ODEs, first of we apply the Torricelli’s Law, that law says that under the influence of gravity the velocity of outflowing water can be modeled as

https://steemitimages.com/DQmduh3NCButA9AhV54UNcxzdkeqsaxtWAUzYGxgmWXjpCU/image.png

http://quicklatex.com/cache3/7b/ql_74a7fa624960c4ad8ee4f07f5fcb797b_l3.png

Where “g” is the acceleration of gravity = 9,81 m/s² = 35316 m/min² for this example. Outflowing volume and its variation can be expressed as the Area of the hole (A) multiplied by the velocity “v” during a variation of time “t”

http://quicklatex.com/cache3/16/ql_d88b14ef38e74ec1c790d9b14e752a16_l3.png

Similarly, the variation of volume in the tank can be modeled as the Area of the tank multiplied by the change in height “h”

http://quicklatex.com/cache3/75/ql_6ece0ec339cbc711acd1fd20585c8f75_l3.png

Equating both formulas of ∆V

http://quicklatex.com/cache3/fe/ql_da44c566481975f166be6cbb691445fe_l3.png

Let´s make the interval of study ∆ approach to a value of 0 so we can get the ODE of this model; this technique is standard in ODE modeling to obtain derivatives.

http://quicklatex.com/cache3/50/ql_d374b92219adb0bdb2d8b4ceef733c50_l3.png

Substituting the Torricelli’s Law in the last expression

http://quicklatex.com/cache3/78/ql_fcfbff4289400014860ed6f444a83f78_l3.png

We can now begin to separate the variables and find the solution to the ODE, but first we should note that the square root of 2g times 0,6 equals to 159,46

http://quicklatex.com/cache3/aa/ql_1ca48e0a4d38b6083b30d5c66da5aaaa_l3.png

http://quicklatex.com/cache3/12/ql_82c96de6cea2d018d97417b12319d012_l3.png

For our problem we have that

http://quicklatex.com/cache3/e6/ql_c5f61a829fea6d67e76f6353b37f30e6_l3.png

Substituting in the former formula we get

http://quicklatex.com/cache3/01/ql_962b9d0a6c9122e24064b664bda59501_l3.png

With the initial value of h(t₀) = 2,25 m or 225 cm we find “c”

http://quicklatex.com/cache3/8b/ql_db8e044a1a564582af13e33ba82a0b8b_l3.png

And the solution is

http://quicklatex.com/cache3/30/ql_57044239c7f50eb806a4ede0e4211d30_l3.png

https://steemitimages.com/DQmeoSNr9GW7Ng19qjduS9pzNvXcr4UHjccrE1YoysmB6Qc/Vaciado.jpeg

About making this post

I wanted to make this post because I find the ODEs fascinating, as everything we use in our efforts to make reason out of the nature.I learned everything here thanks to the Tenth Edition of ''Advanced Engineering Mathematics'' by Erwin Kreyzig, which is an amazing book, the examples were picked, adapted and given another explanation with the help of QuickLatex for the formulas, AutoCAD for the drawings and the mobile version of Mathlab graphic calculator for the graphics. I truly hope you enjoyed this post as much as I did making it.