Abstract

The potential to obtain defective or damaged items with non-defective commodities is common to experience at the production unit or when shipping products from one layer to another. This research focuses on the faulty things that retailers receive from suppliers. The retailer has set a restriction on the percentage of defective things, and the retailer receives a discount on the cost of purchasing defective items. The proposed inventory system handles the uncertainty in inventory costs and also considers the demand and deterioration of items with prioritized maximum product life. This work minimizes total inventory cost when demand rate as a function of reliability and power pattern of time under a crisp and triangular neutrosophic environment. The inventory system for degrading items considers the predictability and power pattern of time with a reasonable payment delay. The interest charges are applied only after a specific permissible time limit in the proposed inventory system. The neutrosophic number that defines three different kinds of membership functions representing the truth, hesitation, and falseness is applied in the inventory model in handling the uncertainty of the cost pattern. The proposed inventory model is investigated using a particle swarm optimization algorithm, and the results are validated using a numerical example and a sensitivity analysis for various parameters.

1. Introduction

In dealing with inventory models, researchers need to give more importance to some key factors, such as deterioration, demand, and reliability. It is obvious that deterioration is a time-dependent factor, and furthermore, it increases with time, resulting in a decrease in the demand for the commodity. When products are stored in warehouses, there are plenty of problems, such as reliability, volume, and deterioration of items, which are different for every commodity. Consequently, the holding costs have a major impact on the value of the quantity stored. Deterioration is a progressive degradation of stored products that results in their devaluation. The rate of deterioration is shallow in durable products such as iron products, steels, toys, consumer electronics, furniture, tools, jewelry, automobiles, sports equipment, and so on. The deterioration rate encounters rapid changes in semi-durable products such as food items, drugs, clothing, cosmetics, and so on. Therefore, the study of the deterioration of items in the inventory systems plays a major role due to different deterioration patterns in the economic order quantity (EOQ) model. Tadikamalla [1] introduced an EOQ model using a gamma distribution to represent the constant, increasing, and decreasing rates of deterioration over time. Alshanbari et al. [2] proposed an inventory model for deteriorating items with two-parameter Weibull distribution. Wang and Lin [3] developed the optimal replenishment strategy with deterioration, market demand, and price changes. An inventory model for deteriorating items with maximum lifetime was proposed by Wang et al. [4]. Skouri et al. [5] proposed an inventory model with shortages, ramp-type demand rate, and time-dependent deterioration rate. Defective items are usually observed in the inventory management system, including the possibility of being defective during shipment from the warehouse to the supplier or retailer’s warehouse. Several researchers [6] considered defective rate of items in the inventory system to develop their models. This article presents an inventory model that handles the rate of deterioration in such a way that the maximum useful life of the product is prioritized, and this study also considers the defective rate.

Demand is critical to inventory management, and future inventory cannot be forecast without considering demand. For each item, the demand rate is different. The demand rate for some items rises at the start of the cycle. The demand rate for some items remains constant throughout the cycle, whereas the demand rate for others increases at the end. The demand rate for cooked items such as bread, sweets, cakes, and so on increases at the beginning of the cycle because customers love just-made goods. Due to the expiration date, the demand rate decreases at the end of the cycle for goods such as fish, vegetables, fruits, and so on. On the contrary, others have an increasing demand rate at the end of the cycle, such as household goods such as oil, sugar, milk, and so on, and the demand rate is constant for furniture, electrical goods, etc. Khedlekar and Sukhla [7] developed a dynamic pricing model for logarithmic demand. Smaila and Chukwu [8] introduced an EOQ model with quadratic demand trends and quasi-partial backlogging. Dutta Choudhury et al. [9] presented an inventory model by taking two-component demand. Prasad and Mukherjee [10] developed an inventory model on stock and time-dependent demand. Wu et al. [11] developed an inventory policy for trapezoidal-type demand patterns and maximum lifetime under trade credits. Mahapatra et al. [12] presented an EPQ model with demand and reliability-dependent unit production cost under limited available intuitionistic fuzzy type storage space. Shaikh and Mishra [13] developed an EOQ model for price-sensitive quadratic demand and inflationary conditions. Delay in payment is a common phenomenon in the inventory management system. A delay in payment is offered by the supplier, where the retailer’s purchase cost is paid at a later date without an interest charge. Liao et al. [14] proposed an inventory model with items that deteriorate under inflation and allow a delay in payment. Several researchers [1518] developed their inventory models under the condition of delay in payment. Some researchers used other payment schemes to develop their models. Teng et al. [19] presented an inventory model under a progressive payment strategy. Teng et al. [20] presented an EOQ model under a trade credit financing scheme. This study of inventory management allows a delay in payment for two different situations based on time.

