Next Article in Journal
Third-Party Reverse Logistics Selection: A Literature Review
Previous Article in Journal
Designing an Intelligent Scoring System for Crediting Manufacturers and Importers of Goods in Industry 4.0
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stock Levels and Repair Sourcing in a Periodic Review Exchangeable Item Repair System

Department of Industrial Engineering, Jerusalem College of Tehnology, Jerusalem 9116001, Israel
Logistics 2024, 8(2), 34; https://doi.org/10.3390/logistics8020034
Submission received: 1 January 2024 / Revised: 27 February 2024 / Accepted: 4 March 2024 / Published: 22 March 2024
(This article belongs to the Section Artificial Intelligence, Logistics Analytics, and Automation)

Abstract

:
Background: Exchangeable item repair systems are inventory systems. A nonfunctional item is exchanged for a functional item and returns to the system after being repaired. In our periodic review setting, repair is performed either in-house or outsourced. When repair is in-house, a repaired item is returned to stock regardless of the repair status of the other items in its order. In contrast, with outsourced repair, the entire order must be repaired for it to return to stock. Methods: We develop formulas for the window fill rate (probability for a customer to be served within a given time window) to measure the system’s performance and compute it for each repair model. The cost of outsourcing is the difference between the number of spares needed to maintain a target performance level when repair is internal and when it is outsourced. Results and Conclusions: In our numerical example, we show that the window fill rate in both models is S-shaped in the number of spares and show how the graph shifts to the right when customer tolerance decreases and order cycle time increases. Further, we show that the cost of outsourcing is increasing with customer tolerance and with the target performance level.

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.

2. Literature Review

2.1. Exchangeable-Item Inventory Systems

The subject of our study lies within the research field of exchangeable-item (also called recoverable-item, repairable-item) inventory management. This is a vast field of research that is described in many professional books and is a main topic in top operations research, management, and logistics journals [9,10,11]. In the classical model, customers arrive at the inventory system (i.e., a logistics warehouse) with a failed item and exchange it for a serviceable item from the system’s stock. The failed item itself is repaired on- or off-site, after which it is returned to the system’s stock. Since repair is time-consuming, spare items are placed in stock in order to increase the number of items available for arriving customers. We follow two assumptions followed by most studies in the field. The first is that customers are served according to a first-come–first-serve policy (FCFS). Models that pivot away from the FCFS service rules assume customer priority (e.g., [12,13]). The second standard assumption is that there is ample repair capacity. This then allows us to assume that repair times are i.i.d. (see, in contrast, ref. [14], for a limited capacity model).
The review policy in a single location can be grouped into two main approaches. The continuous review—typically an ( S , S 1 ) model—is one of the assumptions of Sherbrooke’s seminal paper [15] and more recent studies such as [5,16]. The second group—which this paper belongs to—comprises models with periodic review policies. We now focus our attention to studies belonging to this group.
Periodic review policies may be further grouped to fixed order size policies and the “up-to" policies. In the former, notated by ( t , s , q ) , inventory is reviewed every t periods, and if it falls below s, then an order of q (or its multiples) is issued. Our model belongs to the latter group in which order sizes are designed to raise the current inventory level up to a target level. This group can be further split into ( r , s , S ) and ( r , S ) policies. With ( r , s , S ) , inventory is reviewed every r periods, and if the order falls below the threshold s, then an order is issued so that the inventory level will reach S (e.g., [17,18]). Our model assumes an ( r , S ) policy such that the inventory is reviewed every r periods and upon review, an order is issued so that the inventory level will reach S. Recent studies based on this approach include [19,20,21].
To complete this review, we note that certain operations may be carried out continuously while others are executed periodically. Recent examples include [22], a continuous ( s , S ) model whose parameters are updated periodically, and [23] that analyzes a network of remote locations in which the inventory review is continuous, but order delivery is constrained by periodical shipments.

2.2. Performance Measures

Among a logistic system’s performance measures, the expected waiting time and the fill rate are probably the most prominent. Our research is related to the second main performance measure, the fill rate, which is defined as the portion of demand that is satisfied upon arrival [24]. In contrast to the expected waiting time, the fill rate does not penalize the system for the length of the customer’s weight but on whether the customer had to wait. See [25] [Section 2.1] for a brief review of fill-rate-associated performance measures.
The notion of a tolerable wait, also termed as “reasonable wait” or “acceptable wait” [26], was extensively researched in service-related fields such as medical service [27], transit [28], supplier selection [29], and call centers [30]. In contrast, within inventory sciences, there is surprisingly little research that considers such willingness to wait. This lacuna is clearly unjustified since many inventory contracts actually do state time windows for demand to be supplied. As there are no losses to the inventory systems as long as they meet these time windows, these contracts effectively imply that customers have a tolerable wait. In a similar vein, ref. [31] (p. 744) complained about the failure of inventory models “to capture the time-based aspects of service agreements as they are actually written”.
The term “window fill rate” denotes the probability to wait less than a predetermined window of time. Mathematically, the window fill rate can be derived by assigning the tolerable wait in the waiting time distribution. Ref. [5] characterized the functional form of the window fill rate and used this to develop an efficient spares allocation algorithm that optimizes the window fill rate of a multiple-location inventory system. In a two-echelon system, optimizing the window fill rate requires varying degrees of resource pooling [32] and this was later applied to the battery swapping problem [33,34].
In a periodic review setting—the subject of this paper—there have been very few studies that consider the window fill rate [35]. In the ( r , S ) model, the assumption that repair times are independent stochastic poses a considerable mathematical challenge since to determine accurately when a customer is expected to leave the system, the delivery status of all past orders must be considered. The literature offers solutions to this problem. The first is to assume that orders do not crossover [36,37]. This assumption effectively contradicts the independence of the lead times and [38] shows, using real data, that this results in a poor approximation. The second approach is to overcome the challenge of the evaluation of the exact formula by using simulation for key random variables of the window fill rate [38]. In both studies, however, the lead times of the orders was given exogenously and therefore the probability for order delivery was independent of the order size. In contrast, our model mimics outsourced repair in a more realistic fashion by assuming that item repair times are independent and therefore order repair times do depend on the order size.

