Using Infinite-server Resource Queue with Splitting of Requests for Modeling Two-channel Data Transmission

New Radio Access Technology 3GPP New Radio has become the fundamental wireless technology in the fifth-generation networks, which allows us to achieve high data rates due to the ability to work in the millimeter-wave band. But the key feature and the main problem of 5G New Radio networks is that people themselves, cars, buildings, etc. are signal blockers, while the base stations of the fourth generation networks have widescreen broadcasting and such small obstacles do not cause loss of connection. Service providers and mobile operators are already testing the proposed technology. In this connection, the scientific community has the task of analyzing the performance of these systems and increasing it in the future. Currently, there are known studies of “basic” mathematical models of such networks. By this term, we mean models built in the simplest possible assumptions. However, due to the justified necessity of introducing new technology into the daily lives of subscribers, service providers pose the scientific community with the task of analyzing the effectiveness of the most appropriate mathematical models. For example, a technology of splitting transmitted data into two streams using as 5G and both 4G transmission technologies is considered now by 3GPP Project Coordination Group. The paper is devoted to such a problem. We consider a mathematical model of the message transmitting with the implementation of the splitting function in the communication networks of New Radio technology in the form of a resource queueing system with a renewal arrival process and non-exponential service. For this problem, an approximation of a stationary two-dimensional probability distribution of the number of occupied resources in parallel service units is obtained. It is shown that this approximation coincides with the Gaussian distribution, and its area of applicability is shown.


