A Non-Markovian Multistage Batch Arrival Queue with Breakdown and

The present investigation deals with analysis of non-Markovian queueing model with multistage of services. When the server is unavailable during the systembreakdown (or) vacation periods, we consider reneging to prevail. Supplementary variable techniques have been adopted to obtain steady state system length distributions. The numerical illustrations are provided to validate the tractability of performance measures as far as computational aspect is concerned. Numerical results in the form of graphical representation are also presented. Practical large scale industry applications are described to justify our model.


Introduction
Vacation queueing models play a major role in manufacturing and production, computer and communication, and service and distribution systems.Many models for customer's impatience in queueing systems have been studied in the past, and the source of impatience has always been considered to be either a long wait already experienced at a queue or a long wait anticipated by a customer upon arrival.Queueing models enable organizations to implement critical production strategies or tactics aimed at reducing costs on increasing revenues.By examining the methods of operations research and especially queueing theory, new models can be developed and existing ones can be extended.
One of the earliest works on balking and reneging was by Haight [1,2] and Barrer [3], which was the first to introduce reneging in which they studied deterministic reneging with single server Markovian arrival and service rates.Montazer-Haghighi et al. [4] studied a Markovian mutiserver queueing system with balking and reneging.A two-stage batch arrival queueing system where customers receive a batch service in the first and individual service in the second stage was studied by Doshi [5] in the past.Furthermore [6], Chodhury examined an M [] /G/1 queueing system with a set-up period and a vacation period.This paper deals with an M  /G/1 queueing system with a vacation period which comprises an idle period and a random set-up period.The server is turned off each time when the system becomes empty.At this point of time the idle period starts.As soon as a customer or a batch of customers arrive, the setup of the service facility begins which is needed before starting each busy period.In this paper the steady-state behaviour of the queue size distributions at stationary (random) point of time and at departure point of time is studied.Also, explicit expressions for the system state probabilities and some performance measures of this queueing system are derived analytically.Finally, the probability generating function of the additional queue size distribution due to the vacation period as the limiting behaviour of the M  /M/1 type queueing system is derived.
Madan [7] investigated a batch arrival queueing system, where the server provides two stages of heterogeneous service with a modified Bernoulli schedule under N-policy.The server remains idle till the queue size becomes N (≥1).As soon as the queue size becomes at least N, the server instantly starts working and provides two stages of service in succession to each customer, that is, the first stage service followed by the second stage service.However, after the second stage service, the server may take a vacation or decide to stay in the system to provide service to the next customer, if any.The queue size distributions at a random epoch as well as a departure epoch under the steady-state conditions were derived.
The M/G/1 queue with impatient customers was studied by Bae et al. [8]; the complete formula of the limiting distribution of the virtual waiting time was derived explicitly.The expected busy period of the queue was also obtained by using a martingale argument.
In 2002, Krishna Kumar et al. [9] studied an M/G/1 retrial queueing system with two-phase service and preemptive resume.In this paper, for an arbitrarily distributed retrial time distribution, the necessary and sufficient condition for the system stability is obtained, assuming that only the customer at the head of the orbit has priority access to the server.The steady-state distributions of the server state and the number of customers in the orbit are obtained along with other performance measures.The effects of various parameters on the system performance are analysed numerically.A general decomposition law for this retrial queueing system is established.
Furthermore Madan et al. [10] studied the server vacations based on Bernoulli schedules and a single vacation policy in 2003.In this model, a single server queue with optional server vacations based on exhaustive service is studied.Unlike other vacation policies, the assumption that, only at the completion of service of the last customer in the system, the server has the option to take a vacation or to remain idle in the system waiting for the next customer to arrive is considered.The service times of the customers have been assumed to be deterministic and vacations are phase type exponential.Steady-state results for the probability generating functions of the queue length, the expected number of customers in the queue, and the expected waiting time of the customer are derived explicitly.
In addition, Choudhury [11] investigated a batch arrival queueing system with an additional service channel.This paper deals with an M  /G/1 queueing system with two phases of heterogeneous service under N-policy, where the server remains idle till the queue size becomes N (⩾1).As soon as the queue size becomes N, the server immediately starts first "essential service" for all the units.After completion of the essential service of a unit, it may leave the system with probability (1 − ) or may immediately go for a second phase of service in an additional service channel with probability  (0 ⩽  ⩽ 1).For this model, the queue size distribution at a random epoch as well as at a departure epoch is obtained.This is a generalization of recent papers considered by Medhi and Lee et al. (queueing systems).
In 2004 Madan and Abu Al-Rub [12] investigated on a single server queue with optional phase type server vacations based on exhaustive deterministic service and a single vacation policy.Furthermore Madan and Choudhury [13] proposed a queueing system with restricted admissibility of arriving batches.In this model a batch arrival queue with a Bernoulli vacation schedule is considered.After completion of the service the server either goes for a vacation in random length with probability  (0 ≤  ≤ 1) or may continue to serve the next unit, if any, with probability (1 − ) under a restricted admissibility policy of arriving batches.Unlike the usual batch arrival queueing system the restricted admissibility policy differs during a busy period and a vacation period and hence all arriving batches are not allowed to join the system at all time.Steady-state queue size distribution at a random point of time as well as a departure epoch is derived.Moreover this paper attempts to unify several classes of related batch arrival queueing system.
Choudhury and Madan [14] studied a two-stage batch arrival queueing system with a modified Bernoulli schedule vacation under N-policy.In this paper a batch arrival queueing system, where the server provides two stages of heterogeneous service with a modified Bernoulli schedule under N-policy, is considered.The server remains idle till the queue size becomes N (≥1).As soon as the queue size becomes at least N, the server instantly starts working and provides two stages of service in succession to each customer, that is, the first stage service followed by the second stage service.However, after the second stage service, the server may take a vacation or decide to stay in the system to provide service to the next customer, if any.The queue size distribution at a random epoch as well as a departure epoch under the steady-state conditions is derived.Further, the existence of the stochastic decomposition property to show that the departure point queue size distribution of this model can be decomposed into the distributions of three independent random variables is demonstrated.Finally, a simple procedure to obtain optimal stationary operating policy under a suitable linear cost structure is developed.
Furthermore Chang and Takine [15] studied factorization and stochastic decomposition properties in bulk queues with generalized vacation.This paper considers a class of stationary batch-arrival, bulk-service queues with generalized vacations.The system consists of a single server and a waiting room of infinite capacity.Arrivals of customers follow a batch Markovian arrival process.The server is unavailable for occasional intervals of time called vacations, and, when it is available, customers are served in groups of fixed size B. For this class of queues, the vector probability generating function of the stationary queue length distribution is factored into two terms, one of which is the vector probability generating function of the conditional queue length distribution given that the server is on vacation is shown.The special case of batch Poisson arrivals is examined and a new stochastic decomposition formula is derived for the stationary queue length distribution.
Altman and Yechiali [16] considered the systems with servers vacations where customers' impatience is due to an absentee of servers upon arrival.Such a model, representing frequent behavior by waiting customers in service systems, had never been treated before in the literature.
Altman and Yechiali [17] studied a system that is operating as an M/M/∞ queue.However, when it becomes empty, it is assigned to perform another task, the duration U of which is random.Customers arriving while the system is unavailable for service (i.e., occupied with a U-task) become impatient: each individual activates an "impatience timer" having random duration T such that if the system does not become available by the time the timer expires, the customer leaves the system never to return.When the system completes a U-task and there are waiting customers, each one is taken immediately into service.Explicit expressions for the corresponding mean queue sizes were obtained.
In 2012, Kumar and Sharma [18] analyzed a single server queue with general service time distribution, random system breakdowns, and Bernoulli schedule server vacations where, after a service completion, the server may decide to leave the system with probability  or to continue serving customers with probability 1 − .It is assumed that the customers arrive to the system in batches of variable size but are served one by one.If the system breaks down, it enters a repair process immediately.It is assumed that the repair time has general distribution, while the vacation time has exponential distribution.The purpose is to find the steady-state results in explicit and closed form in terms of the probabilitygenerating functions for the number of customers in the queue, the average number of customers, and the average waiting time in the queue.
Also Maraghi et al. [19] studied an M  /G/1 queue with random breakdowns and Bernoulli schedule server vacations where, after a service completion, the server may decide to leave the system with probability p or to continue serving customers with probability 1−.It is assumed that customers arrive to the system in batches of variable size but are served one by one.
Kumar and Sharma [20] analysed the concept of customer balking and reneging has been exploited to a great extent in the recent past by the queuing modelers.According to this model, a reneged customer can be convinced in many cases by employing certain convincing mechanism to stay in the queue for completion of his service.Thus, a reneged customer can be retained in the queuing system with some probability (say, ) and it may leave the queue without receiving service with probability  (= 1 − ).This process was referred to as customer retention.The effect of probability of retention on the average system size had been studied.
Recently in 2013, Baruah et al. [21] studied a two-stage queuing model where the server provides two stages of service one by one in succession.They considered reneging to occur when the server is unavailable during the system breakdown or vacation periods.
Kumar and Sharma [22] in their queueing model with retention of reneged customers and balking studied a finite capacity and the retention of reneged customers.The customer's impatience may occur in the form of balking also.In this paper, the work of Kumar and Sharma (2013) is extended by including balking to take into consideration the broader perspective of customer's impatience.The model is solved iteratively to obtain the steady-state probabilities of system size.Some useful performance measures are derived, and the effect of the probability of retaining the reneged customers on various performance measures is studied numerically.
1.1.Model Description.In today's competitive world, customer's impatience has become a serious problem.Firms are employing a number of strategies to retain their customers.To the best of our knowledge, studies on multistage batch arrival of service in reneging do not exist.We extend and develop the above models on two stages of service to multistage of services.In particular, a two-stage model of Monita Baruah, Kailash C. Madan, and Tillal Eldabi in the year 2013 motivated us for developing this new model which we consider here.Their model on two stages of service is extended to multistage of service where the arrivals occur in batches.Breakdown and reneging play a vital role in our model.Two stages of service will be benefited only for small scale industries, factories of low production, limited network, and so forth, whereas multistages of services render a huge amount of beneficial factors in large scale industries and in multitier applications in e-commerce in particular.
This results in an introduction of a new model: a non-Markovian multistage batch arrival queue with breakdown and reneging.In this model once the customers enter they will be provided with the service in all the stages of service department one by one in succession.Once the service gets completed, customer departs from the system.Once the multistage service of a unit is completed, the server is assumed to take vacation with probability  or may continue to offer service with probability (1 − ).As soon as the vacation period of the server ends, it joins the system to continue service of the waiting customers.During system breakdown (or) vacation periods, there will be unavailability of service; in that case reneging occurs.Reneging is assumed to follow exponential distribution.Arrival follows a Poisson process and service time follows general distribution.We study the changes in the probability of idle time and in the traffic intensity in some particular cases like (1) if there is no vacation and (2) if the service time and vacation time are exponentially distributed.Using supplementary variable technique we derive the steady-state probabilities and performance measures.
This paper is organized as follows.Applications of model is presented in Section 2, the mathematical description of the model is given in Section 3, the definitions and notations are described in Section 4, the equations governing the system are given in detail using birth and death process in Section 5, the queue size distribution at random epoch is highlighted in Section 6, and the average queue size and average waiting time are given in Section 7. Sections 6 and 7 gives a comprehensive idea about the model in which step by step method is followed for better understanding which helps in any extension of this model.The particular cases are derived in Section 8 in which the model is seen in various aspects or situations.For justification of the model, the numerical illustrations are given in Section 9, and, at last, the conclusion is presented in Section 10.In this section future scope of this model is discussed.

