Determination of Optimal Shipping Quantity for Perishable Goods under Probabilistic Supply

The problems of shipping and transporting perishable goods are commonly considered in the literature as significant topics, but rarely did researchers adopt a probabilistic point of view in their models. It is common in SCM environments that the participating entities’ behaviors are random and unpredictable and so can only be modeled in a probabilistic way. In this paper, we consider the shipping problem of determining the optimal quantity of perishable products with a limited time to be stored in the warehouse.The optimal quantity minimizes the overall operational costs including those of inventory and shipping. We develop a mathematical model, from which the probability distribution function, mean, and variance of the length of the build-up period are derived and we establish a cost function for determining the optimal shipping value.


Introduction
The shipment of goods, in general, aims for cost-effectiveness by loading as many items as possible in a restricted space of transportation.Another goal is the on-time delivery of specific amounts of goods.However, shipping and transporting of perishable goods such as food, biological materials (blood), and hazardous products (bacterial and fungal items) are quite different from shipping inert, nonperishable goods.Shipment of perishable items must involve consideration of the timeline of expiration.Thus, it also includes the goal of shipping and delivering the perishable items before expiration.
When several suppliers providing perishable products participate in supply chain management (SCM) circumstances, there can exist a trade-off between cost of shipping and risk of deterioration of the perishable goods.Specifically, suppose a typical SCM situation in which many suppliers use a warehouse to store their goods until they are delivered to another place, usually in batches.To reduce shipping costs, the amount of goods in shipment should be as large as feasible.It is therefore desirable to wait until enough products have accumulated in the warehouse to ship them together.However, this approach inevitably results in some products that arrive earlier at the warehouse deteriorating before being shipped.If the shipping quantity is large, shipping costs can be greatly reduced; however, the possibility of product deterioration will be much higher because of the time lag.Therefore, it is vital to determine a proper quantity of shipping that balances the two costs and minimizes the overall cost.In this study, we address this problem and study ways to find optimal shipping amounts of perishable goods under probabilistic supply in SCM.
There have been many studies on perishable goods in supply chain management (SCM), involving various perspectives, such as technical, heuristic, analytical, and strategic points of view.Grunow and Piramuthu [1] showed that, under some conditions, the incorporation of radio frequency identification (RFID) can benefit distributors, retailers, and consumers.Recently, Todorovic et al. [2] studied the use of RFID tags.These studies mainly focused on the usage of RFID technology for the logistics of perishable products, not on the mathematical analysis.
van Donselaar et al. [3] worked with two Dutch supermarket chains and improved the automated store ordering (ASO) systems by considering perishable items.They classified the items in the supermarkets and noted that different items require different inventory control rules, designing an inventory control strategy for perishable items.Ignaciuk [4] considered the problem of establishing an efficient control strategy for production-inventory systems in which the stock 2 Mathematical Problems in Engineering replenishment process is unreliable.Considering supply chain strategies for perishable goods, Blackburn and Scudder [5] used melons and sweet corn as two examples and showed that an appropriate model to minimize loss in the supply chain is a hybrid of a responsive model from postharvest to cooling.In these studies, they introduced several inventory strategies for the perishable products but did not consider optimal shipping quantities for them.
A set of agricultural suppliers with low demand that could reduce long-haul transportation costs by consolidating their products was considered in Nguyen et al. [6].They suggested a look-ahead heuristic algorithm taking advantage of economies of scale by aiming to ship larger quantities.Zanoni and Zavanella [7] examined the challenge of shipping a set of perishable products from a single vender to a common buyer.They developed a mixed integer linear programming model incorporating a well-known heuristic algorithm with some modification.Thangam [8] considered optimal price discounting and lot-sizing policies for perishable items in an SCM environment.In his research, an economic order quantity-(EOQ-) based model with perishable items was developed under an advance payment (AP) scheme and a two-echelon trade credit option.In these researches, they considered the control of perishable products by using the heuristic algorithms but failed to find any mathematical relationship between length of build-up period for shipping and expiration of products.Ali et al. [9] dealt with a supply chain of a food in a different perspective; that is, they showed that supply chain integration can mitigate the halal food integrity risk.
Hsieh and Dye [10] suggested a deterministic, stockdependent inventory model for deteriorating items and provided procedures for determining the maximum total profit per unit time.They assumed the storage space for deteriorating items is finite.Gumasta et al. [11] developed a transportation model mapped onto an inventory model with time-varying demand and two types of perishable goods.The model was developed to maximize the revenue and minimize transportation and inventory costs.They consider the relationships between the transportation problem and the inventory model assuming the fixed shortage cost, and selling price of products decreasing with time.
In these studies, the problems of shipping and transporting perishable goods are commonly considered as significant topics, but rarely did researchers adopt a probabilistic point of view in their models.It is common in SCM environments that the participating entities' behaviors are random and unpredictable and so can only be modeled in a probabilistic way.In this paper, we consider the shipping problem of determining the optimal quantity of perishable products with a limited time to be stored in the warehouse.The optimal quantity minimizes the overall operational costs including those of inventory and shipping.We develop a mathematical model by adopting an approach similar to the queuing system with what we call "impatient" products under the -policy, which is described in the following sections.
The paper is organized as follows.Section 2 explains the mathematical model considered in this paper, and Section 3 derives the probability distributions of the build-up period conditioned on the quantity of shipping, with which mean and variance of the period are calculated.Section 4 establishes a cost function that includes inventory holding cost and one-time shipping cost as a function of shipping quantity.This cost function can be used to determine the optimal quantity of one-time shipping.A summary and conclusion are provided in Section 6.

Model
We consider the problem of determining the optimal shipping quantity of perishable goods that minimizes the overall inventory and shipping costs and model it similar to a queuing system.Perishable goods arriving at a warehouse are treated as impatient customers arriving at a service station, so that shipment of these goods is a service provided to the waiting customers in the system, as depicted in Figure 1.
In our model, perishable goods arrive at the warehouse according to the Poisson process with a known positive constant rate  and are stored for later shipment with many other products.Shipping is performed right after the quantity of goods accumulated in the warehouse reaches a predetermined level, say , and we call the storage duration the build-up period.The shipment time for the perishable goods arriving during the build-up period, however, is not certain because it depends on the arrival of goods, which is probabilistic.In this period, stored products might go bad and will be unusable if they have to wait in the warehouse for longer than a fixed time .Those products must then be discarded, resulting in inventory loss.
Figure 2 illustrates a typical sample path of the number of products in the warehouse, assuming  = 9.At point ', the warehouse is empty, and a new build-up period starts.During the build-up period, new products are held in the warehouse until shipment.At point A, the first product is removed from the warehouse because its waiting time exceeded the expiration date .The second product is removed from the warehouse at point B for the same reason.At point C, the number of stored products reaches  before the third item's waiting time reaches , so the items are shipped.In this example, a total of  + 2 = 11 products have arrived at the warehouse, and two products have deteriorated and been scrapped at points A and B. The length of the build-up period is the time it takes to accumulate  products in the warehouse and depends on the arrivals and .

Duration of the Build-Up Period
In this section, we develop a mathematical model representing the number of products stored in the warehouse, from which the probability distribution and length of buildup period are derived.We introduce the following random variables for use in our model.Let  0 be the arrival epoch of the first product after the system is empty, and let   ( = 1, 2, . ..) be the time interval between the th and  + 1th arrival epochs.  's are assumed to be independent of each other and exponentially distributed with a common mean of 1/.Additionally, let   be the random variable representing the total elapsed time, starting at , until the number of  products stored in the warehouse reaches  for the first time, and let  *  () denote the Laplace transform (LT) of the probability density function of   .Then, , the length of the build-up period, is the sum of  0 and   0 , which are independent of each other.
We introduce the concept of a leading product, which is the one whose waiting time is the longest among all products in the warehouse.The first product arrives at the warehouse and becomes the leading product.Products continue to arrive at the warehouse, accumulating until a certain inventory threshold for shipping is reached.If time  has passed but the shipment has not yet occurred, the first product expires and is discarded; the second product then becomes the leading product.Such successive transfers of the leading product position continue until the build-up period ends.
We determine the probability distribution of   0 by conditioning it on two mutually exclusive events based on the status of the first leading products as follows:  1 : event in which the leading product has waited for time  but shipping has not yet occurred because there are still fewer than  products in the warehouse.When  1 occurs, the leading product is discarded, and the next product becomes the leading product.We assume  ≥ 2 because  1 = 0 if  = 1. 0 : complementary event of  1 .The event  0 occurs when  − 1 or more products arrive at the system during time  after the leading product's arrival.
For notational convenience, we introduce the following notations: where   is the probability that, at most,  products arrive during time .Then, the probabilities of events  0 and  1 are given as follows: Pr Because the leading product is removed from the system without being shipped with probability  −2 , the number of leading products that leave before shipping begins is geometrically distributed with an expected value of We observe some properties of   0 under the conditions of either  1 or  0 .When  0 occurs, the following holds.
Property 1.For a given  0 ,   0 is the sum of  − 1 independent interarrivals; that is, For  1 , the leading product is changed, and the following holds.Property 2 holds because the second product is now the new leading product whose arrival epoch becomes the new starting point.Properties 3, 5, and 6 are obvious when the arrivals follow a Poisson process and the interarrival times are exponential.Therefore, the probabilistic behavior repeats itself immediately after each arrival.Finally, Property 4 holds because, as long as  ≥ 2, the first arrival epoch after the leading product is independent of event  1 .

Property 2. 𝑋 𝐴
Then, by the above properties, we have the following: ⋅ Pr ( 0 ) (by Property 6) . ( Because  1 +  2 + ⋅ ⋅ ⋅ +  −1 is a truncated Erlang random variable when event  0 occurs, we obtain the following: Then, LT  *  0 () of   0 is given by (5) Mean value and variance of   0 are as follows: In order to obtain the probability distribution of , we divide the build-up period into two periods: unsuccessful build-up period and successful build-up period.The former is the sum of interarrival times of leading products, and the latter is the duration from the arrival of the last leading product to the start of the busy period, as shown in Figure 3.
Therefore,   0 in the above equations can be expressed as a sum of two random variables,   0 ,fail , representing the unsuccessful build-up period, and   0 ,success , representing the successful build-up period.As indicated, E(  0 ,fail ) is the expected number of discarded products multiplied by the expected interarrival time.Then, we have the following: Now, LT,  * (); the expected value, E(); and the variance, Var(), of the build-up period  =  0 +   0 ,fail +   0 ,success are as follows:

Cost Function for the Optimal Value of 𝑁
In this section, we introduce a cost function to determine the optimal value of  by minimizing the total cost per unit time during a cycle.Let   represent the shipping cost incurred whenever the products are shipped after the number of waiting products reaches .Also, let  ℎ represent the inventory holding cost per product per unit time in the warehouse.A larger  value incurs a smaller shipping cost per product with a larger holding cost per unit time, and vice versa.Therefore, the optimal value of  should establish a balance between these two costs in order to minimize the long-run expected total cost per unit time.
For the total holding cost during the build-up period, we consider the two periods   0 ,fail and   0 ,success , introduced in the previous section.Figure 3 illustrates the cumulative numbers of arrivals (Arr()) and departures (Dep()), respectively, up to time  during the build-up period.The area between Arr( 0 +   0 ) and Dep( 0 +   0 ) represents the total cumulative waiting time of the products during the build-up period.
This area can be divided into two parts: the area of the perished products (corresponding to the shaded area in Figure 3) and the area of the products existing in the warehouse when shipping begins (corresponding to the dotted area in Figure 3).The first part is simply  multiplied by the expected number of perished products  −2 /(1− −2 ) because all perished products have waited in the warehouse for exactly time .The second area is determined as follows: given   0 = , the  − 2 arrival epochs are uniformly distributed on [0, ], and if we denote   ( = 0, 1, . . .,  − 1) as the th arrival epoch after the previous leading product, then Therefore, the total holding cost ( ℎ build-up ) during the buildup period is given as Note that when  → ∞, products do not expire.In this case,  1 = 0; therefore,   =   = 0, for all  and .Then,  * () = (/( + ))  and  ℎ build-up =  ℎ ⋅ ( − 1)/(2), which coincide with the result obtained for the M/G/1 queueing system under -policy without perishable products [12].Now, the long-run total expected cost per unit time, TC(),  ≥ 2, is given as The behavior of the function TC() with regard to  is not easily determined because of its complicated structure.However, the global point minimizing TC() can be obtained numerically.In practical situations, if the empirical range  1 ≤  ≤  2 is given, then the following can be used to identify the set  (if  = 0, then let  = { 1 ,  2 }), and TC() is compared to each  ∈  in order to obtain the local optimal value of threshold  within the range:

Numerical Illustration
We present some numerical examples to illustrate the analysis in this section.Figure 4 shows the change of long-run expected total cost TC() per unit time as a function of  when   / ℎ = 100,  = 8, and  = 3.By plugging the numbers into (1) and (11), we obtain the curve of total cost per unit time, from which the optimal  * = 29 and the optimal total cost TC( * ) = 46.42 can be calculated.More results of the numerical examples are provided in Table 1 under various parameter values.Specifically, two different values of cost ratio   / ℎ = 10, 100, two different values of perishable products' endurable waiting time  = 1, 3, and two different Poisson arrival rates  = 5,8 are considered.
From Table 1, we can see that when the ratio   / ℎ is bigger, the optimal  * also becomes larger.This is because less shipping is preferable in order to reduce the shipping cost.If we compare the cases  = 1 and  = 3, it can be seen that the perishable products' endurable waiting time does not give much effect on  * when   / ℎ is small.That is, when   / ℎ = 10, then  * = 8 (case  = 5) and 10 (case  = 8) under  = 1 and these optimal values become  * = 10 and 13, respectively, as  increases to 3.But when   / ℎ = 100, remarkable increases in  * are shown if  changes from 1 to 3, that is, 14 to 23, and 19 to 29.Therefore, we conclude that the optimal  * becomes more sensitive to the perishable products' endurable waiting time  as the ratio   / ℎ gets bigger.Lastly, we can see by comparing the cases of  = 5 (medium traffic) and  = 8 (heavy traffic) that optimal  * should be obviously bigger as the products arrive more often.

Conclusion
In this paper, we considered a mathematical model to derive the optimal shipping quantity of perishable goods in SCM circumstances.We modeled the system using an approach similar to the queuing system under -policy with impatient customers.All products that wait for longer than a fixed time expire and are removed from the warehouse as loss.
In this model, we derive the probability distribution function, mean, and variance of the length of the buildup period and establish a cost function for determining the optimal shipping value .Because an analytical solution for the optimal value of  is not provided, extensive numerical experiments are necessary to reveal the relations among parameters; such experiments are currently being performed by the authors.

Figure 1 :Figure 2 :
Figure 1: A simplified comparison of a supply chain with a queuing system.

Figure 3 :
Figure 3: Unsuccessful build-up period and successful build-up period.

Table 1 :
Optimal  * for various parameter values.