1 Introduction
This work focus on the modeling and analysis of a twouser buffered ALOHA type multiple access system by taking into account the presence of signals, which either delete a packet or trigger the instantaneous movement of a packet between the user buffers. More precisely, this work introduce the socalled called Random Access Gnetwork (RAGNetwork) to model virus attacks, as well as probabilistic load balancing in highly interacting multiple access communication.
The exact modelling of the system we consider, is in terms of a two dimensional random walk in the quarter plane, which posses a partial spatial homogeneity, due to the high level of interaction among the queues. The investigation of interacting queueing system is a challenging task with severe mathematical difficulties, which have been initially reported in the seminal work in [1], and further in [2, 3]. They have also received considerable attention due to their applicability in random access networks [4, 5, 6, 7, 8, 9] and in shared processor systems [10, 11, 12, 13]. On the other hand, Gnetworks (or queueing networks with negative customers, signals, triggers, etc.) introduced in [14], are characterized by the following feature: in addition to the conventional customers, negative customers or signals arriving to a nonempty queue remove an amount of work from the queue or transfer it to another queue [15]. The analysis of this versatile class of networks has significantly enriched queueing theory as well as contributed to the development of real applications in fields such as computers, communications, manufacturing, energy as well as selfaware networks.
The goal of this work is to provide a general framework in introducing and analysing multiple access systems with the feature of negative and triggering signals. Our ultimate goal is to provide a general model to describe virus attacks or other malfunctions, and load balancing in multiple access networks as well as to study their impact on the overall system performance.
1.1 Related Work
GNetworks:
Motivated by neural network modelling
[16], a novel stochastic network, called Gnetwork or queueing network with signals was introduced as a unifying model for neural and queuing networks. In contrast to traditional queueing networks where (positive) customers obey the specified service and routing disciplines determined by the network dynamics, there is another type of customers with the effect of signal that interact upon arrival at a queue with the queue or with the backlogged customers.GNetworks [14] establish a versatile class of queueing networks with a computationally efficient product form solution, which had been proved to exist by using new techniques from the theory of fixed point equation [17]. In its simplest version, a signal arriving at a non empty queue forces a positive customer to leave the network immediately [14]. Since their introduction, Gnetworks have been extensively studied covering several extensions such as triggered movement, which redirect customers among the queues [15]; catastrophes or batch service [18], adders [19]; multiple classes of positive customers and signals [20], statedependent service disciplines [21, 22, 23], tandem networks [24, 25], deletion of a random amount of work [26, 27], retrials [28, 29] (not exhaustive list). For a complete bibliography see [30, 31, 32].
Gnetworks have been shown to be a diverse application tool to analyse and optimise the effects of dynamic load balancing in large scale networks [33] as well as in Gene Regulatory Networks [34, 35]. A recent application of GNetworks is to the modelling of systems which operate with intermittent sources of energy, known as Energy Packet Networks [36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46]. Other applications include their of approximating continuous and bounded realvalued functions [47], which serves as the foundation for learning algorithms [48]
and Deep Learning
[49]. They have been also used for modelling natural neuronal networks
[50], in image processing applications [51, 52], and as a tool for predicting the toxicity of chemical compounds [53].In the field of computer network performance, the RNN has been used to build distributed controllers for quality of service routing in packet networks [54, 55] and in the design of Software Defined Network controllers [56]. Realtime optimised task allocation algorithms in Cloud systems [57] have also been built and tested. Recent applications has addressed the use of the RNN to detect attacks on Internet of Things (IoT) gateways [58].
Stability and Delay analysis in Random Access Networks:
The simple and decentralized nature of ALOHA protocol [59] made it very popular in multiple access communication systems. The ever increasing need for massive uncoordinated access has increased the interest on random access protocols [60, 61], which remain an active research area with challenging open problems even for very simple networks [62, 63].
In spite of its simple operation, stable throughput and delay analysis are quite challenging due to the high level of interaction among queues. This is due to the fact that each transmitted source interferes with other users with nonempty queues, i.e., the successful transmission probability is a function of the status of the other users. As a consequence, the departure rate of a queue can be computed only when we know the stationary distribution of the joint queue length process
[64]. This is the reason why the vast majority of previous works has focused on smallsized networks and only bounds or approximations are known for the networks with larger number of sources [65, 64, 5, 66, 67]. Delay analysis of random access networks is even more challenging both by mathematical and the application point of view. The high level of interaction among queues is the reason for the limited number of analytical results. For the sake of clarity we mention the works in [67, 68, 4, 69, 70] (not exhaustive list).Recently, the authors have performed considerable contribution both in the investigation of stable throughput region and in the delay analysis by considering sophisticated queueaware transmission protocols in modern random access schemes [7, 8], in IoT networks [71, 72, 73, 74] and in networklevel cooperative wireless networks [75, 9].
1.2 Contributions
In this work, we consider a twouser slotted ALOHA multiple access channel with collisions, i.e., a transmission is successful if and only if a single user transmits. Each user has external bursty arrivals that are stored in its infinitesized capacity queue and accesses the medium in a random access manner. The major contribution of this work relies on the fact that we consider for the first time in the related literature the concept of signals in random access networks, which introduces the socalled Random Access GNetwork.
In particular, we consider both negative signals, which delete a packet from a user’s queue, and signals that trigger the instantaneous movement of packets among user queues. The concept of negative signals can model virus attacks in such networks or other possible malfunctions, whereas, triggering signals can be used to model load balancing schemes. To the best of our knowledge this variation of random access has not been reported so far. For such a network, we investigate stability conditions, the stable throughput region and the queueing delay.
1.2.1 Stability and Stable Throughput Analysis
The presence of negative signals affects the stable throughput region, which now becomes a subset of the stability region, since packets can be dropped before being transmitted to the destination thus, these packets they do not contribute to the achievable throughput. Note that such property has never been reported in the literature of random access networks so far.
When characterizing the stability and/or stable throughput in such a network, we have to cope with the problem of highly interaction among queues. In RAG network, the level of interaction is even higher due to the presence of signals that interact upon arrivals both with the queue and with the backlogged packets. We show that the stochastic dominant technique [64]
is still an efficient tool to bypass the problem of interaction among queues. Although stability conditions can be also derived by well known methods from two dimensional Markov chains
[76, 77], the importance of stochastic dominant technique is more apparent when we characterize the stable throughput region, where we show that it is a subset of the stability region.1.2.2 Queueing Analysis
Queueing theoretic analysis of interacting queues is a quite challenging task. The increased level of interaction due to the presence of signals will further complicate the analysis. In this work, we present a compact mathematical analysis, and we provide exact expressions for the probability generating function (pgf) of the joint stationary queue length distribution of user queues by solving a nonhomogeneous Riemann boundary value problem [78]. Furthermore, we also provide a computationally efficient alternative way to obtain explicit bounds for basic performance metrics without calling for advanced mathematical concepts of boundary value problems.
The rest of the paper is summarized as follows. In Section 2 we describe in detail the mathematical model, while in Section 3 we provide the stability and the stable throughput region. The fundamental functional equation along with some basic preliminary results regarding its investigation are presented in Section 4. In Section 5 we provide expressions for the pgf of the joint stationary queue length distribution at user queues in terms of a solution of a Riemann boundary value problem, while in Section 6 we provide exact bounds for the expected number of backlogged packets at each user queue, without using the theory of boundary value problems. Numerical results are given in Section 7 and show insights in system performance. Some conclusions and future directions are presented in Section 8.
2 The Mathematical Model
We consider an ALOHAtype random access network (RAN) consisting of two users, say , , that communicate with a common destination node. Each user is equipped with infinite capacity buffers, in which it stores the arriving and backlogged packets. The time is slotted and packets have equal length. The transmission time of a packet corresponds to a single slot.
At the beginning of a slot, user , if it is notempty, transmits a packet with probability , or remains silent with probability . We consider a collision based channel model, and thus, if both users attempt to transmit simultaneously, a collision occurs and both transmissions fail. In such case, the packets have to be retransmitted in a latter slot.
Contrary to the traditional RAN, in this work we introduce a generalized RAN, called the Random Access GNetwork (RAGN) by including the concept of signals. In particular, at the beginning of a slot, signals are generated in user with probability . If a signal is generated, then with probability it deletes a packet from the buffer of , while with probability it triggers the instantaneously movement of a packet from the buffer of to the buffer of the other user. A signal arriving in an empty buffer has no effect. In this work, we assume that if a signal and a packet transmission occur simultaneously, the signal occurs first. In case of a signal generation, a temporarily network malfunction occurs, and both users remain silent during the slot. Moreover, packet arrivals are scheduled at the end of the slot, early departure late arrival model. Denote by the number of packets arriving in , in the time interval , with , . Denote the joint pgf of and by , , , .
Let be the queue length at the buffer of at the beginning of time slot . Under the above system model, the twodimensional process is a Markov chain. Clearly, under usual assumptions is irreducible and aperiodic. Then,