Neutrosophic numbers explain the impreciseness of the systems. In real life, most of the parameters are uncertain, so in this situation, neutrosophic numbers play a crucial role in overcoming uncertainty. Mullai and Surya [21] developed a price-break EOQ model with neutrosophic demand and purchasing cost as triangular neutrosophic numbers. Mariagrazia et al. [22] proposed a supplier selection methodology under uncertainty. Ge and Zhang [23] presented an inventory model under a fuzzy uncertain environment. De et al. [24] developed an inventory model for the non-random uncertain environment using the neutrosophic fuzzy approach. Mariagrazia et al. [25] introduced the fuzzy technique for the supply chain network with quantity discounts. Das and Tripathy [26] investigated some properties of neutrosophic multiset topology. Bonilla-Enriquez et al. [27] proposed a supply chain model taking into account uncertain demand. The present paper considers the cost pattern of the inventory system as a triangular neutrosophic number.

Generally, inventory problems are solved by traditional direct optimization methods, but one of the shortcomings is that they are very often stuck to local optima. To avoid several shortcomings of global optimization, particle swarm optimization (PSO) is very useful in solving problems in the field of inventory control. Biuki et al. [28] presented an inventory problem in optimizing through two hybrid metaheuristics as parallel and series combinations of genetic algorithm and PSO. Alejo-Reyes et al. [29] introduced an inventory model for supplier selection and order quantity allocation by using metaheuristic algorithms. Rau et al. [30] proposed a multiobjective green cyclic inventory routing problem via the discrete multiswarm PSO method. Al-Khazraji et al. [31] applied multiobjective PSO to optimize production inventory control systems. Patne et al. [32] presented a closed-loop supply chain network configuration model using game-theoretic PSO. Dabiri et al. [33] presented a bi-objective inventory routing problem with a step cost function using multiobjective PSO. Manatkar et al. [34] presented an integrated inventory distribution optimization model for multiple products using a novel hybrid multiobjective self-learning PSO method. Srinivasan et al. [35] applied PSO to optimize a mathematical model with defective goods. A multiitem EPQ model with a production capacity restriction using the PSO algorithm was developed by Pirayesh and Poormoaied [36]. Araya-Sassi et al. [37] developed an inventory location problem by using the Lagrangian relaxation algorithm. This article considers the PSO algorithm to find the optimal solution under the neutrosophic uncertainty-based cost metrics of the presented inventory system. A comparison of the contributions is presented in Table 1 to determine the most significant contribution of this work. The contribution of this study is to offer the best cycle duration and reduce the total cost of the inventory system under triangular neutrosophic cost pattern when the customer demand depends on time and reliability. To the best of our knowledge, this work addresses all of the themes stated in Table 1 simultaneously first time, which have never been examined together in the literature before.

Simultaneous evaluation of the preceding assumptions allows us to represent a wide range of real-life situations, resulting in a more realistic inventory model. The remainder of this article is organized as follows. The inventory problem is mathematically formulated in Section 2. Section 3 develops the model formulation in a neutrosophic environment. Section 4 contains the proposed inventory model using the PSO algorithm. In Section 5, a numerical example is explained, and then a sensitivity analysis of the optimal inventory policy for the system input parameters is presented, along with some significant managerial insights obtained from the results. Section 6 concludes with some findings and research directions for the future.

