Next Article in Journal
Sample Size Determination for Two-Stage Multiple Comparisons for Exponential Location Parameters with the Average
Previous Article in Journal
Backpropagated Neural Network Modeling for the Non-Fourier Thermal Analysis of a Moving Plate
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Smart Production System with Random Imperfect Process, Partial Backordering, and Deterioration in an Inflationary Environment

1
Department of Mathematics, Vardhaman College, Bijnor 246701, UP, India
2
Department of Mathematics, Maharaj Singh College, Saharanpur 247001, UP, India
3
Department of Mathematics, DN College, Meerut 250002, UP, India
4
Department of Industrial Engineering, Yonsei University, 50 Yonsei-ro, Sinchon-dong, Seodaemun-gu, Seoul 03722, Republic of Korea
5
Center for Transdisciplinary Research (CFTR), Saveetha Dental College, Saveetha Institute of Medical and Technical Sciences, Saveetha University, 162, Poonamallee High Road, Velappanchavadi, Chennai 600077, TN, India
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(2), 440; https://doi.org/10.3390/math11020440
Submission received: 16 October 2022 / Revised: 4 December 2022 / Accepted: 7 January 2023 / Published: 13 January 2023
(This article belongs to the Section Engineering Mathematics)

Abstract

:
In today’s digital age, industrial methods are shifting away from humans and toward machines. We choose automated systems for various jobs related to production systems, such as screening, manufacturing. A smart manufacturing system is one in which machines take the place of humans. Under the influence of inflation, this study proposes a smart production-inventory model with partial backlogging, and an imperfect manufacturing process where the deterioration rate is constant. Every production system, in reality, has a random defect rate. A screening procedure is required due to the manufacture of some defective items, which is carried out by machine, i.e., by an automated system. Carbon is released during the manufacturing process due to actions such as holding deterioration. As a result, carbon emissions are taken into account in the current study. The goal of this study is to reduce total inventory costs as much as possible. To demonstrate the proposed model’s practical application, many numerical examples and sensitivity assessments with graphs are provided.

1. Introduction

Carbon emissions are a major contributor to global warming. Today, carbon emissions are a serious global issue, and many governments are taking considerable steps to reduce them. During the 1990s, several international meetings and conferences were held to try to discover a solution to this problem. In 1997, an international agreement known as the Kyoto Protocol was signed to reduce greenhouse gases and carbon dioxide in the atmosphere. The carbon Cap and Trade program, which is particularly useful in cutting carbon emissions, is one of the most effective methods of the Kyoto Protocol. This system sets a carbon emissions cap for all industries and allows them to buy and sell carbon emission rights within that cap. Industry (www.epa.gov/ghgemissions, accessed on 1 March 2022), transportation, and electricity production are the key sources of carbon emissions, according to the United States (US) Environmental Protection Agency. In 2019, the industry contributed to 23% of total greenhouse gas emissions, transportation accounted for 29%, and power generation accounted for 25% of all greenhouse gas emissions in the US (Figure 1).
Carbon emissions are caused by automobiles, industries, power generation, and so other industries. Industries are directly linked to various economies from an economic standpoint. As a result, while some international policies and agreements can help reduce carbon emissions, a complete removal is not conceivable. At this time and for profitability, a smart production system is required. Because of the high level of market competition nowadays, every firm strives to produce error-free products to improve their brand value, thus inspection is necessary for the production system. Automation is a smart method of maximizing the use of machine resources. Machine inspection, rather than human inspection, ensure errorless items and transforms the manufacturing system into a smart manufacturing system.
Because of increased concerns among stakeholders such as international organizations, governments, supply chain partners, and consumers about the environmental impact due to the manufacturing, shipping, storing, reworking, and disposal of products, many companies are investing in carbon emission reduction. To address these environmental concerns, businesses have implemented or are considering investments in carbon emissions reduction in a variety of ways, including the use of energy-efficient equipment and facilities, renewable energy resources, product and packaging redesign, and process optimization. All over the world, different nations opt for the carbon tax and cap mechanism to achieve this; the term cap is used to describe a limit on the total quantity of carbon that can be emitted by a firm in the system.
To achieve the task of meeting sustainable development goals, knowing the actual picture of the financial position of the business is very important. In developing countries such as Indonesia, India, Iraq, the role of inflation is very crucial. In general, inflation refers to a rise in the price level of an economy over time (or price inflation). When the price level rises, each unit of currency buys fewer goods and services. In other ways, it can says that inflation reduces the purchasing power of a currency. The following sections are included in the research work: Section 1 contains the introduction. Section 2 includes the literature survey and research gaps. Section 3 consists of the assumptions and notation, and in Section 4, the formulation of the mathematical model is described. The solution methodology is discussed in Section 5. Section 6 contains numerical analysis and Section 7 consists of marginal insight, and finally in Section 8, concluding remarks with future scope are given. The flow chart is given below in Figure 2.

2. Literature Review

2.1. Inventory Models with Carbon Emissions

Global warming is mostly caused by carbon emissions. Some governments are focusing their efforts on reducing carbon emissions right now. The main causes of greenhouse gas emissions that contribute to global warming are human activities such as burning fossil fuels for industrial purposes, generating electricity, and driving automobiles. From a global economic standpoint, the industry is intrinsically related to the economy. Hua et al. [1] examined carbon emissions and how corporations manage carbon emissions through carbon emission trading schemes. Carbon emission metrics were linked to several decision variables in a model proposed by Benjaafar et al. [2]. They looked at how emissions control policies affected prices and emissions. The carbon cap-and-trade mechanism was considered by Lou et al. [3] in the study of a two-stage supply chain policy, and in addition to this, investment was made for opting for green technology. In order to maximize profit and curb carbon emissions, Datta [4] designed an inventory model by considering green technology and carbon tax regulations. An economic production quantity (EPQ) model was designed by Daryanto et al. [5], incorporating deterioration as well as a percentage of defective items in the model. They used the concept of a secondary market for defective things with a discount price to see the environmental issues of carbon emission costs, and carbon cap and tax policies are included in its decision model. Jauhari et al. [6] investigated a system of a closed supply chain in which they took a single producer, single retailer, and single collector. Green technologies, faulty production, and remanufacturing procedures were all taken into consideration. Mishra et al. [7] explored a carbon pricing and cap mechanism for a sustainable supply chain management. They looked at energy efficient products with a green technology investment. Yadav et al. [8] examined the reduction of carbon emissions through the selection of items with cross-price elasticity of demand to form a sustainable supply chain. Sarkar et al. [9] discussed carbon tax policy in a sustainable managerial decision-making problem for a substitutable product in a dual channel.

2.2. Inventory Models with Backlogging

Backlog refers to the inability to meet demand quickly from stock. The consumer is expected to wait until the demand is met, which may take some time. In this period, some customers may choose other sources to fulfill demand, so there is a partial need for the present scenario. Considering time-varying backlogging, Sarkar et al. [10] calcalculated inventory within a smart production model considering Stackelrberg game policy. Wee et al. [11] considered two partial backordering costs. They considered linear and fixed backordering costs. Considering stock-and-price-dependent demand, Mishra et al. [12] investigated a controllable deterioration rate inventory model with shortages. Two cases of backordering were discussed: (i) partial backordering and (ii) full backordering. Preservation technology is also being considered. A supply chain model with reverse logistics was developed by Singh and Rani [13] with partially backordered items. They also considered inflation, defective production, and remanufacturing. A model in which demand is taken as stock-and-price-dependent with time-varying holding costs, partial backlog, and a quantity discount was described by Palanivel and Suganya [14]. Further, a concept of partial backlogging was introduced by Kumar et al. [15] in two warehouse inventory models with stock-dependent demand.

2.3. Inventory Models with Random Imperfect Production Rate

It is common knowledge that any imperfect real-world production system has a random imperfect rate. This research takes a step in that approach by extending an inventory model to include random faulty rates. A supply chain with a random demand rate for different warehouse models were explored by Sarkar et al. [16], considering a single-stage production process. Jawla and Singh [17] gave a multi-item and imperfect-item economic production quantity model. They considered a multi-setup for production and rework, and also discussed a learning environment with preservation technology. Khara et al. [18] designed an inventory model considering the influence of imperfect manufacturing on cost. Moreover, demand was taken as reliability-dependent. Considering imperfect production systems, Al-Salamah [19] developed an EPQ model with adjustable rework rates that are both synchronous and asynchronous. In asynchronous production, defective things are held until the entire production lot is finished, but in synchronous production, as soon as the defective goods are generated, they are reworked. With a multiple transportation facility, Sarkar et al. [20] explored a supply chain model for biodegradable products. In a similar direction, Gupta et al. [21] discussed multiple distribution network with integrated policy. Lin [22] developed an EPQ model with uncertain demand, an imperfect rework process, and backlogging. He assumed that imperfect items were stocked separately and then reworked.

2.4. Inventory Models with Inflation

