Research Article  Open Access
D. Kumara Swamy, K. Phaneendra, Y. N. Reddy, "Solution of Singularly Perturbed DifferentialDifference Equations with Mixed Shifts Using Galerkin Method with Exponential Fitting", Chinese Journal of Mathematics, vol. 2016, Article ID 1935853, 10 pages, 2016. https://doi.org/10.1155/2016/1935853
Solution of Singularly Perturbed DifferentialDifference Equations with Mixed Shifts Using Galerkin Method with Exponential Fitting
Abstract
Galerkin method is presented to solve singularly perturbed differentialdifference equations with delay and advanced shifts using fitting factor. In the numerical treatment of such type of problems, Taylor’s approximation is used to tackle the terms containing small shifts. A fitting factor in the Galerkin scheme is introduced which takes care of the rapid changes that occur in the boundary layer. This fitting factor is obtained from the asymptotic solution of singular perturbations. Thomas algorithm is used to solve the tridiagonal system of the fitted Galerkin method. The method is analysed for convergence. Several numerical examples are solved and compared to demonstrate the applicability of the method. Graphs are plotted for the solutions of these problems to illustrate the effect of small shifts on the boundary layer solution.
1. Introduction
Singularly perturbed differentialdifference equations (SPDDEs) arise very frequently in the mathematical modelling of real life situations in science and engineering [1–3]. In the mathematical modelling of a physical system as in control theory, the presence of small time parasitic parameters like moments of inertia, resistances, inductances, and capacitances increases the order and stiffness of these systems. The suppression of these small constants results in the reduction of the order of the system. Such systems are termed as singular perturbation systems and when these systems take into account the past history as well as the present state of the physical system then they are called singularly perturbed delay differential equations. Delay differential equations arise in firstexit time problems in neurobiology and in mathematical formulation of various practical phenomena in biosciences. A differentialdifference equation with the presence of shift terms induces large amplitudes and exhibits oscillations, resonance, turning point behaviour, and boundary and interior layers. Hence, to control such behaviour, we need some simple and efficient numerical techniques.
Lange and Miura [3–7] published a series of papers extending the method of matched asymptotic expansions initially developed for ordinary differential equations to obtain approximate solution of singularly perturbed differentialdifference equations.
Numerical analysis of singularly perturbed differentialdifference turning point problems was initiated by Kadalbajoo and Sharma. In a series of papers, [8–10], they gave many robust numerical techniques for the solution of such type of problems. Kadalbajoo and Sharma [8] elucidate a numerical method to solve boundary value problems for singularly perturbed differentialdifference equation with mixed shifts. Kadalbajoo and Sharma [9] proposed a numerical method to solve boundary value problems for a singularly perturbed differentialdifference equation of a mixed type, that is, which contains both types of terms having negative shifts as well as positive shifts, and considered the case in which the solution of the problem exhibits rapid oscillations. Kadalbajoo and Sharma [10] described a numerical approach based on finite difference method to solve a mathematical model arising from a model of neuronal variability. Kadalbajoo and Kumar [11] used Bspline collocation method with fitted mesh for the solution of singularly perturbed differentialdifference equations with small delay.
Patidar and Sharma [12] combined fittedoperator methods with Micken’s nonstandard finite difference techniques for the numerical approximations of singularly perturbed linear delay differential equations. Kadalbajoo et al. [13] derived uniformly convergent fitted methods for the solution of singularly perturbed differentialdifference equation (SPDDE). Kumar and Sharma [14] presented a numerical scheme based on Bspline collocation to approximate the solution of boundary value problems for singularly perturbed differentialdifference equations with delay and advance.
With this motivation, an exponentially fitting factor is introduced in Galerkin method for the solution of singularly perturbed differentialdifference equation with delay and advanced parameters. In Section 2, description of the problem is given. In Section 3, numerical scheme for the solution of the problem is presented and Section 4 deals with convergence analysis of the proposed scheme. To demonstrate the efficiency of the proposed method, numerical experiments are carried out for several test problems and the results are given in Section 5. Finally the conclusions are given in the last section.
2. Description of the Problem
Consider a linear singularly perturbed differentialdifference equation of the following form:on (0, 1), under the boundary conditionsHere is a small parameter such that , and are smooth functions, and are, respectively, the delay (negative shift) and the advance (positive shift) parameters. If , the solution of (1) with (2) exhibits layer at the left end of the interval and if , the layer exists at the right end of the interval.
Since the solution of boundary value problem equations (1) and (2) is sufficiently differentiable, we expand the terms and using Taylor series; we getUsing (3) and (4) in (1), we getEquation (5) is an asymptotically equivalent secondorder singular perturbation problem of (1) with boundary conditionsSince and , the transition from (1) to (5) is admitted. This replacement is significant from the computational point of view. For more details on the validity of this transition, one can refer El’sgol’ts and Norkin [15]. Thus, the solution of (5) provides a good approximation to the solution of (1).
Here,
3. Numerical Scheme
3.1. LeftEnd Boundary Layer Problem
Let be a decomposition of the considered interval into equal intervals with constant mesh length . Then we have the nodes , for . Assume that , and are sufficiently continuously differentiable functions in . If in where is a positive constant, (5) has a unique solution which, in general, displays a boundary layer of width at .
Lemma 1 (Doolan et al. [16] and O’Malley [17]). Let be the zerothorder asymptotic approximation to the solution of (5), where represents the zerothorder approximate outer solution (i.e., the solution of the reduced problem of (5)) and represents the zerothorder approximate solution in the boundary layer region of (5).
Then for a fixed positive integer ,
Proof. Let be the solution of the reduced problem of (5)and is the solution of the boundary value problem (cf. O’Malley [17])From the theory of singular perturbation, the zerothorder asymptotic approximation to the solution of (4) is (cf. O’Malley [17])As we are considering the differential equations on sufficiently small subintervals, the coefficients could be assumed to be locally constant. Hence,So, at the nodal points, we havethat is,Thereforewhere .
Now, we consider the difference scheme [18] by Galerkin method as follows:
Select a set of basis functions , which will define an interpolation scheme for the approximate solution over a grid of points . For simplicity, we use piecewise Lagrange polynomials of first degree as the basis functions. These interpolating polynomials arein local element coordinates .
The nodal values of the approximate solution at the interior nodes are determined using this basis. The given boundary conditions determine the value of at the end nodes and . The Galerkin method is now employed to obtain the integral equations; we havewhich is an integral equation .
Since is sum of piecewise linear Lagrange polynomials, the secondorder derivatives appearing in (17) vanish except at the element boundaries , where they become infinite.
By integration by parts, (18) becomesUsing the substitution of trial function into the integral equation (19), we havefor .
It can be observed that all quantities on the right side of (20) can be computed from known boundary data to obtain equations in the unknown values at the interior nodes.
The integrals in (20) can be solved by taking advantage of local coordinate system.
Sincewe have, by simple integration,By assuming , and as constants, the integral equation (20) gives, for a typical internal node ,Equation (23), when rearranged, gives the following system of difference equations:Now, introduce a fitting factor in the Galerkin scheme as follows:for with Here is a fitting factor which is to be determined in such a way that the solution of (25) converges uniformly to the solution of (5). Multiplying (25) by and taking the limit as (in [16]), we get Now, approximating the solution y(x) by zerothorder asymptotic approximation and using Lemma 1, we haveUsing the above equations in (27), we getFrom (25), we havefor .
Equation (30) can be written as a threeterm recurrence relation as follows:whereThe tridiagonal system equation (31) is solved using Thomas algorithm.
3.2. RightEnd Layer Problems
We now discuss the method for singularly perturbed twopoint boundary value problems with rightend boundary layer of the underlying interval. Assume that , and are sufficiently continuously differentiable functions in . Furthermore, assume that in , where is a negative constant. Under these assumptions, (5) has a unique solution which, in general, displays a boundary layer of width at .
Lemma 2. Let be the zerothorder asymptotic approximation to the solution of (5), where represents the zerothorder approximate outer solution and represents the zerothorder approximate solution in the boundary layer region.
Then for a fixed positive integer ,
Proof. The proof is based on asymptotic analysis (Doolan et al. [16] and O’Malley [17]) and is similar to the proof of Lemma 1.
Applying the same procedure as in Section 3 and using Lemma 2, we get the tridiagonal system equation (20) with fitting factor as
4. Convergence Analysis
Writing the tridiagonal system equation (31) in matrixvector form, we getin which , is a tridiagonal matrix of order , withand is a column vector with , where with local truncation errorand
We also havewhere denotes the actual solution and is the local truncation error.
From (35) and (38), we getThus, the error equation iswhere .
Let the th row elements sum of matrix be ; then we haveWe can choose sufficiently small so that the matrix is irreducible and monotone. It follows that exists and its elements are nonnegative.
Hence, from (40), we getAlso from the theory of matrices we havewhere is element of the matrix for some between 1 and .
Therefore,where . We define and .
From (37), (40), (43), and (45), we getwhich impliesTherefore,that is, our method reduces to a firstorder convergent for uniform mesh.
5. Numerical Examples
To demonstrate the applicability of the method, we have applied the method on four boundary value problems. These examples have been chosen because they have been widely discussed in literature and exact solutions are available for comparison.
The exact solution of the boundary value problem under the boundary conditionsiswhere , , ,
Example 1. Consider the model boundary value problem of the type given by (1)(2) having the boundary layer at the left endwith boundary conditions , and , .
The maximum absolute errors are given in Tables 1 and 2 for different values of the delay and advanced parameters with perturbation parameter. The effect of the small parameters on the boundary layer solutions is shown in Figures 1 and 2.


