Construction and Analysis of Queuing and Reliability Models Using Random Graphs

: In this paper, the use of the construction of random processes on graphs allows us to expand the models of the theory of queuing and reliability by constructing. These problems are important because the emphasis on the legal component largely determines functioning of these models. The considered models are reliability and queuing. Reliability models arranged according to the modular principle and reliability networks in the form of planar graphs. The queuing models considered here are queuing networks with multi server nodes and failures, changing the parameters of the queuing system in a random environment with absorbing states, and the process of growth of a random network. This is determined by the possibility of using, as traditional probability methods, mathematical logic theorems, geometric images of a queuing network, dual graphs to planar graphs, and a solution to the Dirichlet problem.


Introduction
Random processes on graphs play an important role in applied probability theory (see for example [1][2][3][4][5] and references therein). The study of these processes, being closely related to various applications, includes consideration of a whole range of different issues: combinatorial probability, discrete random sets, applied discrete mathematics, random lattices, etc. All these branches of mathematics and, in particular, probability theory are gradually being transformed into instruments of computer science.
The use of the construction of random processes on graphs allows us to expand the models of the theory of queuing and reliability by constructing reliability models arranged according to the modular principle; queuing networks, in the nodes of which there are multi-channel systems with failures; reliability networks in the form of planar graphs; the random process of changing the parameters of the queuing system in the presence of absorbing states; combinations of special functions for studying the process of growth of a random network. This is determined by the possibility of using, along with traditional methods of probability theory, mathematical logic theorems, geometric image of functioning of a queuing network, the use of dual graphs to study reliability networks in the form of planar graphs, the use of methods for solving the Dirichlet problem to calculate the limit distribution of a random process at absorption points, and the use of special functions to study the convergence rate to the limit value of the average degree of vertices of a growing random network.
The theory of random graphs began in 1950s from several papers by Erdos and Renyi. In [6], along with purely geometric and topological properties of graphs, the dynamics of processes on graphs was studied. Apparently, the most striking continuation of them are works on modelling of random growing networks [7]. As applications of random processes on graphs, along with lattice physical models, we can now specify the questions of traffic prediction [8], machine learning [9], contact processes on graphs [10], computational biology [11] and etc.
In this paper, a stochastic model of reliability system arranged by the modular principle is constructed. Asymptotic formula for probability of incoherence in planar graphs with high reliable edges is obtained. Product theorem for queuing networks with multi server nodes and failures is proved. Limit distribution of customers numbers in one server queuing system M|M|1|∞ in random environment, depending on initial environment state, is calculated. The convergence of the pre-limit distributions to the limit distributions in Barabasi-Albert randomly growing network is proved.
It is shown that random processes on graphs may be used to construct and to analyse a model of a random network arranged by the modular principle. The choice of this task is due to the fact that recently in computer science and informatics there have been many papers devoted to systems built on the modular principle. Therefore, it became necessary to build a fairly general model of such systems within the framework of logical probabilistic modelling. In this paper, this model is constructed with elements of mathematical logic and graph theory. Using the well-known theorems of mathematical logic (see for example [12]), it is shown that this model reduces to calculating the reliability of a random two-pole belonging to the class of recursively definable parallel-sequential two-poles. As a result, the algorithm for calculating it becomes linear in terms of the number of edges of the two-pole.
Another application of the suggested approach is the construction of asymptotic formulas for calculating the probability of incoherence of a planar graph with highly reliable edges. It is shown that the solution of this problem is based on the transition from the original planar graph to a dual graph, the vertices of which are the faces, and the edges are the edges separating these faces (see for example [13]). This circumstance makes it possible to significantly reduce the calculation of the asymptotics of the probability of incoherence, since it is sufficient to calculate the degrees (from two to five) of the adjacency matrices of a dual graph. Thus, a model of a planar graph with highly reliable edges is chosen and the asymptotics of the probability of its incoherence is constructed through the minimum length of cuts and the number of cuts of the minimum length. In turn, the determination of the length and number of such cuts is made using the construction of a dual graph and calculating the number of cycles of minimum length in it. This allows us to proceed to the construction of graphs laid on two-dimensional manifolds without crossing edges at their internal points.
Several problems of queuing theory are considered, the reduction of which to random processes on graphs makes it possible to significantly simplify their consideration. The first of these tasks is the study of an open queuing network of the Jackson type, in which the network nodes are multi-channel queuing systems with failures. These networks are widely used in communication theory and in data transmission networks [14]. For such a study, it is possible to construct protocols for the functioning of these networks using graphs of transient intensities and obtain an analogue of the product theorem. In Jackson's theorem [15], the stationary distribution of the numbers of customers in network nodes is decomposed into the product of the same distributions but in isolated nodes. In accordance with this theorem, the graph of transient intensities is decomposed into a union of complete (consisting of a number of vertices equal to the number of network nodes) sub graphs that can have common vertices, but do not have common edges. Thus, this network is built using a protocol, that allows or prohibits the arrival of input customers to a particular node. The protocol itself is described by a weighted digraph of transient intensities, consisting of sub-graphs (in the case of two nodes in the network-of a triangular form) that do not have common edges. By analogy with the Jackson theorem, a distribution is constructed, that satisfies the stationary Kolmogorov-Chapman equations within individual sub-graphs. This choice of protocol makes the analysed queuing network manageable, i.e., it allows convenient calculations of the stationary distribution. The problem of calculating the limit distribution in a queuing system located in a random environment, in which the system parameters are described by a certain Markov process, is considered. Usually, such limit distributions are searched under the conditions of ergodicity of the Markov process, describing this system. However, if the randomly changing parameters of the system obey the Markov process with absorption in some states, the ergodicity condition is violated and the limit distribution becomes dependent on the initial state of the process, describing the dynamics of the parameters of the queuing system. This dependence is determined by the solution of the discrete analogue of the Dirichlet problem [16] in calculation of a probability of ruin the player (see for example [17]), which allows us to calculate the probability of absorption depending on the initial state of the medium parameters. Thus, it is possible to construct a model of a queuing system, the marginal distribution of which depends on the initial state of the environment and, therefore, this model is not ergodic. That allows us to set and to solve the tasks of statistical assessment of the state of the environment based on observations of the number of customers in the queuing system. The considered queuing model is the simplest, it admits numerous generalizations and is calculated by solving the Dirichlet problem, defined by the graph of transient intensities of the state of the medium with absorbing states.
The last task in the work is to study the Barabasi-Albert and Dorogovtsev's models [7,18] of a randomly growing network. One of the most convenient methods for calculating the limit distribution of the degree of a network node is the continuous approximation [18,19]. This approximation is based on the hypothesis about the asymptotic behavior of the specified distribution, when the number of model steps tends to infinity. In [20], the problem of catastrophic cascade of failures in interdependent networks is considered. So, it is necessary to consider and to analyse the problem of a strict formal description of a statistical ensemble of random networks with a given the distribution of nodes by the number of connections. The absence of a rigorous mathematical proof of the convergence of the pre-limit distributions to the limit distributions, calculated using the continuum approximation, makes the results obtained vulnerable. In an accordance with the problem put in [20] to justify the continuum approximation, used in calculating the average degree of network vertices, an exact asymptotic formula is constructed for the rate of convergence of the average degree of network vertices to a limit value at a time tending to infinity. These models give sufficiently simple description of the Internet network [21], using a power distribution of the degree of network nodes. The solution of this problem is based on combinatorial formulas expressed in terms of gamma functions.
The substantiation of the continuum approach for calculating limit distributions in the main models of growing networks (exponential model, the Barabasi-Albert model and the Dorogovtsev model) is represented on an example of the Barabasi-Albert model. To substantiate the continuum approach, the exact asymptotic's of the differences between the limit and pre-limit distributions are constructed for degrees of nodes when the number of steps of the model tends to infinity. The results obtained are based on recurrent relations and on asymptotic expansions of the differences between the limit and pre-limit distributions.

