Abstract

Nowadays, it is well known that global warming is a great hazard to the planet, and the carbon emissions are a principal source of global warming. For this reason, the customers have become more environment and quality conscious than before, and as a result, they request the firms to be ecofriendly. In this context, it is desirable that companies develop and implement inventory models which consider sustainability issues. Furthermore, the companies face problems of shortages and setting prices in order to persist in a competitive and challenging business. Besides, there exists a kind of items different than the traditional products that it is necessary to feed them until a target weight is reached in order to slaughter and sell to customers. These are named as growing items. In this sense, this research work proposes an inventory model for growing items with imperfect quality when the demand is price sensitive under carbon emissions and shortages. The shortages are fully backordered. The demand is price sensitive according to a polynomial function. The proposed inventory model determines jointly the optimal policy for the selling price of perfect-quality growing items, the order quantity, and the backordering quantity which maximize the expected total profit per unit of time. Some numerical examples are resolved in order to illustrate the use and the applicability of the inventory model. Finally, a sensitivity analysis is conducted and some managerial insights are given.

1. Introduction

Inventory management is a valued function for all firms. The main purpose of the inventory management department is to control the materials from the acquisition to sale, taking decisions of how much and when to buy the items avoiding excess or unplanned stockout.

Since the introduction of the classical economic order quantity (EOQ) inventory model proposed by Harris [1], industries have shown great interest in the development of inventory control models with new features. For example, Ghosh et al. [2] built an inventory model for a deteriorating production system when there are random machine breakdowns considering that there exists a selling price discount. Manna et al. [3] studied the impacts of inspection errors on an imperfect production-inventory model when the demand rate depends on both selling-price discount and warranty time. Khan et al. [4] and Shaikh et al. [5] introduced a two-warehouse inventory model for instantaneously deteriorating goods considering an advance payment agreement. In the same research direction, Khan et al. [6] investigated the effects of noninstantaneous deterioration on the ordering decisions for a two-warehouse inventory system taking into consideration an advance payment scheme. Shaikh et al. [7] developed an EOQ inventory model for deteriorating goods with stock-dependent demand when there exists a price discount facility.

It is well known that the EOQ determines the optimal order quantity to minimize total inventory costs. Although this inventory model has the groundwork for inventory systems, it has unrealistic assumptions. For that reason, many researchers have incorporated different characteristics in order to model real world applications such as imperfect quality, carbon emissions reduction, growing items, shortages, and pricing, among others.

Harris’s [1] inventory model assumes that all items are of perfect quality. However, this assumption not always occurs because the majority of the manufacturing processes generate a portion of imperfect quality items. So, Salameh and Jaber [8] formulated an EOQ inventory model considering that the ordered lot contains imperfect-quality items and these items are sold as a single batch at a lower price. Cárdenas-Barrón [9] identified and corrected an error in the inventory model presented by Salameh and Jaber [8]. Wee et al. [10] extended Salameh and Jaber’s [8] inventory model by incorporating shortages with full backordering. Afterwards, Maddah and Jaber [11] detected and rectified one flaw in Salameh and Jaber’s [8] inventory model. The flaw is in the manner of determining the expected total profit per unit of time . Salameh and Jaber [8] utilized the equation which is not exact. The research work of Wee et al. [10] contains the same flaw. It is important to remark that the process of producing the profit is a renewal process with renewal points at order placement periods. Therefore, in order to calculate the expected total profit per unit of time, it is needed to apply the well-known elementary renewal reward theorem given by Ross [12]. For this reason, the expected total profit per unit of time must be computed with the following equation: . Chang and Ho [13] revisited and rescued Wee et al.’s [10] inventory model by applying the well-known elementary renewal reward theorem to get the expected total profit per unit time. Cárdenas-Barrón [14] derived closed-form expressions to determining the optimal solution to an EOQ inventory model considering items with imperfect quality. Ghiami and Beullens [15] developed a production-inventory model for a deteriorating product with partial backordering using a cash-flow net present value analysis. They generated managerial insights related if it is convenient financially to have planned shortages with partial backordering and lost sales. Zhou et al. [16] developed an economic-order quantity model considering shortages, imperfect quality, and inspection errors. Taleizadeh et al. [17] studied an inventory model when the lot contains imperfect-quality items, and these cannot be replaced with perfect-quality ones immediately, but these are repaired and then they are sold at lower price. In this line of research, Rezaei [18] introduced an inventory model considering imperfect-quality items. Sampling inspection plans are designed so the buyer would be able to decide regarding the next step.

