International Journal of Differential Equations

Volume 2012 (2012), Article ID 572723, 12 pages

doi:10.1155/2012/572723

Research Article

Numerical Integration of a Class of Singularly Perturbed Delay Differential Equations with Small Shift

Gemechis File and Y. N. Reddy

Department of Mathematics, National Institute of Technology, Warangal 506 004, India

Received 22 May 2012; Accepted 1 October 2012

Academic Editor: Samir H. Saker

Copyright © 2012 Gemechis File and Y. N. Reddy. 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.

Abstract.

We have presented a numerical integration method to solve a class of singularly perturbed delay differential equations with small shift. First, we have replaced the second-order singularly perturbed delay differential equation by an asymptotically equivalent first-order delay differential equation. Then, Simpson’s rule and linear interpolation are employed to get the three-term recurrence relation which is solved easily by discrete invariant imbedding algorithm. The method is demonstrated by implementing it on several linear and nonlinear model examples by taking various values for the delay parameter and the perturbation parameter .

1. Introduction

The singularly perturbed delay differential equations with small shift arise very frequently in the modeling of various physical and biological phenomena, for example, micro scale heat transfer [1], hydrodynamics of liquid helium [2], second-sound theory [3], thermoelasticity [4], diffusion in polymers [5], reaction-diffusion equations [6], stability [7], control of chaotic systems [8], a variety of models for physiological processes or diseases [9] and so forth. Hence in the recent times, many researchers have been trying to develop numerical methods for solving these problems. Amiraliyev and Cimen [10] presented numerical method comprising a fitted difference scheme on a uniform mesh to solve second-order delay differential equations. Lange and Miura [11, 12] gave an asymptotic approach for a class of boundary-value problems for linear second-order differential-difference equations. Kadalbajoo and Sharma [13–15] presented numerical approaches to solve singularly perturbed differential-difference equations, which contains negative shift in the convention term (i.e., in the derivative term). Lange and Miura [16] considered the boundary value problem for a singularly perturbed nonlinear differential difference equation with shift and discussed the existence and uniqueness of their solutions. Furthermore, Kadalbajoo and Sharma [17] have discussed the numerical solution of the singularly perturbed nonlinear differential equations with small negative shifts.

In this paper, we have presented a numerical integration method for solving a class of singularly perturbed delay differential equations with small shift. First, the second-order singularly perturbed delay differential equation is replaced by an asymptotically equivalent first-order delay differential equation. Then we employed Simpson’s rule and linear interpolation to get three-term recurrence relation which is solved easily by discrete invariant imbedding algorithm. The method is demonstrated by implementing it on several linear and nonlinear model examples by taking various values for the delay and perturbation parameters.

2. Description of the Method

Consider a class of singularly perturbed boundary value problems of the following form:

with the interval and boundary conditions

where is small parameter, , and is also a small shifting parameter, ; , and are bounded continuous functions in , and are finite constants. Further, we assume that throughout the interval , where is positive constant. This assumption merely implies that the boundary layer will be in the neighborhood of .

By using Taylor series expansion in the neighborhood of the point , we have

and consequently, (2.1) is replaced by the following first-order differential equation:

where

The transition from (2.1) to (2.4) is admitted, because of the condition that is small, . This replacement is significant from the computational point of view. Further details on the validity of this transition can be found in [18].

Now we divide the interval into equal subintervals of mesh size so that  ,  .

Integrating (2.4) with respect to from to for , we get

where ,  ,  ,  ,  .

By using Simpson’s rule to evaluate the integral in (2.6), we get

By the means of Taylor series expansion and then by approximating by linear interpolation, we get

In similar way,

Hence, by making use of (2.8a)–(2.8e) in (2.7) we obtain

To make (2.9) a three-term recurrence relation, we can express in terms of , and using Hermite’s interpolation as follows:

In view of (2.4) and (2.10), we get

By making use of (2.8a)–(2.8e) in (2.11) and finite difference approximations, we get

Finally, making use of (2.12) in (2.9) and rearranging as three-term recurrence relation, we get

for , where

This tridiagonal system is solved by using method of discrete invariant imbedding algorithm which is described in the next section.

3. Discrete Invariant Imbedding Algorithm