Logical-Probabilistic Models of Complex Systems Constructed on the Modular Principle and Their Reliability
I.A. Riabinin introduced a concept of logical-probabilistic modelling as an extension of inductive logic and substitution rules of boolean arguments in functions of the logic algebra by their truth probabilities (see for example [22]). However, recently, the concept of systems built on the modular principle has appeared in computer science (see, for example [23]) with manifold applications to machine tools and robots, industrial processes and etc. However, in computer science there is not mathematical investigation of this concept. For this purpose, the following mathematical model is constructed in this section. Suppose that the system consists of n elements, the state of which is determined by boolean variables x k , k = 1, . . . , n. If x k = 1, then element k is in a working state, otherwise it fails. The state of an entire system is determined by a boolean function A(x 1 , . . . , x n ), where independent random variables x k satisfies the condition P(x k = 1) = p k , k = 1, . . . , n. Then, the reliability of the system is determined by the expression P(A(x 1 , . . . , x n ) = 1).
Boolean function A satisfies the natural monotonicity condition consequently [12] it belongs to the following recursively defined class of boolean functions A : Let us now define a class of boolean functions, that characterize systems according to the modular principle, by a recursively definable class of boolean functions A * ⊂ A : So, boolean function A(x 1 , . . . , x n ) ∈ A * may be depicted as in Figure 1. From the definition of the class A * , the equalities follow Using these formulae, it is possible to calculate P(A(x 1 , . . . , x n ) = 1) in n steps as reliabilities of parallel and/or serial connections of two-poles.

