All Applications For 1st and Higher Order ODEs - Ch1
All Applications For 1st and Higher Order ODEs - Ch1
All Applications For 1st and Higher Order ODEs - Ch1
Al-Ne’aimi 1
dv dv k
m mg kv vg
dt dt m
Solution:
dv k
The velocity v satisfies the equation vg
dt m
where, g is the gravitational constant and k is the constant of proportionality.
Letting b k / m, you can separate variables to obtain
dv 1
dv ( g bv) dt g bv dt
b
ln g bv t C1
ln g bv bt bC1 g bv C e bt
Because the object was dropped, v 0 when t 0; thus g C , and it follows that
bt g ge bt mg
bv g g e v (1 e kt / m ).
b k
Chapter 1: Ordinary Differential Equations ODEs 2
Solution:
Since the object is falling, the resistance is in the upward (positive) direction. Hence,
mv mg kv 2 …………………….……… ……..(1)
where k is a constant. Since the magnitude of the resistance is 8 N when v 2 m/s,
k ( 2 )2 8
so k 2 N-s2/m2. Since m 10 and g 9.8, Eq.(1) becomes
10 v 98 2v 2 2( v 2 49 ) ……....…………..(2)
If vo 7 then v 7 for all t 0. If vo 7 separate the variables to obtain
1 1
v …………..…………………………..(3)
v 2 49 5
1 1
Using partial fraction expansion v
v 2 49 ( v 7 )( v 7 )
1 1 1
……..(4)
14 v 7 v 7
Substituting (4) into (3) yields
1 1 1 1
v
14 v 7 v 7 5
So
1 1
14
v 7 v 7 v
5
Integrating this yields
ln v 7 ln v 7 14t / 5 k
Therefore
v7
e k e14t / 5
v7
v7
C e14t / 5 …………..………………...………..(5)
v7
C e 14t / 5
v 7 ……..…….…………..………..(6)
C e 14t / 5
v0 7
C
v0 7
v0 ( 1 e 14t / 5 ) 7( 1 e 14t / 5 )
v 7
v0 ( 1 e 14t / 5 ) 7( 1 e 14t / 5 )
Consider a tank which initially holds V0 gal of brine that contains a Ib of salt.
Another solution, containing b Ib of salt per gallon, is poured into the tank at the
rate of e gal/min while simultaneously, the well-stirred solution leaves the tank
at the rate of f gal/min (see Figure below).
f
Or y y be
V0 et ft
Solution:
The data given in this example are V0 = 1000, a = 100, b = 5, and e = f = 10.
10
y y 5 10
1000 10t 10t
10 gal/min.
d
dt ( y e 0.01t ) 50 e 0.01t C
50 0.01t
y e 0.01t e C y 5000 C e 0.01t
0.01
Solution:
Let y = number of gallons of alcohol in the tank at any time t. You know that
y 5 when t 0 . Because the number of gallons of solution in the tank at any
time is 50 t , and the tank loses 5 gallons of solution per minute, it must lose
5 gallons of alcohol per minute.
y
50 t
Furthermore, because the tank is gaining 2 gallons of alcohol per minute, the rate
of change of alcohol in the tank is given by
dy 5 dy 5
2 y y 2 (which is linear D.E.)
dt 50 t dt 50 t
5
dt
e e 50 t
P ( t ) dt
Thus,
1
e 5 ln ( 50 t )
( 50 t )5
and, the general solution is
y 2 1
dt C
( 50 t )
5
( 50 t )5
2 ( 50 t ) 4
50 t
y C (50 t )5
2
5
20 50 t 50 t
@ y( 0 ) 5 : C y 20
505 2 50
Finally, when t 10, the amount of alcohol in the tank is
5
50 10 50 10
y 20 13.45 gal
2 50
which represents a solution containing 33.6% alcohol.
Volume of water leaving the tap exit = Volume of water LOSS in the tank
d 2
2 g h( t ) t D h( t )
2
4 4
h( t ) d2
h( t )1 / 2 2 2 g
t D
And since t 0 the 1st. order differential equation for the draining of a water
tank is expressed as:
dh( t ) d2
2 g 2 h( t )
t D
This Equation can be solved by separating the variables and integrating as:
dh( t ) d2
2 g 2 dt
h( t ) D
1 / 2 d2
h dh 2 g 2 dt C
D
d2
2h 1/ 2
2g 2 t C C 2 ho (from initial condition)
D
Thus, the general solution is given by:
2
g d 2
h( t ) t h0
2 D 2
To determine the time required to drain the tank at any time, te ?
2
g d2
Take h( te ) 0 : 0 te h0
D2
2
D2 2h0
te .... seconds
d2 g
Solution:
2
D2 2h0 1 21
te 35.89 seconds
d2 g 1 32.2
12
The outflowing water has velocity: 𝑣(𝑡) = 0.6 √2𝑔 ℎ(𝑡) (𝑻𝒐𝒓𝒓𝒊𝒄𝒆𝒍𝒍𝒊’𝒔 𝒍𝒂𝒘)
where h(t): water height above the hole at time t , and g: acceleration of gravity.
The contraction factor 0.6 is introduced because the stream has a smaller cross
section than the area of the hole.
Solution:
Decrease of water volume in the tank = Volume of water outflow in a time interval
−𝐴𝑡 𝛥ℎ = 𝐴ℎ 𝑣 𝛥𝑡
where, 𝐴𝑡 , 𝐴ℎ = Cross-sectional area of tank and the hole
Δh = Decrease of water height in the tank
Δt = a short time interval, Δt → 0.
𝛥ℎ 𝐴ℎ
=− 𝑣
𝛥𝑡 𝐴𝑡
𝑑ℎ 𝐴ℎ
= − 0.6 √2 (980) ℎ
𝑑𝑡 𝐴𝑡
𝑑ℎ 𝐴ℎ
= −26.56 √ℎ
𝑑𝑡 𝐴𝑡
𝑑ℎ 𝐴ℎ
∫ = −26.56 ∫ 𝑑𝑡
√ℎ 𝐴𝑡
𝐴ℎ
2 √ℎ = −26.56 𝑡 +𝑐 𝐷𝑖𝑣𝑖𝑑𝑖𝑛𝑔 𝑏𝑦 2 𝑎𝑛𝑑 𝑠𝑞𝑢𝑎𝑟𝑖𝑛𝑔
𝐴𝑡
2
𝐴ℎ 2 𝜋 0.52
ℎ = (−13.28 𝑡 + 𝑐) = (−13.28 2 𝑡 + 𝑐) = (−0.000332 𝑡 + 𝑐)2
𝐴𝑡 𝜋 100
Tank empty: At h = 0 , t =?
0 = (−0.000332 𝑡 + 15)2
15
𝑡= = 45181 𝑠𝑒𝑐 = 12.6 ℎ𝑟
0.000332
Solution:
The 1st. order differential equation for heat transfer in solids in the 𝒙-direction is
derived from the Fourier Law of heat conduction as:
dT ( x )
Q q A k A
dx
where, Q total amount of heat flow,
q heat flux - a sense of the intensity of heat flow Q in a solid,
A the area to which heat flows,
k thermal conductivity constant,
dT ( x ) difference of temperatures between left and right surfaces, and
x = distance or thickness of heat flow
The total heat flow Q per unit time t ( Q / t ) in the rod is given by the heat source
to the left end, i.e. 10 kW [From heat flux ( Q / t ) q A 10 kW].
dT ( x ) Q 10
Or
83.33 0C / m
dx kA 100( 1200 10 )6
Solution:
If te required time for the solid to drop its temperature from 80oC to
8oC, we should have: T ( te ) 8 5 75 e 0.0004 te
Solution:
Setting T ( 0 ) T0 yields C T0 Tm , so
T Tm (T0 Tm ) e kt ………..………..(2)
Now, from data given in this example: T0 = 400 and Tm = 25, Eq.(2) becomes
T 25 375 e kt ………..………………..(3)
7
From the stated condition T ( 4 ) 200 200 25 375 e 4k e 4k
15
Taking logarithms and solving for k yields
1 7 1 15
k ln ln
4 15 4 7
Substituting this into (3) yields: T 25 375 e ( t / 4 ) ln 15 / 7
and therefore the temperature of the insulator after 8 minutes is
2
2 ln 15 / 7 7
T ( 8 ) 25 375 e 25 375 107 o C
15
Mechanical Vibrations
Consider the motion of an object with mass m
at the end of a spring that is either vertical (as in
Fig.1) or horizontal on a level surface (as in
Fig.2).
d2x
m Fgrav Fspr Fdamp F (t ) ……………..……..(1)
2
dt
where,
Fgrav is the downward weight of the mass, equal to: Fgrav mg
where, m mass of the object, and g acceleration = 9.8 m/sec2 in MKS
units, g 980 cm/sec2 in CGS units, and g 32 ft/sec2 in British units.
Fspr is modeled by Hooke’s law as: Fspr k ( x xo )
where, k is a positive constant called the spring constant, xo and x the
displacements of the spring without and with the mass attached.
Fdamp c dx , where c 0 is the so-called damping coefficient.
dt
d2x dx
m c kx ( kxo mg) F ( t )
dt 2 dt
Neglecting the term ( kxo mg ) which belongs to the mass rest position from
Eq.(1), the motion of the spring-mass-system is governed by a 2nd. order
nonhomogeneous linear differential equation written as:
d2x dx
m 2
c kx F ( t ) …………………………………..(2)
dt dt
The solution of Eq.(2) depends on whether the vibrating system is free or forced
and damped or undamped, in addition to degrees-of-freedom, and initial
conditions.
The system is free if it is unforced; that is, F ( t ) 0 ; otherwise it is forced. The
system is undamped if c 0 (equivalently 0 ); otherwise, it is damped.
Therefore, four possible cases can be considered as follows:
d2x
m kx 0 n2 x 0 ……………………..(3)
x
2
dt
where, A > 0 is its amplitude and in 0, 2π is its phase angle. Upon comparing
Eq.(3.3) with Eq.(3.2) we see that
c2
A c12 c22 and tan1
c1
Solution:
(ii) The amplitude of the motion A is 5/2 cm, the phase angle is 3 / 2 , the
natural frequency n is 14 1/sec, and the period T 2 / n is / 7 sec.
(iii) The +ve. time at which the mass first returns to its rest position is t / 14.
(b) Forced Udamped Vibrations [In this case c 0 and under harmonic loading
F ( t ) F cos( t ) ] Eq.(2) reduces to:
d2x F
m 2
kx F cos(t ) x
n2 x cos(t ) .……….…..(4)
dt m
F
A , and B0
m(n2 2 )
F ……………………….……….…..(4.2)
Hence x p (t) cos(t )
2
m ( n 2 )
For n :
For n 0 :
Two very important phenomena occur when the driving frequency becomes close
to the system’s natural frequency: resonance and beating.
In case of resonance, the particular solution for Eq.(4) using the modification rule
becomes:
x p ( t ) t ( A cos n t B sin n t )
By substituting this into the Eq.(4) and equating coefficients; gives A 0 and
B F /( 2 mn ) . Hence
F ……………………………….……….…..(4.6)
x p (t ) t sin(nt )
2 mn
We see that because of the factor t, the amplitude becomes larger and larger. So
practically speaking, systems with very little damping may undergo larger
vibrations that can destroy the system (see Fig.4 shown below).
x p (t)
F
2
cos(t ) cos(n t) ( n ).
m ( n 2 )
2F ….……….…..(4.7)
x p (t) sin n t sin n t
2
m ( n 2 ) 2 2
but negligible damping. Any horizontal motion of the system should be ignored.
In normal operation the machine is subjected to a vertical force F ( t ) Fo sin t
where the amplitude Fo is 2500 N. Calculate the response amplitude and the
force transmitted to the foundations when the driving frequency is (a) 20 Hz and
(b) 2 Hz.
Solution:
Fo 2500 2500
Po 1116 N
1 r2 1 1.8 2 2.24
d2x dx
m c kx 0 2 n x n2 x 0
x …..……..…..(5)
2 dt
dt
knowing that:
the natural frequency: n k / m 0
the damping ratio: c / 2 km 0
the damped natural frequency: d n 1 2 (only relevant for 0 1) .
c k
r2 r 0 OR r 2 2 nr n
2
0
m m
c c 2 4 mk
with roots, r1 OR r1 n n 2 1
2 2m 2
The nature of the solution depends upon whether is less than one, equal to one,
or greater than one.
CASE I: c 2 4 mk 0 OR ( 1) (Overdamping)
In this case, r1 and r2 are distinct real roots and the solution is given by:
x (t) c1 e r1t c 2 e r2 t
…..…..(5.1)
( n n 2 1 )t ( n n 2 1 )t
OR x(t) c1 e c2 e
In this case, oscillations do not occur due to strong damping force (high-viscosity
oil or grease) compared with a weak spring or small mass.
x (t) (c1 c 2 t ) e ( c / 2 m ) t
…………………….…………....(5.2)
n t
OR x (t) ( c1 c 2 t ) e
In this case, the damping is just sufficient to suppress vibrations, and any decrease
in the viscosity of the fluid leads to vibrations of the following case III.
CASE III: c 2 4 mk 0 OR ( 1) (Underdamping)
Here the roots are complex:
c 4 mk c 2
r1 i where
2 2m 2m
Solution:
x (t) c1 e 4 t c 2 e 16 t
x(t) 0.05 (e 4 t e 16 t )
F (t )
cx kx F cos(t ) x
mx 2 n x n2 x ………..(6)
m
All possible cases for the homogeneous solution xh (t ) are already given before
in Section (c) for FREE DAMPED VIBRATIONS.
Damping ratio: c c c c
, 2 n
cc 2 mn 2 km m
frequencies ratio:
r
n
2
m ( n 2 ) ( 2 2 ) 2 4 2 2 2 F
AF n n
A
m 2 ( n
2
2 ) 2 2 ( 2 m n ) 2
2
( n 2 ) m
(F / k )
X ……………………..……….…..(6.3)
2 2 2
( 1 r ) ( 2 r )
OR X st M
F
where, st deflection under the static force F and M magnification factor.
k
(F / k )
OR X A2 B 2
( 1 r 2 ) 2 ( 2 r ) 2
and the phase angle is:
1 2 mn 1 2 n r
2
φ tan tan
m 2 m 2 2 2
n n
2r
φ tan1
2
1 r
B 2r
OR φ tan1 tan1 …………………..……….…..(6.4)
A 2
1 r
Solution:
cx kx F cos ( t )
mx
Hence
d
dr
( 1 r 2 )2 ( 2r )2 0
d
Or ( 1 r 4 2 r 2 4 2 r 2 ) 4 r 3 4 r 8 2 r 0
dr
Or r 0 Or r 2 1 2 2 Or r 1 2 2
For r 1 2 2
(F / k)
X max
( 1 1 2 2 )2 4 2 ( 1 2 2 )
(F / k) (F / k)
4 4 4 2 8 4 4 2 4 4
(F / k)
2 1 2
(F / k) 100
So given X max 5 10 3 5 10 3 ..…(1)
2 1 2
k 2 1 2
10 10000 2
r 1 2 2 1 2 2 …..…(2)
n k / 100 k
100 10000 2
Or k
3 1 2 2
5 10 2 1 2
1 4
Or
2 ( 1 2 ) ( 1 2 2 )2
1 4 4 4 2 4 ( 2 4 ) 0
( 4 4 ) 4 ( 4 4 ) 2 1 0
Or
101.4091 ( 101.4091)2 4( 101.4091)( 1 )
2
9.9603 10 3 or 0.0998
2 101.4091
100.7 10 3 N / m 100.7 kN / m
100.7 103 c
c 2 0.0998 100 633.396 N.s/m where 2 n
100 m
Solution:
Buoyancy force, FB = – Weight of displaced water, W
FB = Mass × Acceleration = 𝑚𝑦”
W = Volume of displaced water × Specific weight of water = – 𝜋 𝑟 2 𝑦 𝛾
𝒎𝒚” = – 𝝅 𝒓𝟐 𝒚 𝜸
where y = displacement of cylinder from its
equilibrium position.
√𝛑 𝒓𝟐 𝜸/𝒎 = 𝝎𝟎 hence 𝒚” + 𝝎𝟎 𝟐 𝒚 = 𝟎
𝟐𝛑
cos and sin have the period = = 2 sec 𝜔0 = π
𝝎𝟎
𝛾
𝜔0 2 = π 𝑟 2
𝑚
2 2
9800 (N/m3 )
𝑚 = π 0.3 (m ) 2 = 281 kg
𝜋 (N/kg m)
𝑊 = 𝑚 𝑔 = 281 × 9.81 = 2754 N.
Solution:
𝒚” + 𝝎𝟎 𝟐 𝒚 = 𝟎
Weight of vibrating water = 𝜋 𝑟 2 2𝑦 𝛾
𝛾
𝜔0 2 = 2π 𝑟 2
𝑚
9800
𝜔0 2 = 2π 0.012 = 2.5
1
1 𝜔
Frequency = = 0 = 0.4 sec −1
𝑝𝑒𝑟𝑖𝑜𝑑 2π
Solution:
2 2
𝑻𝒂𝒏𝒌 𝟏: 𝑦1′ = 𝑦2 − 𝑦
100 100 1
𝐷𝑦1 + 0.02 𝑦1 − 0.02 𝑦2 = 0
(𝐷 + 0.02) 𝑦1 − 0.02 𝑦2 = 0 … … … . . (1)
2 2
𝑻𝒂𝒏𝒌 𝟐: 𝑦2′ = 𝑦1 − 𝑦
100 100 2
𝐷𝑦2 − 0.02 𝑦1 + 0.02 𝑦2 = 0
−0.02 𝑦1 + (𝐷 + 0.02) 𝑦2 = 0 … … … . (2)