2.3. Outsourcing Repair

The decision whether to outsource repair or not depends on multiple considerations, including financial, environmental, informational, and logistical [29,39]. Relevant papers considering logistical aspects of this decision include [40] whose model outsources the repair of non-conforming items. Similarly, ref. [41] set repair to be performed either externally or internally depending on whether the “internal repair shop has the authorization (from the original equipment manufacturer) and the capability to repair the part (p. 412).” Finally, ref. [42] characterized outsourced repair as one that does not require the system to acquire repair infrastructure and therefore outsourced repair costs are only variable costs.
The closest study to our proposal with respect to repair outsourcing is [43] that compared between in-house and outsourced recycling (comparable to the concept of repair in our proposal) in a continuous ( r , Q ) model (i.e., when stock level reaches r order quantity Q). In addition to their model being continuous whereas the system in our proposal is under periodic review, their objective was minimizing costs while our paper focuses on the window fill rate as the system’s objective.

3. IR Model: In-House Repair

3.1. Model Setup

Consider a warehouse into which each customer arrives with a single failed item and expects to exchange it for a functional one. Customers’ arrivals are assumed to be distributed Poisson with rate λ , and customers are served according to a first-come–first-serve order. The warehouse is adjacent to a repair workshop so that any nonfunctional item is repaired in-house. The repair workshop has ample identical and independent workstations so that repair times are i.i.d. When repair is outsourced, it is usually the case that the repair contractor will receive a batch of failed items, repair them, and return them as a batch. This is not the case with on-site repair. Here, it is reasonable to assume that individual items from the same repair order will be returned to stock once they are repaired.
Repair orders are managed periodically, and are issued every r units of time. We denote the i’th order, i = ( 1 , 2 , . . . ) , as the order issued at time i r . Let Jane be a random customer that arrives at time i r + t with a failed item. That is, Jane arrived t units of time into the i + 1 ’s cycle. The objective of the following derivations is to compute the probability that within w units of time from Jane’s arrival (i.e., i r + t + w ), Jane has received a functional item. To this end, we track the demand and supply of functional items and observe that Jane will be served by i r + t + w if and only if at that date there is enough supply to meet the demand up to and including her arrival. The supply is S plus all the parts that were repaired on time, whereby on time we mean i r + t + w . Let q denote the number of complete cycles in the tolerable wait, and therefore ( i + q ) r < i r + t + w < ( i + q + 2 ) r . That is, i r + t + w falls either in the i + q + 1 ’th cycle or the i + q + 2 ’th cycle, depending on when Jane arrived during the i + 1 ’th cycle.
Accordingly, to correctly account for all the orders that may contribute to the supply, we distinguish between early and late arrival in the cycle. If t is sufficiently late in the cycle such that i r + t + w > ( i + q + 1 ) r , then the order issued at ( i + q + 1 ) r possibly contributes to the supply, since it may have returned before i r + t + w . In contrast, if t is sufficiently early in the cycle such that i r + t + w ( i + q + 1 ) r , then that order (i.e., the order issued at ( i + q + 1 ) r ) need not be considered for the supply. In Figure 1, we plot Jane’s arrival time (the arrow’s tail) that determines the demand, and the time by which items must return to count for the supply (the arrowhead).
Since late arrival in the cycle requires considering an additional order, it is convenient to define a binary variable a to indicate whether t is late in the cycle. Formally, a is defined as follows:
a = 0 if t + w < ( q + 1 ) r , 1 if t + w ( q + 1 ) r .

3.2. Model Notation

We summarize the model notation:
  • r: the order cycle time.
  • w: the tolerable wait.
  • i + 1 : the arrival cycle; i.e., Jane arrived during the cycle between i r and ( i + 1 ) r .
  • t: the cycle arrival time; i.e., Jane arrived t units of time after the beginning of the cycle.
  • a: A boolean variable indicating if t is late in the cycle.
  • S: the number of spares in the warehouse.
  • λ : items arrive at the warehouse following a Poisson process with rate λ .
  • L ( x ) : item repair times are i.i.d and their cumulative probability distribution function is L ( x ) , and L ( x ) = 0 if x 0 . We assume that L ( x ) has finite support; that is, there exists x * > 0 such that L ( x * ) = 1 .
  • L ¯ ( x ) : = 1 L ( x ) .
  • q : = w r is the number of complete cycles in the tolerable wait; that is, q r w ( q + 1 ) r .
  • A ( x , y ) is customer demand (i.e., the number of nonfunctional items) between x and y and that was repaired by time s + w .
  • B ( x , y ) is customer demand between x and y and that was not repaired by time s + w .
  • N ( x , y ) is customer demand that arrived between x and y (whether it was repaired or not by time s + w ). Note that by definition,
    N ( x , y ) : = A ( x , y ) + B ( x , y ) .
  • For single cycle demand, we use the notation N k : = N ( k r r , k r ) , A k = A ( k r r , k r ) , and B k : = B ( k r r , k r ) .

3.3. Demand

The demand for functional items is equal to the number of customers that arrived until s (since each customer’s demand is a single item) plus the item needed by Jane, A ( 0 , i r + t ) + B ( 0 , i r + t ) + 1 . It is more convenient to break A ( 0 , i r + t ) and B ( 0 , i r + t ) into the series of cycle demands that comprise it. That is, A ( 0 , i r + t ) = k = 1 i A k + A ( i r , i r + t ) and B ( 0 , i r + t ) : = k = 1 i B k + B ( i r , i r + t ) . Therefore, the demand is given by
D e m a n d = k = 1 i ( A k + B k ) + A ( i r , i r + t ) + B ( i r , i r + t ) + 1 .

3.4. Supply