If , then , .

If , ,

If , ,

If , ,
3 Stability Conditions and Stable Throughput Region
From the previous section, it is clear that for a packet located at the head of the queue, there are three options depending if a signal will be generated or not. When a signal is not generated, the packet it will be transmitted successfully to the destination if the node will attempt to transmit and collision will not happen. If a signal is generated, there are two options, either the packet will be dropped from the system or it will be transferred to the other queue. Thus, it is important to emphasize that the stability region is different than the stable throughput region. In fact, the stable throughput region is a subset of the stability region as we will see in this section.
The service probability for is denoted by , , which denotes the probability that a packet will be either transmitted successfully to the destination or it will be relocated to the other queue. The expression for is given by
(3.1) 
Then we define the probability that a packet will be removed from a queue, either it will dropped, or transmitted successfully, or transferred to the other queue. For the users , , we denote this probability . The expression for is given by
(3.2) 
which is equal to
(3.3) 
Similarly, we can write the expression for . Clearly, and depend on the state of the second queue and and on the state of the first queue. Thus, the queues are coupled. We will bypass this difficulty by applying the stochastic dominance technique to obtain the exact stability region and the stable throughput region.
Definition 1.
Denote by the length of queue at the beginning of time slot . The queue is said to be stable if and .
Although we will not make explicit use of this definition we use its corollary consequence which is Loynes’ theorem [80] that states that if the arrival and service processes of a queue are strictly jointly stationary and the average arrival rate is less than the average service rate, then the queue is stable. If the average arrival rate is greater than the average service rate, then the queue is unstable and the value of
approaches infinity almost surely. The stability region of the system is defined as the set of arrival rate vectors
for which the queues in the system are stable. In the previous definition, denotes the total arrival rate at the queue .Before proceeding with the derivation of the stability region and the stable throughput region, we need to calculate the internal arrival probability, for each queue , stands for internal. We have that
(3.4) 
and
(3.5) 
The following theorem provides the stability region for the considered system.
Theorem 3.1 (Stability Region).
The stability region for a fixed transmission probability vector is given by where
(3.6)  
(3.7)  
Proof.
To determine the stability region of our system we apply the stochastic dominance technique [64], i.e. we construct hypothetical dominant systems, in which a node transmits dummy packets for the packet queue that is empty, while for the nonempty queue transmits according to its traffic. Under this approach, we consider the , and dominant systems. In the dominant system, whenever the queue , empties, it continues transmitting a dummy packet. Thus, in , node never empties, and hence, node sees a constant probability that a packet will be removed from its queue, while that probability for node depends on the state of node , i.e., empty or not. We proceed with dominant system . The probability of the first node is given by (3.3). The is given by
(3.8) 
Since in queue never empties, , we have . We can obtain that the total arrival rate at is . By applying Loyne’s criterion [80], the second node is stable if and only if . The probability that the second node is empty and is given by
(3.9) 
Then, we can obtain the internal arrival probability at is given by
(3.10) 
After replacing into (3.3), and applying Loynes criterion [80] we can obtain the stability condition for the first node. Then, we have the stability region given by (3.6). Similarly, we can obtain the stability region for the second dominant system , the proof is omitted. For a detailed treatment of dominant systems please refer to [64].
An important observation made in [64] is that the stability conditions obtained by the stochastic dominance technique are not only sufficient but also necessary for the stability of the original system. The indistinguishability argument [64] applies here as well. Based on the construction of the dominant system, we can see that the queue sizes in the dominant system are always greater than those in the original system, provided they are both initialized to the same value and the arrivals are identical in both systems. Therefore, given , if for some , the queue at the first user is stable in the dominant system, then the corresponding queue in the original system must be stable. Conversely, if for some in the dominant system, the queue at the first node saturates, then it will not transmit dummy packets, and as long as the first user has a packet to transmit, the behavior of the dominant system is identical to that of the original system since dummy packet transmissions are eliminated as we approach the stability boundary. Therefore, the original and the dominant system are indistinguishable at the boundary points. ∎
Theorem 3.2 (Stable Throughput Region).
The stable throughput region for a fixed transmission probability vector is given by where
(3.11)  
(3.12)  
Proof.
Following the same methodology of dominant systems in the previous theorem, we construct two hypothetical systems, and . In , transmits dummy packets when its queue is empty, thus , we have . In the proof of Theorem 3.1, we obtain and is given by (3.9), and the internal arrival probability for is given by (3.10). After replacing (3.9) into (3.1) and applying Loyne’s theorem [80] we have that the stable throughput region obtained from is given by (3.11). Similarly we obtain which is given by (3.12). ∎
Remark 1.
4 The Functional Equation and Preliminary Analysis
In the following, we apply the generating function approach and obtain a fundamental functional equation, which is the key element of our analysis. In order to proceed with the investigation of the functional equation we also need some crucial preliminary results that are also given below. To proceed, the evolution among queues implies
(4.1) 
we obtain the following functional equation:
(4.2) 
where
The function , , is called the kernel of the functional equation (4.2), and plays a crucial role in its solution procedure.
The definition of , implies that for fixed , with , is regular in , with , continuous in , with , and similarly with , interchanged, and the coefficients of , in the series expansion of are nonnegative, with . Then, it also follows that , are both regular in , with , continuous in , with , and the coefficients in the series expansion of are all nonnegative. Similar result holds for .
Note that some interesting relations can be directly derived by the functional equation (4.2). More precisely, setting , rearrange its terms and then taking the limit , and vice versa we come up with the following “conservation of flow relations”:
(4.3) 
Note that (4.3) has a clear probabilistic interpretation. In particular, they equate the flow of packets into the buffer of , with the flow of jobs out of the buffer of , , respectively. Just for the sake of clarity, note that the left hand side in the first of (4.3) is composed of the mean number of external arrivals per slot in the buffer of (i.e. ), while the second term corresponds to the arrivals per slot due to the presence of signals in that trigger the instantaneous transfer of a packet in . Note that such an arrival may occur only in case is not empty, i.e., with probability , which is the fraction of time user is not empty.
4.1 On the kernel
To proceed with the analysis of functional equations of type (4.2), it is required to focus on the investigation of suitable representations of the zerotuples of its kernel . For a complete discussion the interested reader is referred to [3, 81, 2]^{1}^{1}1Note that in case of Bernoulli arrivals, is a quadratic polynomial with respect to , .
Clearly, the definition of , implies that for fixed , with , is regular in , with , continuous in , with , and similarly with , interchanged. In the following we assume symmetrical characterization of the zero tuples of the kernel [77, 3], by letting
and thus,
Then, for , , ,
Since is regular in for , with fixed , , and continuous for , Rouché’s theorem implies that has exactly two roots in . According to Theorem 2.1, pp. 6566, [77], is such that , as , . Put,
In order to proceed, we must show that the contours , are both simple and smooth. To show this in the general case, we need some extra assumptions (see [3], Sections II.3.1II.3.3, pp. 151163) that will further complicate the analysis and worse the readability of the paper.
With that in mind, we will focus on the symmetrical case (see next Section). For such a case Theorem 2.1, pp. 6566 in [77] implies that:

, are both simple and smooth.

, , , , where (resp. ) denotes the interior (resp. the exterior) of , , and , .

(resp. ), then (resp. ), where denotes the complex conjugate of .

the set of the zero tuples of the kernel , and when traverses counterclockwise, , traverses clockwise.
Therefore, the solution of (4.2) is now reduced to the construction of functions , , both regular in , which satisfy
(4.4) 
for all . This problem will be transformed in the next section into a Riemann boundary value problem.
5 The Symmetrical System: Exact Analysis
As for the symmetrical system we assume , , , , , . Note that in such a case, , and . Moreover, for the stability conditions (see Theorem 3.1 or Theorem A.1) the following conditions should hold:
(5.1) 
Set,
We proceed as in [3], sections II.2.4, II.2.13, II.2.16; see also Theorem 3.2.1 in [77] or Theorem 1.1 in [81]. In particular, for the contours , , and due to the symmetry of our system the following hold:

, , ,

and functions
where,

, , ,

is regular and univalent for .

is regular and univalent for .

Then, theorem 3.2.3, p. 123 in [77] states that , , and , , are such that
(5.2) 
where , , , be the unique solution of
(5.3) 
For such a case (resp. ) represents a conformal mapping from (resp. ) onto (resp. ), while for every , with , is a zero tuple of the kernel .
Since for , is such that , (4.4) implies after some algebra that
(5.4) 
where,
Note that due to the regularity of , , and of , ,

, , is regular for , and continuous for .

, , is regular for , and continuous for .
Note also that is well defined due to the fact that , , so that , due to the regularity of , and the maximum modulus theorem [82] (similar result holds for ), and .
From the above discussion, , , can be obtained as a solution of a nonhomogeneous Riemann boundary value problem, with boundary condition given by (5.4). Its solution is solely based on the index of defined as
Note that , , while , posses a continuous derivative along and consequently, they satisfy the Ho1der condition on . Since , are simple contours and (resp. ) traverses counterclockwise (resp clockwise), we have , . Moreover,
Therefore, , and the solution of the Riemann boundary value problem is given by,
Comments
There are no comments yet.