If money was fully neutral, inflation would have no influence on the real economy. Nevertheless, full neutrality is a concept that is not universally accepted, so consideration of inflation is necessary. The impacts of the time value of money and inflation on inventories were explored by Chandra and Michael [23]. For external and internal costs, they assumed two different inflation rates. The effect of inflation was explored by Yang et al. [24] on an inventory model where partial backlogging, stock-dependent demand, and deterioration were considered. Palanivel and Uthayakumar [25] investigated the influence of the time value of money and inflation on an EPQ model with probabilistic deterioration, partial backordering, and variable production. They talked about three different kinds of continuous probabilistic deterioration functions. A green supply chain (GSC) model was developed by Rani et al. [26] and explored the effect of inflation on it. They considered that for remanufacturing products, a secondary market is available. Huang et al. [27] created a model for perishable food with replenishment policies under inflation. Under inflation, an economic order quantity (EOQ) model was designed by Alamri et al. [28] for deteriorating items with carbon emissions. Padiyar et al. [29] discussed the inventory model in an inflationary environment with an imperfect production process and preservation technology (Table 1).

2.5. Research Gap and Research Contributions

From the literature, it is observed that some researchers considered the reworking process, volume flexibility, shortages, inflation, random imperfect processes, or carbon emission independence, or considered two factors or three factors at the same time. As per our knowledge, no researcher has analyzed all these factors together while modelling the production-inventory model. The current work focuses on the following:
  • To make goodwill in the market, screening items is necessary before satisfying the demand of customers. An attempt is made to develop a smart production system by considering that the screening process is performed with the help of an automated machine system.
  • Due to the production process, due to storage, and other different operational activities, carbon emission is unavoidable. To make the model environment sustainable, carbon emissions due to different activities associated with manufacturing systems have been considered.
  • Further, obtaining the goal of sustainability is achieved with help of reworking and proper waste management.
  • The benefits of a flexible manufacturing system are explored by considering the production rate-dependent production cost.
  • The fluctuating nature of the market is absorbed into the model by considering inflation.
  • The model is analyzed under the effect of partial backlogging.

3. Assumptions, Notation, and Problem Description

For modelling purposes, the required assumptions and notation are presented in this section.

3.1. Assumptions

  • Nowadays, production systems have gradually shifted from manual to the machine. Thus, in the current study, we consider a smart production system with an automated screening process. Moreover, we have considered a single-type product inventory.
  • To avoid overstocking and understocking in the production system, it is very crucial to tune up the production rate with the demand for the product. Thus, the production rate p is variable (due to the smart production system).
  • Demand rate d is constant and known.
  • Due to the presence of inventory in the stock, deterioration is unavoidable. Thus, it is important to consider its effect on the model as it has some economical value. Thus, in the current study, the rate of deterioration of inventory is considered as the constant parameter θ, where 0 < θ < 1 .
  • Practically, it is observed that the production cost consists of different components such as raw material cost, tool cost, labor cost, energy cost, etc., and depends on the production rate. Thus, the production cost is considered as follows:
    c p = C 0 + C 1 p + C 2 / p ,
    where C 0 , C 1 , C 2 are raw material cost, labor/energy charges, and tool/die cost, respectively.
  • Shortages are permitted and partially backlogged, i.e., a portion of the demand is backlogged and has been taken as B(t) = e ε t , where t is the time of waiting and ε > 0 is the backlogging parameter.
  • The rate of imperfect production is random and follows the known probability distribution.
  • A portion of the defective items will be reworkable, while the remainder will be scrapped immediately if they cannot be reworked successfully.
  • To absorb the market disturbance, the whole of the study is analyzed under the effect of inflation.
  • Due to the different operational activities associated with the production system, the emission of carbon is unavoidable. Thus, to control carbon emissions, a regulatory mechanism is considered.

3.2. Notation

Decision Variables
ParametersUnitExpression
p units/unit timeProduction rate
t 1 yearsTime at which production starts
t 3 yearsTime at which production stop
Model Parameters
t 2 yearsTime at which backordered cleared partially
t 4 yearsRework process starting time
t 5 yearsRework process ending time
t 6 yearsWhole length of a complete cycle
dunits/unit timeDemand rate
θunits/unit timeDeterioration rate
p r units/unit timeRework process rate
runitDiscount rate
iunitRate of inflation
ε   >0 Shape parameter of backorder rate
s p USD per setupSetup cost for production
s r USD per setupSetup cost for reworking station
c p USD/unitProduction cost
c θ USD/unitDeterioration cost
c r USD/unitCost of rework
c b USD/unit/unit timeBackorder cost
c s USD/unitScrapping cost
c L USD/unitLost sale cost
C 0 USD/unitRaw material cost
C 1 USD/unitLabor/energy charges
C 2 USD/unitTool/die cost
h s USD/unit/unit timeServiceable items holding cost
h r USD/unit/unit timeReworkable items holding cost
e p kg/unitCarbon emissions by production
e s p kg/setupCarbon emissions by setup process of production
e s r kg/setupCarbon emissions by setup process of reworking station
e h kg/unitCarbon emissions by holding the items in a warehouse
e r kg/unitCarbon emissions due to reworking
e θ kg/unitCarbon emissions by deterioration
e s kg/unitCarbon emissions by scrapping
δ USD/kgCarbon tax per kg
ZUSDTotal cost
Z T USD/cycleTotal cost per cycle
I i t unitServiceable inventory level at time t, i = 1, 2, 3, 4, 5, 6
I j t unitReworkable inventory level at time t, j = 7, 8, 9
x p -Defective proportion in production
x r -Defective proportion in rework
E x p -Expected value of defective proportion in production
E x r -Expected value of defective proportion in rework

3.3. Problem Descriptions

In this section, a cleaner production inventory model is developed under the effect of inflation. Waste management and reworking are considered in the model to reduce the waste from the system as much as possible. Because of the environmentally friendly nature of the decision maker, government regulation is imposed on carbon emissions. The inventory situation of the problem is described in Figure 3 and Figure 4. Figure 3 represents the inventory situation of serviceable items, while Figure 4 represents the inventory situation of non-serviceable items. From Figure 3, it is observed that the production process started at the time t 1 and over the interval [ t 1 ,   t 2 ], backlogging was covered partially, and in the interval [ t 2 , t 3 ], the inventory level increased due to production, demand, and deterioration. The production process stopped at the time t 3 , and in the interval [ t 3 , t 4 ] , the inventory level decreased due to demand and deterioration. Rework started at the time t 4 , and in the interval [ t 4 , t 5 ], the inventory level increased due to rework, demand, and deterioration, and decreased in the interval due to demand and deterioration.
From Figure 4, it is observed that over the interval [ t 1 ,   t 3 ], the inventory level increased due to the random imperfectness of the production process, and in the interval [ t 3 , t 4 ] , the inventory level decreased due to the deterioration of imperfect items, and over the interval [ t 4 , t 5 ] , the inventory level decreased due to deterioration and rework, and reached zero at the time t 5 .

4. Mathematical Model of Smart Production System