Applications of Model
With the development of computer and communication networks, queueing systems and networks have been identified as a powerful analysis and design tool for various applications.Our model plays a major role in multitier applications in e-commerce, motor production process, car production process, and textile industry.Some multistage real life applications in large scale manufacturing industry are described in the following sections which come across the necessity of our new model.

Three-Tier
Application in e-Commerce.Consider a threetier application, as shown in Figure 1.The first two tiers are clustered web server and clustered java application, respectively, and the third tier is a database used in ecommerce applications.Each tier is considered to be a stage.Each incoming application is subjected to admission control at the sentry to ensure that the contracted performance guarantees are met.The excess sessions are turned away during overloads.The load balancer efficiently distributes the incoming requests.Each input data is considered to be a customer.The server maintenance period is considered to be a vacation period, and this application can be extended to multistage of services.

Electric Motor Production Industry
Multistage in Motor Production Process (See Figure 2).In the electric motor manufacturing industries, multistage production process is implemented (Figure 2).The stator body is considered as the arrival customer.In the first stage, the stator body is wounded by the coil; in the second stage, the rotor part is fixed in the stator; in the third stage, the cooling fan is fixed in one end of the rotor; in the fourth stage, the end cover is fixed in the other end of the rotor; in the fifth stage, the entire motor is painted; and, finally, in the sixth stage, a complete electric motor is produced.We consider the stator body as a batch arrival customer.In the production sector, vacation denotes obtaining new raw material, maintenance of the tools, stock verification, and so forth.Here, the machine gets idle owing to breakdown, and the machine idle period is considered as the breakdown period.

