Source author record

Bert Zwart

Bert Zwart appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

35works
11topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

35 published item(s)

preprint2026arXiv

Heavy tails in dynamic flow networks: Universal explanation of their emergence

Overload-induced cascading failures can cause extreme disruptions in a wide range of networked systems, such as power grids, transportation networks, or financial systems. Empirical studies across domains report that the size of such disruptions often follows a Pareto- or heavy-tailed distribution. While many models reproduce this scaling behavior, they are either tailored to specific domains or based on simplified mechanisms that overlook key aspects of overload cascading behavior. Hence, a general understanding of the mechanisms driving scale-free behavior in these settings remains incomplete. In this paper, we develop a universal and analytically tractable model of overload cascading failures on flow networks, offering a new perspective on how Pareto-tailed disruptions emerge across networks. Our framework shows, under mild assumptions, that heavy-tailed disruptions can arise naturally from Pareto-tailed external inputs, and it establishes a transformation law linking the input and output tail exponents. We further identify broad conditions under which the resulting cascade cost exhibits a heavy-tailed distribution and show that the mechanism is robust across several domains, including power transmission, traffic networks, and processing systems. Our results provide a unified explanation for the emergence of scale-free failures in overload-driven systems and connect previously disparate, application-specific models under a unified framework.

preprint2022arXiv

Persistence of heavy-tailed sample averages: principle of infinitely many big jumps

We consider the sample average of a centered random walk in $\mathbb{R}^d$ with regularly varying step size distribution. For the first exit time from a compact convex set $A$ not containing the origin, we show that its tail is of lognormal type. Moreover, we show that the typical way for a large exit time to occur is by having a number of jumps growing logarithmically in the scaling parameter.

preprint2022arXiv

Tail Asymptotics for the Delay in a Brownian Fork-Join Queue

In this paper, we study the tail behavior of $\max_{i\leq N}\sup_{s>0}\left(W_i(s)+W_A(s)-βs\right)$ as $N\to\infty$, with $(W_i,i\leq N)$ i.i.d. Brownian motions and $W_A$ an independent Brownian motion. This random variable can be seen as the maximum of $N$ mutually dependent Brownian queues, which in turn can be interpreted as the backlog in a Brownian fork-join queue. In previous work, we have shown that this random variable centers around $\frac{σ^2}{2β}\log N$. Here, we analyze the rare-event that this random variable reaches the value $(\frac{σ^2}{2β}+a)\log N$, with $a>0$. It turns out that its probability behaves roughly as a power law with $N$, where the exponent depends on $a$. However, there are three regimes, around a critical point $a^{\star}$; namely, $0<a<a^{\star}$, $a=a^{\star}$, and $a>a^{\star}$. The latter regime exhibits a form of asymptotic independence, while the first regime reveals highly irregular behavior with a clear dependence structure among the $N$ suprema, with a nontrivial transition at $a=a^{\star}$.

preprint2020arXiv

A Fluid Model of an Electric Vehicle Charging Network

We develop and analyze a measure-valued fluid model keeping track of parking and charging requirements of electric vehicles in a local distribution grid. We show how this model arises as an accumulation point of an appropriately scaled sequence of stochastic network models. The invariant point of the fluid model encodes the electrical characteristics of the network and the stochastic behavior of its users, and it is characterized, when it exists, by the solution of a so-called Alternating Current Optimal Power Flow (ACOPF) problem.

preprint2020arXiv

Analyzing large frequency disruptions in power systems using large deviations theory

We propose a method for determining the most likely cause, in terms of conventional generator outages and renewable fluctuations, of power system frequency reaching a predetermined level that is deemed unacceptable to the system operator. Our parsimonious model of system frequency incorporates primary and secondary control mechanisms, and supposes that conventional outages occur according to a Poisson process and renewable fluctuations follow a diffusion process. We utilize a large deviations theory based approach that outputs the most likely cause of a large excursion of frequency from its desired level. These results yield the insight that current levels of renewable power generation do not significantly increase system vulnerability in terms of frequency deviations relative to conventional failures. However, for a large range of model parameters it is possible that such vulnerabilities may arise as renewable penetration increases.

preprint2020arXiv

Computing first passage times for Markov-modulated fluid models using numerical PDE problem solvers

