Advanced
SPLINE DIFFERENCE SCHEME FOR TWO-PARAMETER SINGULARLY PERTURBED PARTIAL DIFFERENTIAL EQUATIONS
SPLINE DIFFERENCE SCHEME FOR TWO-PARAMETER SINGULARLY PERTURBED PARTIAL DIFFERENTIAL EQUATIONS
Journal of Applied Mathematics & Informatics. 2014. Jan, 32(1_2): 185-201
Copyright © 2014, Korean Society of Computational and Applied Mathematics
  • Received : January 30, 2013
  • Accepted : May 27, 2013
  • Published : January 28, 2014
Download
PDF
e-PUB
PubReader
PPT
Export by style
Share
Article
Author
Metrics
Cited by
TagCloud
About the Authors
W. K. ZAHRA
M. S. EL-AZAB
ASHRAF M. EL MHLAWY

Abstract
In this paper, we construct a numerical method to solve singularly perturbed one-dimensional parabolic convection–diffusion problems. We use Euler method with uniform step size for temporal discretization and exponential-spline scheme on spatial uniform mesh of Shishkin type for full discretization. We show that the resulting method is uniformly convergent with respect to diffusion parameter. An extensive amount of analysis has been carried out to prove the uniform convergence with respect to the singular perturbation parameter. The obtained numerical results show that the method is efficient, stable and reliable for solving convection–diffusion problem accurately even involving diffusion parameter. AMS Mathematics Subject Classification : 65H05.
Keywords
1. Introduction
In the present paper, we are concerned with the numerical approximation of the partial differential equation:
PPT Slide
Lager Image
Here QT = Ω × I , Ω ≡ ( a , b ), ∂Ω ≡ { a , b } , I ≡ (0, T ), and a and b are real positive constants. We consider equation (1) associated with the Dirichlet boundary conditions
PPT Slide
Lager Image
and the initial condition
PPT Slide
Lager Image
where u ( x , t ) is the dependent variable, x and t are the independent variables and f , p , q and u 0 are given functions of their arguments. It is usually assumed that 0 < εc ≪ 1 and 0 < εd ≪ 1.
This problem encloses both the reaction–diffusion problem when εc = 0 and the convection–diffusion problem when εc = 1. These problems occur naturally in various fields of science and engineering, such as, heat transfer with large Peclet numbers, combustion, nuclear engineering, control theory, elasticity, fluid mechanics, aerodynamics, quantum mechanics, optimal control, chemicalreactor theory, convection-diffusion process and geophysics; see [9 , 10 , 18] .
There has been a lot of effort in developing numerical methods for the solution of singular perturbation problems that are uniformly convergent. Clavero et al. [1] and Kadalbajoo et al. [4] gives a uniform convergent numerical method with respect to the diffusion parameter to solve the one-dimensional time–dependent convection–diffusion problem. They used the implicit Euler method for the time discretization and the simple upwind finite difference scheme on a Shishkin mesh for the spatial discretization. Ramos [7] presented an exponentially fitted method for singularly perturbed, one-dimensional (convection–diffusion) parabolic problems, and showed its uniform convergence in the perturbation parameter. Surla and Jerkovic [13] considered a singularly perturbed boundary value problem using a spline collocation method. Sakai and Usmani [11 , 12] gave a new concept of B-spline in terms of hyperbolic and trigonometric splines which are different from the earlier ones. It is proved that the hyperbolic and trigonometric B-splines are characterized by a convolution of some special exponential functions and a characteristic function on the interval [0,1]. Sharma and Kaushik [14] presented a numerical approximation of a singularly perturbed time delayed convection diffusion problem on a rectangular domain. Finally, Zahra and El Mhlawy [18] developed a numerical method to solve two-parameter singularly perturbed semi-Linear boundary value problems. This method is based on exponential spline with a piecewise uniform Shishkin mesh which is shown to be uniformly convergent independent of mesh parameters and perturbation parameter εc and εd .
Our contribution in this paper is to introduce a full discretization scheme based on Rothe’s approximation in time discretization and on the exponential spline approximation in the spatial discretization for solving singularly perturbed, one-dimensional time-dependent convection–diffusion problem associated with homogeneous Dirichlet boundary conditions. It starts with the discretization in time by the 2-point Euler backward finite difference method in the time variable. After that we deal with the exponential spline method for the solution of the time discretized problem to compute the unknown function and the obtained system of algebraic equations is solved by iterative methods. Moreover, we obtain error estimates for the approximation process. The proposed technique is programmed using Matlab 7 , and used in solving the proposed problem.
The rest of the paper is organized as follows: In section 2, we give some notations, assumptions and definitions, In §3, a time discretization scheme for the continuous problem is proposed in addition to the derivation of some a priori estimates. In §4, we prove the convergence results and error estimates of the semidiscretization scheme described in §3. We describe in §5 a full discretization scheme based on Rothe- exponential spline methods and the error estimates of the full discretized solution are derived. Finally, numerical results and discussions are presented and comparisons are made with other solutions in §6.
2. Notations, assumptions and definitions
In the sequel, we will denote by (. , .) either the standard inner product in L 2 = L 2 (Ω) or the pairing between
PPT Slide
Lager Image
and V H −1 (see e.g. [11] ). we use the symbols |.|, ∥.∥ and ∥.∥ as the norms in L 2 (Ω), V , V , respectively. By
PPT Slide
Lager Image
we mean the strong and weak convergence. Also, we introduce some notations concerning the time discretization of our problem.
PPT Slide
Lager Image
for any given family
PPT Slide
Lager Image
The following elementary relations will be used in the following analysis:
PPT Slide
Lager Image
and Young’s inequality
PPT Slide
Lager Image
where ε is a small constant. We will assume, throughout this work, the following hypotheses on the given data.
( H1 ) The function f : Ω × I R is Lipschitz continuous in the sense of
PPT Slide
Lager Image
with f (0, 0) = 0.
( H2 ) The functions p ( x ) and q ( x ) are smooth functions and satisfy
PPT Slide
Lager Image
( H3 )
PPT Slide
Lager Image
Under these assumptions, we can define the weak solution of problem (1)-(2).
Definition 2.1. The measurable function u C ( I ; L 2 (Ω))∩ L 2 ( I ; V ) with ∂ tu L 2 ( I ; V ) and u ( x , 0) = u 0 ( x ) is said to be a weak (variational) solution of (1)−(3) if and only if the integral identity
PPT Slide
Lager Image
holds for all φ V and a.e. t I .
Remark. It is evident that for any functions α , β H 1 (Ω) and γ L 2 (Ω)
PPT Slide
Lager Image
3. The temporal discretization scheme, A prior estimate
Our main goal is to approximate (1)−(3) from a numerical point of view and to prove its convergence. The scheme starts with the discretization in time by the 2-point Euler backward finite difference method in the time variable, so that the problem is converted into a linear system of differential equations that easily solved numerically at each subsequent time level. Let m be a positive integer and subdivide the time interval I by the points tj , where tj = , τ = T/m , j = 0, 1, · · · , m . The suggested discretization scheme of problem (8) consists of the following problem (in the weak sense):
PPT Slide
Lager Image
PPT Slide
Lager Image
where ∂ x denotes to the derivative with respect to x and fj = f ( x , tj ).
The existence of a weak solution wj V is proved in [2] . By means of wj , ( j = 0, 1, · · · , m ) determined by the proposed scheme 10 -11 in each time step tj , we introduce the following piecewise linear functions (Rothe functions)
PPT Slide
Lager Image
and the corresponding step function
PPT Slide
Lager Image
Using the notation of Rothe function and its corresponding step function, a piecewise constant interpolation of equation 11 over I for each φ V yields
PPT Slide
Lager Image
where t j − 1 t tj and
PPT Slide
Lager Image
with
PPT Slide
Lager Image
In order to show the stability of the discrete solution and prove the convergence results, we shall derive some a priori estimates.
Lemma 3.1. Under the assumptions (H1)-(H3), there exists a positive constant C such that for any s ,
PPT Slide
Lager Image
Proof . Let us choose φ = τ δwj in 11 and summing over j from 1 to s , we obtain
PPT Slide
Lager Image
Taking into consideration (4), the second and the fourth terms are estimated by
PPT Slide
Lager Image
PPT Slide
Lager Image
With the aid of 11, Younge’s inequality (5) and (9), we estimate the third term by
PPT Slide
Lager Image
The last term is estimated by the use of (6) as follows
PPT Slide
Lager Image
Collecting (16)-(20), using ( H3 ), and choosing ε sufficiently small, we get
PPT Slide
Lager Image
Applying Gronwall’s inequality, we conclude the proof.
Lemma 3.2. Uniformly with respect to n one has
PPT Slide
Lager Image
PPT Slide
Lager Image
Proof . The estimates (22) and (23) are a consequence of (15). By the use of the identity
PPT Slide
Lager Image
the estimates (22) 3 , and (23) are a consequence of (15) 2 , (15) 3 and the definitions of wm and
PPT Slide
Lager Image
and thus the proof completes.
4. Convergence results
This section is devoted to proving the convergence of the proposed scheme and estimating its accuracy. Before we are able to prove convergence we need to prove the compactness of
PPT Slide
Lager Image
in L 2 ( I ; L 2 (Ω)) which is a consequence of the following assertion.
Lemma 4.1. The estimate
PPT Slide
Lager Image
holds uniformly for 0 < z < z 0 and n .
Proof . [3] We sum up (12) for i = s +1, · · · , s + k considering φ = ws + k ws ) τ . Then we sum it up for s = 1, · · · , m k and obtain the estimate
PPT Slide
Lager Image
Hence for z ≤ ( k + 1) τ we conclude the desired estimate.
Theorem 4.1. There exists u L 2 ( I ; L 2 (Ω)) ∩ H 1 ( I ; L 2 (Ω)) such that
PPT Slide
Lager Image
(in the sense of subsequences). Moreover, we have
PPT Slide
Lager Image
Proof . The estimate (22) 2 implies
PPT Slide
Lager Image
Hence, from lemma 4.1,
PPT Slide
Lager Image
is compact in L 2 ( I ; L 2 (Ω)) because of Kolmogorov’s compactness argument. So we can conclude that wn u and
PPT Slide
Lager Image
in L 2 ( I ; L 2 (Ω)) and also pointwise in Q . Also by the fact that
PPT Slide
Lager Image
we obtain that
PPT Slide
Lager Image
is weakly convergent in L 2 ( I ; V ); i.e.
PPT Slide
Lager Image
and this implies that
PPT Slide
Lager Image
It remains to prove that ∂ twm = ∂ tu . For each t I , by lemma 3.3, ∂ twm is uniformly bounded in the reflexive Banach space L 2 ( I ; L 2 (Ω)) and hence has a subsequence which converges weakly to an element ξ V (Eberlin-Smulian theorem [6] ). Thus
PPT Slide
Lager Image
Using Fubini theorem, we get
PPT Slide
Lager Image
Taking the limit as n → ∞, we obtain
PPT Slide
Lager Image
which implies
PPT Slide
Lager Image
Therefore, we have ∂ tu = ξ . Using the above discussions and the fact that and hence
PPT Slide
Lager Image
and passing with n → ∞ in (14) the proof is complete.
To obtain an error estimate, let us start by introducing the following additional notation
PPT Slide
Lager Image
We now take the difference between (8) and (14) we easily obtain the equality for every φ V
PPT Slide
Lager Image
Take φ = eu and write the new equation as I + II + III + IV = V . To begin with, we split the parabolic term into two parts I = I 1 + I 2 , where
PPT Slide
Lager Image
PPT Slide
Lager Image
The second term of equation (37) is bounded by
PPT Slide
Lager Image
The third and fourth terms of equation (37) are bounded by
PPT Slide
Lager Image
PPT Slide
Lager Image
In view of ( H1 ) and Younge’s inequality we can estimate the right hand by
PPT Slide
Lager Image
Collecting all the previous bounds, choosing ξ , ε , and κ sufficiently small and applying Gronwall lemma we obtain (28).
5. A Full discretization scheme
In this section, we use the exponential spline for the solution of discretization of problem (1-3) which is written in the form (denoting y ( x ) ≡ wj ( x )):
PPT Slide
Lager Image
with boundary conditions
PPT Slide
Lager Image
in the time step tj , j = 1, 2, · · · , n , wj −1 is known and we have to solve the two parameter singularly perturbed boundary value problem
PPT Slide
Lager Image
associated with the boundary condition (43). Here we use G ( x ) and H ( x ) to denote ( q ( x ) + 1/ τ ) and ( fj + wj−1 / τ ) respectively. To this end we consider a uniform mesh Δ with nodal point xi on [ a , b ] such that Δ : a = x 1 < x 2 < . . . < xm −1 < xm = b where
PPT Slide
Lager Image
Let y ( x ) be the exact solution of the problem (44)-(43) and Si be an approximation to yi = y ( xi ) obtained by the segment Qi ( x ) passing through the points ( xi , Si ) and ( xi +1 , Si +1 ). Each mixed spline segment Qi ( x ) has the form, (for more details see [15 - 18] ):
PPT Slide
Lager Image
where ai , bi , ci and di are constants and k is a free parameter.
Following [18] , we get the following relation for i = 1, 2, · · · , m − 1,
PPT Slide
Lager Image
where:
θ = kh , α = (sinh( θ )− θ )/ θ 2 sinh( θ ) and β = (2 θ cosh( θ )−2 sinh( θ ))/ θ 2 sinh( θ ), when k → 0 that θ → 0 then
PPT Slide
Lager Image
and the relation defined by (47) reduces to an ordinary cubic spline relation:
PPT Slide
Lager Image
at the point xi the proposed singularly perturbed problem may discretized by:
PPT Slide
Lager Image
where,
PPT Slide
Lager Image
Gi = G ( xi ), Hi = H ( xi ) and pi = p ( xi ).
Substituting (49) into (47), we get the following linear equation for every i = 1, 2, · · · , m − 1
PPT Slide
Lager Image
Equation (50) gives n −1 linear algebraic equations in n −1 unknowns Si ,where Ai = α pi −1 β pi −3 α pi +1 , Di = 3 α pi −1 + β pi α pi +1 , Ei = 4 α ( pi +1 pi −1 ).
Now we investigate the convergence analysis of the suggested algorithm. The exponential spline solution of (1)-(3) is based on the linear equation given by (50) Let Y = ( yi ), S = ( si ), C = ( ci ), T = ( ti ) and E = ( ei ) = Y S be m − 1 dimensional column vectors. It is easily seen that the system of equation (50) gives n − 1 linear algebraic equations in the n − 1 unknowns Si , i = 1, 2, · · · , m − 1. Then we can rewrite equation (50) in the standard matrix equation as:
PPT Slide
Lager Image
Now the matrix B may be written as
PPT Slide
Lager Image
where
PPT Slide
Lager Image
and
PPT Slide
Lager Image
We have the local truncation error are
PPT Slide
Lager Image
From Eq. (52), then the row sum of the matrix B satisfies: for | pi | ≤ P ≥ 0,
PPT Slide
Lager Image
for small values of εd and εc then the matrix B is irreducible and monotone and it follows that B −1 exist thus the problem (51) has a unique solution.
Lemma 5.1 ( [5] ). For any 0 < d < 1, we have up to a certain order q that it depends on the smoothness of the data
PPT Slide
Lager Image
- Shishkin Mesh Strategy
We form the piecewise-uniform grid in such a way that more points are generated in the boundary layer regions than outside of it. We divide the interval
PPT Slide
Lager Image
into three subintervals:
PPT Slide
Lager Image
where transition parameters are given by
PPT Slide
Lager Image
where
PPT Slide
Lager Image
the real solutions of characteristic equation
PPT Slide
Lager Image
with n to be the number of subdivision points of the interval [0,1] and we place n /4, n /2 and n /4 mesh points, respectively, in [0, γ 1 ], [ γ 1 , 1− γ 2 ] and [1− γ 2 , 1]. Denote the step sizes in each subinterval by
PPT Slide
Lager Image
and
PPT Slide
Lager Image
respectively.
Where the resulting piecewise-uniform Shishkin mesh may be represented by:
PPT Slide
Lager Image
where n is the number of discretization points and set of mesh points
PPT Slide
Lager Image
with:
PPT Slide
Lager Image
Thus, a uniform mesh is placed on each of these subintervals.
Lemma 5.2 ( [5] ). The solution y ( x ) of (42)-(43) has the representation
PPT Slide
Lager Image
Theorem 5.1. Let S ( x ) be the approximate solution of the solution y ( x ) of (42), then
PPT Slide
Lager Image
where C is a constant independent of εd and εc .
Proof . The estimate is obtained on each subinterval Ω i = [0, 1] separately. Let any function z on Ω i , see [5]
PPT Slide
Lager Image
and so it is obvious that, on Ω i
PPT Slide
Lager Image
Taking maximum norm on both sides, we get
PPT Slide
Lager Image
and by appropriate Taylor expansions it is easy to see that
PPT Slide
Lager Image
From (58) and Lemma 5.1, on Ω i ,
PPT Slide
Lager Image
Also, using Lemma 5.3 and equation (58) on Ω i ,
PPT Slide
Lager Image
Case 1. The argument now depends on whether
PPT Slide
Lager Image
In this case ψ 1 C ln n and ψ 2 C ln n . Then the result follows at once from (58)
PPT Slide
Lager Image
Case 2. When
PPT Slide
Lager Image
ln n and
PPT Slide
Lager Image
ln n . Suppose that i satisfies
PPT Slide
Lager Image
and
PPT Slide
Lager Image
Then
PPT Slide
Lager Image
and
PPT Slide
Lager Image
respectively. Now from (58)
PPT Slide
Lager Image
If i satisfied
PPT Slide
Lager Image
Then γ 1 xi and xi ≤ 1 − γ 2 or γ 2 ≤ 1 − xi and so
PPT Slide
Lager Image
PPT Slide
Lager Image
Using the above values in (59), we get the required result.
We are now in a position to formulate the main contribution of this work by combining (28) and (59) in the following theorem.
Theorem 5.2. Let the assumptions (H1)-(H2) be satisfied . Then uα u in u H 1 ( I ; L 2 (Ω))∩ L 2 ( I ; V ) for α → 0 where u is the variational solution of problem (1)-(3) and uα is constructed from the solution of problem (42)-(43). Moreover, we have the error estimates
PPT Slide
Lager Image
6. Numerical Results
We present a case study to illustrate the performance of the error estimation procedures of section 5. For purposes of comparison, contrast and performance, two cases with known solutions were chosen.
Case study 1: Consider the following singularly perturbed parabolic partial differential equation, see [1] .
ut εduxx + ux = f ( x , t ), T = 1,in the interval [0, 1].
associated with the Dirichlet boundary conditions u (0, t ) = 0, u (1, t ) = 0. and the initial condition u 0 ( x ) = e −1/ ε +(1− e −1/ ε ) x e −(1− x )/ ε , where f ( x , t ) = e t (− c 1 + c 2 (1 − x ) + e −(1− x )/ ε ) . The exact solution is given by u ( x , t ) = e t ( c 1 + c 2 x e −(1− x )/ ε ), where c 1 = e −1/ ε , c 2 = 1 − e −1/ ε .
Case study 2: Consider the following singularly perturbed parabolic partial differential equation.
ut εduxx + εcux = f ( x , t ), T = 1,in the interval [0, 1].
associated with the Dirichlet boundary conditions u (0, t ) = 0, u (1, t ) = 0. and the initial condition
PPT Slide
Lager Image
where
PPT Slide
Lager Image
Maximum pointwise errors for Case study1
PPT Slide
Lager Image
Maximum pointwise errors for Case study1
Maximum pointwise errors and order of convergence for Case study1
PPT Slide
Lager Image
Maximum pointwise errors and order of convergence for Case study1
The exact solution is given by
PPT Slide
Lager Image
where
PPT Slide
Lager Image
All the computations were performed by using MATLAB 7.
Comparison of maximum pointwise errors for Case study1
PPT Slide
Lager Image
Comparison of maximum pointwise errors for Case study1
Maximum errors for Case study2n= 128 andτ= 1.6/(n+ 1)2
PPT Slide
Lager Image
Maximum errors for Case study2 n = 128 and τ = 1.6/(n + 1)2
To compute the experimental rates of convergence Ordn for every fixed ε , we use the rate of convergence from:
PPT Slide
Lager Image
The numerical results displayed in Table 1 - 3 clearly indicate that the proposed method based on an exponential-spline function with Shishkin mesh is ε -uniformly convergent. Numerical solution profiles for case study are given in figure 1 , respectively, for n = 128, τ = 0.02and different values of εd = 1, 2 −6 and 2 −12 . For comparison, our results are better than the numerical solution presented in Clavero et al. [1] . Numerical solution profiles for case study2 are given in table 4 and also figure 2 and 3 for different values of εd and εc .
It is clearly that all solutions are in good harmony with each approximate and exact solution and the convergence obtained is parameter uniform and satisfies the theoretical predictions.
PPT Slide
Lager Image
Exact and approximate solutions profiles of the case study1 for n = 128, τ = 0.02 and εd = 2−6and εd = 2−12 respectively
PPT Slide
Lager Image
Exact and approximate solutions profiles of the case study2 for n = 128, 256, τ = 1.6/(n + 1)2, εc = 2−6, and εd = 2−6.
PPT Slide
Lager Image
Exact and approximate solutions profiles of the case study2 for n = 64, 128, τ = 0.002, εc = 2−6, and εd = 2−6.
7. Conclusion
In this paper, we have investigated the application of Rothe’s approximation in time discretization and the exponential spline approximation in the spatial discretization to solve the convection–diffusion problems. The method is shown to be uniformly convergent i.e., independent of mesh parameters and perturbation parameters εd and εc . It has also been seen that the accuracy in the numerical results for the proposed scheme is comparable to that obtained by the uniformly convergent method with a uniform mesh [1] . It has been found that the proposed algorithm gives highly accurate numerical results and higher order of convergence than an upwind finite difference scheme. The present method provides a very simple and accurate method for solving convection–diffusion problems.
BIO
W. K. Zahra received M.Sc. and Ph.D. from Faculty of Engineering, Tanta University, Tanta, Egypt. He is currently working as an Assistant Professor at Faculty of Engineering, Tanta university,Tanta, Egypt. His research interests include Numerical Methods for Singularly Perturbed Differential Equations.
Department of Physics and Engineering Mathematics, Faculty of Engineering, Tanta Univ., Tanta, Egypt.
e-mail: waheed zahra@yahoo.com, wzahra@f-eng.tanta.edu.eg
M. S. El-Azab received M.Sc. from Faculty of Engineering, AlMansurah Univ., Al Mansurah, Egypt and Ph.D. from Comenius University, Slovakia. He is currently working as a Professor at Faculty of Engineering, AlMansurah Univ., Al Mansurah, Egypt. His research interests include Numerical Methods for Singularly Perturbed Differential Equations.
Department of Mathematics, Faculty of Engineering, AlMansurah Univ., Al Mansurah, Egypt.
Ashraf M. El Mhlawy received M.Sc. from Faculty of Engineering, Tanta University, Tanta, Egypt. His research interests include Numerical Methods for Singularly Perturbed Differential Equations.
Department of Physics and Engineering Mathematics, Faculty of Engineering, Tanta Univ., Tanta, Egypt.
References
Clavero C. , Jorge J.C. , Lisbona F. 2003 Uniformly convergent scheme on a nonuniform mesh for convection-diffusion parabolic problems J. Comput. Appl. Math. 154 415 - 429    DOI : 10.1016/S0377-0427(02)00861-0
Amin N.G. , El-Azab M.S. 2012 A new linear scheme for solving the 1 D Burger’s equation American Academic and Scholarly Research Journal 4 (2)
Jager W. , Kacur J. 1990 ”Approximation of porous medium type systems by nondegenerateelliptic systems” Preprint No. 582 Universitat Heidelberg
Kadalbajoo M.K. , Gupta V. , Awasthi A. 2008 A uniformly convergent B-Spline collocation method on a nonuniform mesh for singularly perturbed one-dimensional time-dependent linear convection-diffusion problem Journal of Computational and Applied Mathematics 220 271 - 289    DOI : 10.1016/j.cam.2007.08.016
Journal of Computational and Applied Mathematics M.K. , Yadaw A.S 2008 B-Spline collocation method for a two-parameter singularly perturbed convection-diffusion boundary value problem. Applied Mathematics and Computation 201 504 - 513    DOI : 10.1016/j.amc.2007.12.038
Kufner A. , John O. , Fucik S. 1997 ”Function Spaces” Nordhoff Leyden
Ramos J.I. 2005 An exponentially fitted method for singularly perturbed, one-dimensional, parabolic problems Appl. Math. Comput. 161 513 - 523    DOI : 10.1016/j.amc.2003.12.046
Rashidinia J. , Ghasemi M. , Jalilian R. 2010 A collocation method for the solution of nonlinearone-dimensional parabolic equations Mathematical Sciences 4 (1) 87 - 104
Rao S.C.S , Kumar M. 2008 Exponential B-spline collocation method for self -adjoint singularly perturbed boundary value problems Applied Numerical Mathematics 58 1572 - 1581    DOI : 10.1016/j.apnum.2007.09.008
Rao S.C.S , Kumar M. 2009 Parameter uniformly exponential spline difference scheme for singularly perturbed semilinear reaction diffusion problems Nonlinear Analysis 71 1579 - 1588    DOI : 10.1016/j.na.2009.01.198
Sakai M. , Usmani R.A. 1986 On exponential splines J. approx. theory 47 122 - 131    DOI : 10.1016/0021-9045(86)90037-7
Sakai M. , Usmani R.A. 1989 A class of simple exponential B-splines and their applications to numerical solution to singular perturbation problems Numer. Math. 55 493 - 500    DOI : 10.1007/BF01398912
Surla K. , Jerkovic V. 1985 Some possibilities of applying spline collocations to singular perturbation problems Numerical Methods and Approximation Theory Novisad II 19 - 25
Sharma M. , Kaushik A. 2012 Convergence Analysis of Weighted Difference Approximations on Piecewise Uniform Grids to a Class of Singularly Perturbed Functional Differential Equations J Optim Theory Appl 155 252 - 272    DOI : 10.1007/s10957-011-9965-5
Zahra W.K. 2009 Exponential spline solutions for a class of two point boundary value problems over a semi-infinite range Numer. Algor. 52 561 - 573    DOI : 10.1007/s11075-009-9298-8
Zahra W.K. 2011 A smooth approximation based on exponential spline solutions for nonlinear fourth order two point boundary value problems Applied Mathematics andComputation 217 8447 - 8457    DOI : 10.1016/j.amc.2011.03.043
Zahra W.K. 2011 Finite-difference technique based on exponential splines for the solutionof obstacle problems International Journal of Computer Mathematics 88 (14) 3046 - 3060    DOI : 10.1080/00207160.2011.573846
Zahra W.K. , El Mhlawy A.M. (2013) Numerical Solution of Two-Parameter Singularly Perturbed Boundary Value Problems via Exponential Spline Journal of King Saud University (Science)    DOI : 10.1016/j.jksus.2013.01.003.