Mathematical Problems in Engineering

Volume 2013, Article ID 852082, 18 pages

http://dx.doi.org/10.1155/2013/852082

## Two-Dimensional Convolution Algorithm for Modelling Multiservice Networks with Overflow Traffic

Poznan University of Technology, Polanka 3, 60-965 Poznan, Poland

Received 10 March 2013; Accepted 7 May 2013

Academic Editor: Jun Jiang

Copyright © 2013 Mariusz Głąbowski 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.

#### Abstract

The present paper proposes a new method for analytical modelling multiservice networks with implemented traffic overflow mechanisms. The basis for the proposed method is a special two-dimensional convolution algorithm that enables determination of the occupancy distribution and the blocking probability in network systems in which traffic streams of individual classes can be serviced by both primary and alternative resources. The algorithm worked out by the authors makes it possible to model systems with any type of traffic offered to primary resources. In order to estimate the accuracy of the proposed method, the analytical results of blocking probabilities in selected networks with traffic overflow have been compared with simulation data.

#### 1. Introduction

The basic structure of telecommunications networks is the hierarchical topology (a.k.a. tree topology). The application of the hierarchical topology is followed by a reduction in financial outlays in constructing networks and makes the effective management of their resources possible. Initially, networks with hierarchical structures were public switched telephone networks (PSTN) [1], based on the application of the circuit switching technology. In the 1950s, the hierarchical structure applied to the PSTN network made it possible to introduce a strategy for diverting traffic via alternative routes [2–5]. To achieve that, the following types of resources were distinguished in the hierarchical structure of a network (defined in systems with circuit switching as *link groups*): (i)direct group: a group of links defined directly between end nodes, that is, ingress and egress nodes, (ii)transit group: a group of links between the end node and a transit node or between transit nodes, (iii)basic group: a transit group that connects any switching node with its superordinate or subordinate node, or one that leads between nodes of the highest hierarchy.

The assumption was at the time that both direct groups and transit groups could be high-usage groups or groups with minimal loss. Basic groups constitute the network backbone of a PSTN telecommunications network—the so-called core network. The use of these groups decreases the overall cost of the network and enables operators to extend their transmission capabilities without a simultaneous increase in the capacities of nodes and with the assumed traffic loss factor retained. Basic groups are generally groups with relatively minimal losses, whereas direct groups are in the main high-usage groups. According to the traffic management strategy for alternative routes, calls arriving at the nodes of telecommunications networks free resources in direct groups, that are primary groups, are searched. Some part of this traffic, the so-called overflow traffic, that cannot be serviced by a primary group due to its occupancy, is offered to alternative groups. Traffic that is not carried by this group is rejected traffic.

The above-discussed method for traffic management via alternative routes was then applied to wireless networks: 2G (e.g., GSM—Global System for Mobile Communications), 3G (e.g. UMTS—Universal Mobile Telecommunication System) and 4G networks (e.g., LTE—Long Term Evolution LTE). Overflow traffic in wireless networks was applied within the framework of the same technology and between networks operating in different technologies (e.g. GSM, UMTS, and LTE) [6–8].

In present-day market place, the growing variety of different access technologies—both wireless and wired networks—is accompanied by a concurrent unification with regard to network layer protocols: the standard for the network layer is the IP protocol (v4, v6). At the same time, the increasing number of network service has eventually led to a necessity to ensure predefined quality of service (QoS) parameters demanded by traffic sources in a network. Defined quality of service parameters in networks with packet switching is guaranteed following creation of virtual channels (with dedicated or shared resources) that transmit packet streams related to a given class (or classes) of service, which leads to a connection-oriented packet-switching architecture communications. The allocation of packet streams to given classes of services is executed (in edge routers) as a result of an appropriate determination of the field of the type (class) of service in the IP header. This identification field can then be used as a DSCP (Differentiated Services Code Point) code point in packet marking in the case of the DiffServ (Differentiated Services) architecture. In the case of the MPLS technology (MultiProtocol Label Switching), the most widely used by operators in backbone networks, the field determining the type of class is mapped into a value of a corresponding label that unequivocally determines the path in the network and/or the EXP (experimental) parameter used for distinguishing quality parameters for packet streams transmitted over the same path. The parameters that can be used for a classification of packet streams are, for example, the interface of the edge router (ingress router), the address of the destination network, and the relationship of the source of traffic stream with a given VPN network (virtual private networks) of the second or the third layer [9–11].

From the traffic engineering perspective, the resources (virtual channels, MPLS paths) that are used to service individual packet streams (traffic) can be treated as the *primary resources*. At the same time, in order to optimize the resources of a network and to secure the demanded quality of service, prevent the network from being overloaded, as well as for network survivability, additional resources (such as additional MPLS paths, additional resources within MPLS path) can be introduced to the connection-oriented packet networks in which packets lost as a result of the service in dedicated primary resources would be transmitted. Such an additional resource (virtual channel, MPLS path) that services blocked traffic in the primary resources can be treated as the *alternative resources* that service overflow traffic.

The major difficulty that occurs during an analysis of systems with traffic overflow is to determine the demanded volume of alternative resources (with low losses). If we assume that a given distribution of time between calls (a call in packet networks can be interpreted as a session, e.g., TCP or UDP session) is offered to the primary resources (e.g., exponential distribution), then traffic that overflows from these resources will be of a different nature [12], because calls from an overflow stream can appear only during the total occupancy of the primary resources. This means that an overflow stream is more “concentrated/dense” within certain time intervals, that is, having a dynamic and “peak” nature as compared to traffic offered to primary resources. If we assume identical value of offered traffic and identical value of the blocking probability, then to execute service for overflow traffic, a greater number of resources is required than that for the service of traffic offered to the primary resources.

In traffic engineering, in works devoted to modelling of systems with traffic overflow, it has been usually assumed that the call stream offered to primary resources is consistent with the Poisson distribution or with binomial distribution [8, 12–19]. The variety of network services to be observed in modern-day networks results, however, in a situation where the previously adopted assumptions do not describe call streams generated in packet networks precisely enough. To make an analytical modelling of the network with any distribution of call streams offered to primary resources possible, this paper proposes a new method for modelling multiservice packet networks with implemented traffic overflow mechanisms. The basis of the proposed method is a special two dimensional convolution algorithm. The worked-out method allows researchers to determine the occupancy distribution and the blocking/loss probability in networks with both primary and alternative resources.

