Abstract
We propose a nonstandard finite difference method to integrate coupled systems of singularly perturbed convection–diffusion equations. We prove that the underlying discrete operator satisfies a stability property in the maximum norm. This fact is used to prove that the proposed method converges uniformly with respect to the singular perturbation parameters and the convergence rate is linear with respect to the step size. Numerical investigations are presented to confirm our theoretical findings. We also show that the proposed approach compares well with existing methods.
Similar content being viewed by others
References
Bellew S, O’Riordan E (2004) A parameter robust numerical method for a system of two singularly perturbed convection-diffusion equations. Appl Num Math 51:171–186
Farrell PA, Hegarty AF, Miller JJH, O’Riordan E, Shishkin GI (2000) Robust computational techniques for boundary layers. CRC Press, Boca Raton
Kadalbajoo MK, Patidar KC (2002) A survey of numerical techniques for solving singularly perturbed ordinary differential equations. Appl Math Comput 130:457–510
Kadalbajoo MK, Gupta V (2010) A brief survey on numerical methods for solving singularly perturbed problems. Appl Math Comput 217:3641–3716
Linß T, Madden N (2004) Accurate solution of a system of coupled singularly perturbed reaction-diffusion equations. Computing 73:121–133
Linß T (2005) On a set of singularly perturbed convection-diffusion equations. J Comput Appl Math 180:173–179
Linß T (2007) Analysis of an upwind finite-difference scheme for a system of coupled singularly perturbed convection-diffusion equations. Computing 79:23–32
Linß T (2009) Analysis of a system of singularly perturbed convection-diffusion equations with strong coupling. SIAM J Numer Anal 47(3):1847–1862
Linß T (2010) Layer-adapted meshes for reaction-convection-diffusion problems. Springer, Berlin
Lubuma JM-S, Patidar KC (2007) Non-standard methods for singularly perturbed problems possessing oscillatory/layer solutions. Appl Math Comput 187:1147–1160
Madden N (2000) Numerical methods for wave-current interactions. National University of Ireland, Cork Ph.D
Madden N, Stynes M (2003) A uniformly convergent numerical method for a coupled system of two singularly perturbed linear reaction-diffsion problems. IMA J Numer Anal 23:627–644
Matthews S (2000) Parameter robust methods for a system of coupled singularly perturbed ordinary differential reaction-diffusion equations. School of Mathematical Sciences, Dublin City University, M.Sc
Matthews S, O’Riordan E, Shishkin GI (2002) A numerical method for a system of singularly perturbed reaction-diffusion equations. J Comput Appl Math 145:151–166
Mickens RE (1989) Exact solutions to a finite difference model of a reaction-advection equation: Implications for numerical analysis. Numer Methods Partial Differ. Equ. 4:313–325
Mickens RE (ed) (2005) Advances in the applications of nonstandard finite difference schemes. World Scientific, Singapore
Miller JJH, O’Riordan E, Shishkin GI (1996) Fitted numerical methods for singular perturbation problems. World Scientific, Singapore
Morton KW (1996) Numerical solution of convection-diffusion problems. Applied mathematics and mathematical computation, vol 12. Chapman and Hall, London
Munyakazi JB, Patidar KC (2010) Limitations of Richardson’s extrapolation for a high order fitted mesh method for self-adjoint singularly perturbed problems. J Appl Math Comput 32:219–236
Munyakazi JB, Patidar KC (2010) Higher order numerical methods for singularly perturbed elliptic problems. Neural Parallel Sci Comput 18:75–88
Munyakazi JB, Patidar KC (2012) Novel fitted operator finite difference methods for singularly perturbed elliptic convection-diffusion problems in two dimensions. J Differ Equ Appl 18(5):799–813
Munyakazi JB, Patidar KC (2012) An efficient numerical method for a system of singularly perturbed reaction-diffusion equations. Neural Parallel Sci Comput 20:173–190
O’Riordan E, Stynes M (2009) Numerical analysis of a strongly coupled system of two singularly perturbed convection-diffusion problems. Adv Comput Math 30:101–121
Patidar KC (2005) On the use of nonstandard finite difference methods. J Differ Equ Appl 11(8):735–758
Prandtl L (1904) Über Flussigkeitsbewegung bei sehr kleiner Reibung. Proc. III Intern. Congr Math, Heildeberg
Roos H-G, Stynes M, Tobiska L (1996) Numerical methods for singularly perturbed differential equations. Springer, Berlin
Roos H-G (2012) Special features of strongly coupled systems of convection-diffusion equations with two small parameters. Appl Math Lett 25:1127–1130
Shishkin GI (1995) Mesh approximation of singularly perturbed boundary-value problems for systems of elliptic and parabolic equation. Comput Maths Math Phys 35(4):429–446
Stephens M, Madden N (2009) A parameter-uniform Schwarz method for a coupled system of reaction-diffusion equations. J Comput Appl Math 230:360–370
Zhongdi CEN (2005) Parameter-uniform finite difference scheme for a system of coupled singularly perturbed convection-diffusion equations. J Syst Sci Complex 18(4):498–510
Acknowledgments
This work was completed while the author was a visiting scholar on the University of Michigan African Presidential Scholars (UMAPS) program. He wishes to thank Professor Robert Krasny for reading this work and suggesting improvements. He is also indebted to Professor Kailash C. Patidar for introducing him to the field of nonstandard finite difference methods.
Author information
Authors and Affiliations
Corresponding author
Additional information
Communicated by Raphaèle Herbin.
Rights and permissions
About this article
Cite this article
Munyakazi, J.B. A uniformly convergent nonstandard finite difference scheme for a system of convection–diffusion equations. Comp. Appl. Math. 34, 1153–1165 (2015). https://doi.org/10.1007/s40314-014-0171-6
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s40314-014-0171-6
Keywords
- System of equations
- Convection–diffusion
- Singular perturbation
- Parameter-uniform
- Nonstandard finite difference schemes
- Convergence analysis
- Error bounds