At the present time, as we know, global warming is a big threat to the world and carbon emissions are a leading cause of global warming. For this reason, the sustainability is an extremely significant concern for all people that desire to maintain the planet healthy. It is well known that industries generate large amounts of carbon emissions. Furthermore, the transportation of goods through the supply chains produces carbon emissions too. Consequently, these damage the environment and cause global warming at the same time. In this context, the governments of the countries have imposed strict policies that limit the amounts of carbon emissions with the aim of having a sustainable environment free of pollution. The regulations related to avoid pollution impose additional costs to firms; thus, it is convenient for these to adapt and obey the regulations. In this direction, several companies have adopted inventory models that consider reducing carbon emissions. For example, Hua et al. [19] investigated how companies manage carbon emissions in an inventory management. Arslan and Turkay [20] integrated sustainability features into the design of an inventory model and proposed different policies that include environmental and social criteria. Battini et al. [21] introduced a sustainable EOQ inventory model integrating factors that have a great impact on the environment. Lin and Sarker [22] presented an inventory model that considers carbon tax policy and imperfect-quality items avoiding shortages and incorporated and evaluated some carbon tax systems. Tiwari et al. [23] investigated an integrated vendor-buyer inventory model for deteriorating goods with imperfect quality, taking into account that all processes in the supply chain cause carbon emissions. The target of their inventory model is to decrease the ecological effects. Kazemi et al. [24] proposed an inventory model with imperfect quality from a sustainable point of view, and the aim is to know the impact of emission costs into the total profit. Modak et al. [25] addressed a manufacturer-retailer supply chain inventory model by considering that there is a cost of GHG emission of the production system and determined the optimal pricing policies. Afterwards, Sinha and Modak [26] presented an economic production quantity (EPQ) inventory model taking into account aspects of carbon emissions and carbon trading. Li and Hai [27] studied an inventory system for a warehouse with multiple retailers considering carbon emissions. Manna et al. [28] constructed a production-inventory model for controlling the GHG emissions when the pollution parameters are fuzzy. Later, Manna et al. [29] studied the effects of carbon emissions on an imperfect production-inventory model. Huang et al. [30] proposed an inventory model for a two-echelon supply chain considering that all operations within the chain produce carbon emissions. Their inventory model helps companies to determine the optimal production quantity and delivery quantity which minimize the total costs under different carbon emission policies. Recently, Mishra et al. [31] revised the standard EOQ inventory model to integrate sustainability matters and developed three inventory models with and without shortage situations. Medina-Santana and Cárdenas-Barrón [32] formulated an inventory model considering a discontinuous transportation cost function and carbon emissions function. The most recent investigation about carbon emissions belongs to the authors Modak and Kelle [33]. They proposed to apply the social work donation as a scheme of corporate social responsibility into a closed-loop supply chain, taking into consideration the carbon emission tax and that demand is uncertain.

In recent years, a new topic, that is, inventory management for growing items, has also been included in the derivation of inventory models. A growing item is a kind of product which is capable of growth during time like farmyard animals, unlike inert products that do not increase in weight during storage. Growing items are incorporated for first time in an EOQ inventory model by Rezaei [34]. Here, the weight increment of items is the main difference between the inventory models proposed by Rezaei [34] and Harris [1]. Rezaei [34] developed an EOQ for growing items for different types of animals, specifically for a kind of poultry. Zhang et al. [35] revisited and extended the research work of Rezaei [34] by creating an inventory model for growing items reducing the carbon emissions in operations. Sebatjane [36] presented three inventory models for growing items. Afterwards, Khalilpourazari and Pasandideh [37] worked with growing items too. They formulated a mathematical model for multi-item economic order quantity subjected to some operational constrains. Nobil et al. [38] derived an EOQ inventory model for growing items when the shortages are allowed, and these are fully backordered. They argued that their inventory model helps poultry industries. Under an environment of a two-echelon supply chain composed by one supplier and one farmer, Malekitabar et al. [39] developed an inventory model for growing mortal items considering both the growth function and the mortality rate for the items. Recently, Sebatjane and Adetunji [40] presented an EOQ inventory model for growing items with imperfect quality. It is assumed that a certain fraction of the growing items has lower quality than desired. They mentioned that it is necessary to define the growth function of the items in order to calculate the feeding cost. In this context, Sebatjane and Adetunji [40] considered three growth functions: logistic, linear, and split linear. In the same year, Sebatjane and Adetunji [41] introduced an economic-order quantity inventory model for growing items with incremental quantity discounts. They proposed an optimal inventory policy to minimize the total inventory cost in both the owned and rented facilities. Sebatjane and Adetunji [42] continued working with growing items, but in this case, they formulated an inventory control model more realistically, considering that the items need to be transformed and packaged into a consumable form before customer demand is met. The next year, Sebatjane and Adetunji [43] built a model for managing inventory in a perishable food products supply chain that begins with farming operations where growing items live and finishes with the consumption of processed inventory. Sebatjane and Adetunji [44] derived an inventory model for a four-echelon supply chain with farming, processing, screening, and retail operations. Gharaei and Almehdawe [45] provided an economic growing quantity (EGQ) inventory model to determine the optimal economic growth and slaughter period and the economic growing quantity to minimize the total cost of the inventory system. Hidayat et al. [46] presented an EOQ inventory model with a capacitated warehouse facility and limited budget for growing products when the seller provides an incremental quantity discounts scheme. Mokhtari et al. [47] addressed a production-inventory model for growing goods which deteriorate through time. Nishandhi [48] studied an EOQ inventory model with budget-capacity constraint for growing items when a portion of the items are of imperfect quality. Pourmohammad-Zia and Karimi [49] determined the optimal replenishment and breeding policies for growing products. Afzal and Alfares [50] and Alfares and Afzal [51] developed EOQ inventory models for growing items considering shortages with full backordering. Recently, Sebatjane and Adetunji [52] created an inventory model for growing items under a three-echelon supply chain environment taking into account farming, processing, and retail echelons. Table 1 shows the inventory models related with growing items.