A popular method to compute first-passage probabilities in continuous-time Markov chains is by numerically inverting their Laplace transforms. Past decades, the scientific computing community has developed excellent numerical methods for solving problems governed by partial differential equations (PDEs), making the availability of a Laplace transform not necessary here for computational purposes. In this study we demonstrate that numerical PDE problem solvers are suitable for computing first passage times, and can be very efficient for this purpose. By doing extensive computational experiments, we show that modern PDE problem solvers can outperform numerical Laplace transform inversion, even if a transform is available. When the Laplace transform is explicit (e.g. does not require the computation of an eigensystem), numerical transform inversion remains the primary method of choice.

preprint2020arXiv

Consistency of the PLFit estimator for power-law data

We prove the consistency of the Power-Law Fit PLFit method proposed by Clauset et al.(2009) to estimate the power-law exponent in data coming from a distribution function with regularly-varying tail. In the complex systems community, PLFit has emerged as the method of choice to estimate the power-law exponent. Yet, its mathematical properties are still poorly understood. The difficulty in PLFit is that it is a minimum-distance estimator. It first chooses a threshold that minimizes the Kolmogorov-Smirnov distance between the data points larger than the threshold and the Pareto tail, and then applies the Hill estimator to this restricted data. Since the number of order statistics used is random, the general theory of consistency of power-law exponents from extreme value theory does not apply. Our proof consists in first showing that the Hill estimator is consistent for general intermediate sequences for the number of order statistics used, even when that number is random. Here, we call a sequence intermediate when it grows to infinity, while remaining much smaller than the sample size. The second, and most involved, step is to prove that the optimizer in PLFit is with high probability an intermediate sequence, unless the distribution has a Pareto tail above a certain value. For the latter special case, we give a separate proof.

preprint2020arXiv

Emergence of scale-free blackout sizes in power grids

We model power grids as graphs with heavy-tailed sinks, which represent demand from cities, and study cascading failures on such graphs. Our analysis links the scale-free nature of blackout sizes to the scale-free nature of city sizes, contrasting previous studies suggesting that this nature is governed by self-organized criticality. Our results are based on a new mathematical framework combining the physics of power flow with rare event analysis for heavy-tailed distributions, and are validated using various synthetic networks and the German transmission grid.

preprint2020arXiv

Optimization of stochastic lossy transport networks and applications to power grids

Motivated by developments in renewable energy and smart grids, we formulate a stylized mathematical model of a transport network with stochastic load fluctuations. Using an affine control rule, we explore the trade-off between the number of controllable resources in a lossy transport network and the performance gain they yield in terms of expected power losses. Our results are explicit and reveal the interaction between the level of flexibility, the intrinsic load uncertainty and the network structure.

preprint2020arXiv

Ranking transmission lines by overload probability using the empirical rate function

We develop a non-parametric procedure for ranking transmission lines in a power system according to the probability that they will overload due to stochastic renewable generation or demand-side load fluctuations, and compare this procedure to several benchmark approaches. Using the IEEE 39-bus test network we provide evidence that our approach, which statistically estimates the rate function for each line, is highly promising relative to alternative methods which count overload events or use incorrect parametric assumptions.

preprint2016arXiv

A Dirichlet Process Characterization of RBM in a Wedge

Reflected Brownian motion (RBM) in a wedge is a 2-dimensional stochastic process Z whose state space in R^2 is given in polar coordinates by S={(r,theta): r >= 0, 0 <= theta <= xi} for some 0 < xi < 2 pi. Let alpha= (theta_1+theta_2)/xi, where -pi/2 < theta_1,theta_2 < pi/2 are the directions of reflection of Z off each of the two edges of the wedge as measured from the corresponding inward facing normal. We prove that in the case of 1 < alpha < 2, RBM in a wedge is a Dirichlet process. Specifically, its unique Doob-Meyer type decomposition is given by Z=X+Y, where X is a two-dimensional Brownian motion and Y is a continuous process of zero energy. Furthermore, we show that for p > alpha , the strong p-variation of the sample paths of Y is finite on compact intervals, and, for 0 < p <= alpha, the strong p-variation of Y is infinite on [0,T] whenever Z has been started from the origin. We also show that on excursion intervals of Z away from the origin, (Z,Y) satisfies the standard Skorokhod problem for X. However, on the entire time horizon (Z,Y) does not satisfy the standard Skorokhod problem for X, but nevertheless we show that it satisfies the extended Skorkohod problem.

preprint2016arXiv