The supply comprises three components:
  • The spares in the warehouse, S.
  • The items brought by other customers that were repaired (and therefore returned to the warehouse’s stock) by i r + t + w . As can be seen in Figure 1, when Jane arrived early in the cycle, then the latest order that could be repaired on time had to be issued by ( i + q ) r . In contrast, when Jane arrived late in the cycle, then the order issued at ( i + q + 1 ) r should also be considered. Thus, this component is given by A 0 , ( i + q + a ) r .
  • The (single) item brought by Jane if it were repaired by i r + t + w . We denote this component by the binary random variable Z, where Z is a random variable that is equal to 1 if her item is repaired on time (i.e., by i r + t + w ) and zero otherwise. Thus, the distribution of Z is
    P r [ Z = 0 ] = L ¯ ( t + w r ) ; P r [ Z = 1 ] = L ( t + w r ) .
    where, recall, we use the notation L ¯ ( x ) : = 1 L ( x ) .
Summing these components, the supply is given by
S u p p l y = S + k = 1 i A k + A ( i r , i r + t ) + A ( i r + t , ( i + q + a ) r ) + Z .

3.5. Non-Stationary Window Fill Rate

The non-stationary window fill rate is the probability for the supply to be as great as the demand. Therefore, by (3) and (5),
F N S ( s , w , S ) = P r [ D e m a n d S u p p l y ] = P r k = 1 i B k + B ( i r , i r + t ) + 1 S + A ( i r + t , ( i + q + a ) r ) + Z .
Conditioning on Z and rearranging terms results in
F N S ( s , w , S ) = P r [ Z = 0 ] P r k = 1 i B k + B ( i r , i r + t ) A ( i r + t , ( i + q + a ) r ) S 1 + P r [ Z = 1 ] P r k = 1 i B k + B ( i r , i r + t ) A ( s , ( i + q + a ) r ) S = P r [ Z = 0 ] P r [ M N S S 1 ] + P r [ Z = 1 ] P r [ M N S S ] = P r [ Z = 1 ] P r [ M N S = S ] + P r [ M N S S 1 ] ,
where
M N S : = k = 1 i B k + B ( i r , i r + t ) A ( i r + t , ( i + q + a ) r ) .
To proceed, we need to determine the distribution of M N S . This will be performed separately for each of its components.
  • k = 1 i B k : Each B k is a Poisson process with rate λ r multiplied by the probability that the k’th order was not repaired. Thus,
    B k P o i s s o n λ r L ¯ ( i r + t + w ( k + 1 ) r ) = P o i s s o n λ r L ¯ ( t + w + ( i k 1 ) r ) .
  • B ( i r , s ) : The items arriving during the interval ( i r , i r + t ) are part of the order issued at ( i + 1 ) r . Therefore, the probability for their repair on time is the probability for repair during the interval ( ( i + 1 ) r , i r + t + w ) , which is L ¯ ( t + w r ) . Thus,
    B ( i r , s ) P o i s s o n λ t L ¯ ( t + w r ) .
  • A ( i r + t , ( i + q + a ) r ) : The items that arrived during the interval ( i r + t , ( i + 1 ) r ) were sent for repair at ( i + 1 ) r and were repaired with probability L ( i r + t + w ( i + 1 ) r ) . The remaining items were sent for repair at times k r , k = i + 2 , . . . , i + q + a and their probabilities to be repaired, depending on k, are L ( i r + t + w k r ) . Therefore,
    A i r + t , ( i + q + a ) r P o i s s o n λ ( r t ) L ( t + w r ) + λ r k = i + 2 i + q + a L ( t + w + ( i k ) r ) P o i s s o n λ r k = 1 q + a L ( t + w k r ) λ t L ( t + w r ) .
    Note that if the parameter is nonpositive, such as when q + a = 0 , then A i r + t , ( i + q + a ) r = 0 .
The Poisson random variables in (8) are in different arrival intervals and are therefore independent. The sum of independent Poisson random variables is itself a Poisson random variable with a rate equal to the sum of the rates. Consequently, M N S is the difference between two Poisson random variables, where the positive Poisson is with rate k = 0 i 1 λ r L ¯ ( t + w + ( i k 1 ) r ) + λ t L ¯ ( t + w r ) and the negative Poisson is with rate k = 1 q + a λ r L ( t + w + ( i k ) r ) λ t L ( t + w r ) . The difference between Poisson random variable is a Skellam random variable [6]. Consequently,
M N S S k e l l a m λ r k = 0 i 1 L ¯ ( t + w + k r ) + t L ¯ ( t + w r ) , λ r k = 1 q + a L ( t + w k r ) t L ( t + w r ) .
Here, too, note that when q + a = 0 , then the negative Poisson vanishes and M N S P o i s s o n λ t + r k = 0 i 1 L ¯ ( t + w + k r ) . To obtain the steady state window fill rate, we need to find the mean window fill rate along an order cycle and let i . Thus,
F ( S ) = lim i 1 r 0 r F N S ( s , S ) d t .
In M N S , the value of a depends on the value of the variable of the integration t. When t [ 0 , ( q + 1 ) r w ) , then a = 0 and when t [ ( q + 1 ) r w , r ) , then a = 1 . Therefore, by (7),
F ( S ) = 1 r 0 ( q + 1 ) r w P r [ Z = 1 ] P r [ M 0 = S ] + P r [ M 0 S 1 ] d t + ( q + 1 ) r w r P r [ Z = 1 ] P r [ M 1 = S ] + P r [ M 1 S 1 ] d t ,
where letting i in (12) results in
M 0 S k e l l a m λ r k = 0 L ¯ ( t + w + k r ) + t L ¯ ( t + w r ) , λ r k = 1 q L ( t + w k r ) t L ( t + w r ) ,
M 1 S k e l l a m λ r k = 0 L ¯ ( t + w + k r ) + t L ¯ ( t + w r ) , λ r k = 1 q + 1 L ( t + w k r ) t L ( t + w r ) .
We note that when q = 0 (i.e., w < r ), then M 0 ’s distribution can be simplified to M 0 P o i s s o n λ t + r k = 0 L ¯ ( t + w + k r ) .