The organizations need to determine the optimal selling prices for the products in order to encourage that the customers buy more and more products, and therefore, these organizations can survive in the competitive business environment. In this context, the academicians and researchers are modelling the demand rate as dependent of price. For example, Khan et al. [53] analyzed two supply chain inventory systems when the demand rate depends on price. Khan et al. [54] proposed an inventory model for deteriorating merchandises when the demand is price sensitive, and there exists a discount policy according all-units arrangement. Khan et al. [55] developed inventory models for perishable products when the demand rate is dependent on both price and advertisement. On the one hand, Khan et al. [56] examined the effects of an advance payment with discount facility on ordering decisions for perishable items taking into account that the demand rate is both price and stock dependent. On the other hand, Panda et al. [57] dealt with a two-warehouse system for deteriorating items when the demand rate is dependent on both price and stock; however, the price is not optimized. Sinha et al. [58] formulated an entropic-order quantity inventory model when the demand of the product is dependent on selling price and there is an inspection process to split the imperfect-quality products from the perfect ones. Modak [59] and Modak and Kelle [60] introduced omni- and dual-channel supply chain models, respectively, by considering that there exists a price- and delivery-time-sensitive stochastic demand.

Nowadays, the firms face problems of shortages and setting the selling price in order to survive in a competitive business, which is becoming more challenging day by day. Moreover, the customers are more quality and environmental conscious than before, and therefore, they demand products of perfect quality and request firms to minimize the carbon emissions to the environment. However, the firms always face problems of process quality and contaminate the planet. Besides, there exists a type of goods different than the conventional products which requires feeding until a target weight is reached in order to slaughter and sell to customers. These are named as growing items (chickens, cows, pigs, goats, fish, shrimps, etc.). In this direction, this research work develops an inventory model for growing items with imperfect quality when the demand is price sensitive under carbon emissions and shortages. The shortages are fully backordered. The demand is price sensitive according to a polynomial function. The proposed inventory model determines jointly the optimal policy for the selling price of perfect-quality growing items, the order quantity, and backordering quantity which maximize the total profit per unit of time.

The remaining parts of this research work are described as follows: Section 2 introduces the notation defining symbols for the parameters, decision variables, decision-dependent variables, and functions. Section 3 states the assumptions under which the inventory model is built. Section 4 formulates the inventory model with shortages for growing items with imperfect quality when the demand is price sensitive considering carbon emissions. Section 5 determines some theoretical results, develops the solution procedure to determine the optimal solution to the inventory problem, and identifies some inventory models as special cases. Section 6 solves some numerical examples. Section 7 provides a sensitivity analysis and some managerial insights. Section 8 gives some conclusions and summaries several research points that can be addressed in the near future.

2. Notations

With the aim of having a standard notation for the inventory models with growing items, the nomenclature of Sebatjane and Adetunji [40] is used and extra symbols are defined too. Therefore, the following notations are adopted in order to develop the inventory model for growing items with imperfect quality, carbon emissions, and planned shortages.