Automobile Industry
Multistage in Car Production Process (See Figure 3).In automobile industry, car production process is considered as the multistage production process (Figure 3).Here, we consider the steel plate as the batch arrival customer.In the first stage of production, a roll of steel plate is cut according to the shape of the car parts in the second stage, and steel plates are welded; in the third stage, all the car parts are painted in the painting section; in the fourth stage, the engine assembly is mounted; in the fifth stage, extra fittings are assembled according to the customer requirements; in the sixth stage, quality control inspection is conducted; and, finally, a complete car is produced.We consider the steel plate as the batch arrival customer.In the production sector, vacation denotes obtaining new raw material, maintenance of the tools and machines, stock verification, and so forth.Here, machine gets idle owing to breakdown, and the machine idle period is considered as the breakdown period.

Textile Yarn Production Industry
Multistage in Production Process (See Figure 4).In the textile industry, yarn production is considered as the multistage production process (Figure 4).Here, a bunch of cotton is considered as the arrival customer.In the first stage, a bunch of cotton is inserted into the blow room machine; in the second stage, the refined cotton is inserted into the carding machine; in the third stage, the produced cotton sliver is inserted into the drawing machine and is again processed; in the fourth stage, the processed cotton sliver is inserted into the simplex machine; in the fifth stage, the refined cotton sliver is inserted into the spinning machine; and in the sixth stage, the cotton yarn output from spinning machine is inserted into the cone winding machine, and the entire cotton yarn is wounded as a cone.These six stages of production are required to convert a bunch of cotton into cotton yarn.In the production sector, vacation denotes obtaining a new raw material, maintenance of the tool and machine, stock verification, and so forth.Here, machine gets idle owing to breakdown, and the machine idle period is considered as the breakdown period.

