Next Article in Journal
On the Solvability of a Mixed Problem for a High-Order Partial Differential Equation with Fractional Derivatives with Respect to Time, with Laplace Operators with Spatial Variables and Nonlocal Boundary Conditions in Sobolev Classes
Next Article in Special Issue
A Decision Support System for Dynamic Job-Shop Scheduling Using Real-Time Data with Simulation
Previous Article in Journal
Regular Ordered Ternary Semigroups in Terms of Bipolar Fuzzy Ideals
Previous Article in Special Issue
Effects of Preservation Technology Investment on Waste Generation in a Two-Echelon Supply Chain Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Low Carbon Supply Chain Coordination for Imperfect Quality Deteriorating Items

1
Department of Industrial and Systems Engineering, Chung Yuan Christian University, 200 Chung-Pei Rd., 32023 Chung-li, Taoyuan, Taiwan
2
Department of Industrial Engineering, Universitas Atma Jaya Yogyakarta, Jl. Babarsari 43, 55281 Yogyakarta, Indonesia
3
Department of Industrial Engineering, Petra Christian University, Surabaya, East Java 60236, Indonesia
*
Author to whom correspondence should be addressed.
Mathematics 2019, 7(3), 234; https://doi.org/10.3390/math7030234
Submission received: 19 January 2019 / Revised: 18 February 2019 / Accepted: 20 February 2019 / Published: 5 March 2019

Abstract

:
Nowadays, many countries have implemented carbon pricing policies. Hence, the industry adapts to this policy while striving for its main goal of maximizing financial benefits. Here, we study a single manufacturer–retailer inventory decision considering carbon emission cost and item deterioration for an imperfect production system. This study examines two models considering two cases of quality inspection. The first is when the buyer performs the quality inspection, and the second is when the quality inspection becomes the vendor’s responsibility so that no defective products are passed to the buyer. Carbon emission costs are incorporated under a carbon tax policy, and we consider the carbon footprint from transporting and warehousing the items. The objective is to jointly optimize the delivery quantity and number of deliveries per production cycle that minimize the total cost and reduce the total carbon emissions. This study provides solution procedures to solve the models, as well as two numerical examples.

1. Introduction

Supply chain coordination has a favorable effect on inventory replenishment decisions. Supply chain coordination can be realized through information sharing and joint decision-making. Coordination brings many advantages such as lower inventory-related costs and quality improvement [1]. This study considers supply chain management coordination and examines its effect on both economic and environmental performance. This study proposes supply chain inventory models that consider carbon emission costs and the existence of defective items under different inspection coordination mechanisms. Further, the models also consider the effect of item deterioration. In real life, many inventory items deteriorate over time due to spoilage, physical depletion, or obsolescence.
Due to increasing pressure from legislation, customers, and other organizations, business and industry are striving for more eco-friendly operation. The production, distribution, consumption, and other post-consumption processes of a product are sources of carbon emission. Therefore, the concept of a low-carbon supply chain has gained massive interest among researchers and industry practitioners [2,3]. The objective is to control and reduce CO2 emissions (the major part of greenhouse gas emission) from the supply chain. Recently, Kazemi et al. [4] considered the effect of carbon emissions on several economic order quantity (EOQ) models. Sarkar et al. [5] considered warehouse emissions in the EOQ model with a rework for the defective items. The model also considered partial backorder and multi-trade-credit-period. Taleizadeh et al. [6] proposed economic production quantity (EPQ) models that considered carbon emissions. Recently, Sarkar et al. [7,8] and Daryanto and Wee [9] incorporated a carbon tax in a supply chain total cost model. Wahab et al. [10], Jauhari et al. [11], Sarkar et al. [12], Jauhari [13], Gautam and Khanna [14], and Tiwari et al. [15] incorporated both carbon emissions and imperfect quality in a low-carbon supply chain model. The quality inspection is performed by the buyer, and the defective products are sent back to the vendor or sold into the secondary market at a discounted price.
Table 1 illustrates the research gap by comparing this paper with the existing literature. This study focuses on supply chain inventory models for a system that contains imperfect quality items. The decisions of the supply chain dealing with the defective items in the imperfect production processes affect carbon emissions, because defective item processing also adds to the total emissions. Moreover, the loss due to imperfect quality and deterioration also forces the manufacturer to produce more products to satisfy customer demand per period, resulting in the increase in carbon emission from production, holding, and distribution. The objectives of these studies are to simultaneously minimize the total cost and reduce carbon emissions. This paper also contributes to low-carbon supply chain models by considering two cases of quality inspection. In the first case, the buyer performs the quality inspection, and in the second case, the quality inspection becomes the vendor’s responsibility. The first model extends the studies of Wahab et al. [10], Jauhari et al. [11], Sarkar et al. [12], Jauhari [13], and Gautam and Khanna [14] to consider the effect of deterioration. In addition to the fixed and variable inspection costs, the model also extends Tiwari et al.’s [15] model by introducing weight and distance-dependent transportation cost and emission variables. The second model extends the first model by introducing an inspection option to prevent defective products from being shipped to the buyer. This model reduces the expected total costs and emission costs of the supply chain.
This study incorporates carbon emissions, item deterioration, and defective percentage to guide the supply chain managers to make the inventory decisions on the delivery size and the number of deliveries per cycle. This introduction section is followed by reviews of previous related studies in Section 2. Then, Section 3 defines the problem, assumptions, and notations in this study. Section 4 presents two mathematical models. Section 5 provides two numerical examples and the sensitivity analysis to find some insights from the proposed models. In the end, Section 6 summarizes the findings and discusses some opportunities for further research.

2. Literature Review

This study incorporates both economic and emission costs in a two-echelon supply chain production–inventory model assuming that defective products exist in each delivered lot. This section presents the existing literature that supports this study.

2.1. Imperfect Quality Inventory Model

In many industries, production systems are imperfect, producing a certain percentage of defective products. Rosenblatt and Lee [16] and Porteus [17] studied the relationship between the optimal lot size and quality performance. Rosenblatt and Lee [16] studied the optimal production cycle through considering the proportion of defective items, while Porteus [17] related the model to the opportunity for quality improvement and setup cost reduction through investment. Salameh and Jaber [18] incorporated defective items into the EPQ model and considered the screening time and cost. Many other researchers have continued the research on the EOQ and EPQ models with imperfect quality. Those researchers assume that at the end of the screening period, or the end of the cycle, the defective products will be sold at a lower price.
Other researchers bring the effect of imperfect quality items into the integrated vendor–buyer or multi-echelon inventory model. Huang [19] considered imperfect quality and assumed that the vendor provides a product warranty for the defective items. The buyer conducts a 100% inspection, and at the end, the vendor treats the faulty items. Goyal et al. [20] extended the previous model with a single production and multiple deliveries containing defective products. They assumed that the buyer sells the defective items at a discounted price. Figure 1 illustrates the scenario. Wee et al. [21] considered imperfect quality, shortage backorder, and item deterioration in an integrated production–inventory model. Lee and Kim [22] also developed an integrated production–inventory model of imperfect quality deteriorating items.
Bazan et al. [23] studied the effect of imperfect quality in different vendor–buyer inventory models. In their study, the vendor performs the inspection and considers one of three possible decisions regarding the defective items: (1) scrap off, (2) salvage at a discounted price, and (3) rework. Figure 2 illustrates the scenario when the vendor performs the inspection. Sarkar et al. [24] studied the integrated inventory model with two-stage inspection by the vendor considering the rework process and variable transportation cost. Yu and Hsu [25] developed a production–inventory model in which the defective items are returned to the vendor immediately for rework.

2.2. Low-Carbon Supply Chain Management

Research on low-carbon supply chain management has increased rapidly in recent years and has been marked by a surge in the amount of literature in this area. Much of this research revealed that supply chain collaboration and the adjustment on operational decisions could reduce carbon emissions without significantly increasing their costs. Wahab et al. [10] studied the optimal shipment size and number of shipments for a two-echelon supply chain with carbon emission cost from transporting the inventory. The emissions are affected by the distance traveled, vehicle fuel efficiency, and the actual shipment weight. Benjaafar et al. [26] modified the traditional supply chain model by associating the carbon footprint from placing an order to the supplier, production setup, production process, and inventory holding. Fahimnia et al. [27] studied the impact of carbon pricing on a closed-loop supply chain through a case study.
Bozorgi et al. [28] considered carbon emissions from transporting and storing cold items that require temperature-controlled trucks and freezers. Bozorgi [29] extended the previous model considering multi-product cold items under limited capacity. Hariga et al. [30] incorporated carbon emissions from transporting and storing the cold items in a three-echelon supply chain. Ghosh et al. [31] considered carbon emissions from production, inventory holding, and transportation in a vendor–buyer supply chain under a single setup and multiple deliveries policy. Toptal and Çetinkaya [32] studied the effect of supply chain coordination and carbon emissions on vendor–buyer order quantity under lot-for-lot delivery. Bouchery et al. [33] examined different supply chain coordination configurations considering carbon emissions under limited vehicle capacity. Dwicahyani et al. [34] incorporated carbon emission costs, energy cost, and waste disposal for a two-echelon supply chain with remanufacturing. Li et al. [35] considered joint carbon tax and cap-and-trade policies for a two-echelon supply chain production–distribution model with transportation outsourcing. Wangsa [36] incorporated the government’s penalties and incentive policies to reduce carbon emissions.
Anvar et al. [37] considered emissions from transportation and inventory-holding activities in a one-supplier multi-retailer supply chain. Hariga et al. [38] considered carbon tax and carbon cap policies for a two-echelon supply chain with vendor-managed consignment inventory partnership. The model incorporated emissions from the ordering process, production setup, and holding the inventory. Ji et al. [39] considered a carbon reduction investment from the supplier to get higher customer demand. Wang and Ye [40] compared the effect of considering carbon emissions on just-in-time and economic order quantity decisions for two-echelon supply chain inventory models. Ghosh et al. [41] considered a carbon tax regulation to minimize the total expected cost of a supply chain under stochastic demand and shortage backorder. Ma et al. [42] considered the effect of the carbon tax scheme between suppliers and buyer for production, procurement, and pricing decisions. Darom et al. [43] developed a manufacturer–retailer inventory model considering disruption risks and recovery with safety stock and the effect of carbon emission costs. The model considered carbon emissions from the transportation activities for a better recovery plan. Huang et al. [44] studied inventory and pricing decisions considering carbon emission, production disruption, and controllable deterioration using preservation technology. Recently, Daryanto et al. [45] proposed a low-carbon three-echelon supply chain inventory model considering item deterioration. Kundu and Chakrabarti [46] developed a low-carbon supply chain inventory model taking into account the effect of inflation and the time-value of money. Other researchers incorporated carbon emissions in supplier selection and order allocation [47,48,49].

