Fluid Queue Driven by an M / M / 1 Queue Subject to Bernoulli-Schedule-Controlled Vacation and Vacation Interruption

This paper deals with the stationary analysis of a fluid queue driven by anM/M/1 queueing model subject to Bernoulli-ScheduleControlled Vacation and Vacation Interruption. The model under consideration can be viewed as a quasi-birth and death process. The governing system of differential difference equations is solved using matrix-geometric method in the Laplacian domain. The resulting solutions are then inverted to obtain an explicit expression for the joint steady state probabilities of the content of the buffer and the state of the background queueingmodel. Numerical illustrations are added to depict the convergence of the stationary buffer content distribution to one subject to suitable stability conditions.


Introduction
In many real time situations, the server in the background queueing model may become unavailable for a random period of time to perform a secondary task, when there are no customers in the waiting line at the service completion epoch.Such period of server absence is termed as server vacation.Queueing models subject to various vacation policies are of interest to researchers in recent times owing to their widespread applicability.
There are different types of vacation queueing systems.In the single vacation scheme, the server takes a vacation of some random duration when the queue is empty.At the end of the vacation, the server returns to the queue.The server resumes service if there is at least one customer waiting upon his return from vacation.However, if the queue is empty on the server's return, the server waits to complete a busy period.In the multiple vacation scheme, if the server returns from a vacation and finds the queue empty, he immediately commences another vacation.If there is at least one waiting customer, then he will commence the service.Queueing models subject to single or multiple exponential vacation are apt to model many practical scenarios [1][2][3].However, a better modeling assumption would be to assume that the server works at a slower rate during vacation periods in comparison to that of a regular working period.Such models are classified as queues subject to working vacations [4][5][6].In addition, the server can stop the vacation once some indices of the system, such as the number of customers, achieve a certain value in the vacation period.Certainly, it is possible for the server to take an interrupted vacation, so we call this policy vacation interruption.Li and Tian [7] studied the single server queueing model with working vacation and vacation interruption.The modulating queueing model considered in this paper is an //1 queue wherein the server is subject to regular vacation with probability  or working vacation with probability 1−.Further, the vacation duration of the server during working vacation epoch may be interrupted due to vacation interruption.
Fluid queues have become a fascinating area of research in recent years due to their widespread applicability in computer and communication systems [8,9], manufacturing systems [10], and so forth.A stochastic fluid flow system is an inputoutput model where the input is modeled as a continuous 2 Advances in Operations Research fluid that enters and leaves the storage device called a buffer, according to randomly varying rates.They are appropriate in a situation wherein the arrival is comprised of a discrete unit, but the interarrival time between successive arrivals is negligible.Therefore, the arrivals can be approximated by a continuous flow of fluid as individuals units have less impact on the performance of the system.In these models, a fluid buffer is either filled or depleted or both at rates determined by the current state of the background queueing model.Markov modulated fluid queues are a particular class of fluid models useful for modeling many physical phenomena and they often allow tractable analysis.In addition, fluid models are quite useful as approximate models for certain queueing and inventory systems where the flow consists of discrete entities, but the behavior of the individual is not important to identify the performance analysis.Certain interesting real world applications of Markov Modulated Fluid Flow models can be found in [11][12][13][14].Besides, fluid queues also have successful applications in the field of congestion control [15] and risk processes [16].More recently, Bosman and Nunez-Queija [17] considered a tandem fluid queue model to evaluate the performance of streaming media over an unreliable network.
For example, consider a production inventory model operating in a stochastic environment.The inventory level increases when the production rate exceeds the demand rate and decreases otherwise.The inventory level under continuous review can be viewed as a fluid process that fluctuates according to the evolution of the underlying background environment.For example, consider a machine shop with a single server.When the server is busy, items are produced continuously at a rate  and if he is idle, there is no production.However, for all practical reasons, the server might either take a vacation of some random duration with probability  or decide to provide service at a reduced rate with probability 1 − .Further, by offering service at a reduced rate, the server may continue to do so with probability  or due to certain unforeseen reasons, like a sudden increase in the demand, interrupt the vacation with probability 1 − , and continue the busy period.The demands are assumed to vary from time to time at the rate   independent of the state of the server.The level of inventory thus oscillates between  −   and −  depending on the busy or idle state of the server.Such scenario can be modeled as a fluid queue driven by an //1 queue subject to Bernoulli-Schedule-Controlled Vacation and Vacation Interruption.
The stationary analysis of fluid queueing models in a stochastic environment has been discussed by many authors.Fluid queues driven by an //1 queueing model are extensively studied in the literature.Various techniques have been employed by researchers to obtain the stationary buffer content distribution.To mention a few, Adan and Resing [18] analyze the buffer content distribution by viewing the arrival process as an alternating renewal process and Virtamo and Norros [19] provide the buffer content distribution by finding the spectrum of the eigenvalue equation and explicit expressions for the corresponding eigenvectors in terms of Chebyshev polynomials of the second kind.Sericola and Tuffin [20] express the stationary distribution of the buffer occupancy in terms of a sequence of recursively defined polynomials.Parthasarathy et al. [21] present an explicit expression for the buffer content distribution in terms of modified Bessel function of the first kind using continued fraction methodology.
Furthermore, fluid models driven by an //1 queue subject to various vacation strategies were analyzed in steady state by Mao et al. [22] and Wang et al. [23].The work was further extended to the stationary analysis of fluid queues driven by an //1 queue with multiple exponential vacation and  policy [24].Fluid model driven by an //1 queue with working vacations and vacation interruption was studied by Xu et al. [25].However, in most of the literature relating to fluid queues driven by vacation queueing models, the buffer content distribution is expressed in the Laplace domain.More recently, Vijayashree and Anjuka [26,27] presented an explicit expression for the buffer content distribution of a fluid queueing model modulated by an //1 queue subject to catastrophes and subsequent repair and a fluid queue driven by an / 2 /1 queueing model, respectively.Also, Ammar [28] derives an explicit expression for the fluid queue driven by an //1 queue with multiple exponential vacation using generating function methodology.
This paper presents an analytical solution for the fluid queue driven by an //1 queue subject to Bernoulli-Schedule-Controlled Vacation and Vacation Interruption in stationary regime.When the background queueing model is empty, the server will take either an ordinary vacation with probability  or a working vacation with probability 1 − .If the system is in a working vacation, upon completion of service, the server either ends the vacation and enters a regular busy period with probability 1 −  or continues the vacation with probability .It is assumed that the fluid in the buffer content increases at a constant rate, , when there is one or more customers in the background queueing model and decreases at a constant rate,  0 , when the queue is empty.The system of equations governing the process is modeled in terms of quasi-birth and death process and solved using matrix-geometric method.The stationary distribution of the buffer content is thereby obtained in the Laplace domain and hence inverted to obtain explicit expressions for the joint steady state probabilities of the state of the background queueing model and content of the buffer.Closed form expressions help to gain a deeper insight into the model.As a special case, when  = 0,  = 0, and  = 1 the theoretical results, so obtained, are seen to coincide with the existing results of Xu et al. [25] and Wang et al. [23], respectively.
The rest of the paper is organized as follows: Section 2 gives a brief description of the background queueing model under consideration.Section 3 presents the system of differential difference equations that governs the fluid queueing model under steady state subject to suitable stability conditions.Section 4 gives detailed derivations of the closed form expressions for the joint steady state probabilities of the state of the background queueing model and also the buffer content distributions.Section 5 presents the numerical illustration of the stationary buffer content distribution for suitable choice of the parameter values.

