Mathematics and Statistics 9(2): 166-171, 2021 http://www.hrpub.org
DOI: 10.13189/ms.2021.090211
Solving One-Dimensional Porous Medium Equation
Using Unconditionally Stable Half-Sweep Finite
Difference and SOR Method
Jackel Vui Lung Chew
1,*
, Jumat Sulaiman
2
, Andang Sunarto
3
1
Faculty of Computing and Informatics, Universiti Malaysia Sabah Labuan International Campus, 87000, Labuan F.T., Malaysia
2
Faculty of Science and Natural Resources, Universiti Malaysia Sabah, 88400 Kota Kinabalu, Sabah, Malaysia
3
IAIN Bengkulu, Jl Raden Fatah Kota Bengkulu, Bengkulu, Indonesia
Received December 19, 2020; Revised March 15, 2021; Accepted March 24, 2021
Cite This Paper in the following Citation Styles
(a): [1] Jackel Vui Lung Chew, Jumat Sulaiman, Andang Sunarto , "Solving One-Dimensional Porous Medium Equation
Using Unconditionally Stable Half-Sweep Finite Difference and SOR Method," Mathematics and Statistics, Vol. 9, No. 2,
pp. 166 - 171, 2021. DOI: 10.13189/ms.2021.090211.
(b): Jackel Vui Lung Chew, Jumat Sulaiman, Andang Sunarto (2021). Solving One-Dimensional Porous Medium
Equation Using Unconditionally Stable Half-Sweep Finite Difference and SOR Method. Mathematics and Statistics, 9(2),
166 - 171. DOI: 10.13189/ms.2021.090211.
Copyright©2021 by authors, all rights reserved. Authors agree that this article remains permanently open access under
the terms of the Creative Commons Attribution License 4.0 International License
Abstract A porous medium equation is a nonlinear
parabolic partial differential equation that presents many
physical occurrences. The solutions of the porous medium
equation are important to facilitate the investigation on
nonlinear processes involving fluid flow, heat transfer,
diffusion of gas-particles or population dynamics. As part
of the development of a family of efficient iterative
methods to solve the porous medium equation, the
Half-Sweep technique has been adopted. Prior works in the
existing literature on the application of Half-Sweep to
successfully approximate the solutions of several types of
mathematical problems are the underlying motivation of
this research. This work aims to solve the one-dimensional
porous medium equation efficiently by incorporating the
Half-Sweep technique in the formulation of an
unconditionally-stable implicit finite difference scheme.
The noticeable unique property of Half-Sweep is its ability
to secure a low computational complexity in computing
numerical solutions. This work involves the application of
the Half-Sweep finite difference scheme on the general
porous medium equation, until the formulation of a
nonlinear approximation function. The Newton method is
used to linearize the formulated Half-Sweep finite
difference approximation, so that the linear system in the
form of a matrix can be constructed. Next, the Successive
Over Relaxation method with a single parameter was
applied to efficiently solve the generated linear system per
time step. Next, to evaluate the efficiency of the developed
method, deemed as the Half-Sweep Newton Successive
Over Relaxation (HSNSOR) method, the criteria such as
the number of iterations, the program execution time and
the magnitude of absolute errors were investigated.
According to the numerical results, the numerical solutions
obtained by the HSNSOR are as accurate as those of the
Half-Sweep Newton Gauss-Seidel (HSNGS), which is
under the same family of Half-Sweep iterations, and the
benchmark, Newton-Gauss-Seidel (NGS) method. The
improvement in the numerical results produced by the
HSNSOR is significant, and requires a lesser number of
iterations and a shorter program execution time, as
compared to the HSNGS and NGS methods.
Keywords One-Dimensional Porous Medium
Equation, Half-Sweep, Finite Difference Method, Newton,
Successive Over Relaxation, Iterative Method
1. Introduction
Porous medium equation (PME) is a nonlinear parabolic
partial differential equation that exists in many nonlinear