3.6. The Fill Rate

The fill rate is obtained by setting w = 0 in F ( S ) . The first integral in (14) vanishes, and setting w = 0 in (4) reveals that P [ Z = 1 ] = 0 . Thus, the fill rate is given by
F ( S ) = 1 r 0 r P r [ M 1 S 1 ] d t ,
where setting w = 0 in (16) renders the negative Poisson in the Skellam random variable to zero, and therefore,
M 1 P o i s s o n λ r k = 0 L ¯ ( t + k r ) + t L ¯ ( t r ) .
The first difference of the window fill rate Δ F ( S ) : = F ( S + 1 ) F ( S ) .
Δ F ( S ) = 1 r 0 r P r [ M 1 S ] P r [ M 1 S 1 ] d t = 1 r 0 r P r [ M 1 = S ] d t .
This has both computational and theoretical implications.

3.7. Computational Considerations

In (15) and (16), the expressions P r [ M 0 S 1 ] and P r [ M 1 S 1 ] are time-consuming (computation-wise) since they require computing the Skellam cumulative distribution. It is therefore more efficient to compute F ( 0 ) and then construct F ( s ) , s = 1 , . . . , S using the fact that Δ F ( s ) depends only on the probability mass function and not the cumulative distribution function. The Skellam distribution is closely tied to the Bessel function of first order [6] and therefore can be easily computed using standard software that store this function (e.g., Excel, Matlab).

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 L ( x ) , 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), N ( x , y ) = A ( x , y ) + B ( x , y ) , then (3) can be written as follows:
D e m a n d = j = 1 i N j + N ( i r , i r + t ) + 1 .
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., S u p p l y ) follows a Poisson process. In contrast, in the OR model, deliveries are of repaired orders. As a result, if the order was delivered, then A j = N j and B j = 0 , and if not, then A j = 0 and B j = N j . Consequently, S u p p l y 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 1 j as follows:
1 j = 0 if order j was not delivered on time . 1 if order j was delivered on time .
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 i r + t + w . if Jane arrived early in the cycle, then the orders relevant to the supply are all the orders issued until date ( i + q ) r , whereas if Jane arrived late in the cycle, then the order issued at ( i + q + 1 ) r must also be considered. Notice that when the i + 1 ’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
S u p p l y = S + j = 1 i N j 1 j + 0 if q = 0 and a = 0 , ( N i + 1 + 1 ) 1 i + 1 + j = i + 2 i + q + a N j 1 j otherwise .
where, recall, a = 1 when Jane’s arrival is late in the cycle and a = 0 when it is early in the cycle. The net demand ( N e t D e m ) is given by the demand minus the supply as follows:
N e t D e m = D e m a n d S u p p l y = j = 1 i N j + N ( i r , i r + t ) + 1 S j = 1 i N j 1 j 0 if q , a = 0 , ( N i + 1 + 1 ) 1 i + 1 + j = i + 2 i + q + a N j 1 j otherwise , = S + j = 1 i N j 1 ¯ j + 0 if q , a = 0 , N ( i r , i r + t ) + 1 1 ¯ i + 1 N ( i r + t , i r + r ) 1 i + 1 j = i + 2 i + q + a N j 1 j otherwise .
In the first row above, the i + 1 ’th order includes the additional item brought by Jane. In the second row, we use the notation 1 ¯ j : = 1 1 j and the fact that N i + 1 = N ( i r , i r + t ) + N ( i r + t , i r + r ) . The non-stationary window fill rate is the probability that the net demand is non-positive as follows:
F N S ( t ) = P r [ N e t D e m 0 ] .

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 i r + t + w . The probability that an arbitrary item from order j was repaired by i r + t + w is
p j : = L ( i r + t + w r j ) = L ( t + w + ( i j ) r ) .
The probability that the order is repaired (i.e., delivered) on time is itself a random variable, P j = p j N j , where, recall, N j P o i s s o n { λ r } . A notable exception is the order containing Jane’s arrival for which P j = p j N j + 1 since we must also include her item. Thus, when j i + 1 (i.e., j is not the period in which Jane arrived),
E [ P j ] = E p j N j = n = 0 n = P r [ N j = n ] p j n = e λ r n = 0 n = ( λ r p j ) n n ! = e λ r e λ r p j = e λ r 1 L ( t + w + ( i j ) r ) .
When j = i + 1 , then we must also ensure that Jane’s item has also been repaired. Therefore,
E [ P i + 1 ] = E p i + 1 N i + 1 + 1 = n = 0 n = P r [ N i + 1 = n ] p i + 1 n + 1 = L ( t + w r ) e λ r 1 L ( t + w r ) .

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 i r + t + w . 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 E [ P j ] 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 L ( · ) has a finite support, and therefore there exists a K such that for all x > r K , L ( x ) = 1 . Consequently, we can disregard dated orders that have been surely supplied since for these orders, the supply cancels with the demand (i.e., 1 ¯ j = 0 ). 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 ( a = 0 ), then we need to consider the orders between i K + q + 1 and i + q , whereas if t is late in the cycle, the relevant orders are the orders between i K + q + 2 and i + q + 1 (see Figure 2). Thus, we can rewrite N e t D e m , (23), as the following:
N e t D e m = S + j = i K + q + 1 + a i N j 1 ¯ j + 0 if q , a = 0 , N ( i r , i r + t ) + 1 1 ¯ i + 1 N ( i r + t , i r + r ) 1 i + 1 j = i + 2 i + q + a N j 1 j otherwise .
We note that if the tolerable wait is sufficiently long such that q K , then the window fill rate is 1 since order i + 1 was delivered on time, implying supply is certainly greater than the demand. We therefore proceed under the assumption that q < K . In (28), there are at most K + 1 Poisson random variables that need to be simulated. These are, N i + q K + 1 + a , . . . , N i , N ( i r , i r + t ) , N ( i r + t , i r + r ) , N i + 2 , . . . , N i + q , N i + q + a . Their distributions are N ( i r , i r + t ) P o i s s o n { λ t } , N ( i r + t , i r + r ) P o i s s o n { λ ( r t ) } and the other N j , j i + 1 are distributed N j P o i s s o n { λ r } .
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., i > K ), 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 i + 1 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 F N S ( t ) as follows:
F N S ( t ) = P r j = K 1 N j 1 ¯ j S if q , a = 0 , P r j = K + q + a 1 N j 1 ¯ j + N ( r , r + t ) + 1 1 ¯ 0 N ( r + t , 0 ) 1 0 j = 1 q 1 + a N j 1 j S otherwise .

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 d t = r / T . The arrival cycle itself, i.e., [ r , 0 ) , is simulated as a sum of T independent random variables, denoted by N 0 1 , . . . N 0 T , where N 0 j represents the demand along r + ( j 1 ) r T , r + j r T and therefore each N 0 j is distributed P o i s s o n { λ r / T } . Consequently, the demand components of Jane’s arrival cycle can be stated as follows:
N 0 = j = 1 r / d t = T N 0 j ; N ( r , r + t ) = j = 1 t / d t N 0 j ; N ( r + t , 0 ) = j = t / d t + 1 r / d t N 0 j
In summary, the window fill rate for a customer that arrived t units of time into a cycle is F N S ( t ) = P r [ N e t D e m 0 ] where
N e t D e m : = j = K + q + a 1 N j 1 ¯ j + 1 ¯ 0 1 + j = 1 t / d t N 0 j 1 0 j = t / d t + 1 r / d t N 0 j j = 1 q 1 + a N j 1 j S .
The following algorithm computes the stationary window fill rate, W F R = t = 0 r F N S ( t ) d t . 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 n q K , n q K + 1 , . . . , n q of the periodic demands N q K , N q K + 1 , . . . N q , respectively. N 0 is unique in that it is constructed by its subsegments. That is, the realization n 0 is the sum of n 0 1 , n 0 2 , . . . n 0 T , each of which are the realizations of N 0 1 , N 0 2 , . . . , N 0 T , respectively. For each simulation thread and arrival time t, the integral of F N S ( t ) is computed.
Algorithm 1: Simulation Algorithm
  • Initialization:
    (a)
    T; d t : = r T         /* integration resolution parameter */
    (b)
    Set M         /* the number of simulation threads */
    (c)
    W F R : = 0         /* W F R is the algorithm’s output */
  • Computation: For m = 1 , 2 , . . . , M do
    (a)
    Initialization:
    • F m : = 0         /* The WFR for a single realization */
    • Let n q K , n q K + 1 , . . . , n 1 , n 1 , n 2 , . . . , n q be each a realization of P o i s s o n ( λ r }
    • Let n 0 1 , n 0 2 , . . . , n 0 T be each a realization of P o i s s o n ( λ d t }
    • n 0 : = j = 1 T n 0 j , n : = ( n q K , n q K + 1 , . . . , n q )         /* Note: n includes n 0 */
    (b)
    For t : = d t , 2 d t , . . . , T d t do         /* Single Realization Computation */
    • i : = t / d t ; n 0 E : = j = 1 i n 0 j ; n 0 L : = j = i + 1 T n 0 j
    • F m : = F m + F N S ( t , n , n 0 E , n 0 L ) · d t         /* F N S ( t , n , n 0 E , n 0 L ) is defined below */
    (c)
    W F R : = W F R + F m / M
  • Output: W F R

