1. Introduction
Queueing inventory models have been extensively analyzed since 1992. Very few of these discuss multi-commodity systems in randomly changing environments. Queueing systems which evolve under influences from external sources have, for a long time, inspired interest. In real life situations, inventory systems are often subject to randomly changing exogenous environment conditions that affect the demand for the product, the supply, and the cost structure. The area of queues in random environments is today a field of active research in applied probability. Queueing systems with correlated arrival flow of customers give adequate mathematical models for different real world systems including computer and telecommunication systems, and network protocols [
1]. The following papers are relevant to the present paper only in that the authors consider multi-commodity inventory systems without any specified main and optional items. FaizAl-Khayyal et al. [
2] consider a multi-commodity network model in maritime routing and scheduling. They tried to optimize the cost and the quantity of each commodity with constrained production rates, consumption rates, and storage capacities in each port. This paper addresses the common problem faced in the maritime transportation of petrochemical products. The problem is formulated as a mixed-integer non-linear programming problem.
Two-echelon, multi-commodity supply chain network design problem is considered by Hannan et al. [
3]. The authors formulate the problem as a mixed-integer programming model in deterministic, single-period, and multi-commodity contexts. They develop a heuristic solution procedure based on Lagrangian relaxation due to the complexity of the problem. The problem deals with developing an optimum strategy in locating and sizing factories and warehouses and minimizing the total cost of the system which includes the costs associated with production, storage, transportation, and lead-times of commodities.
Jin QIN et al. [
4] discuss an optimization problem in multi-commodity logistic network design. The problem is formulated as a non-linear mixed-integer programming model based on assumed normally distributed stochastic demands of the retailers. In this work, the strategic decisions regarding inventory controls and facility locations are incorporated simultaneously. The authors also developed a combined simulated annealing (CSA) algorithm to solve the problem. Optimization of the cost function involving costs associated with location, inventory, and transportation is also done.
Ronald et al. [
5] examine a multi-commodity logistics network design problem with simultaneous emphasis on establishing ideal location facilities and distribution of commodities to ensure minimization of costs together with improved services. Factors such as the location of facilities and warehouses, storage capacity of the warehouses, and transportation routes involving all these are considered in the optimization problem. They later designed, tested, and compared a genetic algorithm and a specific problem heuristic on several realistic scenarios.
Claudio et al. [
6] discuss multi-commodity inventory location models where inventory control policies are reviewed periodically and continuously under modular stochastic capacity constraints. The model is formulated as a mixed-integer non-linear programming model. The logistic problem of supplying certain commodities to the warehouses, which acts as customer service centers, from a single factory is under consideration. An objective function with factors associated with the selection of warehouses, type, and quantity of commodities to be assigned, type of customers to be served, etc., is optimized. They have developed a Lagrangian heuristic to obtain a feasible integer solution at each iteration of the subgradient method.
Ali et al. [
7] examine a dynamic multi-commodity inventory and facility location problem in steel supply chain networks. Demand is assumed to be stochastic with normal distribution. In this model, the authors suggested a potential production capacity with emergency and shared safety stocks. The authors have presented a mixed-integer non-linear programming model and a mixed-integer linear programming model in this paper. The paper focuses on the strategic and tactical design of steel supply chain networks.
In Shajin et al. [
8], the authors discuss a single server multi-commodity queueing inventory system with one essential and m optional items. They were the first to introduce the concept of optional items for sales/service. Customer arrival follows the Markovian arrival process, service completion with respect to the essential inventory follows Phase-type distribution and that with respect to optional inventories follows an exponential distribution. In this model, immediately after the service of an essential item, the customer either leaves the system with probability
p or with probability
the customer goes for optional item(s). The system is assumed to be idle either in the absence of an essential item or when there is no customer in the system. Each customer is allowed to purchase only one unit of the essential item, whereas more than one type of optional items can be purchased with an imposed restriction of, at most, one unit per item. The stability condition is obtained by using the well-known fact that the left drift rate should be less than that of the right drift. Optimization of the control variables with respect to the cost function is also done numerically.
The following papers deal with queueing inventory systems influenced by randomly changing environments. In Song et al. [
9], the authors consider an inventory model where the rate of demand is dependent on the environment variables. These variables can be anything, such as different stages in the life cycle of the particular inventory or changes in various factors linked with the economy, etc. They not only derived basic characteristics of the optimal policies but also observed the influence of various patterns in problem data on optimal policies and developed algorithms for computing optimal policies
Özekici et al. [
10] describe inventory models with unreliable suppliers in randomly changing environments. The environment change follows a Markov chain. The dependence of the stock-flow equations of the system on random environments is represented by a two-dimensional stochastic process. Under specified conditions, they have derived an optimality condition for the base-stock policy and (
s,
S) policy. Computational issues and some extensions are also determined.
A single item inventory model which is observed periodically in a randomly changing environment is considered in Erdem et al. [
11]. All the model parameters are dependent on a time-homogenous Markov chain environment. The replenishment quantity is minimum{ Order quantity, Vendors capacity}. The problem is analyzed in single, multiple, and infinite periods. In all these cases, the authors prove that the optimal base-stock level depends on the state of the environment. Comparisons of the results with the case when the replenishment quantity equals the quantity ordered is also done.
Perry et al. [
12] discuss production-inventory models with an unreliable facility operating in a two-state random environment. The system is characterized by a production machine. The production can even be stopped purposefully when there is a limited stocking capacity. When the machine is in ON period, the input into the buffer is assumed to be continuous and uniform until the threshold is reached, whereas the output from the buffer follows a compound Poisson process during OFF periods. Two different models are discussed and the factors controlling OFF periods are determined.
A continuous review (
s,
S) inventory system in a randomly changing environment is discussed in Feldman et al. [
13] and its steady-state distribution obtained. The demand process is an environment dependent compound Poisson process when the environment is in a fixed state during an interval of time. The environmental process follows a continuous-time Markov process.
Kalpakam et al. [
14] consider a lost sales (
s,
S) inventory system in a random environment. No backlog is allowed. The demand and supply rates are influenced by the environment process which is a finite irreducible Markov chain in continuous time. They have obtained the transform solution of the inventory level distribution and also an efficient algorithm to evaluate the long run system state is provided. Moreover, transient and limiting values of the mean reorder and shortage rates are also obtained. Goh et al. [
15] discuss price-dependent inventory models with discount offers at random times. The offer is accepted when the inventory position is lower than a threshold level. Three different pricing policies are considered in which demand is induced by the retailer’s price variation. They have obtained expressions for optimal order quantities, prices, and profits under the assumptions of constant demand rates.
Highlights of this paper are:
It considers multi-commodity inventory with positive service time [
16] in finite number of randomly changing environments;
The first paper to introduce optional items for service in random environments;
Except for one item (essential), all others are optional;
The customer demand process follows Markovian arrival process (MAP);
The environment change process follows marked Markovian environment arrival process (MMEAP)[n] of order n;
Service time of customers, being served with the essential inventory follows phase type distribution and that w.r.t optional item(s) follows exponential distribution (depending on the environment). The latter has a parameter, depending on the specific item(s) demanded by the customer.
The rest of the paper is organized as follows. The mathematical formulation of the model including the stability condition and the steady-state probability vector is described in detail in
Section 2.
Section 3 deals with some system performance measures and in
Section 4, the construction of the cost function for optimizing the system control variables is discussed. Numerical illustrations and the numerical analysis of the cost function are discussed in
Section 5.
Section 6 gives the conclusion followed by references.
Notations and abbreviations used:
ordering policy: An inventory policy which says that when the inventory level falls below a certain minimum number s, the order for replenishment is made to restore the inventory to a maximum number S;
e = Column vector of appropriate order with all its entries as s;
= Matrix of appropriate order with all its entries as 0;
= Identity matrix of order k;
= element of the matrix G;
Continuous time Markov chain;
Level Independent Quasi-Birth and Death process;
= Markovian arrival process;
= Marked Markovian environment arrival process (with n distinct environments);
= The Kronecker product of two given matrices and , given by of order ;
= The Kronecker sum of two square matrices B and C of orders m and n, respectively, given by ;
= The Correlation coefficient;
= Markovian arrival process with positive ;
= Markovian arrival process with negative ;
: Phase type distribution with the initial probability vector and the transition generator matrix T.
2. Mathematical Formulation
Consider a single server multi-commodity queueing inventory system with one essential and
m optional inventories in
n random environments. Only one environment will be in operation at any given time. The arrival of customers follows Markovian arrival process (
) with representation (
,
), where each
for
is of order
. The generator matrix of the underlying
(
) on the state space {
} is given by
. These matrices,
and
, are of the form
where,
for
. Thus,
gives the transition rate from
state to
state through an arrival, while
, gives the transition from
state to
state without an arrival. Note that the transition rate between the
states, given by
occurs only with an arrival. Let
be the steady-state probability vector of
H. Then,
satisfy
= 0 and
= 1. The fundamental rate
of this
is given by
=
which gives the expected number of arrivals per unit of time. The coefficient of variation
of intervals between arrivals is calculated as
and coefficient of correlation
of intervals between successive arrivals is given as
.
There are n environments that occur randomly and the occurrence of the environments follows marked Markovian environment arrival process () with representation (,,,⋯,), where each for is of order . As stated earlier, only one environment will be in operation at any given time. The change in environment is directed by the stochastic process which is an irreducible continuous time Markov chain with the state space . The sojourn time of this chain in the state is exponentially distributed with parameter . When the sojourn time in the state v expires, the process jumps to the state without any change in the environment with probability where . On the other hand, the process jumps to the state with the arrival of environment with probability where and .
The behavior of the
is completely characterized by the matrices
defined by
The matrix
represents the generator of the process
.
Service time of those customers who are served with the essential item is phase type distributed with representation of order . This service time is the time until the undergoing Markov chain () with a finite state space {} reaches the absorbing state . = () gives the initial probability of starting in any of the states. T is the generator matrix that gives transition rates within the states {}. The absorption rates from the individual transient states {} to the absorption state is given by . = gives the mean service of the customer.
Service time of those customers who are served with the optional items are environment dependent and they are exponentially distributed with parameter , where and i ∈ {}. It is important to note that no order preference has been given to any element, i.e., , and so on where each ∈ {} with .
In this model, a customer is allowed to demand exactly one unit of the essential inventory where as more than one type of optional inventory can be demanded by a customer with an imposed restriction of, at most, one item from each optional inventories. Service rates of the optional items are assumed to be environment dependent. The optional item is served in the environment with probability , similarly the and optional inventories are served in the environment with probability , and so on. If the demanded optional inventory is not available, the customer is expected to quit the system after acquiring the essential item together with those available optional inventories. The server is assumed to be in the idle state in the absence of customers, as well as essential inventories. Essential and optional inventories have exponentially distributed positive lead time with parameters and for , respectively. The essential inventories are under the control policy whereas the environment dependent optional inventories are under the (,) for and control policies in the environment.
At any given time
t, let
,
,
,
,
,
, and
denote, respectively, number of customers in the system, status of environment, number of essential inventory items, number of
optional inventory items for
, service phase, environment phase and arrival phase of the customers. The status of the server at any given time
t is defined as,
Let be the collection of all the permitted combinations of different optional inventories and let denotes the server status, in general, for the combined service of u optional items, for . Thus, the process = is a which is a level independent quasi-birth and death process with state space as follows
for
for
for
for if where then for and
for if where then for and for
for if where then for and for for .
The infinitesimal generator
of the system is of the form
Matrices and are of order and , respectively, their entries are due to the arrival of customers following MAP with representation . Matrices and are of order and , respectively, their entries are due to the service of essential inventories following phase type distribution with representation and also due to the environment dependent, exponentially distributed service of optional inventories. Matrices and are square matrices of order a and c, respectively, their entries includes the replenishment rates of the inventories in addition to the negative sign of sum of other entries of the same row found in , , , and , where , and .
When then .
When then for
When then for
When then for and so on.
When then
In order to have a better understanding of the system, a detailed illustration of the model has been provided in
Appendix A by fixing the number of optional items
and the number of environments
. All the transitions and resultant component matrices are shown clearly in the
Appendix A.
2.1. Stability Condition
Let be the steady-state probability vector of where for
Refer to
Appendix A for the component matrix representations From (1),
where,
Solving Equations (2)–(5), we get
where,
The only unknown probability vector
is obtained from the normalizing condition
Theorem 2.1. The necessary and sufficient condition for the stability of queuing inventory system under study is Proof. The queueing system with the generator
under study is stable if, and only if,
Refer to
Appendix A for the component matrix representations.
Using Equations (6)–(9) together with the matrices
and
we get
where
for
where
for
Let
Then, by (10) we get the stated result. □
2.2. Steady State Probability Vector
Let
denote the steady state probability vector of the generator
. Then, we have
Partitioning
as
, from (23) we get
By assuming the stability condition, we see that
is obtained as (see [
17])
where
R is the minimal non-negative solution of the matrix quadratic equation
The boundary conditions are given by
From Equation (23) we get,
and by the normalizing condition in
, we get
where