KEMBAR78
Holistic Numerical Methods | PDF | Mathematics Of Computing | Subtraction
Scribd
Upload
513 views6 pages

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
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
513 views6 pages

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
We take content rights seriously. If you suspect this is your content, claim it here.
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