We now describe the Thomas algorithm which is also called discrete invariant imbedding [19] to solve the three-term recurrence relation:

Let us set a difference relation of the form

where and are to be determined.

From (3.2), we have

Substituting (3.3) in (3.1), we have

By comparing (3.2) and (3.4), we get the recurrence relations

To solve these recurrence relations for , we need the initial conditions for and . If we choose , then we get . With these initial values, we compute and for from (3.5) and (3.6) in forward process and then obtain in the backward process from (3.2).

The conditions for the discrete invariant imbedding algorithm to be stable are (see [18–21])

In our method, one can easily show that if the assumptions , and hold, then the above conditions (3.7) hold, and thus the discrete invariant imbedding algorithm is stable.

4. Numerical Experiments

To demonstrate the applicability of the method, we have implemented it on two linear and two nonlinear problems with left-end boundary layers. Computational results are compared with exact solutions wherever exact solutions are available. When exact solution is not available, we have tested the effect of small delay parameter on solution of the problem for different values of of .

4.1. Linear Problems

Example 4.1. Consider an example of singularly perturbed delay differential equation with left layer:

The exact solution is given by

where and .

The computational results are presented in Tables 1, 2, 3, and 4 for and 0.0001 for different values of .

Table 1: Numerical results of Example 4.1 for , .

Numerical solution
Exact solution
Absolute error

0.00
1.0000000
1.0000000
0.000E + 00

0.01
0.3724909
0.3719167
5.743E − 04

0.02
0.3753635
0.3756417
2.781E − 04

0.03
0.3791343
0.3794135
2.792E − 04

0.04
0.3829441
0.3832233
2.791E − 04

0.06
0.3906790
0.3909578
2.788E − 04

0.08
0.3985702
0.3988485
2.784E − 04

0.20
0.4493791
0.4496520
2.730E − 04

0.50
0.6065730
0.6068032
2.302E − 04

0.60
0.6703575
0.6705610
2.035E − 04

0.90
0.9048500
0.9049187
6.870E − 05

1.00
1.0000000
1.0000000
0.000E + 00

Table 2: Numerical results of Example 4.1 for , .

Numerical solution
Exact solution
Absolute error

0.00
1.0000000
1.0000000
0.000E + 00

0.02
0.3785059
0.3753847
3.121E − 03

0.03
0.3786281
0.3791566
5.284E − 04

0.04
0.3827057
0.3829664
2.607E − 04

0.05
0.3865343
0.3868145
2.802E − 04

0.06
0.3904227
0.3907013
2.786E − 04

0.08
0.3983141
0.3985924
2.782E − 04

0.20
0.4491281
0.4494008
2.728E − 04

0.40
0.5486277
0.5488775
2.498E − 04

0.60
0.6701703
0.6703736
2.033E − 04

0.90
0.9047868
0.9048555
6.870E − 05

1.00
1.0000000
1.0000000
0.000E + 00

Table 3: Numerical results of Example 4.1 for , .

Numerical solution
Exact solution
Absolute error

0.00
1.0000000
1.0000000
0.000E + 00

0.01
0.3737204
0.3716098
2.111E − 03

0.02
0.3754769
0.3753442
1.327E − 04

0.03
0.3792426
0.3791161
1.264E − 04

0.04
0.3830523
0.3829260
1.264E − 04

0.08
0.3986781
0.3985520
1.261E − 04

0.20
0.4494850
0.4493613
1.237E − 04

0.40
0.5489547
0.5488413
1.134E − 04

0.50
0.6066623
0.6065580
1.043E − 04

0.80
0.8188018
0.8187455
5.630E − 05

0.90
0.9048767
0.9048455
3.120E − 05

1.00
1.0000000
1.0000000
0.000E + 00

Table 4: Numerical results of Example 4.1 for , .

Numerical solution
Exact solution
Absolute error

0.00
1.0000000
1.0000000
0.000E + 00

0.02
0.3754557
0.3753185
1.372E − 04

0.03
0.3792183
0.3790904
1.279E − 04

0.04
0.3830281
0.3829003
1.279E − 04

0.06
0.3907630
0.3906352
1.278E − 04