Fluid Limit of a PS-queue with Multistage Service

The PS-model treated in this paper is motivated by freelance job websites where multiple freelancers compete for a single job. In the context of such websites, multistage service of a job means collection of applications from multiple freelancers. Under Markovian stochastic assumptions, we develop fluid limit approximations for the PS-model in overload. Based on this approximation, we estimate what proportion of freelancers get the jobs they apply for. In addition, the PS-model studied here is an instant of PS with routing and impatience, for which no Lyapunov function is known, and we suggest some partial solutions.

preprint2016arXiv

Importance sampling of heavy-tailed iterated random functions

We consider a stochastic recurrence equation of the form $Z_{n+1} = A_{n+1} Z_n+B_{n+1}$, where $\mathbb{E}[\log A_1]<0$, $\mathbb{E}[\log^+ B_1]<\infty$ and $\{(A_n,B_n)\}_{n\in\mathbb{N}}$ is an i.i.d. sequence of positive random vectors. The stationary distribution of this Markov chain can be represented as the distribution of the random variable $Z \triangleq \sum_{n=0}^\infty B_{n+1}\prod_{k=1}^nA_k$. Such random variables can be found in the analysis of probabilistic algorithms or financial mathematics, where $Z$ would be called a stochastic perpetuity. If one interprets $-\log A_n$ as the interest rate at time $n$, then $Z$ is the present value of a bond that generates $B_n$ unit of money at each time point $n$. We are interested in estimating the probability of the rare event $\{Z>x\}$, when $x$ is large; we provide a consistent simulation estimator using state-dependent importance sampling for the case, where $\log A_1$ is heavy-tailed and the so-called Cramér condition is not satisfied. Our algorithm leads to an estimator for $P(Z>x)$. We show that under natural conditions, our estimator is strongly efficient. Furthermore, we extend our method to the case, where $\{Z_n\}_{n\in\mathbb{N}}$ is defined via the recursive formula $Z_{n+1}=Ψ_{n+1}(Z_n)$ and $\{Ψ_n\}_{n\in\mathbb{N}}$ is a sequence of i.i.d. random Lipschitz functions.

preprint2016arXiv

Line failure probability bounds for power grids

We develop upper bounds for line failure probabilities in power grids, under the DC approximation and assuming Gaussian noise for the power injections. Our upper bounds are explicit, and lead to characterization of safe operational capacity regions that are convex and polyhedral, making our tools compatible with existing planning methods. Our probabilistic bounds are derived through the use of powerful concentration inequalities.

preprint2016arXiv

Minimizing heat loss in DC networks using batteries

Electricity transmission networks dissipate a non-negligible fraction of the power they transport due to the heat loss in the transmission lines. In this work we explore how the transport of energy can be more efficient by adding to the network multiple batteries that can coordinate their operations. Such batteries can both charge using the current excess in the network or discharge to meet the network current demand. Either way, the presence of batteries in the network can be leveraged to mitigate the intrinsic uncertainty in the power generation and demand and, hence, transport the energy more efficiently through the network. We consider a resistive DC network with stochastic external current injections or consumptions and show how the expected total heat loss depends on the network structure and on the batteries operations. Furthermore, in the case where the external currents are modeled by Ornstein-Uhlenbeck processes, we derive the dynamical optimal control for the batteries over a finite time interval.

preprint2016arXiv

Transient error approximation in a Lévy queue

Motivated by a capacity allocation problem within a finite planning period, we conduct a transient analysis of a single-server queue with Lévy input. From a cost minimization perspective, we investigate the error induced by using stationary congestion measures as opposed to time-dependent measures. Invoking recent results from fluctuation theory of Lévy processes, we derive a refined cost function, that accounts for transient effects. This leads to a corrected capacity allocation rule for the transient single-server queue. Extensive numerical experiments indicate that the cost reductions achieved by this correction can by significant.

preprint2015arXiv

Insensitivity of Proportional Fairness in Critically Loaded Bandwidth Sharing Networks

Proportional fairness is a popular service allocation mechanism to describe and analyze the performance of data networks at flow level. Recently, several authors have shown that the invariant distribution of such networks admits a product form distribution under critical loading. Assuming exponential job size distributions, they leave the case of general job size distributions as an open question. In this paper we show the conjecture holds for a dense class of distributions. This yields a key example of a stochastic network in which the heavy traffic limit has an invariant distribution that does not depend on second moments. Our analysis relies on a uniform convergence result for a fluid model which may be of independent interest.