Introduction
The interest in the study of resource queueing systems is determined by the possibility of their application in the modeling of modern technical devices, next-generation data transmission networks, and computer systems, including cloud computing systems.The theory of stochastic processes, queueing theory, and the theory of teletraffic are served as the basis for the development of both theoretical and applied research in the field of data transmission and processing in information and communication systems.Various mathematical models and methods can be used to analyze the performance indicators of such systems.This includes, for example, queueing systems with repeated calls (retrial queues) Artalejoa and Gomez-Corral (2008); Falin and Templeton (1997), resource queueing systems Naumov et al. (2016); Samouylov et al. (2016); Tikhonenko (2010) and models of stochastic geometry Benkhelifa et al. (2020); Chu et al. (2020).
The growing popularity of wireless network services necessitates the creation of new approaches to assessing the quality of services of communication operators Buturlin et al. (2012); Galinina et al. (2014).In these networks, each active connection requires a certain amount of radio resources (call, message, video) provided to the request at the time of its receipt and should be released at the end of the connection.The amount of required resources is determined by a predefined probability distribution, which can take into account the features of various radio resource allocation schemes when analyzing the performance of wireless networks Galinina et al. (2014); Gudkova et al. (2014).Such systems can be modeled using resource queueing systems Lisovskaya et al. (2019); Naumov et al. (2016); Sopin et al. (2017); Tikhonenko (2010), where each arrival, in addition to one server, takes a certain amount of resource (deterministic or random, discrete or continuous).While the customers are being serviced, it occupies both the server and the resources, at the end of the service, the customers leave the system and frees up the occupied server and resources.Modeling of wireless communication networks using resource queues is represented by a large number of publications.However, most of the results were obtained under simplifying assumptions: deterministic resource requests, exponentially distributed service time, Poisson arrival process, simple configuration of a queueing system, which is associated with the complexity of constructing the corresponding stochastic processes (see references Basharin et al. (2009); Lisovskaya et al. (2019); Naumov and Samouylov (2017) and reviews inside them).
In addition, resource queues allow us to model any features of the distribution of resources in modern wireless networks.The inclusion of signals that trigger the redistribution of resources makes it possible to take into account the mobility of user devices Ageev et al. (2018); Samouylov et al. (2018); Sopin et al. (2017), as well as to model Network Slicing Guan et al. (2018); Song et al. (2019).Multicast technology offers a possible solution to the problem of transferring the same data to a number of devices, which leads to a significant improvement in the spectral efficiency and throughput of a wireless network.Since the same frequency band is used for several devices in the multicast mode, the data transfer speed can reach higher values in comparison to the unicast mode, where for each device only a small separate frequency band is allocated Araniti et al. (2015).Simultaneous servicing at several stations is modeled using resource queues with parallel servicing, and by using multi-stage (tandem) resource queues for sequential servicing at several stations Lisovskaya et al. (2017); Galileyskaya et al. (2019).Also, there is a possibility to take into account the heterogeneity of requests Lisovskaya et al. (2019) and their requirements for various resources Lisovskaya et al. (2019).
According to Cisco (2018) it is expected that by 2022 the number of devices which will connect to the wireless access network will reach 28 billion, this means about 3.5 devices per person as average.According to the same data, it is expected that monthly traffic on all devices will increase 2-5 times depending on the type of device: Smartphone, Tablet, Laptop or PC, Ultra High Definition TV.Moreover, the main part of the traffic will be watching video content, and a 12-fold increase will affect VR/AR traffic.It is worth noting that for such traffic, the amount of allocated radio resource does not affect the duration of the sessions, but affects the quality of the video image (4K, Immersive 360 video, VR, High Dynamic Range, High Frame Rate, 8K).
The 3GPP Project Coordination Group is considering now simultaneously supporting active links between multiple Radio Access Technologies to ensure session continuity in fifth-generation wireless networks.This is due to data transmission over the 5G network is very sensitive to interference on the line of sight, and the movement of the receiver can lead to additional problems in data transmission.One of the solutions may be a splitting of transmitted data into two streams one of which uses 5G New Radio and another uses 4G Long-Term Evolution technology.
Due to the necessity of introducing new technology into the daily life of subscribers, service providers pose the scientific community with the task of analyzing the effectiveness of the most appropriate mathematical models.A detailed review of modern research in the field of communications, modern computer networks, and information systems is reflected in Dudin et al. (2020).
This paper proposes a mathematical model of data transmitting with a splitting transmitted data into two independent streams.We represent the model in the form of a resource queueing system with the splitting of requests for their parallel servicing in two separate units, each of which contains an unlimited number of servers.Arrivals are modeled as a renewal process with a given distribution of interarrival time.The service discipline is as follows: when a request arrives, it is split into two parts and they are processed in parallel units of servers (each part in a single unit).The parts of the request capture (allocate) some amount of resources in their servicing units.This amount we model as random variables given by distribution functions that may be different for different servicing units.Service durations of the parts are stochastically independent of each other and of the amount of allocated resources, and they are determined by given distributions with finite first and second moments.After the service is completed for any part, it leaves the system immediately independently of another part.
The paper studies the two-dimensional stochastic process of the total amount of resources captured in each unit of servers.The complexity of the study of resource systems is due to the fact that at present there is no universal approach to solve the problems of this type.We use asymptotic methods Moiseev andNazarov (2014, 2016) for the study which give asymptotic expressions acceptable for practical use for the desired characteristics of the system in cases where an exact analysis is impossible.In the paper, we use an asymptotic condition of equivalent growth of service times, the practical meaning of which is that the average service speed is much less than the intensity of requests Bushkova et al. (2019).The paper continues studies started in Bushkova et al. (2019) where a similar model with Poisson arrivals was considered.
Briefly about the content of the paper.In Sect.2, we describe a mathematical model of the system under study, introduce main notations, and describe the problem under study.Section 3 contains a description of the dynamic screening method adapted to the analysis of resource queues with two units of servers that are used to solve the problem.The balance equation is formulated in this section too.In Sect.4, we apply an asymptotic analysis procedure to solve the balance equation.As a result, we obtain Gaussian approximation for the stochastic process under study (subsection 4.3).Numerical analysis and example are presented in Sect. 5. Using numerical experiments and simulation, we estimate the applicability area of the obtained approximation.Then, using technique Moiseev and Nazarov (2016), we obtain an estimation of the amount of resources that is enough that the probability of losses is equal to or less than a given value for a similar system with a limited reserve of resources.

