1. Introduction
Periodic review models are employed in many logistical systems in which continuous review is restrictively expensive or infeasible for lack of necessary resources [
1]. When such review is conducted and some items are found to be nonfunctional, depending on the system, replacement parts may be ordered or the faulty units themselves may be repaired. In the former—the order model—the supplier of these parts typically has these items readily available and therefore once an order is issued for a batch of items, the entire batch is shipped to the ordering facility. It is generally assumed that order shipping times are i.i.d. and therefore each order may take a different time to arrive and possibly, orders may even crossover. However, within an order, all the items will arrive at the same time. In the latter setting—the repair model—items must be repaired at a repair facility with the standard assumption that item repair times are i.i.d. The key difference between the two models is what element forms the basic stochastic variable of the model—the order (as in the order model) or the item (as in the repair model).
Repair models are an important component of sustainable logistics, hence their importance. Examples abound in the practice of operations and logistics. Closed-loop supply chains encourage repair to recycle materials and resources used in the manufacturing process (e.g., [
2]). Other examples include hybrid manufacturing–remanufacturing and green products remanufacturing [
3,
4].
This paper analyses a repair model under two inherently different settings, in-house repair (IR model) and outsourced repair (OR model). In the first, items are repaired on-site and therefore once each item is repaired, it can immediately be returned to stock. In contrast, when repair is outsourced, the repair facility may be reluctant to ship each item individually, and therefore only once all the items that compose a single order are repaired, then the order is shipped to the ordering facility.
Many inventory models use the fill rate to measure their performance. In many settings of practical interest, however, the fill rate is not a reasonable approximation of the customer–inventory-system relationship. The reason for this is that the fill rate tacitly assumes that the reputation cost to the system due to customer waiting begins immediately when the customer arrives. However, it is often that firms are obliged, whether by government regulation or contractual commitment, to reduce the waiting time to be below a predetermined time threshold. From the customers’ standpoint, too, there is usually a certain tolerable or acceptable period of waiting, which may depend on their level of patience or expectation. In these cases, the time from which the firm suffers reputation costs begins only after this tolerable wait has passed [
5]. Accordingly, the performance measure we use for our facility is the window fill rate, which is the probability that each item is replaced within a predetermined window of time.
One obvious way to improve performance is keeping a large stock of spare parts. Since this may be a very expensive solution, it is of special interest to managers to determine how spares affect the window fill rate, thereby keeping the minimum inventory to meet their target window fill rate. We use our analysis of the window fill rate of both models to examine the cost of outsourcing in terms of additional inventory. This is performed by comparing between the number of spares needed in each model (i.e., the IR and OR models) to meet a required window fill rate level. Prima facie, the OR model requires more spares than the IR model. This is because, all else being equal, the OR model requires waiting for all the items in an order to be repaired before they are returned to the inventory, whereas in the IR models, each item returns to stock as soon as it is repaired.
We begin by deriving the stationary window fill rate for the IR model by tracking the demand and supply for spares and computing the probability that within
w units of time from a random customer’s arrival there will be a functional item available for the customer. Since this entails finding the difference between supply and demand, each of which are Poisson random variables, our formula makes repeated use of Skellam random variables [
6], whose distribution is nowadays ubiquitously available.
Next, we derive the stationary window fill rate for the OR model. Here, however, there are two major complications. First, in contrast to the IR model, the delivery times of repaired items in the OR model do not follow a Poisson process, since the delivery times of items in a particular order are dependent on each other. To overcome this, we track each order’s status, that is, whether it was delivered or not, and for each state in the delivery state space, we derive the probability of the state and the expected window fill rate. The weighted average of these expected window fill rates is the system’s window fill rate. The second complication in the OR model is with respect to the evaluation of the window fill rate’s formula, which requires conditioning on the values of numerous random variables and therefore cannot be executed within a reasonable time and accuracy. This is overcome by identifying the components of the formula whose evaluation can be expedited through simulation. This hybrid approach of using simulations and analytical evaluation allows us to evaluate the window fill rate accurately (within less than 1% margin of error) within reasonable time.
We complete our mathematical analysis with a numerical illustration that illustrates functional form of the window fill rate. In both models, the window fill rate is S-shaped in the number of spares and the implication of this to managers is discussed. The difference between the graphs of the IR and OR model reveals that all else being equal, a shift from internal repair to external repair results in a lower window fill rate. Thus, to maintain the same level of performance, more spares must be procured when outsourcing repair compared to in-house repair. We estimate this cost in our numerical example for different levels of tolerable waiting and target performance.
Our study, therefore, ignores managerial and operational costs of in-house and outsourced repair. See, in contrast, [
7] that optimizes the stock management and ignores the inventory cost and [
8] that considers maintenance and inventory jointly. Therefore, the practical decision of whether to repair in-house or outsource repair should weigh possible managerial advantages of outsourcing against the cost of additional spares needed to maintain the required system performance.
The contributions of our model are therefore both theoretical and applicable in their nature. We contribute to theory by developing a formula for the window fill rate in the IR model and an efficient algorithm for the estimation of the window fill rate in the OR model. Our contribution to the practice of inventory modeling is threefold. We provide a tool for managers to quantify how many spares are needed to meet their performance objectives. Second, it allows managers to assess the cost—in terms of spares—of changing contract terms such as service times and performance guarantees. Third, it allows managers to be informed about inventory requirements when deciding whether to outsource their repair operations.
The remainder of this paper is as follows: In the next section, we review the literature focusing on inventory systems and performance measures related to our study. In
Section 3 and
Section 4, we develop the formula and algorithm for the window fill rate in the IR and OR model, respectively.
Section 5 is a numerical illustration that demonstrates the cost of outsourcing in terms of spare capacity and
Section 6 concludes and details future paths for research.
4. OR Model: Outsourced Repair
In this model, instead of repairing the failed items in-house, the warehouse ships the failed items that arrived in the recent cycle to an external repair contractor. Since the repair is outsourced, the warehouse must wait for the contractor to ship back the repaired items. However, a contractor will be reluctant to ship each item separately, and will ship all the items of the orders together. Thus, although each item repair time is independent with cumulative distribution , their delivery does not follow the same distribution. Furthermore, we assume that there are a number of such contractors so that it is possible that orders may crossover. Finally, to simplify the comparison between the IR and OR models, we assume here that shipping times are deterministic. Reducing the shipping time from the window fill rate is equivalent to zero shipping times and therefore, without loss of generality, we assume shipping times to be zero.
The OR model does not differ from the in-house repair model with respect to the demand, and therefore (
3) applies in this model, too. Since by (
2),
, then (
3) can be written as follows:
In contrast to the demand, there is a marked difference between the two models with respect to the supply. The key difference stems from the different definitions of what constitutes a delivery. In the IR model, deliveries are of single repaired items and this implies that the arrivals of items delivered on time (i.e.,
) follows a Poisson process. In contrast, in the OR model, deliveries are of repaired orders. As a result, if the order was delivered, then
and
, and if not, then
and
. Consequently,
cannot be formulated as a Poisson process and the analytical techniques used in the IR model fail. To continue the analysis, therefore, we must track the delivery state of each order. To that end, we define the stochastic binary variable
as follows:
Similarly to the analysis in the IR model, the supply comprises the
S spares in the system plus items that were brought by customers and were repaired until
. if Jane arrived early in the cycle, then the orders relevant to the supply are all the orders issued until date
, whereas if Jane arrived late in the cycle, then the order issued at
must also be considered. Notice that when the
’th order is part of the supply, then the (one) item brought by Jane must be considered with it. Thus, the supply is given by
where, recall,
when Jane’s arrival is late in the cycle and
when it is early in the cycle. The net
demand (
) is given by the demand minus the supply as follows:
In the first row above, the
’th order includes the additional item brought by Jane. In the second row, we use the notation
and the fact that
. The non-stationary window fill rate is the probability that the net demand is non-positive as follows:
4.1. The Probability for Delivery on Time
To continue the derivation, we begin with determining the probability that order
j is delivered on time, i.e., that all the items that compose it have been repaired by
. The probability that an arbitrary item from order
j was repaired by
is
The probability that the
order is repaired (i.e., delivered) on time is itself a random variable,
, where, recall,
. A notable exception is the order containing Jane’s arrival for which
since we must also include her item. Thus, when
(i.e.,
j is not the period in which Jane arrived),
When
, then we must also ensure that Jane’s item has also been repaired. Therefore,
4.2. Evaluating the Window Fill Rate
To evaluate (
24), it is necessary to track the delivery status of each order. In the literature review, we describe two ways how this was performed in similar models. One way is to assume that orders do not crossover—an assumption that effectively contradicts the independence of the lead times—and therefore is only an approximation. The other approach is to account for each combination of arrival status of the orders. In both ways, however, the lead times of the orders was given exogenously and therefore the probability for order delivery was independent of the order size and depended only on the time between the order’s issuance and
. In contrast, our model attempts to mimic outsourced repair in a more realistic fashion and therefore order repair times depend on the order size.
One may be tempted to use the window fill rate formula of these models and use the expected probabilities for repair
as the probability that order
j was repaired on time. However, this approach generally results in a very poor approximation since the real probability for repair is strongly dependent on the order size. To obtain an accurate evaluation, one must condition on the value of the order size of each of the orders that comprise (
24). Unfortunately, such an evaluation is extremely cumbersome and accuracy must be severely compromised if it is to be executed in reasonable computer running time. In what follows, we overcome this challenge by integrating simulation into the evaluation of the window fill rate.
4.3. Simulation
Simulation is commonly used to solve logistic models [
44,
45]. One of the main drawbacks of simulation is running time, a problem that is further exacerbated when the system is complex or when the parameters of the system change repeatedly and the simulation must be frequently re-executed. Therefore, in our proposition, we use a minimalist approach in which the simulation is incorporated into the evaluation narrowly, just enough to overcome the mathematical complexity of the problem, and the remainder of the evaluation makes use of (
24), the window fill rate formula (see a similar approach in [
16]). Other methods to overcome these accuracy and running time tradeoffs include constructing artificial neural networks [
46], and employing various heuristics such as interior point optimization, genetic searching, and particle swarm optimization [
47].
Recall, it is assumed that
has a finite support, and therefore there exists a
K such that for all
,
. Consequently, we can disregard dated orders that have been surely supplied since for these orders, the supply cancels with the demand (i.e.,
). The range of orders that are relevant to the analysis depend on whether
t is early or late in the cycle, denoted by
a, defined in (
1). Specifically, if
t is early in the cycle (
), then we need to consider the orders between
and
, whereas if
t is late in the cycle, the relevant orders are the orders between
and
(see
Figure 2). Thus, we can rewrite
, (
23), as the following:
We note that if the tolerable wait is sufficiently long such that
, then the window fill rate is 1 since order
was delivered on time, implying supply is certainly greater than the demand. We therefore proceed under the assumption that
. In (
28), there are at most
Poisson random variables that need to be simulated. These are,
. Their distributions are
,
and the other
,
are distributed
.
Similarly to the IR model, the steady state window fill rate satisfies formula (
13); namely, we need to remove the dependency on
i and take the mean along a cycle. Note that as long as
i is sufficiently large (i.e.,
), then (
28) is in fact independent of
i, since at most
K periods prior to Jane’s arrival are relevant to the analysis and all the probabilities are independent of the value of
i. It remains, therefore, to take the mean of the non-stationary window fill rate along the arrival cycle. This is performed in the algorithm described below.
4.3.1. Re-Indexing the Orders
To describe the simulation algorithm, it is more convenient to re-index the periods such that the arrival cycle, over which we will take the mean, is indexed as cycle zero. That is, order
is now relabeled as order zero. See
Figure 2 in which we plot the relevant orders (i.e., orders that are not surely cancelled out) re-indexed. Thus, with some abuse of notation (we allow the Poisson random variables to take negative times), we can rewrite
as follows:
4.3.2. Main Simulation Algorithm
We define
T as the number of segments used to evaluate the integral (
13). Since each cycle is of length
r, we have that
. The arrival cycle itself, i.e.,
, is simulated as a sum of
T independent random variables, denoted by
, where
represents the demand along
and therefore each
is distributed
. Consequently, the demand components of Jane’s arrival cycle can be stated as follows:
In summary, the window fill rate for a customer that arrived
t units of time into a cycle is
where
The following algorithm computes the stationary window fill rate,
. It simulates
M simulation threads and computes the window fill rate for each thread. The output of the algorithm is the expected window fill rate. Each simulation thread comprises the realizations
of the periodic demands
, respectively.
is unique in that it is constructed by its subsegments. That is, the realization
is the sum of
, each of which are the realizations of
, respectively. For each simulation thread and arrival time
t, the integral of
is computed.
Algorithm 1: Simulation Algorithm |
|
4.3.3. Evaluating the Non-Stationary Window Fill Rate
It remains to describe how
is defined. To this end, we note that for a given demand realization and the arrival time
t, the distributions of
(
29) are known and are given by
The
algorithm receives a single demand realization and arrival time in the cycle,
t, and returns the value of (
29), i.e.,
. It considers all the possible states of the demand (i.e., the arrival status of each of the relevant orders). Since there are
relevant orders, the state space is
, represented here by the numbers
. The binary representation of
k dictates the status of each of the orders. For example, if
, and
, then
. When
t is early in the cycle, the
firstK (of
relevant) periods are considered since the latest one is irrelevant because it neither belongs to the supply nor the demand. In contrast, when
t is late in the cycle, the
lastK periods are considered since the earliest period is irrelevant since it has surely arrived and therefore the supply and demand of that period cancel each other.
For each demand state
k, we compute the probability of this state (
) using (
26) and (
27) and the net demand (
) using the left side of the inequality in (
29). The sum of the probabilities whose state meets the condition
is the value of (
29) and is the output of the algorithm.
Algorithm 2: Algorithm |
|
5. The Capacity Cost of Outsourcing Repair
In this section, we demonstrate the functional form of the window fill rate for both models and use these graphs to determine the number of spares needed to meet the target performance level for each model. In what follows, we use the following baseline values and change at most one parameter at a time to highlight its effect. The order cycle time days, the tolerable wait days, the demand arrivals rate items per day, and the item repair times are distributed uniformly days.
In
Table 1, we describe the result of the baseline case for different levels of spares for the IR and OR model. For the OR model, a single simulation comprised
simulation threads. To test the accuracy of the simulation, for each value of spares, the simulation was repeated 30 times. We reported the mean and standard deviation. The largest standard error is 0.8% and is obtained at
for which the window fill rate is 72.2%.
In
Figure 3 and
Figure 4, we show the window fill rate and the IR and OR models, respectively. In each figure, we depict the window fill rate for different order cycle times. The window fill rate is generally S-shaped, initially convex and then concave. This feature is exploited by [
5] to develop a fast optimization algorithm to allocate spares in a multiple-location setting, and can be employed in both the IR and OR models, too.
Figure 3 and
Figure 4 also highlight the profound effect of the periodic review. While shorter review periods are costly in terms of order costs, they result in a higher window fill rate and therefore a need for less spare items. This, too, has been demonstrated in other models such as [
37] and should be always considered by managers when they consider reducing order frequency to ordering costs. For example, if managers wish to maintain a performance level of an 80% window fill rate, then by
Figure 3, in the IR model, if
, then only 9 spares are needed whereas if
and
, then 14 and 17 spares are needed, respectively. For the OR model, by
Figure 4, the numbers of spares needed are 17, 22, and 27 when
, and 10, respectively.
A comparison between the graph of the OR and IR models allows us to determine the cost of outsourcing repairs (compared to in-house repair) in terms of spares needed to maintain a required level of performance. This cost can be deduced from
Figure 5 that plots the window fill rates of the IR and OR models against the number of spares when the tolerable wait is
(left panel) and
(right panel).
In
Table 2, we list the numbers of spares that are required to meet common industry performance target levels;
, and
, for different values of tolerable waiting. The column in bold in the table is the difference between the two models and represents the cost that we are seeking to determine. For example, when
and the required performance level is
, then the OR model dictates that 22 spares must be procured whereas the IR model requires only 14. Therefore, if managers consider outsourcing repair, they must also consider that they will need to augment their stock with eight more spares.
Managerial Implications
Our numerical example illustrates a number of implications that are useful for managers. First, the functional form of the window fill rate in both models is generally S-shaped with the number of spares. That is, when spares are sufficiently low or sufficiently high, adding one more spare offers only a slight improvement. In the former, this is because the system is so short on spares that one more item still leaves it lagging behind. In the latter, the system is so well-off that there is hardly a need for improvement. It is therefore only when the system is somewhere in-between that any additional spare makes an impressional change since it enables the system to serve a large portion of customers on time, something it would not be able to conduct otherwise. Therefore, managers should put checks on the impulse to procure additional spares whenever possible. This should be performed only if the effect of the additional spares will indeed justify their cost, something that heavily depends on the current number of spares in the system.
Second, our sensitivity analysis demonstrates that the window fill rate shifts to the right when the order cycle time increases (
Figure 3 and
Figure 4) and when the tolerable wait decreases (
Figure 5). In other words, if managers have a shortage of spare parts, they must renegotiate contracts for longer service times. Alternatively, managers can offer contracts with shorter service times at a higher premium. The necessary changes to the contract (service times, contract price) can be estimated using our model.
One more implication is with respect to the decision to outsoure repair. As shown in
Figure 5, the decision to outsource repairs results in a lower window fill rate. To compensate for this, more spares must be purchased. The cost—in terms of spares—for the baseline case of this decision is demonstrated in
Table 2. The most important managerial lesson, therefore, is that while outsourcing has many advantages, managers must keep in mind costs associated with inventory requirements. Our model serves as a tool to make an informed decision about these costs. Furthermore, our baseline example shows that these costs depend very much on the model variables. For example, in the baseline case, the cost of outsourcing is increasing with the tolerable wait and the target performance level.
6. Conclusions
In this paper, we develop the window fill rate formula when item repair is conducted in-house and when repair is outsourced. The difference between the two models is that with in-house repair, items are repaired and returned independently of the other items in the order and therefore an item that was “lucky" to be repaired quickly will return sooner to the system’s stock whilst the other items in its order are still being repaired. In contrast, when repair is outsourced, then it is assumed that when an order for repair is shipped to the external repair firm, then no item will return separately from the other items in the order. Therefore, the system’s stock depends on the longest repair time within each order. Consequently, all else being equal, outsourced repair will result in a lower window fill rate compared to in-house repair.
The derivation of the window fill rate in both models requires tracking the supply of items when the customer demands an item (i.e., w units of time after the customer’s arrival) and the demand of items up to and including the customer’s arrival. In our Poisson demand model, the derivation of the in-house repair formula is quite straightforward since the independence of the delivery times of repaired items implies that these deliveries behave as a Poisson process, too. Therefore, the window fill rate in the in-house repair case can be reduced to a difference between two Poisson processes (i.e., a Skellam random variable) and can be easily computed. In contrast, the derivation of the window fill rate under outsourcing is computationally challenging since items of the same order depend on each other for their delivery.
We overcome this challenge by developing a Monte-Carlo simulation algorithm that considers all the possible delivery states of orders that are relevant to the supply–demand equation. The efficiency of the simulation algorithm stems from the fact that we do not simulate the entire logistical process but instead simulate only the demand arrivals.
We complete our analysis with a numerical illustration that demonstrates the functional form of the window fill rate as a function of the number of spares. The cost of outsourcing repair (in terms of required spares) is discerned by comparing the window fill rate graphs of both models. A manager of a logistic system that considers outsourcing repair must therefore consider how many more spares they will need to procure to maintain the same window fill rate level.
While our model assumes a Poisson demand, our mathematical approach can be used in a similar manner when arrivals follow a compound Poisson process. Unfortunately, for demand that does not follow a Poisson process, deriving a formula for the window fill rate may prove to be too mathematically challenging. Notwithstanding, this study serves even these cases by drawing managers’ attention to the need of adjusting stock levels when conducting organizational changes such as outsourcing repair or establishing in-house repair.
Our model may be extended in a number of ways. As previously stated, our model requires Poisson arrivals and therefore additional research is needed to develop efficient algorithms for estimating the window fill rate in non-Poisson arrival cases. Furthermore, we focus on a single warehouse, whereas large inventory systems employ a network of warehouses with central control of inventory. These systems may be organized in multiple echelons and may provide emergency replenishment from nearby sites. Future work can build on our approach to analyze these complex systems and estimate the optimal number of spares and their allocation.