0.08
0.3986540
0.3985264
1.276E − 04

0.20
0.4494613
0.4493362
1.251E − 04

0.40
0.5489329
0.5488182
1.146E − 04

0.60
0.6704187
0.6703254
9.330E − 05

0.70
0.7409000
0.7408227
7.730E − 05

0.90
0.9048707
0.9048392
3.150E − 05

1.00
1.0000000
1.0000000
0.000E + 00

Example 4.2. Now we consider an example of variable coefficient singularly perturbed delay differential equation with left layer:

For which the exact solution is not known. This example is considered to show the effect of the small shift on the boundary layer solution.

The computational results are presented in Tables 5 and 6 for and 0.0001 for different values of .

Table 5: Numerical results of Example 4.2 for , , and different values of .

Numerical solutions

= 0.0001
= 0.0003
= 0.0006
= 0.0008

0.00
1.0000000
1.0000000
1.0000000
1.0000000

0.02
0.2608070
0.2574476
0.2531103
0.2507717

0.04
0.2666194
0.2664989
0.2663023
0.2661506

0.05
0.2696995
0.2695796
0.2694007
0.2692834

0.06
0.2728263
0.2727058
0.2725251
0.2724043

0.08
0.2792234
0.2791020
0.2789198
0.2787981

0.20
0.3220214
0.3218944
0.3217039
0.3215766

0.40
0.4142966
0.4141648
0.4139674
0.4138352

0.60
0.5434980
0.5433738
0.5431879
0.5430634

0.80
0.7285067
0.7284169
0.7282822
0.7281920

0.90
0.8508641
0.8508093
0.8507276
0.8506728

1.00
1.0000000
1.0000000
1.0000000
1.0000000

Table 6: Numerical results of Example 4.2 for , , and different values of .

Numerical solutions

= 0.00001
= 0.00003
= 0.00006
= 0.00008

0.00
1.0000000
1.0000000
1.0000000
1.0000000

0.01
0.4289609
0.4278078
0.4260698
0.4249051

0.02
0.2612188
0.2608663
0.2603437
0.2599993

0.03
0.2637076
0.2636943
0.2636772
0.2636673

0.04
0.2667407
0.2667287
0.2667109
0.2666990

0.06
0.2729484
0.2729363
0.2729184
0.2729064

0.09
0.2826190
0.2826068
0.2825886
0.2825764

0.20
0.3221486
0.3221358
0.3221169
0.3221042

0.60
0.5436167
0.5436044
0.5435862
0.5435733

0.70
0.6275683
0.6275574
0.6275409
0.6275294

0.90
0.8509144
0.8509088
0.8509008
0.8508952

1.00
1.0000000
1.0000000
1.0000000
1.0000000

4.2. Nonlinear Problems

Nonlinear problems are linearized by the quasilinearization process. Then we have applied the present method.

Example 4.3. Consider a singularly perturbed nonlinear delay differential equation:

under the interval and boundary conditions

The exact solution is not known.

The computational results are presented in Tables 7 and 8 for for different values of .

Table 7: Numerical results of Example 4.3 for , , and different values of .

Numerical solutions

= 0.0001
= 0.0003
= 0.0006
= 0.0008

0.00
1.0000000
1.0000000
1.0000000
1.0000000

0.01
0.3724909
0.3596174
0.3392884
0.3250014

0.03
0.3791343
0.3790570
0.3788695
0.3786281

0.04
0.3829441
0.3828712
0.3827668
0.3827057

0.05
0.3867922
0.3867193
0.3866106
0.3865343

0.06
0.3906790
0.3906061
0.3904977
0.3904227

0.08
0.3985702
0.3984973
0.3983891
0.3983141

0.20
0.4493791
0.4493076
0.4492016
0.4491281

0.40
0.5488575
0.5487921
0.5486948
0.5486277

0.60
0.6703575
0.6703041
0.6702249
0.6701703

0.90
0.9048500
0.9048320
0.9048053
0.9047868

1.00
1.0000000
1.0000000
1.0000000
1.0000000

Table 8: Numerical results of Example 4.3 for , , and different values of .

Numerical solutions

= 0.00001
= 0.00003
= 0.00006
= 0.00008

