B-Spline Collocation Method for Solving Singularly Perturbed Boundary Value Problems


We use fifth order B-spline functions to construct the numerical method for solving singularly perturbed boundary value problems. We use B-spline collocation method, which leads to a tri-diagonal linear system. The accuracy of the proposed method is demonstrated by test problems. The numerical results are found in good agreement with exact solutions.

Share and Cite:

Lin, B. (2016) B-Spline Collocation Method for Solving Singularly Perturbed Boundary Value Problems. Journal of Applied Mathematics and Physics, 4, 1699-1704. doi: 10.4236/jamp.2016.49178.

1. Introduction

Consider following singularly perturbed boundary value problem


with boundary conditions


B-spline functions are useful wavelet basis functions; the stiffness matrix is sparse when it is used as trial functions. B-splines were introduced by Schoenberg in 1946 [5] . Up to now, B-spline approximation method for numerical solutions has been researched by various researchers [6] - [8] .

2. Description of the B-Spline Collocation Method

The expression of fifth order B-spline function is as follows:


The fifth order B-spline function is used to calculate in this work and possesses the following characters: piecewise smooth, compact support, Symmetry, rapidly decaying, differentiability, linear combination.

The region [a,b] is partitioned into uniformly sized finite elements of length h by the knots such that with, ,. Let be fifth order B-spline function with knots at the points,. The set of splines forms a basis for functions defined over [a,b].

In the proposed algorithm, The fifth order B-spline function is used as a single mother wavelet, i.e. and dilation and translation of mother wavelet functions can construct any function of.



So the global approximation to the function can be written in terms of the B-spline as follows

, (5)

where, are unknown real coefficients.

Using the fifth order B-spline function and the approximate solution Equation (5), the nodal values, and at the node are given in terms of element parameters by




where the symbols and denote first and second differentiation with respect to x, respectively.

Substituting Equations (6)-(8) into Equation (1) and Equation (2), we can obtain following linear equations






It is easily seen that the matrix B is strictly diagonally dominant and hence nonsingular. Since B is nonsingular, we can solve the system for . Hence the method of collocation using the fifth order B-spline function as a basis function applied to the singularly perturbed boundary value problem has a unique solution given by Equation (5).

3. Numerical Results

In the section, we illustrate the numerical techniques discussed in the previous section by the following problems.

Example 1. Consider the convention-dominated equation:


with boundary conditions:,



The exact solution is given by



Comparison of the numerical results and point-wise errors is given in Table 1.

It observed that

1) when h decreases (i.e. collocation number increases) for fixed the point-wise errors decrease;

2) when decreases for fixed h the point-wise errors increase;

3) when, the errors are very large.

Example 2. Solve the following non-homogeneous equation:


with boundary conditions

, ,.

The analytical solution is given by




And and are the real solutions of the characteristic equation.

Approximation solutions for different values of and for fixed p are given in Figure 1. It observed that

Table 1. Example 1. Comparison of results and point-wise errors.

Figure 1. Approximation solutions of example 2 for different values of epsilon g and for fixed p.

1) when and, the approximation solutions are in good agreement with exact solution; 2) when and, and the errors are very large; 3) when decreases for fixed p the width of boundary layer becomes small and wave shape change more and more stiff at and.

4. Conclusion

The numerical results show clearly the effect of on the boundary layer and the B-spline collocation method solving singular boundary value problems is relatively simple to collocate the solution at the mesh points. It is applicable technique and approximates the exact solution very well.


The authors would like to thank the editor and the reviewers for their valuable comments and suggestions to improve the results of this paper. This work was supported by the Natural Science Foundation of Guangdong (No. 2015A030313827).

Conflicts of Interest

The authors declare no conflicts of interest.


[1] Evrenosoglu, M. and Somali, S. (2008) Least Squares Methods for Solving Singularly Perturbed Two-Point Boundary Value Problems Using Bezier Control Points. Applied Mathematics Letters, 21, 1029-1032.
[2] Lin, B., Li, K.T. and Cheng, Z.X. (2009) B-Spline Solution of a Singularly Perturbed Boundary Value Problem Arising in Biology. Chaos, Solitons & Fractals, 42, 2934-2948.
[3] Fan, X.H. and Tian, L.X. (2005) The Existence of Solitary Waves of Singularly Perturbed mKdV-KS Equation. Chaos, Solitons & Fractals, 26, 1111-1118.
[4] Bigge, J. and Bohl, E. (1985) Deformations of the Bifurcation Diagram Due to Discretization. Mathematics of Computation, 45, 393-403.
[5] De Boor, C. (1978) A Practical Guide to Splines. Springer-Verlag, Berlin.
[6] Siddiqi, S.S. and Akram, G. (2007) Sextic Spline Solutions of Fifth Order Boundary Value Problems. Applied Mathematics Letters, 20, 591-597.
[7] Caglar, H.N. and Caglar, S.H. (1997) The Numerical Solution of Fifth Order Boundary Value Problems with Sixth Degree B-Spline Functions. Applied Mathematics Letters, 12, 25-30. http://dx.doi.org/10.1016/S0893-9659(99)00052-X
[8] Caglar, H., Özer, M. and Caglar, N. (2008) The Numerical Solution of the One-Dimensional Heat Equation by Using Third Degree B-Spline Functions. Chaos, Solitons & Fractals, 38, 1197-1201.

Copyright © 2024 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.