4.3.3. Evaluating the Non-Stationary Window Fill Rate

It remains to describe how F N S ( t , n , n 0 E , n 0 L ) is defined. To this end, we note that for a given demand realization and the arrival time t, the distributions of 1 j (29) are known and are given by
1 j = 1 w . p . p j n j = L t + w ( j + 1 ) r ) n j , 0 otherwise , if j 0 , 1 w . p . p j n j + 1 = L t + w r ) n 0 + 1 , 0 otherwise , if j = 0 .
The F N S algorithm receives a single demand realization and arrival time in the cycle, t, and returns the value of (29), i.e., P r [ N e t D e m 0 ] . It considers all the possible states of the demand (i.e., the arrival status of each of the relevant orders). Since there are K + 1 relevant orders, the state space is 2 K + 1 , represented here by the numbers k = 0 , 1 , . . . , 2 K + 1 . The binary representation of k dictates the status of each of the orders. For example, if K = 4 , q = 1 , and k = 28 = 11100 2 , then ( 1 3 , 1 2 , 1 1 , 1 0 , 1 1 ) = ( 1 , 1 , 1 , 0 , 0 ) . When t is early in the cycle, the firstK (of K + 1 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 ( N e t D e m ) using the left side of the inequality in (29). The sum of the probabilities whose state meets the condition N e t D e m 0 is the value of (29) and is the output of the algorithm.
Algorithm 2:  F N S Algorithm
  • Input:
    (a)
    Demand realization: n q K , n q K + 1 , . . . , n q
    (b)
    Arrival period demand breakup: n 0 E , n 0 L         /* n 0 E + n 0 L n 0 */
    (c)
    Arrival time: t
  • Initialization:
    (a)
    If t + w < ( q + 1 ) r do: a : = 0 Else do: a : = 1
    (b)
    F N S = 0         /* F N S is the algorithm’s output */
  • Computation: For k = 0 , 1 , . . . , 2 K + 1 1 do         /* k is a state */
    (a)
    For j = 1 . . . K + 1 do: s j K + q 1 := jth digit of binary representation of k                        /* s i is order i’s delivery status */
    (b)
    For j = q K + a , q K + a + 1 , . . . , q 1 + a do:
    • If j 0 do:
      π j : = s j L t + w ( j + 1 ) r ) n j + ( 1 s j ) 1 L t + w ( j + 1 ) r ) n j
    • Else do:
      π 0 : = s 0 L t + w r n 0 + 1 + ( 1 s 0 ) 1 L t + w r n 0 + 1
    (c)
    π : = j = q K + a q 1 + a π j         /* π is the probability of the state. */
    (d)
    If q + a = 0 do:
    N e t D e m : = j = K 1 n j ( 1 s j )
    (e)
    Else do:
    N e t D e m : = j = q K + a 1 n j ( 1 s j ) + ( n 0 E + 1 ) ( 1 s 0 ) n 0 L s 0 j = 1 q 1 + a n j s j S
    (f)
    If N e t D e m 0 do: F N S : = F N S + π
  • Output: F N S

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 r = 7 days, the tolerable wait w = 5 days, the demand arrivals rate λ = 2 items per day, and the item repair times are distributed uniformly U ( 0 , 10 ) 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 M = 50 , 000 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 S = 20 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 r = 4 , then only 9 spares are needed whereas if r = 7 and r = 10 , 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 r = 4 , 7 , 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 w = 5 (left panel) and w = 10 (right panel).
In Table 2, we list the numbers of spares that are required to meet common industry performance target levels; W F R t = 80 % , 90 % , and 95 % , 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 w = 5 and the required performance level is W F R t = 80 % , 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.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Wensing, T. Periodic Review Inventory Systems: Lecture Notes in Economics and Mathematical Systems; Springer: New York, NY, USA, 2011; Volume 651. [Google Scholar]
  2. Ullah, M.; Asghar, I.; Zahid, M.; Omair, M.; AlArjani, A.; Sarkar, B. Ramification of remanufacturing in a sustainable three-echelon closed-loop supply chain management for returnable products. J. Clean. Prod. 2021, 290, 125609. [Google Scholar] [CrossRef]
  3. Ullah, M.; Sarkar, B. Recovery-channel selection in a hybrid manufacturing-remanufacturing production model with RFID and product quality. Int. J. Prod. Econ. 2020, 219, 360–374. [Google Scholar] [CrossRef]
  4. Sarkar, B.; Ullah, M.; Sarkar, M. Environmental and economic sustainability through innovative green products by remanufacturing. J. Clean. Prod. 2022, 332, 129813. [Google Scholar] [CrossRef]
  5. Dreyfuss, M.; Giat, Y. Optimal spares allocation to an exchangeable-item repair system with tolerable wait. Eur. J. Oper. Res. 2017, 261, 584–594. [Google Scholar] [CrossRef]
  6. Skellam, J.G. The frequency distribution of the difference between two Poisson variates belonging to different populations. J. R. Stat. Soc. Ser. Stat. Soc. 1946, 109, 296. [Google Scholar] [CrossRef]
  7. Vasilev, J.; Milkova, T. Optimisation Models for Inventory Management with Limited Number of Stock Items. Logistics 2022, 6, 54. [Google Scholar] [CrossRef]
  8. Scarf, P.; Syntetos, A.; Teunter, R. Joint maintenance and spare-parts inventory models: A review and discussion of practical stock-keeping rules. Ima J. Manag. Math. 2024, 35, 83–109. [Google Scholar] [CrossRef]
  9. Sherbrooke, C.C. Optimal Inventory Modeling of Systems: Multi-Echelon Techniques; Kluwer: Dordrecht, The Netherlands, 2004; Volume 72. [Google Scholar]
  10. Muckstadt, J.A. Analysis and Algorithms for Service Parts Supply Chains; Springer Science & Business Media: New York, NY, USA, 2005. [Google Scholar]
  11. Van Houtum, G.J.; Kranenburg, B. Spare Parts Inventory Control under System Availability Constraints; Springer: New York, NY, USA, 2015; Volume 227. [Google Scholar]
  12. Basten, R.J.; Ryan, J.K. The value of maintenance delay flexibility for improved spare parts inventory management. Eur. J. Oper. Res. 2019, 278, 646–657. [Google Scholar] [CrossRef]
  13. Kouki, C.; Larsen, C. Rationing policies in a spare parts inventory system with customers differentiation. Int. J. Prod. Res. 2021, 59, 6270–6290. [Google Scholar] [CrossRef]
  14. Chen, Y.K.; Gao, Q.; Su, X.B.; Fang, S.; Guo, C.M. Research on optimization of spare parts inventory policy considering maintenance priority. Int. J. Syst. Assur. Eng. Manag. 2018, 9, 1336–1345. [Google Scholar] [CrossRef]
  15. Sherbrooke, C.C. METRIC: A multi-echelon technique for recoverable item control. Oper. Res. 1968, 16, 122–141. [Google Scholar] [CrossRef]
  16. Dreyfuss, M.; Giat, Y.; Stulman, A. An analytical approach to determine the window fill rate in a repair shop with cannibalization. Comput. Oper. Res. 2018, 98, 13–23. [Google Scholar] [CrossRef]
  17. Van Donselaar, K.; Broekmeulen, R.; de Kok, T. Heuristics for setting reorder levels in periodic review inventory systems with an aggregate service constraint. Int. J. Prod. Econ. 2021, 237, 108137. [Google Scholar] [CrossRef]
  18. Gutierrez, M.; Rivera, F.A. Undershoot and order quantity probability distributions in periodic review, reorder point, order-up-to-level inventory systems with continuous demand. Appl. Math. Model. 2021, 91, 791–814. [Google Scholar] [CrossRef]
  19. Dehghani, M.; Abbasi, B.; Oliveira, F. Proactive transshipment in the blood supply chain: A stochastic programming approach. Omega 2021, 98, 102112. [Google Scholar] [CrossRef]
  20. Poormoaied, S. Inventory decision in a periodic review inventory model with two complementary products. Ann. Oper. Res. 2022, 315, 1937–1970. [Google Scholar] [CrossRef]
  21. Yıldız, B.; Sütçü, M. A variant SDDP approach for periodic-review approximately optimal pricing of a slow-moving a item in a duopoly under price protection with end-of-life return and retail fixed markdown policy. Expert Syst. Appl. 2023, 212, 118801. [Google Scholar] [CrossRef]
  22. Rodrigues, L.R.; Yoneyama, T. A spare parts inventory control model based on Prognostics and Health monitoring data under a fill rate constraint. Comput. Ind. Eng. 2020, 148, 106724. [Google Scholar] [CrossRef]
  23. Giat, Y. Allocation of scarce resources in a network with periodic deliveries and customer tolerable wait. Comput. Ind. Eng. 2021, 159, 107462. [Google Scholar] [CrossRef]
  24. Baloch, G.; Gzara, F. Inventory planning for self-serve pharmacy kiosks: A fill rate maximization approach. Comput. Ind. Eng. 2024, 188, 109836. [Google Scholar] [CrossRef]
  25. Lamghari-Idrissi, D.; Basten, R.; van Houtum, G.J. Spare parts inventory control under a fixed-term contract with a long-down constraint. Int. J. Prod. Econ. 2020, 219, 123–137. [Google Scholar] [CrossRef]
  26. Durrande-Moreau, A. Waiting for service: Ten years of empirical research. Int. J. Serv. Ind. Manag. 1999, 10, 171–194. [Google Scholar] [CrossRef]
  27. Atkinson, M.; Ntontis, E.; Neville, F.; Reicher, S. “I’ll wait for the English one”: COVID-19 vaccine country of origin, national identity, and their effects on vaccine perceptions and uptake willingness. Soc. Personal. Psychol. Compass 2023, 17, e12837. [Google Scholar] [CrossRef]
  28. Drabicki, A.; Cats, O.; Kucharski, R.; Fonzone, A.; Szarata, A. Should I stay or should I board? Willingness to wait with real-time crowding information in urban public transport. Res. Transp. Bus. Manag. 2023, 47, 100963. [Google Scholar] [CrossRef]
  29. Feng, B.; Fan, Z.P.; Li, Y. A decision method for supplier selection in multi-service outsourcing. Int. J. Prod. Econ. 2011, 132, 240–250. [Google Scholar] [CrossRef]
  30. Hathaway, B.A.; Emadi, S.M.; Deshpande, V. Personalized priority policies in call centers using past customer interaction information. Manag. Sci. 2022, 68, 2806–2823. [Google Scholar] [CrossRef]
  31. Caggiano, K.E.; Jackson, P.L.; Muckstadt, J.A.; Rappold, J.A. Efficient computation of time-based customer service levels in a multi-item, multi-echelon supply chain: A practical approach for inventory optimization. Eur. J. Oper. Res. 2009, 199, 744–749. [Google Scholar] [CrossRef]
  32. Dreyfuss, M.; Giat, Y. Optimal allocation of spares to maximize the window fill rate in a two-echelon exchangeable-item repair system. Eur. J. Oper. Res. 2018, 270, 1053–1062. [Google Scholar] [CrossRef]
  33. Dreyfuss, M.; Giat, Y. The window fill rate with nonzero assembly times: Application to a battery swapping network. In Proceedings of the Operations Research and Enterprise Systems: 6th International Conference, ICORES 2017, Porto, Portugal, 23–25 February 2017; Revised Selected Papers 6. Springer: Berlin/Heidelberg, Germany, 2018; pp. 42–62. [Google Scholar]
  34. Yang, J.; Liu, W.; Ma, K.; Yue, Z.; Zhu, A.; Guo, S. An optimal battery allocation model for battery swapping station of electric vehicles. Energy 2023, 272, 127109. [Google Scholar] [CrossRef]
  35. Xu, Y.; Serel, D.A.; Bisi, A.; Dada, M. Setting fulfillment-time guarantees for accepting customer orders in a periodic-review base-stock inventory system. IISE Trans. 2023, 1–14. [Google Scholar] [CrossRef]
  36. Van der Heijden, M.C.; De Kok, A. Customer waiting times in an (R, S) inventory system with compound Poisson demand. Z. FÜR Oper. Res. 1992, 36, 315–332. [Google Scholar]
  37. Dreyfuss, M.; Giat, Y. Allocating spares to maximize the window fill rate in a periodic review inventory system. Int. J. Prod. Econ. 2019, 214, 151–162. [Google Scholar] [CrossRef]
  38. Giat, Y.; Dreyfuss, M. The window fill rate in a periodic review inventory system with order crossover. Int. J. Prod. Res. 2021, 59, 6791–6808. [Google Scholar] [CrossRef]
  39. Bouhnik, D.; Giat, Y. ISO 9001 as a tool for improving knowledge management in business ecosystems. Int. J. Knowl. Based Dev. 2015, 6, 261–272. [Google Scholar] [CrossRef]
  40. Öztürk, H. A deterministic economic production quantity model with partial backordering and outsourcing repair of non-conforming products. Int. J. Syst. Sci. Oper. Logist. 2022, 9, 355–383. [Google Scholar] [CrossRef]
  41. Driessen, M.; Arts, J.; van Houtum, G.J.; Rustenburg, J.W.; Huisman, B. Maintenance spare parts planning and control: A framework for control and agenda for future research. Prod. Plan. Control. 2015, 26, 407–426. [Google Scholar] [CrossRef]
  42. Basten, R.J.; van der Heijden, M.C.; Schutten, J.M. Practical extensions to a minimum cost flow model for level of repair analysis. Eur. J. Oper. Res. 2011, 211, 333–342. [Google Scholar] [CrossRef]
  43. Hallak, B.K.; Nasr, W.W.; Jaber, M.Y. Re-ordering policies for inventory systems with recyclable items and stochastic demand–Outsourcing vs. in-house recycling. Omega 2021, 105, 102514. [Google Scholar] [CrossRef]
  44. Benmoussa, O. Improving Replenishment Flows Using Simulation Results: A Case Study. Logistics 2022, 6, 34. [Google Scholar] [CrossRef]
  45. Papaleonidas, C.; Androulakis, E.; Lyridis, D.V. A Simulation-based planning tool for floating storage and regasification units. Logistics 2020, 4, 31. [Google Scholar] [CrossRef]
  46. Sarkar, A.; Guchhait, R.; Sarkar, B. Application of the artificial neural network with multithreading within an inventory model under uncertainty and inflation. Int. J. Fuzzy Syst. 2022, 24, 2318–2332. [Google Scholar] [CrossRef]
  47. Sarkar, B.; Tayyab, M.; Kim, N.; Habib, M.S. Optimal production delivery policies for supplier and manufacturer in a constrained closed-loop supply chain for returnable transport packaging through metaheuristic approach. Comput. Ind. Eng. 2019, 135, 987–1003. [Google Scholar] [CrossRef]