The further part of the paper is organized as follows. Section 2 presents an overview of research studies in modelling multiservice networks with traffic overflow. Section 3 proposes a new model of the overflow system with multirate traffic. In Section 4, the results of the analytical calculations are compared with the results of the simulation. Section 5 sums up the paper’s main conclusions.

#### 2. Related Works

In traffic theory both single-service systems [12, 14, 16] and multiservice systems [17, 20] have been considered. In the case of single-service systems, all calls that are serviced by a telecommunications system always demand the same amount of resources. In the case of multiservice systems, in which multirate models are used for modelling, offered traffic is a mixture of different classes of calls, each demanding an integer multiplicity of a certain unit of the resources, called basic bandwidth unit (BBU) [21–28]. In the case of constant bit rate sources, resources are expressed in bit rates of original traffic streams. In the case of variable bit rate sources prior the BBU determination, the so-called *equivalent bandwidth* [27] for particular traffic packet streams is calculated (Methods for determination of equivalent bandwidth depend on such parameters such as admissible packet loss rate, admissible latency, link capacity, the average and the maximum value of bit rate, the nature of packet streams (e.g., self-similar streams), and the type of the network. Algorithms for determination of the equivalent bandwidth for defined types of the network and services are proposed, for example, in [27, 29–36].)

For single-rate traffic, two basic methods for determination of blocking probabilities in alternative groups (groups to which traffic overflows from other groups) have been worked-out: the equivalent random technique (ERT) method [12, 14] and the Fredericks-Hayward method [16]. The latter method has become the basis for a generalization for multirate overflow traffic [17]. In this method [17], probabilities of subsequent states were determined in a recurrent way using the generalization of Kaufman-Roberts formula [37, 38].

A different technique for modelling networks with traffic overflow is proposed in [39]. The basis for the method proposed in [39] is the Erlang’s Ideal Grading model with multirate traffic [40, 41]. The objective of the proposed method is to simplify the process of the determination of the occupancy distribution in systems with traffic overflow as it does not require calculations for the parameters of overflow traffic. Its accuracy is comparable to the accuracy of the method [17]. The basic limitation of this method is that only one class of calls is offered to each of the primary resources.

In calculations of multiservice systems with traffic overflow, the so-called convolution algorithms [20, 22, 42, 43] can also be used. The advantage of these algorithms is that they offer a possibility to model systems with any streams of offered traffic. References [22, 42, 43] propose convolution models of systems with reservation. In [44], the authors of the paper describe, for the first time, the problem of modelling systems with overflow traffic with the application of the convolution mechanism. The presented method is characterized by high accuracy, though the scope of its application is limited to systems in which each of the primary resources services directly calls of only one traffic class. This paper proposes, using the concept idea presented in [44], a new, generalized method that makes it possible to determine characteristics of multiservice systems with traffic overflow in which each primary group is offered a number of traffic streams. To improve the order of computational complexity of the new method, the simultaneous convolution operation (proposed in [44]) has been replaced by a two-dimensional convolution operation worked-out for the purposes of this paper (Section 3).

#### 3. Modelling of Multiservice Systems with Traffic Overflow

##### 3.1. Basic Assumptions

In order to present the basic assumptions of the proposed method for modelling systems with overflow traffic, let us consider an overflow system that consists of primary groups that belong to the set and one alternative group (Figure 1). Each primary group () is offered traffic classes from the set with cardinality . Traffic of class , the calls of which demand BBUs, is offered to the primary group . Let denote a set of all offered classes: A primary group with the capacity services traffic streams with the intensity . Discarded traffic that is not serviced in these groups forms an overflow traffic stream. This traffic is then offered to the alternative group with the capacity . Hence, calls of classes of the set can occupy at the maximum BBUs in the system: The calls of class can occupy at the maximum BBUs in the system: The total capacity of the system is

##### 3.2. Occupancy Distribution

Firstly, let us consider the occupancy distribution in the basic subsystem of the considered overflow system presented in Figure 1, that is, a system with one primary group and one alternative group for calls of classes of the set . A single th element of the distribution , that is, , denotes the probability of busy BBUs. The maximum number of BBUs that can be occupied in such a system is (equation (2)).

By expanding our analysis of the system into a subsystem that is composed of two primary groups “” and “”, we can determine the occupancy distribution for two sets of service classes and . The single element of the occupancy distribution in the system with two primary groups and can be obtained by a convolution operation of distributions and : where and are determined by (2), in which determines a new length of the distribution equal to The parameter in (5) is the normalization coefficient. Notice that according to the convolution operation defined by (5), the distribution is shortened to the length that is lower than . After shortening of the distribution, it is necessary then to normalize it. The normalization coefficient ensures the sum of all elements of the distribution to be equal to one:

Let us notice that the expression (5) does not precisely reflect the operation of the overflow system since the convolution operation does not take into consideration states that appear as a result of the termination of a service of certain calls. States that occur directly after the termination of the service of a call in the primary group, in instances where all BBUs of the primary group and some (or all) BBUs of the alternative group were busy before the termination of this service, are not taken into consideration properly. To present this phenomenon, let us consider a simple example of an overflow system that is composed of three primary groups , , and with the capacities , , and BBUs. Each group is offered one traffic class with demands equal to 1 BBU. Traffic coming from the three groups overflows to alternative group with the capacity BBUs. The total capacity of this system is then equal to BBUs. Let us consider now the combination . The adopted notation is such that calls of the first class, offered to the first group, occupy 4 BBUs and calls of the second class, offered to the second group, occupy 8 BBUs, whereas no call that was offered to the third group is currently being serviced. The combination is feasible because calls of the first class occupy the entire primary group 1 (4 BBUs). Then, calls of class can occupy 8 BBUs (entire primary group 2 (4 BBUs) and the whole alternative group (4 BBUs)) because calls of class are not serviced at all. The considered combination is shown in Figure 2(a).