The Model Is Based on the Following Assumptions. (a)
Customers arrive in batches following a compound Poisson process with a rate of arrival .Let    ( = 1, 2, 3, . . ., ) be the first-order probability that customers arrive at the system in batches of size  at the system in a short interval of time (,  + ), where 0 ≤   ≤ 1 and ∑ ∞ =1   = 1.(b) The server provides multistage of heterogeneous services one after the other in succession.An arrival batch receives the service offered at multistage in succession, defined as the first stage, second stage, and so forth, respectively.The service discipline is assumed to be on a first come first served basis (FCFS).Let us assume that the service time   ( = 1 to ) of the th stage service follows a general probability distribution with a distribution function   (  ), where   (  ) is the probability density function and (   ) is the th moment of the service time, with  = 1, 2, . . ., . Let , where  = 1, 2, . . ., . ( (c) Once the multistage service of a unit is completed, the server is assumed to take vacation with probability  or may continue to offer service with probability (1 − ).As soon as the vacation period of the server ends, it joins the system to continue the service of the waiting customers.
Let us assume the vacation time to be a random variable following general probability law with distribution  (V) and density function by  (V).(  ) denotes the th moment ( = 1, 2, . ..) of the vacation time.Here, let us consider that () is the conditional probability of a vacation period during the interval (,  + ), given that elapsed time is , which can be given as Thus, (d) Customers arriving for service may become impatient and renege after joining during vacations and breakdown periods.Reneging is assumed to follow exponential distribution with parameter .Thus,  () =  − ,  > 0. Thus  is the probability that a customer can renege during a short interval of time (,  + ).
The system may fail or be subjected to breakdown at random.The customer receiving service during breakdown returns back to the head of the queue.Let us assume that the interval between breakdowns occurs according to a Poisson process with a mean rate of breakdown  > 0.

Definitions and Notations
We assume that steady state exists and define  , () as the probability that there are n (≥1) customers in the system including one customer in type  service,  = 1, 2, . . ., , and elapsed service time is . , = ∫ ∞ 0  , ()  is the corresponding steady-state probability irrespective of elapsed time .  () is the probability that there are n (≥0) customers in the queue and server is on vacation and elapsed vacation time is .  = ∫ ∞ 0   ()  is the corresponding steady-state probability irrespective of elapsed vacation time .
Q is the steady-state probability of the server being idle as the server takes vacation.
The probability generating functions are defined as

Equations Governing the System
We frame the difference equations with Poisson input related to birth-death process.
Let  be the size of customers at time .Let   () be the probability of "n" customers in the system.Let (1)   be the average arrival rate when "" customers are in the system; (2)   the average service rate when "" customers are in the system; (3)   (+Δ) the probability of "" customers at (+Δ); (4)  the probability of breakdown arrival; (5)   the probability of batch arrival.
Note.LHS of (14) indicates that there are n customers in the system and the service is about to start in stage 1.RHS of ( 14) indicates the various situations in which the system has n customers and the service is about to start in stage 1.
The first term of RHS of ( 14) indicates that the system is in idle (), an arrival of batch consisting of  customers occurs, and 1st stage service is about to start.
The second term indicates that there are  + 1 customers in the system including one customer whose service is completed in Nth stage ( ≥ 2) with probability   () remaining  customers are there in the system.1st stage of service is about to start.Moreover the server is not going for vacation; instead stays in the system with probability 1 − .
The third term indicates the repair process gets completed with probability Υ ().So one customer just enters into the 1st stage service.Now there are  customers in the system.
The fourth term indicates that vacation gets completed with probability  ().There are  customers in the system.1st stage of service is about to start.
The combination of all the four terms gives the LHS of (14).
The integral in the RHS of (15) states that, after completion of 1st stage of service, there are  customers in the system to whom 2nd stage of service is about to start which gives the LHS of (15).
Similarly (16) indicates that the Nth stage of service is about to start.
RHS of ( 17) explains the fact that, after completion of vacation with probability , the server has just returned to the system consisting of n customers and the 1st stage of service is about to start.This leads to the LHS of (17).

Queue Size Distribution at Random Epoch
We use supplementary variable technique to get the following comprehensive idea about the model.

Average Queue Size and Average Waiting Time
As   = (/)   () | =1, with the mean queue size is of the 0/0 form, the L'Hopital's rule is applied twice to obtain where the primes and double primes denote the first and second derivatives, respectively: By utilizing (58) in (57), we can obtain the mean length of the queue size   while the mean waiting time of the queue   can be derived by using   =   /.Furthermore, we can also determine  =   + , the mean queue size of the system, and  = /, the mean waiting time of the system.

Particular Cases
Case 1 (no server vacations).In this case, the server has no vacation.Hence,  () = 0. Thus  = 0 in (45)-(47).Our For Nth stage, The probability of idle time is Furthermore by considering  = 0 in (58) we can obtain the mean queue size   and mean waiting time   .
Case 2 (exponential service time and exponential vacation time).In this case, we assume that the service time for the multistage of service with service rate   > 0,  = 1 to , is exponentially distributed.Furthermore, the repair time and vacation time are exponentially distributed with a repair rate Υ > 0 and vacation time  > 0. Thus, where ℎ 1 =  −  () + ; ℎ 2 =  −  () +  − (/).
By utilizing the abovementioned relations in (44)-( 46), we can obtain For Nth stage, where Therefore, the probability that the server is providing service in the first, second, and Nth stages at a random point of time can be given as presented in (67), (68), and (69), respectively.Consider the following: Furthermore, the probability that server is under repairs at random point of time is The probability that server is on vacation at random point of time is given by The probability that the server is idle but available in the system is given by Similarly the mean queue size and mean waiting time can be derived by determining   (1),   (1),   (1), and   (1) and utilizing them in (57).Consider the following: 5,6,7, and 8 clearly show that owing to the breakdown of the system and reneging, the proportion of the idle time of the server increases and the utilization factor decreases.Furthermore, Figures 9, 10, 11, and 12 demonstrate the effect of reneging and breakdown on the mean queue size   and the mean waiting time   and evidently indicate that, as the breakdown occurs as well as customers reneging from the queue, the average length of the queue decreases and the average waiting time decreases.

Conclusion
The present study clearly analyzed the multistage batch arrival queue with reneging during vacation and breakdown periods.The steady-state solutions are obtained by using supplementary variable technique and the mean queue length and mean waiting time are calculated.Furthermore, some particular cases are also discussed and numerical illustrations are presented.It can be concluded that, with the increase in the value of the parameter of customer's impatience reneging  for varying values of the breakdown parameter , the  utilization factor  decreases, whereas the probability of the server idle time  increases.In this model, due to reneging and breadown the mean queue size   and the mean waiting time   decteases and the average waiting time decreases.
For the future study compulsory vacation can be included in the present model and the effect of reneging during breakdown can be obtained.The model presented in this study can be utilized for communication networks and largescale manufacturing industries.

Figure 8 :
Figure 8: Effect of reneging  = 5, 9, 10, and 12 and breakdown at  = 4 on the proportion of idle time  and utilization factor .The figure gives the fact that as idle time of the server  increases, it leads to a decrease in utilization factor .

Figure 9 :
Figure 9: Effect of  = 5 and  = 1, 2, 3, and 4 on the mean queue size   and mean waiting time   .The figure explains that the average length of the queue   decreases owing to customers reneging from the queue, so the average waiting time   also decreases.

Figure 10 :Figure 11 :Figure 12 :
Figure 10: Effect of  = 9 and  = 1, 2, 3, and 4 on the mean queue size   and mean waiting time   .The graph indicates that as the average length of the queue   decreases, owing to customers reneging from the queue, the average waiting time   also decreases.

Table 1 :
Effect of reneging and breakdown.
Q Reneging