0.00
1.0000000
1.0000000
1.0000000
1.0000000

0.01
0.3737204
0.3724622
0.3705641
0.3692932

0.03
0.3792426
0.3792360
0.3792249
0.3792183

0.04
0.3830523
0.3830458
0.3830347
0.3830281

0.05
0.3869004
0.3868939
0.3868828
0.3868762

0.08
0.3986781
0.3986716
0.3986605
0.3986540

0.20
0.4494850
0.4494785
0.4494676
0.4494613

0.40
0.5489547
0.5489486
0.5489386
0.5489329

0.50
0.6066623
0.6066567
0.6066476
0.6066423

0.70
0.7409147
0.7409106
0.7409040
0.7409000

0.90
0.9048767
0.9048751
0.9048723
0.9048707

1.00
1.0000000
1.0000000
1.0000000
1.0000000

Example 4.4. Consider an example of singularly perturbed nonlinear delay differential equation:

under the interval and boundary conditions

The exact solution is not known.

The computational results are presented in Tables 9 and 10 for and 0.001 for different values of  .

Table 9: Numerical results of Example 4.4 for , , and different values of .

Numerical solutions

= 0.0001
= 0.0003
= 0.0006
= 0.0008

0.00
0.0000000
0.0000000
0.0000000
0.0000000

0.02
−0.2107353
−0.2105360
−0.2098469
−0.2090868

0.04
−0.2053359
−0.2053088
−0.2052677
−0.2052346

0.05
−0.2026495
−0.2026230
−0.2025854
−0.2025608

0.06
−0.1999765
−0.1999504
−0.1999132
−0.1998883

0.08
−0.1946704
−0.1946451
−0.1946091
−0.1945850

0.10
−0.1894171
−0.1893927
−0.1893577
−0.1893344

0.30
−0.1396748
−0.1396577
−0.1396331
−0.1396166

0.60
−0.0737937
−0.0737854
−0.0737733
−0.0737652

0.80
−0.0350537
−0.0350500
−0.0350445
−0.0350408

0.90
−0.0170888
−0.0170870
−0.0170844
−0.0170827

1.00
0.0000000
0.0000000
0.0000000
0.0000000

Table 10: Numerical results of Example 4.4 for , , and different values of .

Numerical Solutions

= 0.00001
= 0.00003
= 0.00006
= 0.00008

0.00
0.0000000
0.0000000
0.0000000
0.0000000

0.02
−0.2107172
−0.2107119
−0.2106993
−0.2106895

0.03
−0.2080040
−0.2080019
−0.2079975
−0.2079955

0.04
−0.2053046
−0.2053025
−0.2052982
−0.2052962

0.05
−0.2026187
−0.2026166
−0.2026124
−0.2026104

0.08
−0.1946411
−0.1946391
−0.1946350
−0.1946331

0.10
−0.1893887
−0.1893868
−0.1893829
−0.1893810

0.30
−0.1396550
−0.1396536
−0.1396509
−0.1396495

0.60
−0.0737842
−0.0737834
−0.0737821
−0.0737814

0.80
−0.0350494
−0.0350491
−0.0350485
−0.0350481

0.90
−0.0170868
−0.0170866
−0.0170863
−0.0170862

1.00
0.0000000
0.0000000
0.0000000
0.0000000

5. Discussions and Conclusions

We have presented a numerical integration method to solve singularly perturbed delay differential equations. The scheme is repeated for different choices of the delay parameter, , and perturbation parameter, . The choice of is not unique but can assume any number of values satisfying the condition with and is not too large Lange and Miura [12]. To demonstrate the efficiency of the method, we have implemented it on two linear and two nonlinear model examples with the boundary layer on the left for different values of and . From the computational results, it is observed that the proposed method approximates the exact solution very well (see Tables 1–4), and the small shift, , affects the boundary layer solutions. That is, as increases, the size/thickness of the left boundary layer decreases (see Tables 5–10). This method does not depend on asymptotic expansion as well as on the matching of the coefficients. Thus, we have devised an alternative technique of solving boundary value problems for singularly perturbed delay differential equations, which is easily implemented on computer and is also practical.

References