Mathematical Model
Consider a queueing system with two service units each of which contains an unlimited number of homogeneous service devices (Fig. 1).These units correspond to different transmission environments (channels or technologies).
Arrivals are determined as a renewal process, interarrival periods have cumulative distribution function A(z).Arrived request is split into two parts: the first part goes to the first service unit and the second one goes to the second service unit.Each part occupies any free server in their unit, where they are serviced for a random times with probability distribution functions B k (x) (we denote the number of parts and respective service unit by k = 1, 2 ).Service times correspond to the part's transmission duration.
When entering into a service unit, a part of the request occupies some resources of the amount v k > 0 , k = 1, 2 distributed with the probability distribution function G k (y) , k = 1, 2 .Occupied resources correspond to channel resources captured by the transmission session (e.g., radio frequencies, bandwidths, etc.)After the service of a part of the request is completed, this part leaves the system and free resource occupied by it independently from another part of the request.
Suppose, functions A(z), B k (x) and G k (y), k = 1, 2 have finite first and second order moments, and the function A(z) satisfies the condition o(A(Δt)) = o(Δt).
In the paper, we consider the problem of analyzing the total amount of occupied resources in both channels.Problem of parts synchronization for a single request or analyzing and optimizing total transmission time are out of scope.We use an infinite-server model with unlimited resources assuming the channels have great enough bandwidth.Due to the resulting distributions will be Gaussian (see Sect. 4.3), we can use the hyper-ellipsoid method presented in Moiseev and Nazarov (2016) to estimate the number of servers and reserve amount of resources that will be enough for the systems with finite limitations.
Let us denote the total amount of resources occupied in the k -th service unit at time instant t by V k (t) .The goal of the study is to obtain the steady-state probability distribution of Fig. 1 Resource queueing system with splitting of requests two-dimensional stochastic process {V 1 (t), V 2 (t)} .This process is not Markovian and classic methods of study do not allow to solve the problem.Due to this, we use the dynamic screening method Moiseev and Nazarov (2016) adapted for resource systems Lisovskaya et al. (2016) combined with the splitting of requests.

Dynamic Screening Method and Balance Equation
The dynamic screening method allows us to transform the original problem of studying of a queueing system by the analysis of some stochastic (screened) process which is constructed in the following way.Let us draw three parallel time axes.We number them from 0 to 2 (Fig. 2).Axes with numbers 1 and 2 will correspond to the service units with the same numbers, and the axis with number 0 is used to display arrival epochs.
Let some moment of time T in the future be fixed.We denote the probability that a part of a request arrived in the system at the time instant t ≤ T will be served at the time instant T by S k (t) , where k = 1, 2 is the number of part of the request (or the number of service unit where it is served).These probabilities we call as screening probabilities for the 1st or for the 2nd axis.We will mark epochs on these axes if original arrival was screened on them with probabilities S 1 (t) and S 2 (t) respectively.A sequence of these epochs on each axis we call as screened flow or process.Obviously, that original arrival do not generate an epoch in k-th screened flow with probability 1 − S k (t) , k = 1, 2 .This case corresponds to the situation when a part leaves the k-th unit of servers before the time instant T. Generally, the S k (t) functions also depend on the fixed moment T, i.e. S k (T, t) , but in the future we will omit T to simplify the notation.
Let us denote the number of requests serviced in k-th unit at time instant t by i k (t) , and the total number of epochs screened on k-th axis before time instant t by n k (t) .Suppose that at some initial instant in time t 0 < T the system is empty.Then, according to the method of constructing the screened processes, the following equality holds Moiseev and Nazarov (2016): for k = 1, 2 and m = 0, 1, … This is the main idea of the dynamic screening method: we analyze two-dimensional process {n 1 (t), n 2 (t)} and, substituting t = T , we obtain a result for the process {i 1 (t), i 2 (t)} at time instant T.
For resource queueing systems, we denote the amount of resources allocated by the original request arrived at time instant t and screened on the k-th axis by W k (t) , k = 1, 2 .Then the probability distribution of the amount of resources allocated in the system at time instant T coincides with the probability distribution in the screened process Lisovskaya et al. (2016): So, we can analyze process {W 1 (t), W 2 (t)} instead of process {V 1 (t), V 2 (t)} and obtain a result for the last one by substituting t = T into the result for process {W 1 (t), W 2 (t)}.
Process {W 1 (t), W 2 (t)} is not Markovian, therefore we consider a three-dimensional Markov process {z(t), W 1 (t), W 2 (t)} , where z(t) is the length of the interval from time instant t until the next request arrives in the system (residual time for the arrival process).We denote its probability distribution as follows: The following events may happen in interval (t, t + Δt]: -no new customers enter into the system with probability P(z + Δt, w 1 , w 2 , t) − P(Δt, w 1 , w 2 , t); -exactly one customer enters into the system and will be in both service units in moment T, and the length of the interarrival period until the next customer's arrival is less than -exactly one customer enters into the system and will be in only first service unit in moment T, and the length of the interarrival period until the next customer's arrival is -exactly one customer enters into the system and will be in only second service unit in moment T, and the length of the interarrival period until the next customer's arrival is less than z -with probability A(z)(1 − S 1 (t))S 2 (t) -exactly one customer enters into the system and will finish service on both service units before moment T, and the length of the interarrival period until the next customer's arrival is less than z -with probability A(z)P(Δt, w 1 , w 2 , t)(1 − S 1 (t))(1 − S 2 (t)); -more than one customer enters into the system with probability o(Δt).
For the three-dimensional Markov process, we apply the formula of total probability and write down the equality (1) (2) whose solution is defined as follows: Let us introduce partial characteristic functions where j = √ −1.

