A robust fitted finite difference method for semi-linear two-parameter singularly perturbed PDEs

Document Type : Research Article

Authors

1 Department of Mathematics, College of Natural and Computational Science, Wolkite University, Wolkite, Ethiopia

2 Department of Mathematics and Applied Mathematics University of Western Cape

3 Department of Applied Mathematics Adama Science and Technology University

4 Department of Mathematics, Arsi University, Asela, Ethiopia

5 Department of Mathematics, Oda Bultum University, Chiro, Ethiopia

Abstract

In this article, a new numerical approach is developed for nonlinear two-parameter singularly perturbed initial-boundary value problems. The implicit backward Euler discretization for the time derivative and the fitted operator technique in the spatial domain are employed. Newton's quasilinearization technique is applied to the nonlinear terms. An investigation of parameter-uniform error estimates shows that the developed approach is first-order accurate in both time and space. However, a temporal mesh refinement technique is introduced to improve the order of accuracy to two. Two examples are provided and implemented in Python to validate the applicability of the method, and the results are displayed in tables and graphs.

Keywords

Main Subjects


[1] K. Bansal, P. Rai, K.K. Sharma, Numerical treatment for the class of time dependent singularly perturbed parabolic problems with general shift arguments, Differ. Equ. Dyn. Syst. 25 (2017) 327–346.
[2] R.E. Bellman, R.E. Kalaba, Quasilinearization and Nonlinear Boundary-Value Problems, Elsevier, 1965.
[3] P. Chakravarthy, R.N. Rao, A modified numerov method for solving singularly perturbed differential-difference equations arising in science and engineering, Results Phys. 2 (2012) 100–103.
[4] A. Das, L. Govindarao, J. Mohapatra, A second order weighted monotone numerical scheme for time-delayed parabolic initial-boundary-value problem involving a small parameter, Int. J. Math. Model. Numer. Optim. 12 (2022) 233–251.
[5] P. Deepika, A. Das, A robust numerical scheme via grid equidistribution for singularly perturbed delay partial differential equations arising in control theory, Int. J. Appl. Comput. Math. 10 (2024) 72.
[6] E.P. Doolan, J.J. Miller, W.H. Schilders, Uniform Numerical Methods for Problems with Initial and Boundary Layers, Boole Press, 1980.
[7] G.F. Duressa, T.B. Mekonnen, An exponentially fitted method for two parameter singularly perturbed parabolic boundary value problems, Commun. Korean Math. Soc 38 (2023) 299–318.
[8] S. Elango, L. Govindarao, R. Vadivel, A comparative study on numerical methods for fredholm integro-differential equations of convection-diffusion problem with integral boundary conditions, Appl. Numer. Math. 207 (2025) 323–338.
[9] S. Elango, L. Govindarao, M. Awadalla, H. Zaway, Efficient numerical methods for reaction diffusion problems governed by singularly perturbed fredholm integro-differential equations, Mathematics 13 (2025) 2227–7390.
[10] L. Govindarao, A. Das, A second-order fractional step method for two-dimensional delay parabolic partial differential equations with a small parameter, Palest. J. Math. 11 (2022) 96–111.
[11] L. Govindarao, H. Ramos, S. Elango, Numerical scheme for singularly perturbed fredholm integro differential equations with non-local boundary conditions, Comput. Appl. Math. 43 (2024) 126.
[12] J.L. Gracia, E. O'Riordan, M. Pickett, A parameter robust second order numerical method for a singularly perturbed two-parameter problem, Appl. Numer. Math. 56 (2006) 962–980.
[13] V. Gupta, M.K. Kadalbajoo, A layer adaptive b-spline collocation method for singularly perturbed one-dimensional parabolic problem with a boundary turning point, Numer. Methods Partial Differ. Equ. 27 (2011) 1143–1164.
[14] P.W. Hemker, G.I. Shishkin, L.P. Shishkina, ε-uniform schemes with high-order time-accuracy for parabolic singular perturbation problems, IMA J. Numer. Anal. 20 (2000) 99–121.
[15] M.K. Kadalbajoo, A. Awasthi, A parameter uniform difference scheme for singularly perturbed parabolic problem in one space dimension, Appl. Math. Comput. 183 (2006) 42–60.
[16] M.K. Kadalbajoo, V. Gupta, Numerical solution of singularly perturbed convection-diffusion problem using parameter uniform b-spline collocation method, J. Math. Anal. Appl. 355 (2009) 439–452.
[17] M.K. Kadalbajoo, V. Gupta, A. Awasthi, A uniformly convergent b-spline collocation method on a nonuniform mesh for singularly perturbed one-dimensional time-dependent linear convection diffusion problem, J. Comput. Appl. Math. 220 (2008) 271–289.
[18] M.K. Kadalbajoo, A.S. Yadaw, B-spline collocation method for a two-parameter singularly perturbed convection-diffusion boundary value problems, Appl. Math. Comput. 201 (2008) 504–513.
[19] B. Lin, K. Li, Z. Cheng, B-spline solution of a singularly perturbed boundary value problem arising in biology, Chaos Solit. Fractals 42 (2009) 2934–2948.
[20] T. Linß, H.G. Roos, Analysis of a finite-difference scheme for a singularly perturbed problem with two small parameters, J. Math. Anal. Appl. 289 (2004) 355–366.
[21] R.E. Mickens, Advances in the Applications of Nonstandard Finite Difference Schemes, World Scientific, 2005.
[22] J. Mohapatra, L. Govindarao, S. Priyadarshana, A splitting based higher-order numerical scheme for 2d time-dependent singularly perturbed reaction-diffusion problems, J. Supercomput. 81 (2025) 203.
[23] M.A. Mohye, J.B. Munyakazi, T.G. Dinka, A second-order numerical method for two-parameter singularly perturbed time-delay parabolic problems, J. Math. Model. 11 (2023) 745–746.
[24] J.B. Munyakazi, K.C. Patidar, A fitted numerical method for singularly perturbed parabolic reaction-diffusion problems, Comput. Appl. Math. 32 (2013) 509–519.
[25] N.T. Negero, A parameter-uniform efficient numerical scheme for singularly perturbed time-delay parabolic problems with two small parameters, Partial Differ. Equ. Appl. Math. 7 (2023) 100518.
[26] M.J. Ng-Stynes, E. O'Riordan, M. Stynes, Numerical methods for time-dependent convection diffusion equations, J. Comput. Appl. Math. 21 (1988) 289–310.
[27] S. Priyadarshana, J. Mohapatra, An efficient fractional step numerical algorithm for time-delayed singularly perturbed 2d convection-diffusion–reaction problem with two small parameters, Numer. Algorithms 97 (2024) 687–726.
[28] S. Priyadarshana, Monotone hybrid numerical method for singularly perturbed time-lagged semilinear parabolic problems, Nat. Acad. Sci. Letters 46 (2023) 347–350.
[29] S. Priyadarshana, J. Mohapatra, An efficient computational technique for time dependent semilinear parabolic problems involving two small parameters, J. Appl. Math. Comput. 69 (2023) 3721–3754.
[30] S. Priyadarshana, J. Mohapatra, Weighted variable based numerical scheme for time-lagged semilinear parabolic problems including small parameter, J. Appl. Math. Comput. 69 (2023) 2439–2463.
[31] S. Priyadarshana, J. Mohapatra, S.R. Pattanaik, Parameter uniform optimal order numerical approximations for time-delayed parabolic convection diffusion problems involving two small parameters, Comput. Appl. Math. 41 (2022) 233.
[32] S. Priyadarshana, J. Mohapatra, S.R. Pattanaik, An improved time accurate numerical estimation for singularly perturbed semilinear parabolic differential equations with small space shifts and a large time lag, Math. Comput. Simul. 214 (2023) 183–203.
[33] C.S. Rao, M. Kumar, Exponential b-spline collocation method for self-adjoint singularly perturbed boundary value problems, Appl. Numer. Math. 58 (2008) 1572–1581.
[34] G.I. Shishkin, On finite difference fitted schemes for singularly perturbed boundary value problems with a parabolic boundary layer, J. Math. Anal. Appl. 208 (1997) 181–204.
[35] M. Udupa, S. Saha, A. Das, Computational analysis of blood flow through stenosed artery under the influence of body acceleration using shishkin mesh, Int. J. Appl. Comput. Math. 8 (2022) 107.
[36] A. Wineman, Nonlinear viscoelastic membranes, Comput. Math. Appl. 53 (2007) 168–181.
[37] W.K. Zahra, A.M. ElMhlawy, Numerical solution of two-parameter singularly perturbed boundary value problems via exponential spline, J. King Saud Univ. Sci. 25 (2013) 201–208.