preprint2014arXiv

Corrected phase-type approximations for the workload of the MAP/G/1 queue with heavy-tailed service times

In many applications, significant correlations between arrivals of load-generating events make the numerical evaluation of the load of a system a challenging problem. Here, we construct very accurate approximations of the workload distribution of the MAP/G/1 queue that capture the tail behavior of the exact workload distribution and provide a small relative error. Motivated by statistical analysis, we assume that the service times are a mixture of a phase-type and a heavy-tailed distribution. With the aid of perturbation analysis, we derive our approximations as a sum of the workload distribution of the MAP/PH/1 queue and a heavy-tailed component that depends on the perturbation parameter. We refer to our approximations as corrected phase-type approximations, and we exhibit their performance with a numerical study.

preprint2014arXiv

Corrected phase-type approximations of heavy-tailed queueing models in a Markovian environment

Significant correlations between arrivals of load-generating events make the numerical evaluation of the workload of a system a challenging problem. In this paper, we construct highly accurate approximations of the workload distribution of the MAP/G/1 queue that capture the tail behavior of the exact workload distribution and provide a bounded relative error. Motivated by statistical analysis, we consider the service times as a mixture of a phase-type and a heavy-tailed distribution. With the aid of perturbation analysis, we derive our approximations as a sum of the workload distribution of the MAP/PH/1 queue and a heavy-tailed component that depends on the perturbation parameter. We refer to our approximations as corrected phase-type approximations, and we exhibit their performance with a numerical study.

preprint2014arXiv

Corrected phase-type approximations of heavy-tailed queueing models in a Markovian environment

We develop accurate approximations of the delay distribution of the MArP/G/1 queue that cap- ture the exact tail behavior and provide bounded relative errors. Motivated by statistical analysis, we consider the service times as a mixture of a phase-type and a heavy-tailed distribution. With the aid of perturbation analysis, we derive corrected phase-type approximations as a sum of the delay in an MArP/PH/1 queue and a heavy-tailed component depending on the perturbation parameter. We exhibit their performance with numerical examples.

preprint2014arXiv

Corrected phase-type approximations of heavy-tailed risk models using perturbation analysis

Numerical evaluation of performance measures in heavy-tailed risk models is an important and challenging problem. In this paper, we construct very accurate approximations of such performance measures that provide small absolute and relative errors. Motivated by statistical analysis, we assume that the claim sizes are a mixture of a phase-type and a heavy-tailed distribution and with the aid of perturbation analysis we derive a series expansion for the performance measure under consideration. Our proposed approximations consist of the first two terms of this series expansion, where the first term is a phase-type approximation of our measure. We refer to our approximations collectively as corrected phase-type approximations. We show that the corrected phase-type approximations exhibit a nice behavior both in finite and infinite time horizon, and we check their accuracy through numerical experiments.

preprint2014arXiv

On the accuracy of phase-type approximations of heavy-tailed risk models

Numerical evaluation of ruin probabilities in the classical risk model is an important problem. If claim sizes are heavy-tailed, then such evaluations are challenging. To overcome this, an attractive way is to approximate the claim sizes with a phase-type distribution. What is not clear though is how many phases are enough in order to achieve a specific accuracy in the approximation of the ruin probability. The goals of this paper are to investigate the number of phases required so that we can achieve a pre-specified accuracy for the ruin probability and to provide error bounds. Also, in the special case of a completely monotone claim size distribution we develop an algorithm to estimate the ruin probability by approximating the excess claim size distribution with a hyperexponential one. Finally, we compare our approximation with the heavy traffic and heavy tail approximations.

preprint2014arXiv

Separation of timescales in a two-layered network

We investigate a computer network consisting of two layers occurring in, for example, application servers. The first layer incorporates the arrival of jobs at a network of multi-server nodes, which we model as a many-server Jackson network. At the second layer, active servers at these nodes act now as customers who are served by a common CPU. Our main result shows a separation of time scales in heavy traffic: the main source of randomness occurs at the (aggregate) CPU layer; the interactions between different types of nodes at the other layer is shown to converge to a fixed point at a faster time scale; this also yields a state-space collapse property. Apart from these fundamental insights, we also obtain an explicit approximation for the joint law of the number of jobs in the system, which is provably accurate for heavily loaded systems and performs numerically well for moderately loaded systems. The obtained results for the model under consideration can be applied to thread-pool dimensioning in application servers, while the technique seems applicable to other layered systems too.