Note that the method for determination of the occupancy distribution in the system with traffic overflow presented in [44] assumes that after termination of each of the calls serviced in the primary group, calls are transferred from the alternative group to the appropriate primary group provided that the latter has free resources. This means that after a disconnection of a call of class 2 in the primary group 2 in the considered state , the system transits to state (Figure 2(b)) in which immediately a transfer of the serviced connection of class 2 in the alternative group to the primary group is executed (Figure 2(c)). Consequently, the combination means the occupancy of 4 BBUs by calls of class 1 in the primary group 1 and 7 BBUs by calls of class 2; while all 4 BBUs are busy in the primary group, 2 and 3 BBUs are busy in the alternative group. Hence, the convolution (5) determines the occupancy distribution in the overflow system that includes a transfer of connections from alternative group to primary group.

Subsequently, let us try to determine, for the considered exemplary system, the aggregated distribution for the case when the system simultaneously services calls that are offered to all of the three primary groups. This distribution cannot be determined in a direct way since not all of the combinations included and taken into consideration in the distribution (5) are allowable. Consider an exemplary combination in which calls offered to primary group 1 occupy BBUs, calls offered to primary group 2 occupy BBUs, and calls offered to primary group 3 occupy BBUs. In total, there are BBUs occupied for such a combination, less than the total capacity of the system . Note, however, that according to the combination under consideration calls offered initially to primary group occupy BBUs in the alternative group, whereas calls offered initially to primary group occupy BBU. This means that calls offered initially to primary group and “occupy” more resources of the alternative group than the number of resources available in the group. Therefore, while determining probabilities for each of the states, it is necessary to omit those combinations that lead to an occupancy state in the alternative group higher than its capacity.

The method that assumes call transfers between the alternative and the primary group and eliminates all forbidden states in which the number of busy BBUs in the alternative group exceeds the capacity of this group, the so-called simultaneous convolution operation of many distributions, is proposed in [44]. The method is characterized by high accuracy, while its limitation involves the assumption of only one call class that is offered to each of the primary groups. In Section 3.3, a new method will be proposed. The method is based on the so-called two-dimensional convolution operation that, due to a lower order of complexity, can make it possible to model more effectively systems with overflow traffic in which the primary group services single classes of calls. This will be followed by a presentation of two generalized methods (continuous method—Section 3.5.1—and discrete method—Section 3.5.2) that enable modelling systems in which each primary group can be offered many classes of calls.

##### 3.3. Two-Dimensional Distribution

Let state denote a state in which there are occupied resources of the alternative group, with the total occupancy of the system equal to . The two-dimensional occupancy distribution in this system will be denoted by the symbol . In this distribution a single element determines the state probability .

According to the adopted notation in the paper, the first dimension determines the occupancy state of the alternative group, whereas the second dimension determines the occupancy state of the whole system. To justify such an approach let us consider an exemplary relation between the one-dimensional distribution and the two-dimensional distribution . Both distributions describe state probabilities for the system with traffic overflow that is composed of the primary group with the capacity and the alternative group with the capacity . The primary group is offered calls from the set . The relations between the elements of both distributions can be written in the following way: Figure 3 shows the interpretation of (8) for the system composed of the primary group with the number , with the capacity and the alternative group with the capacity . The primary group is offered calls from the set . The figure shows only those states——of the distribution that are permitted by the system. Thus, the states in which a larger amount of resources of the alternative group than the capacity of this group would be occupied (), as well as the states in which more resources of the alternative group than the total amount of occupied resources would be occupied (), are omitted.

##### 3.4. Modelling of Systems in Which the Primary Group Is Offered a Single Class of Calls

Let us recall that the method proposed in [44] is limited to the instances of modelling of systems in which the primary group is offered only one class of calls. The method [44] can be optimized effectively with the application of the two-dimensional distribution presented in the further part of this subsection. The input data that enable the determination of the two-dimensional distribution are then the occupancy distributions of the so-called subsystems. A subsystem is composed of a set of a certain number of primary groups and an alternative group. In this paper, the basic subsystem will be a system that is composed of the primary group and an alternative group. In our further considerations, the assumption that the system with traffic overflow is approximated by a system with calls transfer still holds good.

###### 3.4.1. Determination of the Two-Dimensional Occupancy Distribution of the Basic Subsystem

Using the two-dimensional distribution , it is possible to describe the characteristics of the occupancy for the basic subsystem . A single element of this distribution determines the occupancy probability of BBUs in the alternative group and BBUs in the whole system by calls of the classes from the set . In the case of the basic subsystem the set is a one-element set. To improve the readability of the notation, we adopt that the symbols and will always denote occupancy distributions for one class of calls exclusively.

In the basic subsystem , which is offered one call of class , the distribution is the equivalent of the distribution . In a similar way, the distribution is the equivalent of the distribution . A single element of the two-dimensional distribution can be determined on the basis of a single element of the one-dimensional distribution for the group with the capacity . For this purpose, we adopt that for the occupancy states , where denotes the total number of busy (occupied) resources, the relation is fulfilled. The parameter denotes the number of busy resources of the alternative group. The adopted assumption is fulfilled for the systems with calls’ transfer and allows for the occurrence of a state in which the resources of the alternative group are busy, while part of the resources of the primary group has not been entirely used. Unused resources of the primary group result from the necessity of servicing a single call requiring within one group.

Thus according to adopted assumptions, the relation between the expressions of both distributions (one-dimensional and two-dimensional) can then be written in the following way:

Following the assumption for transferring the calls, it is possible to determine in an easy way the two-dimensional distribution of the occupancy in the basic system to which one class of calls is offered. The following subsection will present a description of the method for a determination of the two-dimensional occupancy distribution for multiservice systems with traffic overflow. The method (Section 3.4.2) is universal and does not depend on the number of traffic classes offered to basic subsystems. The distributions determined in the following subsection will provide a starting point for the algorithms for modelling systems in which many classes of calls are offered to the primary group.

###### 3.4.2. Two-Dimensional Convolution Operation