2. Mathematical Modeling of the Inventory System

Based on the discussion in Introduction, the following notations are considered to present the inventory system.

2.1. Notations

: ordering cost per order ($/order).: purchasing cost per item ($/unit).: selling price ($/unit).: holding cost per item ($/unit/unit time).: rate of interest charged per year in stocks by suppliers.: rate of interest earned by investment per year.: deterioration rate of items.: maximum life time in years of item, .: supplier permissible delay period.: cycle time in per cycle.: percentage of non-defective items .: reliability of item .: reduction percentage of purchasing cost of defective items.: average demand per cycle .: total demand per cycle.: demand pattern index .

2.2. Problem Definition

This paper considers mathematical modeling of the fact that, in reality, the rate of deterioration increases over time, where the maximum lifetime of an item is known to the suppliers. The shortages are not allowed for the proposed inventory system through the maximum life span is considered for the item in each replenishment cycle. Demand is the main factor in inventory modeling, which depends on several factors, but reliability is one of the main characteristics. Therefore, the demand rate of the commodities for this study depends on the reliability and the power pattern of time, i.e., the demand rate in a period of time is defined as , where n seems to be the demand pattern index, and , and it is worth noting that if , a larger portion of demand occurs towards the beginning of the period. If , the demand rate remains constant throughout the inventory cycle. However, if , a large portion of the demand comes near the end of the cycle, as shown in Figure 1(a).

The deterioration rate of the inventory system increases with time. It depends on the maximum shelf life of the product in the inventory system, so the deterioration rate of this model is  = . The supplier offers a reduction rate on the purchase cost of damaged or defective products of this model to the buyer. The inventory system considers a permissible delay in payments, where the credit period is less than or equal to during which the buyer does not need to pay any interest to the retailer, where the interest is applied beyond the period . The costs related to the inventory system are imprecise. Therefore, to represent the impreciseness of the inventory costs, we consider the inventory cost parameters as triangular neutrosophic numbers instead of a fixed crisp value. The retailer separates defective and non-defective items by screening items received from the supplier. The retailer fixed the rate of the range of non-defective items which is , which means that the retailer does not accept stock that contains 10% or above defective items. The retailer also receives a discount on the purchase cost of defective items.

2.3. Mathematical Formulation

Based on these considerations, the inventory system (shown in Figure 1(b)) can be represented in the time interval , using the differential equation as given below:subject to the boundary conditions: , where is the initial inventory for each cycle.

Lemma 1. The deterioration rateis decreasing with maximum life time of the product and increases with time.

Proof. The deterioration rate of the proposed inventory system is , where .
for .
Therefore, the deterioration rate decreases with the maximum shelf life of the product.
for .
Therefore, the deterioration rate decreases with the maximum shelf life of the product and increases with time.
Hence the proof.

Lemma 2. The demand rate D(t) is decreasing with time for, increasing with time for, and constant with time for.

Proof. The demand rate of the given model is . For time t, , for and for all t, which shows that demand rate increasing with time.
For the time t, , for and for all t, which shows that demand rate is constant with time.
For the time t, , for and for all t, which shows that demand rate is decreasing with time.
Hence the proof.
Using these boundary conditions, the solution of differential (1) of the proposed inventory model during the time interval is as follows:Again using boundary condition , we getNow, using (2) and (3), the holding cost (HC), the purchasing cost (PC), and the ordering cost (OC) of the proposed inventory model can be obtained as given in the succeeding equations:The purchase cost of the inventory system consider per item and required to purchase the item for inventory.The cost of ordering is fixed for each cycle of unit inventory model.The supplier gives a discount to the purchasing cost of damaged or defective items of the inventory. The discount amount obtained by purchasing defective items isThe retailers of this inventory system allow for a conditional delay in payment for the buyer. Therefore, the two cases are arrived at such that the delay period is greater than the cycle time and delay period is less than the cycle time. Now, we discuss the interests due to delayed payments for both cases to obtain the total inventory costs of the following two cases.