Example 2. Consider the boundary value problem having the boundary layer at the left endwith boundary conditions , and .
The maximum absolute errors are given in Table 3 for different values of with the delay and advance parameter values . The effect of the small parameters on the boundary layer solutions is shown in Figures 3 and 4.

Example 3. Consider the boundary value problem of the type given by (1)(2) having the boundary layer at the right end with boundary conditions , and .
The maximum absolute errors are given in Table 4 with for different values of the delay and advance parameters. The effect of the small parameters on the boundary layer solutions is shown in Figures 5 and 6.

Example 4. Consider the boundary value problem having the boundary layer at the left end with boundary conditions , and .
The maximum absolute errors are given in Table 5 and Table 6 for different values of delay and advanced parameters with perturbation parameter. The effect of the small parameters on the boundary layer solutions is shown in Figures 7 and 8.


6. Discussion and Conclusion
An exponentially fitted Galerkin method has been presented for solving singularly perturbed differentialdifference equations with delay as well as advance parameters. To demonstrate the applicability of the method, three examples with leftend and one with rightend boundary layer have been solved for different values of the delay, advance, and perturbation parameters. The numerical results are taken by using MATLAB coding and solutions have been compared with the exact solutions and maximum absolute errors are presented in tables. To show the efficiency of the method, we have compared results of the proposed scheme with the results of Kadalbajoo and Sharma [10]. The rate of convergence in the examples is given in Table 7. It is observed that the present method approximates the exact solution very well for which other classical finite difference methods fail to give good results. The effect of the delay and advance parameters on the solutions has also been investigated and presented by using graphs. When the solution of the boundary value problem exhibits layer behaviour on the left side, the effect of delay or advance on the solution in the boundary layer region is negligible while, in the outer region, it is considerable and the change in the advance affects the solution in similar fashion; that is, the width of the layer in the outer region increases as the advance parameter increases and decreases as the delay increases, but reversely (Figures 1–4). When the solution of the boundary value problem exhibits layer behaviour on the right side, the changes in delay or advance affect the solution in boundary layer region as well as outer region. The thickness of the layer increases as the size of the delay increases while it decreases as the size of the advance increases (Figures 5 and 6).