Disconnection Probability of Planar Graph with High Reliable Edges
The next model of reliability theory based on graph-theoretic representation consists of obtaining an asymptotic formula for the probability of incoherence of a planar graph with highly reliable edges. Assume that S is an undirected connected planar graph. It means that the graph S is placed on the plane and its edges do not intersect in internal points. Suppose that the graph S edges are highly reliable and fail independently with the probability h → 0. A cut of the graph S consists of edges whose deletion causes the graph to become disconnected. Further, we consider only cuts with minimal length (number of edges).
Calculate an asymptotic of a disconnection probability P h for such graph. From the Burtin-Pittel formula (see, for example [24]) we have: where m is (minimal) cut length and c m is the number of cuts with length m. In [13], it is proved that there is one-to-one correspondence between all cuts (with minimal length) in planar graph and all cycles with minimal length in dual graph S * .
Let us denote m * the minimum vertices degree in S, then m ≤ m * and from Euler's formula [13] we have m * ≤ 5, where an equality may be reached.
Since we are interested in cuts with minimal length and, therefore, dual cycles with minimal length, then we shall limit ourselves to enumerating simple (without selfintersections) cycles of length 2, 3, 4, 5 in the dual graph. Denote a ij a number of edges contained in the intersection of the faces z i ∩ z j , i = j, a ii = 0, in the planar graph S.
Elements of l-th degree A l , l > 1, of the matrix A are denoted by a (l) ij . Using wellknown formulas of the graph theory [25] we can calculate the number c k of simple cycles of length k, k = 2, 3, 4, 5 in the dual graph: Then, the following equalities are valid: m = min(k : 2 ≤ k ≤ 5 : c k > 0). It follows from these formulas that the algorithm defined by them for calculating the asymptotic constants m, c m in the case of a planar graph has a complexity of no more than cubic.
In the examples below, it is easy to see that this complexity can be significantly reduced. It should be noted that it is possible not only to calculate the constants m, c m , but also to list the cycles of the minimum length in the graph S * and the corresponding sections interpreted by them as the bottlenecks of the graph S using computer algebra algorithms (see, for example [26]). It is worthwhile to note that the idea of using dual variables is also widely used in deterministic optimization problems (see for example [27,28]). These results may be generalized for graphs laid on two-dimensional surfaces.
Numerical example. Consider the graph in Figure 2. In this graph m = 3, c 3 = 22. The disconnection probability in this graph is calculated by the Monte-Carlo method with 10 6 realizations (designated P MC ) and by asymptotic . Obtained numerical results are represented in Table 1.

