1. Introduction
In this research, we consider a single-server FIFO (First In First Out) queueing system , with a superposition of n independent renewal stationary inputs (components) corresponding to customers of different classes. The superposition may include component processes with general (component-dependent) interarrival times. The service time distributions are also assumed to be general and class-dependent. (In what follows, we use the terms ‘distribution’ and ‘distribution function’ as synonymous terms.)
The superposed input process system has been intensively studied since the 1970s (see, for example, [
1,
2,
3,
4,
5,
6,
7], among others). It is well known that, unless the components are Poissonian, the superposed input process system is not (classically) regenerative. On the other hand, the regeneration property is quite important for both the stability and performance analysis of the queueing system. Indeed, it admits accurate estimation of stationary performance metrics on the basis of regenerative simulation [
8,
9,
10,
11], which allows one to use the independence of regeneration cycles (of a basic process) in the framework of the regenerative variant of the Central Limit Theorem (CLT) in order to estimate the required metrics with a given precision. (In
Section 5, we discuss this topic in more detail.) What is especially important is that the values of a queueing process, within a regeneration cycle, are typically highly correlated, and a direct application of the classical variant of the CLT is not justified.
In general, the superposed process generated by independent renewal processes can be transformed, under mild assumptions, into a regenerative process with the so-called
one-dependent regeneration cycles (see [
6,
12]). This construction, which is based on splitting and coupling procedures, is rather complicated and hardly applied in practice to evaluate the required performance measures by simulation. (We will discuss it in
Section 3.) Moreover, to apply this approach for accurate estimation, we must overcome an even bigger problem; we need to construct regeneration events of the queueing process (not only of the input process), which, in turn, requires the synchronization of a regeneration point of the input process with an empty state of the system. For this reason, we instead will base our analysis on the construction of the artificial regeneration obtained by
exponential splitting. (For a deeper understanding of the splitting method, we refer readers to the fundamental monograph [
13].) A key property opening the possibility of this construction in the initially non-regenerative process is that the interarrival distributions of the component processes are assumed to be
heavy-tailed. Thus, in this research, we focused on the case where the superposition consists of
Poisson inputs and
inputs with Pareto interarrival times. The splitting approach, in keeping the distributional properties of the original process, allows one to use regenerative simulation for accurate estimations of the required stationary performance metrics. To the best of our knowledge, this approach is realized for the first time in systems with a superposed input process, and this finding is the main contribution of this research.
Our second main purpose is to study the basic system by means of another
classically regenerative system with a single renewal input process. It is not a new approach as, for instance, it has been realized in a series of papers, such as [
3,
5,
7]. In these works, the so-called stationary-interval method (based on the Palm approach) was used to approximate the distribution of the interarrival times in the superposition process by a stationary distribution with independent identically distributed (iid) intervals. We emphasize that such a replacement is indeed an approximation because the successive intervals in the superposed process are dependent in general. The latter issue has been analyzed in [
1], in which the correlation structure of the successive interevent times was considered in detail. Another motivation of the approximation, as we mentioned, is that the target processes in the basic system (for example, workload or queue size) are not regenerative; however, a regenerative structure is highly desirable both to establish the stability of the system and to accurately estimate its stationary performance metrics. An issue in the study of queueing systems with superposed input (with independent renewal components) concerns the approximation of such a process by a Poisson process. It is a well-known problem going back to A. Khintchine [
14]. More exactly, when the number of summands (components) increases, while the (properly defined) traffic intensity
remains bounded (and not close to 1), then the superposed process approaches a Poisson process, with the rate being the sum of the rates of the summands. On the other hand, it has been shown in [
2] that, as the number
n of the summands increases, this approximation works well only if the product
. The latter result indicates that the Poissonian approximation is not suitable when the number of components is small and/or the traffic intensity
is close to 1, that is, the system is close to ‘saturation’. The saturated system (in the so-called
heavy-traffic regime) with the superposed input process is studied in the paper [
15]. Notice also that the basic results for ordinary renewal process, obtained originally by W. Smith [
16] and D. Cox [
17], have been extended to the superposition of independent renewal processes in [
4].
Thus, our second important contribution is the detailed analysis of the approximating
-type system when the interarrival time distribution is a mixture containing Pareto distributions, and when the service time distribution is also a mixture of distributions that describe the service times of the customers in the component processes. Indeed, we show that, in the approximating system, the stationary workload distribution is analytically available when the service times of the different classes of customers have a common exponential distribution (
Section 4.2); this is also a new contribution of this research. Moreover, we compare the accuracy of the approximation in terms of the Kolmogorov distance between the (empirical) distributions of the stationary waiting times in the original and approximated systems (
Section 6.1). In particular, in spite of the presence of heavy-tailed components in interarrival distribution, the simulation results show quite acceptable distance values (less than
) for most of the considered cases (when the traffic intensity
), and it is consistent with the results of the abovementioned works [
3,
5,
7]. Hence, the proposed approximation may be acceptable if the splitting construction turns out to be too complicated, e.g., if the number of component processes is very large. Another possible application of the approximating system is as a base for appropriate lower/upper bounds for the stationary workload process, exploiting the known monotonicity properties of the classical queueing system under some mild assumptions. (In this regard, see [
18] and recent papers of one of the authors [
19,
20] on the monotonicity of systems with mixed service times.)
As mentioned above, the approximation of the superposed input process by a renewal process is indeed a mixed distribution, consisting of an appropriately weighted sum of the stationary distributions of the remaining renewal times (integrated-tail distributions) in the component stationary renewal processes. This distribution, based on the Palm probability concept, is well known (see, for instance, [
7,
17]) and is the exact one-dimensional distribution of the interarrival times, which, however, does not capture possible dependencies between the neighboring intervals. For this reason, reducing the superposed input process to an ordinary renewal process is indeed an approximation. We note that, in general, a dependence between the adjacent intervals cannot be ignored in the careful study of the superposed input, as shown in [
5].
The rest of the paper is organized as follows. In
Section 2, we give some important preliminary results from the theory of renewal and point processes with elements of the Palm theory. In
Section 3, we outline the construction of
one-dependent regeneration using the approach from [
12] and highlight the difficulties arising in its practical realization in simulation. Nevertheless, we use some basic ideas of this approach to construct regeneration events by the exponential splitting of heavy-tailed interarrival times.
In
Section 4, we examine an approximating
system with iid interarrival times, generated by the Palm distribution of the stationary intervals in the original superposed input process. We then focus on special cases where we find explicit expressions for the Laplace–Stieltjes transform (LST) of the stationary workload. In
Section 5, we introduce exponential splitting and explain how it is used to obtain classical regenerations in the superposed process generated by Pareto components. Next, in
Section 6, we analyze our theoretical findings through discrete-event simulations. Specifically, we first evaluate the accuracy of the input process approximation in terms of Kolmogorov distance for different system parameter settings. Then, we apply exponential splitting to construct a confidence estimate of the empirical mean of the stationary waiting time in the original system. Finally,
Section 7 concludes the paper by providing an overview of its main contributions.
2. Preliminary Results from the Theory of Renewal and Point Processes
In this section, we present theoretical results from the theory of renewal and point processes [
21,
22]. These results are then utilized in the analysis of queueing systems with the superposed input process generated by independent renewal processes. Among the most important sources on renewal theory, besides the pioneering work [
23], we also mention the recognized monographs [
24,
25].
To start, let us construct a stationary renewal process. To do this, consider iid non-negative random variables (r.v.s)
and an independent r.v.
. We denote by
F the distribution function of an r.v.
X, which is a stochastic copy of any
. It is assumed that the renewal rate
. Let
be the distribution of the first r.v.
. Also, denote by
the renewal instants, where
. Hence,
represents the number of renewal events that occur in the time interval
, with
. The process
is a time-stationary renewal process [
22,
26]. Specifically, the stationary renewal process satisfies
and the distribution of the increments
does not depend on
for each fixed
There are other equivalent statements to identify a stationary renewal process (see, for example, Chapter 2 in [
25]). It is also known that (
1) represents the distribution of the stationary remaining (and attained) renewal time, and it is referred to as the integrated-tail distribution [
12].
In the following analysis, we will consider n independent stationary renewal processes that model the arrivals of different classes of customers. Then, the renewal points correspond to the arrival instants of class-i customers from process i, and let denote the generic interarrival time (with a rate of ) in process i, where . These processes will compose a merged superposed input process. To provide a more detailed description, let be the overall set of renewal points for all n processes, and let be the ordered sequence of these renewal instants. In other words, is the first smallest element of , is the second smallest element, and so on.
Assuming, for simplicity, that there are no batch arrivals, it follows that the sequence
satisfies
It is clear that, unlike a renewal process, the distances between events (arrivals) in the superposed process are, in general, dependent.
Since the component renewal processes are assumed to be time-stationary, then the superposed process, defined (in evident notation) as
is also time-stationary (counting), but not renewal in general. (A detailed proof can be found, for example, in [
21], Chapter 1.)
To explain this process in more detail, we recall some findings from the theory of stationary point processes, while maintaining the previously established notation. Following [
21,
25], we consider a set of points (occurrence times)
on the positive real line that satisfy (
2), with the condition that
is an occurrence time. This process is assumed to be interval-stationary and is referred to as the stationary Palm process. We also consider a counting process
, which records the number of points in the interval
, formally defined as
with the assumption that
. Let
denote the interval between the
kth and
st points, where
. These intervals are identically distributed with distribution function
F. This new process is also referred to as the stationary Palm process [
7,
21,
22].
In the following discussion, we will distinguish between the Palm process (
4) and the related time-stationary (also counting) process
. Note that an interval sequence
(forming the Palm process) is stationary if the joint distribution of the
k-tuple
is independent of the positive integer
h for all integers
and all non-negative
. It is worth pointing out that the Palm process
and the associated sequence
are usually not both stationary at the same time [
7]. However, there is a one-to-one correspondence between the stationary counting process
and the stationary sequence
based on the Palm theory [
21]. To establish this correspondence, let us introduce the
Palm function , which represents the probability of having
k points (
) in a time interval of length
t provided that the beginning of the interval is an occurrence point. Hence, Palm functions describe the dynamics of the Palm process (
4). The distribution of the process
associated with the Palm process (
4) can be expressed via Palm functions by the following relations [
7,
14]:
implying, in particular, the following key connection:
where, as already stated,
F is the distribution function of the intervals in the stationary Palm process. (In what follows, for any distribution function
F, we denote by
its tail.) Note that, by construction, the stationary process
satisfies both required conditions (see [
7]):
Recall that
represents the
ith component (stationary renewal process) and we will supply with index
i the corresponding quantities in process
i. Note that, by (
1), the first interval in the process
i has distribution
and rate
. Let
It is known [
7,
14,
21] that the Palm functions
of the superposed stationary Palm process
are expressed via the original time-stationary processes
and their stationary Palm versions
as follows:
where the summation is taken over all different partitions of
such that
. In particular, by (
8), the tail of the distribution
F satisfies [
7]:
It is easy to verify that
, as required. We emphasize that the Palm process
is interval-stationary, while the process
is time-stationary. Furthermore, since, for a single component process with interevent distribution
F and rate
, we have
Equation (
5) implies that
, satisfying (
6). Additionally, expression (
10) has a straightforward intuitive interpretation.
Since a renewal process is a particular case of the point process, it is worth mentioning that
F represents the exact marginal (one-dimensional) distribution of the interarrival time in the superposed Palm process
. However, this distribution does not capture the dependence between adjacent intervals. As a result, expression (
10) only provides an approximation of the original superimposed process, embedded at the instances of the points’ occurrences.
Therefore, it is crucial to examine the accuracy of this approximation using various metrics. Specifically, as detailed in
Section 6, the approximation (
10) of the original
system yields an acceptable Kolmogorov distance between the empirical distributions of the stationary waiting times when
.
Remark 1. Representations (10), (9) hold for any superposition in which the components , are time-stationary processes, not necessary renewal. In this general case, the intensities are defined as [7]. 3. Splitting of a Superposed Process
In general, for superposed processes, it is only possible to construct the so-called one-dependent regenerations. In this section, we will briefly describe this construction following [
12]. Roughly speaking, in a one-dependent regenerative process, two consecutive regeneration cycles are dependent, but the cycles separated by at least one cycle are independent. However, as we will discuss below, this approach is challenging to apply in simulations. Therefore, in
Section 5, we will introduce exponential splitting, which can be effectively applied when the renewal intervals follow heavy-tailed distributions [
22].
Let us describe the construction of one-dependent regenerations for the superposition of two renewal processes, with inter-renewal time distributions
. Firstly, we assume that both these processes are positive recurrent, i.e., the mean inter-renewal time is finite for both processes. Furthermore, let us assume that
satisfies the following condition:
for some positive constants
, and
c. This means that
is a particular case of a spread-out distribution
F. (In the general case, some convolution of
F with itself must satisfy (
11); see [
22].)
Denote by
the (right-continuous) remaining renewal time at instant
t in the
ith renewal process,
. Then, under assumption (
11), the distribution of
has the following smoothness property:
for each Borel set
B on
, where
is the uniform measure on
,
, and
C is some constant. (Condition (
12) is sometimes called minorization.) That is, for large enough
t, the distribution of
has a uniform component. Let
and denote by
the renewal points of the process
. Let
where
C satisfies (
12). Also, denote by
the iid r.v.s with uniform distribution
(in what follows, we will denote it as
.) Also, let
be iid
Bernoulli r.v.s such that
. Take
where the r.v.
has the rest distribution (see (
12))
Continuing as in (
13), define recursively
where
are independent of all
. Eventually, let
Then, the sequence
and moreover, at each instant
, the remaining times vector has distribution
which is independent of the pre-history prior to instant
. (Here,
denotes the distribution generated by the marginal distributions
and
.) Thus,
are one-dependent regeneration points for the superposed remaining time process
, and, according to (
14)–(
16), they constitute a subsequence of the renewal points of input 2. The extension of this construction to arbitrary
n renewal processes is straightforward [
12].
Remark 2. It is worth pointing out that, under classical regeneration, the distribution of in (17) would be independent of the pre-history prior to instant . The primary challenge in practically implementing regenerations in the superposed process is that the constants in Equation (
12) are not explicitly defined. An even more difficult problem is the need to synchronize these regeneration points with empty states of the queueing system to obtain regenerations of the entire system. However, in some cases, particularly when the inter-renewal intervals follow heavy-tailed distributions, this difficulty can be overcome. In this case, as we will show in
Section 5, by utilizing the aforementioned construction of one-dependent regenerations, the splitting procedure can be effectively applied to construct the classical regenerations of the superposed process.