Dear Label_Name,

Greetings! I hope you are well.

I am writing on behalf of an international publishing house, Science and Education Publishing, USA (http://www.sciepub.com).

Science and Education Publishing (SciEP) is a multidisciplinary open access publisher of Journals covering the fields of Science, Technology and Medicine. All SciEP publications are subject to high-quality peer review, editorial and production processes and are freely available online.

Welcome to the April, 2013 edition of SciEP newsletter, providing monthly summaries of articles on science & education related topics of interest to scientist and professionals. The mission of our journals is to provide an open access forum for publishing new findings on issues pertinent to scientific and academic field. We also invite you to submit original research papers, review articles, short reports, or case studies are considered for publication in our Multidisciplinary Journals.

American Journal of Public Health Research
DOI: 10.12691/ajphr-1-1-3

The Implementation of a Tele-Homecare System with Service Oriented Architecture and HL7 Message Transmission Standard
Ching-Sung Wang

American Journal of Medical and Biological Research
DOI: 10.12691/ajmbr-1-1-2

The Influence of Ethanol on Pyruvate Kinases Activity in Vivo, in Vitro, in Silico
Lelevich Sergey Vladimirovich, Khrustalev Vladislav Victorovich, Barkovsky Eugene Victorovich, Shedogubova Tatyana Aleksandrovna

American Journal of Medical Sciences and Medicine
DOI: 10.12691/ajmsm-1-1-1

Clinical and Histological Features in Wegener Granulomatosis
Li Cavoli Gioacchino, Passantino Rita, Bono Luisa, Tortorici Calogera, Ferrantelli Angelo, Giammarresi Carlo, Rotolo Ugo

American Journal of Infectious Diseases and Microbiology
DOI: 10.12691/ajidm-1-1-2

The Curse of Lymphatic Filariasis: Would the Continual Use of Diethylcarbamazine Eliminate this Scourge in Papua New Guinea?
Anne Lanham, Lillian Mwanri

We have released issues of some journals on February 2013.
American Journal of Public Health Research
American Journal of Medical Sciences and Medicine
American Journal of Medical and Biological Research
American Journal of Pharmacological Sciences
American Journal of Clinical Medicine Research
American Journal of Nursing Research
American Journal of Cancer Prevention
American Journal of Medicine Studies
American Journal of Infectious Diseases and Microbiology
American Journal of Cardiovascular Disease Research
American Journal of Biomedical Research
American Journal of Mechanical Engineering
American Journal of Microbiological Research
Journal of Preventive Medicine
American Journal of Food and Nutrition
American Journal of Hypertension Research
American Journal of Sports Science and Medicine
Journal of Cancer Research and Treatment
American Journal of Zoological Research
Journal of Biomedical Engineering and Technology
American Journal of Microbiological Research
Biomedical Science and Engineering

We are pleased to announce that submissions are now being accepted for peer review for all the journals across all our subject areas in SciEP.

Please submit your manuscript via our online submission system:
http://www.sciepub.com/author/signin

To unsubscribe from our mailing list, please click on the following link: unsubscribe

I am looking forward to hearing from you.

Kind regards,

David Rain
Email:join@sciepub.com
Science and Education Publishing (http://www.sciepub.com)

Solution of Singularly Perturbed Differential Difference Equations Using Higher Order Finite Differences

Lakshmi Sirisha, Y. N. Reddy*

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

Abstract: In this paper, we discuss the solution of singularly perturbed differential-difference equations exhibiting dual layer using the higher order finite differences. First, the second order singularly perturbed differential-difference equations is replaced by an asymptotically equivalent second order singular perturbed ordinary differential equation. Then, fourth order stable finite difference scheme is applied to get a three term recurrence relation which is easily solved by Thomas algorithm. Some numerical examples have been solved to validate the computational efficiency of the proposed numerical scheme. To analyze the effect of the parameters on the solution, the numerical solution has also been plotted using graphs. The error bound and convergence of the method have also been established.

Keywords: differential-difference equations, delay parameter, advance parameter, dual layer

Received January 11, 2015; Revised February 10, 2015; Accepted February 13, 2015

Citation: Lakshmi Sirisha, Y. N. Reddy, Solution of Singularly Perturbed Differential Difference Equations Using Higher Order Finite Differences. American Journal of Numerical Analysis, 3 (1). 8-17. January 2015. doi:10.12691/ajna-3-1-2

1. Introduction

A class of functional differential equations which have the characteristics of both classes i.e., delay and/or advanced and singularly perturbed behaviour are known as Singularly Perturbed Differential-Difference Equations (SPDDEs). These differential-difference equation models have richer mathematical framework for the analysis of bio system dynamics, compared with ordinary differential equations. Also, they exhibit better consistency with the nature of the underlying process and predictive results. SPDDEs provide mathematical models in optics, physiology, mechanics, epidemiology, neural networks as well as many other applications Ref. [1,2,3,4]. As a result, numerical treatment of SPDDEs has received a great deal of attention. Stein [5] gave stochastic effects due to neuronal excitation with help of a model that represents a differential-difference equation. Lange and Miura [6,7,8,9,10] published a series of papers to study the numerical solution of SPDDEs. They extended the matched asymptotic expansions, initially developed for ordinary differential equations, to obtain approximate solutions to differential difference equations. The same authors [7] extended the problems involving boundary and interior layer phenomenon, rapid oscillations and resonance behaviour in [6] and [10] to the problems having solutions which exhibit turning point behaviour i.e., transition regions between rapid oscillations and exponential behaviour. Kadalbajoo and Sharma [11] devised a finite difference scheme to obtain the solution for singularly perturbed differential difference equations of mixed type i.e., which contain negative as well as positive shifts. The same authors [12] described a numerical approach based on finite difference method to solve a mathematical model arising from neuronal variability. Sharma and Kaushik [13] discussed two approaches, namely, fitted operator and fitted mesh method to solve SPDDEs with small delay as well as advanced, with layer behaviour. Amiraliyev and Cimen [14] used an exponentially fitted difference scheme on uniform mesh, accomplished by the method of integral identities, for linear second order delay differential equations. Prathima and Sharma [15] presented a numerical method to solve SPDDEs with negative shift and isolated turning point at x=0. Chakravarthy and Rao [16] presented a modified fourth order Numerov method for solving SPDDEs of mixed type.

With this motivation, in this paper, we discuss the solution of singularly perturbed differential-difference equations exhibiting dual layer using the higher order finite differences. First, the second order singularly perturbed differential-difference equations is replaced by an asymptotically equivalent second order singular perturbed ordinary differential equation. Then, fourth order stable finite difference scheme is applied to get a three term recurrence relation which is easily solved by Thomas algorithm. Some numerical examples have been solved to validate the computational efficiency of the proposed numerical scheme. To analyze the effect of the parameters on the solution, the numerical solution has also been plotted using graphs. The error bound and convergence of the method have also been established.

2. Description of the Method

Consider the singularly perturbed differential difference equation of mixed type, i.e., with delay as well as advance parameters of the form:

[picture] | (1)

and subject to the interval and boundary conditions

[picture] | (2)

[picture] | (3)

where are bounded and continuously differentiable functions on (0,1), is the singular perturbation parameter; and and are the delay and the advance parameters respectively.

Using Taylor series expansion in the neighborhood of the point x, we have

[picture] | (4)

[picture] | (5)

Substituting (4) and (5) into Eq. (1), we obtain an asymptotically equivalent singularly perturbed boundary value problem of the form:

[picture] | (6)

[picture] | (7)

[picture] | (8)

where

and .

The transition from Eq.(1) to Eq. (6) is admitted, because of the condition that and are sufficiently small. Further details on the validity of this transition can be found in Elsgolt’s and Norkin [17]. If on the interval , then the solution of Eq. (1) exhibits boundary layers at both ends of the interval , whereas it exhibits oscillatory behaviour for . Here, we have considered the first case where the solutions of the problem exhibit dual layers. Dividing the interval [0, 1] into N equal parts with constant mesh length h, with mesh points:, we have .Assuming that has continuous fourth derivatives on and by making use of ’s expansion, we have the fourth order central differences:

[picture] | (9)

[picture] | (10)

where and for .

Substituting (9) and (10) into Eq. (6) we can write the central difference approximation of Eq.(6) in the form that includes all the error terms as follows:

[picture] | (11)

where .

Now, from Eq.(6) we have

[picture] | (12)

Differentiating both sides of (12), we get:

[picture] | (13)

[picture] | (14)

Using (13) and (14) into (11) we obtain

[picture] | (15)

Now, approximating the converted error term, which has a stabilizing effect, in (15) by using the central difference formulas (9) and (10) for and we obtain

[picture] | (16)

Where

is the truncation error and .

Rearranging Eq. (16) we obtain the three term recurrence relation of the form

[picture] | (17)

for i = 1 (1) N-1.

where

This gives us the tri-diagonal system which can easily be solved by Thomas Algorithm.

3. Thomas Algorithm

We give a brief description of Thomas algorithm for solving the tri-diagonal system. Consider the system:

[picture] | (18)

subject to the boundary conditions

[picture] | (19)

[picture] | (20)

We set

[picture] | (21)

for ,

where and which are to be determined.

From (21), we have:

[picture] | (22)

By substituting (22) in (18) and comparing with (21) we get the recurrence relations:

[picture] | (23)

[picture] | (24)

To solve these recurrence relations for , we need the initial conditions for and . For this we take. We choose so that the value of . With these initial values, we compute and for from (23) and (24) in forward process, and then obtain in the backward process from (20) and (21).

4. Error Analysis

In this section we discuss the error analysis of the method. Writing the tri-diagonal system (18) in matrix-vector form, we get

[picture] | (26)

where, is a tri diagonal matrix of order N-1, with

is a column vector with , where with local truncation error

[picture] | (27)

where

We also have

[picture] | (28)

denotes the actual solution and

is the local truncation error.

From (26) and (28), we get

[picture] | (29)

Thus, we obtain the error equation

[picture] | (30)

where .

Let be the sum of elements of the ith row of A, then we have

is a column vector with , where with local truncation error

Since and , for sufficiently small the matrix A is irreducible and monotone Ref. [18]. Then it follows that exists and its elements are non negative.

Hence, from Eq. (30), we get

[picture] | (31)

and

[picture] | (32)

Let be the element of . Since , from the theory of matrices we have

[picture] | (33)

Therefore,

[picture] | (34)

We define and .

From (23), (27), (28) and (30), we get

which implies

[picture] | (35)

where .

Therefore,

Hence, the method gives a fourth order convergence for uniform mesh as given in Ref. [19].

5. Numerical Examples

To validate the computational efficiency of the scheme, we have applied it to two cases of the problem (1)-(2,3). These cases have been chosen because they have been widely discussed in literature and because approximate solutions are available for comparison.

and subject to the conditions

The exact solution of such boundary value problems having constant coefficients (i.e. and are constants) is given by:

[picture] | (36)

where

Example 1: Consider the model problem given by equations (1)-(3) with

The exact solution of the problem is given by Eq. (36). The numerical results are given in Table 1-Table 4 for =0.1, 0.01 and different values of the delay and advance parameters.

Table 1. Numerical solution of Example 1 for ε=0.1; δ= 0.07 and N=100

Table 2. Numerical solution of Example 1for ε=0.1; δ= 0.05 and N=100

Table 3. Numerical solution of Example 1 for ε=0.01; δ= 0.005 and N=100

Table 4. Numerical solution of Example 1 for ε=0.01; δ= 0.007 and N=100

Example 2: Consider the model problem given by equations (1)-(3) with

The exact solution of the problem is given by Eq. (36). The numerical results are given in Table 5- Table 8 for =0.1, 0.01 and different values of the delay and advance parameters.

Table 5. Numerical solution of Example 2 for ε=0.1; δ= 0.07 and N=100

Table 6. Numerical solution of Example 2 for ε=0.1; η=0.05 and N=100

Table 7. Numerical solution of Example 2 for ε=0.01; η=0.007 and N=100

Table 8. Numerical solution of Example 2 for ε=0.01; δ=0.005 and N=100

Figure 1. Numerical solution of Example 1 for ε=0.1 and δ=0.07

6. Discussion and Conclusion

In this paper, we have discussed the solution of singularly perturbed differential-difference equations exhibiting dual layer using the higher order finite differences. First, the second order singularly perturbed differential-difference equations is replaced by an asymptotically equivalent second order singular perturbed ordinary differential equation. Then, fourth order stable finite difference scheme is applied to get a three term recurrence relation which is easily solved by Thomas algorithm. The error bound and convergence of the method have also been established. Some model examples (i.e., examples in which both the negative and positive shifts are non zero) have been solved to demonstrate the applicability of the proposed approach, by taking various values of the delay, advanced and perturbation parameters. These cases have been chosen because they have been widely discussed in literature Ref. [20,21,22] and also because the exact solutions are available for comparison. The numerical solutions are presented in tables and compared with the exact solution. To analyze the effect of the parameters on the solution, the numerical solution has also been plotted using graphs. From the graphs, it is observed that by varying either or the thickness of the left boundary layer decreases and that of the right boundary layer increases. The effect of negative shift is more dominant than the positive shift.

Figure 2. Numerical solution of Example 1 for ε=0.1 and η=0.05

Figure 3. Numerical solution of Example 1 for ε=0.01 and δ= 0.007

Figure 4. Numerical solution of Example 2 for ε=0.1 and δ=0.07

Figure 5. Numerical solution of Example 2 for ε=0.1 and η=0.05

Figure 6. Numerical solution of Example 2 for ε=0.01 and η=0.007

References

[1]   | Bellman and R. K. L. Cooke, Differential-Difference Equations, Academic Press, , 1963.
[2]   | Derstine, M. W., Gibbs, H. M., Hopf, F. A. and Kalplan, D. L.: Bifurcation gap in hybrid optical system, Phys. Rev. A26(1982)3720-3722.
[3]   | Kolmanovskii, V. and Myshkis, A.:Applied Theory of Functional Differential Equations, Volume 85 of Mathematics and its Applications(Soviet Series), Kluwer Academic Publishers Group, Dordrecht, 1992.
[4]   | H. C. Tuckwell and W. Richter, Neuronal inter-spike time distributions and the estimation of neuro-physiological and neuro-anatomical parameters, J. Theor. Biol., 71 (1978) 167-183.
[5]   | Stein R.B., A theoretical analysis of neuronal variability, Biophys. J., 5 (1965) 173-194.
[6]   | C.G. Lange R.M. Miura, Singular perturbation analysis of boundary-value problems for differential-difference equations. II. Rapid oscillations and resonances, SIAM J. Appl. Math., 45 (1985) 687-707.
[7]   | C.G. Lange R.M. Miura, Singular perturbation analysis of boundary- value problems for differential-difference equations. III. Turning point problems, SIAM J. Appl. Math., 45 (1985) 708-734.
[8]   | C.G. Lange R.M. Miura, Singular perturbation analysis of boundary-value problems for differential-difference equations. V. Small shifts with layer behaviour, SIAM J. Appl. Math., 54 (1994) 249-272.
[9]   | C.G. Lange R.M. Miura, Singular perturbation analysis of boundary-value problems for differential-difference equations. VI. Small shifts with rapid oscillations, SIAM J. Appl. Math., 54 (1994) 273-283.
[10]   | C.G. Lange R.M. Miura, Singular perturbation analysis of boundary-value problems for differential difference equations, SIAM J. Appl. Math., 42 (1982) 502-531.
[11]   | M. K. Kadalbajoo and K. K. Sharma, Numerical analysis of boundary-value problems for singularly perturbed differential-difference equations with small shifts of mixed type, J. Optim. Theory Appl., 115 (1) (2002) 145-163.
[12]   | M. K. Kadalbajoo and K. K. Sharma, Numerical treatment of a mathematical model arising from a model of neuronal variability, J. Math. Anal. Appl., 307 (2005) 606-627.
[13]   | Sharma, K. K. and Kaushik, A.: A solution of the discrepancy occurs due to using the fitted mesh approach rather than to the fitted operator for solving singularly perturbed differential equations, Appl. Math. Comput., 181(2006) 756-766.
[14]   | Amiraliyeva, I. G., and Cimen, E: Numerical method for a singularly perturbed convection-diffusion problem with delay, Applied Mathematics and Computation, 216(2010) 2351-2359.
[15]   | Pratima, R. and Sharma, K. K.: Numerical method for singularly perturbed differential-difference equations with turning point, International Journal of Pure and Applied Mathematics, 73(4) (2011) 451-470.
[16]   | Chakarvarthy, P. P. and Rao R. N.: A modified Numerov method for solving singularly perturbed differential-difference equations arising in science and engineering, Results in Physics, 2 (2012) 100-103.
[17]   | L. E. El’sgolts and S. B. Norkin, Introduction to the Theory and Applications of Differential Equations with Deviating Arguments, Academic Press, New York, 1973.
[18]   | R. K. Mohanty and N. Jha, A class of variable mesh spline in compression methods for singularly perturbed two point singular boundary value problems, Applied Mathematics and Computation, 168 (2005) 704-716.
[19]   | J. Y. Choo and D. H.Schultz: Stable higher order methods for differential equations with small coefficients for the second order terms, Journal of Computers and Math. with Applications, 25 (1993), 105-123.
[20]   | Pratima, R. and Sharma, K. K.: Numerical study of singularly perturbed differential difference equation arising in the modelling of neuronal variability, Computers and Mathematics with Applications, 63 (2012) 118-132.
[21]   | Rao, R. N. and Chakravarthy, P. P.: An initial value technique for singularly perturbed differential-difference equations with a small negative shift, J. Appl. Math. & Informatics, 31 (2013) 131-145.
[22]   | Gemechis FD and Reddy Y.N.: A Non-Asymptotic Method for Singularly Perturbed Delay Differential Equations, Journal of Applied Mathematics and Informatics, 32:1(2014) 39-53.