2.1. Parameters
: Scale parameter for the price-dependent demand: Sensitivity parameter for the price-dependent demand: Demand power index: Selling price of imperfect items (currency symbol/unit of weight): Setup cost (currency symbol/cycle): Holding cost (currency symbol/unit of weight/unit of time): Shortage cost (currency symbol/unit of weight/unit of time): Feeding cost (currency symbol/unit of weight): Purchasing cost (currency symbol/unit of weight): Inspection cost (currency symbol/unit of weight): Carbon tax rate (currency symbol/amount of carbon emissions): Carbon emissions cost (currency symbol): Amount of carbon emissions produced during the setup process (unit of weight/unit of time): Amount of carbon emissions caused by holding items in the warehouse (unit of weight/unit of time): Amount of carbon emissions generated during the feeding period (unit of weight/unit of time): Amount of carbon emissions made during the purchasing activity (unit of weight/unit of time): Amount of carbon emissions created during the inspection process (unit of weight/unit of time): Inspection rate (unit of weight/unit of time): Asymptotic weight of each item (unit of weight): Integration constant (numeric value): Growth rate (numeric value/unit of time): Percentage of slaughtered items that are of imperfect quality : Expected value of the percentage of imperfect items : Expected value of the percentage of perfect items : Weight of a newborn item (unit of weight): Target weight of a grown item (unit of weight): Weight of an item at time (unit of weight): Growing period (unit of time): Inspection period for the backordering quantity (unit of time): Inspection period for units of weight (unit of time): Consumption period of perfect items after inspection time (unit of time): Shortages accumulation period (unit of time)Decision variables:: Order quantity of newborn items (units): Backordering quantity (unit of weight): Selling price of perfect items (currency symbol/unit of weight)Decision-dependent variables:: Cycle time (unit of time): Total weight at the beginning of the growing period, (unit of weight): Total weight at the end of the growing period , (unit of weight)Functions:: Price-dependent demand function (unit of weight/unit of time): Growth function: Probability density function of the percentage of imperfect items: Total profit (currency symbol/unit of time)

3. Assumptions

The inventory model is based on the following assumptions:(1)The planning horizon is infinite and a single kind of items is purchased, and these are capable of growing before the slaughter process.(2)The shortages are permitted, and these are completely backordered.(3)The items are slaughtered and are immediately inspected in order to sell them to consumers. Firstly, the backordering quantity is inspected in order to cover the shortages of the previous cycle.(4)There exists an inspection process that is 100% effective.(5)A random percentage of the slaughtered items is of imperfect quality.(6)Imperfect-quality items are not reworked or replaced.(7)All imperfect-quality items are salvaged and sold as a single lot at the end of the inspection process.(8)The feeding cost for growing the items is directly related to weight gained by these.(9)The holding cost to keep a weight unit of the slaughtered item in storage is incurred during both the inspection process and the consumption period.(10)The demand rate is a polynomial function of selling price of the perfect-quality items. It is as follows: .(11)The selling price of perfect-quality items is optimized, and it must be greater than that of the imperfect-quality items.(12)Carbon emissions are taken into account, and these occur in all operations of the inventory system, except in the shortage period.

4. Inventory Model Development

The inventory model for growing items with imperfect quality, carbon emissions, and shortages is depicted in Figure 1. Consider a situation where a company orders newborn growing items from an outside supplier at the beginning of the growing period . Each newborn growing item has an initial weight of . In this moment, the total initial weight of the inventory is . The growing items are fed, and eventually, they grow through time until an objective weight of is attained. Then, the growing items are slaughtered at the end of the growing time . At this point, the final weight of the inventory is , and this total weight contains a percentage of imperfect items. The portion of imperfect items is a random variable with a known probability density distribution , and its expected value is . The shortages are permitted, and these are fully backordered. Therefore, immediately, the inspection process starts to screen the items to complete the backordering quantity at a rate of units of weight per unit of time during the period in order to satisfy immediate shortages from the previous cycle. So, at the end of the inspection period , the inventory model diminishes vertically by units of weight. It is important to remark that the items continue to be inspected at the same rate during the period till the total weight is screened. The length of the inspection time is . It is worth mentioning that, for the duration of , the on-hand inventory declines by both removing the imperfect items and current demand rate. At the end of period , the imperfect items are salvaged and sold as a single lot with a less price. Consequently, the on-hand inventory drops vertically by units of weight. On the other hand, in the course of , the on-hand inventory decreases due to the current demand rate. The inventory level continues gradually consuming until it reaches zero at the end of period and the shortages period starts. As a final point, during , the shortages are accumulated at the current demand rate which are eventually satisfied in the next cycle. Without loss of generality and with the purpose to make the mathematical expressions more tractable, the cycle time is determined as follows: .

Throughout , the items are growing according to a logistic growth function which relates the weight of items with time using three input parameters. These input parameters are the asymptotic weight of the items, the integration constant, and the growth rate, which are represented by , , and , respectively. Thus, the logistic growth function of the items is mathematically expressed by

