`Discrete Dynamics in Nature and SocietyVolume 2011, Article ID 628369, 12 pageshttp://dx.doi.org/10.1155/2011/628369`
Research Article

## Dynamics of a Discrete Internet Congestion Control Model

School of Mathematics and Computer Science, Fujian Normal University, Fuzhou 350007, China

Received 24 July 2011; Accepted 1 September 2011

Academic Editor: Carlo Piccardi

Copyright © 2011 Yingguo Li. 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 consider a discrete Internet model with a single link accessed by a single source, which responds to congestion signals from the network. Firstly, the stability of the equilibria of the system is investigated by analyzing the characteristic equation. By choosing the time delay as a bifurcation parameter, we prove that Neimark-Sacker bifurcations occur when the delay passes a sequence of critical values. Then, the explicit algorithm for determining the direction of the Neimark-Sacker bifurcations and the stability of the bifurcating periodic solutions is derived. Finally, some numerical simulations are given to verify the theoretical analysis.

#### 1. Introduction

Congestion control algorithms and active queue management (AQM) for Internet have been the focus of intense research since the seminal work of Kelly et al. [1]. The primal algorithms they have introduced are based on explicit feedback functions, where the sum of user utilities is maximized within the capacity constraints of the links. In primal algorithms, the users adapt the source rates dynamically based on the congestion signal generated by the link, and the links select a static law to determine the link prices directly from the arrival rates at the links [2]. In [3], a stability condition was provided for a single proportionally fair congestion controller with delayed feedback. Since then, this result was extended to networks in [4, 5] and for different classes of controllers in [6].

In this paper, we consider an Internet model with a single link and single source, which can be formulated as a congestion control system with feedback delay. The model is described by [3, 7, 8] where is the sending rate of the source at time , is a positive gain parameter, is the sum of forward and return delays, is a target (set point), and the congestion indication function is increasing, nonnegative, and not identically zero, which can be viewed as the probability that a packet at the resource receives a feedback congestion indication signal. System (1.1) corresponds to a rate control algorithm for the user, which can be used to adjust the sending rate so that the expected number of signals received by the user will tend to target . A great deal of research has been devoted to the asymptotic behaviour, periodic solutions, global stability, and other properties of this model, for which, we refer to [3, 79].

But considering the need of scientific computation and real-time simulation, our interest is focused on the behaviors of discrete dynamics system corresponding to (1.1). Many authors considered the numerical approximation of a scalar delay differential equation by using different numerical methods, such as -method, Euler method, and Trapezoidal method (see [1015]). In this paper, we use the forward Euler scheme to make the discretization for system (1.1).

Moreover, it is also of interest to find what will happen when the congestion control system loses stability. The purpose of this paper is to discuss this version as a discrete dynamical system by using Neimark-Sacker bifurcation theory of discrete systems. We not only investigate the stability of the fixed point and the existence of the Neimark-Sacker bifurcations, but also the stability and direction of the Neimark-Sacker bifurcation of the discrete system.

The paper is organized as follows: in Section 2, we analyze the distribution of the characteristic equation associated with the discrete model and obtain the existence of the local Neimark-Sacker bifurcation. In Section 3, the direction and stability of closed invariant curve from the Neimark-Sacker bifurcation of the discrete delay model are determined by using the theories of discrete systems in [16]. In the final section, some computer simulations are performed to illustrate the analytical results found.

#### 2. Stability Analysis

Let . Then (1.1) can be rewritten as We consider sufficient small step size of the form , where . The Euler method applied to this equation yields the delay difference equation where is an approximate value to . Let be a positive equilibrium of (2.2), then we have Set , then there follows: By introducing a new variable , we can rewrite (2.4) in the form where , and Clearly the origin is a fixed point of (2.5), and the linear part of (2.5) is where The characteristic equation of is given by It is well known that the stability of the zero equilibrium solution of (2.5) depends on the distribution of the zeros of the roots of (2.9). In order to proof the existence of the local Neimark-Sacker bifurcation at fixed point, we need some lemmas as follows.

Lemma 2.1. There exists a such that for all roots of (2.9) have modulus less than one.

Proof. When , (2.9) becomes The equation has, at , an -fold root and a simple root .
Consider the root such that . This root depends continuously on and is a differential function of . From (2.9), we have Since is increasing and nonnegative, we have So with the increasing of , cannot cross . Consequently, all roots of (2.9) lie in the unit circle for sufficiently small positive , and the existence of the follows.

