Next Article in Journal
Changes in Energy Consumption in Agriculture in the EU Countries
Previous Article in Journal
EMC Component Modeling and System-Level Simulations of Power Converters: AC Motor Drives
Previous Article in Special Issue
Effect of Renewable Energy to Reduce Carbon Emissions under a Flexible Production System: A Step Toward Sustainability
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Impact of Energy and Carbon Emission of a Supply Chain Management with Two-Level Trade-Credit Policy

1
Department of Mathematics, Inderprastha Engg. College, Ghaziabad 201010, India
2
Department of Mathematics, Chaudhary Charan Singh University, Meerut 250001, India
3
Department of Mathematics, Vardhaman College, MJP Rohilkhand University, Uttar Pradesh 246701, India
4
Department of Industrial Engineering, Yonsei University, 50 Yonsei-ro, Sinchon-dong, Seodaemun-gu, Seoul 03722, Korea
5
Information Technology Research Center, Chung-Ang University, Seoul 06974, Korea
*
Author to whom correspondence should be addressed.
Energies 2021, 14(6), 1569; https://doi.org/10.3390/en14061569
Submission received: 20 December 2020 / Revised: 17 February 2021 / Accepted: 25 February 2021 / Published: 12 March 2021

Abstract

:
Supply chain management aims to integrate environmental thinking with efficient energy consumption into supply chain management. It includes a flexible manufacturing process, more product delivery to customers, optimum energy consumption, and reduced waste. The manufacturing process can be made more flexible through volume agility. In this scenario, production cannot be constant, and with the concept of volume agility, production is taken as a decision variable under the effect of optimum energy consumption. Considering a two-echelon supply chain, we consider a producer and supplier with two-level-trade-credit policies (TLTCP) with the optimum consumption. To reduce the integrated total inventory cost, we believe that demand is a function of the credit period and selling price. The cost function is analyzed, either with the credit period dependent demand rate or with the selling price dependent demand rate through the numerical examples under energy costs. Energy and carbon emission costs are introduced in setup/ordering cost, holding cost, and item cost for producer and supplier. The effect of inflation on the total cost cannot be ignored; this model is being developed for deteriorating items with the simultaneous impact of volume agility, energy, carbon emission cost, and two-level-trade-credit policies with inflation. This supply chain model was solved analytically and obtained the optimum decision variables in a quasi-closed form solution. An illustrative theorem is being utilized to analyze the optimum result for all the decision parameters. The convexity of the objective function is being obtained analytically as well as graphically. Finally, numerical examples and sensitivity analysis are employed to illustrate the present study and with managerial insights.

1. Introduction

The efficient use of energy can make a supply chain more efficient. Supply chain management (SCM) is characterized by complementing the traditional economic (cost) focus on environmental (emissions) considerations only without the concept of carbon emission. This research gap for the effect of energy and carbon emission can be conducted in this study. It is challenging to reduce carbon emission and energy consumption in traditional production, but energy consumption can be reduced by controlling the production rate. Due to a lack of coordination in the supply chain, the impact of deterioration from one stage to other increases. Effects of carbon emission in a global sustainable supply chain play an important role. Carbon emission and energy cost may happen in many different production stages while holding the produced stock, placing the order, and item cost.
Consequently, due to the deterioration of products, waste is generated in every supply chain stage. It poses a severe risk to nature and human health. It introduces global warming to the environment. Global warming is increasing due to the greenhouse gas emissions (GHGs), such as carbon dioxide (CO2), methane, and nitrous oxide. Cárdenas-Barrón and Sana [1] introduced a model to show the payment policies and reduce costs without considering waste items. Taleizadeh et al. [2] introduced a reworked strategy within a production model. They proved that multiple shipments could reduce the waste quickly compared to a single shipment with a large lot. Taleizadeh et al. [3] developed an improved strategy to repair defective products to reduce waste, but they did not consider controllable production rates. They just thought of a fixed production rate. Chen et al. [4] formulated a model on management of warehouses with low carbon under the effect of the cap-and-trade policy. This model proved that by utilizing cap and trade policy within production management, the carbon emission could be controlled, which is significant support for an environmental issue. Sarkar et al. [5] formulated a model based on the effect of multi-delay-in-payments and variable carbon emission cost. The single-setup-multiple-delivery policy is also discussed for the reduction of waste production. Khanna et al. [6] introduced a two-level-trade-credit policy (TLTCP) in an integrated model. However, they did not consider the reduction of those defective items by any strategy. Aljazzar et al. [7] formulated a delay-in-payment policy model to decrease carbon emission from the integrated model. Many researchers have tried to manage and reduce the cost of transportation and the cost of carbon emission. Tiwari et al. [8] established a multi-item sustainable green production process with the effect of partial backordering and payment delay policy. Sarkar [9] developed a model for defective goods with a multi-stage production system. Tayyab and Sarkar [10] investigated a model on optimal quantity under a multi-stage production with a random defective rate without any waste reduction policy or cost reduction policy. Kim and Sarkar [11] analyzed a model for upgrading the quality of items with a multi-stage cleaner production process. In this model, lead time-dependent ordering cost is assumed. Ahmed and Sarkar [12] formulated a model on sustainable supply chain management (SSCM) for a second-generation biofuel on carbon emission in an integrated model. Sarkar and Chung [13] developed a flexible work-in-process production system in SCM under quality improvement.
Agile manufacturing has also been applied to the economic production quantity (EPQ) model. The controllable production rate can be easily used through agile manufacturing to increase or decrease the production rate based on the management’s necessity or based on the system required. It indicates that the production rate must vary with the situation, whatever is necessary at the right time with the right judgment. This production rate is assumed as a decision variable. Sana and Chaudhuri [14] formulated a model using stock-dependent demand with volume flexibility for deteriorating items. According to the authors, the unit production cost is a function of the finite production rate, treated as a decision variable. Sana and Chaudhuri [15] extended this model by using a time-dependent demand rate. Physical delay of the stocked item over time is taken into account, and inventory shortages are allowed. Singh et al. [16] formulated a supply chain model with variable holding costs for a flexible manufacturing system. Singh et al. [17] inspected an ordering model with volume agility, variable demand, and inflation. Dem and Singh [18] studied an ordering system with multiple items and greening under the flexible production rate assumption. Sangal and Gupta [19] formulated a two-warehouse vendor-supplier framework with volume agility. Singh and Gupta [20] discussed the theory based on the supply chain model with error in quality inspection with selling price dependent demand rate under the effect of volume agility. Gautam et al. [21] examined a model on sustainable production policies under the impact of volume agility, preservation technology, and price-reliant demand.
An EPQ model depends on the supposition that as the supplier receives the producer’s goods, he must pay the goods’ cost. Nevertheless, in today’s scenario, it is very general to observe that the producer will allocate a specific time interval for paying the total price of items that the producer owes to the supplier for the goods, known as the trade-credit period (TCP). Generally, interest is not charged if it is reimbursed during the time given by the producer. There are two benefits of this credit-period, which the producer offers.
  • This price reduction policy for green products attracts new customers.
  • It should cause a decline in sales since it takes time to profit from this delay period more often, and some customers will pay more rapidly.
In supply chain management, TLTCP plays an important role. In this credit policy, the producer offers a credit period to the supplier, and the supplier also provides some credit period to his customers. Ouyang et al. [22] analyzed a strategy with TCP and non-instantaneous delivery. Jaggi et al. [23] proposed a model with TLTCP to assume that the demand rate is a function of that period given for paying the money. Sarkar [24] proposed a model that is based on TCP. Here deterioration rate depends on time, but he did not consider any idea regarding agile manufacturing or its effect. The author just developed a delay-period concept on the basic inventory model and did not find any waste reduction in the inventory system. Khanra et al. [25] extended this idea with the effect of shortages and time-varying demand rates. Here, the impact of TCP is also discussed. Chung [26] developed a production model based on TLTCP instead of the basic inventory model with limited storage capacity in an integrated system. Sarkar et al. [27] finally established a model with a price-discount offer, finite replenishment rate, and credit-period policy. Sarkar et al. [28] formulated a model for defective units. In this model, variable lead time and TCP are also discussed without reducing those faulty products. Sarkar et al. [29] analyzed the TCP policy’s effect on outcomes with a fixed lifetime. Here deterioration rate is taken as a variable. Yadav et al. [30] explained a retailer’s optimal policy under inflation in a fuzzy environment with TCP. The TLTCP is considered, and the demand rate depends on the selling price and credit period. An integrated total cost function is evaluated in three different cases and analyzes varying cost function parameters’ effects.
The concept of energy and carbon emission cost with agile manufacturing for product management is the significant study gap in inventory models. The planned model fulfills this study gap within an integrated model. Thus, TLTCP is utilized for this purpose. The energy and carbon emission costs can be reduced, production can be decided as per the situation, and the total cost can be globally minimized.

