Research Article  Open Access
Musa Çakır, Gabil M. Amiraliyev, "A Numerical Method for a Singularly Perturbed ThreePoint Boundary Value Problem", Journal of Applied Mathematics, vol. 2010, Article ID 495184, 17 pages, 2010. https://doi.org/10.1155/2010/495184
A Numerical Method for a Singularly Perturbed ThreePoint Boundary Value Problem
Abstract
The purpose of this paper is to present a uniform finite difference method for numerical solution of nonlinear singularly perturbed convectiondiffusion problem with nonlocal and third type boundary conditions. The numerical method is constructed on piecewise uniform Shishkin type mesh. The method is shown to be convergent, uniformly in the diffusion parameter , of first order in the discrete maximum norm. Some numerical experiments illustrate in practice the result of convergence proved theoretically.
1. Introduction
This paper is concerned with uniform numerical method for the singularly perturbed semilinear boundaryvalue problem (BVP): where is a small positive parameter, the functions , and are sufficiently smooth on , and , respectively, and furthermore The solution generally has boundary layers near and .
Singularly perturbed differential equations are characterized by the presence of a small parameter multiplying the highestorder derivatives. Such problems arise in many areas of applied mathematics. Among these are the NavierStokes equations of fluid flow at high Reynolds number, mathematical models of liquid crystal materials and chemical reactions, control theory, reactiondiffusion processes, quantum mechanics, and electrical networks. The solutions of singularly perturbed differential equations typically have steep gradients, in thin regions of the domain, whose magnitude depends inversely on some positive power of . Such regions are called either interior or boundary layers, depending on whether their location is in the interior or at the boundary of the domain. An overview of some existence and uniqueness results and applications of singularly perturbed equations can be found in [1–6].
It is known that these problems depend on a small positive parameter in such away that the solution exhibits a multiscale character; that is, there are thin transition layers where the solution varies rapidly, while away from layers it behaves regularly and varies slowly. The treatment of singularly perturbed problems presents severe difficulties that have to be addressed to ensure accurate numerical solutions. Therefore it is important to develop suitable numerical methods for solving these problems, whose accuracy does not depend on the value of parameter , that is, methods that are convergent uniformly. These include fitted finitedifference methods, finite element methods using special elements such as exponential elements, and methods which use a priori refined or special piecewise uniform grids which condense in the boundary layers in a special manner. One of the simplest ways to derive parameteruniform methods consists of using a class of special piecewise uniform meshes, such as Shishkin meshes (see [3, 6–15] for the motivation for this type of mesh), which are constructed a priori and depend on the parameter , the problem data, and the number of corresponding mesh points. The various approaches to numerical solution of differential equations with stepwise continuous solutions can be found in [2, 3, 6].
There is also an increasing interest in the application of Shishkin meshes to singularly perturbed convectiondiffusion problems (see [16, 17] and references cited therein). However, much of the Shishkin mesh literature is concerned with linear or quasilinear singularly perturbed twopoint problems with firstorder reduced equations. In [18] has been obtained a result uniform for the twopoint boundary value problem of (1.1), by using a fitted operator method on uniform meshes.
In the present paper, we analyse a fitted finitedifference scheme on a Shishkin type mesh for the numerical solution of the semilinear nonlocal boundary value problem (1.1)–(1.3). The origin of the fitted finite difference method can be traced to [19]; for subsequent work on fitted operator method and its applications, see [2, 3]. Nonlocal boundary value problems have also been studied extensively in the literature. For a discussion of existence and uniqueness results and for applications of nonlocal problems see [20–26] and the references cited therein. Some approaches to approximating this type of problem have also been considered [20, 21, 26–28]. However, the algorithms developed in the papers cited above are mainly concerned with regular cases (i.e., when boundary layers are absent). In [27] has been studied the fitted difference schemes on an equidistant mesh for the numerical solution of the linear threepoint reactiondiffusion problem.
The numerical method presented here comprises a fitted difference scheme on a piecewise uniform mesh. We have derived this approach on the basis of the method of integral identities using interpolating quadrature rules with the weight and remainder terms in integral form. This results in a local truncation error containing only firstorder derivatives of exact solution and hence facilitates examination of the convergence. A summary of paper is as follows. Section 2 contains results concerning the exact solutions of problem (1.1)–(1.3). In Section 3, we describe the finitedifference discretization and construct a piecewise uniform mesh, which is fitted to the boundary layers. In Section 4, we present the error analysis for the approximate solution. Uniform convergence is proved in the discrete maximum norm. In the following section numerical results are presented, which are in agreement with the theoretical results. The approach to the construction of the discrete problem and the error analysis for the approximate solution are similar to those in [18, 27–29].
2. Continuous Solution
In this section, we give uniform bounds for the solution of the BVP (1.1)–(1.3), which will be used to analyze properties of appropriate difference problem.
Lemma 2.1. Let . Then the solution of problem (1.1)–(1.3) satisfies the inequalities where with providing that is bounded for and .
Proof. We rewrite the problem (1.1)–(1.3) as
where
Here we use the Maximum Principle: let and be the differential operators in (2.5)(2.6) and . If and for all, then for all . Then, from (2.5)(2.6) we have the following inequality:
Next, from boundary condition (2.7), we get
By setting the value in the inequality (2.9), we obtain
From (2.10) and (2.11), then we have
which along with (2.9) leads to (2.1).
After establishing (2.1) the further part of the proof is almost identical to that of [28].
3. Discretization and Mesh
Let be any nonuniform mesh on : and . For each we set the stepsize . Before describing our numerical method, we introduce some notation for the mesh functions. For any mesh function defined on we use
The difference scheme we will construct follows from the identity with the basis functions having the form where and , respectively, are the solutions of the following problems: The functions and can be explicitly expressed as follows: The coefficient in (3.3) is given by Rearranging (3.3) gives with As consistent with [26, 27], we obtain the precise relation: where
It then follows from (3.8) that where and after a simple calculation
To define an approximation for the boundary condition (1.2), we proceed our discretization process by with Here the function is the solution of the following problem: In the analogous way, as in construction of (3.12), we obtain where
Now, it remains to define an approximation for the second boundary condition (1.3). Let be the mesh point nearest to . Then where being intermediate point between and .
Based on (3.12), (3.18), and (3.22), we propose the following difference scheme for approximating (1.1) and (1.3): where and are defined by (3.14), (3.19), and (3.20), respectively.
The difference scheme (3.24)–(3.26) in order to be uniform convergent, we will use the Shishkin mesh. For a divisible by 4 positive integer , we divide each of the intervals and into equidistant subintervals and also into equidistant subintervals, where the transition points and , which separate the fine and coarse portions of the mesh, are obtained by taking where and are given in Lemma 2.1. In practice one usually has ; so the mesh is fine on , and coarse on . Hence, if we denote by , and the stepsize in , and , respectively, we have and so In the rest of the paper we only consider this mesh.
We note that on this mesh the coefficient which is defined by (3.14) simplifies to For the evaluation of the rest values and we will use the form (3.14).
4. Error Analysis
Let . Then for the error of the difference scheme (3.24)–(3.26) we get where the truncation errors , and are defined by (3.9), (3.21), and (3.23), respectively.
Lemma 4.1. The solution of problem (4.1)–(4.3) satisfies
Proof. The problem (4.1)–(4.3) can be rewritten as
where
intermediate points called for by the mean value theorem.
Since the discrete maximum principle is valid here, we have the proof of (4.4) by analogy with the proof of Lemma 2.1.
Lemma 4.2. Under the above assumptions of Section 1 and Lemma 2.1, for the error functions , and , the following estimates hold:
Proof. From explicit expression (3.9) for , on an arbitrary mesh we have
This inequality together with (2.3) enables us to write
in which
We consider first the case , and so , and . Hereby, from (4.9) we can write
Hence
We now consider the case and , and so and estimate on , and separately. In the layer region , the inequality (4.9) reduces to
Hence
The same estimate is obtained in the layer region in the similar manner.
It remains to estimate for . In this case we are able to write (4.9) as
Since , it follows that
Also, if we rewrite the mesh points in the form , evidently
The last two inequalities together with (4.15) give the bound
It remains to estimate for the mesh points and . For the mesh point , inequality (4.9) reduces to
Since
it then follows that
The same estimate is obtained for in the similar manner.
The same estimate is valid when only one of the values and is equal to .
Next, we estimate the remainder term . From the explicit expression (3.21), taking into consideration that and , we obtain
Hence
Finally, we estimate the remainder term . From the expression (3.23) we obtain
where we have assumed that is on the lefthand side of (if is on right side of , the integral will be over ()). In the same manner as above we therefore obtain from here that
Thus Lemma 4.2 is proved.
Combining the two previous lemmas gives us the following convergence result.
Theorem 4.3. Let be the solution of (1.1)–(1.3) and the solution (3.24)–(3.26). Then
5. Numerical Results
In this section, we present some numerical results which illustrate the present method.
Example 5.1. Consider the test problem:
The exact solution of our test problem is unknown. Therefore, we use a doublemesh method [2] to estimate the errors and compute the experimental rates of convergence in our computed solutions; that is, we compare the computed solution with the solution on a mesh that is twice as fine (for details see [13, 28]). The error estimates obtained in this way are denoted by
where is the approximate solution of the respective method on the mesh
with
The convergence rates are
Approximations to the uniform rates of convergence are estimated from
The corresponding uniform convergence rates are computed using the formula
To solve the nonlinear problem (3.24)–(3.26) we use the following iteration technique:
To solve (5.8)–(5.10), we take the initial approximation as and the stopping criterion is
The computed maximum pointwise errors and , and the orders of uniform convergence for different values of and , based on the doublemesh principle are presented in Tables 1 and 2. The results established here are that the discrete solution is uniformly convergent with respect to the perturbation parameter and the errors are uniformly convergent with rates of almost unity as predicted by our theoretical analysis.