Asymptotic Analysis
Since it is not possible to find an explicit form of a solution to the problem ( 5)-( 6), we will seek the solution under the asymptotic condition of equivalent growing of service times in the units of servers.This asymptotic condition means proportional growth of the average service times in both service units and it is taken from practice.For example, during transmitting high-quality streaming video, resources are used to provide quality of the content, and service time corresponds to the duration of transmitting process.So, in such a case, average service times in both transmission channels will be proportional each other.

Asymptotic Analysis of the First Order
we can write the asymptotic condition of equivalent growing of service times in the units Theorem 1 The first-order asymptotic characteristic function Proof Let us perform the following substitutions in Eqs. ( 5)-( 6): we derive the equation with initial condition Performing the limit transition → 0 in Eq. ( 9) and using the expansions and we obtain Suppose function F 1 (z, x 1 , x 2 , , ) is differentiable in = 0 .Then, taking into account (4), we will find it in the expansion form (the Maclaurin series on powers of ) where Φ 1 (x 1 , x 2 , ) is some function which satisfies initial condition Then let us substitute (11) into (9): where denoting r � (0) = .Let us, make in (12) limit transition z → ∞: Dividing the obtained equation by and tend → 0 , we derive the following differential equation for function Φ 1 (x 1 , x 2 , ): which solution is as follows: Substituting ( 13) into ( 11) and performing substitutions inverse to (8), we obtain The theorem is proved.
Remark Function F 1 (z, x 1 , x 2 , , ) must satisfy the condition of existing finite derivative in = 0 to we be able to use expansion (11).But the problem of checking such property or trying to establish its meaningful for specific distribution/characteristic functions seem very hard to be solved and lay out of the topic of the paper.Due to this, we rely on the fact that we obtain an expression for the characteristic function at the end of the derivations, and on the comparison of this analytical result with a simulation one presented in Sect. 5.The same applies to the differentiability of function F 2 (z, x 1 , x 2 , , ) in the next section.
Letting z → ∞ in ( 14), we obtain the complete characteristic function of the stochastic pro- cess {V 1 (t), V 2 (t)} , and setting t = T and letting t 0 → −∞ , we obtain the characteristic function of the stationary probability distribution of the two-dimensional process under study: The obtained approximation determines only the means of resource amounts.To construct a qualitatively more accurate approximation, we perform a second-order asymptotic analysis.

Asymptotic Analysis of the Second Order
Let us prove the following statement.

Theorem 2
The second-order asymptotic characteristic function h (2) (z, u 1 , u 2 , t) ≈ h(z, u 1 , u 2 , t) of the probability distribution of the process {z(t), W 1 (t), W 2 (t)} under the condition of equivalently growing service times has the form , a and 2 are the mean and variance of random variable given by distribution function A(x).
Proof Taking into account expression (7) for function h (1) (z, u 1 , u 2 , t) , we make the follow- ing substitution in (3): We obtain Let us make the following substitutions in this equation: We obtain equation ( 15)