Research Gap

At the producer’s end, different production rates are considered to formulate the model more effective in the literature. Agile manufacturing received little attention in the literature at the producer’s end in the SCM. If the producer will decide the production rate as per the market situation, this affects the whole supply chain to reduce the holding cost and minimize the total cost. Here, to increase the demand, the producer offers some credit period to the supplier, and the supplier also provides some credit period to the customer. In all the above papers, the effect of TLTCP with an integrated model is unobserved. Therefore, a new agile manufacturing model has been developed to determine the supplier’s optimal TCP when producer and supplier propose a time-period to motivate their customers. In this research, the demand rate is a function of the credit period offered by the supplier and selling price, which is not considered with agile manufacturing in the above literature review. This study highlights agile manufacturing and TLTCP in the supply chain.
Further, it develops an assessment of the energy and carbon emission in setup/ordering cost, holding cost, and item cost. In this model, there is an analysis of total integrated cost with different cases of demand rate by numerical examples. Here production rate, credit period offered by the supplier, and complete cycle length are taken as decision variables. From the above literature review, the research gap is shown in contribution to Table 1.
The rest of this paper is ordered as follows: Section 2 shows the notations and assumptions. Section 3 shows the material and methods for producer and supplier. Illustrative examples, which explain the numerical verifications, sensitivity analysis, and observations, are provided in Section 4. Discussion and managerial insights represent the outcomes of this study in Section 5. Conclusion and suggestions for further research shown in Section 6 and Appendix A.

2. Problem Description, Notations, and Assumptions

2.1. Problem Description

A supply chain coordination is studied to minimize the energy and carbon emission cost with agile manufacturing. In agile manufacturing, per unit production cost is different from the raw material, and production rate is assumed as a decision variable. Here the demand rate varies with the credit period and selling price. The carbon emission and energy cost are found during setup/ordering cost, item cost, and holding cost. This research aims to control the energy and carbon emission cost, minimize the total cost, and maintain an improved strategy towards the integrated model. In this study, TLTCP plays a vital role in motivating the customers and increasing the demand.

2.2. Notation

2.2.1. Decision Variables

  • T: Producer’s cycle length;
  • N: Credit period offered by supplier;
  • P: Production rate for the producer.

2.2.2. Producer’s Parameters

  • C1p: Producer’s stock holding cost ($/unit/week);
  • Cp: Producer’s item cost ($/unit/week);
  • C2p: Producer’s setup cost ($/set up);
  • C’1p: Cost which include energy and carbon emission due to holding items ($/unit/week);
  • C’p: Cost which include energy and carbon emission due to deterioration ($/unit/week);
  • C’2p: Cost which include energy and carbon emission due to setup ($/setup);
  • φ(P): Production cost.

2.2.3. Supplier’s Parameters

  • M: Credit period offered by producer;
  • ie: Interest rate earned by supplier;
  • ip: Interest rate payable by supplier;
  • C1s: Holding cost of supplier ($/unit/week);
  • Cs: Item cost of supplier ($/unit);
  • C2s: Ordering cost ($/order);
  • C1s: Cost which include energy and carbon emission due to holding items ($/unit/week);
  • C’2s: Cost which include energy and carbon emission due to ordering items ($/order);
  • C’s: Cost which include carbon emission and energy due to deterioration ($/unit/week);
  • c: Purchasing cost ($/unit);
  • r: Inflation rate;
  • n: Number of cycles;
  • T2: Per delivery cycle length; T2= T/n.

2.3. Assumptions

  • As for the environment requirements, this paper assumed the carbon emission and energy cost in holding cost, item cost, and setup/ordering costs for producer and supplier (Ahmed and Sarkar [12]).
  • Due to agile manufacturing, a controllable production rate controls the excess holding or shortage any time during the whole cycle-length. There is no shortage, since the production rate is considered a decision variable (Sana and Chaudhuri [2]).
  • Here, the demand rate is variable, and demand varies with the supplier’s credit period and selling price. Therefore, D = W − (Wf)(1z)N + α pτ, where W and f are the maximum and initial demand of products (unit/week), respectively, z is the demand saturation rate (unit/week). α and τ are the scaling parameter and price elasticity parameter of the selling price p, respectively (Jaggi et al. [13]).
  • The deterioration rate (0 ≤ θ ≤ 1) exists and is very small.
  • In this paper, TLTCP is applied. In this policy, the producer gives a fixed time M to the supplier, and the supplier also proposes a period N to their clients for the settlement of the account (Jaggi et al. [13]).
  • A complete model is studied in the inflationary environment (Yadav et al. [20]).
  • Due to agile manufacturing, production cannot be constant. That is why production cost is given below (Sana and Chaudhuri [3]):
