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,QuasilinearizationandNonlinearBoundary-ValueProblems,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 per
turbed 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] Gracia, E. O’Riordan, M. Pickett, A parameter robust second order numerical method for a singu
larly 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 Part. Differ.
Equ. 27 (2011) 1143–1164.
[14] Hemker, G. Shishkin, L. 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 prob
lem 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 per
turbed 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. Equat. 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 semi
linear 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 semi
linear 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 ap
proximations for time-delayed parabolic convection diffusion problems involving two small param
eters, 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, Numericalsolution of two-parameter singularly perturbed boundary
value problems via exponential spline, J. King Saud Univ. Sci. 25 (2013) 201–208