Example 5.2. Consider the test problem: The exact solution of our test problem is unknown. In the same manner as above we solve this problem.
Acknowledgment
The authors wish to thank the referees for their suggestions and comments which helped improve the quality of manuscript.
References
 E. M. Dejager and J. Furu, The Theory of Singular Perturbations, Elsevier Science, 1996.
 E. P. Doolan, J. J. H. Miller, and W. H. A. Schilders, Uniform Numerical Methods for Problems with Initial and Boundary Layers, Boole Press, Dublin, Ireland, 1980. View at: Zentralblatt MATH  MathSciNet
 P. A. Farrell, A. F. Hegarty, J. J. H. Miller, E. O'Riordan, and G. I. Shishkin, Robust Computational Techniques for Boundary Layers, vol. 16 of Applied Mathematics, Chapman & Hall/CRC, Boca Raton, Fla, USA, 2000. View at: Zentralblatt MATH  MathSciNet
 A. H. Nayfeh, Introduction to Perturbation Techniques, John Wiley & Sons, New York, NY, USA, 1993. View at: Zentralblatt MATH  MathSciNet
 R. E. O'Malley Jr., Singular Perturbation Methods for Ordinary Differential Equations, vol. 89 of Applied Mathematical Sciences, Springer, New York, NY, USA, 1991. View at: Zentralblatt MATH  MathSciNet
 H.G. Roos, M. Stynes, and L. Tobiska, Robust Numerical Methods for Singularly Perturbed Differential Equations, vol. 24 of Springer Series in Computational Mathematics, Springer, Berlin, Germany, 2nd edition, 2008. View at: Zentralblatt MATH  MathSciNet
 P. W. Hemker, A numerical study of stiff twopoint boundary value problems, Ph.D. thesis, Mathematical Center MCT 80, Amsterdam, The Netherlands, 1977.
 C. Clavero, J. L. Gracia, and F. Lisbona, “High order methods on Shishkin meshes for singular perturbation problems of convectiondiffusion type,” Numerical Algorithms, vol. 22, no. 1, pp. 73–97, 1999. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 P. A. Farrell, E. O'Riordan, and G. I. Shishkin, “A class of singularly perturbed semilinear differential equations with interior layers,” Mathematics of Computation, vol. 74, no. 252, pp. 1759–1776, 2005. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J. L. Gracia, E. O'Riordan, and M. L. Pickett, “A parameter robust second order numerical method for a singularly perturbed twoparameter problem,” Applied Numerical Mathematics, vol. 56, no. 7, pp. 962–980, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 E. O'Riordan, M. L. Pickett, and G. I. Shishkin, “Parameteruniform finite difference schemes for singularly perturbed parabolic diffusionconvectionreaction problems,” Mathematics of Computation, vol. 75, no. 255, pp. 1135–1154, 2006. View at: Publisher Site  Google Scholar  MathSciNet
 H.G. Roos, “Ten ways to generate the Il'in and related schemes,” Journal of Computational and Applied Mathematics, vol. 53, no. 1, pp. 43–59, 1994. View at: Publisher Site  Google Scholar  MathSciNet
 T. Linß and M. Stynes, “A hybrid difference scheme on a Shishkin mesh for linear convectiondiffusion problems,” Applied Numerical Mathematics, vol. 31, no. 3, pp. 255–270, 1999. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 I. A. Savin, “On the rate of convergence, uniform with respect to a small parameter, of a difference scheme for an ordinary differential equation,” Computational Mathematics and Mathematical Physics, vol. 35, no. 11, pp. 1417–1422, 1995. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 G. F. Sun and M. Stynes, “A uniformly convergent method for a singularly perturbed semilinear reactiondiffusion problem with multiple solutions,” Mathematics of Computation, vol. 65, no. 215, pp. 1085–1109, 1996. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 T. Linß, “Layeradapted meshes for convectiondiffusion problems,” Computer Methods in Applied Mechanics and Engineering, vol. 192, no. 910, pp. 1061–1105, 2003. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 T. Linß and H.G. Roos, “Analysis of a finitedifference scheme for a singularly perturbed problem with two small parameters,” Journal of Mathematical Analysis and Applications, vol. 289, no. 2, pp. 355–366, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 G. Amiraliyev and M. Çakır, “A uniformly convergent difference scheme for a singularly perturbed problem with convective term and zeroth order reduced equation,” International Journal of Applied Mathematics, vol. 2, no. 12, pp. 1407–1419, 2000. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 A. M. Il'in, “A difference scheme for a differential equation with a small parameter multiplying the highest derivative,” Mathematical Notes of the Academy of Science of the USSR, vol. 6, pp. 237–248, 1969. View at: Google Scholar  MathSciNet
 R. Čiegis, “The numerical of singularly perturbed non local problem,” Lietuvas Matematica Rink, vol. 28, pp. 144–152, 1988 (Russian). View at: Google Scholar
 R. Čiegis, “On the difference schemes for problems with nonlocal boundary conditions,” Informatica (Lietuva), vol. 2, no. 2, pp. 155–170, 1991. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 C. P. Gupta and S. I. Trofimchuk, “A sharper condition for the solvability of a threepoint second order boundary value problem,” Journal of Mathematical Analysis and Applications, vol. 205, no. 2, pp. 586–597, 1997. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 B. Liu, “Positive solutions of secondorder threepoint boundary value problems with change of sign,” Computers & Mathematics with Applications, vol. 47, no. 89, pp. 1351–1361, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 R. Ma, “Positive solutions for nonhomogeneous mpoint boundary value problems,” Computers & Mathematics with Applications, vol. 47, no. 45, pp. 689–698, 2004. View at: Publisher Site  Google Scholar  MathSciNet
 A. M. Nahushev, “On non local boundary value problems,” Differential Equations, vol. 21, pp. 92–101, 1985 (Russian). View at: Google Scholar
 M. P. Sapagovas and R. Yu. Čiegis, “Numerical solution of some nonlocal problems,” Lietuvas Matematica Rink, vol. 27, no. 2, pp. 348–356, 1987 (Russian). View at: Google Scholar  MathSciNet
 G. M. Amiraliyev and M. Çakır, “Numerical solution of the singularly perturbed problem with nonlocal boundary condition,” Applied Mathematics and Mechanics, vol. 23, no. 7, pp. 755–764, 2002. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. Çakır and G. M. Amiraliyev, “Numerical solution of a singularly perturbed threepoint boundary value problem,” International Journal of Computer Mathematics, vol. 84, no. 10, pp. 1465–1481, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 G. M. Amiraliyev and H. Duru, “A uniformly convergent difference method for the periodical boundary value problem,” Computers & Mathematics with Applications, vol. 46, no. 56, pp. 695–703, 2003. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
Copyright
Copyright © 2010 Musa Çakır and Gabil M. Amiraliyev. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.