As it was mentioned above, the growing items are slaughtered when their weight attains the objective weight of which occurs at the end of the growth period . Hence,

The duration of the growth period is calculated by solving equation (2) for . Thus,

The inventory must be inspected before being sold in order to avoid to vend imperfect quality as good ones; therefore, firstly, the backordering quantity of the previous cycle must be screened at an inspection rate . Therefore, the duration of the inspection period is computed as follows:

After the inspection of the backordering quantity is conducted, the screening process continues until the total weight is screened due to the fact that there are pending units of weight to be inspected. This is performed during the inspection time which is given as follows:

After the inspection time , the imperfect-quality items are withdrawn from the storage and sold. Now, the on-hand inventory contains only perfect-quality items, and these are consumed during . The time is determined as follows:

The shortages accumulation period is obtained with

With the aim to have more manageable equations and without loss of generality, it is defined that the expected cycle time is computed by the sum of , and . From equations (5)–(7), the expected cycle time is calculated as

The main goal of the inventory model for growing items with imperfect quality, carbon emissions, and shortages is to determine the optimal values for the selling price, order quantity, and backordering quantity that maximize the total profit of the company. The total profit is obtained by the difference between the total revenue and the total cost of the inventory system. The details of how to calculate each component of the total revenue and the total cost are given below.

4.1. Expected Revenue per Period

Due to the fact that the company vends both perfect- and imperfect-quality items, for that reason, the expected total revenue is calculated by the sum of total sales of imperfect and perfect items. The perfect-quality items are sold at a price of per unit of weight, while at the end of the inspection period , the imperfect-quality items are vended as a single batch with a less price of per unit of weight. As a result, the expected total revenue per period is given by

4.2. Total Cost per Period

The expected total cost per period includes the purchasing cost, setup cost, feeding cost, expected holding cost, inspection cost, backordering cost, and expected carbon emissions cost. One has

The objective of this inventory model is to maximize the expected total profit which is obtained by subtracting the expected total cost per period from the expected total revenue per period. Then,

The following section provides a detailed discussion of the calculations of the aforesaid costs.

4.3. Purchasing Cost per Period

At the beginning of each period, the company buys newborn items at a cost of per unit of weight; each one with a weight of . Thus, the purchasing cost per period is obtained with

4.4. Setup Cost per Period

At the commencement of each period, a setup cost of is carried out. So, the setup cost per period is determined as

4.5. Feeding Cost per Period

The growing items are fed during , and a feeding cost is incurred by the company at per unit of weight. As the growing items become older and bigger, they need more food. For that reason, the quantity of food consumed by the items depends on the age (weight) of the items according to the growth function . The feeding cost per period is calculated with

4.6. Expected Holding Cost per Period

The company incurs in a holding cost for maintaining in good conditions the items in storage during the cycle time . The expected holding cost is determined by multiplying the holding cost per unit of weight per unit of time by the inventory carried out during . The whole inventory held is calculated as the sum of the areas (see Figure 5 in Appendix A). So, the expected holding cost is given by

The detailed derivation of the expected holding cost is given in Appendix A.

4.7. Inspection Cost per Period

Throughout , a 100% inspection process is performed at a rate of with the aim of splitting the perfect-quality items from imperfect-quality items. The company incurs in an inspection cost of per unit of weight. In consequence, the inspection cost per period is expressed as

4.8. Shortage Cost per Period

During the course of , the shortages are accumulated till these reach the backordering quantity of units of weight. The company has a backordering cost for the management of the accumulation of shortages. The backordering cost is computed by multiplying shortage cost per unit of weight per unit of time by area (see Figure 5 in Appendix A). Therefore, the backordering cost is calculated as

The detailed derivation of the backordering cost is specified in Appendix B.

There are several research works that suppose that the carbon emissions generated by the companies are due to the transportation and warehousing activities. But, in fact, there are a lot of operations that emit carbon emissions. For example, the operations involved in the following process also cause carbon emissions: purchasing, setup, feeding, holding inventory, and inspection, among others. It is important to mention that some growing items per se generate carbon emissions (cows, goats, pigs, etc.). The carbon emissions caused by these processes are determined as follows:

4.9. Carbon Emissions Caused by the Purchasing Action

4.10. Carbon Emissions Produced by the Setup Activity

4.11. Carbon Emissions Generated during the Feeding Process

4.12. Carbon Emissions Created in Holding Inventory’s Operations

4.13. Carbon Emissions Made by the Inspection Process

Carbon tax is one of the well-known mechanisms imposed by government regulations as a penalty. This means that the companies need to pay a tax on the amount of carbon emissions. As a result, the carbon emissions cost per period that the company must pay is