D. Y. Tzou, Micro-to-Macroscale Heat Transfer, Taylor and Francis, Washington, DC, USA, 1997.

D. D. Joseph and L. Preziosi, “Heat waves,” Reviews of Modern Physics, vol. 61, no. 1, pp. 41–73, 1989.

D. D. Joseph and L. Preziosi, “Addendum to the paper heat waves,” Reviews of Modern Physics, vol. 62, no. 2, pp. 375–391, 1990.

M. A. Ezzat, M. I. Othman, and A. M. S. El-Karamany, “State space approach to two-dimensional generalized thermo-viscoelasticity with two relaxation times,” International Journal of Engineering Science, vol. 40, no. 11, pp. 1251–1274, 2002.

Q. Liu, X. Wang, and D. De Kee, “Mass transport through swelling membranes,” International Journal of Engineering Science, vol. 43, pp. 1464–1470, 2005.

M. Bestehorn and E. V. Grigorieva, “Formation and propagation of localized states in extended systems,” Annalen der Physik, vol. 13, no. 7-8, pp. 423–431, 2004.

T. A. Burton, “Fixed points, stability, and exact linearization,” Nonlinear Analysis, vol. 61, no. 5, pp. 857–870, 2005.

X. Liao, “Hopf and resonant codimension two bifurcation in van der Pol equation with two time delays,” Chaos, Solitons & Fractals, vol. 23, no. 3, pp. 857–871, 2005.

M. C. Mackey and L. Glass, “Oscillations and chaos in physiological control systems,” Science, vol. 197, pp. 287–289, 1977.

G. M. Amiraliyev and E. Cimen, “Numerical method for a singularly perturbed convection-diffusion problem with delay,” Applied Mathematics and Computation, vol. 216, no. 8, pp. 2351–2359, 2010.

C. G. Lange and R. M. Miura, “Singular perturbation analysis of boundary value problems for differential-difference equations. V. Small shifts with layer behavior,” SIAM Journal on Applied Mathematics, vol. 54, no. 1, pp. 249–272, 1994.

C. G. Lange and R. M. Miura, “Singular perturbation analysis of boundary value problems for differential-difference equations. VI. Small shifts with rapid oscillations,” SIAM Journal on Applied Mathematics, vol. 54, no. 1, pp. 273–283, 1994.

M. K. Kadalbajoo and K. K. Sharma, “Numerical treatment of boundary value problems for second order singularly perturbed delay differential equations,” Computational & Applied Mathematics, vol. 24, no. 2, pp. 151–172, 2005.

M. K. Kadalbajoo and K. K. Sharma, “A numerical method based on finite difference for boundary value problems for singularly perturbed delay differential equations,” Applied Mathematics and Computation, vol. 197, no. 2, pp. 692–707, 2008.

M. K. Kadalbajoo and K. K. Sharma, “Numerical analysis of singularly perturbed delay differential equations with layer behavior,” Applied Mathematics and Computation, vol. 157, no. 1, pp. 11–28, 2004.

C. G. Lange and R. M. Miura, “Singular perturbation analysis of boundary value problems for differential-difference equations. IV. A nonlinear example with layer behavior,” Studies in Applied Mathematics, vol. 84, no. 3, pp. 231–273, 1991.

M. K. Kadalbajoo and K. K. Sharma, “Numerical treatment for singularly perturbed nonlinear differential difference equations with negative shift,” Nonlinear Analysis, vol. 63, pp. 1909–1924, 2005.

L. E. Elsgolt’s and S. B. Norkin, Introduction to the Theory and Applications of Differential Equations with Deviating Arguments, Academic Press, New York, NY, USA, 1973.

E. Angel and R. Bellman, Dynamic Programming and Partial Differential Equations, Academic Press, New York, NY, USA, 1972.

V. Y. Glizer, “Asymptotic analysis and solution of a finite-horizon control problem for singularly-perturbed linear systems with small state delay,” Journal of Optimization Theory and Applications, vol. 117, no. 2, pp. 295–325, 2003.

M. K. Kadalbajoo and Y. N. Reddy, “A non asymptotic method for general linear singular perturbation problems,” Journal of Optimization Theory and Applications, vol. 55, pp. 256–269, 1986.