Solving one-dimensional Porous Medium Equation using unconditionally stable Half-Sweep finite difference and SOR method
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...
| Main Authors: | , , |
|---|---|
| Format: | Article |
| Language: | English English |
| Published: |
Horizon Research Publishing
2021
|
| Subjects: | |
| Online Access: | https://eprints.ums.edu.my/id/eprint/37195/1/ABSTRACT.pdf https://eprints.ums.edu.my/id/eprint/37195/2/FULL%20TEXT.pdf |
| _version_ | 1825715343711010816 |
|---|---|
| author | Jackel Vui Lung Chew Jumat Sulaiman Andang Sunarto |
| author_facet | Jackel Vui Lung Chew Jumat Sulaiman Andang Sunarto |
| author_sort | Jackel Vui Lung Chew |
| collection | UMS |
| description | 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. |
| first_indexed | 2024-03-06T03:25:25Z |
| format | Article |
| id | ums.eprints-37195 |
| institution | Universiti Malaysia Sabah |
| language | English English |
| last_indexed | 2024-03-06T03:25:25Z |
| publishDate | 2021 |
| publisher | Horizon Research Publishing |
| record_format | dspace |
| spelling | ums.eprints-371952023-09-19T01:59:49Z https://eprints.ums.edu.my/id/eprint/37195/ Solving one-dimensional Porous Medium Equation using unconditionally stable Half-Sweep finite difference and SOR method Jackel Vui Lung Chew Jumat Sulaiman Andang Sunarto TA1-2040 Engineering (General). Civil engineering (General) TA401-492 Materials of engineering and construction. Mechanics of materials 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. Horizon Research Publishing 2021-12 Article NonPeerReviewed text en https://eprints.ums.edu.my/id/eprint/37195/1/ABSTRACT.pdf text en https://eprints.ums.edu.my/id/eprint/37195/2/FULL%20TEXT.pdf Jackel Vui Lung Chew and Jumat Sulaiman and Andang Sunarto (2021) Solving one-dimensional Porous Medium Equation using unconditionally stable Half-Sweep finite difference and SOR method. Mathematics and Statistics, 9. pp. 166-171. ISSN 2332-2071 https://doi.org/10.13189/ms.2021.090211 |
| spellingShingle | TA1-2040 Engineering (General). Civil engineering (General) TA401-492 Materials of engineering and construction. Mechanics of materials Jackel Vui Lung Chew Jumat Sulaiman Andang Sunarto Solving one-dimensional Porous Medium Equation using unconditionally stable Half-Sweep finite difference and SOR method |
| title | Solving one-dimensional Porous Medium Equation using unconditionally stable Half-Sweep finite difference and SOR method |
| title_full | Solving one-dimensional Porous Medium Equation using unconditionally stable Half-Sweep finite difference and SOR method |
| title_fullStr | Solving one-dimensional Porous Medium Equation using unconditionally stable Half-Sweep finite difference and SOR method |
| title_full_unstemmed | Solving one-dimensional Porous Medium Equation using unconditionally stable Half-Sweep finite difference and SOR method |
| title_short | Solving one-dimensional Porous Medium Equation using unconditionally stable Half-Sweep finite difference and SOR method |
| title_sort | solving one dimensional porous medium equation using unconditionally stable half sweep finite difference and sor method |
| topic | TA1-2040 Engineering (General). Civil engineering (General) TA401-492 Materials of engineering and construction. Mechanics of materials |
| url | https://eprints.ums.edu.my/id/eprint/37195/1/ABSTRACT.pdf https://eprints.ums.edu.my/id/eprint/37195/2/FULL%20TEXT.pdf |
| work_keys_str_mv | AT jackelvuilungchew solvingonedimensionalporousmediumequationusingunconditionallystablehalfsweepfinitedifferenceandsormethod AT jumatsulaiman solvingonedimensionalporousmediumequationusingunconditionallystablehalfsweepfinitedifferenceandsormethod AT andangsunarto solvingonedimensionalporousmediumequationusingunconditionallystablehalfsweepfinitedifferenceandsormethod |