A Neimark-Sacker bifurcation occurs when a complex conjugate pair of eigenvalues of cross the unit circle as varies. We have to find values of such that there are roots on the unit circle. Denote the roots on the unit circle by , . Since we are dealing with complex roots of a real polynomial, we only need to look for . For , is a root of (2.9) if and only if Separating the real part and imaginary part from (2.13), there are So From (2.14) and (2.15), we obtain that that is, then there exists a sequence of where denotes the greatest integer function. It is clear that there exists a sequence of the time delay parameters satisfying (2.14) and (2.15), according to .

From (2.15), we have Differentiating (2.16) with respect to , we obtain that This gives Therefore,

Lemma 2.2. Let be a root of (2.9) near satisfying and . Then

Proof. From (2.11), we have From (2.14) and (2.15), we obtain that It is easy to see that
Using (2.25)–(2.27), we have This completes the proof.

Lemma 2.3. Equation (2.9) has a pair of simple roots on the unit circle when . Furthermore, if , then all the roots of (2.9) have modulus less than one; If , then all roots of (2.9) except have modulus less than one. But if , for , (2.9) has , roots have modulus more than one.

Proof. From (2.18)–(2.20), we have that (2.9) has toots if and only if and given in (2.19) and (2.20). By Lemmas 2.1 and 2.2, we know that if , then all roots of (2.9) have modulus less than one; if , then all roots of (2.9) except have modulus less than one; furthermore, by the Rouchs theorem (Dieudonné [17, Theorem 9.17.4]), the statement on the number of eigenvalues with modulus more than one is as follows.

Lemma 2.3 immediately leads to stability of the zero solution of (2.4), and equivalently, of the positive fixed point of (2.2). So we have the following results on stability and bifurcation in system (2.2).

Theorem 2.4. There exists a sequence of values of the time delay parameter such that the positive fixed point of (2.2) is asymptotically stable for and unstable for . Equation (2.2) undergoes a Neimark-Sacker bifurcation at the when , where satisfies (2.20).

#### 3. Direction and Stability of the Neimark-Sacker Bifurcation in Discrete Model

In the previous section, we obtain the conditions under which a family of periodic solutions bifurcate from the steady state at the critical value . Without loss of generality, denote the critical value by . For continuous systems, as pointed out in Hassard et al. [18], it is interesting and important to determine the direction, stability, and the period of these periodic solutions bifurcating from the steady state. In this section, following the idea of Hassard et al., we shall study the direction, stability, and the period of the bifurcating periodic solution when in the discrete Internet congestion control model. The method we used is based on the theories of discrete system by Kuznetsov [16].

Set . Then is a Neimark-Sacker bifurcation value for (2.4). Rewrite (2.4) as So system (2.5) is turned into where Let be an eigenvector of corresponding to , then We also introduce an adjoint eigenvector having the properties and satisfying the normalization , where .

Lemma 3.1 (see [19]). Define a vector-valued function by If is an eigenvalue of , then .

In view of Lemma 3.1, we have

Lemma 3.2. Suppose that is the eigenvector of corresponding to eigenvalue , and . Then where and

Proof. Assign satisfies with , then the following identities hold Let , then From normalization and computation, (3.10) is held.

Let be characteristic polynomial of and . Following the algorithms in [16] and using a computation process similar to that in [19], we can compute an expression for the critical coefficient , where By (3.4), (3.8), and Lemma 3.2, we get Substituting these into (3.13), we have

Lemma 3.3 (see [20]). Given the map (2.5) and assume that(1), where and ;(2) for ;(3).Then an invariant closed curve, topologically equivalent to a circle, for map (2.5) exists for in a one-side neighborhood of . The radius of the invariant curve grows like . One of the four cases below applies:(1). The origin is asymptotically stable for and unstable for . An attracting invariant closed curve exists for .(2). The origin is asymptotically stable for and unstable for . An repelling invariant closed curve exists for .(3). The origin is asymptotically stable for and unstable for . An attracting invariant closed curve exists for .(4). The origin is asymptotically stable for and unstable for . An repelling invariant closed curve exists for .

From the discussion in Section 2, we know that ; therefore, by Lemma 3.3, we have the following result.

Theorem 3.4. For (2.2), is asymptotically stable for and unstable for . An attracting (repelling) invariant closed curve exists for if .

#### 4. Computer Simulation

In this section, we will confirm our theoretical analysis by numerical simulation. We choose as in [8] and set . Then is the Neimark-Sacker bifurcation value.

Figures 1 and 2 are about delay difference equation (2.2) when step size .

Figure 1: (a) The equilibrium of (2.2) is asymptotically stable for , (b) Phase plot of (2.2) on the plane for .
Figure 2: (a) Waveform plot of (2.2) for , (b) A bifurcating periodic solution appears on the plane for .