Figure 1. Time axes of demand and supply when arrival is early and late in the cycle. D e m a n d is all the customers that arrived until the arrow’s tail plus the item required by Jane. S u p p l y is S plus all the items that were repaired until the arrowhead.
Figure 1. Time axes of demand and supply when arrival is early and late in the cycle. D e m a n d is all the customers that arrived until the arrow’s tail plus the item required by Jane. S u p p l y is S plus all the items that were repaired until the arrowhead.
Logistics 08 00034 g001
Figure 2. The relevant orders after re-indexing. When Jane arrived early in the cycle, then any order prior to order q K is surely delivered on time and is therefore irrelevant. In contrast, when Jane arrived late in the cycle, then also order q K is surely delivered on time and is irrelevant.
Figure 2. The relevant orders after re-indexing. When Jane arrived early in the cycle, then any order prior to order q K is surely delivered on time and is therefore irrelevant. In contrast, when Jane arrived late in the cycle, then also order q K is surely delivered on time and is irrelevant.
Logistics 08 00034 g002
Figure 3. The window fill rate as a function of the spares when repair is in-house (IR model) for different values of cycle time. The other parameters’ values are when λ = 2 , L ( x ) U ( 0 , 10 ) , and w = 5 . The dotted line expresses the number of spares needed to maintain an 80% window fill rate.
Figure 3. The window fill rate as a function of the spares when repair is in-house (IR model) for different values of cycle time. The other parameters’ values are when λ = 2 , L ( x ) U ( 0 , 10 ) , and w = 5 . The dotted line expresses the number of spares needed to maintain an 80% window fill rate.
Logistics 08 00034 g003
Figure 4. The window fill rate as a function of the spares when repair is outsourced (OR model) for different values of cycle time. The other parameters’ values are when λ = 2 , L ( x ) U ( 0 , 10 ) , and w = 5 . The dotted line expresses the number of spares needed to maintain an 80% window fill rate.
Figure 4. The window fill rate as a function of the spares when repair is outsourced (OR model) for different values of cycle time. The other parameters’ values are when λ = 2 , L ( x ) U ( 0 , 10 ) , and w = 5 . The dotted line expresses the number of spares needed to maintain an 80% window fill rate.
Logistics 08 00034 g004
Figure 5. The window fill rate as a function of the spares when repair is in-house (IR model). In the left panel, the tolerable wait is w = 5 and in the right panel, the tolerable wait is w = 7 . The other parameters’ values are λ = 2 , L ( x ) U ( 0 , 10 ) , and r = 7 .
Figure 5. The window fill rate as a function of the spares when repair is in-house (IR model). In the left panel, the tolerable wait is w = 5 and in the right panel, the tolerable wait is w = 7 . The other parameters’ values are λ = 2 , L ( x ) U ( 0 , 10 ) , and r = 7 .
Logistics 08 00034 g005
Table 1. The window fill rate of the IR and OR model in the baseline case for different spare levels.
Table 1. The window fill rate of the IR and OR model in the baseline case for different spare levels.
Model S = 0S = 5S = 10S = 15S = 20S = 25S = 30
IR 0.3%14.1%54.4%86.5%98.3%99.9%100.0%
ORMean0.0%1.0%13.4%42.5%72.2%91.6%98.4%
SD0.00%0.05%0.43%0.74%0.77%0.25%0.07%
Notes: For the IR model, the window fill rate is reported. For the OR model, the mean and standard deviation (SD) of 30 simulation executions of the window fill rate are reported.
Table 2. The number of spares required to meet the target window fill rate in each model and the cost of outsourcing.
Table 2. The number of spares required to meet the target window fill rate in each model and the cost of outsourcing.
Tolerable WaitSpares WFR t = 80 % WFR t = 90 % WFR t = 95 %
w = 2 OR model233234
IR model202729
Cost355
w = 5 OR model222527
IR model141618
Cost899
w = 8 OR model161821
IR model457
Cost121314
Notes: Cost is the difference between the number of spares in the OR model and the number of spares in the IR model.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Giat, Y. Stock Levels and Repair Sourcing in a Periodic Review Exchangeable Item Repair System. Logistics 2024, 8, 34. https://doi.org/10.3390/logistics8020034

AMA Style

Giat Y. Stock Levels and Repair Sourcing in a Periodic Review Exchangeable Item Repair System. Logistics. 2024; 8(2):34. https://doi.org/10.3390/logistics8020034

Chicago/Turabian Style

Giat, Yahel. 2024. "Stock Levels and Repair Sourcing in a Periodic Review Exchangeable Item Repair System" Logistics 8, no. 2: 34. https://doi.org/10.3390/logistics8020034

APA Style

Giat, Y. (2024). Stock Levels and Repair Sourcing in a Periodic Review Exchangeable Item Repair System. Logistics, 8(2), 34. https://doi.org/10.3390/logistics8020034

Article Metrics

Back to TopTop