Case 1. Delay period is greater than the cycle time .
Here, the retailer earns interest per cycle at a return rate ; if , then the annual interest earned per cycle is given byThe total cost of inventory per unit time in situation is obtained as follows:The total inventory cost of the proposed inventory model under the delay period is greater than the cycle time which is to be optimized with respect to the optimal cycle using PSO.

Case 2. Delay period is less than the cycle time .
The interest charged (IC) by the retailer per cycle is obtained byThe interest earned during the time 0 to is given byIn situation , the total cost of inventory per unit of time is obtained as

3. Proposed Inventory Model under Neutrosophic Environment

Neutrosophic number deals with uncertainties better than fuzzy and intuitionistic fuzzy numbers, since neutrosophic numbers contain truth, hesitant, and falsity membership functions to deal with all types of uncertainties of the parameters. In a real market, the cost parameters are uncertain and contain a dilemma in the decision maker’s mind. Thus, this study attempts to manifest the inventory system by introducing a neutrosophic set to represent the different inventory costs and rates of the proposed inventory system. Thus, observe the effect of uncertainty by comparing it with the crisp model. In the proposed inventory model, this paper has considered holding cost , purchase cost , ordering cost , and selling price (S) as triangular neutrosophic numbers [49]. The representation of triangular neutrosophic numbers , , , is as follows: , , , .

Now using the removal area technique for the de-neutrosophic of the neutrosophic numbers to find the optimum inventory costs, are obtained as follows: , , , and .

To obtain the total cost in the neutrosophic domain , substituting the values of de-neutrosophic values , , , and into the total cost of both cases in (9) and (12), we get

Next section finds the optimal value of equations (9) and (12)–(14) through the PSO technique.

4. Optimization of Proposed Inventory Model Using PSO

Particle swarm optimization is a nature-influenced optimization technique to solve the most complex optimization problems. The PSO developed by Kennedy and Eberhart [50] is influenced by the movement and intelligence of swarms. The swarms are considered as vector points in the domain space where the optimum value of a given objective function lies. There are four important vectors: (i) x-vector, which records the current position, (ii) p-vector that records the personal best, (iii) v- vector, which controls the velocity of the moving particle at each instance, and (iv) g-vector, which records the direction of the particle towards the global best position. The design of the algorithm is that the particles at each iteration move from its current position towards a new position by traveling to some extent parallel to its velocity vector, then to some extent towards its personal best, and to some extent towards the global best. The new position of the particle is the sum of these three vectors. The global best position is found by taking the best of all the personal best positions of all particles. The attributes are updated and proceeded for the next iteration. The process is repeated for a suitable number of iterations to achieve the best solution.

Four essential parameters control the movement of particles; a parameter to represent the coefficient of inertia “w” along with a damping inertial coefficient “wamp,” and two constants, representing the acceleration of individual swarms, and representing the social acceleration. These parameters are varied according to the choice of the optimization problem. To find the optimal total cost of the proposed inventory model, the algorithm of the variant of the PSO technique is presented as follows (see Algorithm 1).

(1)Read: Swarm Population (SwarmP), Maximum Iteration (MaxI), particle(i).Velocity , particle(i).Position , particle(i).Best.Position , Global Best.Position , particle(i).Cost , particle(i).Best.Cost , Global Best.Cost ,
(2)Set the objective function:
(a); (Optimum for ),
(b); (Optimum for ).
(3)Set the parameters of PSO for Inventory system:
(a), , , , .
(b) Set coefficient of inertia, ; the damping ratio of inertial coefficient , , and .
(c) Set inventory parameters as .
(4)Define the structure of void particles:
(a) Set , , , , and .
(b) Create vector of void particles of size equal to SwarmP.
(c) Set initial  =  .
(5)Initialize the population member: For to SwarmP initiate and .
(6)Estimate and evaluate:
(a) Set , find , and .
(b) then .
(c).
(7)Optimize and update for every iteration:
(a) For to MaxI, and for to SwarmP update the velocity and position.
(b).
(c).
(d) If , then set ; .
(e) If , then set ;
(8)Set up the global best cost:
(a).
(b) Set and continue the iteration by going back to 6.
(9)Stop