preprint2014arXiv

Time-dependent behaviour of an alternating service queue

We consider a model describing the waiting time of a server alternating between two service points. This model is described by a Lindley-type equation. We are interested in the time-dependent behaviour of this system and derive explicit expressions for its time-dependent waiting-time distribution, the correlation between waiting times, and the distribution of the cycle length. Since our model is closely related to Lindley's recursion, we compare our results to those derived for Lindley's recursion.

preprint2013arXiv

Fluid Approximation of a Call Center Model with Redials and Reconnects

In many call centers, callers may call multiple times. Some of the calls are re-attempts after abandonments (redials), and some are re-attempts after connected calls (reconnects). The combination of redials and reconnects has not been considered when making staffing decisions, while ignoring them will inevitably lead to under- or overestimation of call volumes, which results in improper and hence costly staffing decisions. Motivated by this, in this paper we study call centers where customers can abandon, and abandoned customers may redial, and when a customer finishes his conversation with an agent, he may reconnect. We use a fluid model to derive first order approximations for the number of customers in the redial and reconnect orbits in the heavy traffic. We show that the fluid limit of such a model is the unique solution to a system of three differential equations. Furthermore, we use the fluid limit to calculate the expected total arrival rate, which is then given as an input to the Erlang A model for the purpose of calculating service levels and abandonment rates. The performance of such a procedure is validated in the case of single intervals as well as multiple intervals with changing parameters.

preprint2013arXiv

Fluid Limits for Bandwidth-Sharing Networks with Rate Constraints

Bandwidth-sharing networks as introduced by Massoulié & Roberts (1998) model the dynamic interaction among an evolving population of elastic flows competing for several links. With policies based on optimization procedures, such models are of interest both from a Queueing Theory and Operations Research perspective. In the present paper, we focus on bandwidth-sharing networks with capacities and arrival rates of a large order of magnitude compared to transfer rates of individual flows. This regime is standard in practice. In particular, we extend previous work by Reed & Zwart (2010) on fluid approximations for such networks: we allow interarrival times, flow sizes and patient times (i.e. abandonment times measured from the arrival epochs) to be generally distributed, rather than exponentially distributed. We also develop polynomial-time computable fixed-point approximations for stationary distributions of bandwidth-sharing networks, and suggest new techniques for deriving these types of results.

preprint2013arXiv

Random Fluid Limit of an Overloaded Polling Model

In the present paper, we study the evolution of an overloaded cyclic polling model that starts empty. Exploiting a connection with multitype branching processes, we derive fluid asymptotics for the joint queue length process. Under passage to the fluid dynamics, the server switches between the queues infinitely many times in any finite time interval causing frequent oscillatory behavior of the fluid limit in the neighborhood of zero. Moreover, the fluid limit is random. Additionally, we suggest a method that establishes finiteness of moments of the busy period in an M/G/1 queue.

preprint2013arXiv

Scaling limits via excursion theory: Interplay between Crump-Mode-Jagers branching processes and processor-sharing queues

We study the convergence of the $M/G/1$ processor-sharing, queue length process in the heavy traffic regime, in the finite variance case. To do so, we combine results pertaining to Lévy processes, branching processes and queuing theory. These results yield the convergence of long excursions of the queue length processes, toward excursions obtained from those of some reflected Brownian motion with drift, after taking the image of their local time process by the Lamperti transformation. We also show, via excursion theoretic arguments, that this entails the convergence of the entire processes to some (other) reflected Brownian motion with drift. Along the way, we prove various invariance principles for homogeneous, binary Crump-Mode-Jagers processes. In the last section we discuss potential implications of the state space collapse property, well known in the queuing literature, to branching processes.

preprint2012arXiv

Efficient Rare-event Simulation for Perpetuities

We consider perpetuities of the form D = B_1 exp(Y_1) + B_2 exp(Y_1+Y_2) + ... where the Y_j's and B_j's might be i.i.d. or jointly driven by a suitable Markov chain. We assume that the Y_j's satisfy the so-called Cramer condition with associated root theta_{ast} in (0,infty) and that the tails of the B_j's are appropriately behaved so that D is regularly varying with index theta_{ast}. We illustrate by means of an example that the natural state-independent importance sampling estimator obtained by exponentially tilting the Y_j's according to theta_{ast} fails to provide an efficient estimator (in the sense of appropriately controlling the relative mean squared error as the tail probability of interest gets smaller). Then, we construct estimators based on state-dependent importance sampling that are rigorously shown to be efficient.