Model Description
Consider an //1 queueing model with infinite capacity.Let the customers arrive according to a Poisson process with parameter .The server provides service to the arriving customers according to an exponential distribution with parameter .When the system becomes empty, the server begins a vacation of random length and takes an ordinary vacation with probability  or a working vacation with probability 1 − , where 0 ≤  ≤ 1.In an ordinary vacation, the server will stop working even if there are new arrivals during the vacation period.In a working vacation, customers are served at a lower rate  V < .Further, in a working vacation, it is assumed that, at the instants of service completion, either the vacation is interrupted and the server resumes to a regular busy period with probability 1 −  or the server continues the vacation with probability .When the vacation period ends and the system is nonempty, a new busy period starts.The ordinary vacation times and the working vacation times are also assumed to be exponentially distributed with parameters  and  V , respectively.Let () denote the number of customers at time .Define 1, if the server is in an ordinary vacation at time , 2, if the server is in a regular busy period at time . ( It is well known that the process {((), ()),  ≥ 0} is a quasi-birth and death (QBD) process with state space given by Ω = {((0, 0) ∪ (0, 1)) ∪ (, ) ,  = 1, 2, . . .,  = 0, 1, 2} .
The state transition diagram of the background queueing model is given in Figure 1.Let represent the steady state probabilities of the background queueing model.Further, let  0 = ( 0,0  0,1 ) and   = ( ,0  ,1  ,2 ) for  ≥ 1.Then, the stationary probability vector is denoted by It is readily seen that the system of equations governing the background queueing model under steady state can be written in the form of matrix as where  1 = (1, 1)  ,  2 = (1, 1, 1)  , and  = ( ).
Note that The steady state probabilities of the single server queueing models with Poisson arrival and exponentially distributed service times subject to Bernoulli-Schedule-Controlled Vacation and Vacation Interruption were studied by Zhang and Shi [29].

Analysis of Fluid Queue
This section deals with the stationary analysis of a fluid queue modulated by an //1 queueing model subject to Bernoulli-Schedule-Controlled Vacation and Vacation Interruption.Let () be the content of the buffer at time .Furthermore, it is assumed that the content of the buffer increases at the rate , when there are customers in the background Variations of the content of the buffer queueing model, while the buffer content decreases at the rate  0 , when the system is empty.The rate at which the content of the buffer varies with time is given by