4.14. Expected Total Profit Function

The expected total profit per period is calculated in the following manner: the expected total revenue per period minus the expected total cost per period. Basically, the expected total profit is formulated by substituting equations (9) and (12) to (17) and equation (23) into equation (11). Then, the expected total profit per unit time per period is determined as . Thus,

Considering that the demand rate depends on selling price with a polynomial function , the expected total profit per unit of time becomes

Then, the objective is to maximize the expected total profit per unit of time. Therefore, the optimization problem is formulated as follows:

It is important to remark that the abovementioned maximization formulation is a nonlinear optimization problem.

5. Solution Procedure

Firstly, the unconstrained optimization problem is considered in order to obtain some theoretical results. Secondly, an algorithm for finding the optimal solution to the constrained problem given by equation (26) is developed.

5.1. Theoretical Results

The aim is to find the optimal values for order quantity , backordering quantity , and selling price that maximize the expected total profit per unit of time. As the expected total profit per unit of time function is continuous and twice differentiable with respect to the three decision variables on the interval , there exists a global maximum on that interval.

The necessary conditions which must be satisfied for the optimality of solution that maximizes the expected total profit per unit of time function are , , and .

The first partial derivative of given in equation (25) with respect to is expressed as follows:

The first partial derivative of with respect to is given by

The first partial derivative of with respect to is

By solving equations (27) and (28) simultaneously for the decision variables and , one has

The sufficient conditions for the optimality of the solution that maximize are given in Appendix C.

5.2. Algorithm for Finding the Optimal Solution

The following Algorithm 1 is proposed taking into account the theoretical results and constrains.

Step 1. Provide the input parameters of the inventory system.
Step 2. Compute the selling price , the order quantity , and the backordering quantity by solving simultaneously equations (29)–(31).
Step 3. If the optimality conditions are satisfied, then go to Step 4. Else, go to Step 8.
Step 4. If , then go to Step 7. Else, go to Step 5.
Step 5. If , then set , determine the order quantity with equation (30) and the backordering quantity with equation (31), and go to Step 6. Else, set and calculate the order quantity with equation (30) and the backordering quantity with equation (31) and go to Step 6.
Step 6. Calculate the expected total profit per unit of time with equation (25).
Step 7. Report the solution: and .
Step 8. Stop.
5.3. Special Cases

The proposed inventory model developed in this research work is a general inventory model because it contains the following inventory models as special cases: Sebatjane and Adetunji’s [40] improved inventory model, Maddah and Jaber [11], Chang and Ho [13], Shih [61], Silver [62], and the traditional EOQ inventory models with and without shortages.(i)When , it means that the product has a linear price-dependent demand . It is expressed as follows:(ii)When shortages and carbon tax rate are not considered and the demand is not dependent on price (i.e., and ), then an improved version of Sebatjane and Adetunji’s [40] inventory model is obtained. It is shown as follows:(iii)When shortages are permitted, carbon tax rate is not involved, the demand is not dependent on price, and the products are not of the type of growing items, and thus, the feeding process is not incorporated (i.e., , and ), and then, the inventory model of Chang and Ho [13] is determined. It is given by(iv)When shortages are not allowed, carbon tax rate is not taken account, the demand is not dependent on price, and the products are not of the type of growing items (, and ); then, Maddah and Jaber’s [11] inventory model is derived. It is presented as follows:(v)When shortages are not tolerable, carbon tax rate is not considered, the demand does not depend on price, the goods are not of the type of growing items, and inspection rate is sufficiently large (i.e., and ); then, it is converted to the inventory model of Shih [61] and Silver [62]. It is expressed as(vi)When shortages are permissible, carbon tax rate is not incorporated, the demand is not dependent on price, the products are not of the type of growing items, and the products are of perfect quality (i.e., , and ); then, the traditional inventory model with shortages and with full backordering is found. It is written as follows:(vii)When shortages are not allowable, carbon tax rate is not taken into consideration, the demand does not depend on price, the goods are not of type of growing items, and these are of good quality (i.e., , and ); then, it is transformed to the traditional inventory model without shortages proposed by Harris [1]. It is

It is important to remark that the inventory model developed in this research does not reduce to Salameh and Jaber’s [8] and Wee et al.’s [10] inventory models because those inventory models contain a flaw that was identified and corrected by Maddah and Jaber [11] and Chang and Ho [13], respectively. Basically, Maddah and Jaber [11] and Chang and Ho [13] proposed improved versions for Salameh and Jaber’s [8] and Wee et al.’s [10] inventory models, respectively. For this reason, the proposed inventory model developed in this research work, in fact, converges to Maddah and Jaber’s [11] and Chang and Ho’s [13] inventory models as it was mentioned above.