Consider now a method for a determination of the two-dimensional occupancy distribution for a multiservice system with traffic overflow. Each expression of this distribution determines the occupancy probability BBU of the alternative group and BBU of all groups (including the alternative group). The two-dimensional convolution operation makes it possible to effectively determine the distribution for the whole system with traffic overflow on the basis of the occupancy distributions of basic subsystems with traffic overflow. The distribution can be a distribution for a subsystem in which the group is offered one class of calls (its determination is described in Section 3.4.1) or a subsystem in which the primary group is offered a number of call classes (two different methods for a determination of these distributions are presented in Section 3.5.1 and Section 3.5.2, resp.).

In order to determine the occupancy distribution for a system composed of primary groups that belong to the set , all distributions () are to be aggregated. The distribution , which determines state probabilities for a subsystem composed of primary groups and the alternative group, is determined on the basis of the convolution operation formulated as follows: for and . The symbol in formula (10) is the normalization coefficient that is determined as follows:

Let us interpret (10) with the example of a system that is composed of primary groups with the capacity BBUs each and an alternative group with the capacity BBU. We will be considering combinations of an overflow system composed of three primary groups and one alternative group with the following capacities: . Let us consider then a combination for the respective subsystems. The considered combination describes a state in which calls of the classes from the set occupy BBUs in a primary group and do not occupy resources in the alternative group of the subsystem. Calls of the classes from the set occupy BBUs in the primary group and BBUs in the alternative group of the subsystem, whereas calls of the classes from the set occupy BBUs in the primary group and BBU in the alternative group of the subsystem. Thus, the total number of busy BBUs in the alternative group is equal to , which means that this combination is not applicable. Let us analyse now a method for determination of the probability of an occurrence of such a combination using the two-dimensional convolution operation.

Note that after the classes from the set and the classes from the set are aggregated, the probabilities of the occurrence of the state combination will be added up to the probability of the permitted state . While proceeding with further aggregation of the distribution of the classes from the set with the distribution of the classes from the set , the probability of occurrence of the state combination should be added to the state probability . This state, however, is not a permitted state in view of the convolution operation determined in (10), which allows us to eliminate properly nonpermitted states in the process of the determination of the two-dimensional distribution for the system with traffic overflow.

Note that the convolution operation increases the length of the distribution in the dimension that determines the total number of busy BBUs, while the dimension that defines the capacity of the alternative group remains without any changes.

The introduced notation, in which one dimension determines the occupancy of the alternative group while the other determines the occupancy of primary groups and the alternative group, makes it possible to directly convolute two-dimensional distributions. The order of complexity of a single convolution operation of two-dimensional distributions is equal to , whereas the order of computational complexity of the algorithm for primary groups is , which is a major advantage over the method discussed in [44], where the order of complexity is .

###### 3.4.3. Blocking Probability

On the basis of the convolution operation (10) it is possible to determine the aggregated occupancy distributions for all groups except the group : Distribution (12) will allow us to work out a method for determination of the blocking probability for calls of traffic classes offered to group , when the alternative group is offered overflow traffic from all primary groups.

For this purpose, let us consider now the feasible combinations of the distributions and . Let the set be a set of all permitted occupancy combinations for the basic subsystem with a subsystem that is composed of primary groups that belong to the set . Let us define the possible combinations of the occupancy of the basic subsystem with the subsystem that is composed of the remaining primary groups and the alternative group: The condition is self-explanatory since the total number of BBUs serviced in the system cannot exceed the capacity of the system. The condition in formula (13) means that the number of all busy BBUs in the alternative group will always be equal to or at the maximum, if . The condition means that the total number of busy BBUs occupied by calls offered to the group cannot be higher than the availability of the basic subsystem or than the total number of busy BBUs in the system (if ). The condition , in turn, expresses the fact the number of busy BBUs in the alternative group occupied by calls that overflow from the group will never be higher than the total number of the resources that occupy calls offered to the group or than the total number of occupied resources of the alternative group, if .

In the same way we will define the set that is a subset of the set and that determines the blocking state for the class . The basis for the determination of this set is the lack of free BBUs in both the alternative and the primary group . Therefore we get

The probability of the combination from the set , in which calls of all classes occupy BBUs in the primary groups and in the alternative group, while calls of the classes from the set occupy BBUs in the primary group and the alternative group ( BBUs in the alternative group), can be determined on the basis of the distributions and : where is the normalization coefficient equal to The blocking probability for the class can be determined on the basis of the sum of the probabilities of blocking combinations from the set :

###### 3.4.4. Algorithm for Modelling Systems with Overflow for Primary Groups with Single Classes of Calls

The algorithm that allows us to determine the blocking probability in networks with traffic overflow in which each primary group is offered a single class of calls can be written in the form of the following four stages: (1)determination of the one-dimensional distribution for a single class offered to the primary group on the basis of models worked-out for single-service systems, (2)determination of two-dimensional distributions for all primary groups on the basis of one-dimensional distributions (formula (9)), (3)determination, for every primary group , of aggregated distributions of all primary groups except group (formula (12)), (4)determination of the blocking probability (formula (17)) for each traffic class offered to each primary group .

##### 3.5. Modelling of Systems in Which the Primary Group Is Offered a Number of Call Classes

###### 3.5.1. Continuous Method

Consider now an overflow system in which primary groups are offered many classes of calls (multiple classes of calls). According to the adopted notation, the primary group is offered call classes from the set . In the first proposed method, that is, the continuous method, we adopt the assumption that a single call can be serviced simultaneously by the resources of the primary as well as the alternative group. The occupancy distribution in the considered basic subsystem is determined on the basis of a modified convolution operation of two-dimensional distributions.

Consider a system in which the primary group has the capacity . Traffic from this group overflows to the alternative group with the capacity . Assume that the sets and are separable subsets of the set . The primary group is offered mixtures of call classes that belong to the sets and for which the distributions and have been determined.

For any group in the overflow system with the assumed calls’ transfer, the two-dimensional convolution operation makes it possible to determine the elements of the distribution on the basis of the summation of probabilities (combinations) of states that satisfy the following condition: Condition (18) determines all combinations for states whose probabilities can be added up to the probability of the aggregated state . For the states that satisfy the condition (18), the following relations then ensue:

Let us determine now for the state a set of all state combinations for which the condition (18) is fulfilled, that is, for which the value of deficit (described in the next paragraph) is equal to 0:

The system also offers such combinations of states that do not satisfy the condition (18). These combinations will be permitted in the system when part of the serviced calls is transferred to the alternative group. Then, the relation (19) that determines the total occupancy in the system is still retained, while the relation (20) is not satisfied because the occupancy of the alternative group changes. The failure of the relation (20) results from the so-called deficit in the resources of the primary group. This deficit can be determined in the following way: The interpretation of the deficit of the primary resources is shown in Figure 4. Figure 4(a) presents the combination , when the resources of the primary and the alternative groups are occupied only by calls from the set and when the resources are occupied only by calls from the set . Figure 4(b) shows the attempt to aggregate both combinations. It is clear to see the deficit in the resources that is to be “transferred” to the resources of the alternative group (Figure 4(c)). Note that the amount of all resources of the alternative group after the transfer of the deficit resources taken into consideration is equal to Taking into account (19), which is always true (for and for ), we can express (23) in the following form: Now, we are in a position to define the set of all permitted combinations ; for which the deficit value is higher than zero (the condition (18) has not been fulfilled): Taking into account (22), the set (25) can be rewritten as follows:

On the basis of the sets and it is possible to determine the sum of the probabilities of all combinations for the state , that is, the probability for this state:

Equation (27) can be rewritten in the following way: The transformation of (27) into (28) is presented in Appendix A. Using (28) we are in a position to aggregate subsequent distributions in the primary group. It is thus a convenient notation for a construction of the computational algorithm.

###### 3.5.2. Discrete Method

The method for determination of the occupancy distribution in the basic subsystem with multiple call classes, that is, a system that is composed of the primary group and the alternative group presented in Section 3.5.1, adopts the assumption that the transfer of calls results in a total exploitation of the primary group. Such an approach assumes a possibility of simultaneous occupation of the resources of the primary and the alternative group by one call in the stage of the aggregation of classes within the basic subsystem and, in consequence, leads to inaccurate and imprecise determination of the deficit in resources.

Let us consider now a modified method for determination of the two-dimensional occupancy distribution for the basic subsystem on the basis of the occupancy distributions and in which a possibility of a simultaneous servicing of calls by the resources of the primary and the alternative groups is excluded. Note that, on the basis of the given combination , it is not possible to unequivocally determine the history of call acceptance and call service. It is not possible then to determine the amount of resources of the primary group that cannot be used due to the lack of the possibility of a division of the serviced call between the primary group and the alternative group. The value of the parameter depends not only on the combination itself but also on the history of call acceptance and call service.

As previously mentioned, let us assume that the basic subsystem is offered mixtures of call classes that belong to the set and . The sets and are separate subsets of the set traffic classes.

The number of occupied resources in the alternative group can be determined in the same way as in Section 3.5.1 but with the number of unused resources in the primary group taken additionally into account. If the total number of occupied resources is lower than or equal to the capacity of the primary group, then there is no deficit in resources and all resources of the primary group can be used. In this situation, the condition (18) is also fulfilled. If the number of occupied resources exceeds the capacity of the primary group and the condition (18) is at the same time fulfilled, then there is no deficit in the primary group. When this is the case, the number of busy resources of the alternative group is higher by the number of unused resources of the primary group . In the case when the condition (18) is not fulfilled, the number of busy resources of the primary group is increased by the parameters and . We can thus write

Consider now the way the parameter that determines the number of unused resources of the primary group is estimated. Assuming that there is a possibility of calls’ transfer, there can be at the maximum unused BBU or unused BBU in the primary group, depending on which call class has been admitted for service as the last one. The upper boundary for the value of the parameter is BBU, that is, the number of free resources in both groups (primary and alternative). Let denote the maximum number of unused resources in the primary group, on condition that a call that belongs to the set has been accepted as the last one:

Let denote the probability of BBU unused in the primary group in state , on condition that the last accepted call was a call of the class that belongs to the set . A determination of this probability is fairly complex. Let us assume then that the distributions and are approximated by the uniform distribution:

Let denote the probability that the state has been reached when the last accepted call was a call of the class that belonged to the set , and denotes the probability that the state has been reached when the last accepted call was a call of the class that belonged to the set . The proposed method makes the assumption that these probabilities are directly proportional to the number of resources that are occupied by calls that belong to the sets and :

Let us define the convolution of the distributions and . To determine the state probability , let us consider the three intervals defined in (29). The first interval includes the states for . The second and the third intervals include states in which the total number of occupied resources is higher than the capacity of the primary group . Each interval includes the combinations of states denoted as .

Let us consider the method for determination of the probability of occurrence of the states for as the first. The probability of the occurrence of the state is equal to the sum of the probability of occurrence for all combinations for this particular state, defined in (19) and (29): For the remaining two intervals defined in (29), to determine the state probability , for , it is necessary to take into consideration additional parameters such as , (formula (31)) and and (formula (32)). In addition, it is necessary to define the set of all possible combinations for the second interval of (29) and the set of all possible combinations for the third interval.

Definitions of the sets and depend on the last admitted call. If the last admitted call is a call of the classes from the set , then for the state we denote these sets as and . Similarly, if the last admitted call was a call of the class that belonged to the set , then the sets will be denoted as and . The sets , , and define all the possible combinations as well as the number of unused resources for the state . In addition, the sets and define the deficit of the resources of the primary group. As a consequence of our considerations we can write

In line with the adopted notation, the state probability (for ) can be written as follows: Equations (33) and (35) can be combined and eventually rewritten in the following form: The way (33) is transformed to the form (36) is described in Appendix B.

#### 4. Numerical Results