preprint2012arXiv

Fluid Limits for an ALOHA-type Model with Impatient Customers

Random multiple-access protocols of type ALOHA are used to regulate networks with a star configuration where client nodes talk to the hub node at the same frequency (finding a wide range of applications among telecommunication systems, including mobile telephone networks and WiFi networks). Such protocols control who talks at what time sharing the common idea "try to send your data and, if your message collides with another transmission, try resending later". In the present paper, we consider a time-slotted ALOHA model where users are allowed to renege before transmission completion. We focus on the scenario that leads to overload in the absence of impatience. Under mild assumptions, we show that the fluid (or law-of-large-numbers) limit of the system workload coincides a.s. with the unique solution to a certain integral equation. We also demonstrate that the fluid limits for distinct initial conditions converge to the same value as time tends to infinity.

preprint2011arXiv

A Lévy input fluid queue with input and workload regulation

We consider a queuing model with the workload evolving between consecutive i.i.d.\ exponential timers $\{e_q^{(i)}\}_{i=1,2,...}$ according to a spectrally positive Lévy process $Y_i(t)$ that is reflected at zero, and where the environment $i$ equals 0 or 1. When the exponential clock $e_q^{(i)}$ ends, the workload, as well as the Lévy input process, are modified; this modification may depend on the current value of the workload, the maximum and the minimum workload observed during the previous cycle, and the environment $i$ of the Lévy input process itself during the previous cycle. We analyse the steady-state workload distribution for this model. The main theme of the analysis is the systematic application of non-trivial functionals, derived within the framework of fluctuation theory of Lévy processes, to workload and queuing models.

preprint2011arXiv

Convergence of the all-time supremum of a Lévy process in the heavy-traffic regime

In this paper we derive a technique of obtaining limit theorems for suprema of Lévy processes from their random walk counterparts. For each $a>0$, let $\{Y^{(a)}_n:n\ge 1\}$ be a sequence of independent and identically distributed random variables and $\{X^{(a)}_t:t\ge 0\}$ be a Lévy processes such that $X_1^{(a)}\stackrel{d}{=} Y_1^{(a)}$, $\mathbb E X_1^{(a)}<0$ and $\mathbb E X_1^{(a)}\uparrow0$ as $a\downarrow0$. Let $S^{(a)}_n=\sum_{k=1}^n Y^{(a)}_k$. Then, under some mild assumptions, $Δ(a)\max_{n\ge 0} S_n^{(a)}\stackrel{d}{\to} R\iffΔ(a)\sup_{t\ge 0} X^{(a)}_t\stackrel{d}{\to} R$, for some random variable $R$ and some function $Δ(\cdot)$. We utilize this result to present a number of limit theorems for suprema of Lévy processes in the heavy-traffic regime.

preprint2011arXiv

Diffusion limits of limited processor sharing queues

We consider a processor sharing queue where the number of jobs served at any time is limited to $K$, with the excess jobs waiting in a buffer. We use random counting measures on the positive axis to model this system. The limit of this measure-valued process is obtained under diffusion scaling and heavy traffic conditions. As a consequence, the limit of the system size process is proved to be a piece-wise reflected Brownian motion.

preprint2004arXiv

Exact asymptotics for fluid queues fed by multiple heavy-tailed on-off flows

We consider a fluid queue fed by multiple On-Off flows with heavy-tailed (regularly varying) On periods. Under fairly mild assumptions, we prove that the workload distribution is asymptotically equivalent to that in a reduced system. The reduced system consists of a ``dominant'' subset of the flows, with the original service rate subtracted by the mean rate of the other flows. We describe how a dominant set may be determined from a simple knapsack formulation. The dominant set consists of a ``minimally critical'' set of On-Off flows with regularly varying On periods. In case the dominant set contains just a single On-Off flow, the exact asymptotics for the reduced system follow from known results. For the case of several On-Off flows, we exploit a powerful intuitive argument to obtain the exact asymptotics. Combined with the reduced-load equivalence, the results for the reduced system provide a characterization of the tail of the workload distribution for a wide range of traffic scenarios.