with initial condition
Let us make here asymptotic transition → 0 , denoting lim →0 F 2 (z, x 1 , x 2 , , ) = F 2 (z, x 1 , x 2 , ) , we derive Suppose function F 2 (z, x 1 , x 2 , , ) has finite derivatives of the second and greater orders in = 0 .Then we can represent this function in the expansion form where Φ 2 (x 1 , x 2 , ) and f 2 (z) are some functions.The first one satisfies the initial condition Expansion ( 20) is the Maclaurin series on powers of where the second term is chosen in the specific form to avoid cumbersome derivations that the authors have made before their presentation in the paper, including the property f 2 (∞) = 0.
Substituting (20) into (18), dividing on and setting → 0 , we derive the following equation for function f 2 (z): Let us substitute expansion (20) into (18) and make asymptotic transition z → ∞ , we derive the following equation for function Φ 2 (x 1 , x 2 , ): Its solution is as follows: Substituting this expression into (20) and making substitutions inverse to (17), we obtain:

Main result
Let us set in (15) z → ∞ , t = T to satisfy (1), and t 0 → −∞ to reach steady-state regime, then we obtain the following expression for characteristic function of joint distribution of resources allocated in both units of servers in steady-state regime: where So, we can see that the probability distribution of the two-dimensional process {V 1 (t), V 1 (t)} is asymptotically Gaussian with mean vector

and covariance matrix
We may notice that in the case of the Poisson arrival process we obtain exact matching with a result obtained in Bushkova et al. (2019).

Numerical Analysis and Example
Since we have no exact results to compare obtained approximation with, we use simulation and perform analysis of the approximation accuracy and applicability area in a numerical form.Consider the queueing system described in Sect. 2 with the following parameters.Interarrival time has a gamma distribution with shape and rate parameters both equal to 0.5.So, the intensity of the arrival process 1 is equal to 1, and we obtain 2 = 1 too.
Let us represent asymptotic parameter in the form = 1∕N , where N be some scale parameter.This is made just for convenience of presentation.
Service times have gamma distributions with shape and rate parameters equal to 0.5 and 0.5 for the first unit of servers, and 1.5 and 0.75 for the second one.So, the average ser- vice time b 1 is equal to N for the first unit of servers and b 2 = 2N for the second one.
We choose gamma distributions for the probability distribution of resources allocated by a single request.Their shape and rate parameters are equal both to 2.5 for the first service unit, and 0.4, 0.8 for the second one.So, the average amount of resources allocated by a request in the first unit of servers ( a 1 ) is equal to 1 and it is equal to 0.5 for the second unit ( a 2 ).
So, according to (23), we obtain the average amount of resources allocated in each unit of servers in a steady-state regime equal to N.
Varying scale (asymptotic) parameter N, we can obtain an applicability area for obtained results.To estimate the accuracy of results, we perform simulations and calculate error in the form of Kolmogorov distance where F(y 1 , y 2 ) is empiric cumulative distribution function constructed on the base of sim- ulations results and G(y 1 , y 2 ) is two-dimensional Gaussian cumulative distribution function with means vector (23) and covariance matrix (24).
Obtained results are presented in the first row of Table 1.The same analysis also is done for various types of probability distributions for allocated resources with parameters that give a 1 = 1 and a 2 = 0.5 as for the case of a gamma distribution.In Table 1, you can find the results for Uniform, Weibull, Log-normal, and Pareto distributions.As we can notice, for all cases, while N is growing, Kolmogorov distance d is decreasing, therefore, results ( 22)-( 24) become more precise.For example, if we choose d ≤ 0.05 as enough error for the results accuracy, then we can conclude that results ( 22)-( 24) are applicable for N ≥ 40 ( ≤ 0.025 ).Figures 3 and 4 show the probability density functions of the Gaussian approximation ( 22)-( 24) in comparison with distribution built on the base of simulations result for the first and the second service units and various values of scale parameter N. As we can see, while N is increasing, the approximation and actual distribution get closer to each other.
In practice, we usually have no infinite reserve of resources to be used, but we can use the results from Sect.4.3 in the case of the system with limited resources.Consider the problem of design such transmission system and the task of choosing maximum values (limits) of resources provided in each channel.On one hand, these limit values should be minimal as possible to make a cost of the system very low.On the other hand, these values should be enough to avoid losing of requests due to the absence of resources for their servicing.We call values of the total amount in resources that satisfy both criteria as optimal and denote them by a opt k where k is the number of the channel.Evaluation of these values is based on the Gaussian distribution ( 22)-( 24) and technique explained in Moiseev and Nazarov (2016).Consider the probability that a request arrived in the system will be rejected due to there is no enough amount of resources for its transmission in any of the channels because the most resources are occupied by other requests.Let us estimate the total amount of resources in each channel that will be enough to make this probability be not more than some given value p loss .According to Moiseev and Nazarov (2016), due to the joint probability distribution of occupied resources is two-dimensional Gaussian, we can evaluate the optimal values of the total resources in each channel using the following expressions: Here r is the so-called hyper-ellipsoid radius that depends on p loss and can be evaluated using the method described in Moiseev and Nazarov (2016).
For example, consider the system described above with gamma distribution of resource allocation.We choose N = 40 to be sure that approximations ( 22)-( 24) are accurate ( 25) enough.Then we obtain the following mean vector and covariance matrix for occupied resources in both channels: Applying formula (25), if we choose value p loss = 10 −3 then r ≈ 3.717 and we calculate (rounding up): a opt 1 = 72 , a opt 2 = 74 .These values are the total amounts of resources in the corresponding channels that will be enough to a new request arrived in the system in any time with probability 1 − p loss finds enough resources in both channels to be transmitted in parts via the channels.We may interpret these amounts as the number of radio transmission bands used in the channels of the system.It is interesting that if we choose p loss = 10 −6 decreasing losing probability by a thousand, we obtain r ≈ 5.256 and a