6. Numerical Examples

This section presents and solves some numerical examples in order to illustrate the applicability of the proposed inventory model.

Example 1. This example considers the input parameters of Sebatjane and Adetunji [40]. To solve this inventory problem, it is stated , and . The values for the data are , , , ZAR/g, ZAR/g, ZAR/cycle, ZAR/g/year, ZAR/g, ZAR/g, ZAR/g, g/year, g, /year, g, and g, and the percentage of imperfect growing items follows a uniform distribution with the probability density function which is given as follows:Considering ,Then, and are computed as follows:As it was mentioned before, the proposed inventory model generates an improved version of Sebatjane and Adetunji’s [40] inventory model. Additionally, the proposed inventory model also can be used for optimizing the order quantity and backordering quantity when the selling price is given. Obviously, it also optimizes the three decision variables . Therefore, the numerical example is solved for different values of for the four inventory models: (I) Sebatjane and Adetunji [40] (original version), (II) Sebatjane and Adetunji [40] (improved version), (III) the proposed inventory model when selling price is given, and (IV) the proposed inventory model when the selling price is optimized. Table 2 shows the comparison of the optimal solutions among the four inventory models.
From Table 2, it is observed that the expected total profit of both the proposed inventory models (see columns III and IV) is greater than that of Sebatjane and Adetunji’s [40] inventory models (see columns I and II).

Example 2. This example uses the input parameters of Wee et al. [10] and Chang and Ho [13]. To solve this inventory system, the following is established: , and . The data are as follows: , , , , $/unit, $/unit, $/cycle, $/unit/year, $/unit/year $/unit, $/unit, units/year, and . The optimal solution is as follows: $/year, units, and units. This solution is the same as in the work of Chang and Ho [13].

Example 3. Now, let us consider some dataset of Sebatjane and Adetunji [40] which is ZAR/g, ZAR/cycle, ZAR/g, ZAR/g, ZAR/g, g/year, g, , /year, g, g, and . Here, the holding cost is ZAR/g/year. For the implementation of the proposed inventory model, additional information is required. These data are related to the type of the demand function dependent on the price, the backordering cost, and carbon emissions. Therefore, it is necessary to state the following data: Assume that the demand rate of the growing items follows a polynomial function given by with values of , , and . The backordering cost is 0.1 ZAR/g/year. The relevant input parameters related to carbon emissions are 0.0045 ZAR/tons, 2000 tons/year, 0.2 tons/year, 0.65 tons/year, 0.375 tons/year, and 0.005 tons/year. By applying the proposed algorithm, the optimal solution is obtained: ZAR/year, units of newborn growing items, g, and ZAR/g. Figures 24 demonstrate graphically the concavity property of the expected total profit with respect to pairs of decision variables and ; and ; and and , respectively.

7. Sensitivity Analysis

This section provides a sensitivity analysis in order to investigate the effects of changing the input parameters on the expected total profit per unit of time , order quantity , backordering quantity , and selling price of the growing items. The sensitivity analysis is performed taking into consideration the data of Example 3 when only one parameter changes at a time and other parameters are kept at their original values. Specifically, the sensitivity analysis is performed to study the effect of scale, sensitivity and power index demand parameters, selling price and expected value of a percent of poorer-quality growing items, setup cost, holding cost, backordering cost, feeding cost, purchasing cost, inspection cost, carbon tax rate, and the amount of carbon emissions caused by the operations of the processes in ordering, holding, feeding, purchasing, and screening. The results of the sensitivity analysis are given in Tables 35. Based on the behavioral changes as reflected in Tables 35, the following observations are made:(1)The value of is highly sensitive to the demand parameters and less sensitive to other parameters. On the one hand, the higher the value of the scale parameter of demand , the higher the value of due to the fact that demand increases; therefore, the sales increase, and this leads to high profits. For this reason, it is suggested for the manager to apply advertising actions in order to boost the demand. On the other hand, when the sensitivity parameter of the demand grows, the declines because this parameter has a negative impact on the demand, making it drop. For this reason, it is advisable for the decision maker to reduce the selling price. However, a higher value of induces higher values of but a lower value for . This means that as the setup cost increases, the company must raise selling price .(2)The value of is more sensitive to the parameters and less sensitive to other parameters. The higher the value of the parameter , the smaller the value of .(3)The value of is more sensitive to the parameters and less sensitive to other parameters. The higher the value of , the smaller the value of .(4)The value of is more sensitive to the parameters and less sensitive to other parameters. The higher the value of the parameter , the higher the value of . This means that as demand scale parameter increases, the company must raise selling price which will also directly impact positively the expected total profit. The higher the value of the parameters , the smaller the value of .(5)When the expected value of the defective growing items decreases, the expected total profit also decreases .(6)Changes on carbon emission parameters have a regular influence on the expected total profit .