3. Problem Definition, Assumption, and Notations

3.1. Problem Definition

This study considers a manufacturer–retailer supply chain that produces one type of item sold solely through one channel. The retailer orders n deliveries of equal lot size (Q) per cycle. The manufacturer implements single-setup multiple-deliveries (SSMD). Hence, it produces nQ units of item per production cycle. This study develops two models considering two cases of quality inspection. (1) In the first, the buyer performs a complete quality inspection process. (2) In the second, the vendor performs the quality inspection so that no defective products are passed to the buyer. The defective items are sold at a discounted price with no additional cost in both scenarios. Both the vendor and the buyer consider the carbon emission costs in their decision to comply with the carbon tax regulation. Model 1 is an extension of Tiwari et al.’s [15] model by introducing weight and distance-dependent transportation costs in addition to the fixed and variable inspection costs. Later, the second model extends the first model by studying another inspection option to reduce the expected total costs and emissions.

3.2. Assumption and Notation

This study explores real-life problems of cost minimization and carbon emission reduction under certain assumptions of a controlled situation. The assumptions are listed below, while the notations are presented in Table 2.
  • The retailer’s demand rate and the manufacturer’s production rate are known and constant.
  • The manufacturer implements a single-setup multiple-deliveries (SSMD) policy. Based on the retailer’s order, the manufacturer produces nQ units of item per production cycle to reduce the setup time and cost. Then, it delivers the item in an equal lot sizes and constant time intervals [50].
  • The replenishment is instantaneous.
  • The items deteriorate in the manufacturer and retailer’s inventory. The deterioration rate for both the manufacturer and retailer are equal and constant.
  • The defective percentage, u, has a uniform distribution where 0 ≤ α < β < 1.
  • Good products are always available during the quality inspection as x > D.
  • The retailer (in Model 1) and the manufacturer (in Model 2) perform a 100% quality inspection to ensure an excellent service.
  • The fixed inspection cost per cycle is constant, whether performed by the buyer or the manufacturer.
  • Carbon emissions come from the fuel and electricity consumption during transporting and holding the inventory.
  • Shortage is not considered.
  • The additional fuel consumption is a linear function of truckloads. Figure 3 illustrates the linear fuel consumption model, which is similar to that of Hariga et al. [30] with an example of the dataset from Volvo Corporation [51].

4. Model Development

This section provides model development for two inspection cases.

4.1. Model Development with Retailer Inspection

This sub-section presents model development when quality inspection becomes the retailer’s responsibility. The model is adapted from Tiwari et al. [15]. However, this study considers weight and distance-dependent transportation cost.
The inventory level of the manufacturer and the retailer is illustrated in Figure 4. In one production cycle, the manufacturer produces PT1 units of the item, and delivers all the produced items to the retailer n times with a constant lot size Q. Right after receiving each lot, the retailer starts the quality inspection that ends at Ti. At Ti, uQ units of defective products will be removed from the inventory. During the period [0, T/n], the retailer’s inventory decreases due to demand.

4.1.1. Retailer Cost and Emission

As c is the retailer’s ordering cost, the ordering cost per year is given by c/T (Lee and Kim [22], Yang and Wee [52]). After a lot arrived, the 100% quality inspection starts and then finishes at Ti. As there are fixed inspection costs per delivery, ic, and unit inspection costs, uc (Sarkar et al. [12]), the inspection cost per year is given by:
i c n T + u c Q n T = n T ( i c + u c Q )
Figure 5 illustrates the retailer’s on-hand inventory per delivery cycle. The inventory level, considering deterioration and imperfect quality, has been studied by researchers such as Jaggi et al. [53] and results in the following equations:
I d ( t ) = Q e θ   t + D θ ( e θ   t 1 )   ,   0   t   Q / x
I d ( t ) = Q e θ   t + D θ ( e θ   t 1 ) u Q   , Q / x   t   T / n
At t = T/n, Id(T/n) = 0, therefore equation (3) becomes:
Q e θ T / n + D θ ( e θ T / n 1 ) u Q = 0
Solving the equation for Q, one has:
Q = D ( e θ T / n 1 ) θ ( e θ T / n u ) = D ( e θ T / n 1 ) θ ( u e θ T / n )   which   is   = D ( e θ T / n 1 ) θ ( 1 u e θ T / n )   ( Jaggi   et   al .   [ 53 ] )
Hence:
Q = I d ( 0 ) = D ( e θ T / n 1 ) θ ( 1 u e θ T / n )
Further, the on-hand inventory per year for the retailer is:
n T [ 0 Q / x I d ( t ) d t + Q / x T / n I d ( t ) d t ]
n T [ Q θ ( 1 e θ Q x ) D θ 2 ( Q θ x + e θ Q x 1 ) 1 θ ( e θ T n e θ Q x ) ( Q + D θ ) ( T n Q x ) ( u Q + D θ ) ]
We assumed that storing the items requires electrical energy with a certain amount of carbon footprint. Therefore, substituting Equations (4) with (5) and considering the retailer’s holding costs and emissions, the inventory holding and emission cost per year becomes:
( h d + w e ) n T [ D θ 2 ( u e θ T n 1 ) 2 ( u D x ( 1 + e 2 θ T n 2 e θ T n ) + u ( 1 e 2 θ T n + u e 2 θ T n u e θ T n ) + θ T n ( u 1 + u e θ T n u 2 e θ T n ) + e θ T n 1 ) ]
Considering the expected probability value of the defective products (E[u]), the average warehouse energy consumption per unit product (ec), and the average emission from electricity generation (Ee), from Equation (6), the retailer’s expected carbon footprint (ETEd) per year of holding the inventory is:
E T E d = ( e c E e ) n T [ D θ 2 ( E [ u ] e θ T n 1 ) 2 ( E [ u ] D x ( 1 + e 2 θ T n 2 e θ T n ) + E [ u ] ( 1 e 2 θ T n + E [ u ] e 2 θ T n E [ u ] e θ T n ) + θ T n ( E [ u ] 1 + E [ u ] e θ T n E [ u ] 2 e θ T n ) + e θ T n 1 ) ]
The retailer’s deteriorating cost per year is:
d d n T ( Q u Q D T n )   =   d d ( ( 1 u ) n T D ( e θ T / n 1 ) θ ( 1 u e θ T / n ) D )
The retailer‘s total cost is the sum of the ordering, inspection, deteriorating, inventory holding, and emission costs. Therefore, considering the probability of the defective products, the expected total cost per year is:
E T C d = c T + n T ( i c + u c D ( e θ T n 1 ) θ ( 1 E [ u ] e θ T n ) ) + d d ( ( 1 E [ u ] ) n T D ( e θ T / n 1 ) θ ( 1 E [ u ] e θ T / n ) D ) + ( h d + w e ) n D T θ 2 ( E [ u ] e θ T n 1 ) 2 ( E [ u ] D x ( 1 + e 2 θ T n 2 e θ T n ) + E [ u ] ( 1 e 2 θ T n + E [ u ] e 2 θ T n E [ u ] e θ T n ) + θ T n ( E [ u ] 1 + E [ u ] e θ T n E [ u ] 2 e θ T n ) + e θ T n 1 )

4.1.2. Manufacturer Cost and Emission