Conclusions
According to Cisco, it is expected that by 2022 the number of devices connected to wireless networks will reach about 3.5 devices per person.Video content will be the main part of the traffic.The amount of allocated radio resources for such traffic does not affect the duration of the sessions but affects the quality of the video image.The 3GPP Project Coordination Group is considering now supporting active links between multiple Radio Access Technologies simultaneously to ensure session continuity in fifth-generation wireless a opt p loss 10 -1 10 -2 10 -3 10 -4 10 -12 10 -11 10 -10 10 -9 10 -8 10 -7 10 -6 10 -5 a 1 opt a 2 opt Fig. 5 Optimal amounts of resources a opt 1 and a opt 2 for channels of a finite-resource system when various values of losing probability p loss is choosing networks.One of the solutions may be a splitting of transmitted data into two streams one of which uses 5G New Radio and another uses 4G Long-Term Evolution technology.
In the paper, we have considered the problem of analyzing the total amount of resources allocated in both channels.To do this, we have used an infinite-server resource queueing model with the splitting of incoming requests into two service units (transmitting channels) that have unlimited resources assuming the channels have great enough bandwidth.We have found the solution under the asymptotic condition of equivalent growth of service times in the units which has a practical interpretation: during transmitting high-quality streaming traffic, resources are used to provide quality of the content, and service time corresponds to the duration of transmitting process.We have found the solution in the form of the distribution approximation and then we have analyzed the applicability area of it.Due to the resulting distribution (approximation) is two-dimensional Gaussian, we have shown for example how to estimate the reserve amount of resources for the system with finite limitations of resources.
The approach presented in the paper may be applied to study resource queues with the splitting of requests with the arrival processes of other various types (Markovian arrival process, semi-Markov process, a mixture of arrival processes, etc.).And also to obtain the main characteristics of splitting of transmitted data into three and more streams (for example, 5G New Radio, 4G Long-Term Evolution and Wi-Fi technology).But we note that the formulas' dimensional will be extra-large.

Fig. 2
Fig. 2 Constructing of the screened processes k (y k ) is the average amount of resources allocated by one request in the k-th service unit, k = 1, 2.

Fig. 4
Fig. 4 Probability density functions of Gaussian approximation (dashed line) and distribution of simulations result (solid line) for the second service unit for various values of scale parameter N which means that we need enlarge amount of available radio resources only by 20% approximately.Figure5depicts changing of a opt 1 and a opt 2 while p loss decreases.

Table 1
Kolmogorov distance for various values of scale parameter N and various distributions of allocated resources Fig. 3 Probability density functions of Gaussian approximation (dashed line) and distribution of simulations result (solid line) for the first service unit for various values of scale parameter N