States of the background queueing model
where  0 < 0 and  > 0. Figure 2 depicts the interaction between the buffer content process and the background queueing model.It is seen that the content of the infinite capacity buffer decreases at the rate  0 < 0 when the background queueing model is empty with no waiting customers and it increases at the rate  when () ̸ = 0 irrespective of the states of ().

Solution Methodology
The stationary distributions of the fluid process play a vital role as they give more information relating to quantities of interest for practical applications like tail probabilities, expected buffer content, traffic intensity, expected delay, and sojourn time.This section presents explicit expressions for the joint steady state probabilities of the background queueing model and the content of the buffer in terms of modified Bessel function of the first kind.The governing system of equations in the Laplace domain is expressed as a system of matrix equations.The minimal nonnegative solution of the matrix quadratic equation is determined.The stationary joint probability distributions are expressed in terms of this minimal nonnegative solution and are shown to satisfy the governing system of matrix equations.This section presents an explicit analytical solution to the governing system of equations represented by (11).In this sequel, define  0 () = ( 0,0 () ,  0,1 ()) ,   () = ( ,0 () ,  ,1 () ,  ,2 ())  = 1, 2, . . . .
Our objective is to solve the above system of matrix difference equations to obtain explicit expressions for the stationary probability distribution and hence determine the stationary buffer content distribution.Towards that end, we present a lemma below followed by a theorem.) , where ) .

Advances in Operations Research
The solutions of ( 24) and ( 25) are given by Let () ( 1 ()) and () ( 1 ()) denote the negative (positive) roots of  11 () and  33 (), respectively.Considering the root that lies inside the unit circle, () and () are taken for further analysis.Further () satisfies the following relations: From ( 26), we get Substituting for  11 () and  33 () into ( 27) yields Using the relation given by (30) in the above leads to Again substituting for  22 () and  33 () into (28) leads to Using the relation given by (30) in the above equation yields This completes the proof.Note that   () for  = 1, 2, 3, . . .can be simplified as Also (0) = , where  is given by with  = / and  = ((+

Theorem 2. The stationary joint probability distributions of the content of the buffer and the state of the background queueing model in Laplace domain are given by
where  = ( 1 0 0 0 1 0 ) .
Having determined all the joint steady state probabilities in the Laplace domain, we now present the explicit analytical solution by inverting using transform techniques.From (41), we get which can be written as With  = 2√/ and  = 2√ V /, inversion of (46), (47), and (50) yields where Thus all the joint steady state probabilities of the state of the system and the content of the buffer are explicitly obtained in terms of modified Bessel function of the first kind.
Buffer Content Distribution.The stationary buffer content distribution of the fluid model under consideration is given by Taking Laplace transform of the above equation yields In matrix notation, the above equation can be rewritten as Then, which on inversion leads to where  0,0 () and  0,1 () are given by (51) and (52), respectively.Thus all the joint steady state probabilities and the buffer content distribution of the fluid queue driven by an //1 queue subject to Bernoulli-Schedule-Controlled Vacation and Vacation Interruption are explicitly obtained under steady state.Explicit analytical expressions help the practitioner to better understand the behavior of any quantity of interest, like the mean buffer content, for varying values of the parameters involved in the model.

Numerical Illustrations
This section illustrates the variation of the buffer content distribution against the content of the buffer for  = 1,  = 2,  V = 1.1,  = 0.9,  V = 0.6,  = 0.5,  = 0.5,  = 1 and varying values of  0 .The choice of  0 is relatively high as compared to  because of our assumptions that  happens when the background queueing model is nonempty and  0 happens otherwise.To compensate for the rarity in the occurrence of  0 , it is assumed to be larger.Figure 3 depicts the behavior of the buffer content distribution, (), against  for the above choice of the parameter values with  0 = −450.9.For this choice of the parameters, it is seen that  = −140.86< 0. Therefore, the stability condition is satisfied.It is seen that () increases  (68) Furthermore, as the value of  0 greatly affects buffer content distribution, its variation against  for a different value of  0 is presented in Table 1.

Conclusion
Markov Modulated Fluid Flows (MMFF) are a class of fluid models wherein the rates at which the content of the fluid buffer varies are modulated by the Markov process evolving in the background.This paper studies a fluid model driven by an //1 queue subject to Bernoulli-Schedule-Controlled Vacation and Vacation Interruption.The study of such models provides greater flexibility to the design and control of input and output rates of fluid flow thereby adapting the fluid models to wider application background.
The governing system of infinite differential difference equations is explicitly solved using Laplace transform and matrixgeometric methodology.Most of the existing results in the literature pertaining to MMFF have presented the solution to the buffer content distribution in the Laplace domain.However, closed form analytical solutions help to gain a deeper insight into the model and other related performance measures.The current findings can be thought of as one of the key contributions to the theoretical development of MMFF rather than the practical context.The theoretical results so obtained are verified with the existing results in the literature as a special case.The variations of the stationary buffer content distribution against the content of the buffer for varying values of  0 are numerically illustrated.

Figure 2 :
Figure 2: Interaction between fluid model and the background queueing model.

Table 1 :
Convergence of stationary buffer content distribution for varying values of  0 .