Scribd
Upload
351 views

Holistic Numerical Methods

This document provides an example of using the 4th order Runge-Kutta method to solve an ordinary differential equation modeling the concentration of salt in a tank over time. The method is applied with different step sizes to calculate the salt concentration at 3 minutes. Tables and figures show that decreasing the step size improves the accuracy of the numerical solution compared to the exact solution. The document also compares the 4th order Runge-Kutta method to lower order Runge-Kutta methods.

Uploaded by

Ohol Rohan Bhaskar
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
351 views

Holistic Numerical Methods

This document provides an example of using the 4th order Runge-Kutta method to solve an ordinary differential equation modeling the concentration of salt in a tank over time. The method is applied with different step sizes to calculate the salt concentration at 3 minutes. Tables and figures show that decreasing the step size improves the accuracy of the numerical solution compared to the exact solution. The document also compares the 4th order Runge-Kutta method to lower order Runge-Kutta methods.

Uploaded by

Ohol Rohan Bhaskar
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 6

Chapter 08.

04
Runge-Kutta 4th Order Method for Ordinary
Differential Equations-More Examples
Chemical Engineering
Example 1
The concentration of salt x in a home made soap maker is given as a function of time by
dx
 37.5  3.5 x
dt
At the initial time, t  0 , the salt concentration in the tank is 50 g/L Using Runge-Kutta 4th
order method and a step size of, h  1.5 min , what is the salt concentration after 3 minutes?
Solution
dx
 37.5  3.5 x
dt
f t , x  37.5  3.5x
xi 1  xi  k1  2k 2  2k 3  k 4 h
1
6
For i  0 , t 0  0 , x0  50
k1  f t 0 , x0 
 f 0,50
 37.5  3.5(50)
 137.5
 1 1 
k 2  f  t 0  h, x0  k1 h 
 2 2 
 
 f  0  1.5,50   137.51.5 
1 1
 2 2 
 f 0.75,53.125
 37.5  3.5(53.125)
 223.44
 1 1 
k 3  f  t 0  h, x0  k 2 h 
 2 2 
 
 f  0  1.5,50  223.441.5 
1 1
 2 2 

08.04.1
08.04.2 Chapter 08.04

 f 0.75,217.58
 37.5  3.5(217.58)
 724.02
k 4  f t 0  h, x0  k 3 h
 f 0  1.5,50   724.031.5
 f 1.5,1036.0
 37.5  3.5(1036.0)
 3663.6
1
x1  x0  (k1  2k 2  2k 3  k 4 )h
6
 50   137.5  2223.44  2 724.02  3663.6 1.5
1
6
 50  2525.0 1.5
1
6
 681.24 g/L
x1 is the approximate concentration of salt at
t  t1  t 0  h  0  1.5  1.5
x1.5  x1  681.24 g/L
For i  1, t1  1.5, x1  681.24
k1  f t1 , x1 
 f 1.5,681.24
 37.5  3.5(681.24)
 2346.8
 1 1 
k 2  f  t1  h, x1  k1 h 
 2 2 
 
 f 1.5  1.5, 681.24   2346.81.5 
1 1
 2 2 
 f 2.25,1078.9
 37.5  3.5(1078.9)
 3813.6
 1 1 
k 3  f  t1  h, x1  k 2 h 
 2 2 
 
 f 1.5  1.5, 681.24  3813.61.5 
1 1
 2 2 
 f 2.25, 3541.4
 37.5  3.5(3541.4)
 12358

k 4  f t1  h, x1  k 3 h 
Runge-Kutta 4th order Method for ODE-More Examples: Chemical Engineering 08.04.3

 f 1.5  1.5, 681.24   123581.5


 f 3,  17855
 37.5  3.5(17855)
 62530
1
x 2  x1  (k1  2k 2  2k 3  k 4 )h
6
 681.24   2346.8  23813.6  2 12358  625301.5
1
6
 681.24  430961.5
1
6
 11455 g/L
x2 is the approximate concentration of salt at
t 2  t1  h = 1.5  1.5  3 min
x3  x2  11455 g/L

The exact solution of the ordinary differential equation is given by


xt   10.714  39.286e 3.5t
The solution to this nonlinear equation at t  3 min is
x3  10.715 g/L

Figure 1 compares the exact solution with the numerical solution using Runge-Kutta 4th order
method using different step sizes.
08.04.4 Chapter 08.04

Figure 1 Comparison of Runge-Kutta 4th order method with exact solution for
different step sizes.

Table 1 and Figure 2 show the effect of step size on the value of the calculated temperature at
t  3 min.

Table 1 Value of concentration of salt at 3 minutes for different step sizes.

Step size, h x3 Et |t | %


3 14120 14109 131680
1.5 11455 11444 106800
0.75 25.559 14.843 138.53
0.375 10.717 0.0014969 0.013969
0.1875 10.715 0.00031657 0.0029544
Runge-Kutta 4th order Method for ODE-More Examples: Chemical Engineering 08.04.5

Figure 2 Effect of step size in Runge-Kutta 4th order method.

In Figure 3, we are comparing the exact results with Euler’s method (Runge-Kutta 1st order
method), Heun’s method (Runge-Kutta 2nd order method) and Runge-Kutta 4th order method.
08.04.6 Chapter 08.04

Figure 3 Comparison of Runge-Kutta methods of 1st, 2nd, and 4th order.

You might also like