Competing Interests
The authors declare that they have no competing interests.
References
 M. W. Derstine, H. M. Gibbs, F. A. Hopf, and D. L. Kaplan, “Bifurcation gap in a hybrid optically bistable system,” Physical Review A, vol. 26, no. 6, pp. 3720–3722, 1982. View at: Publisher Site  Google Scholar
 V. Y. Glizer, “Asymptotic analysis and solution of a finitehorizon ${H}_{\infty}$ control problem for singularlyperturbed linear systems with small state delay,” Journal of Optimization Theory and Applications, vol. 117, no. 2, pp. 295–325, 2003. View at: Publisher Site  Google Scholar
 C. G. Lange and R. M. Miura, “Singular perturbation analysis of boundaryvalue problems for differentialdifference equations. II. Rapid oscillations and resonances,” SIAM Journal on Applied Mathematics, vol. 45, no. 5, pp. 687–707, 1985. View at: Publisher Site  Google Scholar
 C. G. Lange and R. M. Miura, “Singular perturbation analysis of boundary value problems for differentialdifference equations III. Turning point problems,” SIAM Journal on Applied Mathematics, vol. 45, no. 5, pp. 708–734, 1985. View at: Publisher Site  Google Scholar  MathSciNet
 C. G. Lange and R. M. Miura, “Singular perturbation analysis of boundary value problems for differentialdifference equations. V. Small shifts with layer behavior,” SIAM Journal on Applied Mathematics, vol. 54, no. 1, pp. 249–272, 1994. View at: Publisher Site  Google Scholar  MathSciNet
 C. G. Lange and R. M. Miura, “Singular perturbation analysis of boundaryvalue problems for differentialdifference equations. VI. Small shifts with rapid oscillations,” SIAM Journal on Applied Mathematics, vol. 54, no. 1, pp. 273–283, 1994. View at: Publisher Site  Google Scholar  MathSciNet
 C. G. Lange and R. M. Miura, “Singular perturbation analysis of boundary value problems for differentialdifference equations,” SIAM Journal on Applied Mathematics, vol. 42, no. 3, pp. 502–531, 1982. View at: Publisher Site  Google Scholar  MathSciNet
 M. K. Kadalbajoo and K. K. Sharma, “Numerical analysis of boundaryvalue problems for singularlyperturbed differentialdifference equations with small shifts of mixed type,” Journal of Optimization Theory and Applications, vol. 115, no. 1, pp. 145–163, 2002. View at: Publisher Site  Google Scholar  MathSciNet
 M. K. Kadalbajoo and K. K. Sharma, “Numerical analysis of boundaryvalue problems for singularly perturbed differentialdifference equations: small shifts of mixed type with rapid oscillations,” Communications in Numerical Methods in Engineering, vol. 20, no. 3, pp. 167–182, 2004. View at: Publisher Site  Google Scholar  MathSciNet
 M. K. Kadalbajoo and K. K. Sharma, “Numerical treatment of a mathematical model arising from a model of neuronal variability,” Journal of Mathematical Analysis and Applications, vol. 307, no. 2, pp. 606–627, 2005. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 M. K. Kadalbajoo and D. Kumar, “Fitted mesh Bspline collocation method for singularly perturbed differentialdifference equations with small delay,” Applied Mathematics and Computation, vol. 204, no. 1, pp. 90–98, 2008. View at: Publisher Site  Google Scholar  MathSciNet
 K. C. Patidar and K. K. Sharma, “Uniformly convergent nonstandard finite difference methods for singularly perturbed differentialdifference equations with delay and advance,” International Journal for Numerical Methods in Engineering, vol. 66, no. 2, pp. 272–296, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. K. Kadalbajoo, K. C. Patidar, and K. K. Sharma, “Uniformly convergent fitted methods for the numerical solution of the problems arising from singularly perturbed general DDEs,” Applied Mathematics and Computation, vol. 182, no. 1, pp. 119–139, 2006. View at: Publisher Site  Google Scholar  MathSciNet
 V. Kumar and K. K. Sharma, “An optimized Bspline method for solving singularly perturbed differential difference equations with delay as well as advance,” Neural, Parallel & Scientific Computations, vol. 16, no. 3, pp. 371–385, 2008. View at: Google Scholar  MathSciNet
 L. E. El'sgol'ts and S. B. Norkin, Introduction to the Theory and Application of Differential Equations with Deviating Arguments, Academic Press, 1973.
 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.
 R. E. O'Malley, Introduction to Singular Perturbations, Academic Press, New York, NY, USA, 1974.
 D. J. Evans, “CageGalerkin method for solving the 2point boundary value problem,” International Journal of Computer Mathematics, vol. 43, no. 34, pp. 219–230, 1992. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2016 D. Kumara Swamy et al. 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.