To confirm the accuracy and effectiveness of the proposed convolution method, the results of the blocking probability in the overflow systems obtained on the basis of the analytical calculations, both for the continuous method (Section 3.5.1), and the discrete method (Section 3.5.2), were compared with the results of the simulation experiments. Additionally, in the case of call streams generated by infinite number of traffic sources, the results obtained on the basis of a modified Hayward method [17] are presented. The characteristics presented in the graphs show the dependence between the blocking probability of each of the offered classes in the whole system and offered traffic , offered to a single BBU of the primary group:

The event-oriented discrete simulation method was used for simulation [45, 46]. Each simulation experiment consisted of 10 series. The length of each series was determined by the number of lost calls (at least 200000 lost calls for each traffic class). The 99% confidence intervals were calculated according to the Student-Fisher distribution.

The study was carried out for the telecommunications systems defined in Table 1. The systems are offered three types of Erlang (Poisson call streams), Engset (binomial call streams), and Pascal (negative binomial call streams) traffic streams [20, 47]. The selected types of traffic cover three different types of the dependence between the mean arrival rates of calls and the occupancy state of the system: the mean arrival rate of new calls does not depend on the occupancy state of the system (Erlang traffic), the mean arrival rate of new calls decreases with the increase in the occupancy state of the system (Engset traffic), and the mean arrival rate of new calls increases with the increase in the occupancy state of the system (Pascal traffic).

The graphs presented in Figures 5, 6, 7, 8, and 9 showing the blocking probability, determined on the basis of the two-dimensional convolution algorithm, both continuous and discrete, as well as the results of the simulation experiments, are presented for each of the systems under consideration. In addition, for the systems and , the results obtained on the basis of the generalized Hayward method are included. The confidential interval in the graphs is too small to be visible due to the length and the number of series in the simulation.

We can notice that both convolution algorithms, continuous and discrete, ensure higher accuracy than the modified Hayward method. Additionally, we can observe, that the continuous algorithm lowers probabilities of states that belong to the blocking area for classes that demand the highest number of BBUs in relation to the remaining classes. This phenomenon can be justified in the following way. The continuous method assumes a possibility of servicing a single call by the primary and the alternative group. In this way, the model assumes a better use of resources than it is actually the case in a real system. As a consequence of the adoption of this assumption, states in which some resources of primary groups would not be busy do not occur. This leads to an increase in the number of free resources of the alternative group and, in consequence, to an increase in the number of serviced calls that require a higher number of BBUs. At the same time, elimination of states in which the primary group would have free resources leads to an increase in the blocking state for calls that demand a lower number of resources. As a result, the probabilities determined for classes that demand a higher number of BBUs are lower than those obtained in real systems, while blocking probabilities for classes that demand a low number of BBUs are overestimated. The highest accuracy of calculation process is ensured by the discrete algorithm.

#### 5. Conclusions

The paper proposes new methods for modelling of network systems with traffic overflow. The methods are based on the proposed two-dimensional convolution algorithm and can be applied to determine the parameters of the overflow system exclusively on the basis of the knowledge of the parameters of offered traffic and the volume of the resources. Additionally, the method makes it possible to determine the blocking probability in systems that are offered traffic streams of any type, which is an important advantage of the proposed method. The proposed methods differ considerably in the way two-dimensional distributions that describe the occupancy of the primary group and the alternative group are determined.

The first method for determination of the two-dimensional convolution, the so-called continuous method, gives the same results as the simultaneous convolution method [44]. The proposed algorithm offers, however, a lower order of computational complexity and can be used for modelling systems in which the primary group is offered multiple call classes. The other method, the so-called discrete method for determination of the two-dimensional convolution, gives more accurate results with the same order of computational complexity retained.

It should be emphasized that the methods for modelling multirate systems with traffic overflow developed in [17] aimed at determining the blocking probability only in alternative resources (groups). Moreover, all hitherto known methods for modelling systems with both single-service overflow traffic (ERT method and Hayward method), and with multiservice overflow traffic [17], require the overflow traffic parameters to be first determined in the process of a determination of losses in the overflow system. Consequently, these methods are limited to modelling only systems with Erlang and Engset traffic streams, whereas the method proposed in the present article makes it possible to model systems with multiservice overflow traffic generated by sources of any type.

#### Appendices

#### A.

Let us consider again (27): Equation (A.1) can be expressed in the form of the following functional dependency: where the parameter determines the product convolution depending on four parameters: the occupancy of the alternative group , and the occupancy of the whole system , by the calls of the classes from sets and :

Taking into account the conditions determining set in definition (21), we get Therefore the convolution defined by Function (A.3) can be presented in the form of a double sum with regard to the two parameters and only:

By analyzing the conditions defining the set in Definition (25), Function can be made dependent on the three parameters: The notation (A.6) allows us to express the second sum in (A.2) as a triple sum: Observe that the parameter can, at the maximum, take on the value equal to defined by formula (24). This occurs when the total resources serviced in an alternative group result from the deficit (i.e., when ).

In line with (A.5) and (A.7), the sums and can be expressed as follows: Observe that the case for the first equation in formula (A.8) can be taken into consideration in the second equation in formula (A.8) as a result of the exchange of the lower bound of the sum to . Eventually, (A.8) can be written in the following way:

According to (27) and (A.9), we get

#### B.

Let us consider again (35): Equation (B.1) will be written in a more convenient way by replacing the product with the function . Then, to this modified formula we substitute (32) and we get Now we transform (B.2) into a form that enables its simple application in engineering calculations. Observe that the determinacy interval for the parameter , that is, the number of unused resources in the primary group, has been determined in the same way in both the probability definitions , (formula (31)) and in the definitions of the sets , , , (formula (34)). Hence, the determinacy interval of the parameter can be removed from the definition of the appropriate sets and with no influence upon the distribution (B.2). Thus, the sets , can be redefined in the following way: Similarly, the sets , can be rewritten as follows: The definitions for the sets and , according to (B.3) and (B.5), make it possible to simplify formula (B.2):

In (B.6), the resulting sum has been divided into two sums over the separable sets and . Let us consider first, then, the first sum in formula (B.6). Notice that according to the definition of the set (B.3), the pair of the parameters and unequivocally determines one and only one value : Taking into account all the conditions determining set in Definition (B.3), we are in a position to write the following equation:

Therefore, the convolution defined by Function (B.8), that is, the first sum in formula (B.6), can be presented in the form of double sum with regard to two parameters , only: Let us proceed to a transformation of the second sum of (B.6). By analyzing the conditions determining the set in Definition (B.5) we can make Function dependent on the three parameters , , and : The notation (B.10) makes it possible to present the second sum in formula (B.6) in the form of the following triple sum: Now, by adding the expressions (B.9) and (B.11) we can eventually reduce formula (B.6) to the form used in the expression (36) for : Observe that in formula (B.12) the value of the sum for corresponds to the value of the expression (B.9).

#### References

- L. Harte and A. Ofrane,
*Telecom Systems: Pstn, Pbx, Datacom, Ip Telephony, Iptv, Wireless and Billing*, Althos, 2006. - R. J. Harris, “Concepts of optimality in alternate routing networks,” in
*Proceedings of 7th International Teletraffic Congress*, 1973. - Y. Rapp, “The use of computers for network planning,”
*Telecommunication Journal*, vol. 12, 1971. View at Google Scholar - B. Wallstrom, “Methods for optimazing alternative routing networks,”
*Ericsson Technics*, no. 1, 1969. View at Google Scholar - I. Weber, “Some traffic characteristics of communication networks with automatic alternate routing,”
*Bell System Technical Journal*, vol. 41, no. 2, pp. 769–786, 1962. View at Google Scholar - M. Glabowski, M. Stasiak, A. Wisniewski, and P. Zwierzykowski, “Uplink blocking probability calculation for cellular systems with WCDMA radio interface, finite source population and differently loaded neighbouring cells,” in
*Proceedings of the Asia-Pacific Conference on Communications*, pp. 138–142, Perth, Australia, 2005. - H. Holma and A. Toskala,
*WCDMA for UMTS: HSPA Evolution and LTE*, John Wiley & Sons, New York, NY, USA, 2007. - M. Stasiak, M. Glabowski, A. Wiśniewski, and P. Zwierzykowski,
*Modeling and Dimensioning of Mobile Networks*, Wiley, 2011. - L. De Ghein,
*MPLS Fundamentals*, Cisco Press, 2007. - J. Evans and C. Filsfils,
*Deploying IP and MPLS QoS for Multiservice Networks: Theory & Practice*, The Morgan Kaufmann Series in Networking, Elsevier Science, 2010. - C. Lewis and S. Pickavance,
*Selecting Mpls Vpn Services*, Networking Technology Series, Cisco Press, 2006. - R. I. Wilkinson, “Theories for toll traffic engineering in the USA,”
*Bell System Technical Journal*, vol. 35, no. 2, pp. 421–514, 1956. View at Google Scholar · View at MathSciNet - M. Asaduzzaman, T. J. Chaussalet, and N. J. Robertson, “A loss network model with overflow for capacity planning of a neonatal unit,”
*Annals of Operations Research*, vol. 178, pp. 67–76, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - G. Bretschneider, “Die Berechnung von Leitungsgruppen für berfließenden Verkehr in Fernsprechwählanlagen,”
*Nachrichtentechnische Zeitung*, vol. 11, pp. 533–540, 1956. View at Google Scholar - L. E. N. Delbrouck, “On the steady-state distribution in a service facility carrying mixtures of traffic
with different peakedness factors and capacity requirements,”
*IEEE Transactions on Communications*, vol. 31, no. 11, pp. 1209–1211, 1983. View at Google Scholar · View at Scopus - A. Fredericks, “Congestion in blocking systems—a simple approximation technique,”
*The Bell System technical journal*, vol. 59, no. 6, pp. 805–828, 1980. View at Google Scholar · View at Scopus - M. Glabowski, K. Kubasik, and M. Stasiak, “Modeling of systems with overflow multi-rate traffic,”
*Telecommunication Systems*, vol. 37, no. 1–3, pp. 85–96, 2008. View at Publisher · View at Google Scholar · View at Scopus - Q. Huang, K. T. Ko, and V. B. Iversen, “An approximation method for multiservice loss performance in hierarchical networks,” in
*Managing Traffic Performance in Converged Networks, Proceedings of the 20th International Teletraffic Congress (ITC 07)*, L. Mason, T. Drwiega, and J. Yan, Eds., vol. 4516 of*Lecture Notes in Computer Science*, pp. 901–912, Springer, 2007. - J. S. Kaufman and K. M. Rege, “Blocking in a shared resource environment with batched Poisson arrival processes,”
*Performance Evaluation*, vol. 24, no. 4, pp. 249–263, 1996. View at Publisher · View at Google Scholar · View at Scopus - V. Iversen,
*Teletraffic Engineering Handbook*, ITU-D, Study Group 2, Question 16/2, ITC, Geneva, Switzerland, 2005. - M. Glabowski, “Modelling of state-dependent multi-rate systems carrying BPP traffic,”
*Annals of Telecommunications*, vol. 63, no. 7-8, pp. 393–407, 2008. View at Publisher · View at Google Scholar - M. Glabowski, A. Kaliszan, and M. Stasiak, “Modeling product-form state-dependent systems with BPP traffic,”
*Performance Evaluation*, vol. 67, no. 3, pp. 174–197, 2010. View at Publisher · View at Google Scholar · View at Scopus - M. Glabowski and M. Stasiak, “Point-to-point blocking probability in switching networks with reservation,”
*Annales des Télécommunications*, vol. 57, no. 7-8, pp. 798–831, 2002. View at Google Scholar - M. Glabowski and M. Stasiak, “Multi-rate model of the group of separated transmission links of various capacities,” in
*Telecommunications and Networking, Proceedings of the 11th International Conference on Telecommunications (ICT '04)*, J. N. de Souza, P. Dini, and P. Lorenz, Eds., vol. 3124 of*Lecture Notes in Computer Science*, pp. 1101–1106, Springer, Fortaleza, Brazil, August 2004. - M. Glabowski, M. Stasiak, and P. Zwierzykowski, “Modelling of virtual-circuit switching nodes with multicast connections,”
*European Transactions on Telecommunications*, vol. 20, no. 2, pp. 123–137, 2009. View at Publisher · View at Google Scholar · View at Scopus - D. Parniewicz, M. Stasiak, and P. Zwierzykowski, “Multicast connections in mobile networks with embedded threshold mechanism,” in
*Computer Networks, Proceedings of the 18th Conference on Computer Networks (CN '11)*, A. Kwiecien, P. Gaj, and P. Stera, Eds., vol. 160 of*Communications in Computer and Information Science*, pp. 407–416, Springer, Ustron, Poland, June 2011. - J. Roberts, V. Mocci, and I. Virtamo, “Broadband network teletraffic,”
*Final Report of Action COST*242, Commission of the European Communities, Springer, Berlin, Germany, 1996. View at Google Scholar - M. Stasiak, J. Wiewiora, P. Zwierzykowski, and D. Parniewicz, “Analytical model of traffic compression in the UMTS network,” in
*Computer Performance Engineering, Proceedings of the 6th European Performance Engineering Workshop*, J. T. Bradley, Ed., vol. 5652 of*Lecture Notes in Computer Science*, pp. 79–93, Springer, London, UK, July 2009. - S. Bodamer and J. Charzinski, “Evaluation of effective bandwidth schemes for self-similar traffic,” in
*Proceedings of the 13th ITC Specialist Seminar on IP Traffic Measurement, Modeling and Management*, pp. 21.1–21.10, Monterey, Calif, USA, September 2000. - J. Charzinski, “Internet traffic measurment and characterisation results,” in
*Proceedings of the 13th Internetional Symposium on Services and Local Access (ISSLS '00)*, Stockholm, Sweden, June 2000. - N. L. S. Fonseca, G. S. Mayor, and C. A. V. Neto, “On the equivalent bandwidth of self-similar sources,”
*ACM Transactions on Modeling and Computer Simulation*, vol. 10, no. 2, pp. 104–124, 2000. View at Publisher · View at Google Scholar · View at Scopus - R. Gibbens and Y. The, “Critical time and space scales for statistical multiplexing,” in
*Proceedings of the 16th International Teletraffic Congress (ITC '99)*, pp. 87–96, Edinburgh, UK, June 1999. - J. Y. Hui, “Resource allocation in broadband networks,”
*Journal on Selected Areas in Communications*, vol. 6, no. 9, pp. 1598–1608, 1988. View at Google Scholar - F. P. Kelly, “Effective bandwidths at multi-class queues,”
*Queueing Systems*, vol. 9, no. 1-2, pp. 5–15, 1991. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - T. Neame, M. Zuckerman, and R. Addie, “A paractical approach for multimedia traffic modeling,” in
*Proceedings of the 5th International Conference on Broadband Communications*, pp. 73–82, Hong Kong, November 1999. - I. Norros, “On the use of fractional Brownian motion in the theory of connectionless networks,”
*IEEE Journal on Selected Areas in Communications*, vol. 13, no. 6, pp. 953–962, 1995. View at Publisher · View at Google Scholar · View at Scopus - J. S. Kaufman, “Blocking in a shared resource environment,”
*IEEE Transactions on Communications*, vol. 29, no. 10, pp. 1474–1481, 1981. View at Google Scholar · View at Scopus - J. Roberts, “A service system with heterogeneous user requirements—application to multiservice
telecommunications systems,” in
*Proceedings of the the International Conference on Performance of Data Communications Systems and Their Applications*, G. Pujolle, Ed., pp. 423–431, North-Holland, Paris, France, September 1981. - M. Glabowski, S. Hanczewski, and M. Stasiak, “Erlang’s ideal grading in diffserv modelling,” in
*Proceedings of IEEE Africon 2011*, pp. 1–6, (IEEE, Livingstone, Zambia, September 2011. View at Publisher · View at Google Scholar - M. Glabowski, S. Hanczewski, M. Stasiak, and J. Weissenberg, “Modeling erlang’s ideal grading with multi-rate BPP traffic,”
*Mathematical Problems in Engineering*, vol. 2012, Article ID 456910, 35 pages, 2012. View at Publisher · View at Google Scholar - M. Stasiak and S. Hanczewski, “Approximation for multi-service systems with reservation by systems with limited-availability,” in
*Computer Performance Engineering, Proceedings of the 5th European Performance Engineering Workshop (EPEW '08)*, N. Thomas and C. Juiz, Eds., vol. 5261 of*Lecture Notes in Computer Science*, pp. 257–267, Springer, Palma de Mallorca, Spain, September 2008. - M. Glabowski, A. Kaliszan, and M. Stasiak, “Asymmetric convolution algorithm for blocking probability calculation in full-availability group with bandwidth reservation,”
*IET Circuits, Devices and Systems*, vol. 2, no. 1, pp. 87–94, 2008. View at Publisher · View at Google Scholar · View at Scopus - M. Glabowski, A. Kaliszan, and M. Stasiak, “On the application of the asymmetric convolution algorithm in modeling of full-availability group with bandwidth reservation,” in
*Managing Traffic Performance in Converged Networks, Proceedings of the 20th International Teletraffic Congress (ITC '07)*, L. Mason, T. Drwiega, and J. Yan, Eds., vol. 4516 of*Lecture Notes in Computer Science*, pp. 878–889, Springer, Ottawa, Canada, June 2007. - M. Glabowski and A. Kaliszan, “Convolution algorithm for overflow calculation in integrated services networks,” in
*Proceedings of the 17th Asia-Pacific Conference on Communication (APCC '11)*, pp. 428–433, IEEE, Kota Kinabalu, Malaysia, October 2011. - M. Glabowski and A. Kaliszan, “Simulator of full-availability group with bandwidth reservation and
multi-rate bernoulli-poisson-pascal traffic streams,” in
*Proceedings of Eurocon 2007*, pp. 2271–2277, Warszawa, Poland, September 2007. - J. Tyszer,
*Object-Oriented Computer Simulation of Discrete-Event Systems*, Kluwer Academic Publishers, 1999. - H. Akimaru and K. Kawashima,
*Teletraffic: Theory and Application*, Springer, New York, NY, USA, 1999.