With the information of Tables 35 is constructed Table 6. Table 6 summarizes the sensitivity analysis study in a visual manner which is more helpful for the decision makers due to that these can more easily observe how changes in parameters affect to the decisions variables and the expected total profit per unit of time in order to take the best decisions based on numerical facts.

8. Conclusions and Future Research

This paper develops an inventory model for growing items with imperfect quality and carbon emissions in which the demand rate is price sensitive according to a polynomial function. Furthermore, the shortages are permissible and these are fully backordered. To optimize the expected total profit per unit of time, some theoretical results are derived, and with these results, an effective and efficient solution procedure is developed. It is important to state that the algorithm generates the optimal solution to order quantity, backordering quantity, and selling price which maximizes the expected total profit per unit of time. The effects of carbon tax and the amount of carbon emissions are studied. Some important managerial insights are obtained from numerical examples and sensitivity analysis. The results show that carbon tax has a regular role on the reduction of carbon emissions. It is important to pay attention to quality and ensure that the percentage of defective items is kept to a minimum. Furthermore, if the setup cost increases, the company must raise selling price , and when the demand elasticity parameter increases, the company must raise selling price which also directly impacts the total profit. Furthermore, the proposed inventory model is a generalized inventory model due the fact that several previously published inventory models are particular cases. It was found that the policy that permits shortages with full backordering is more economical that the one that avoids shortages. Moreover, it was also found that the proposed inventory model outperforms the inventory model of Sebatjane and Adetunji [40].

There are several possible extensions of the proposed inventory model that can be explored. These extensions constitute future research endeavors in the inventory management of imperfect growing items. For instance, consider to investigate the effect that only a percentage of defective items can be reworked, and the others must be eliminated immediately. Other research studies that can be conducted are to include new aspects such as stock-dependent demand, nonlinear holding cost, vendor-managed inventory (VMI) with consignment stock (CS), inflation, volume discounts, deterioration, trade credit, supply chain environment, and a vendor-buyer inventory model with multiple shipments, partial backordering, advertising, and multiple products subject to constrains such as space, budget, and time. These are, among others, some interesting and challenging subjects of ongoing future investigation that academicians and researchers would like to study in the future.

Appendix

A. Determination of the Expected Holding Cost

The holding cost is calculated by the multiplication of the unit holding cost and the total inventory accumulated during , and this is the area above level zero in Figure 5. This area is computed as the sum of the areas .

Considering Figure 5, the five areas are defined as follows.

The area is obtained as follows:

The area is found as follows:

The area is determined in the following manner:

The area is defined as follows:

The area is given by

The total inventory accumulated during is computed as follows:

Substituting the corresponding areas,

Simplifying,

Taking the expected value of the total inventory accumulated , the following expression is obtained:

Rearranging terms,

Finally, the expected holding cost is given by

B. Determination of the Backordering Cost

The backordering cost is determined by the multiplication of the unit backordering cost and the total shortages accumulated during , and this is the area shown below level zero in Figure 5. This area is obtained as follows:

Thus, the backordering cost () is given by

C. Sufficient Conditions for the Optimality

For the sake of brevity, only the final equations are provided for the direct and cross second-order partial derivatives.

The direct second-order partial derivatives are given as follows.

The second-order partial derivative of equation (25) with respect to is as follows:

The second-order partial derivative of equation (25) with respect to is expressed as

The second-order partial derivative of equation (25) with respect to is given as follows:

The second-order cross partial derivatives are given as follows.

The second-order cross partial derivative of equation (25) is presented as follows:

The second-order cross partial derivative of equation (25) is written in the following manner:

The second-order cross partial derivative of equation (25) is shown in the following way:

The second-order cross partial derivative of equation (25) is given as follows:

The second-order cross partial derivative of equation (25) is given as follows:

The second-order cross partial derivative of equation (25) is

Optimality is given as follows.

To prove the optimality of three decision variables, it is required to construct the Hessian matrix , which is given as follows:

The sufficient conditions for optimality through the Hessian matrix are as follows:

, and , where , , and are the following matrixes:

Thus,

The abovementioned analysis proves that the function is strictly concave and shows that the Hessian is negative-definite. Therefore, the optimal solution for the decision variables exists and maximizes the expected total profit.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research was supported by the Tecnológico de Monterrey Research Group in Optimization and Data Science (0822B01006).