ϕ ( P ) = ( μ + ( k / P ) + s P + β ( P P c )   . g ( P P c ) ) ,   where   g   ( P P c ) = { 1 0 P > P c P P c .
Here,
  • µ is the material cost.
  • The second term is ( k / P ) is associated with labor and energy cost. As the production rate (P) increases, this cost decreases.
  • The term (sP) represents a tool or die costs.
  • The last term represents a critical value of the production rate (Pc) for the machine. The formed stock is to be defective for a very high production rate (P > Pc).

3. Material and Methods

In this section, we proposed a two-echelon integrated model for a producer and a supplier. Here the producer produces the goods with a production rate (P). These goods are transferred to the supplier in the number of shipments. The demand rate depends on the credit-period and selling price. The credit period is given to the supplier by the producer to return the number of goods. The retailer also allows a credit period to its customers.

3.1. Producer’s Model

The producer’s production model is formulated under agile manufacturing. The production process starts at t = 0 time and goes up to t = T1 time, where the stock level goes to its highest level. At this time, the inventory level goes up because of production and declines because of deterioration and demand. At time T1, production stops, and inventory level decreases because of demand and deterioration rate. The level of the stock goes down to zero at t = T. For the duration of [0, T], the stock system’s behavior is exhibited in Figure 1.
It is observed that the rate of inventory level increases due to the production rate and decreases due to demand and deterioration rate. The inventory position can be represented as follows:
d I p 1 ( t ) d t = θ I p 1 ( t ) D + P ,   0 t T 1 ,
At the time duration, T 1 t T , we observe that the inventory level rate decreases due to the demand and deterioration rate. The inventory position can be represented as follows:
d I p 2 ( t ) d t = θ I p 2 ( t ) D ,   T 1 t T ,
On solving the above equations by using boundary conditions I p 1 ( 0 ) = 0 and I p 2 ( T ) = 0 , we have, inventory level during the production time at any time t is as follows:
I p 1 ( t ) = ( P D θ ) ( 1 e θ t ) ,
Inventory level during the non-production time at any time t is as follows:
I p 2 ( t ) = ( D θ ) ( e θ ( T t ) 1 ) ,
With the initial condition Ip1(T1) = Im, the maximum inventory after production is
I m = ( P D θ ) ( 1 e θ T 1 ) ,
Using the equation of continuity of inventory I p 1 ( T 1 ) = I p 2 ( T 1 ) , we have
( P D θ ) ( 1 e θ T 1 ) = ( D θ ) ( e θ ( T T 1 ) 1 ) .
This gives the products production time as
T 1 = 1 θ log ( D e θ T + P D P ) .
The producer bears the material setup cost, holding cost, item cost, and production cost. The setup cost, holding cost, and item cost include energy and carbon emission cost. The following components give the total cost function:

3.1.1. Holding Cost

For maintaining the EURO-14, the industry must take care of the environmental cost and the extra amount of investment. Therefore, the holding cost consists of two components. One is associated due to the holding of the products as C1p, and the other part is due to the energy consumption and carbon emission as C1p. Finally, the total holding cost under the effect of the environment is
H C p = ( C 1 p + C 1 p ) T [ 0 T 1 I p 1 ( t ) e r t d t + T 1 T I p 2 ( t ) e r t d t ] , = ( C 1 p + C 1 p ) T [ ( P D θ ) { θ ( θ + r ) e r T 1 + r e ( θ + r ) T 1 r ( θ + r ) } + D θ ( e r T e r T 1 r + e θ T ( θ + r ) T 1 e r T θ + r ) ] .

3.1.2. Setup Cost

The producer’s setup cost consists of two components one is associated due to the setup of the products as C2p, and the other component is the energy consumption and carbon emission as C2p. Finally, the total setup cost under the effect of the environment is
S C p = ( C 2 p + C 2 p ) T .

3.1.3. Item Cost

The producer’s item cost consists of two components. One is associated with the deterioration and items sold as Cp and other components due to the energy consumption and carbon emission as Cp. Finally, the total item cost under the effect of the environment is
I C p = ( C p + C p ) T 0 T 1 P e r t d t = ( C p + C p ) T P ( 1 e r T 1 ) .

3.1.4. Production Cost

The producer’s production cost includes labor cost, material cost, tool/die cost, and some excess labor and energy cost.
P C = ϕ ( P ) [ 0 T 1 P e r t d t ] = ( μ + ( k / P ) + s P + β ( P P c ) g ( P P c ) ) ( ( 1 e r T 1 ) P r ) .

3.1.5. Producer’s Total Cost

The producer’s total cost can be represented by:
T C p = P C + D C p + S C p + H C p , T C p = ( μ + k P + s P + g ( P P c ) β ( P P c ) ) ( ( 1 e r T 1 ) P T r ) + ( C 2 p + C 2 p ) T + ( C p + C p ) T P ( 1 e r T 1 ) + ( C 1 p + C 1 p ) T [ ( P D θ ) { θ ( θ + r ) e r T 1 + r e ( θ + r ) T 1 r ( θ + r ) } + D θ ( e r T e r T 1 r + e θ T ( θ + r ) T 1 e r T θ + r ) ] .

3.2. Supplier’s Model for Selling Strategy of Products

The supplier’s inventory storage and demand pattern can be represented as in Figure 2. At the starting of the period, Is(t) is stock received from the producer. It is supposed that the supplier has a credit-period which the producer offers, and the supplier also gives a time interval to his clients. The level of stock Is(t) continuously goes down to meet the demand due to the deterioration and demand, and inventory level goes down to zero at t = T2. Here, the supplier supplies the inventory to their customers in n different orders.
During the time interval, 0 t T 2 , the inventory level decreases due to the demand and deterioration rate. Hence, the variation in the level of inventory versus time t is given below:
d I s ( t ) d t = θ I s ( t ) D 0 t T 2 ,
With the inventory status I s ( T 2 ) = 0, current inventory at any time t in the period 0 t T 2 is as follows:
I s ( t ) = ( D / θ ) ( e θ ( T 2 t ) 1 ) .
The supplier bears the ordering cost, holding cost, item cost, interest payable cost, and interest earned. The ordering cost, holding cost, and item cost include energy and carbon emission cost. The following components give the total cost function.

3.2.1. Holding Cost for Items

Here holding cost consists of two components. One is associated with the holding of the products as C1s, and the other component is the energy consumption and carbon emission as C1s. Finally, the total holding cost under the effect of the environment is
H C s = n ( C 1 s + C 1 s ) T 0 T 2 e r t I s ( t ) d t = n ( C 1 s + C 1 s ) D T θ ( e θ T 2 e r T 2 θ + r + e r T 2 1 r ) .

3.2.2. Ordering Cost for Items

The supplier’s ordering cost consists of two components. One is associated with the ordering of the products as C2s, and the other component is the energy consumption and carbon emission as C2s. Finally, the total ordering cost under the effect of the environment is
S C s = n ( C 2 s + C 2 s ) T .

3.2.3. Item Cost

The supplier’s item cost consists of two components. One is associated due to the deterioration and item sold as Cs, and the other component is due to the energy consumption and carbon emission as Cs. Finally, the total item cost under the effect of the environment is
I C s = n ( C s + C s ) T I s ( 0 ) = n ( C s + C s ) T ( W     ( W f ) ( 1 z ) N +   τ   p ρ θ ) ( e θ ( T 2 ) 1 ) .
There are different possibilities of permissible delay payment period. Variation of revenue versus time is shown in Figure 3a–c.

3.2.4. Supplier’s Total Cost

The supplier’s total cost can be represented by:
T C s = H C s + I C s + S C s , = n ( C 1 s + C 1 s ) D T θ ( e θ T 2 e r T 2 θ + r + e r T 2 1 r ) + n ( C s + C s ) T ( D θ ) ( e θ ( T 2 ) 1 ) + n ( C 2 s + C 2 s ) T

3.2.5. Supplier’s Total Cost in Two Different Credit Policy

Case I-When NMT2 + N

Here, TCP offered by the supplier to his clients is significantly less than the TCP proposed by the producer to the supplier. After getting the revenue, the supplier earns interest on the average sales revenue in the time interval (M–N). The finances are arranged to make the payment to the producer at time M. It is shown in Figure 3a.
Interest earned for the products can be expressed as:
I E 1 = n i e p T N M D e r t d t = n i e p D ( e r N e r M ) r T ,
Interest payable for the products can be defined as:
I P 1 = n i p c T M N + T e t r I s ( t ) d t = n i p c D T θ [ e θ T 2 ( e ( r + θ ) M e ( r + θ ) ( T 2 + N ) ) θ + r + ( e r ( T 2 + N ) e r M ) r ] ,
Supplier’s total cost TCs1 for the inventory per unit time can be expressed as:
T C s 1 = T C s + I P 1 I E 1 , = n ( C 1 s + C 1 s ) D T θ ( e θ T 2 e r T 2 θ + r + e r T 2 1 r ) + n ( C s + C s ) T ( D θ ) ( e θ ( T 2 ) 1 ) + n ( C 2 s + C 2 s ) T + n i p c D T θ [ e θ T 2 ( e r M θ M e ( T 2 + N ) ( r + θ ) ) θ + r + ( e r T 2 r N e r M ) r ] n i e p D ( e r N e r M ) r . T .

Case II-NT2 + NM

Here, TCP offered by the producer to the supplier is more than the cycle length of the supplier and the TCP provided by the supplier to his customers, and the supplier earns interest on the received average sales revenue during (N, T2 + N) and on total sales revenue for (M–(T2 + N)), but there is no interest payable by the supplier. It is represented in Figure 3b.
Interest earned for the inventory can be express as:
I E 2 = n i e p T [ N T 2 + N D e r t 2 d t 2 + D T 2 T 2 + N M e r t 2 d t 2 ] , = n i e p T [ D e r N r ( 1 e r T 2 r ) + D T 2 ( e r ( T 2 + N ) r e r M r ) ] ,
Supplier’s total cost TCs2 of inventory product per cycle can be express as:
T C s 2 = T C s I E 2 , = n ( C 1 s + C 1 s ) D T θ ( e θ T 2 e r T 2 θ + r + e r T 2 1 r ) + n ( C s + C s ) T ( D θ ) ( e θ ( T 2 ) 1 ) + n ( C 2 s + C 2 s ) T n i e p T [ D e r N r ( 1 e r T 2 r ) + D T 2 ( e r ( T 2 + N ) r e r M r ) ] .

Case III-MNT2 + N

Here, TCP offered by the supplier to his customers is more than that period provided by the producer to the supplier, and the supplier has no interest, but pays interest on the full order of products for a time (N–M) and average product held during the cycle T2. It is represented in Figure 3c.
Interest payable for a product is expressed as:
I P 2 = n i p c T [ D T 2 ( N M ) + N T 2 + N D e r t 3 d t 3 ] = n i p c D T [ T 2 ( N M ) + e r N r ( 1 e r T 2 r ) ] ,
Supplier’s total cost TCs3 for product per unit time can be express as:
T C s 3 = H C s + I C s + S C s + I P 2 , = n ( C 1 s + C 1 s ) D T θ ( e θ T 2 e r T 2 θ + r + e r T 2 1 r ) + n ( C s + C s ) T ( D θ ) ( e θ ( T 2 ) 1 ) + n ( C 2 s + C 2 s ) T + n i p c D T [ T 2 ( N M ) + e r N r ( 1 e r T 2 r ) ] ,
Therefore, the total cost is
T C s c = { T C s 1 i f N     M     T 2 + N   T C s 2 i f N     T 2 + N     M     T C s 3 i f M     N     T 2 + N         .
where T C s 1 , T C s 2 , and T C s 3 are given by Equations (20), (22) and (24) for products.
Finally, the system’s integrated cost comprises the producer’s cost (TCp) and the supplier’s cost (TCsc). The relevant cost for the whole model is different in various cases.

3.3. Integrated Cost Function for Two-Echelon Supply Chain Management

Three cases of the joint total cost are formed because of different situations of the supply chain, and can be developed as:
T C S 1 ( P ,   N ,   T ) = T C p + T C s 1 , N M T 2 + N , T C S 2 ( P ,   N ,   T ) = T C p + T C s 2 , N   T 2 + N M , T C S 3 ( P ,   N ,   T ) = T C p + T C s 3 , M N T 2 + N .
This theory develops a production model for items with a specific credit period for the supplier. The inventory cost model is derived from the policy of optimal production and credit-period.

3.3.1. The Total Integrated Cost in Case I

T C S 1 ( P ,   N ,   T ) = b ( P ( 1 e r T 1 ) T r ) + ( C 1 p + C 1 p ) T [ ( P D θ ) { θ ( θ + r ) e r T 1 + r e ( θ + r ) T 1 r ( θ + r ) } + D θ ( e r T e r T 1 r + e θ T ( θ + r ) T 1 e r T θ + r ) ] + ( C 2 p + C 2 p ) T + ( C p + C p ) T P ( 1 e r T 1 ) + n ( C 1 s + C 1 s ) D T θ ( e θ T 2 e r T 2 θ + r + e r T 2 1 r ) + n ( C s + C s ) T ( D θ ) ( e θ ( T 2 ) 1 ) + n ( C 2 s + C 2 s ) T + n i p c D T θ [ e θ T 2 ( e r M θ M e ( T 2 + N ) ( r + θ ) ) θ + r + ( e r T 2 r N e r M ) r ] n i e p D ( e r N e r M ) r T .

3.3.2. The Total Integrated Cost in Case II

T C S 2 ( P ,   N ,   T ) = b ( P ( 1 e r T 1 ) T r ) + ( C 1 p + C 1 p ) T [ ( P D θ ) { θ ( θ + r ) e r T 1 + r e ( θ + r ) T 1 r ( θ + r ) } + D θ ( e r T e r T 1 r + e θ T ( θ + r ) T 1 e r T θ + r ) ] + ( C 2 p + C 2 p ) T + ( C p + C p ) T P ( 1 e r T 1 ) + n ( C 1 s + C 1 s ) D T θ ( e θ T 2 e r T 2 θ + r + e r T 2 1 r ) + n ( C s + C s ) T ( D θ ) ( e θ ( T 2 ) 1 ) + n ( C 2 s + C 2 s ) T n i e p D T r [ e r N ( 1 ( e r T 2 r ) ) + T 2 ( e r ( T 2 + N ) e r M ) ] .

3.3.3. The Integrated Total Cost in Case III

T C S 3 ( P ,   N ,   T ) = b ( P ( 1 e r T 1 ) T r ) + ( C 1 p + C 1 p ) T [ ( P D θ ) { θ ( θ + r ) e r T 1 + r e ( θ + r ) T 1 r ( θ + r ) } + D e r T 1 θ ( e r T e r T 1 r + e θ T ( θ + r ) T 1 e r T θ + r ) ] + ( C 2 p + C 2 p ) T + ( C p + C p ) T P ( 1 e r T 1 ) + n ( C 1 s + C 1 s ) D T θ ( e θ T 2 e r T 2 θ + r + e r T 2 1 r ) + n ( C s + C s ) T ( D θ ) ( e θ ( T 2 ) 1 ) + n ( C 2 s + C 2 s ) T + n i p c D T [ T 2 ( N M ) + e r N r ( 1 e r T 2 r ) ] .

3.4. Solution Methodology

Kuhn-Tucker’s method is the best approach to prove the convexity of this model. From the necessary conditions of optimization, one can obtain:
T C S 1 P = ( C p + C p ) T ( ( 1 m ) m r D ( e θ T 1 ) P θ ) + P ( 1 m ) r T ( k P 2 + s + β ) + ( C 1 p + C 1 p ) T { ( l e θ T m ) D 2 ( e θ T 1 ) P 2 θ 2 + l r + θ m ( r + θ ) r θ ( θ + r ) + ( P D ) D ( e θ T 1 ) ( l m ) θ P 2 θ } + ( 1 m ) b r T m b D ( e θ T 1 ) T P θ = 0 ,
T C S 1 N = D e r N n p i e T + n i p c D K 4 T θ + K n p i e D 1 r T ( C p + C p ) D l 2 r T θ ( C s + C s ) ( e T θ n 1 ) n D 1 T θ n i p c D 1 T θ ( K 2 r + e T θ / n K 1 r + θ ) + n ( C 1 s + C 1 s ) D 1 T θ ( e r ( T / n ) 1 r + K 3 r + θ ) D b l 2 T θ + ( C 1 p + C 1 p ) T [ ( e r T m r + e r T + l e T θ r + θ ) D 1 θ + ( r l + θ m ( r + θ ) ) D 1 r ( r + θ ) θ + D θ ( l 2 P θ + e T θ l 1 P θ ) + ( P D ) ( l 1 l 2 ) P θ 2 ] = 0 ,
T C S 1 T = ( C p + C p ) P ( 1 m ) T 2 ( C 2 p + C 2 p ) T 2 n ( C 2 s + C 2 s ) T 2 + ( C s + C s ) e θ T / n D T + K n p i e D r T 2 + ( C p + C p ) r D m e T θ T n ( C s + C s ) ( e T θ n 1 ) D T 2 θ n i p c D T 2 θ ( K 2 r + e T θ / n K 1 r + θ ) + n ( C 1 s + C 1 s ) D T 2 θ ( e r ( T / n ) 1 r + e T θ / n e r ( T / n ) r + θ ) + n i p c D T θ ( K 4 n + θ e θ T / n K 1 n ( r + θ ) ) + n ( C 1 s + C 1 s ) D T θ ( e r ( T / n ) n r e r ( T / n ) + θ e θ ( T / n ) n ( r + θ ) ) ( C 1 p + C 1 p ) T 2 { D θ ( e r T m r + e r T + l e T θ r + θ ) + ( P D ) ( l m ) r θ } + ( C 1 p + C 1 p ) T 2 { ( P D ) D r e T θ ( l m ) P r θ + D θ ( e r T + l e T θ D P + r e r T + m e T θ r + θ ) } = 0 .
The sufficient condition must be satisfied to prove the globality of the optimal solution. Now, to acquire the optimal solution, a theorem is formulated as follows:
Theorem 1.
If ξ 1 > 0 , ξ 1 ξ 2 ξ 4 2 > 0 , ( ξ 1 ξ 2 ξ 3 ξ 1 ξ 6 2 ξ 3 ξ 4 2 + 2 ξ 4 ξ 5 ξ 6 ξ 2 ξ 5 2 )   > 0 , then TCS1 (P, N, T) will be minimum at point (P*, N*, T*)
Proof. 
To find the sufficient condition of the global optimality, taking second-order derivatives of Equations (29)–(31) versus P, T, and N, one can find:
2 T C S 1 N 2 = D e r N n r p i e T + n i p c D ( r K 4 θ e ( r θ ) ( N + T / n ) + T θ / n ) T θ 2 D 1 e r N n p i e T 2 n i p c D 1 K 4 T θ + n K i e p D 1 T r 2 ( C p + C p ) l 2 r T θ 2 ( C s + C s ) ( e T θ n 1 ) n D 1 T θ 2 n i p c D 1 T θ ( K 2 r + e T θ / n K 1 r + θ ) + 2 n ( C 1 s + C 1 s ) D 1 T θ ( e r ( T / n ) 1 r + K 3 r + θ ) 2 l 2 b T θ + ( C 1 p + C 1 p ) T [ ( e r T m r + e r T + l e T θ r + θ ) 2 D 1 θ + 2 ( r l + θ m ( r + θ ) ) D 1 r ( r + θ ) θ + D θ ( 2 l 2 P θ + 2 e T θ l 1 P θ ) + ( P D ) ( l 1 l 2 ) P θ 2 ] =   ξ 1 ( s a y ) ,
2 T C S 1 P 2 = 2 ( 1 m ) ( s + β ) r T m D ( e θ T 1 ) P θ T ( 2 ( k P 2 + s + β ) + ( C p + C p ) ( e θ T 1 ) r 2 D P 2 ) + ( C 1 p + C 1 p ) T { D θ ( D 2 ( e θ T 1 ) 2 ( l ( r + θ ) e θ T m r ) P 4 θ 2 + 2 ( m l e θ T ) D ( e θ T 1 ) P 3 θ ) + ( P D ) D ( e θ T 1 ) θ ( D ( e θ T 1 ) ( l ( r + θ ) m r ) P 4 θ 2 2 ( l m ) P 3 θ ) + 2 ( l m ) D ( e θ T 1 ) P 2 θ 2 } m r b D 2 ( e θ T 1 ) 2 P 3 θ 2 T = ξ 2 ( s a y ) ,
2 T C S 1 T 2 = 2 ( C p + C p ) P ( 1 m ) T 3 + 2 ( C 2 p + C 2 p ) T 3 + 2 n ( C 2 s + C 2 s ) T 3 2 ( C s + C s ) e θ T / n D T 2 2 K n p i e D r T 3 ( C p + C p ) r D m e T θ T 2 + 2 n ( C s + C s ) ( e T θ n 1 ) D T 3 θ + 2 ( C s + C s ) n ( e θ T / n 1 ) D T 3 θ + θ ( C s + C s ) e θ T / n D n T + ( C p + C p ) r D m e T θ T ( θ r D e T θ T ) + 2 n i p c D T 3 θ ( e r ( N + T / n ) e r M r + e T θ / n K 1 r + θ ) + 2 n ( C 1 s + C 1 s ) D T 3 θ ( e r ( T / n ) 1 r + e T θ / n e r ( T / n ) r + θ ) 2 n i p c D T 2 θ ( K 4 n + θ e T θ / n K 1 n ( r + θ ) ) + n i p c D T θ ( r e r ( N + T / n ) + θ e ( r θ ) ( N + T / n ) + T θ / n n 2 + θ 2 e θ T / n K 1 n 2 ( r + θ ) r e ( r θ ) ( N + T / n ) + T θ / n n 2 ) 2 n ( C 1 s + C 1 s ) D T 2 θ ( e r ( T / n ) n r e r ( T / n ) + θ e θ ( T / n ) n ( r + θ ) ) + n ( C 1 s + C 1 s ) D T θ ( r e r ( T / n ) n 2 + r 2 e r ( T / n ) + θ 2 e θ ( T / n ) n 2 ( r + θ ) ) + 2 ( C 1 p + C 1 p ) T 3 { D θ ( e r T m r + e r T + l e T θ r + θ ) + ( P D ) ( l m ) r θ } 2 ( C 1 p + C 1 p ) T 2 { ( P D ) D e T θ ( l m ) P θ + D θ ( e r T + l e T θ D P + r e r T + l e T θ ( θ e T θ D ( r + θ ) ) / P r + θ ) } + ( C 1 p + C 1 p ) T { ( P D ) P θ ( e T θ m D ( θ r e T θ D P ) + l e T θ D P ( θ + ( r + θ ) e T θ D P ) ) + D θ ( r 2 e r T + m r e T θ D P ( θ r e T θ D P ) + 1 r + θ ( r 2 e r T l e 2 T θ D θ ( r + θ ) P + e T θ l ( θ ( r + θ ) e T θ D P ) 2 ) ) } + 2 P b ( 1 m ) r T 3 2 e T θ m D b T 2 + e T θ m D b T ( θ r e T θ D P ) = ξ 3 ( s a y ) ,
2 T C S 1 N T = K n p i e T 2 ( D e r N + K D 1 r ) + i p c T θ ( n K 4 D T r K 4 D + n D 1 T ( K 2 r + K 1 e T θ / n r + θ ) n D 1 ( K 4 + K 1 θ e T θ / n r + θ ) ) D 1 e T θ ( r m ( C p + C p ) + ( C s + C s ) ) T + n ( C s + C s ) D 1 ( e T θ n 1 ) T 2 θ + ( C p + C p ) r l 2 T θ ( 1 T + D r 2 e T θ P ) + n ( C 1 s + C 1 s ) D 1 T θ ( e r ( T / n ) n + r e r ( T / n ) + θ e θ ( T / n ) n ( r + θ ) ( e r ( T / n ) 1 ) r K 3 r + θ ) m b e T θ D 1 T + b l 2 T θ ( 1 T + r e T θ P ) ( C 1 p + C 1 p ) T 2 { D 1 θ ( e r T m r + e r T + l e T θ r + θ ) + D 1 ( l r + θ m r m θ ) r θ ( r + θ ) D ( l 2 + l 1 e T θ ) P θ 2 + ( P D ) ( l 1 l 2 ) θ } + ( C 1 p + C 1 p ) T 2 { D 1 e T θ ( l 1 + l 2 ) θ D 1 θ ( e r T + m D e T θ P ) ( P D ) D 1 e T θ ( l + m ) θ P + D θ ( m e T θ D 1 P + l e T θ D 1 P ) } = ξ 5 ( s a y ) ,
2 T C S 1 N P = r 2 D l 2 ( b + ( C p + C p ) ) P θ 2 l 2 r ( k P 2 + s + β ) + ( C 1 p + C 1 p ) T { D ( l 1 e T θ l 2 ) P 2 θ 2 + D ( l 1 l 2 ) θ + D θ ( ( l 1 e T θ + l 2 ) P 2 θ + ( e T θ 1 ) D ( l 1 e T θ ( r + θ ) r l 2 ) P 2 θ 2 ) + ( P D ) P 2 θ 2 ( l 1 + l 2 + ( e T θ 1 ) D θ ( r l 2 + ( r + θ ) l 1 ) ) + ( l 1 l 2 ) P θ 2 } = ξ 6 ( s a y ) .
To minimize the total integrated cost function of inventory management, the principle minors must be positive in sign. Hence, the conditions are made this way: (i) Satisfy; then, the total costs function contains a global minimum at the decision variables’ optimum value. This model is tested through numerical experiments, and at the optimal points, the optimality is tested.
Similarly, we can check for Case II and Case III. □

4. Numerical Examples

This section used data to demonstrate the application of the model. The integrated model consists of improvement over procurement, manufacturing, and distribution. The model is formulated to examine the effect of controllable production rate and permissible delay in payment policy to minimize energy and carbon emission costs. The model is reviewed by the following numerical example. The numerical data is considered from Jaggi et al. [5] and Sarkar and Chung [26].

4.1. Numerical Example 1

This numerical considers the TCP offered by the supplier to his clients is significantly less than the TCP proposed by the producer to the supplier. The data values are:
Let Cp = $6.7/unit/week, C1p = $2.8/unit/week, C2p = $990/set up, C’1p = $0.2/unit/week, C’p = $0.3/unit/week, C’2p = $10/set up, θ = 0.02, r = 0.01, W = 100 units/week, f = 20 units/week, z = 0.2, τ = 0.03, α = 0.11, p = $5/unit, c = $11/unit, M = 12week, Pc = 220 units, μ = $90/week, k = $3500/week, s = $0.01/week, β = $0.04/week, Cs = $5.5/unit/week, C1s = $2.7/unit/week, C2s = $298/order, C’1s = $0.3/unit/week, C’s = $0.5/unit/week, C’2s = $2/order, n = 3, r = 0.01, ie = 0.15/year, ip = 0.10/year.
The best possible value of TCS1 is $11,356.2, P is 182.436 units, N is 8.995 weeks, and T is 10.1445 weeks. T1 and T2 can be found out from the Equation (7) and T2 = T/n. Complete numerical is analyzed with the help of Mathematica software.
Using the above values of all parameters and optimal values of decision parameters, we can validate our numerical through the Hessian matrix.
2 ( T C S ) T 2 = 31.5807 > 0 , | ξ 3 ξ 4 ξ 4 ξ 2 | = | 31.581 1.439 1.439 0.006105 | = 2.123 > 0 a n d | | ξ 3 ξ 4 ξ 5 ξ 4 ξ 2 ξ 6 ξ 5 ξ 6 ξ 1 | | = | 31.581 1.439 25.517 1.439 0.006105 0.3925 25.517 0.3925 110.635 | = 506.112 > 0
Therefore, at the critical point, the total integrated cost TCS1 is minimum.

4.1.1. Nature of Total Cost Concerning Two Decision Variables and the Third One Is Fixed

Here, graphical representations of the convexity of the integrated cost function with respect to the pair of decision variables are shown in Figure 4a–c. The impact of the producer’s cycle length and production rate on total cost is illustrated in Figure 4a. Here, credit period N is fixed on its optimal value. The total cost function shows the convexity for two decision variables. Figure 4b shows the convexity of the objective function with respect to production rate P and credit period N when the producer’s cycle length T is fixed. Figure 4c provides the impact of variation in the producer’s cycle length T and credit period N on the complete system’s total cost. Here, production rate is fixed. It shows the convexity of the cost function.

4.1.2. Effect of the Various Parameters on the Decision Variables and the Total Integrated Cost

Table 2 shows the impact of the number of shipments on the decision variables and the total integrated cost. As the number of shipments increases, the producer’s cycle length and the integrated cost increases, and the production rate and credit period decrease. It also shows the effect of different demand rates on the decision variables and the integrated cost function. As the demand rate is credit period dependent, total cost and production rate increase, and there is a slight decrement in the producer’s cycle length. As the demand rate is selling price dependent, then the producer’s cycle length and integrated cost increases, and the production rate and credit period decrease. An increase in the value of all energy and carbon emission parameters increases the integrated total cost function.
Figure 4d shows the impact of deterioration rate θ on an integrated cost function; as the deterioration rate increases, the total cost increases accordingly.

4.2. Numerical Example 2

This numerical considers the TCP offered by the producer to the supplier is more than the cycle length of the supplier and the TCP provided by the supplier to his customers. The data values are:
Let Cp = $4.5/unit/week, C1p = $2.1/unit/week, C2p = $1000/set up, C1p = $0.9/unit/week, C’p = $1.5/unit/week, C2p = $30/set up, θ = 0.03, r = 0.01, W = 100 units/week, f = 20 units/week, z = 0.2, τ = 0.03, α = 0.11, p = $6/unit, c = $10/unit, M = 30week, Pc = 220 units, μ = $90/week, k = $3200/week, s = $0.02/week, β = $0.03/week, Cs = $4.5/unit/week, C1s = $2.1/unit/week, C2s = $290/order, C1s= $0.9/unit/week, C’s = $0.4/unit/week, C2s = $11/order, n = 2, r = 0.01, ie = 0.15/year, ip = 0.10/year.
The best possible value of TCS2 is $146,438, P is 12,849.9 units, N is 13.0415 weeks, and T is 22.9585 weeks. T1 and T2 values can be found out from the Equation (7) and T2 = T/n.
By using the above values of all parameters and optimal values of decision parameters, we can validate our numerical through Hessian matrix.
ξ 3 = 2173.69 > 0 , | ξ 3 ξ 4 ξ 4 ξ 2 | = | 2173.69 209 209 201.526 | = 450280 > 0 a n d | ξ 3 ξ 4 ξ 5 ξ 4 ξ 2 ξ 6 ξ 5 ξ 6 ξ 1 | = | 2173.69 209 16.7189 209 201.526 30.24 16.7189 30.24 2238.89 | = 1076410000 > 0
Therefore, at the critical point, the total integrated cost TCS2 is minimum.

4.2.1. Nature of Total Cost with Respect to Two Decision Variables and the Third One Is Fixed

Here, graphical representations of the convexity of the integrated cost function with respect to the pair of decision variables are shown in Figure 5a–c. The impact of the producer’s cycle length and production rate on the total cost is illustrated in Figure 5a, here credit period N is fixed. It shows the convexity for two decision variables. Figure 5b shows the convexity of the objective function for production rate P and credit period N when the producer’s cycle length T is fixed. Figure 5c provides the impact of variation in the producer’s cycle length T and credit period N on the integrated cost. This shows the convexity with respect to the pair of decision variables.

4.2.2. Effect of the Various Parameters on the Decision Variables and the Total Integrated Cost

Table 3 shows the impact of the number of shipments, different demand rates, and cost parameters of energy and carbon emission on the decision variables and the integrated cost. As the number of shipments increases, the producer’s cycle length and the integrated cost increases, and the production rate and credit period decrease. As the demand rate is credit period dependent, total cost and production rate increase, and there is a slight decrement in the producer’s cycle length. As the demand rate is selling price-dependent, the producer’s cycle length and integrated cost increases, and the production rate and credit period decrease. An increase in the value of all energy and carbon emission parameters increases the integrated total cost function.

4.2.3. Effect of the Deterioration Rate and Selling Price on the Total Integrated Cost

Figure 5d shows that as the deterioration rate increases, the total cost increases accordingly. Figure 5e shows that as the selling price increases, demand decreases, holding cost increases, and the total cost increases.

5. Discussion and Managerial Insights

  • The result presents the strategy and critical implication for producers in managing inventory systems under agile manufacturing to incorporate energy and carbon emissions costs. The producer should also focus on the deterioration and demand rates, as it is susceptible to the total cost. Since the deterioration rate increases, the total cost of the model increases. In numerical 1 and 2, it has been proved that if the demand rate changes, then the total cost is very much affected, which is shown in Table 2 and Table 3, respectively. Here the total cost is analyzed in three cases. One is demand depends on selling price and credit period. Second, demand depends on the credit period only, and the third one is that demand depends on the price of sale only. It is observed that the total cost is minimum in case one. Thus, it is suggested to the manager to balance the different studies regarding the demand parameter.
  • The producer can adjust the production rate during production to meet the customer’s demand since it is a decision variable. The producer must be vigilant about the demands of the market. This model gives the fundamental insight for planning his credit period as per the production to fulfill all the requirements.
  • The inventory model can benefit from creating decisions to progress the inventory system by incorporating energy and carbon emission costs. The decision-maker should take care of energy and carbon emission costs while modeling as different rules impose regularity authorities.
  • To acquire more demand, TLTCP is used, due to which, both producer and supplier can earn more interest. Different cases of credit periods are shown through numerical examples. The results revealed that the TCP offered to the customer has a good effect on demand, reducing the total integrated cost. Offering TCP is the supply chain’s motivational tool, so the decision-maker should take care of this issue.

6. Conclusions

Energy and carbon emission is inevitable in every manufacturing system. Machining operations, chemical reactions, storage, and transportation are the main activities that emit carbon significantly. In the present paper, the author analyzes a production-inventory model under energy and carbon emission. The study formulated a two-level integrated model containing a producer and a supplier. The production process is assumed to be agile to depict a more realistic production model. Here production rate is considered to be a decision variable. This model has been developed based on the inventory management theory and examined how all optimal decision variables and the supply chain’s total cost are affected by critical parameters. Our computational results supply managerial insights to the production and marketing managers to help plan the supply chain’s successful coordination. The problem addressed in this model is applied to a wide array of production practices and motivated by real business concerns. Future research can be extended with a fuzzy environment for a multi-echelon integrated model with additive or flexible manufacturing to make zero waste for saving the environment.

Author Contributions

Conceptualization, V.; methodology, B.S. and S.R.S.; software, V.; validation, B.S.; formal analysis, D.Y. and M.S.; investigation, B.S. and D.Y.; resources, B.S. and M.S.; data curation, B.S. and V.; writing—original draft preparation, V.; writing—review and editing, B.S. and V.; visualization, V. and S.R.S.; supervision, B.S. All authors have read and agreed to the published version of the manuscript.

Funding

There is no external funding for this study.

Institutional Review Board Statement

Not Applicable.

Informed Consent Statement

Not Applicable.

Data Availability Statement

Not Applicable.

Acknowledgments

The authors are grateful to the Editor of the journal and reviewers of the paper.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

b = k P + s P + β ( P P c ) + μ , l 1 = l ( e T θ 1 ) D 1 , l 2 = m ( e T θ 1 ) D 1 , K = e r N e r M , l = e ( e T θ 1 ) ( r + θ ) D P θ , m = e ( e T θ 1 ) r D P θ ,
K 1 = ( e ( θ + r )     M e ( N + T / n ) ( θ + r ) ) , K 2 = ( e r M e ( N + T / n ) r ) , K 3 = ( e θ T / n e r T / n ) , K 4 = ( e r ( N + T / n ) + e ( N + T / n ) ( θ + r ) + θ T / n ) , D 1 = ( W f ) ( 1 z ) N log ( 1 z ) ,
l 1 = 1 2 P D T 2 θ 2 2 P 2 + D 2 T 2 θ 2 D P T 2 θ 2 , l 2 = ( 2 D 2 T θ 2 2 D P T θ 2 ) ( 2 P D T 2 θ 2 ) 4 ( P 2 + D 2 T 2 θ 2 D P T 2 θ 2 ) 3 / 2 + D T θ 2 P 2 + D 2 T 2 θ 2 D P T 2 θ 2 ,
l 3 = r l 1 P θ + r T 1 P , l 4 = ( 2 P D T 2 θ 2 ) 2 4 ( P 2 + D 2 T 2 θ 2 D P T 2 θ 2 ) 3 / 2 1 P 2 + D 2 T 2 θ 2 D P T 2 θ 2 , l 6 = r l 4 P θ + 2 r l 1 P 2 θ 2 r T 1 P 2 , l 5 = ( D θ l ) P θ [ 1 + ( D θ l ) 2 D 2 T θ 2 D P T θ 2 ] , l 7 = r l P θ + θ ( 1 l P θ ) .

References

  1. Cárdenas-Barrón, L.E.; Sana, S.S. A production-inventory model for a two-echelon supply chain when demand is dependent on sales teams׳ initiatives. Int. J. Prod. Econ. 2014, 155, 249–258. [Google Scholar] [CrossRef]
  2. Taleizadeh, A.A.; Kalantari, S.S.; Cárdenas-Barrón, L.E. Pricing and lot sizing for an EPQ inventory model with rework and multiple shipments. TOP 2016, 24, 143–155. [Google Scholar] [CrossRef]
  3. Taleizadeh, A.A.; Khanbaglo, M.P.S.; Cárdenas-Barrón, L.E. An EOQ inventory model with partial backordering and reparation of imperfect products. Int. J. Prod. Econ. 2016, 182, 418–434. [Google Scholar] [CrossRef]
  4. Chen, X.; Wang, X.; Kumar, V.; Kumar, N. Low carbon warehouse management under cap-and-trade policy. J. Clean. Prod. 2016, 139, 894–904. [Google Scholar] [CrossRef] [Green Version]
  5. 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]
  6. Khanna, A.; Kishore, A.; Sarkar, B.; Jaggi, C. Supply chain with customer-based two-level credit policies under an imperfect quality environment. Mathematics 2018, 6, 299. [Google Scholar] [CrossRef] [Green Version]
  7. Aljazzar, S.M.; Gurtu, A.; Jaber, M.Y. Delay-in-payment—A strategy to reduce carbon emissions from supply chains. J. Clean. Prod. 2018, 170, 636–644. [Google Scholar] [CrossRef]
  8. Tiwari, S.; Ahmed, W.; Sarkar, B. Multi-item sustainable green production system under trade-credit and partial backordering. J. Clean. Prod. 2018, 204, 82–95. [Google Scholar] [CrossRef]
  9. Sarkar, B. Mathematical and analytical approach for the management of defective items in a multi-stage production system. J. Clean. Prod. 2019, 218, 896–919. [Google Scholar] [CrossRef]
  10. Tayyab, M.; Sarkar, B. Optimal batch quantity in a cleaner multi-stage lean production system with random defective rate. J. Clean. Prod. 2016, 139, 922–934. [Google Scholar] [CrossRef]
  11. Kim, M.S.; Sarkar, B. Multi-stage cleaner production process with quality improvement and lead time dependent ordering cost. J. Clean. Prod. 2017, 144, 572–590. [Google Scholar] [CrossRef]
  12. Ahmed, W.; Sarkar, B. Impact of carbon emissions in a sustainable supply chain management for a second generation biofuel. J. Clean. Prod. 2018, 186, 807–820. [Google Scholar] [CrossRef]
  13. Sarkar, M.; Chung, B.D. Flexible work-in-process production system in supply chain management under quality improvement. Int. J. Prod. Res. 2019, 58, 3821–3838. [Google Scholar] [CrossRef]
  14. Sana, S.; Chaudhuri, K.S. On a volume flexible stock dependent inventory model. Adv. Model. Optim. 2003, 5, 197–210. [Google Scholar]
  15. Sana, S.; Chaudhuri, K.S. On a volume flexible production policy for a deteriorating item with time dependent demand and shortage. Adv. Model. Optim. 2004, 6, 57–73. [Google Scholar]
  16. Singh, S.R.; Singh, A.P.; Bhatia, D. A supply chain model with variable holding cost for flexible manufacturing system. Int. J. Oper. Res. Optim. 2010, 1, 107–120. [Google Scholar]
  17. Singh, S.R.; Gupta, V.; Bansal, P. EOQ model with volume agility, variable demand rate, Weibull deterioration rate and inflation. Int. J. Comput. Appl. 2013, 72, 1–6. [Google Scholar]
  18. Dem, H.; Singh, S.R. Joint replenishment modelling of a multi-item system with greening policy and volume flexibility. Int. J. Oper. Res. 2015, 22, 148–166. [Google Scholar] [CrossRef]
  19. Sangal, I.; Gupta, V. Vendor-supplier cooperative inventory model with two warehouse and variable demand rate under the environment of volume agility. In Proceedings of the Fifth International Conference on Soft Computing for Problem Solving; Springer: Singapore, 2016; pp. 109–123. [Google Scholar]
  20. Singh, S.R.; Gupta, V. Vendor–buyer model with error in quality inspection and selling price dependent demand rate under the effect of volume agility. Int. J. Oper. Quant. Manag. 2016, 22, 357–371. [Google Scholar]
  21. Gautam, P.; Kamna, K.M.; Jaggi, C.K. Sustainable production policies under the effect of volume agility, preservation technology and price- reliant demand. Yugosl. J. Oper. Res. 2020, 30, 307–324. [Google Scholar] [CrossRef]
  22. Ouyang, L.Y.; Teng, J.T.; Chuang, K.W.; Chuang, B.R. An optimal inventory policy with non instantaneous receipt under trade-credit. Int. J. Prod. Econ. 2005, 98, 290–300. [Google Scholar] [CrossRef]
  23. Jaggi, C.K.; Goyal, S.K.; Goel, S.K. Retailer’s optimal replenishment decisions with credit-linked demand under permissible delay in payments. Eur. J. Oper. Res. 2008, 190, 130–135. [Google Scholar] [CrossRef]
  24. Sarkar, B. An EOQ model with delay in payments and time varying deterioration rate. Math. Comput. Mod. 2012, 55, 367–377. [Google Scholar] [CrossRef]
  25. Khanra, S.; Mandal, B.; Sarkar, B. An inventory model with time dependent demand and shortages under trade-credit policy. Econ. Model. 2013, 35, 349–355. [Google Scholar] [CrossRef]
  26. Chung, K.J. The EPQ model under conditions of two levels of trade-credit and limited storage capacity in supply chain management. Int. J. Syst. Sci. 2013, 44, 1675–1691. [Google Scholar] [CrossRef]
  27. Sarkar, B.; Sana, S.S.; Chaudhuri, K. An inventory model with finite replenishment rate, trade-credit policy and price-discount offer. J. Ind. Eng. 2013, 2013, 672504. [Google Scholar] [CrossRef]
  28. Sarkar, B.; Gupta, H.; Chaudhuri, K.; Goyal, S.K. An integrated inventory model with variable lead time, defective units and delay in payments. Appl. Math. Comput. 2014, 237, 650–658. [Google Scholar] [CrossRef]
  29. Sarkar, B.; Saren, B.; Cárdenas-Barrón, L.E. An inventory model with trade-credit policy and variable deterioration for fixed lifetime products. Ann. Oper. Res. 2015, 229, 677–702. [Google Scholar] [CrossRef]
  30. Yadav, D.; Singh, S.R.; Kumari, R. Retailer’s optimal policy under inflation in fuzzy environment with trade credit. Int. J. Syst. Sci. 2015, 42, 754–762. [Google Scholar] [CrossRef]
Figure 1. Behavior of producer’s inventory level versus time.
Figure 1. Behavior of producer’s inventory level versus time.
Energies 14 01569 g001
Figure 2. The behavior of the supplier’s inventory level versus time.
Figure 2. The behavior of the supplier’s inventory level versus time.
Energies 14 01569 g002
Figure 3. (a) Revenue level versus time. (b) Revenue level versus time for the accumulated inventory. (c) Revenue level versus time.
Figure 3. (a) Revenue level versus time. (b) Revenue level versus time for the accumulated inventory. (c) Revenue level versus time.
Energies 14 01569 g003
Figure 4. (a) Nature of total cost with respect to producer’s cycle length and production rate. (b) Nature of total cost with respect to credit period and production rate. (c) Nature of total cost with respect to producer’s cycle length and credit period. (d) Integrated cost for the products versus deterioration rate.
Figure 4. (a) Nature of total cost with respect to producer’s cycle length and production rate. (b) Nature of total cost with respect to credit period and production rate. (c) Nature of total cost with respect to producer’s cycle length and credit period. (d) Integrated cost for the products versus deterioration rate.
Energies 14 01569 g004
Figure 5. (a) Nature of total cost with respect to producer’s cycle length and production rate. (b) Nature of total cost with respect to credit period and production rate. (c) Nature of total cost with respect to producer’s cycle length and credit period. (d) Integrated cost for the products versus the deterioration rate. (e) Integrated cost for the products versus selling price.
Figure 5. (a) Nature of total cost with respect to producer’s cycle length and production rate. (b) Nature of total cost with respect to credit period and production rate. (c) Nature of total cost with respect to producer’s cycle length and credit period. (d) Integrated cost for the products versus the deterioration rate. (e) Integrated cost for the products versus selling price.
Energies 14 01569 g005
Table 1. Contribution from various researchers.
Table 1. Contribution from various researchers.
Researcher (s)Model TypeManufacturing TypeFinancial PolicyDemand PatternEnvironmental Effect
Cárdenas-Barrón [1]EPQTraditionalTrade-creditVariableNot applicable
Taleizadeh et al. [2,3]EOQNot applicableNot applicableConstantNot applicable
Sarkar et al. [5]SSCMTraditionalTrade-creditConstantCarbon emission
Tiwari et al. [8]EPQTraditionalTrade-creditConstantNot applicable
Sarkar [9]SCMTraditionalNot applicableConstantNot applicable
Kim and Sarkar [11]EPQTraditionalNot applicableConstantCarbon emission
Ahmed and Sarkar [12]SSCMTraditionalNot applicableConstantCarbon emission
Sana and Chaudhuri [14]EPQAgileNot applicableVariableNot applicable
Ouyang et al. [22]EOQNot applicableTrade-creditConstantNot applicable
Jaggi et al. [23]Integrated inventoryTraditionalTrade-creditConstantNot applicable
Chung [26]EPQTraditionalTrade-creditConstantNot applicable
This modelGSCMAgileTrade-creditVariableCarbon emission
Table 2. Sensitivity analysis of key parameters of the model.
Table 2. Sensitivity analysis of key parameters of the model.
Parameter ChangesT (week)N (week)PTCS2 ($)
n28.16899.4597199.61110,867.2
310.14458.994182.43511,356.2
412.15368.7306164.18311,760.5
514.4838.5827142.53212,103.6
618.6378.6091105.99812,127.4
Demand rateD = W − (Wf)(1z)N10.13428.996182.62912,155.3
D = W + α p−τ18.79610.4741984.626514,451.6
D = W − (Wf)(1z)N ++ α p−τ10.14458.995182.43611,356.2
C1p0.210.14458.995182.43611,356.2
0.39.9208.950182.02311,374.1
0.49.74218.912181.62911,396.3
0.59.53248.88181.14711,405.2
0.69.24568.85180.84311,422.8
Cp0.110.23129.120182.3611,356.0
0.210.14458.995182.43611,356.1
0.310.01238.812182.14811,356.2
0.49.82458.723181.96211,356.3
0.59.73108.621181.82311,356.4
C2p810.14218.992182.42111,356.0
1010.14458.995182.43611,356.2
1210.14628.999182.44611,356.4
1410.14839.002182.45911,356.6
1610.14949.014182.46311,356.8
Table 3. Sensitivity analysis of key parameters of the model.
Table 3. Sensitivity analysis of key parameters of the model.
Parameter ChangesT (week)N (week)PTCS2 ($)
n222.958513.041512,849.9146,438
328.295112.50411,061.8153,318
431.598512.216811,052.3170,615
533.956112.032411,497.2190,263
636.178511.880911,650.9206,448
Demand rateD = W − (Wf)(1z)N22.870613.05212,991.3147,304
D = W + α pτ6.17274.841238.7168,514
D = W − (Wf)(1z)N ++ α pτ22.958513.041512,849.9146,438
C1p0.623.710113.64112,906.1146,369
0.723.524513.61012,884.3126,401
0.823.012413.37012,860.1146,438
0.922.955913.04112,849.9146,461
1.022.515712.88012,812.3146,492
Cp1.323.341013.4412,852.3146,425
1.423.012013.2812,851.2146,435
1.522.951113.0412,849.9146,438
1.622.814712.8412,847.0146,440
1.722.631212.5112,846.6146,443
C2p2622.95.312.31512,818.4146,420
2822.957412.71212,824.1146,426
3022.958513.04112,849.9146,438
3222.959113.51212,855.6146,445
3422.960014.14312,901.2146,450
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Vandana; Singh, S.R.; Yadav, D.; Sarkar, B.; Sarkar, M. Impact of Energy and Carbon Emission of a Supply Chain Management with Two-Level Trade-Credit Policy. Energies 2021, 14, 1569. https://doi.org/10.3390/en14061569

AMA Style

Vandana, Singh SR, Yadav D, Sarkar B, Sarkar M. Impact of Energy and Carbon Emission of a Supply Chain Management with Two-Level Trade-Credit Policy. Energies. 2021; 14(6):1569. https://doi.org/10.3390/en14061569

Chicago/Turabian Style

Vandana, S. R. Singh, Dharmendra Yadav, Biswajit Sarkar, and Mitali Sarkar. 2021. "Impact of Energy and Carbon Emission of a Supply Chain Management with Two-Level Trade-Credit Policy" Energies 14, no. 6: 1569. https://doi.org/10.3390/en14061569

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