After the arrival of the retailer’s order, the manufacturer starts the production of nQ units of the item at a production rate P. Since s is the setup cost per production cycle, the manufacturer’s setup cost per year is s/T.
The first delivery occurs as soon as the quantity is met. The following deliveries occur at T/n intervals. The transportation cost belongs to the manufacturer and consist of fixed and variable costs. Swenseth and Godfrey [54], Nie et al. [55], Rahman et al. [56], and Wangsa and Wee [57] incorporated the variable transportation cost, which is affected by the shipping distance and truckloads. The manufacturer’s transportation cost per delivery is given by:
t f + 2 d c 1 t v + d Q w c 2 t v
The first element is the fixed transportation setup cost. The second element calculates the transportation cost of an empty truck. As the truck goes from the manufacturer to the retailer and then goes back, the distance is multiplied by two. Then, the transportation cost for the truckload is calculated, which depends on the delivery distance and quantity, product weight, additional fuel consumption per ton per km, and the fuel price. Substituting Equation (4) to (10), the manufacturer’s transportation cost per year is given by:
n T ( t f + 2 d c 1 t v + d D ( e θ T / n 1 ) θ ( 1 u e θ T / n ) w c 2 t v )
Wahab et al. [10] identified that the emissions from transportation were affected by the delivery distance, actual shipment weight, fuel consumption per km, and CO2 emissions per liter of fuel. Therefore, the amount of the manufacturer’s carbon emission per year as the result of transportation activity can be derived as follows:
n T ( 2 d c 1 + d Q w c 2 ) F e = n T ( 2 d c 1 + d D ( e θ T / n 1 ) θ ( 1 u e θ T / n ) w c 2 ) F e
As depicted in Figure 6, Lee and Kim [22] studied a similar inventory model for the manufacturer using Yang and Wee’s [52] approach. Both production and consumption occur during T1, while only consumption occurs during T2. Hsu [58] suggested some revision considering the effect of defective products. Due to defective products and the retailer’s quality inspection, the effective demand rate for the manufacturer becomes D/(1-u). Therefore, the inventory functions are as follows:
I p 1 ( t 1 ) = P ( D / ( 1 u ) ) θ ( 1 e θ   t 1 ) ,   0   t 1   T 1
I p 2 ( t 2 ) = ( D / ( 1 u ) ) θ ( e θ   ( T 2 t 2 ) 1 ) , 0     t 2   T 2
From the boundary condition Ip1(T1) = Ip2(0) and following Misra’s [59] approximation:
( P D 1 u ) T 1 ( 1 1 2 θ T 1 ) = ( D 1 u ) T 2 ( 1 + 1 2 θ T 2 )
Hence, one has:
T 1 D ( 1 u ) P D T 2 ( 1 + 1 2 θ T 2 )
T T 2 ( 1 u ) P D ( ( 1 u ) P + 1 2 D θ T 2 )
From Yang and Wee [52], the manufacturer’s inventory per cycle is:
0 T 1 I p 1 ( t 1 ) d t 1 + 0 T 2 I p 2 ( t 2 ) d t 2 n 0 T / n I d ( t ) d t
Hence, the manufacturer’s holding cost per year is:
h p T   [ P ( D / ( 1 u ) ) θ T 1 + P ( D / ( 1 u ) ) θ 2 ( e θ   T 1 1 ) ( D / ( 1 u ) ) T 2 θ ( D / ( 1 u ) ) θ 2 ( 1 e θ   T 2 ) n ( D θ 2 ( u e θ T n 1 ) 2 ( u D x ( 1 + e 2 θ T n 2 e θ T n ) + u ( 1 e 2 θ T n + u e 2 θ T n u e θ T n ) + θ T n ( u 1 + u e θ T n u 2 e θ T n ) + e θ T n 1 ) ) ]
The amount of a manufacturer’s carbon emissions per year as the result of warehousing activity can be derived as follows:
e c E e T [ P ( D / ( 1 u ) ) θ T 1 + P ( D / ( 1 u ) ) θ 2 ( e θ   T 1 1 ) ( D / ( 1 u ) ) T 2 θ ( D / ( 1 u ) ) θ 2 ( 1 e θ   T 2 ) n ( D θ 2 ( u e θ T n 1 ) 2 ( u D x ( 1 + e 2 θ T n 2 e θ T n ) + u ( 1 e 2 θ T n + u e 2 θ T n u e θ T n ) + θ T n ( u 1 + u e θ T n u 2 e θ T n ) + e θ T n 1 ) ) ]
The manufacturer’s carbon emissions per year come from Equations (12) and (19). Therefore, the manufacturer’s carbon emission cost is:
n T ( 2 d e 1 + d D ( e θ T / n 1 ) θ ( 1 u e θ T / n ) e 2 ) + w e T [ P ( D / ( 1 u ) ) θ T 1 + P ( D / ( 1 u ) ) θ 2 ( e θ   T 1 1 )   ( D / ( 1 u ) ) T 2 θ ( D / ( 1 u ) ) θ 2 ( 1 e θ   T 2 ) n ( D θ 2 ( u e θ T n 1 ) 2 ( u D x ( 1 + e 2 θ T n 2 e θ T n ) + u ( 1 e 2 θ T n + u e 2 θ T n u e θ T n ) + θ T n ( u 1 + u e θ T n u 2 e θ T n ) + e θ T n 1 ) ) ]
The loss due to deterioration in the manufacturer’s inventory is the total production during the period T1, minus the total delivered products to the retailer’s inventory. Therefore, the manufacturer’s deteriorating cost per year is:
d p T ( P T 1 n ( D ( e θ T / n 1 ) θ ( 1 u e θ T / n ) ) )
From Equations (11), (18), (20), and (21), as well as the setup cost, and considering the probability of the defective products, the manufacturer‘s expected total cost per year is:
E T C p = s T + n T ( t f + 2 d c 1 t v + d D ( e θ T / n 1 ) θ ( 1 E [ u ] e θ T / n ) w c 2 t v ) + ( h p + w e ) T [ P ( D / ( 1 E [ u ] ) ) θ T 1 + P ( D / ( 1 E [ u ] ) ) θ 2 ( e θ   T 1 1 )   ( D / ( 1 E [ u ] ) ) T 2 θ ( D / ( 1 E [ u ] ) ) θ 2 ( 1 e θ   T 2 ) n ( D θ 2 ( E [ u ] e θ T n 1 ) 2 ( e θ T n 1 + E [ u ] D x ( 1 + e 2 θ T n 2 e θ T n ) + E [ u ] ( 1 e 2 θ T n + E [ u ] e 2 θ T n E [ u ] e θ T n ) + θ T n ( E [ u ] 1 + E [ u ] e θ T n E [ u ] 2 e θ T n ) ) ] + n T ( 2 d e 1 + d D ( e θ T / n 1 ) θ ( 1 E [ u ] e θ T / n ) e 2 ) +   d p T ( P T 1 n ( D ( e θ T / n 1 ) θ ( 1 u e θ T / n ) ) )

4.1.3. The Integrated Manufacturer and Retailer Cost Function

In an integrated decision, the manufacturer and the retailer jointly specify n, which minimizes the expected total cost (ETC). The ETC is the sum of ETCd and ETCp in Equations (9) and (22).
Using Taylor’s series expansion for a small value of θT/n, θT1, and θT2, we can solve the cost function by assuming ex as 1 + x + x2/2 + x3/6. Furthermore, the expected total emissions (ETE) per year from the manufacturer and the retailer can be derived from Equations (7), (12), and (19).

4.1.4. Methodology and Solution Search

The objective is to determine the optimal number of deliveries (n*) that minimize the expected total cost function (ETC). The value of n* and the respective T, T1, and T2 will lead us to the optimal delivery quantity Q* and production quantity R. The proposed procedure to derive the positive integer decision variable n is adapted from Tiwari et al. [15] and Yang and Wee [52] as follows:
Step 1.
Substitute the T1 and T functions in Equations (15) and (16) into ETC;
Step 2.
Input all the known parameters;
Step 3.
Set n = 1;
Step 4.
Derive the partial derivative of ETC with respect to T2 and set it to zero. Solve the equation to find the value of T2;
Step 5.
Use the known n and T2 to find the value of T1 and T using Equations (15) and (16).
Step 6.
Derive the corresponding ETC;
Step 7.
If ETC(n) > ETC(n − 1), then n* = n − 1 and go to Step 8; otherwise, set n = n + 1 and go back to Step 4;
Step 8.
Use n* and the corresponding T* to find Q* from Equation (4) and calculate R = PT1*.

4.2. Model Development with Manufacturer Inspection

This sub-section presents the model development when quality inspection becomes the manufacturer’s responsibility. The purpose of such a policy is to prevent transporting defective products. The manufacturer performs a quality inspection of all the produced products and keeps the defective products separately until the end of the production period T1. The defective products will be sold at a discounted price to the secondary market. The inventory level for both the manufacturer (Ip) and the retailer (Id), including the manufacturer’s inventory of the defective products (Ipd), is illustrated in Figure 7. In this model, the Ipd is accumulated during the T1 period. Besides, Id decreases solely due to demand.

4.2.1. Retailer Cost and Emission

The retailer‘s total cost is the sum of the ordering, inventory holding, deteriorating, and emission costs. The ordering cost per year is given by c/T, which is similar to Model 1 (Section 4.1.1). During the T/n period, Id decreases due to demand and deterioration. From Yang and Wee [52], the inventory function is as follows:
I d ( t ) = D θ ( e θ ( T n   t )   ) , 0     t   T / n
and
Q = I d ( 0 ) = D θ ( e θ T n 1 )
and the holding cost per year is:
h d n T ( 0 T / n I d ( t ) d t ) = h d n T ( D θ ( 1 θ ( e θ T n 1 ) T n ) )
Further, the deteriorating cost per year is:
d d n T ( Q D T n ) = d d n T ( D θ ( e θ T n 1 ) D T n )
E T E d = e c E e ( n T D θ ( 1 θ ( e θ T n 1 ) T n ) )
Equation (27) shows the retailer’s carbon emission (ETEd) from holding the inventory. Therefore, the retailer’s carbon emission cost per year is:
w e ( n T D θ ( 1 θ ( e θ T n 1 ) T n ) )
The retailer‘s expected total cost per year (ETCd) becomes:
E T C d = c T + ( h d + w e ) n T ( D θ ( 1 θ ( e θ T n 1 ) T n ) ) + d d n T ( D θ ( e θ T n 1 ) D T n )

4.2.2. Manufacturer Cost and Emission

Due to some percentage of defective products, the production rate of the perfect product is (1 − u)P. The manufacturer’s setup cost per year is s/T. In this model, the manufacturer will have an additional inspection cost. Since the total number of products being produced per production cycle is PT1, considering a fixed inspection cost per cycle (ic) and unit inspection cost (uc), the manufacturer’s inspection cost per year is:
i c T + u c P T 1 T
The manufacturer’s transportation function is similar to Equation (10); therefore, the transportation cost and emissions per year become:
n T ( t f + 2 d c 1 t v + d D θ ( e θ T n 1 ) w c 2 t v )
n T ( 2 d c 1 + d D θ ( e θ T n 1 ) w c 2 ) F e
From Figure 6, the inventory differential equations are:
d I p 1 ( t 1 ) = ( ( 1 u ) P D ) d t 1 θ I p 1 ( t 1 ) d t 1 ,   0   t 1   T 1
d I p 2 ( t 2 ) = D d t 2 θ I p 2 ( t 2 ) d t 2 ,   0   t 2   T 2
For the boundary condition for t1 = 0, I1 (0) = 0 and for t2 = 0, I2 (0) = Io and for t2 = T2, I2 (T2) = 0, the manufacturer’s inventory functions for the good products are:
I p 1 ( t 1 ) = ( 1 u ) P D θ ( 1 e θ   t 1 ) ,   0   t 1   T 1
I p 2 ( t 2 ) = D θ ( e θ   ( T 2 t 2 ) 1 ) , 0     t 2   T 2
From the boundary condition Ip1(T1) = Ip2(0), we have the following equation:
( ( 1 u ) P D ) θ ( 1 e θ   T 1 ) = D θ ( e θ T 2 1 )
From Taylor’s series expansion and the assumption of 𝜃T << 1, following Misra’s [59] approximation, one has:
( ( 1 u ) P D ) T 1 ( 1 1 2 θ T 1 ) = D T 2 ( 1 + 1 2 θ T 2 )
T 1 D ( 1 u ) P D T 2 ( 1 + 1 2 θ T 2 )
T T 2 ( 1 u ) P D ( ( 1 u ) P + 1 2 D θ T 2 )
Therefore, the manufacturer’s inventory for good products becomes:
0 T 1 ( 1 u ) P D θ ( 1 e θ   t 1 ) d t 1 + 0 T 2 D θ ( e θ   ( T 2 t 2 ) 1 ) d t 2 n [ D θ ( 1 θ ( e θ T n 1 ) T n ) ]
Besides, there is an inventory of defective products. From Figure 6, the inventory differential equation for the defective products is:
d I p d ( t 1 ) = u P d t 1 θ I p d ( t 1 ) d t 1 ,   0   t 1   T 1
For the boundary condition for t1 = 0, I1(0) = 0, the manufacturer’s inventory function for the defective products is:
I p d ( t 1 ) = u P θ ( 1 e θ   t 1 ) ,   0   t 1   T 1
Therefore, the manufacturer’s inventory of the defective products becomes:
0 T 1 u P θ ( 1 e θ   t 1 ) d t 1
Hence, the manufacturer’s holding cost per year is:
h p T ( ( 1 u ) P D θ T 1 + ( 1 u ) P D θ 2 ( e θ T 1 1 ) D T 2 θ D θ 2 ( 1 e θ T 2 ) n ( D θ ( 1 θ ( e θ T n 1 ) T n ) ) + u P T 1 θ + u P θ 2 ( e θ T 1 1 ) )
Therefore, based on Equations (32) and (40), the manufacturer’s carbon emission cost and the total expected carbon emissions per year can be calculated as follows:
n T ( 2 d e 1 + d D θ ( e θ T n 1 ) e 2 )    + w e T ( ( 1 u ) P D θ T 1 + ( 1 u ) P D θ 2 ( e θ T 1 1 ) D T 2 θ    − D θ 2 ( 1 e θ T 2 ) n ( D θ ( 1 θ ( e θ T n 1 ) T n ) ) + u P T 1 θ    + u P θ 2 ( e θ T 1 1 ) )
E T E p = n T ( 2 d c 1 + d D θ ( e θ T n 1 ) w c 2 ) F e + e c E e T ( ( 1 u ) P D θ T 1 + ( 1 u ) P D θ 2 ( e θ T 1 1 ) D T 2 θ D θ 2 ( 1 e θ T 2 ) n ( D θ ( 1 θ ( e θ T n 1 ) T n ) ) + u P T 1 θ + u P θ 2 ( e θ T 1 1 ) )
The number of deteriorated items in the manufacturer’s inventory is the total production during the period T1, minus the total products delivered to the buyer and the inventory of the defective products. Therefore, the manufacturer’s deteriorating cost per year is:
d p T ( ( 1 u ) P T 1 n ( D θ ( e θ T n 1 ) ) + ( u P T 1 u P θ ( 1 e θ   T 1 ) ) )
Considering the additional inspection cost and the probability of the defective products, the manufacturer‘s expected total cost per year is:
E T C p = s T + i c T + u c P T 1 T + n T ( t f + 2 d c 1 t v + d D θ ( e θ T n 1 ) w c 2 t v ) + n T ( 2 d e 1 + d D θ ( e θ T n 1 ) e 2 ) + ( h p + w e T ( ( 1 E [ u ] ) P D θ T 1 + ( 1 E [ u ] ) P D θ 2 ( e θ T 1 1 ) D T 2 θ D θ 2 ( 1 e θ T 2 ) n ( D θ ( 1 θ ( e θ T n 1 ) T n ) ) + E [ u ] P T 1 θ + E [ u ] P θ 2 ( e θ T 1 1 ) ) + d p T ( ( 1 E [ u ] ) P T 1 n ( D θ ( e θ T n 1 ) ) + ( E [ u ] P T 1 E [ u ] P θ ( 1 e θ T 1 ) ) )

4.2.3. The Integrated Manufacturer and Retailer Cost Function

The ETC of the integrated system is the sum of Equations (29) and (44). Using Taylor’s series expansion for a small value of 𝜃T/n, 𝜃T1, and 𝜃T2, we can solve the cost function by assuming ex as 1 + x + x2/2 + x3/6. Furthermore, the ETE can be derived from Equation (27) and Equation (42).

4.2.4. Methodology and Solution Search

Similar to Model 1, the objective is to determine the optimal number of deliveries (n*) that minimize the expected total cost function ETC. The proposed procedure to search for the optimum solution is as follows:
Step 1.
Substitute the T1 and T functions in Equations (36) and (37) into ETC;
Step 2.
Input all the known parameters;
Step 3.
Set n = 1;
Step 4.
Derive the partial derivative of ETC with respect to T2 and set it to zero. Solve the equation to find the value of T2;
Step 5.
Use the known n and T2 to find the value of T1 and T using Equations (36) and (37).
Step 6.
Derive the corresponding ETC;
Step 7.
If ETC(n) > ETC(n-1) then n* = n − 1 and go to step 8, otherwise set n = n + 1 and back to Step 4;
Step 8.
Use n* and the corresponding T* to find Q* from Equation (24) and calculate R = PT1*.

5. Numerical Example and Management Insights

5.1. Numerical Example 1

The values of the parameters are considered by adopting data from Yang and Wee [52], Hariga et al. [30], and Tiwari et al. [15] as P = 2,000,000 units/year, D = 500,000 units/year, x = 1,725,000 unit/year, ic = $500/delivery, uc = $0.5/unit, c = $2,000/order, s = $100,000/setup, hd = $60/unit/year, hp = $40/unit/year, dd = $600/unit, dp = $400/unit, 𝜃 = 0.1, d = 100 km, tf = $1000/delivery, tv = $0.75/liter, w = 0.01 ton/unit, c1 = 27 L/100 km, c2 = 0.57 L/100 km/ton truckload, ec = 1.44 kWh/unit/year, Tx = $75/tonCO2, Fe = 2.6 × 10−3 tonCO2/L (US. EPA [60]), Ee = 0.5 × 10−3 tonCO2/kWh (McCarthy [61]), and u is uniformly distributed in which α = 0 and β = 0.04, with E[u] = 0.02.
The minimum value of joint expected total cost can be obtained at n* = 7 with T2 = 0.0651856, T1 = 0.0223966, and T = 0.0875822, as shown in Table 3. The ETC is $2,834,922/year, which is from Equation (4), the optimum Q is 6,387.7 units. The optimum R is 44,793.2 units, and the ETE is 30.598 tonCO2/year. If the supply chain solely minimizes the total amount of carbon footprint, the decision is to perform a single-setup single-delivery (SSSD) as n = 1 with ETE = 18.460 tonCO2/year (saving 39.7%). However, this situation increases the ETC into $3,366,391 (18.7%). Figure 8 shows the convexity of ETC when n = 7.
When the probability of defective products, unit inspection costs, carbon tax, and variable transport cost are equal to zero (E[u] = uc = Tx = tv = 0), the results are n = 7, T = 0.08791, and ETC = $2,559,246 which are similar to the results of Yang and Wee [52].

5.2. Numerical Example 2

We consider the parameters in numerical example one and solve it using the Model 2 results of the following values. The minimum value of joint expected total cost can be obtained at n* = 9 and T = 0.08869, as shown in Table 4. The Q, R, and ETC are 4,929.6 units, 45,360.7 units, and $2,782,396/year, respectively. The ETE is 33.52 tonCO2/year.
Table 5 provides the cost comparison between the two models for the result of examples one and two. The number of deliveries per cycle (n) is higher in Model 2 (when the manufacturer performs the quality inspection), while the delivery lot size (Q) is lower. All the retailer’s cost components decrease, while all the manufacturer’s cost components increase, except for the setup cost. In total, the ETC of Model 2 is 1.85% lower than the ETC in Model 1. However, the retailer’s total costs were reduced by 54.7%, while the manufacturer’s total costs increased by 15.6%. Considering this situation, cost-saving compensation from the retailer to the manufacturer is an alternative solution so that both parties take advantage of the implementation of the second inspection policy. Based on Goyal [62]:
z = E T C d c a s e 1 E T C c a s e 1
z is the retailer’s cost coefficient. Therefore, ETCd and ETCp after cost-saving compensation are
E T C d a = z E T C c a s e 2
E T C p a = ( 1 z ) E T C c a s e 2
Hence, the retailer and manufacturer total cost per year become $690,574.4 and $2,091,821.6, respectively. Finally, by using this compensation policy, the cost decreases $13,036.7 for the retailer and $39,489.6 for the vendor, or 1.85% for both parties. Table 5 also shows that the ETE of Model 2 is 33.52 tonCO2/year, which is 9.55% higher than the ETE in Model 1. We can obtain both the cost-saving and emissions-reducing objectives from Table 4, as there is a chance to reduce the ETE of Model 2. For n = 7, the ETE is 30.292, and the ETC is $2,781,779 in which now the ETE and the ETC are 1.0% and 1.87% lower than those of Model 1, respectively. Thus, the objectives of cost efficiency and carbon footprint level reduction can be obtained simultaneously. The new comparison between the two models is presented in Table 6. It is also observed that the total delivered products to the retailer in Model 2 after adjustment (Model 2 adj) is less than those in Model 1.
Sensitivity analysis is performed by increasing or decreasing the value of the parameter by ±25% and ±50% from the original values, as shown in Table 7. The results confirm that the second model is superior to the first model in terms of total cost. The number of deliveries per cycle (n) is sensitive to changes in parameters P, D, s, hd, hp, dd, dp, and tf. As the values of parameters P, hp, dp, and tf increase, the smaller the value of n. Contradictory conditions occur for parameters D, s, hd, and dd. The expected total cost is highly sensitive to the changes in parameters P, D, θ, s, uc, hd, hp, dd, dp, and tf, and almost insensitive to the changes in other parameters.
It is observed that when the deterioration rate (θ) increases, the expected total cost and the number of deliveries increase, but the delivery quantity decreases. When the probability of defective products (u) increases, the expected total cost increases very slightly, especially when the inspection is performed by the vendor. When the carbon tax (Tx) increases, the number of deliveries remains stable. Otherwise, the delivery quantity and expected total cost increase are very small.

6. Conclusions and Future Research

This study considers a two-echelon supply chain consisting of a manufacturer and a retailer where the production activities are resulting in a certain percentage of defective products. The supply chain entities are willing to reduce their environmental impact by coordinating the delivery quantity and number of deliveries per cycle. The effect of carbon emissions, item deterioration, and two choices of inspection are examined. The models are illustrated with two numerical examples, and the results give some insights. This study is an initial exploratory study that attempts to provide a mathematical solution for a controlled situation; it may be applied to handle larger problems of cost minimization and carbon emission reduction in the future.
From the research finding, it is observed that the numbers of delivered products from the manufacturer are less when the inspection is performed by the vendor. As a result, the total cost of the supply chain is less, because the total inventory-holding cost and the total deteriorating cost are decreasing. However, the vendor’s total cost becomes higher when it performs the inspection. Therefore, the retailer needs to compensate a certain amount of cost-saving to the manufacturer so that both parties take advantage.
The research finding also revealed that although the total cost is less when the inspection is performed by the vendor, it does not guarantee a reduction in emissions. However, both the cost-saving and emission-reducing objectives can still be obtained simultaneously by reducing the level of cost savings. In this situation, there is a tradeoff between cost savings and reduction in carbon emissions.
Although this study addresses some practical aspects of a supply chain scenario to deal with lower carbon emissions, the scope has a wide opportunity to be extended. Applying the approach in a three-echelon supply chain or more is one opportunity. Future works can consider the possibility of reworking the defective products, the capacity of the vehicle and storage facility, and investment to reduce the carbon emissions. This study assumes a 100% inspection by the manufacturer. Future research may assume a sampling inspection by the manufacturer, as well as incorporate the issue of imperfect quality inspection.

Author Contributions

Conceptualization, Y.D.; Methodology, Y.D.; Software, Y.D.; Validation, G.A.W. and H.M.W.; Formal Analysis, Y.D.; Investigation, Y.D.; Resources, H.M.W.; Writing—Original Draft Preparation, Y.D.; Writing—Review & Editing, H.M.W. and G.A.W.; Visualization, Y.D.; Supervision, H.M.W.

Funding

This research received no external funding.

Acknowledgments

The authors express their gratitude to the editor and the two anonymous reviewers for their comments and valuable suggestions to improve this paper. The first author also thanks to the United Board for Christian Higher Education in Asia for financing his study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Glock, C.H. The joint economic lot size problem: A review. Int. J. Prod. Econ. 2012, 135, 671–686. [Google Scholar] [CrossRef]
  2. Luo, Z.; Gunasekaran, A.; Dubey, R.; Childe, S.J.; Papadopoulos, T. Antecedents of low carbon emissions supply chains. Int. J. Clim. Chang. Strateg. Manag. 2017, 9, 707–727. [Google Scholar] [CrossRef] [Green Version]
  3. Das, C.; Jharkharia, S. Low carbon supply chain: A state-of-the-art literature review. J. Manuf. Technol. Manag. 2018, 29, 398–428. [Google Scholar] [CrossRef]
  4. Kazemi, N.; Abdul-Rashid, S.H.; Ghazilla, R.A.R.; Shekarian, E.; Zanoni, S. Economic order quantity models for items with imperfect quality and emission considerations. Int. J. Syst. Sci. Oper. Logist. 2018, 5, 99–115. [Google Scholar] [CrossRef]
  5. Sarkar, B.; Ahmed, W.; Choi, S.B.; Tayyab, M. Sustainable inventory management for environmental impact through partial backordering and multi-trade-credit period. Sustainability 2018, 10, 4761. [Google Scholar] [CrossRef]
  6. Taleizadeh, A.A.; Soleymanfar, V.R.; Govindan, K. Sustainable economic production quantity models for inventory system with shortage. J. Clean. Prod. 2018, 174, 1011–1020. [Google Scholar] [CrossRef]
  7. Sarkar, B.; Ganguly, B.; Sarkar, M.; Pareek, S. Effect of variable transportation and carbon emission in a three-echelon supply chain model. Transp. Res. Part E Logist. Transp. Rev. 2016, 91, 112–128. [Google Scholar] [CrossRef]
  8. Sarkar, B.; Ahmed, W.; Kim, N. Joint effects of variable carbon emission cost and multi-delay-in-payments under single-setup-multiple-delivery policy in a global sustainable supply chain. J. Clean. Prod. 2018, 185, 421–445. [Google Scholar] [CrossRef]
  9. Daryanto, Y.; Wee, H.M. Single vendor-buyer integrated inventory model for deteriorating items considering carbon emission. In Proceedings of the 8th International Conference on Industrial Engineering and Operations Management (IEOM), Bandung, Indonesia, 6–8 March 2018; pp. 544–555. [Google Scholar]
  10. Wahab, M.I.M.; Mamun, S.M.H.; Ongkunarak, P. EOQ models for a coordinated two-level international supply chain considering imperfect items and environmental impact. Int. J. Prod. Econ. 2011, 134, 151–158. [Google Scholar] [CrossRef]
  11. Jauhari, W.A.; Pamuji, A.S.; Rosyidi, C.N. Cooperative inventory model for vendor-buyer system with unequal-sized shipment, defective items and carbon emission cost. Int. J. Logist. Syst. Manag. 2014, 19, 163–186. [Google Scholar] [CrossRef]
  12. Sarkar, B.; Saren, S.; Sarkar, M.; Seo, Y.W. A Stackelberg game approach in an integrated inventory model with carbon-emission and setup cost reduction. Sustainability 2016, 8, 1244. [Google Scholar] [CrossRef]
  13. Jauhari, W.A. A collaborative inventory model for vendor-buyer system with stochastic demand, defective items and carbon emission cost. Int. J. Logist. Syst. Manag. 2018, 29, 241–269. [Google Scholar] [CrossRef]
  14. Gautam, P.; Khanna, A. An imperfect production inventory model with setup cost reduction and carbon emission for an integrated supply chain. Uncertain Supply Chain Manag. 2018, 6, 271–286. [Google Scholar] [CrossRef]
  15. Tiwari, S.; Daryanto, Y.; Wee, H.M. Sustainable inventory management with deteriorating and imperfect quality items considering carbon emissions. J. Clean. Prod. 2018, 192, 281–292. [Google Scholar] [CrossRef]
  16. Rosenblatt, M.J.; Lee, H.L. Economic production cycles with imperfect production processes. IIE Trans. 1986, 18, 48–55. [Google Scholar] [CrossRef]
  17. Porteus, E.L. Optimal lot sizing, process quality improvement and setup cost reduction. Oper. Res. 1986, 34, 137–144. [Google Scholar] [CrossRef]
  18. Salameh, M.K.; Jaber, M.Y. Economic production quantity model for items with imperfect quality. Int. J. Prod. Econ. 2000, 64, 59–64. [Google Scholar] [CrossRef]
  19. Huang, C.K. An integrated vendor-buyer cooperative inventory model for items with imperfect quality. Prod. Plan. Control 2002, 13, 355–361. [Google Scholar] [CrossRef]
  20. Goyal, S.K.; Huang, C.K.; Chen, K.C. A simple integrated production policy of an imperfect item for vendor and buyer. Prod. Plan. Control 2003, 14, 596–602. [Google Scholar] [CrossRef]
  21. Wee, H.M.; Yu, J.C.P.; Wang, K.J. An integrated production-inventory model for deteriorating items with imperfect quality and shortage backordering considerations. In Proceedings of the International Conference on Computational Science and Its Applications (ICCSA), Glasgow, UK, 8–11 May 2006; Springer: Berlin/Heidelberg, Germany, 2006; pp. 885–897. [Google Scholar]
  22. Lee, S.; Kim, D. An optimal policy for a single vendor single-buyer integrated production-distribution model with both deteriorating and defective items. Int. J. Prod. Econ. 2014, 147, 161–170. [Google Scholar] [CrossRef]
  23. Bazan, E.; Jaber, M.Y.; Zanoni, S.; Zavanella, L.E. Vendor managed inventory (VMI) with consignment stock (CS) agreement for a two-level supply chain with an imperfect production process with/without restoration interruptions. Int. J. Prod. Econ. 2014, 157, 289–301. [Google Scholar] [CrossRef]
  24. Sarkar, B.; Shaw, B.K.; Kim, T.; Sarkar, M.; Shin, D. An integrated inventory model with variable transportation cost, two-stage inspection, and defective items. J. Ind. Manag. Optim. 2017, 13, 1975–1990. [Google Scholar] [CrossRef]
  25. Yu, H.F.; Hsu, W.K. An integrated inventory model with immediate return for defective items under unequal-sized shipments. J. Ind. Prod. Eng. 2017, 34, 70–77. [Google Scholar] [CrossRef]
  26. Benjaafar, S.; Li, Y.; Daskin, M. Carbon footprint and the management of supply chains: Insights from simple models. IEEE Trans. Autom. Sci. Eng. 2013, 10, 99–116. [Google Scholar] [CrossRef]
  27. Fahimnia, B.; Sarkis, J.; Dehghanian, F.; Banihashemi, N.; Rahman, S. The impact of carbon pricing on a closed-loop supply chain: An Australian case study. J. Clean. Prod. 2013, 59, 210–225. [Google Scholar] [CrossRef]
  28. Bozorgi, A.; Pazour, J.; Nazzal, D. A new inventory model for cold items that considers costs and emissions. Int. J. Prod. Econ. 2014, 155, 114–125. [Google Scholar] [CrossRef]
  29. Bozorgi, A. Multi-product inventory model for cold items with cost and emission consideration. Int. J. Prod. Econ. 2016, 176, 123–142. [Google Scholar] [CrossRef]
  30. Hariga, M.; As’ad, R.; Shamayleh, A. Integrated economic and environmental models for a multi stage cold supply chain under carbon tax regulation. J. Clean. Prod. 2017, 166, 1357–1371. [Google Scholar] [CrossRef]
  31. Ghosh, A.; Jha, J.K.; Sarmah, S.P. Optimizing a two-echelon serial supply chain with different carbon policies. Int. J. Sustain. Eng. 2016, 9, 363–377. [Google Scholar] [CrossRef]
  32. Toptal, A.; Çetinkaya, B. How supply chain coordination affects the environment: A carbon footprint perspective. Ann. Oper. Res. 2017, 250, 487–519. [Google Scholar] [CrossRef]
  33. Bouchery, Y.; Ghaffari, A.; Jemai, Z.; Tan, T. Impact of coordination on cost and carbon emissions for a two-echelon serial economic order quantity problem. Eur. J. Oper. Res. 2017, 260, 520–533. [Google Scholar] [CrossRef]
  34. Dwicahyani, A.R.; Jauhari, W.A.; Rosyidi, C.N.; Laksono, P.W. Inventory decisions in a two-echelon system with remanufacturing, carbon emission, and energy effects. Cogent Eng. 2017, 4, 1–17. [Google Scholar] [CrossRef]
  35. Li, J.; Su, Q.; Ma, L. Production and transportation outsourcing decisions in the supply chain under single and multiple carbon policies. J. Clean. Prod. 2017, 141, 1109–1122. [Google Scholar] [CrossRef]
  36. Wangsa, I.D. Greenhouse gas penalty and incentive policies for a joint economic lot size model with industrial and transport emissions. Int. J. Ind. Eng. Comput. 2017, 8, 453–480. [Google Scholar]
  37. Anvar, S.H.; Sadegheih, A.; Zad, M.A.V. Carbon emission management for greening supply chains at the operational level. Environ. Eng. Manag. J. 2018, 17, 1337–1347. [Google Scholar]
  38. Hariga, M.; Babekian, S.; Bahroun, Z. Operational and environmental decisions for a two-stage supply chain under vendor managed consignment inventory partnership. Int. J. Prod. Res. 2018. [Google Scholar] [CrossRef]
  39. Ji, S.; Zhao, D.; Peng, X. Joint decisions on emission reduction and inventory replenishment with overconfidence and low-carbon preference. Sustainability 2018, 10, 1119. [Google Scholar]
  40. Wang, S.; Ye, B. A comparison between just-in-time and economic order quantity models with carbon emissions. J. Clean. Prod. 2018, 187, 662–671. [Google Scholar] [CrossRef]
  41. Ghosh, A.; Sarmah, S.P.; Jha, J.K. Collaborative model for a two-echelon supply chain with uncertain demand under carbon tax policy. Sādhanā 2018, 43, 144. [Google Scholar] [CrossRef]
  42. Ma, X.; Ho, W.; Ji, P.; Talluri, S. Coordinated pricing analysis with the carbon tax scheme in a supply chain. Decis. Sci. 2018, 49, 863–900. [Google Scholar] [CrossRef]
  43. Darom, N.A.; Hishamuddin, H.; Ramli, R.; Nopiah, Z.M. An inventory model of supply chain disruption recovery with safety stock and carbon emission consideration. J. Clean. Prod. 2018, 197, 1011–1021. [Google Scholar] [CrossRef]
  44. Huang, H.; He, Y.; Li, D. Pricing and inventory decisions in the food supply chain with production disruption and controllable deterioration. J. Clean. Prod. 2018, 180, 280–296. [Google Scholar] [CrossRef]
  45. Daryanto, Y.; Wee, H.M.; Astanti, R.D. Three-echelon supply chain model considering carbon emission and item deterioration. Transp. Res. Part E Logist. Transp. Rev. 2019, 122, 368–383. [Google Scholar] [CrossRef]
  46. Kundu, S.; Chakrabarti, T. A fuzzy rough integrated multi-stage supply chain inventory model with carbon emissions under inflation and time-value of money. Int. J. Math. Oper. Res. 2019, 14, 123–145. [Google Scholar] [CrossRef]
  47. Shalke, P.N.; Paydar, M.M.; Hajiaghaei-Keshteli, M. Sustainable supplier selection and order allocation through quantity discounts. Int. J. Manag. Sci. Eng. Manag. 2018, 13, 20–32. [Google Scholar]
  48. Moheb-Alizadeh, H.; Handfield, R. An integrated chance-constrained stochastic model for efficient and sustainable supplier selection and order allocation. Int. J. Prod. Res. 2018, 56, 6890–6916. [Google Scholar] [CrossRef]
  49. Moheb-Alizadeh, H.; Handfield, R. Sustainable supplier selection and order allocation: A novel multi-objective programming model with a hybrid solution approach. Comput. Ind. Eng. 2019, 129, 192–209. [Google Scholar] [CrossRef]
  50. Cao, W.; Hu, Y.; Li, C.; Wang, X. Single setup multiple delivery model of JIT system. Int. J. Adv. Manuf. Technol. 2007, 33, 1222–1228. [Google Scholar] [CrossRef]
  51. Volvo Truck Corporation. Emissions from Volvo’s Truck. Issue 3. 09 March 2018. Available online: http://www.volvotrucks.com/content/dam/volvo/volvo-trucks/markets/global/pdf/our-trucks/Emis_eng_10110_14001.pdf (accessed on 25 September 2018).
  52. Yang, P.C.; Wee, H.M. Economic ordering policy of deteriorated item for vendor and buyer: An integrated approach. Prod. Plan. Control 2000, 11, 474–480. [Google Scholar] [CrossRef]
  53. Jaggi, C.K.; Goel, S.K.; Mittal, M. Economic order quantity model for deteriorating items with imperfect quality and permissible delay in payment. Int. J. Ind. Eng. Comput. 2011, 2, 237–248. [Google Scholar] [CrossRef]
  54. Swenseth, S.R.; Godfrey, M.R. Incorporating transportation costs into inventory replenishment decisions. Int. J. Prod. Econo. 2002, 77, 113–130. [Google Scholar] [CrossRef]
  55. Nie, L.; Xu, X.; Zhan, D. Incorporating transportation costs into JIT lot splitting decisions for coordinated supply chains. J. Adv. Manuf. Syst. 2006, 5, 111–121. [Google Scholar] [CrossRef]
  56. Rahman, M.N.A.; Leuveano, R.A.C.; bin Jafar, F.A.; Saleh, C.; Deros, B.M.; Mahmood, W.M.F.W.; Mahmood, W.H.W. Incorporating logistic costs into a single vendor-buyer JELS model. Appl. Math. Model. 2016, 40, 10809–10819. [Google Scholar] [CrossRef]
  57. Wangsa, I.D.; Wee, H.M. An integrated vendor-buyer inventory model with transportation cost and stochastic demand. Int. J. Syst. Sci. Oper. Logist. 2018, 5, 295–309. [Google Scholar]
  58. Hsu, L.F. Erratum to: An optimal policy for a single-vendor single-buyer integrated production-distribution model with both deteriorating and defective items. [Int. J. Prod. Econ. 147 (2014) 161–170]. Int. J. Prod. Econ. 2016, 178, 187–188. [Google Scholar] [CrossRef]
  59. Misra, R.B. Optimum production lot size model for a system with deteriorating inventory. Int. J. Prod. Res. 1975, 13, 495–505. [Google Scholar] [CrossRef]
  60. The United States Environmental Protection Agency (US. EPA). Emission Facts: Average Carbon Dioxide Emissions Resulting from Gasoline and Diesel Fuel. February 2005. Available online: https://nepis.epa.gov/ (accessed on 25 September 2018).
  61. McCarthy, J.E. EPA Standards for Greenhouse Gas Emissions from Power Plants: Many Questions, Some Answers. CRS Report for Congress, 7–5700. 2013. Available online: http://nationalaglawcenter.org/wp-content/uploads/assets/crs/R43127.pdf (accessed on 25 September 2018).
  62. Goyal, S.K. An integrated inventory model for a single supplier-single customer problem. Int. J. Prod. Res. 1977, 15, 107–111. [Google Scholar] [CrossRef]
Figure 1. Illustration of an integrated inventory model for imperfect quality items where defective products are sold at a discounted price.
Figure 1. Illustration of an integrated inventory model for imperfect quality items where defective products are sold at a discounted price.
Mathematics 07 00234 g001
Figure 2. Illustration of an integrated inventory model for imperfect quality items with vendor’s inspection.
Figure 2. Illustration of an integrated inventory model for imperfect quality items with vendor’s inspection.
Mathematics 07 00234 g002
Figure 3. Linear function of average fuel consumption vs. total weight for regional traffic (Data source: Volvo Corporation, 2018).
Figure 3. Linear function of average fuel consumption vs. total weight for regional traffic (Data source: Volvo Corporation, 2018).
Mathematics 07 00234 g003
Figure 4. Manufacturer’s and retailer’s inventory model for constantly deteriorating items with retailer’s inspection for n = 5.
Figure 4. Manufacturer’s and retailer’s inventory model for constantly deteriorating items with retailer’s inspection for n = 5.
Mathematics 07 00234 g004
Figure 5. Retailer’s inventory model with imperfect quality per delivery cycle.
Figure 5. Retailer’s inventory model with imperfect quality per delivery cycle.
Mathematics 07 00234 g005
Figure 6. Ip1(t1) and Ip2(t2) vs. time.
Figure 6. Ip1(t1) and Ip2(t2) vs. time.
Mathematics 07 00234 g006
Figure 7. Manufacturer’s and retailer’s inventory model for constant deteriorating items with manufacturer’s inspection for n = 5.
Figure 7. Manufacturer’s and retailer’s inventory model for constant deteriorating items with manufacturer’s inspection for n = 5.
Mathematics 07 00234 g007
Figure 8. Graphical representation of expected total cost (ETC) for a fixed n in Model 1.
Figure 8. Graphical representation of expected total cost (ETC) for a fixed n in Model 1.
Mathematics 07 00234 g008
Table 1. Gap analysis with existing literature.
Table 1. Gap analysis with existing literature.
AuthorsImperfect QualityDeteriorating ItemVariable Transportation CostCarbon Emission
Vendor’s InspectionBuyer’s Inspection
Huang (2002)
Goyal et al. (2003)
Wee et al. (2006)
Wahab et al. (2011)
Benjaafar et al. (2013)
Lee and Kim (2014)
Bazan et al. (2014)
Bozorgi et al. (2014)
Jauhari et al. (2014)
Bozorgi (2016)
Ghosh et al. (2016)
Sarkar et al. (2016b)
Yu and Hsu (2017)
Sarkar et al. (2017)
Toptal and Çetinkaya (2017)
Bouchery et al. (2017)
Dwicahyani et al. (2017)
Wangsa (2017)
Li et al. (2017)
Anvar et al. (2018)
Hariga et al. (2018)
Ji et al. (2018)
Wang and Ye (2018)
Gosh et al. (2018)
Ma et al. (2018)
Darom et al. (2018)
Jauhari (2018)
Gautam and Khanna (2018)
Wangsa and Wee (2018)
Tiwari et al. (2018)
Kundu and Chakrabarti (2019)
This paper
Table 2. List of notations.
Table 2. List of notations.
SymbolDefinition
Ddemand rate (unit/year);
Pproduction rate (unit/year);
Rproduction quantity; R = PT1;
𝜃deterioration rate; (0 ≤ 𝜃 < 1);
uthe probability of defective products per delivery lot size;
xquality screening rate (unit/year);
icfixed quality inspection cost ($/cycle);
ucunit inspection cost ($/unit);
cretailer’s ordering cost ($/order);
hdretailer’s holding cost ($/unit/year);
ddretailer’s deteriorating cost ($/unit);
smanufacturer’s setup cost ($/order);
hpmanufacturer’s holding cost ($/unit/year);
dpmanufacturer’s deteriorating cost ($/unit);
tfmanufacturer’s fixed transportation cost per delivery ($/delivery);
tvfuel price for manufacturer’s variable transportation cost ($/liter);
ddistance traveled from vendor to buyer (km);
wproduct weight (ton/unit);
c1average vehicle fuel consumption when empty (liter/km);
c2average additional fuel consumption per ton of load (liter/km/ton);
Txcarbon emission tax ($/tonCO2);
Feaverage emissions from fuel combustion (tonCO2/liter);
Eeaverage emissions from electricity generation (tonCO2/kWh);
e1transportation emission cost ($/km); e1 = c1FeTx;
e2average additional transportation emission cost per unit product ($/unit/km); e2 = c2wFeTx;
ecaverage warehouse energy consumption per unit product (kWh/unit/year);
wewarehouse emissions cost per unit product ($/unit/year); we = ecEeTx;
Tcycle length;
T1production period for the manufacturer in each cycle;
T2nonproduction period for the manufacturer in each cycle;
Tiinspection time per delivery for the retailer;
Tbinventory cycle length per delivery for the retailer; Tb = T/n;
Ip(t)manufacturer’s inventory level at time t;
Ipd(t)manufacturer’s inventory for defective products at time t;
Id(t)retailer’s inventory level at time t;
ETCdretailer’s expected total cost per year ($/year);
ETCpmanufacturer’s expected total cost per year ($/year);
ETCjoint expected total cost per year ($/year);
ETEdretailer’s expected total carbon emissions per year (tonCO2/year);
ETEpmanufacturer’s expected total carbon emissions per year (tonCO2/year);
ETEjoint expected total carbon emissions per year (tonCO2/year);
Decision variables
Qdelivery lot size (unit);
nnumber of deliveries per order (positive integer).
*indicates optimal solution
Table 3. Expected total cost for different n in Model 1.
Table 3. Expected total cost for different n in Model 1.
nT2(10−5)T1 (10−5)T (10−5)ETCdETCpETCETE
14960170366632,348,9911,017,4003,366,39118.460
25641193775781,463,3231,568,4933,031,81621.333
35961204880091,122,0121,799,1932,921,20523.501
4616121178278941,5091,930,5092,872,01825.422
5630621678473830,6732,017,7482,848,42127.216
6642222068628756,3702,081,5262,837,89628.935
7 *651822408758703,6112,131,3112,834,92230.598
8660322698872664,6202,172,0672,836,68732.218
9668022958976634,9522,206,6512,841,60333.805
10675123209071611,8892,236,8232,848,71235.360
Table 4. Expected total cost for different n in Model 2.
Table 4. Expected total cost for different n in Model 2.
nT2(10−5)T1(10−5)T (10−5)ETCdETCpETCETE
14977170966862,041,0051,312,7953,353,80018.22
25653194175951,167,5061,844,1933,011,69921.07
3596320488011827,1782,067,7592,894,93723.21
4615121138264644,6842,195,2732,839,95725.12
5628221588441530,6522,280,0882,810,74126.91
6638421938577425,5682,342,1422,794,71128.63
7646722228689395,7142,390,6072,786,32230.29
8653822468784352,2402,430,2972,782,74731.92
9 *660122688869318,4112,463,9852,782,39633.52
10665822888945290,9232,493,3782,784,30135.11
Table 5. Comparison between Model 1 and Model 2.
Table 5. Comparison between Model 1 and Model 2.
Decision Variables and Cost ItemsModel 1Model 2Saving (%)
Number of deliveries per cycle (n*)79
Cycle time (T)0.087580.08869
Delivery lot size (Q); units6387.74929.6
Ordering cost ($)22,835.722,550.71.25
Inspection cost ($)295,230.70100
Inventory holding cost ($)190,027.5147,863.722.2
Deteriorating cost ($)195,346.3147,863.724.3
Emission cost ($)171.0133.122.2
Total retailer’s cost per year ($)703,611.2318,909.154.7
Total retailer’s cost per year after compensation ($) 690,574.41.85
Setup cost ($)1,141,784.61,127,534.21.25
Inspection cost ($)0261,366.3−100
Transportation cost ($)85,289.5107,690.0−26.3
Inventory holding cost ($)539,976.9567,658.3−5.13
Deteriorating cost ($)362,122.2397,345.0−9.73
Emission cost ($)2,138.02390.9−11.8
Total manufacturer’s cost per year ($)2,131,311.22,463,984.8−15.6
Total manufacturer’s cost per year after compensation ($) 2,091,821.61.85
Expected total cost ($)2,834,922.42,782,396.01.85
Expected total emission (tonCO2/year)30.598 33.523−9.56
Table 6. Comparison between Model 1 and the adjusted Model 2.
Table 6. Comparison between Model 1 and the adjusted Model 2.
Decision Variables and Cost ItemsModel 1Model 2 adjSaving (%)
Number of deliveries per cycle (n*)77
Cycle time (T)0.087580.08704
Delivery lot size (Q); units6,387.76221.2
Ordering cost ($)22,835.722,977.4−0.62
Inspection cost ($)295,230.70100
Inventory holding cost ($)190,027.5186,596.11.80
Deteriorating cost ($)195,346.3186,596.14.48
Emission cost ($)171.0167.91.81
Total retailer’s cost per year ($)703,611.2396,337.543.7
Total retailer’s cost per year after compensation ($) 690,428.81.87
Setup cost ($)1,141,784.61,148,869.4−0.62
Inspection cost ($)0261,461.4−100
Transportation cost ($)85,289.585,764.5−0.55
Inventory holding cost ($)539,976.9529,446.71.95
Deteriorating cost ($)362,122.2357,812.11.20
Emission cost ($)2,138.02117.70.95
Total manufacturer’s cost per year ($)2,131,311.22,385,471.8−11.9
Total manufacturer’s cost per year after compensation ($) 2,091,380.51.87
Expected total cost ($)2,834,922.42,781,809.31.87
Expected total emission (tonCO2/year)30.59830.2921.00
Table 7. Sensitivity analysis of the two models.
Table 7. Sensitivity analysis of the two models.
ParameterValue ChangeModel 1Model 2
n *TQ aETC%CTCn *TQ bETC%CTC
P = 2,000,000+50%70.08165949.33,024,2496.6880.08195120.32,966,117.16.60
+25%70.08396117.12,948,6154.0180.08425263.22,892,821.23.97
070.08766387.72,834,922090.08874929.62,782,396.00
−25%80.09596120.62,643,862−6.7490.09575322.02,596,679.1−6.67
−50%80.11477321.02,253,495−20.5110.11535296.52,217,427.4−20.3
D = 500,000+50%80.08167811.63,190,85812.5100.08226168.23,136,392.112.7
+25%80.08416705.83,044,1907.3890.08405833.62,989,178.57.43
070.08766387.72,834,922090.08874929.62,782,396.00
−25%70.09585242.72,548,150−10.180.09624510.92,499,238.8−10.2
−50%60.11004683.32,151,051−24.180.11203501.72,108,977.7−24.2
c = 2000+50%70.08806415.92,846,3150.4090.08914951.62,793,646.20.40
+25%70.08786401.82,840,6250.2090.08894940.62,788,027.50.20
070.08766387.72,834,922090.08874929.62,782,396.00
−25%70.08746373.52,829,207−0.2090.08854918.62,776,752.0−0.20
−50%70.08726359.32,823,479−0.4090.08834907.52,771,095.4−0.41
s = 100,000+50%90.10746090.83,348,78418,1110.10814917.33,292,159.718.3
+25%80.09796250.13,104,5469.51100.09884945.13,049,823.59.61
070.08766387.72,834,922090.08874929.62,782,396.00
−25%60.07606464.72,529,737−10.880.07734834.42,480,006.3−10.9
−50%50.06256383.92,169,316−23.560.06315257.72,122,805.7−23.7
ic = 500+50%70.08836437.12,854,8270.7090.08884935.12,785,213.30.10
+25%70.08796412.42,844,8940.3590.08874932.42,783,805.10.05
070.08766387.72,834,922090.08874929.62,782,396.00
−25%70.08726362.92,824,912−0.3590.08864926.82,780,986.3−0.05
−50%70.08696338.02,814,883−0.7190.08864924.12,779,575.6−0.10
uc = 0.5+50%70.08766387.52,962,5574.5090.08874929.02,910,260.44.60
+25%70.08766387.62,898,7402.2590.08874929.32,846,328.12.30
070.08766387.72,834,922090.08874929.62,782,396.00
−25%70.08766387.82,771,105−2.2590.08874929.92,718,463.9−2.30
–50%70.08766387.92,707,289−4.5090.08874930.22,654,531.7−4.60
hd = 60+50%90.08724947.92,916,2532,87110.08804003.72,848,604.12.38
+25%80.08735571.52,878,4571.54100.08834014.82,817,636.91.27
070.08766387.72,834,922090.08874929.62,782,396.00
−25%60.08827502.82,782,702−1.8470.08856327.82,739,318.4−1.55
−50%50.08939119.22,716,341−4.1850.08898901.62,680,745.8−3.65
hp = 40+50%50.07797955.93,075,1298.4750.07767768.53,033,547.39.03
+25%60.08237002.82,962,6504.5170.08275908.52,915,155.64.77
070.08766387.72,834,922090.08874929.62,782,396.00
−25%80.09405999.52,691,834−5.05100.09514758.02,634,200.0−5.33
−50%90.10205786.32,531,836−10.7130.10514044.52,468,523.7−11.3
dd = 600+50%90.0872494412,918,3182,94110.08804003.72,848,604.12.38
+25%80.08735569.02,879,6191.58100.08834416.42,817,636.91.27
070.08766387.72,834,922090.08874929.62,782,396.00
−25%60.08827507.42,781,144−1.9070.08856327.82,739,318.4−1.55
−50%30.087214852.72,703,928−4.6250.08898901.62,680,745.8−3.65
dp = 400+50%70.08205981.93,010,0616,1860.08066721.32,957,798.86.30
+25%70.08476174.82,923,9243.1470.08395999.62,875,189.03.33
070.08766387.72,834,922090.08874929.62,782,396.00
−25%80.09225884.42,738,923−3.39100.09334668.12,679,457.1−3.70
−50%90.09765538.22,632,492−7.14120.09974155.52,566,509.0−7.76
𝜃 = 0.1+50%70.07945794.43,100,0019.3590.08044470.63,041,924.79.33
+25%70.08326069.22,970,7374.7990.08434683.22,915,364.14.78
070.08766387.72,834,922090.08874929.62,782,396.00
−25%70.09276762.92,691,447−5.0690.09395219.82,641,938.0−5.05
−50%70.09897213.92,538,850−10.480.09925516.82,492,404.9−10.4
tf = 1000+50%60.08747440.52,872,4361.3270.08836308.32,826,288.01.58
+25%70.08836437.12,854,8270.7080.08865542.62,805,413.70.83
070.08766387.72,834,922090.08874929.62,782,396.00
−25%80.08795612.02,814,046−0.74100.08854425.02,756,197.9−0.94
−50%90.08804991.12,790,972−1.55120.08843685.72,725,660.0−2.04
tv = 0.75+50%70.08766391.82,837,6050.0990.08884933.72,785,501.10.11
+25%70.08766389.72,836,2640.0590.08874931.62,783,948.80.06
070.08766387.72,834,922090.08874929.62,782,396.00
−25%70.08766385.72,833,581−0.0590.08864927.62,780,842.7−0.06
−50%70.08756383.62,832,239−0.0990.08864925.52,779,289.1−0.11
d = 100+50%70.08776392.82,838,3090.1290.08884934.72,786,313.00.14
+25%70.08766390.32,836,6160.0690.08874932.22,784,355.00.07
070.08766387.72,834,922090.08874929.62,782,396.00
−25%70.08756385.12,833,229−0.0690.08864927.02,780,436.7−0.07
−50%70.08756382.62,831,535−0.1290.08864924.52,778,476.4−0.14
E[u] = 0.02+50%70.08746437.82,843,9080.3290.08874931.32,784,137.70.063
+25%70.08756412.62,839,4140.1690.08874930.52,783,248.60.031
070.08766387.72,834,922090.08874929.62,782,396.00
−25%70.08776363.12,830,435−0.1690.08874928.72,781,579.0−0.029
−50%70.08786338.72,825,951−0.3290.08874927.72,780,796.5−0.058
w = 0.01+50%70.08766387.72,835,9560.0490.08874929.62,783,408.90.036
+25%70.08766387.72,835.4390.0290.08874929.62,782,902.60.018
070.08766387.72,834,922090.08874929.62,782,396.00
−25%70.08766387.72,834,405−0.0290.08874929.62,781,889.5−0.018
−50%70.08766.387.72,833,889−0.0490.08874929.62,781,383.0−0.036
c1, c2= 0.27, 0.0057+50%70.08766388.72,835,6270.0290.08874930.52,783,244.80.031
+25%70.08766388.22,835,2750.0190.08874930.02,782,838.60.016
070.08766387.72,834,922090.08874929.62,782,396.00
−25%70.08766387.22,834,570−0.0190.08874928.92,782,026.3−0.013
−50%70.08766386.72,834,218−0.0290.08874928.42,781,620.2−0.028
ec = 1.44+50%70.08766386.62,835,3720.01690.08874928.72,782,845.70.016
+25%70.08766387.12,835,1480.00890.08874929.22,782,620.80.008
070.08766387.72,834,922090.08874929.62,782,396.00
−25%70.08766388.32,834,698−0.00890.08874930.02,782,171.2−0.008
−50%70.08766388.82,834,472−0.01690.08874930.52,781,946.2−0.016
Tx = 75+50%70.08766387.72,836,0770.0490.08874929.82,783,658.00.045
+25%70.08766387.72,835,5000.0290.08874929.72,783,027.10.023
070.08766387.72,834,922090.08874929.62,782,396.00
−25%70.08766387.72,834,345−0.0290.08874929.52,781,764.9−0.023
−50%70.08766387.62,833,768−0.0490.08874929.42,781,134.2−0.045

Share and Cite

MDPI and ACS Style

Daryanto, Y.; Wee, H.M.; Widyadana, G.A. Low Carbon Supply Chain Coordination for Imperfect Quality Deteriorating Items. Mathematics 2019, 7, 234. https://doi.org/10.3390/math7030234

AMA Style

Daryanto Y, Wee HM, Widyadana GA. Low Carbon Supply Chain Coordination for Imperfect Quality Deteriorating Items. Mathematics. 2019; 7(3):234. https://doi.org/10.3390/math7030234

Chicago/Turabian Style

Daryanto, Yosef, Hui Ming Wee, and Gede Agus Widyadana. 2019. "Low Carbon Supply Chain Coordination for Imperfect Quality Deteriorating Items" Mathematics 7, no. 3: 234. https://doi.org/10.3390/math7030234

APA Style

Daryanto, Y., Wee, H. M., & Widyadana, G. A. (2019). Low Carbon Supply Chain Coordination for Imperfect Quality Deteriorating Items. Mathematics, 7(3), 234. https://doi.org/10.3390/math7030234

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop