Skip to main content Accessibility help

A Novel Technique for Constructing Difference Schemes for Systems of Singularly Perturbed Equations

  • Po-Wen Hsieh (a1), Yin-Tzer Shih (a2), Suh-Yuh Yang (a3) and Cheng-Shu You (a3)


In this paper, we propose a novel and simple technique to construct effective difference schemes for solving systems of singularly perturbed convection-diffusion-reaction equations, whose solutions may display boundary or interior layers. We illustrate the technique by taking the Il'in-Allen-Southwell scheme for 1-D scalar equations as a basis to derive a formally second-order scheme for 1-D coupled systems and then extend the scheme to 2-D case by employing an alternating direction approach. Numerical examples are given to demonstrate the high performance of the obtained scheme on uniform meshes as well as piecewise-uniform Shishkin meshes.


Corresponding author

*Corresponding author. Email (P.-W. Hsieh), (Y.-T. Shih), (S.-Y. Yang), (C.-S. You)


Hide All
[1]Fletcher, C. A. J., Generating exact solutions of the two-dimensional Burgers' equation, Int. J. Numer. Methods Fluids, 3 (1983), pp. 213216.
[2]Hsieh, P.-W., Shih, Y., and Yang, S.-Y., A tailored finite point method for solving steady MHD duct flow problems with boundary layers, Commun. Comput. Phys., 10 (2011), pp. 161182.
[3]Hsieh, P.-W. and Yang, S.-Y., Two new upwind difference schemes for a coupled system of convection-diffusion equations arising from the steady MHD duct flow problems, J. Comput. Phys., 229 (2010), pp. 92169234.
[4]Hsieh, P.-W., Yang, S.-Y., and You, C.-S., A high-accuracy finite difference scheme for solving reaction-convection-diffusion problems with a small diffusivity, Adv. Appl. Math. Mech., 6 (2014), pp. 637662.
[5]Kelleci, A. and Yildirim, A., An efficient numerical method for solving coupled Burgers' equation by combining homotopy perturbation and Pade techniques, Numer. Methods Partial Differential Eq., 27 (2011), pp. 982995.
[6]Kokotović, P. V., Applications of singular perturbation techniques to control problems, SIAM Rev., 26 (1984), pp. 501550.
[7]Morton, K. W., Numerical Solution of Convection-Diffusion Problems, Chapman & Hall, London, UK, 1996.
[8]O'Riordan, E. and Stynes, M., Numerical analysis of a strongly coupled system of two singularly perturbed convection-diffusion problems, Adv. Comput. Math., 30 (2009), pp. 101121.
[9]O'Riordan, E., Stynes, J., and Stynes, M., A parameter-uniform finite difference method for a coupled system of convection-diffusion two-point boundary value problems, Numer. Math. Theor. Meth. Appl., 1 (2008), pp. 176197.
[10]O'Riordan, E., Stynes, J., and Stynes, M., An iterative numerical algorithm for a strongly coupled system of singularly perturbed convection-diffusion problems, in NAA 2008, LNCS 5434, Springer-Verlag, Berlin Heidelberg, 2009, pp. 104115.
[11]Roos, H.-G., Ten ways to generate the Il'in and related schemes, J. Comput. Appl. Math., 53 (1994), pp. 4359.
[12]Roos, H.-G., Stynes, M., and Tobiska, L., Robust Numerical Methods for Singularly Perturbed Differential Equations, Second Edition, Springer-Verlag, Berlin, 2008.
[13]Sanyasiraju, Y. V. S. S. and Mishra, N., Exponential compact higher order scheme for nonlinear steady convection-diffusion equations, Commun. Comput. Phys., 9 (2011), 897916.
[14]Shishkin, G. I., Mesh approximation of singularly perturbed boundary-value problems for systems of elliptic and parabolic equations, Comput. Maths. Math. Phys., 35 (1995), pp. 429446.
[15]Thomas, G. P., Towards an improved turbulence model for wave-current interactions, in Second Annual Report to EU MAST-III Project “The Kinematics and Dynamics of Wave-Current Interactions,” 1998.
[16]Tian, Z. F. and Dai, S. Q., High-order compact exponential finite difference methods for convection-diffusion type problems, J. Comput. Phys., 220 (2007), pp. 952974.
[17]Zhu, H., Shu, H. and Ding, M., Numerical solutions of two-dimensional Burgers' equations by discrete Adomian decomposition method, Comput. Math. Appl., 60 (2010), pp. 840848.


MSC classification


Full text views

Total number of HTML views: 0
Total number of PDF views: 0 *
Loading metrics...

Abstract views

Total abstract views: 0 *
Loading metrics...

* Views captured on Cambridge Core between <date>. This data will be updated every 24 hours.

Usage data cannot currently be displayed