5. Numerical Solution

Numerical example of the proposed inventory system is presented to support analytical derivation and discussion. Let us consider the proposed inventory system with the following setup: assume that the ordering cost is $750 per order, cost of raw material is $10 per unit, the cost to hold the item is $2.5 per item per unit time, and the selling price of each item is $18. Since the system allows both defective and perfect items, let the percentage of non-defective items be taken as 0.95, while the maximum lifetime of the item is considered as three years. Let the supplier admit a permissible delay of the payment period of two years for which the buyer earns 5% interest, where the interest charged per year in stocks by the supplier is taken as 8% per year. Let us further consider the parameter settings for the PSO:  =  , , , , years, years (Case 1), year (Case 2), , , , , , , , and .

Use the PSO for the inventory system to obtain the optimal solution of the proposed inventory model, as shown in Table 2.

5.1. Numerical Example Using Triangular Neutrosophic Set through PSO

A numerical example has been given to test the proposed inventory system using a triangular neutrosophic set to show the effect of the impreciseness of the cost parameters , , , . Thus, by de-neutrosophication of the triangular neutrosophic parameters, the values of the cost parameters are , , , . For , years, years (Case 1), year (Case 2), , , , , , , , and . The optimal solution in the neutrosophic environment is shown in Table 2.

The PSO is implemented to find the optimal total inventory costs , and with neutrosophic cost parameters along with the parameters of the proposed PSO algorithm , , , . Table 2 shows a higher total cost at because the demand rate increases with time at , so the total quantity is higher and the total cost of inventory systems increases. The total cost is lower at because the demand rate decreases with time at ; therefore, the total quantity is lower, and consequently, the total cost of the inventory systems decreases. Estimating the values of the cost parameters under actual market conditions is a big task for decision makers. The advantage of using neutrosophic numbers is that the decision maker can choose appropriate values for the cost parameters based on real market conditions. Neutrosophic numbers contain truth, hesitance, and falsity membership functions. These membership functions are helpful for establishing cost parameters to overcome impreciseness. Using neutrosophic numbers, the decision maker can choose appropriate values for cost parameters based on real market conditions.

By changing the values of and in the numerical example, the optimal solution was obtained as shown in Table 3. Both total inventory costs and increase when and decrease when . If , then the cost parameters increase and the total cost also increases. If , then the values of the cost parameters decrease and the total cost also decreases. Here, and are the smallest positive constants between (0, 1) for the proposed inventory cost parameters. For all values of and , the optimal inventory costs in the neutrosophic domain are lower than those in the crisp environment.

5.2. Sensitivity Analysis

Based on the numerical example of the inventory system shown in the previous section, a sensitivity analysis is performed by changing the parameters by a specific percentage, depending on the parameter boundary, based on the sensitivity analysis performed. For the sensitivity analysis, consider one parameter at a time and keep the other parameters fixed.

Table 4 shows that total inventory costs and are more sensitive to the purchase cost and moderately sensitive to the holding cost and the ordering cost but less sensitive to the selling price . In both and , the holding cost changes by 20%; if the holding cost increases, then also increases in an analogous manner obtained from Figure 2. The total inventory costs and are less sensitive to the selling price, which can be seen in Figure 2. It is also observed that and are moderately sensitive to the ordering cost as shown in Figure 2. and are highly sensitive to the purchase cost, since the total cost is directly related to the purchase cost. Furthermore, when reliability is high, the price of buying an item dependent on reliability also increases.

