Applications of First Order DE
Applications of First Order DE
A. Population Growth/Decay
Population growth/decay follows the idea that the rate of growth/decay of a certain population is proportional to
the current population, say P. Putting it into equation,
dP
± αP
dt
Equation [1] will be our working equation in solving population growth/decay problems.
Positive sign will denote a growth, negative sign will denote a decay. Also, k has a unit of per unit time (1/sec,
1/year, etc.).
Examples:
1. A current state has a population of 4200. With the rate of 0.1/year, what will be the population after 3 years?
Note: You answers f or populations must be in whole numbers, because there is no person counted as
0.41. You need to round it up to the nearest whole number (i.e. 2.4 will be rounded up to 3).
2. A certain population doubles every 5 years. How long would it take f or the population to be 5 times the
initial population?
3. A town with initial population of 6245 developed a mysterious disease that resulted in deaths of 0.2/year.
The mayor started the head-count and f ound out that there are only 3600 people lef t. How long did the
deaths persisted before the mayor decided to do head count?
4. The current population of a certain state is 32,602. Over the past decade, there have been deaths of
0.23/year due to risks of obesity. Determine the population before the said deaths began.
5. A strain of Staphylococcus aureus doubles every 12 hours. A colony of 200 was cultured in a petri dish for
12 hours. What is the population of the strain after 30 hours?
B. Radioactive Decay
Radioactive decay f ollows the idea that the rate disappearance of a certain substance is proportional to the
current amount, say C. Putting it into equation,
dC
− αC
dt
Equation [1] will be our working equation in solving radioactive decay problems.
k still has a unit of per unit time (1/sec, 1/year, etc.).
• Half-Life
Half -life is the time taken to reduce the amount of a certain substance to half of its original value.
In terms of equation:
1
@ t = t 1 ; C = Co
2 2
Where t is any given time, C is the final amount or concentration and Co is the initial amount/concentration.
Plugging-in the condition:
1
@ t = t 1 ; C = Co
2 2
1 −kt 1
Co = Co e 2
2
1 −kt 1
=e 2
2
Taking the natural logarithm (ln) for both sides:
1
ln ( ) = −kt 1
2 2
1
ln ( ) = − ln2
2
− ln 2 = −kt 1
2
ln2 = kt 1 [𝟐]
2
Equation [2] can be used to solve the rate constant using half-life.
There are cases in which initial and f inal amounts are not given, but rather expressed in terms of %
conversion (conversion is the fraction of reactant that proceeds into product), or in decimal form XA. Putting
it into equation:
remained = initial − converted
C = Co − Co XA
C = Co (1 − XA ) [𝟑]
C = Co e−kt
Co (1 − XA ) = Co e−kt
(1 − XA ) = e−kt
ln(1 − XA ) = −kt [ 𝟒]
Examples:
1. A certain compound has a half-life of 23 days. What percent of the original amount will remain if 4 half-
lives had passed?
Solution:
# of half-lives # of days % remained
1 23 50
2 46 25
3 69 12.5
4 92 6.25
Based on the table above, after 4 half-lives 6.25% of the original amount remains.
Another solution is to use this formula if only the number of half-lives (n) is given:
C 1 n
=( ) [ 𝟓]
Co 2
Substituting n=4:
C 1 4
= ( ) = 𝟎. 𝟎𝟔𝟐𝟓 (𝟔.𝟐𝟓%)
Co 2
2. 30 g of an oxide of sulfur has a half-life of 2.3 years. How much will remain after a decade (t=10
years)?
We will use equation [1] to solve for the final amount, but we will use equation [2] f irst to solve for the rate
constant k.
ln2 = kt 1
2
ln2 = k(2.3 years )
0.3014
k=
yr
C = Co e−kt
0.3014
− (10 yrs)
C = 30e yr Note: k has a unit of per 𝐮𝐧𝐢𝐭 𝐭𝐢𝐦𝐞, it must cancel with the unit of your t.
𝐂 = 𝟏. 𝟒𝟕𝟑 𝐠 → 𝐅𝐢𝐧𝐚𝐥 𝐚𝐧𝐬𝐰𝐞𝐫
3. An isotope of carbon (C-14) has a half-life of 5730 years. A sample of wood was analyzed to contain
only 0.08 % of its carbon content. Estimate the age of the wood.
4. 50 g of an organochlorine compound has a half-life of 100 years. How long would it take to reduce the
amount to 30 g.
5. 50 g of an organochlorine compound has a half-life of 100 years. How long would it take to reduce the
amount by 30 g.
Newton’s Law of Cooling/Heating follows the idea that the rate of change in temperature is proportional to the
dif ference of temperature T and the medium temperature Tm. Putting it into equation,
dT
α (T − Tm)
dt
Examples:
1. An object with temperature of 80°C is placed in a room 30°C. After 2 hours, the object cools to 50°C. What
will be its temperature after 6 hours?
3. An object initially at 23°C is placed in a water bath at 100°C f or sterilization purposes. After 30 minutes, it
reached 73 °C. How long would it take to reach 90°C?
4. A metal at 86°C is placed in an oil bath at 200°C. After 30 minutes, the temperature reached 120°C. What
will be the temperature after 40 minutes?
T = 200 − 114e−0.0118(40)
T = 𝟏𝟐𝟖. 𝟗𝟏°𝐂 → 𝐅𝐢𝐧𝐚𝐥 𝐚𝐧𝐬𝐰𝐞𝐫
5. At 10:00 PM, a body was found with gunshot wounds, in an apartment with room temperature of 30°C. The
criminologist immediately checked the body temperature and was found out to be 33°C. After 30 minutes,
the criminologist checked the body temperature again and it was f ound out to be 31°C. Considering that
the temperature of the human body is 37°C, what time did the victim died?
D. Mixtures
Consider the diagram above, Qin (unit volume per unit time) is the entering volumetric f low rate with
concentration Cin (unit mass per unit volume). This entering stream is poured in a container with volume V and
initial salt mass x o. On the other hand, Qout is the exiting flow rate with concentration Cout. The assumption is
that the solution inside the container is well mixed, that the concentration at any given point is same although
out. The mass of salt after mixing and leaving the container is labeled x.
dx
= rin − rout [𝟏]
dt
Where rin is the mass flow rate of salt entering, obtained from
rin = Qin Cin
The solution will differ from with the flowrates entering. Here are the cases:
For this case, equation [1] simplifies into a variable-separable DE. However, it can be solved also using
linear DE.
For this case, there is accumulation since entering f low rate is greater than exiting flow rate. The f ormula
f or Cout changes due to accumulation.
rout = Qout Cout
x
rout = Qout ( )
Vt
Where:
Vt = V + (Qin − Qout )t
For this case, there is no more accumulation. But since entering flow rate is now less than exiting flow
rate, the container is emptying. The f ormula for Cout is the same with Case 2 but expect the Q in – Qout
term to become negative.
rout = Qout Cout
x
rout = Qout ( )
Vt
Where:
Vt = V + (Qin − Qout )t
Examples:
1. A tank has salt solution flowing into it at 10 L/min with salt concentration 0.2 kg/L. The contents of the tank
are kept thoroughly mixed, and the contents flow out at 10 L/min. Initially, the tank contains 10 kg of salt in
100 L of water. How much salt will there be in the tank after 30 minutes? What is the maximum amount of
salt in the tank?
Solving rin:
L kg kg
rin = 10 x 0.2 = 2
min L min
Solving rout:
L x kg x kg
rout = 10 ( )=
min 100L 10 min
dx x
+ =2
dt 10
t t
xe10 = ∫ 2e10 dt + C
t t
xe10 = 20e10 + C
Now, we need to solve f or the constant of integration. Note that it is stated in the problem that the tank
initially contains 10 kg salt. Thus, putting it into the equation:
t t
xe10 = 20e10 + C, x(0) = 10
10 = 20 + C, C = −10
t t t
xe10 = 20e10 − 10; dividing both sides by e10
t
x = 20 − 10e− 10
Finally, solving for amount of salt in the tank after 30 minutes:
30
x = 20 − 10e− 10
𝐱 = 𝟏𝟗.𝟓𝟎 𝐤𝐠 → 𝐅𝐢𝐧𝐚𝐥 𝐚𝐧𝐬𝐰𝐞𝐫
2. A tank has pure water f lowing into it at 10 L/min. The contents of the tank are kept thoroughly mixed, and
the contents flow out at 10 L/min. Initially, the tank contains 10 kg of salt in 100 L of water. How much salt
will there be in the tank af ter 30 minutes?
Solving rin:
L kg kg
rin = 10 x0 =0 , since it says pure water
min L min
Solving rout:
L x kg x kg
rout = 10 ( )=
min 100L 10 min
dx dt
=−
x 10
dx dt
∫ =∫−
x 10
t
lnx = − + C
10
t t
x = e 10 = Ce− 10 ,x(0) = 10
− +C
C = 10
t
x = 10e− 10
at t = 30 mins.,
3. A tank with a 500 L capacity has salt solution f lowing into it at 18 L/min with salt concentration 0.15 kg/L.
The contents of the tank are kept thoroughly mixed, and the contents flow out at 8 L/min. Initially, the tank
contains 26 kg of salt in 150 L of water.
a. How much salt will there be in the tank after 25 minutes?
Using equation [1]:
dx
= rin − rout
dt
Solving rin:
L kg kg
rin = 18 x 0.15 = 2.7
min L min
Solving rout:
Accumulation = 150 + (18 − 8)t = 150 + 10t
L x kg 4x kg
rout = 8 ( )=
min (150 + 10t)L 75 + 5t min
dx 4x
+ = 2.7
dt 75 + 5t
4 4 ( ) 4
µ = e∫75+5t dt = e 5 ln 75+5t = (75 + 5t) 5
4 4
x(75 + 5t) 5 = ∫ 2.7(75 + 5t) 5 dt + C
4 9
x(75 + 5t) 5 = 1.5(75 + 5t) 5 + C
4
x = 1.5(75 + 5t) + C(75 + 5t) −5
Now, we need to solve for the constant of integration. Note that it is stated in the problem that the tank
initially contains 26 kg salt. Thus, putting it into the equation:
4
x = 1.5(75 + 5t) + C(75 + 5t) −5 ,x(0) = 26
4
26 = 1.5(75) + C(75) − 5, C = 750.8 (26 − 1.5(75))
4
x = 1.5(75 + 5t) + [750.8 (26− 1.5(75))] (75 + 5t) − 5
4. A tank with a 500 L capacity has pure water f lowing into it at 18 L/min. The contents of the tank are kept
thoroughly mixed, and the contents flow out at 8 L/min. Initially, the tank contains 26 kg of salt in 150 L of
water.
a. How much salt will there be in the tank after 25 minutes?
Using equation [1]:
dx
= rin − rout
dt
Solving rin:
L kg kg
rin = 18 x0 =0
min L min
Solving rout:
Accumulation = 150 + (18 − 8)t = 150 + 10t
L x kg 4x kg
rout = 8 ( )=
min (150 + 10t)L 75 + 5t min
dx 4x
=−
dt 75 + 5t
dx 4 dt
=−
x 75 + 5t
dx 4 dt
∫ = −∫
x 75 + 5t
4
lnx = − ln(75 + 5t) + C
5
4
x = C(75 + 5t) −5 ;x(0) = 26
4 4
26 = C(75) −5 ;C = 26(75) 5
4 4
x = 26(75) 5(75 + 5t) −5
At t = 25 mins,
4
4 −
5
x= 26(75) 5(75 + 5(25))
𝐱 = 𝟏𝟏. 𝟖𝟔 𝐤𝐠 → 𝐅𝐢𝐧𝐚𝐥 𝐚𝐧𝐬𝐰𝐞𝐫
5. A tank with a 1000 L capacity has salt solution flowing into it at 23 L/min with salt concentration 0.17 kg/L.
The contents of the tank are kept thoroughly mixed, and the contents flow out at 35 L/min. Initially, the tank
contains 46 kg of salt in 350 L of water.
a. How much salt will there be in the tank after 27 minutes?
Using equation [1]:
dx
= rin − rout
dt
Solving rin:
L kg kg
rin = 23 x 0.17 = 3.91
min L min
Solving rout:
Accumulation = 350 + (23 − 35)t = 350 − 12t
L x kg 30x kg
rout = 30 ( )=
min (350 − 12t) L 350 − 12t min
dx 30x
− = 3.91
dt 350 − 12t
30 5 (350−12t) 5
µ = e∫− 350−12t dt = e− 2 ln = (350 − 12t) −2
5 5
x(350 − 12t) −2 = ∫ 3.91(35 − 12t) −2 dt + C
5 3.91 3
x(350 − 12t) −2 = (35 − 12t) − 2 + C
36
3.91 5
x= (350 − 12t) + C(350 − 12t) 2 ;x(0) = 46
36
3.91 5
46 = (350) + C(350) 2
36
36(46) − 3.91(350)
C= 5
36(350) 2
3.91 36(46) − 3.91(350) 5
x= (350 − 12t) + [ 5
] (350 − 12t) 2; t = 27 mins
36
36(350) 2
3.91 36(46) − 3.91(350) 5
x= (350 − 12(27)) + [ 5
] (350 − 12 ( 27)) 2
36
36(350) 2
𝐱 = 𝟐. 𝟖𝟒 𝐤𝐠 → 𝐅𝐢𝐧𝐚𝐥 𝐚𝐧𝐬𝐰𝐞𝐫
6. A tank with a 1000 L capacity has pure water flowing into it at 23 L/min. The contents of the tank are kept
thoroughly mixed, and the contents flow out at 35 L/min. Initially, the tank contains 46 kg of salt in 350 L of
water.
a. How much salt will there be in the tank after 27 minutes?
Using equation [1]:
dx
= rin − rout
dt
Solving rin:
L kg kg
rin = 23 x0 =0
min L min
Solving rout:
Accumulation = 350 + (23 − 35)t = 350 − 12t
L x kg 30x kg
rout = 30 ( )=
min (350 − 12t) L 350 − 12t min
dx 30
=− dt
x 350 − 12t
dx −30
∫ =∫ dt
x 350 − 12t
5
lnx = ln(350 − 12t) + C
2
5
x = C(350 − 12t) 2 ;x(0) = 46
5 5
46 = C(350) 2; C = 46(350) −2
5 5
x = (46(350) −2 ) (350 − 12t) 2 ;t = 27 mins
5 5
x = (46(350) −2 )(350 − 12(27) )2
𝐱 = 𝟎. 𝟎𝟕 𝐤𝐠 → 𝐅𝐢𝐧𝐚𝐥 𝐚𝐧𝐬𝐰𝐞𝐫
c. When will the salt concentration in the tank reaches 0.02 kg?
5 5
x = (46(350) −2 ) (350 − 12t) 2 ;x = 0.02 kg
5 5
0.02 = (46(350) − 2) (350 − 12t) 2
𝐭 = 𝟐𝟕. 𝟖𝟓 𝐦𝐢𝐧𝐬 → 𝐅𝐢𝐧𝐚𝐥 𝐚𝐧𝐬𝐰𝐞𝐫
APPLICATIONS OF FIRST-ORDER DIFFERENTIAL EQUATION
E. Orthogonal Trajectories
Consider a f amily of curves f(x,y) = C, where C is constant. For this family of curves, we can generate another
f amily of curves g(x,y) = C (or its orthogonal trajectories) in which it crosses the given curves at right angles.
Examples:
The slope of the curve is equal to -y/x. From the def inition,
orthogonal trajectories cross the curves at right angles. Thus, we need to get the negative reciprocal of the
slope and equate it to dy/dx again.
dy x
( ) =−
dx old y
dy y
( ) =
dx new x
𝐱 𝟐 + 𝐲 𝟐 = 𝐂𝐱 → 𝐅𝐢𝐧𝐚𝐥 𝐚𝐧𝐬𝐰𝐞𝐫
dy y
( ) =
dx old x + 1
dy x+1
( ) =−
dx new y
dy x+1
=−
dx y
y dy = −(x + 1) dx
∫ y dy = ∫ −(x + 1) dx
y2 (x + 1) 2
=− +C
2 2
2 2
𝑦 = −(x + 1) + C
(x + 1) 2 + y 2 = C
𝐱 𝟐 + 𝟐𝐱 + 𝐲 𝟐 = 𝐂 → 𝐅𝐢𝐧𝐚𝐥 𝐚𝐧𝐬𝐰𝐞𝐫
RL Circuits
R V R
IeL t = e L t + C; I(0) = 0
R
V
C =−
R
R
t V Rt V
IeL = eL −
R R
V V − Rt
I= − e L
R R
V R
I = (1 − e− L t ) [ 𝟐]
R
Equation [2] can be used to solve the current I at any given time t.
Examples:
Find:
a. the equation for I
Using equation [2]
V R
I = (1 − e− L t )
R
70 30
I = (1 − e− 12t )
30
𝟕 𝟓
𝐈 = (𝟏 − 𝐞− 𝟐𝐭 ) → 𝐅𝐢𝐧𝐚𝐥 𝐚𝐧𝐬𝐰𝐞𝐫
𝟑
dI
VL = L
dt
d 7 5
VL = L ( (1 − e− 2t ))
dt 3
35 5
VL = L ( e− 2t )
6
35 − 5t
VL = ( Le 2 )
6
35 5
VL = ( (12)e− 2t )
6
𝟓
𝐕𝐋 = 𝟕𝟎𝐞−𝟐𝐭 → 𝐅𝐢𝐧𝐚𝐥 𝐚𝐧𝐬𝐰𝐞𝐫
5 1
lne− 2t = ln
2
5
− t = −ln2
2
5
t = ln2
2
2 ln2
t=
5
t = 𝟎. 𝟐𝟕𝟕 𝐬 → 𝐅𝐢𝐧𝐚𝐥 𝐚𝐧𝐬𝐰𝐞𝐫
Find:
a. the equation for I
Using equation [2]
V R
I = (1 − e− L t )
R
43 28
I = (1 − e− 19t )
28
𝟒𝟑 𝟐𝟖
𝐈 = (𝟏 − 𝐞− 𝟏𝟗𝐭 ) → 𝐅𝐢𝐧𝐚𝐥 𝐚𝐧𝐬𝐰𝐞𝐫
𝟐𝟖
dI
VL = L
dt
d 43 28
VL = L ( (1 − e− 19t ))
dt 28
43 28
VL = L ( e− 19t )
19
43 28
VL = ( (19)e− 19t )
19
𝟐𝟖
𝐕𝐋 = 𝟒𝟑𝐞−𝟏𝟗𝐭 → 𝐅𝐢𝐧𝐚𝐥 𝐚𝐧𝐬𝐰𝐞𝐫
28
− t = −ln2
19
28
t = ln2
19
19 ln2
t=
28
t = 𝟎. 𝟒𝟕 𝐬 → 𝐅𝐢𝐧𝐚𝐥 𝐚𝐧𝐬𝐰𝐞𝐫
d. Time constant
Time is constant (τ, or Greek letter tau) is defined as the time at which the system’s response already
reached 63.2% (equivalent to 1/e) of its final value. In the case of RL circuits, τ is equivalent to L/R.
L
τ= [ 𝟑]
R
19
τ= = 𝟎. 𝟔𝟕𝟗 𝐬 → 𝐅𝐢𝐧𝐚𝐥 𝐚𝐧𝐬𝐰𝐞𝐫
28
For this problem, transient current is represented by the equation for I itself, which is:
𝟒𝟑 𝟐𝟖
𝐈 = (𝟏 − 𝐞− 𝟏𝟗𝐭 ) 𝐀 → 𝐅𝐢𝐧𝐚𝐥 𝐚𝐧𝐬𝐰𝐞𝐫
𝟐𝟖
Find:
a. the equation for I
Using equation [2]
V R
I = (1 − e− L t )
R
53 21
I = (1 − e− 16t )
21
𝟓𝟑 𝟐𝟏
𝐈 = (𝟏 − 𝐞− 𝟏𝟔𝐭 ) → 𝐅𝐢𝐧𝐚𝐥 𝐚𝐧𝐬𝐰𝐞𝐫
𝟐𝟏
dI
VL = L
dt
d 53 21
VL = L ( (1 − e− 16t ))
dt 21
53 21
VL = L ( e− 16t )
21
𝟐𝟏
𝐕𝐋 = 𝟓𝟑𝐞−𝟏𝟔𝐭 → 𝐅𝐢𝐧𝐚𝐥 𝐚𝐧𝐬𝐰𝐞𝐫
d. Time constant
Using equation [3]
L
τ=
R
16
τ= = 𝟎. 𝟕𝟔𝟐 𝐬 → 𝐅𝐢𝐧𝐚𝐥 𝐚𝐧𝐬𝐰𝐞𝐫
21
Steady-state current.
53 21
I = (1 − e− 16t )
21
53
I=
21
𝐈 = 𝟐. 𝟓𝟐𝟒 𝐀 → 𝐅𝐢𝐧𝐚𝐥 𝐚𝐧𝐬𝐰𝐞𝐫
RC Circuits
Kirchhoff's voltage law says the total voltages must be zero. Thus, applying this law to a series RC circuit
results in the equation:
1
IR + ∫ I dt − V = 0
C
1
IR + ∫ I dt = V
C
dI 1
R + I = 0 [𝟒]
dt C
Equation [4] can be solved either using separation of variables or linear DE. In this case, we’ll use
separation of variables:
dI 1
R + I=0
dt C
dI 1
R =− I
dt C
dI dt
=−
I RC
dI dt
∫ = ∫−
I RC
t
ln I = − + C1
RC
t
I = e− RC+C1
t V
I = C2 e− RC;I(0) =
R
V
= C2
R
V −t
I = e RC [𝟓]
R
Example:
1. A series RC circuit with R = 8 Ω and C = 0.02 F is connected with a battery of E = 100 V. At t = 0, the
voltage across the capacitor is zero.
Find:
a. the equation for I
Using equation [5]
V t
I = e− RC
R
100 − t
I= (e (0.02)(8))
8
𝐈 = 𝟏𝟐.𝟓(𝐞−𝟔.𝟐𝟓𝐭 ) → 𝐅𝐢𝐧𝐚𝐥 𝐚𝐧𝐬𝐰𝐞𝐫
1
VC = ∫ I dt
C
1
VC = ∫ 12.5(e−6.25t ) 𝑑t
0.02
VC = −100e−6.25t + C1
At t = 0, V = 0
VC = −100e−6.25t + C1 ; V(0) = 0
𝐶1 = 100
VC = 100 − 100e−6.25t
𝐕𝐂 = 𝟏𝟎𝟎− 𝟏𝟎𝟎𝐞−𝟔.𝟐𝟓𝐭 → 𝐅𝐢𝐧𝐚𝐥 𝐚𝐧𝐬𝐰𝐞𝐫
1
e−6.25t =
2
1
lne −6.25t = ln
2
−6.25t = −ln2
6.25t = ln2
ln2
t=
6.25
t = 𝟎. 𝟏𝟏𝟏 𝐬 → 𝐅𝐢𝐧𝐚𝐥 𝐚𝐧𝐬𝐰𝐞𝐫
Example:
1. A series RC circuit with R = 31 Ω and C = 0.0019 F is connected with a battery of E = 100 cos(60t) V. At
t = 0, the charge in the capacitor is -0.039 C.
Find:
a. the equation for q