Construction and Calculation of Queuing Networks with Failures
However, graph-theoretic consideration is convenient to use not only in models of reliability theory. It turns out that it is possible to represent the classical Jackson network with various modifications caused by modern problems of the theory of communication and data transmission networks using original graph-theoretic considerations. This section discusses a queuing network with failures in nodes. Such a network arises in the mathematical theory of communications [14]. Assume that G is an open queuing network with m nodes (one-server or multi-server queuing systems), input intensity λ > 0, and route matrix Here, θ ki is a probability of customer transition from ith node to jth node after its serving in ith node (0th node is an external source or external space). Then, the route matrix Θ satisfies the conditions and so the system of motion equations has a unique solution with λ 1 > 0, . . . , λ m > 0.
Suppose that kth node has N k servers with serving intensities µ k , k = 1, . . . , m. Denote n k (t) a number of customers in kth node at moment t. Then (n 1 (t), . . . , n m (t)), t ≥ 0, is discrete Markov process, describing the network G.
For simplicity and without limiting generality everywhere, we further assume m = 2. It is convenient to divide the construction of a weighted digraph of transition intensities on the set of states of a discrete Markov process (n 1 (t), n 2 (t)) into two stages. First, an undirected graph Γ is constructed (Figure 3), consisting of triangular sub-graphs with vertices (n 1 , n 2 ), (n 1 + 1, n 2 ), (n 1 , n 2 + 1). Then, each edge of such a triangular sub-graph is replaced by two oriented edges directed in opposite sides with corresponding transition intensities A, B, C, D, E, F. The combination of these edges forms a triangular digraph Ψ(n 1 , n 2 ) (Figure 4). Next, each undirected triangular sub-graph in Γ should be replaced by the corresponding digraph of the form Ψ(n 1 , n 2 ). Here, the coefficients A, C characterize the intensities of the arrivals of input customers to the network G nodes, the coefficients B, D, E, F -the intensities of service at the nodes of the digraph Ψ(n 1 , n 2 ), A = λθ 02 , B = µ 2 n 2 θ 20 , C = λθ 01 , D = µ 1 n 1 θ 10 , E = µ 1 n 1 θ 12 , F = µ 2 n 2 θ 21 . ( Thus, digraph Ψ(n 1 , n 2 ) defines the network operation protocol, which includes permission or prohibition for certain transitions of customers. If edge [(n 1 , n 2 ), (n 1 + 1, n 2 )] or edge [(n 1 , n 2 ), (n 1 , n 2 + 1)] belong to the digraph Ψ(n 1 , n 2 ) then there is an allowance for input customer arrival to the first node or to the second node. If edge [(n 1 , n 2 ), (n 1 − 1, n 2 )] or edge [(n 1 , n 2 ), (n 1 , n 2 − 1)] belong to the digraph Ψ(n 1 , n 2 ) then there is an allowance for customer exit after service at the first node or at the second node from the network G. If these edges do not belong to the digraph Ψ(n 1 , n 2 ) then appropriate transitions are prohibited.
Thus, the Markov process (n 1 (t), n 2 (t)) satisfies the stationary Kolmogorov-Chapman equations on the Ψ(n 1 , n 2 ) digraph. However, such different digraphs have not common edges, which means that the distribution P(n 1 , n 2 ) in this form satisfies the stationary Kolmogorov-Chapman equations on the entire graph Γ.
Another generalization is the change in the puncture of the network functioning in accordance with the change of the graph Γ as follows ( Figure 5) with N ≤ min(N 1 , N 2 ). We shall not discuss closed queuing networks in detail and only note that for them a connected graph Γ is constructed from complete digraphs Ψ with vertices (n 1 + 1, n 2 , . . . , n m ), (n 1 , n 2 + 1, n 3 , . . . , n m ), . . . , (n 1 , n 2 , . . . , n m−1 , n m + 1).

Limit Distribution in Queuing System M|M|1|∞ with Randomly Varying Intensity of Input Flow
The graph-theoretic approach used in the previous section can be extended to the calculation of stationary distributions in queuing systems located in a random environment. Moreover, this approach allows us to identify new problems in which the stationary distribution does not obey the ergodicity hypothesis and to calculate the dependence of the stationary distribution on the initial state of the external environment. A large number of papers are devoted to queuing models with randomly varying parameters (see, for example [30,31]) with limit distributions, which do not depend on the initial state. However, there are models in which a random process, determining the parameters, develops on a segment with absorbing ends.
In the queuing system M|M|1|∞ with λ(t) ≡ λ limit distribution of customers numbers Assume that Λ(0), . . . , Λ(N) < µ and denote ρ(0) = Λ(0) It is not difficult to obtain that limit distribution P z 0 (k), k = 0, 1, . . . , of customers numbers k(t) in the system M satisfies the equality, depending on initial state z 0 So, the Markov process (λ(t), k(t)) is not ergodic. These results may be spread onto multi-server systems with and without failures, on systems in which input to service intensities are randomly varying.

Convergence to the Limit Distribution in the Barabasi-Albert Model of a Growing Random Network
This section considers a random growing network. It cannot be attributed either to models of reliability theory or to queuing models. It is of an intermediate nature, because it is used in the description of the Internet and, therefore, is close to the theory of queuing. On the other hand, purely mathematically, this network is close to the problems of reliability theory and has analogues with systems built on the modular principle.
The main purpose of this section is to estimate the rate of convergence to limit distributions in the Barabashi-Albert model of growing random network (Figure 7) [7]. Let us consider the growing random network [7], in which a newly appearing vertex connects to any of the already existing vertices with a probability proportional to the degree of this vertex. At the initial step t = 1, there is a single vertex and a loop associated with it. At the step t ≥ 1, we can determine the degree k of the vertex s, equal to the number of incident edges.
Denote p(k, s, t) the probability that at step t ≥ 1, k edges of an undirected network graph are connected to the vertex s. Formally, we can define p(k, s, t) = 0, s > t, and so the following relations are obtained in [18]: denote P(k, t) = 1 t t ∑ s=1 p(k, s, t).
Denote ∏ t t+1 = 1, from Formula (6) it is not difficult to obtain for t ≥ 1 From here, it is not difficult to obtain

Remark 1.
In the future, we shall need to understand the expression , if A, t are natural numbers. Here Γ(z) is a gamma function.
Proof. First, let A > 0 not be a natural number. Denote hence, the Formula (9) if t = 1 is proved. Suppose that this formula is true for a given t > 1, we prove it for t + 1 : So, the relation (9) is true for all natural t. The relation (10) may be proved similarly. For natural A, the function V(t) in Formulas (9) and (10) is understood in the sense that remark 2 contains.
Proof. For A > 0 it is possible to obtain from Stirling's formula: Indeed, it is well known ( [32], Chapter 1) that Let k = 1, then the relation (12) follows from Formulas (7), (9) and (13) when A = 1/2 : Let us assume that the relation (12) holds for k − 1 > 0, and prove its validity for k. We look for f k (t) in the form f k (t) = a k (t) Therefore, by virtue of Formulas (10) and (11) for A = k/2 a k (t) = 4 k(k + 1) In turn, by the assumption of induction f k−1 (t) ∼ t −3/2 3 √ π , t → ∞, properties of the gamma function and Formula (13) Hence the validity of the asymptotic relation (12) follows for an arbitrary k.
Let us consider a model of Dorogovtsev's growing network, in which an edge from a newly appearing vertex to existing vertices has a probability proportional to the degree of this vertex plus a. Here, a > 0 is a parameter of the model. In [18], the following relations are obtained: From Formula (6), we have By analogy with (7) and (8) it is not difficult to obtain Denote A k = (k + a)/(a + 1), k ≥ 0, Theorem 2. The following relations are fulfilled Proof. By virtue of Formulas (9), (13) and (14) for A = A 0 we have For k > 0, we are looking for f k (t) = P(k, t Then, by virtue of Formula (13) for t → ∞ and A = A k we have In turn, from the assumption of induction f k−1 (t) ∼ C k−1 t −1−A 0 , by virtue of the relation A k > A 0 and Formula (13) Thus, the asymptotic relation (16) is proved for an arbitrary k.

Remark 2.
Dorogovtsev's model [18] is widely used in models of growing random networks. For small a, it gives sufficiently simple description of the Internet network, using a power distribution of the degree of network nodes in which its parameter is close to two [21].

Remark 3.
Similarly, an exponentially growing random network, in which the newly appearing vertex connects to any of the already existing vertices with equal probabilities, may also be studied.

Discussion
This paper presents the main examples known to the author of the use of random processes on graphs for modelling and analysing various modern computer science systems. At the same time, each of the above examples allows for different generalizations. However, these generalizations require a more detailed definition and research using both probabilistic methods and methods of graph theory.
In particular, models of systems built on the modular principle have just begun to penetrate into logical-probabilistic modelling [22]. Their application can expand the possibilities of this method of stochastic systems modelling. Models of reliability systems represented by planar graphs with highly reliable edges are built at the junction of reliability theory and planar graph theory. Such a study allows us to expand the possibilities of logical probabilistic modelling by introducing ideas from graph theory and topology [13] into it. The use of oriented graphs defining protocols in queuing networks makes it possible to obtain formulas for limit distributions in data transmission networks, in particular, describing networks of the latest generations [14]. The calculation of limit distributions in queuing systems placed in a random environment and not obeying the ergodicity hypothesis allows us to expand the capabilities of known methods [30,31] for performing such calculations. The asymptotic formulas of the difference between the prelimit and limit distributions in models of growing random networks are accurate. The asymptotic formulas of the difference between the limit and limit distributions in the models of growing random networks obtained in this paper are accurate. This allows us to strengthen the results of the papers [7,18].

Conclusions
The possibility of strengthening the known results on the theory of reliability, queuing and on the theory of growing random networks is presented. It seems to the author that the continuation of this research can become a source of new and significant results in the listed and already known areas of probabilistic modelling.