In Figure 1, we show the waveform plot and phase plot for (2.2) with initial values for . The equilibrium of (2.2) is asymptotically stable. In Figure 2, we show the waveform plot and phase plot for (2.2) with initial values . The equilibrium of (2.2) is unstable for . When varies and passes through , the equilibrium loses its stability and a periodic solution bifurcates from the equilibrium for . That is the delay difference equation (2.2) which has a Neimark-Sacker bifurcation at .

#### Acknowledgment

The authors are very grateful to the referees for their valuable suggestions, which help to improve the paper significantly.

#### References

1. F. P. Kelly, A. K. Maulloo, and D. Tan, “Rate control for communication networks: Shadow prices, proportional fairness and stability,” Journal of the Operational Research Society, vol. 49, no. 3, pp. 237–252, 1998.
2. T. Alpcan and T. Basar, “A globally stable adaptive congestion control scheme for internet-style networks with delay,” IEEE/ACM Transactions on Networking, vol. 13, no. 6, pp. 1261–1274, 2005.
3. F. P. Kelly, “Models for a self-managed Internet,” Philosophical Transactions of the Royal Society A, vol. 358, no. 1773, pp. 2335–2348, 2000.
4. L. Massoulié, “Stability of distributed congestion control with heterogeneous feedback delays,” IEEE. Transactions on Automatic Control, vol. 47, no. 6, pp. 895–902, 2002.
5. X. F. Wang, G. Chen, and K.-T. Ko, “A stability theorem for internet congestion control,” Systems & Control Letters, vol. 45, no. 2, pp. 81–85, 2002.
6. F. Paganini, J. Doyle, and S. Low, “Scalable laws for stable network congestion control,” in Proceedings of the IEEE Conference on Decision and Control, vol. 1, pp. 185–190, 2001.
7. S. Deb and R. Srikant, “Global stability of congestion controllers for the internet,” in Proceedings of the IEEE Conference on Decision and Control, vol. 4, pp. 3626–3631, 2002.
8. C. Li, G. Chen, X. Liao, and J. Yu, “Hopf bifurcation in an internet congestion control model,” Chaos, Solitons and Fractals, vol. 19, no. 4, pp. 853–862, 2004.
9. G. Raina, “Local bifurcation analysis of some dual congestion control algorithms,” Institute of Electrical and Electronics Engineers. Transactions on Automatic Control, vol. 50, no. 8, pp. 1135–1146, 2005.
10. C. Zhang, Y. Zu, and B. Zheng, “Stability and bifurcation of a discrete red blood cell survival model,” Chaos, Solitons and Fractals, vol. 28, no. 2, pp. 386–394, 2006.
11. Z. He, X. Lai, and A. Hou, “Stability and Neimark-Sacker bifurcation of numerical discretization of delay differential equations,” Chaos, Solitons and Fractals, vol. 41, no. 4, pp. 2010–2017, 2009.
12. H. Su and X. Ding, “Dynamics of a nonstandard finite-difference scheme for Mackey-Glass system,” Journal of Mathematical Analysis and Applications, vol. 344, no. 2, pp. 932–941, 2008.
13. X. Ding and H. Su, “Dynamics of a discretization physiological control system,” Discrete Dynamics in Nature and Society, Article ID 51406, 16 pages, 2007.
14. J. F. Neville and V. Wulf, “Numerical Hopf bifurcation for a class of delay differential equations,” Journal of Computational and Applied Mathematics, vol. 115, no. 1-2, pp. 601–616, 2000.
15. C. Zhang, M. Liu, and B. Zheng, “Hopf bifurcation in numerical approximation of a class delay differential equations,” Applied Mathematics and Computation, vol. 146, no. 2-3, pp. 335–349, 2003.
16. Y. A. Kuznetsov, Elements of Applied Bifurcation Theory, vol. 112 of Applied Mathematical Sciences, Springer, New York, NY, USA, 1995.
17. J. Dieudonné, Foundations of Modern Analysis, Pure and Applied Mathematics, Vol. X, Academic Press, New York, NY, USA, 1960.
18. B. D. Hassard, N. D. Kazarinoff, and Y. H. Wan, Theory and Applications of Hopf Bifurcation, vol. 41 of London Mathematical Society Lecture Note Series, Cambridge University Press, Cambridge, UK, 1981.
19. V. Wulf, Numerical Analysis of Delay Differential Equations Undergoing a Hopf Bifurcation, University of Liverpool, 1999.
20. S. Wiggins, Introduction to Applied Nonlinear Dynamical Systems and Chaos, vol. 2 of Texts in Applied Mathematics, Springer, New York, NY, USA, 1990.