Table 5 shows that the total costs and are highly sensitive to d and moderately sensitive to m, f, and r but less sensitive to , , and . Both total costs and are highly sensitive to d; according to credit concepts, from a practical point of view, if the parameter d increases, then total costs also increase, as shown in Figure 3. Both and are moderately sensitive to m, f, and r, which is graphically shown in Figure 3. The parameter f is increased as well decreased by 5% and 2%, because if a percentage of f increases by more than that value of f exceeds its permissible value, which shows that and decrease a limited amount but not very much, if non-defective items of inventory increase. The total optimal costs and of the proposed inventory model decrease with a lower rate, as shown in Figure 3 if the parameters and increase.

5.3. Managerial Implications

The value of any product is generally time-dependent. The cost of several products will increase as time progresses and the same will decrease for a few specific products. Most food items are perishable for which the value of such products will decrease, prolonging the storage period. Like wine, the cost of a few select products will increase when the storage period increases. This research convinced decision makers that perishable products should be stored in a healthier environment to lengthen their shelf life and reduce their rate of deterioration.

From a managerial standpoint, this study contributes to a better understanding of the reliability of the product. A well-maintained product will have higher demand, resulting in higher profits during sales. Furthermore, the managerial fact is that in order to keep such a product, management may have to incur higher holding costs. This means that we must balance these factors so that our total inventory cost is kept to a minimum.

The efficiency of deteriorating items, the cost patterns, and the reliability factors have been protectively represented in the flow of the inventory system. Therefore, this research can be widely applied in managing the inventory system where the reliability and deteriorating rates are concerned. Consequently, the current inventory model was created to depict these realistic features that may be utilized to forecast the system’s elements.

6. Conclusion and Future Direction

The proposed inventory systems consist of various essential and complex functions to deal with the ordering and storing of products in which there are several uncertainties. Finding a generalized optimum order quantity model with uncertainty is a challenging task. In this article, an optimum order quantity inventory model without shortages in determining the items with reliability, the time-dependent power pattern, and the allowed time delay. This improves the inventory system, addresses the nature of items for time-dependent deterioration factors, and the cost pattern imprecise for a real scenario. The key advantage of this model is that it covers many sorts of product demand rates. The uncertainty in the cost pattern is dealt with the concept of impreciseness using triangular neutrosophic number. This study discussed and formulated the inventory model with reliability-dependent demand and analyzed the solutions by altering various parameter settings. A modified PSO algorithm is presented that fits the inventory model, and the same is implemented to find the optimal total cost in two different scenarios relating to time delay. PSO is utilized in both the crisp and neutrosophic senses, and the neutrosophic environment in which the solution is obtained is analyzed by sensitivity analysis. The proposed inventory system can be developed multiitems with different maximum lifetimes for items. The inventory system under a probabilistic or possibilistic environment with a different form of uncertainty such as intuitionistic fuzzy, interval type-2 fuzzy, Hesitant fuzzy sets, etc. The main limitation of this model is that it cannot allow more than 10% of defective items. Shortages are not allowed in this model, but shortages generally occur in the inventory system. This is the drawback of this model, so it can further develop this model under shortages.

Appendix

Definition A.1. Single-valued neutrosophic set: “a single-valued neutrosophic setof a single-valued independent variable (x) is defined by=, whererepresent the concept of truth, hesitation, and falsity membership functions, respectively.”
Neutro-normal: let us consider three points p, q, and r for which ; then, is defined as neutro-normal.

Definition A.2. Neutro-convex: “a neutrosophic setis called neutro-convex if the following condition holds: (i), (ii), and (iii), where, and.”

Definition A.3. Triangular single-valued neutrosophic number: “a triangular single-valued neutrosophic numberis defined as. Here,is the truth membership function,is the hesitation membership function, and the falsity membership function is, where the membership functions are mathematically and graphically (Figure 4) defined as follows.”De-neutrosophication of triangular single-valued neutrosophic number: use the removal area technique to evaluate the de-neutrosophication value of triangular single-valued neutrosophic number . The de-neutrosophic form of is given by  =  .

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.