Initially, the inventory cycle starts with shortages in the interval [0, t 1 ], and at the time t 1 , a fraction of the shortage is cleared by the backlogging rate B(t). In the interval [ t 1 ,   t 2 ], the inventory grows with (production, demand, and deterioration), in [ t 2 , t 3 ], the inventory decreases by (demand and deterioration), in [ t 3 , t 4 ] with (rework, demand, and deterioration), and in the interval [ t 4 , t 5 ], the inventory decreases (demand and deterioration).
The inventory level declines due to the production, rework, demand, and deterioration rates. The inventory level of serviceable inventories can be represented mathematically with the help of the following differential equations:
d I 1 t d t   = d B t 1 t   = d e ε t 1 t , 0     t     t 1 ;   I 1 0 = 0
d I 2 t d t =   p { 1 E [ x p ] } d , t 1   t   t 2 ,   I 2 t 2 = 0
d I 3 t d t   = p { 1 E [ x p ] } θ I 3 t d , t 2     t     t 3 , I 3 t 2 = 0
d I 4 t d t   = d θ I 4 t , t 3     t     t 4 , I 4 t 4 = 0
d I 5 t d t   = p r { 1 E x r } d θ I 5 t , t 4     t   t 5 , I 5 t 4 = 0
d I 6 t d t = d θ I 6 , t 5     t     t 6 , I 6 t 6 = 0 .
The above-mentioned differential equations’ solutions are expressed as
I 1 t   = d ε e ε t 1   [ 1 e ε t ] , 0     t     t 1
I 2 t = [ p { 1 E [ x p ] } d ] t t 2 , t 1   t   t 2
I 3 t = 1 θ [ 1 e θ t t 2 ] [ p { 1 E [ x p ] } d ] , t 2     t     t 3 ,
I 4 t   = d θ [ e θ ( t 4 t ) 1 ] , t 3     t     t 4 ,
I 5 t = 1 θ [ 1 e θ t t 4 ] p r 1 E x r d , t 4     t   t 5 ,
I 6 t   = d θ [ e θ ( t 6 t ) 1 ] , t 5     t     t 6 .
Similar to the above, the inventory level of reworkable inventories is governed by the following differential equations.
d I 7 t d   t   = p E [ x p ] θ I 7 , t 1   t   t 3 , I 7 t 1 = 0
d I 8 t d t   = θ I 8 , t 3     t     t 4   , I 7 t 3 = I 8 t 3
d I 9 t d t   = p r θ I 9 , t 4     t     t 5 I 9 t 5 = 0 .
The above-mentioned differential equations’ solutions are expressed as
I 7 t = 1 θ [ 1 e θ t t 1 p E [ x p , t 1   t   t 3 ,
I 8 t   = 1 θ e θ t   [ e θ t 3 e θ t 1 p E [ x p , t 3     t     t 4   ,
I 9 t   = p r θ [ e θ ( t 5 t ) 1 ] , t 4     t     t 5 .
Using continuity I 1 t 1 = I 2 t 1 , I 3 ( t 3 ) = I 4 t 3 , I 5 t 5 = I 6 t 5 , and I 8 t 4 = I 9 t 4 and Taylor’s series expansion and omitting higher-order terms, we obtain the values of t 2 , t 4 , t 5 , t 6 in terms of t 1 and t 3 .
Assume   a = p { 1 E [ x p ] }   and   b = p E [ x p ]   c = p r 1 E x r .   Then ,
t 2 = t 1 + d a d ε 1 e ε t 1
t 4 = t 3 + a d d ( t 3 t 2 ) [ 1 θ 2 t 3 t 2 ]
t 6   = t 5 + c d d ( t 5 t 4 ) [ 1 θ 2 t 5 t 4 ]
t 5   = t 4 + b p r ( t 3 t 1 ) [ 1 + θ 2 t 3 + t 1 2 t 4 ] .
Now, different costs can be evaluated step by step as follows:
Production   cost   per   cycle   ( C P ) = c p t 1 t 2 p e R t d t + c p t 2 t 3 p e R t d t   = c p p R e R t 1 e R t 3 .
Rework   cost   per   cycle   ( C R ) = c r t 4 t 5 p r e R t d t = c r p r R e R t 4 e R t 5 .
Setup   cost   per   cycle   ( C S p r ) = s p + s r .
Backorder   cost   ( C B ) = c b [ 0 t 1 I 1 t e R t d t + t 1 t 2 I 2 t e R t   d t ] . = c b   [ d ε e ε t 1 { 1 R ( 1 e R t 1 ) + 1 ε R ( 1 e ε R t 1 ) } + [ p { 1 E [ x p ] } d ] . [ 1 R e R t 1 t 1 t 2 + 1 R 2 ( e R t 1 e R t 2 ) ] ] .
Lost   sale   cost   per   cycle   ( C L ) = c L 0 t 1 ( 1 e ε t 1 t ) d e R t d t = c L d [ 1 R ( 1 e R t 1 e ε t 1 ε R ( e t 1 ε R 1 ) ] .
The total deterioration cost is given by
C θ = c p t 2 t 3 p { 1 E [ x p ] } e R t d t c p t 2 t 4 d e R t d t + c r t 4 t 5 p r 1 E x r e R t d t c r t 4 t 6 d e R t d t   + c p t 1 t 3 p E [ x p ] e R t d t c p t 4 t 5 p r 1 E x r e R t d t = c p p { 1 E [ x p ] } 1 R e R t 2 e R t 3 c p d 1 R [ e R t 2 e R t 4 ] + c r p r 1 E x r   1 R [ e R t 4 e R t 5 ] [ c r d 1 R e R t 4 e R t 6 ] + c p p E [ x p ]   1 R e R t 1 e R t 3 c p p r 1 E x r 1 R e R t 4 e R t 5 .
Scrapping cost per cycle (CS) = c s p r E x   r t 5 t 4   e R t 5 .
Per cycle holding cost is the sum of (i) per cycle holding cost for serviceable inventory and (ii) per cycle holding cost for reworkable inventory. Per cycle holding cost for serviceable inventory
H C s = h s [ t 2 t 3 I 3 t e R t d t + t 3 t 4 I 4 t e R t d t + t 4 t 5 I 5 t e R t d t + t 5 t 6 I 6 t e R t d t ] = h s [ { p { 1 E [ x p ] } d } θ [ 1 R ( e R t 2 e R t 3 ) e θ t 2 θ + R ( e θ + R t 2 e θ + R t 3 ) ] + d θ [ e θ t 4 θ + R ( e θ + R t 3 e θ + R t 4 ) + 1 R ( e R t 4 e R t 3 ) ] + p r 1 E x r d θ [ 1 R ( e R t 4 e R t 5 ) e θ t 4 θ + R ( e θ + R t 4 e θ + R t 5 ) ] + d θ [ e θ t 6 θ + R ( e θ + R t 5 e θ + R t 6 ) + 1 R ( e R t 6   e R t 5 ) ] ] .
Per cycle holding cost for recoverable inventory
H C r = h r [ t 1 t 3 I 7 t e R t d t + t 3 t 4 I 8 t e R t d t + t 4 t 5 I 9 t e R t d t   ] = h r [ p E [ x p ]   θ [ 1 R ( e R t 1 e R t 3 ) e θ t 1 θ + R ( e θ + R t 1 e θ + R t 3 ) ] + 1 θ [ e θ t 3 e θ t 1   p E [ x p 1 θ + R ( e θ + R t 3   e θ + R t 4 ) + p r θ [ e θ t 5 θ + R ( e θ + R t 4   e θ + R t 5 ) + 1 R ( e R t 5 e R t 4 ) ] .
Total inventory cost per cycle
T C 1 = 1 t 6 ( C P + C R + C S p r + C B + C L + C θ + C S + H C s + H C r ) .
Now, due to production, rework, and deterioration, the holding carbon emissions costs can be obtained step by step as follow.
  • Carbon emissions cost due to setup of production and reworking station
    ( C S E   ) = δ   t 6   e s p + e s r
  • Per cycle carbon emissions cost due to production
    ( C E p ) = δ   t 6 t 1 t 3 e p p e R t d t   = δ e   p p   t 6 R [ e R t 1   e R t 3 ]
  • Per cycle carbon emissions cost due to rework
    ( C E r ) = δ   t 6 t 4 t 5 e r p r e R t d t   = δ e r p r   t 6 R [ e R t 4   e R t 5   ]
  • Per cycle carbon emissions cost due to scrapping
    ( C E s   ) = δ   t 6 e s p r E x r t 5 t 4   e R t 5
  • Per cycle carbon emissions cost due to holding inventory
    C E h = δ T e h [ t 2 t 3 I 3 t e R t d t + t 3 t 4 I 4 t e R t d t + t 4 t 5 I 5 t e R t d t + t 5 t 6 I 6 t e R t d t + t 1 t 3 I 7 t e R t d t + t 3 t 4 I 8 t e R t d t + t 4 t 5 I 9 t e R t d t ] = δ   t 6 e h [ { p { 1 E [ x p ] } d } θ [ 1 R ( e R t 2 e R t 3 ) e θ t 2 θ + R ( e θ + R t 2 e θ + R t 3 ) ] + d θ [ e θ t 4 θ + R ( e θ + R t 3 e θ + R t 4 ) + 1 R ( e R t 4 e R t 3 ) ] + p r 1 E x r d θ [ 1 R ( e R t 4 e R t 5 ) e θ t 4 θ + R ( e θ + R t 4 e θ + R t 5 ) ] + d θ [ e θ t 6 θ + R ( e θ + R t 5 e θ + R t 6 ) + 1 R ( e R t 6 e R t 5 ) + p E [ x p ]   θ [ 1 R ( e R t 1 e R t 3 e θ t 1 θ + R ( e θ + R t 1 e θ + R t 3 ) ] + 1 θ [ e θ t 3 e θ t 1   p E [ x p 1 θ + R ( e θ + R t 3 e θ + R t 4 ) + p r θ [ e θ t 5 θ + R ( e θ + R t 4 e θ + R t 5 ) ] 1 R ( e R t 4 e R t 5 ) ] .
Carbon emissions cost per cycle due to deterioration
C E θ = δ   t 6   e θ   [ t 2 t 3 p { 1 E [ x p ] } e R t d t   t 2 t 4 d e R t d t + t 4 t 5 p r 1 E x r e R t d t   t 4 t 6 d e R t d t + t 1 t 3 θ I 7 t e R t d t + t 3 t 4 θ I 8 t e R t d t + t 3 t 4 θ I 9 t e R t d t ] = δ   t 6   e θ   [ p { 1 E [ x p ] } 1 R e R t 2 e R t 3   d 1 R e R t 2 e R t 4 + p r 1 E x r   1 R [ e R t 4 e R t 5 ]   d 1 R e R t 4 e R t 6 + p E [ x p ] [   1 R e R t 1 e R t 3 + e θ t 1 θ + R ( e θ + R t 3 e θ + R t 1 ) + p E [ x p ]   ( e θ t 3 e θ t 1 )   1 θ + R ( e θ + R t 3 e θ + R t 4 ) + p r [ e θ t 5 θ + R ( e θ + R t 4 e θ + R t 5 ) + 1 R ( e R t 5 e R t 4 ) ] ] .
Total carbon emission cost
T C 2   = C S E   + C E p + C E r + C E s   + C E h   + C E θ  
Total   cost   T C = T C 1   + T C 2   .
Now, it is considered that the rate of random defectiveness during production and rework follows two kinds of distribution: (1) uniform distribution and (2) triangular distribution.
Case 1. Rates of random defectiveness follow a uniform distribution.
It is assumed that the rates of random defectiveness x p , x r both follow a uniform distribution. Then, the expected values of x p , x r are given by
E x p = l 1 + m 1 2   ,   l 1   > 0 ,   m 1 > 0   ,   l 1   <   m 1
E x r = l 2 + m 2 2   ,   l 2   > 0 ,   m 2 > 0   ,   l 2   <   m 2 .
After putting the values of E x p , E x r in Equation (25), we can obtain the total cost of the system.
Objective
Min   T C = T C ( p ,   t 1 ,   t 3 )
Subjected   to   p > 0 ,   t 1 > 0 ,   t 3 > 0 .
Case 2. Rates of random defectiveness follow the triangular distribution.
It is assumed that the rates of random defectiveness x p , x r both follow the triangular distribution. The expected values of x p , x r are given by
E x p = l 1 + m 1 + n 1 3 , l 1   <   m 1   and   l 1     n 1   m 1
E x r = l 2 + m 2 + n 2 3 ,   l 2   <   m 2 ,   and   l 2     n 2   m 2 .
Now, after putting the values of E x p , E x r in Equation (25), we can obtain the total cost.
Objective
Min   T C = T C ( p ,   t 1 ,   t 3 )
Subjected   to   p > 0 ,   t 1 > 0 ,   t 3 > 0 .
Theoretical Result: To derive the theoretical results, we apply some of the theorems of Cambini and Martein. According to them, any function can be written as
f ( x ) = F x G x     x   ϵ   R n ,
f(x) is a (strictly) pseudo-convex function, if F(x) is convex and differentiable and G(x) is positive and affine.
Theorem 1.
Cost function TC ( p ,   t 1 ,   t 3 ) is a pseudo-convex function with respect to p ,   t 1 , and t 3 , and hence there exist unique p ,   t 1 , and   t 3   that satisfy the following equations:
T C p = 0 ,   T C t 1 = 0 ,   T C t 3 = 0
such that TC( p ,   t 1 ,   t 3 ) is minimum.
Proof: 
See Appendix A. □

5. Solution Methodology

The objective functions given by Equations (26) and (27) are nonlinear functions of p ,   t 1 ,   and   t 3 . Thus, the closed-form solution is not an easy task. Thus, in order to obtain the optimal solution, the following solution methodology has been adopted.
  • Step 1: The values of all parameters are put in the objective functions in Equations (26) and (27).
  • Step 2: Different derivatives of objective functions p , t 1 , and t 3 are evaluated, and then each expression is put to zero to obtain the stationary points:
    T C p = 0
    T C t 1 = 0
    T C t 3 = 0 .
    On solving the above system of equations, we obtain the values p ,   t 1 ,   t 3 .
  • Step 3: The Hessian matrix is obtained as follows to check the convexity of the objective function.
    H = 2 T C p 2   2 T C p t 1   2 T C p t 3   2 T C t 1 p   2 T C t 1 2   2 T C t 1 t 3   2 T C t 3 p   2 T C t 3 t 1   2 T C t 3 2   .
After that, we calculate the different principal minors of the Hessian matrix as follows:
D 11 = 2 T C p 2 p ,   t 1 ,   t 3 ,   D 22 = 2 T C p 2 2 T C p t 1 2 T C t 1 p 2 T C t 1 2 p ,   t 1 ,   t 3 ,   D 33 = 2 T C p 2   2 T C p t 1   2 T C p t 3   2 T C t 1 p   2 T C t 1 2   2 T C t 1 t 3   2 T C t 3 p   2 T C t 3 t 1   2 T C t 3 2   p ,   t 1 ,   t 3 .
If D11 > 0; D22 > 0; and D33 > 0, i.e., the Hessian matrix is positive-definite, then at the obtained point p ,   t 1 ,   t 3 , the objective function is minimum.
  • Step 4: Thus, p * ,   t 1 * ,   t 3 * is an extreme point and T C p * ,   t 1 * ,   t 3 * is the optimal value of the objective function.

6. Numerical Analysis

The suggested model is illustrated with numerical examples in this section. For this purpose, data have been taken from the literature with appropriate modifications.
Example 1.
When x p , x r follow a uniform distribution.
d = 1200   unit / unit   time ,   p r = 0.05   unit / unit   time ,   c θ = 10   USD / unit ,   c r = 4   USD / unit ;   c b = 32   USD / unit / unit   time ;   c L = 27   USD / unit ,   c s = 3   USD / unit ,   c = 6 ,   s p = 250   USD / setup ,   s = 100   USD / setup ,   h s   = 5   USD / unit / unit   time ,   h r   = 6   USD / unit / unit   time ,   = 0.08 ,   e p = 5   kg / unit ,   e = 2   kg / unit ,   e sp = 100   kg / setup ,   e sr =   50   kg / setup ,   e h = 1   kg / unit ,   e = 10   kg / unit ,   e θ = 3   kg / unit ,   δ = 5   USD / kg ,   l 1 = 0.09 ,   m 1 =   0.11 ,   l 2 = 0.07 ,   m 2 = 0.08 ,   C 0 = 14   USD / unit ,   C 1 = 0.005   USD / unit ;   C 2 = 0.5   USD / unit ;   ε =   0.3 .
Now, on applying the solution methodology described in Section 4, the solution of Equations (28)–(30) is as follows:
p = 2050   unit / unit   time ,   t 1 =   0.47   years ,   t 3 = 1.48   years
Now, we evaluate the different principal minors of the Hessian matrix as follows:
D11 = 0.0114417 > 0, D12 = 27.1434 > 0, D13 = 21,897.9 > 0,
The above result shows that the objective function is convex. Hence, the optimal solution is as follows:
p * = 2050   unit / unit   time , t 1 * =   0.47   years ;   t 3 * = 1.48   years ;   TC * = 85,103.2   USD / cycle .
Example 2.
When x p , x r follow the triangular distribution.
d = 1200   unit / unit   time ,   p r = 0.05   unit / unit   time ,   c θ = 10   USD / unit ,   c r = 4   USD / unit ;   c b = 32   USD / unit / unit   time ;   c L = 27   USD / unit ,   c s = 3   USD / unit ,   c = 6 ,   s p = 250   USD / setup ,   s = 100   USD / setup ,   h s   = 5   USD / unit / unit   time ,   h r   = 6   USD / unit / unit   time ,   = 0.08 ,   e p = 5   kg / unit ,   e = 2   kg / unit ,   e sp = 100   kg / setup ,   e sr =   50   kg / setup ,   e h = 1   kg / unit ,   e = 10   kg / unit ,   e θ = 3   kg / unit ,   δ = 5   USD / kg ,   l 1 = 0.11 ,   m 1 =   0.09 ,   l 2 = 0.080 ,   m 2 = 0.070 ,   C 0 = 14   USD / unit ,   C 1 = 0.005   USD / unit ;   C 2 = 0.5   USD / unit ;   ε =   0.3 .
Now, on applying the solution methodology described in Section 4, the solution of Equations (28)–(30) is as follows:
p = 2068   unit / unit   time , t 1 = years , t 3 years .
Now, we evaluate different principal minors of the Hessian matrix as follows:
D11 = 0.0114417 > 0, D12 = 27.1434 > 0, D13 = 21,897.9 > 0.
The above result shows that the objective function is convex. Hence, the optimal solution is as follows:
p * = 2068   unit / unit   time ,   t 1 * =   0.39   years ;   t 3 * = 1.43   years ;   TC * = 109855   USD / cycle .

7. Sensitivity Analysis

Sensitivity analysis with respect to important parameters provides the decision maker with fundamental information about the model’s behavior (Table 2). This section discusses sensitivity to key parameters and their implications for the model.
The effect of the holding cost for non-serviceable and serviceable items, cost of scrapping, backorder cost, demand, setup cost for both the reworking station and manufacturing station, carbon tax, carbon emissions by production and rework, carbon emission by deterioration, and carbon emission by holding have been observed. The following inferences have been obtained.
(i)
The total cost of the system changes only slightly due to the rise in the holding cost for serviceable products, but there is a significant rise due to a rise in the holding cost for reworkable items.
(ii)
A minor change in the total inventory cost in the system is detected due to increases in the scrapping cost, setup cost for the reworking station, and setup cost for serviceable products.
(iii)
A high inclination in the total inventory cost of the system is observed due to the increase in the customer base.
(iv)
A high change in the total cost of the system is observed due to the change in emission by production, while a slight change is observed due to the change in emission by rework.
(v)
It is observed that holding cost, reworkable cost, scrapping cost, backorder cost, lost-sale cost, and carbon emissions costs are all positively correlated with total cost.

8. Managerial Insights

The current study’s goal is to deliver some constructive and economically useful insights to industry executives. The following are some of the most relevant managerial implications of the current study.
  • To lessen the influence of human mistakes in the screening process, the current study advises that automated screening procedures be established.
  • According to the study, the cost of carbon emissions related to production has a significant impact on the system’s total cost. As a result, managers must consider measures to reduce carbon emissions from manufacturing from both an environmental and economic standpoint.
  • Managers will have an understanding of how rework benefits the environment and how it can be used to meet customer requests as a result of this research.
  • Inflation is extremely sensitive in developing economies such as India, Bangladesh, and others. As a result, managers must consider this element when designing the optimal strategies for the system.

9. Concluding Remarks and Future Extension

Today’s competitive business environment had compelled sectors to make fact-based decisions to survive. Furthermore, organizations suffered financial losses as a result of carbon emissions and human errors. Furthermore, environmental challenges were linked to waste generation and carbon emissions. Taking cognizance of all these, a smart production-inventory model was developed for deteriorating items with an automated screening process. From the environmental point of view, costs due to carbon emissions, reworking, and waste management were considered in the model. Here, randomness was considered in proportion to defectiveness. The developed models were illustrated with numerical examples in both cases of uniform and triangular distributions, demonstrating which was better from an economic standpoint. Total inventory costs were determined to be 29% greater in the case of triangular distribution versus uniform distribution. The cost of carbon emissions had a significant negative influence on overall inventory cost, with an approximately 7% change in inventory cost seen as a result of this cost change. The findings suggested that managers should look for ways to expand their customer base, as the change in demand caused a 21% increase in inventory expenses. According to one study, carbon emissions have a negative impact on both the environment and the economy. Carbon emissions in the manufacturing process raise inventory costs by about 6%. As a result, managers must investigate the possibility of greenness in the manufacturing process. In this modern environment, we prefer machine systems over traditional human systems, and when it comes to reworking systems, we do not require a secondary market to sell damaged products at a discount. As a result, from an industrial standpoint, this study is relevant when considering the reworking process, backlog, and automation policy in screening. In the future, this model can be expanded to include selling price and advertisement-dependent demand to build a customer base, as well as time-varying deterioration, rework [31], and a carbon cap and trade policy for environmental sustainability. Preservation, green technologies, and outsourcing [32] can be utilized to expand this research and limit the number of imperfect objects. This study will aid in the reduction of carbon emissions and economic development analysis [33] by applying green technology within a complex production system [34].

Author Contributions

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

Funding

This research received no external funding.

Institutional Review Board Statement

Not Applicable.

Informed Consent Statement

Not Applicable.

Data Availability Statement

Provided within the paper.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Equation (25) can be written as
T C = 1 t 6 T C 1 + T C 2 .
Let F   =   T C 1 + T C 2 and G = t 6 .
Thus, F can be written as follows:
F = c p p X 1 + c r p r X 2 + ( s p + s r ) c b [ d ε e ε t 1 { 1 R ( 1 e R t 1 ) + 1   ε R ( 1 e ε R t 1 ) } + a d { 1 R e R t 1 t 1 t 2 + 1 R   X 3 } ] + c L d [ 1 R ( 1 e R t 1 e ε t 1 ε   R ( e t 1 ε R 1 ) ] + c p a X 4 + c p d   X 5 + c r b   X 2 c r d   X 6 + c p p E [ x p ]   X 1 c p b X 2 + h s [ a d θ { X 4 Y 2 X 7 } + d θ { Y 4 X 8 + X 9   } + b d θ { X 2 Y 4 X 10 } + d θ   { Y 6 X 11 + X 12   } + h r [   p E [ x p ]   θ { X 1 Y 1 X 13 } + p E [ x p ]   θ X 14 X 8 + p r θ { Y 5 X 10 X 2   } + δ   e s p + e s r + δ e p p X 1 + δ e r p r X 2 + δ e s p r   E x r X 15 + δ e h   [ a d θ { X 4 Y 2 X 7 } + d θ { X 8   Y 4 + X 9 } + b d θ   { X 2 Y 4 X 10 } + d θ   { Y 6 X 11 + X 12   } + p E [ x p ]   θ { X 1 Y 1 X 13 } + p E [ x p ]   θ X 14 X 8 + p r θ { Y 5 X 10 X 2   } ] + δ e θ [ a X 4 d X 5 + b X 2 d X 6 + p   E [ x p ] {   X 1 Y 1 X 13   } ]
where
X 1 = 1 R e R t 1 e R t 3 ;   X 2 = 1 R e R t 4 e R t 5 ;   X 3 = 1 R e R t 1 e R t 2 ;
X 4 = 1 R e R t 2 e R t 3 ;   X 5 = 1 R e R t 2 e R t 4 ;   X 6 = 1 R e R t 4 e R t 6 ] ;
X 7 = 1 θ + R [ e θ + R t 2 e θ + R t 3 ] ;   X 8 = 1 θ + R [ e θ + R t 3 e θ + R t 4 ] ; X 9 = 1 R [ e R t 4 e R t 3 ] ;
X 10 = 1 θ + R [ e θ + R t 4 e θ + R t 5 ] ;   X 11 = 1 θ + R [ ( e θ + R t 5 e θ + R t 6 ] ;
X 12 = 1 R [ e R t 6 e R t 5   ] ;   X 13 = 1 θ + R [ e θ + R t 2 e θ + R t 3 ] ;   X 14 = [ e θ t 3 e θ t 1 ] ;
X 15 = ( t 5 t 4 ) e R t 5 ;   Y 1 = e θ t 1 ; Y 2 = e θ t 2 ; Y 3 = e θ t 3 ;   Y 4 = e θ t 4 ;   Y 5 = e θ t 5 ;   Y 6 = e θ t 6
Now, we take the different derivatives with regard to different decision variables.
2 F p 2 = 2 C 1 X 3 + 2 C 1 ( 1 p E [ x p ] )   X 4 + 2 C 2 / p 3 d   X 5 + 2 C 1 E [ x p ]   X 1 2 C 2 / p 3   p r 1 E [ x r ) X 2 , 2 F p t 1 = 2 F t 1 p = ( C 0 + 2 C 1 p )   X 1 t 1 c b   ( 1 E [ x p ] ) {   1 R e R t 1 t 1 e R t 1 + t 2   e R t 1   + 1 R X 3 t 1 } + ( C 0 + 2 C 1 p )   ( 1 E [ x p ] )   X 4 t 1   ( C 1 C 2 / p 2 )   d X 5 t 1 + ( C 0 + 2 C 1 p )   E [ x p ]   X 1 t 1   ( C 1 C 2 / p 2 )   p r   1 E x r X 2 t 1 + h s [ 1 E [ x p ]   θ X 4 t 1 Y 2 X 7 t 1 Y 2 t 1 X 7 + h r   [ E [ x p ]   θ X 1 t 1 Y 1 X 13 t 1 Y 1 t 1 X 13   + E [ x p ]   θ { X 14 X 8 t 1 + X 14 t 1 X 8 } ] , 2 F p t 3 = 2 F t 3 p = ( C 0 + 2 C 1 p )   ( 1 E [ x p ] )   X 4 t 3     ( C 1 C 2 / p 2 )   d X 5 t 3 + ( C 0 + 2 C 1 p )   E [ x p ]   X 1 t 3   ( C 1 C 2 / p 2 )     p r   1 E x r X 2 t 3 + h s [ 1 E [ x p ]   θ X 4 t 3 Y 2 X 7 t 3 + h r   [ E [ x p ]   θ X 1 t 1 Y 1 X 13 t 3   + E [ x p ]   θ { X 14 X 8 t 3 + X 14 t 3 X 8 } ] , 2 F t 1 2 = c p p X 1 t 1 t 1 + c r p r X 2 t 1 t 1 c b [ d ε { 1 R   ε 2 e ε t 1 + ε + R e ε + R t 1 + 1   ε R (   ε 2 e ε t 1 R 2 e R t 1 ) } + a d   { 1 R R 2 t 1 e R t 1 + R 2 t 2 e R t 1 2 R e R t 1 } + 1 R X 3 t 1 t 1 ] + c L d [ R e R t 1 1   ε R ( R e R t 1 ε 2 e ε t 1 ]   + c p a X 4 t 1 t 1 c p d X 5 t 1 t 1 + c r b X 2 t 1 t 1 + h s [ a d θ { X 4 t 1 t 1 Y 2 X 7 t 1 t 1 2 Y 2 t 1 X 7 t 1 Y 2 t 1 t 1 X 7 } + d θ { Y 4 X 8 t 1 t 1 + 2 Y 4 t 1 X 8 t 1 + Y 4 t 1 t 1 X 8 + X 9 t 1 t 1 } + b d θ { X 2 t 1 t 1 Y 4 X 10 t 1 t 1 2 Y 4 t 1 X 10 t 1 Y 4 t 1 t 1 X 10 } + d θ { Y 6 X 11 t 1 t 1 + 2 Y 6 t 1 X 11 t 1 + Y 6 t 1 t 1 X 11 + X 12 t 1 t 1 } ] + h r [ p E [ x p ]   θ   { X 1 t 1 t 1 Y 1 X 13 t 1 t 1 2 Y 1 t 1 X 13 t 1 Y 1 t 1 t 1 X 13 } + p E [ x p ]   θ   { X 14 t 1 t 1 X 8 + 2 X 8 t 1 X 14 t 1 + X 8 t 1 t 1 X 15 } + p r θ   { Y 5 X 10 t 1 t 1 + 2 Y 1 t 1 X 10 t 1 + Y 5 t 1 t 1 X 10 X 2 t 1 t 1 } + δ e p p X 1 t 1 t 1 + δ e r p r X 2 t 1 t 1 + δ e s p r   E x r   X 15 t 1 t 1 + δ e h [ a d θ   { X 4 t 1 t 1 Y 2 X 7 t 1 t 1 2 Y 2 t 1 X 7 t 1 Y 2 t 1 t 1 X 7 } + + d θ { Y 4 X 8 t 1 t 1 + 2 Y 4 t 1 X 8 t 1 + Y 4 t 1 t 1 X 8 + X 9 t 1 t 1 } + b d θ X 2 t 1 t 1 Y 4 X 10 t 1 t 1 2 Y 4 t 1 X 10 t 1 Y 4 t 1 t 1 X 10 + d θ { Y 6 X 11 t 1 t 1 + 2 Y 6 t 1 X 11 t 1 + Y 6 t 1 t 1 X 11 + X 12 t 1 t 1 } + p E [ x p ]   θ { X 1 t 1 t 1 Y 1 X 13 t 1 t 1 2 Y 1 t 1 X 13 t 1 Y 1 t 1 t 1 X 13 } + p E [ x p ]   θ   { X 14 t 1 t 1 X 8 + 2 X 8 t 1 X 14 t 1 + X 8 t 1 t 1 X 15 } + p r θ   { Y 5 X 10 t 1 t 1 + 2 Y 1 t 1 X 10 t 1 + Y 5 t 1 t 1 X 10 X 2 t 1 t 1 } ] + δ e θ   [ a X 4 t 1 t 1 d X 5 t 1 t 1 + b X 2 t 1 t 1 d X 6 t 1 t 1 ] + p E [ x p ] X 1 t 1 t 1 Y 1 X 13 t 1 t 1 2 Y 1 t 1 X 13 t 1 Y 1 t 1 t 1 X 13
2 F t 1 t 3 = 2 F t 3 t 1 = c p p X 1 t 1 t 3 + c r p r X 2 t 1 t 3 + c p a X 4 t 1 t 3 + c p d X 5 t 1 t 3 + c r b X 2 t 1 t 3 + h s [ a d θ { X 4 t 1 t 3 Y 2 X 7 t 1 t 3 Y 2 t 1 X 7 t 3 Y 2 t 3 X 7 t 1 Y 2 t 1 t 3 X 7 } + d θ   Y 4 X 8 t 1 t 3 + Y 4 t 1 X 8 t 3 + Y 4 t 3 X 8 t 1 + Y 4 t 1 t 3 X 8 + X 9 t 1 t 3 + b d θ { X 2 t 1 t 3 Y 4 X 10 t 1 t 3   Y 4 t 1 X 10 t 3 Y 4 t 3 X 10 t 1 Y 4 t 1 t 3 X 10 } + d θ { Y 6 X 11 t 1 t 3 + Y 6 t 1 X 11 t 3 + Y 6 t 3 X 11 t 1 + Y 6 t 1 t 3 X 11 + X 12 t 1 t 3 } ] +   h r [ p E [ x p ]   θ   { X 1 t 1 t 3 Y 1 X 13 t 1 t 3 Y 1 t 1 X 13 t 3 Y 1 t 3 X 13 t 1 Y 1 t 1 t 3 X 13 } + p E [ x p ]   θ   { X 14 t 1 t 3 X 8 + X 8 t 1 X 14 t 3 + X 8 t 3 X 14 t 1 + X 8 t 1 t 3 X 15 } + p r θ   { Y 5 X 10 t 1 t 3 + Y 1 t 1 X 10 t 3 + Y 1 t 3 X 10 t 1 + Y 5 t 1 t 3 X 10 X 2 t 1 t 3 } + δ e p p X 1 t 1 t 3 + δ e r p r X 2 t 1 t 3   +   δ e s p r   E x r   X 15 t 1 t 3 + δ e h [ a d θ   { X 4 t 1 t 3 Y 2 X 7 t 1 t 3 Y 2 t 1 X 7 t 3 Y 2 t 3 X 7 t 1 Y 2 t 1 t 3 X 7 }   + d θ { Y 4 X 8 t 1 t 3 + Y 4 t 1 X 8 t 3 + Y 4 t 3 X 8 t 1 + Y 4 t 1 t 3 X 8 + X 9 t 1 t 3 } + b d θ X 2 t 1 t 3 Y 4 X 10 t 1 t 3 Y 4 t 1 X 10 t 3 Y 4 t 3 X 10 t 1 Y 4 t 1 t 3 X 10 + d θ { Y 6 X 11 t 1 t 3 + Y 6 t 1 X 11 t 3 + Y 6 t 3 X 11 t 1 + Y 6 t 1 t 3 X 11 + X 12 t 1 t 3 } + p E [ x p ]   θ   { X 1 t 1 t 3 Y 1 X 13 t 1 t 3 Y 1 t 1 X 13 t 3 Y 1 t 3 X 13 t 1 Y 1 t 1 t 3 X 13 }   + p E [ x p ]   θ   { X 14 t 1 t 3 X 8 + X 8 t 1 X 14 t 3 + X 8 t 3 X 14 t 1 + X 8 t 1 t 3 X 15 }   + p r θ   { Y 5 X 10 t 1 t 1 + 2 Y 1 t 1 X 10 t 1 + Y 5 t 1 t 1 X 10 X 2 t 1 t 3 } ] + δ e θ   [ a X 4 t 1 t 3 d X 5 t 1 t 3 + b X 2 t 1 t 3 d X 6 t 1 t 3 ]   + p E [ x p ]   X 1 t 1 t 3 Y 1 X 13 t 1 t 3 Y 1 t 1 X 13 t 3 Y 1 t 3 X 13 t 1 Y 1 t 1 t 3 X 13 2 F t 3 2   = c p a X 4 t 3 t 3 c p d X 5 t 3 t 3 + c r b X 2 t 3 t 3 + h s [ a d θ { X 4 t 3 t 3 Y 2 X 7 t 3 t 3 } + d θ { Y 4 X 8 t 3 t 3 4 + 2 Y 4 t 3 X 8 t 3 + Y 4 t 3 t 3 X 8 } + b d θ { X 2 t 3 t 3 Y 4 X 10 t 3 t 3 2 Y 4 t 3 X 10 t 3 Y 4 t 3 t 3 X 10 } + d θ { Y 6 X 11 t 3 t 3 + 2 Y 6 t 3 X 11 t 3 + Y 6 t 3 t 3 X 11 + X 12 t 3 t 3 } ] + h r [ p E [ x p ]   θ { X 1 t 3 t 3 Y 1 X 13 t 3 t 3 } + p E [ x p ]   θ   { X 14 t 3 t 3 X 8 + 2 X 8 t 3 X 14 t 3 + X 8 t 3 t 3 X 15 } + p r θ   { Y 5 X 10 t 3 t 3 + 2 Y 1 t 3 X 10 t 3 + Y 5 t 3 t 3 X 10 X 2 t 3 t 3 } + δ e p p X 1 t 3 t 3 + δ e r p r X 2 t 3 t 3 + δ e s p r   E x r   X 15 t 3 t 3 + δ e h [ a d θ   { X 4 t 3 t 3 Y 2 X 7 t 3 t 3 } + d θ { Y 4 X 8 t 3 t 3 + 2 Y 4 t 3 X 8 t 3 + Y 4 t 3 t 3 X 8 + X 9 t 3 3 } + b d θ X 2 t 3 t 3 Y 4 X 10 t 3 t 3 2 Y 4 t 3 X 10 t 3 Y 4 t 3 t 3 X 10 + d θ { Y 6 X 11 t 3 t 3 + 2 Y 6 t 3 X 11 t 3 + Y 6 t 3 t 3 X 11 + X 12 t 3 t 3 } + p E [ x p ]   θ { X 1 t 3 t 3 Y 1 X 13 t 3 t 3 } + p E [ x p ]   θ   { X 14 t 3 t 3 X 8 + 2 X 8 t 3 X 14 t 3 + X 8 t 3 t 3 X 15 } + p r θ   { Y 5 X 10 t 3 t 3 + 2 Y 5 t 3 X 10 t 3 + Y 5 t 3 t 3 X 10 X 2 t 3 t 3 } ] + δ e θ   [ a X 4 t 3 t 3 d X 5 t 3 t 3 + b X 2 t 3 t 3 d X 6 t 3 t 3 ] + p E [ x p ]    X 1 t 3 t 3 Y 1 X 13 t 3 t 3
where X 1 t 1 = X 1 t 1 ;   X 2 t 1 = X 2 t 1 ;   X 3 t 1 = X 3 t 1 ;   X 4 t 1 = X 4 t 1 ;   X 5 t 1 = X 5 t 1 ;   X 7 t 1 = X 7 t 1 ;   X 8 t 1 = X 8 t 1 ; X 13 t 1 = X 13 t 1 ; X 14 t 1 = X 14 t 1 ;   X 1 t 3 = X 1 t 3 ;   X 2 t 3 = X 2 t 3 ;   X 3 t 3 = X 3 t 3 ;   X 4 t 3 = X 4 t 3 ;   X 5 t 3 = X 5 t 3 ;   X 7 t 3 = X 7 t 3 ;   X 8 t 3 = X 8 t 3 ;   X 13 t 3 = X 13 t 3 ; X 14 t 3 = X 14 t 3 ; X 1 t 1 t 1 = 2 X 1 t 1 2 ;   X 2 t 1 t 1 = 2 X 2 t 1 2 ;   X 3 t 1 t 1 = 2 X 3 t 1 2 ; X 4 t 1 t 1 = 2 X 4 t 1 2 ;   X 5 t 1 t 1 = 2 X 5 t 1 2 ;   X 6 t 1 t 1 = 2 X 6 t 1 2 ;   X 7 t 1 t 1 = 2 X 7 t 1 2 ;   X 8 t 1 t 1 = 2 X 8 t 1 2 ; X 9 t 1 t 1 = 2 X 9 t 1 2 ;   X 10 t 1 t 1 = 2 X 10 t 1 2 ;   X 11 t 1 t 1 = 2 X 11 t 1 2 ;   X 12 t 1 t 1 = 2 X 12 t 1 2 ;   X 13 t 1 t 1 = 2 X 13 t 1 2 ; X 14 t 1 t 1 = 2 X 14 t 1 2 ;   X 15 t 1 t 1 = 2 X 15 t 1 2 ;   X 1 t 3 t 3 = 2 X 1 t 3 2 ;   X 2 t 3 t 3 = 2 X 2 t 3 2 ;   X 3 t 3 t 3 = 2 X 3 t 3 2 ; X 4 t 3 t 3 = 2 X 4 t 3 2 ;   X 5 t 3 t 3 = 2 X 5 t 3 2 ;   X 6 t 3 t 3 = 2 X 6 t 3 2 ;   X 7 t 3 t 3 = 2 X 7 t 3 2 ;   X 8 t 3 t 3 = 2 X 8 t 3 2 ; X 9 t 3 t 3 = 2 X 9 t 3 2 ;   X 10 t 3 t 3 = 2 X 10 t 3 2 ;   X 11 t 3 t 3 = 2 X 11 t 3 2 ;   X 12 t 3 t 3 = 2 X 12 t 3 2 ;   X 13 t 3 t 3 = 2 X 13 t # 2 ; X 14 t 3 t 3 = 2 X 14 t 3 2 ;   X 15 t 3 t 3 = 2 X 15 t 3 2 ;   X 1 t 1 t 3 = 2 X 1 t 1 t 3 ;   X 2 t 1 t 3 = 2 X 2 t 1 t 3 ;   X 3 t 1 t 3 = 2 X 3 t 1 t 3 ;   X 4 t 1 t 3 = 2 X 4 t 1 t 3 ;   X 5 t 1 t 3 = 2 X 5 t 1 t 3 ;   X 6 t 1 t 3 = 2 X 6 t 1 t 3 ;   X 7 t 1 t 3 = 2 X 7 t 1 t 3 ;   X 8 t 1 t 3 = 2 X 8 t 1 t 3 ;   X 9 t 1 t 3 = 2 X 9 t 1 t 3 ;   X 10 t 1 t 3 = 2 X 10 t 1 t 3 ;   X 11 t 1 t 3 = 2 X 11 t 1 t 3 ;   X 12 t 1 t 3 = 2 X 12 t 1 t 3 ;   X 13 t 1 t 3 = 2 X 13 t 1 t 3 ; X 14 t 1 t 3 = 2 X 14 t 1 t 3 ;   X 15 t 1 t 3 = 2 X 15 t 1 t 3 ;   Y 1 t 1 = Y 1 t 1 ;   Y 2 t 1 = Y 2 t 1 ;   Y 3 t 1 = Y 3 t 1 ;   Y 4 t 1 = Y 4 t 1 ; Y 5 t 1 = Y 5 t 1 ;   Y 6 t 1 = Y 6 t 1 ;   Y 1 t 3 = Y 1 t 3 ;   Y 2 t 3 = Y 2 t 3 ;   Y 3 t 3 = Y 3 t 3 ;   Y 4 t 3 = Y 4 t 3 ;   Y 5 t 3 = Y 5 t 3 ;   Y 6 t 3 = Y 6 t 3 ;   Y 1 t 1 t 3 = 2 Y 1 t 1 t 3 ;   Y 2 t 1 t 3 = 2 Y 2 t 1 t 3 ;   Y 3 t 1 t 3 = 2 Y 3 t 1 t 3 ;   Y 4 t 1 t 3 = 2 Y 4 t 1 t 3 ;   Y 5 t 1 t 3 = 2 Y 5 t 1 t 3 ;   Y 6 t 1 t 3 = 2 Y 6 t 1 t 3 .
The Hessian matrix for F is as follows:
2 F p 2   2 F p t 1   2 F p t 3   2 F t 1 p   2 F t 1 2   2 F t 1 t 3   2 F t 3 p   2 F t 3 t 1   2 F t 3 2   .
Now, the first principal minor is
H 11 = 2 T C p 2 > 0 .
The second principal minor is
H 22 = 2 T C p 2 2 T C p t 1 2 T C t 1 p 2 T C t 1 2 > 0
and the third principal minor is
H 33 = 2 T C p 2   2 T C p t 1   2 T C p t 3   2 T C t 1 p   2 T C t 1 2   2 T C t 1 t 3   2 T C t 3 p   2 T C t 3 t 1   2 T C t 3 2   > 0 .
Thus, the Hessian matrix is positive-definite. Therefore, F is a differentiable as well as a convex function with regard to decision variables. Further, G = t 6 > 0 and a differentiable and affine function. Therefore, F attains the global minimum value at p * , t 1 * , t 3 * , which satisfies the equation T C p = 0 ; T C t 1 = 0 ; T C t 2 = 0 .

References

  1. Hua, G.; Cheng, T.; Wang, S. Managing carbon footprints in inventory management. Int. J. Prod. Econ. 2011, 132, 178–185. [Google Scholar] [CrossRef]
  2. Benjaffar, S.; Li, Y.; Daskin, M. Carbon footprint and the management of supply chain: Insights from simple model. IEEE Trans. Autom. Sci. Eng. 2012, 10, 99–116. [Google Scholar] [CrossRef]
  3. Lou, G.X.; Yang, H.; Zhang, J.Q.; Fan, T.J. Investment strategy of emission reduction technology in a supply chain. Sustainability 2015, 7, 10684–10708. [Google Scholar] [CrossRef] [Green Version]
  4. Datta, T.K. Effect of green technology investment on a production-inventory system with carbon tax. Adv. Oper. Res. 2017, 2017, 4834839. [Google Scholar] [CrossRef] [Green Version]
  5. 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]
  6. Jauhari, W.A.; Adam, N.A.F.P.; Rosyidi, C.N.; Pujawan, I.N.; Shah, N.H. A closed-loop supply chain model with rework waste disposal, and carbon emission. Oper. Res. Perspect. 2020, 7, 100155. [Google Scholar] [CrossRef]
  7. Mishra, M.; Ghosh, S.K.; Sarkar, B. Maintaining energy efficiencies and reducing carbon emissions under a sustainable supply chain management. AIMS Environ. Sci. 2022, 9, 603–635. [Google Scholar] [CrossRef]
  8. Yadav, D.; Singh, R.; Kumar, A.; Sarkar, B. Reduction of pollution through sustainable and flexible production by controlling by-products. J. Environ. Inform. 2022, 40, 106–124. [Google Scholar] [CrossRef]
  9. Sarkar, B.; Kar, S.; Basu, K.; Guchhait, R. A sustainable managerial decision-making problem for a substitutable product in a dual-channel under carbon tax policy. Comput. Ind. Eng. 2022, 172, 108635. [Google Scholar] [CrossRef]
  10. Sarkar, M.; Sarkar, B.; Dolgui, A. An automated smart production with system reliability under a leader-follower strategy of supply chain management. In Advances in Production Management Systems. Smart Manufacturing and Logistics Systems: Turning Ideas into Action. APMS 2022. IFIP Advances in Information and Communication Technology; Kim, D.Y., von Cieminski, G., Romero, D., Eds.; Springer: Cham, Switzerland, 2022; Volume 663, pp. 459–467. [Google Scholar] [CrossRef]
  11. Wee, H.M.; Huang, Y.D.; Wang, W.T.; Cheng, Y.L. An EPQ model with partial backorders considering two backordering costs. Appl. Math. Comput. 2014, 232, 898–907. [Google Scholar] [CrossRef]
  12. Mishra, U.; Cárdenas Barrón, L.E.; Tiwari, S.; Shaikh, A.A.; Trevino-Gazra, G. An inventory model under price and stock dependent demand for controllable deterioration rate with shortages and preservation technology investment. Ann. Oper. Res. 2017, 254, 165–190. [Google Scholar] [CrossRef]
  13. Singh, S.R.; Rani, M. An inventory model for exponential time—Sensitive demand and parabolic time linked holding cost under inflation and shortages with salvage value. PalArchs J. Archaeol. Egypt Egyptol. 2020, 17, 10941–10949. [Google Scholar]
  14. Palanivel, M.; Suganya, M. Partial backlogging inventory model with price and stock level dependent demand, time varying holding cost and quantity discounts. J. Manag. Anal. 2021, 9, 32–59. [Google Scholar] [CrossRef]
  15. Kumar, N.; Dahiya, S.; Kumar, S. Two warehouse inventory model for deteriorating items with fixed shelf-life stock-dependent demand and partial backlogging. J. Math. Comput. Sci. 2022, 12, 85. [Google Scholar] [CrossRef]
  16. Sarkar, B.; Takeyeva, D.; Guchhait, R.; Sarkar, M. Optimized radio-frequency identification system for different warehouse shapes. Knowl. Based Syst. 2022, 258, 109811. [Google Scholar] [CrossRef]
  17. Jawla, P.; Singh, S. Multi-item economic production quantity model for imperfect items with multiple production setups and rework under the effect of preservation technology and learning environment. Int. J. Ind. Eng. Comput. 2016, 7, 703–716. [Google Scholar] [CrossRef]
  18. Khara, B.; Dey, J.K.; Mondal, S.K. An inventory model under development cost-dependent imperfect production and reliability dependent demand. J. Manag. Anal. 2017, 4, 258–275. [Google Scholar] [CrossRef]
  19. Al-Salamah, M. Economic production quantity in an imperfect manufacturing process with synchronous and Asynchronous flexible rework rates. Oper. Res. Perspect. 2019, 6, 100103. [Google Scholar] [CrossRef]
  20. Sarkar, B.; Ganguly, B.; Pareek, S.; Cárdenas-Barrón, L.E. A three-echelon green supply chain management for biodegradable products with three transportation modes. Comp. Indust. Eng. 2022, 174, 108727. [Google Scholar] [CrossRef]
  21. Gupta, S.; Vijaygargy, L.; Sarkar, B. A bi-objective integrated transportation and inventory management under a supply chain network considering multiple distribution networks. RAIRO-Operat. Res. 2022, 56, 3991–4022. [Google Scholar] [CrossRef]
  22. Lin, H.J. An economic production quantity model with backlogging and imperfect rework process for uncertain demand. Int. J. Prod. Res. 2021, 59, 467–482. [Google Scholar] [CrossRef]
  23. Chandra, M.J.; Michael, L.B. The effect of inflation and the time value of money on some inventory system. Int. J. Prod. Res. 2007, 23, 723–730. [Google Scholar] [CrossRef]
  24. Yang, H.L.; Teng, J.T.; Chern, M.S. An inventory model under inflation for deteriorating items with stock-dependent consumption rate and partial backlogging shortages. Int. J. Prod. Econ. 2010, 123, 8–23. [Google Scholar] [CrossRef]
  25. Palanivel, M.; Uthayakumar, R. An EPQ model with variable production, Probabilistic deterioration and partial backlogging under inflation. J. Manag. Anal. 2014, 1, 200–223. [Google Scholar] [CrossRef]
  26. Rani, S.; Ali, R.; Agarwal, A. Green supply chain inventory model for deteriorating items with variable demand under inflation. Int. J. Bus. Forecast. Mark. Intell. 2017, 3, 50–77. [Google Scholar] [CrossRef]
  27. Huang, X.; Yang, S.; Wang, Z. Optimal pricing and replenishment policy for perishable food supply chain under inflation. Comput. Ind. Eng. 2021, 158, 107433. [Google Scholar] [CrossRef]
  28. Alamri, O.A.; Jayaswal, M.K.; Khan, F.A.; Mittal, M. An EOQ Model with Carbon Emissions and Inflation for Deteriorating Imperfect Quality Items under Learning Effect. Sustainability 2022, 14, 1365. [Google Scholar] [CrossRef]
  29. Padiyar, S.V.S.; Vandana, V.; Bhagat, N.; Singh, S.R.; Sarkar, B. Joint replenishment strategy for deteriorating multi-item through multi-echelon supply chain model with imperfect production under imprecise and inflationary environment. RAIRO-Oper. Res. 2022, 56, 3071–3096. [Google Scholar] [CrossRef]
  30. Singh, S.R.; Sharma, S. A partially backlogged supply chain model for deteriorating items under reverse logistics, imperfect production/ remanufacturing and inflation. Int. J. Logist. Syst. Manag. 2019, 33, 221–255. [Google Scholar] [CrossRef]
  31. Rout, C.; Chakraborty, D.; Goswami, A. An EPQ model for deteriorating items with imperfect production, two type of inspections errors and rework under complete backordering. Int. Game Theory Rev. 2020, 22, 2040011. [Google Scholar] [CrossRef]
  32. Bachar, R.K.; Bhuniya, S.; Ghosh, S.K.; Sarkar, B. Sustainable green production model considering variable demand, partial outsourcing, and rework. AIMS Environ. Sci. 2022, 9, 325–353. [Google Scholar] [CrossRef]
  33. Oryani, B.; Moridian, A.; Sarkar, B.; Rezania, S.; Kamyab, H.; Khan, M.K. Assessing the financial resourse curse hypothesis in Iran: The novel dynamic ARDL approach. Resour. Pol. 2022, 78, 102899. [Google Scholar] [CrossRef]
  34. Hota, S.K.; Ghosh, S.K.; Sarkar, B. Involvement of smart technologies in an advanced supply chain management to solve unreliability under distribution robust approach. AIMS Environ. Sci. 2022, 9, 461–492. [Google Scholar] [CrossRef]
Figure 1. Emissions of greenhouse gases in the US by economic sector in 2019 (www.epa.gov/ghgemissions/overview-greenhouse-gases, accessed on 7 July 2021).
Figure 1. Emissions of greenhouse gases in the US by economic sector in 2019 (www.epa.gov/ghgemissions/overview-greenhouse-gases, accessed on 7 July 2021).
Mathematics 11 00440 g001
Figure 2. The flow chart of the current study.
Figure 2. The flow chart of the current study.
Mathematics 11 00440 g002
Figure 3. Level of serviceable inventory in production and rework.
Figure 3. Level of serviceable inventory in production and rework.
Mathematics 11 00440 g003
Figure 4. Level of reworkable inventory in production and rework.
Figure 4. Level of reworkable inventory in production and rework.
Mathematics 11 00440 g004
Table 1. Summary of selected literature and the proposed model.
Table 1. Summary of selected literature and the proposed model.
Author(s)Model TypeCarbon EmissionVolume FlexibilityInflationBackorderImperfect Process
Hua et al. [1]EOQ----
Benjaafar et al. [2]SCM--Full-
Lou et al. [3]EOQ
Datta [4]EPQ---Other
Daryanto et al. [5]SCM---Constant
Jauhari et al. [6]SCM---Constant
Mishra at el. [7]SCM----
Yadav et al. [8]SCM-PartialConstant
Sarkar et al. [9]SCM----
Sarkar et al. [10]SCM----Random
Wee et al. [11]EPQ---Partial-
Mishra et al. [12]EOQ---Both-
Singh and Rani [13]EOQ--Partial-
Kumar [15]EOQ---Partial-
Sarkar et al. [16]SCM-----
Jawla and Singh [17]EPQ---Constant
Khara [18]EPQ----Constant
Al-Salamah [19]EPQ---FullConstant
Sarkar et al. [20]SCM----
Sarkar et al. [21]EOQ-----
Lin [22]EPQ---FullConstant
Chandra and Michael [23]EOQ--Full-
Yang et al. [24]EOQ--Partial-
Palanivel and Uthayakumar [25]EPQ--Partial-
Rani and Ali [26]SCM----
Huang et al. [27]EOQ----
Alamri et al. [28]EOQ--Random
Padiyar et al. [29]SCM---Constant
Singh and Sharma [30]SCM--PartiallyConstant
This paperEPQ
SCM: supply chain model; EPQ: economic production quantity; EOQ: economic order quantity.
Table 2. Sensitive analysis with respect to key parameters.
Table 2. Sensitive analysis with respect to key parameters.
Parameters% ChangeTC% Change in TC
Serviceable items holding cost ( h s )+2085,114.9+0.014
+1085,109+0.007
−1085,097.3−0.007
−2085,091.5−0.014
Reworkable items holding cost ( h r )+2090,048.8+5.81
+1087,576+2.91
−1082,630.4−2.91
−2080,157.6−5.81
Scrapping cost ( c s )+2085,108.4+0.006
+1085,105.8+0.003
−1085,100.6−0.003
−2085,098−0.006
Setup cost for
production (sp)
+2085,132+0.034
+1085,117.6+0.017
−1085,088.8−0.017
−2085,074.4−0.034
Setup cost for
reworking station (sr)
+2085,114.7+0.014
+1085,108.9+0.007
−1085,097.4−0.007
−2085,091.7−0.014
Demand rate ( d )+20103,356+21.45
+1093,765.1+10.18
−1076,909.6−9.63
−2068,973.7−18.95
Backorder cost ( c b )+2085,536.5+0.51
+1085,319.9+0.26
−1084,886−0.26
−2084,669.9−0.51
Carbon tax per cycle ( δ )+2091,084.2+7.03
+1088,093.7+3.51
−1082,112.7−3.51
−2079,122.2−7.03
Carbon emission by production ( e p )+2090,612.4+6.47
+1087,857.8+3.24
−1082,348.6−3.24
−2079,594−6.47
Carbon emission due to re-working ( e r )+2085,305.8+0.24
+1085,204.5+0.12
−1085,001.9−0.12
−2084,900.6−0.24
Carbon emission by holding the items in warehouse ( e h )+2085,185.1+0.096
+1085,144.1+0.048
−1085,062.3−0.048
−2085,021.3−0.096
Carbon emission by deterioration ( e θ )+2085,117.7+0.017
+1085,110.4+0.008
−1085,095.9−0.009
−2085,088.7−0.012
Carbon emission by scrapping ( e s )+2085,189.8+0.10
+1085,146.5+0.05
−1085,059.9−0.05
−2085,016.6−0.10
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yadav, D.; Chand, U.; Goel, R.; Sarkar, B. Smart Production System with Random Imperfect Process, Partial Backordering, and Deterioration in an Inflationary Environment. Mathematics 2023, 11, 440. https://doi.org/10.3390/math11020440

AMA Style

Yadav D, Chand U, Goel R, Sarkar B. Smart Production System with Random Imperfect Process, Partial Backordering, and Deterioration in an Inflationary Environment. Mathematics. 2023; 11(2):440. https://doi.org/10.3390/math11020440

Chicago/Turabian Style

Yadav, Dharmendra, Umesh Chand, Ruchi Goel, and Biswajit Sarkar. 2023. "Smart Production System with Random Imperfect Process, Partial Backordering, and Deterioration in an Inflationary Environment" Mathematics 11, no. 2: 440. https://doi.org/10.3390/math11020440

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