Next Article in Journal
Assessing CO2 Emissions of Online Food Delivery in Japan
Next Article in Special Issue
Research on Sustainable Cooperation Strategies for Cross-Regional Supply Chain Enterprises in Uncertain Environments
Previous Article in Journal
Cenozoic Sedimentary Characteristics of the East China Sea Shelf Basin and an Evaluation of the Suitability of Geological Storage of Carbon Dioxide in the Saline Water Layer
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Economic Growing Quantity Model with Mortality in Newborn Items and Inhibition Cost of Ammonia Production under All-Units Discount Policy

by
Amir Hossein Nobil
,
Erfan Nobil
,
Leopoldo Eduardo Cárdenas-Barrón
*,
Dagoberto Garza-Núñez
,
Gerardo Treviño-Garza
,
Armando Céspedes-Mota
,
Imelda de Jesús Loera-Hernández
and
Neale R. Smith
Tecnológico de Monterrey, School of Engineering and Sciences, E. Garza Sada 2501 Sur, Monterrey 64849, Mexico
*
Author to whom correspondence should be addressed.
Sustainability 2023, 15(10), 8086; https://doi.org/10.3390/su15108086
Submission received: 27 January 2023 / Revised: 19 April 2023 / Accepted: 10 May 2023 / Published: 16 May 2023
(This article belongs to the Special Issue Sustainable Supply Chain and Inventory Management)

Abstract

:
The amount of global ammonia (NH3) emissions is growing continuously, similar to the damage to the environment, particularly humans and animals, caused by those emissions. Various problems derived from pollution by ammonia emissions have attracted increasing attention in recent years. In particular, accumulation of ammonia in poultry farms is a concern for the poultry industry as it can lead to possible damage due to reduced bird performance, damage to the respiratory tract and skin of birds, and thus loss of customers. As birds age, ammonia production increases due to factors such as feeding and mobility, requiring the application of solutions to reduce it such as the use of fans, feed supplements, and temperature adjustments to improve bird health. These solutions impose additional costs on poultry farms to combat ammonia emissions. This study presents a general economic growing quantity (EGQ) model that includes the cost of inhibition of ammonia production during the growing period. In addition, the model is formulated under an all-units discount policy, where the price of newborn items is related to the size of the order purchased from the supplier. Furthermore, the model assumed that some newborn items are dead when the lot is received because of stress experiences and incidents during the catching, loading, transportation, and unloading. Finally, two versions of the proposed general EGQ model are presented: EGQ with no discount and EGQ with known slaughter age.

1. Introduction

The analysis of inventory problems for growing items is a new and significant topic in the inventory management literature. The model of economic order quantity (EOQ) for growing items or economic growing quantity (EGQ) was first presented by Rezaei [1]. EGQ models do not consider the common assumption of lot size problems, in which the weight of the items remains constant until consumption, while the former models assume that the weight increases after the items are received. EGQ models are common and usable in inventory systems such as livestock, poultry, and fish farming. In these cases, a company purchases several newborn animals, such as chicks, calves, or lambs from suppliers, and upon receipt, the farm raises them until they reach proper slaughter weight. The EGQ model balances the amount of the initial order for newborn items and the age of slaughter to maximize profit. The model considers the items’ growth function that shows the relationship between their weight and age. Therefore, the weight of the items is considered, instead of the number of items, in the EGQ model inventory graph.
Zhang et al. [2] presented a second paper on lot sizing problems for growing items in 2016. They considered one of the important aspects of environmental effects for companies, as they modeled a government-imposed carbon production tax over the growth period of the items. Therefore, they added this tax to the total inventory cost and obtained the optimal order quantity and the optimal slaughter age. Carbon emissions that occur from animal meat production are primarily associated with manure emissions and the production of feed for livestock and poultry. Due to these factors, the government sets a tax that emitters must pay for carbon production. Next, De-la-Cruz-Márquez et al. [3] developed an EGQ model with shortage and price-dependent demand, based on the model presented by Zhang et al. [2]. Additionally, Rana et al. [4] proposed an EGQ model for the effects of carbon emissions with a deteriorating process under a permissible delay policy. Concurrently, Gharaei and Almehdawe [5] addressed the impacts of carbon emissions emitted from produced manure, fermentation processes, and transportation in an EGQ model.
Nobil et al. [6] developed an EGQ model considering the admissibility of deficiency. An innovation in their work is the use of a linear approximation function for the growth function, which made the optimal solution easily derived in a few simple steps. After that, Nobil and Taleizadeh [7] presented a practical and realistic mode for mathematical modeling of the EGQ problems. When ordering newborn animals, it must be ensured that the quantity ordered is discrete; therefore, in their model, the corresponding decision variable is modeled as a discrete number.
A common assumption in the EGQ models has been to consider an inventory system for raising only one growing type of animal, until Khalilpourazari and Pasandideh [8] presented an EGQ model for an inventory system with multiple growing items. They assumed that several types of newborn animals could be bought and bred. Additionally, they integrated other assumptions such as storage capacity, budget in hand, and the maximum allowed cost of breeding and maintenance. Finally, the model, due to its structural complexity, was solved with two new meta-heuristic algorithms named the sine cosine crow search algorithm and the water cycle moth–flame optimization algorithm. The vast majority of EGQ models have focused on the farmer level, but Malekitabar et al. [9] considered two levels: farmer and supplier. The supplier raises newborn rainbow trout and then sells them to the farmer. Subsequently, the fish farmer grows the rainbow trout until they reach maturity and slaughter age, finally selling them to the market.
Recently, Gharaei and Almehdawe [10] applied the uniform distribution as the probability density function of survival and mortality over the growth period. Likewise, Sebatjane and Adetunji [11] proposed an EGQ model with the mortality function as a fraction of total items that do not survive the growing period. Another important extension to the EGQ models is the number of poor-quality items produced during the growing period that should be identified and removed. Sebatjane and Adetunji [12] extended the models of Rezaei [1] and Nobil et al. [6]. They considered an inspection time at the end of the breeding period, in which poor-quality items are identified and sold to the market at a lower price. In accordance with a sustainable green breeding policy, Nobil et al. [13] suggested an EGQ model with a practical polynomial function of CO2 emissions, which is dependent both on the age of the animals and on the mortality function. The goal was to choose the best number of chicks to purchase from the supplier at the correct age for slaughter to save on overall costs. To increase the profits of livestock and poultry farms and reduce overall carbon emissions, Khan et al. [14] examined combined prepayment installment, pricing, and replenishment choices for a growing item under cap-and-price, cap-and-trade, and carbon tax legislation.
Under inventory policies, suppliers commonly recommend one of two forms of quantity discounts. The all-units quantity discount results in decreased procuring costs for the whole order if the quantity requested is above a precise quantity named the breakpoint. The incremental-units quantity discount results in decreased procuring cost that only applies to products purchased above the breakpoint. We found only one study in growing items that works with discount policies. Sebatjane and Adetunji [15] applied the incremental discount policy for the purchase of growing items, considering a price that depends on the order size.
The inclusion of the food supply chain echelons is another significant extension in the literature on growing items. Sebatjane and Adetunji [16] addressed an inventory-production problem for growing items in a supply chain with three levels: farming, processing, and retail. The farmer raises newborn animals and then sends them to a processor once the animals reach the slaughter age. The processing plant slaughters the animals and packages them at a specified rate. The processor sends a certain number of equally sized shipments of processed items to a retailer that meets the consumer’s demand. Sebatjane and Adetunji [17] presented an EGQ model for optimizing management decisions in a perishable food items supply chain that starts with farming operations in which live inventory products are raised and ends with the consumption of processed stock. The farming and consumption stages are linked by a processing stage during which livestock or poultry are processed into a consumable product. Consumer demand at the consumption stage is a function of the selling price and the freshness of the processed stock. The farming, processing, and consumption stages comprise the three-tier network model with the objective to maximize the joint supply chain profit. Later, Sebatjane and Adetunji [18] addressed an EGQ model for a three-level supply chain, including a farmer, a processor, and a retailer. They assumed that customer demand is dependent on the inventory level and the product expiration date. In addition, Pourmohammad-Zia et al. [19] evaluated the effects of marketing policies and perishable items in place for a two-echelon supply chain, including a supplier and a retailer. The model focuses on two scenarios: centralized and decentralized supply chain with a profit-sharing contract.
Pourmohammad-Zia et al. [20] also extended an EOQ model for growing items in a three-level supply chain (supplier, producer, and multiple retailers) by establishing a trade-off between cost-efficiency and market coverage. Lately, Sebatjane [21] suggested some inventory management techniques for a multi-level supply chain for perishable foods that includes components that are deteriorating and growing. His study assumed that the farming echelon, where young growing items are bred until they reach an age of maturity, is at the upstream end of the supply chain. The live inventory products are delivered to the processing echelon, which includes activities such as slaughtering, cutting, and packaging.
A summary of the major contributions made by scholars on these topics is presented in Table 1. According to our search, the effects of ammonia on an inventory system have not been considered in any of the past studies on inventory models. Wang et al. [22] reported that the emission of ammonia, which is a hazardous gas that exists in industry and agriculture, has gradually been identified as a major participant in the high levels of haze in many developing and developed countries. Based on the study of Koerkamp [23], absorption of surplus amino acids is mainly excreted as uric acid in poultry, while undigested protein is excreted through the feces. Microbial degradation of undigested protein and uric acid in broiler litters creates ammonia (NH3) as a by-product. Additionally, excreted nitrogen (N) is associated with deteriorating litter quality, which exists in NH3 content and contributes to the spread of hock and foot lesions in broilers. High NH3 concentrations inside the broiler farms, and prolonged exposure to this volatile N compound, negatively impact the performance and health of birds (Miles et al. [24]), as well as the health of farmworkers and customers (Webb et al. [25]).
Wang et al. [26] stated that the deposition of ammonia in the environment can impact plant diversity and productivity, as well as the life cycle of aquatics. Additionally, based on the work of Brink et al. [27], NH3 is a precursor for particulate matter that contributes to lower air quality. Strategies such as nutrition adaptations and ventilation are used in poultry farms to reduce NH3 emissions. Although these strategies avoid emissions, they represent an extra cost for the farms. Therefore, this study presents a general EGQ model with an ammonia production inhibition cost during the breeding period.
Another innovative aspect of this study is the use of a discount policy when purchasing items from suppliers. None of the previous studies, except Sebatjane and Adetunji [15], considered the effects of discounts on inventory control decisions for growing items. Their model investigates the incremental discount policy on an EGQ model for a livestock farm. The proposed model considers the impact of an all-units discount policy on the EGQ problem and investigates the mortality for newborn items purchased from the supplier. Additionally, Yılmaz et al. [28] investigated that stress-causing incidents and practices during the catching, loading, transportation, and unloading may increase the mortality of birds. Hence, this study assumes that a percentage of the newborn items received from the supplier are dead. In summary, a general EGQ model with mortality in newborn items and an inhibition cost of NH3 production under an all-units discount policy is presented. In addition, two specific cases are considered for the general model: EGQ with no discount and EGQ with known slaughter age.
Table 1. A summary of the major contributions made by some scholars in the literature of EGQs.
Table 1. A summary of the major contributions made by some scholars in the literature of EGQs.
PaperYearObjective FunctionSolution MethodGrowth CurveDead Newborn AnimalsNH3 EmissionDiscountTypes
of Items
Rezaei [1]2014PFAnalyticalBW P
Zhang et al. [2]2016CFAnalyticalBW P
Sebatjane and Adetunji [12]2019PFAnalyticalBW/LR P
Sebatjane and Adetunji [15]2019CFAnalyticalBW IDL
Sebatjane and Adetunji [16]2019CFAnalyticalBW L
Khalilpourazari and Pasandideh [8]2019PFMeta-heuristicBW P
Malekitabar et al. [9]2019PFGame theoryBW* F
Nobil et al. [6]2019CFAnalyticalLR P
Nobil and Taleizadeh [7]2019CFClosed formLR P
Sebatjane and Adetunji [17]2020PFAnalyticalBW L
Gharaei and Almehdawe [10]2020CFClosed formLR P
De-la-Cruz-Márquez et al. [29]2021PFAnalyticalBW P
Gharaei and Almehdawe [5]2021CFMeta-heuristicBW P
Pourmohammad-Zia et al. [19]2021PFAnalyticalBW P
Pourmohammad-Zia et al. [20]2021PFGame theoryBW P
Rana et al. [4]2021CFAnalyticalBW P
Sebatjane and Adetunji [19]2021PFAnalyticalBW P
De-la-Cruz-Márquez et al. [3]2022PFAnalyticalBW P
Sebatjane and Adetunji [11]2022CFAnalyticalBW P
Sebatjane [21]2022CFAnalyticalBW P
Nobil et al. [13]2023CFAnalyticalBW P
Khan et al. [14]2023PFAnalyticalBW P
This paper CFAnalyticalBW**ADP
PF: profit function, CF: cost function, BW: biological weight, LR: linear rate, AD: all-units, ID: incremental, P: poultry, L: livestock, F: fish. The “*” means that the research work includes the characteristic.
The rest of this paper is organized as follows: The three proposed models are explained and formulated in Section 2. Next, Section 3 presents the solution procedures for each model, whereas Section 4 uses some numerical examples to validate the proposed algorithms. Finally, Section 5 presents some general conclusions and future research ideas.

2. Proposed Models

This section documents the proposed models, the notations, and the assumptions. Then, the mathematical formulation of the models (a) with no discount, (b) all-units discount, and (c) all-units discount with known slaughter age are explained by considering the nature of the growing items.

2.1. Notations and Assumptions

The indices of the time and discount categories based on a price are defined as t = 1,2 , and i = 1,2 , , J , respectively, and the rest of notation of the model is as follows:
  • Parameters:
J : the number of discount categories based on a price; i = 1,2 , , J
m t : the weight of a broiler at time t (weight unit)
E : the highest length of breeding time (day)
W : the lowest length of the breeding time (day)
D : the rate of demand to weight unit (weight unit/year)
α : the fraction of the purchased chicks from the supplier that are dead or of poor quality.
r : the holding cost percentage rate per year.
u i : the weight volume breakpoint in units for a unit purchasing cost (weight unit)
p i : the purchasing cost per weight unit for price range i (monetary unit/weight unit) if u i y < u i + 1
S : the setup cost per growing cycle (monetary unit/setup)
x : the cost of production (feeding) during growing period per unit item (monetary unit/unit item)
z : the inhibition cost of ammonia production (monetary unit/microliter/day·square meter·weight unit)
f t : the polynomial function of feeding consumption (unit items)
v t : the exponential function of ammonia production (microliter/day·square meter·bird).
k : the growth rate.
A : the asymptotic weight.
b : the integration constant of the growth function.
n : the shape parameter of the growth function.
t d : the demand periods.
Q : the total weight at the slaughter time (weight unit).
T C U : the total cost per year (monetary unit).
T C U i : the total cost of breakpoint i per year (monetary unit).
  • Decision variables:
y : the total number of growth items ordered at the beginning of each cycle (unit items)
t : the slaughter age (day)
Figure 1 illustrates the weight stock level of broilers in the proposed model. In this inventory system, y newborn chicks with an initial weight m 0 are purchased from the supplier. When the chicks are transported from the supplier to the poultry farm, a percentage of them (α) arrive dead. As a result, the number of live chicks that are available at the beginning of the period is equal to 1 α y . The poultry farm raises them for a period until they reach the ideal slaughter weight ( m t ). At the time of slaughter ( t ), all broilers begin processed and shipped to the market for consumption, during the consumption period ( t d ).
There are several works on live bird growth that define different growth functions. However, Richards in 1959 suggested one of the most important growth functions for live bird weight, which is named Richards’ function (Rezaei, [1]). This function is as follows:
m t = A 1 + b e k t 1 n
where m t is the body weight of a live bird at age t .
In the breeding period, birds grow and feed with a function that is proposed by Goliomytis et al. [30]. Their age is represented by a polynomial function. The function is as follows:
f t = a 0 + a 1 t + a 2 t 2 + a 3 t 3
where a 0 , a 1 , a 2 , and a 3 are fixed and known coefficients.
According to Shah et al. [31], the biological process of broiler litter, manure, and protein composition has an important outcome, NH3, emissions in poultry farms. NH3 has a harmful effect that jeopardizes the health of live birds. The ammonia production function based on the age of the birds is taken from the work of Leonard et al. [32]. This function is an exponential curve and is defined as:
v t = θ 0 e θ 1 t
where θ 0 and θ 1 are fixed and known coefficients.

2.2. EGQ with No Discount

Due to customer demand, the model considers that a poultry farm purchases newborn birds (raw material), grows them to meet market slaughter weight, and processes them. The farm needs to cover the total cost of the inventory system, which includes setup cost, purchasing cost, holding cost, feeding cost, and ammonia production inhibition cost.
The cycle begins with cleaning activities to start the breeding process. According to Figure 1, when time is zero, the number of newborn chicks ( y ) is received with initial weight ( m 0 ) from the supplier, so the purchasing cost per growth cycle is equal to p y m 0 , where p is the purchasing cost per weight unit without discount. Additionally, the length of the consumption period is t d , while the average weight of the inventory during the consumption period is 1 α y m t / 2 . The holding cost per cycle is r p t d y m t / 2 , where r p is the annual inventory cost per weight unit with no discount. The feeding function per weight unit, described in Equation (2), depends on the age of the bird. It is noted that x and t are considered as the feeding cost per unit and the length of the breeding period, respectively. The total feeding cost per cycle is x 1 α y 0 t f t d t = x 1 α y 0 t ( a 0 + a 1 t + a 2 t 2 + a 3 t 3 ) d t .
Finally, there is a cost of dealing with ammonia production defined by z 1 α y 0 t b 0 e b 1 t d t . Therefore, the total cost per cycle with no discount is as follows:
T C = S + p y m 0 + r p t d 1 α y m t 2                                                                                                             + x 1 α y 0 t ( a 0 + a 1 t + a 2 t 2 + a 3 t 3 ) d t                                             + z 1 α y 0 t θ 0 e θ 1 t d t
To meet the total cost function per year, T C should be divided by t d = 1 α y m t / D = 1 α y A 1 + b e k t 1 n / D . This means that the total cost function per year without discount is:
T C U = D S 1 α y A 1 + b e k t 1 n + D p m 0 1 α A 1 + b e k t 1 n                                                                             + r p 1 α y A 1 + b e k t 1 n 2 + D z θ 0 e θ 1 t 1 θ 1 A 1 + b e k t 1 n             + D x a 0 t + a 1 t 2 2 + a 2 t 3 3 + a 3 t 4 4 A 1 + b e k t 1 n
where y and t are the decision variables.

2.3. General EGQ

An EGQ model is considered in which the supplier offers an order quantity based on purchase cost per weight unit as given by:
p i = p 1 u 1 y < u 2 p 2 u 2 y < u 3                           p J u J y
where y is the order quantity, u 1 < u 2 < < u J and p 1 > p 2 > > p J .
The total cost with an all-units discount is approximately equal to the case with no discount with only one difference: the purchasing and holding costs are changed with respect to the number of newborn birds that are ordered, as the acquisition price depends on the ordered quantity. Hence, the total cost function per year with an all-units discount ( T C U i ) is defined as follows:
T C U i = D S 1 α y A 1 + b e k t 1 n + D p i m 0 1 α A 1 + b e k t 1 n                                                         + r p i 1 α y A 1 + b e k t 1 n 2 + D z 0 t θ 0 e θ 1 t d t A 1 + b e k t 1 n                                                         + D x 0 t a 0 + a 1 t + a 2 t 2 + a 3 t 3 d t A 1 + b e k t 1 n ; u i y < u i + 1
where y and t are the decision variables and the price ( p i ) depends on the order size of y .

2.4. EGQ with Known Slaughter Age

This model assumes that the slaughter weight of chickens is known due to market behavior with respect to consumption weight through purchasing data. For example, in Asian countries such as Iran, this weight is 1.5 kg, in European countries such as the United Kingdom is about 2.2 kg and in the United States is about 3 kg (Nobil and Taleizadeh, [6]). As a result, once the consumption weight of broilers is known, the slaughter age is obtained according to the biological weight function Equation (1). Thus, slaughter age ( t ) can be considered as an input parameter and not as a decision variable. Therefore, the birds’ weight is easily computed using the slaughter age using Equation (1). In addition, the values of f ( t ) , and v ( t ) are obtained by defining the slaughter age from Equations (2) and (3), respectively. Based on the value of t , the total cost in Equation (7) is adapted as follows:
T C U i = D S 1 α y m t + r p i 1 α y m t 2 + D p i m 0 1 α m t + D z v t m t + D x f t m t ; u i y < u i + 1
where y is the only decision variable, and p i is a dependent variable.

3. Solution Procedures

In this section, the solution procedures of the three proposed models are shown. Due to the nature and assumptions of the models, each one has a specific procedure to obtain their optimal solution.

3.1. EGQ with No Discount

The decision variables used in this model are the slaughter age and the order size. Clearly, based on objective function (5), the purchasing cost is an input parameter with a constant value and does not depend on the order size. As the model is convex (see Appendix A), then the partial derivative of objective function (5) with respect to t and y are set equal to zero, as follows:
T C U t = D S k b e k t A n y α 1 1 + b e k t 1 n + 1                                         D p m 0 k b e k t A n α 1 1 + b e k t 1 n + 1                                                                   + A k p r y b e k t α 1 1 + b e k t 1 n 1 2 n                                                                                             + D z θ 0 e θ 1 t A 1 + b e k t 1 n + D z b k θ 0 e k t e θ 1 t 1 A n θ 1 1 + b e k t 1 n + 1                                 + D x a 0 + a 1 t + a 2 t 2 + a 3 t 3 A 1 + b e k t 1 n                                                                           + D x b k e k t a 0 t + a 1 t 2 2 + a 2 t 3 3 + a 3 t 4 4 A n 1 + b e k t 1 n + 1 = 0
T C U y = D S 1 α y 2 A 1 + b e k t 1 n + r p 1 α A 1 + b e k t 1 n 2 = 0
Based on Equation (10), the value of the order size is equal to:
y = 2 D S r p 1 α 2 A 2 1 + b e k t 2 n
The optimal slaughter time cannot be obtained by a closed-form method from Equation (9), so a numerical method is used to solve it. Therefore, by substituting the order size value from Equation (11) into Equation (9), then the optimal slaughter age is obtained using the Newton–Raphson method by applying the following solution procedure:
  • Step 1. Substituting the order size value of Equation (11) into Equation (9), solve for the value of t from Equation (9) using the Newton–Raphson method.
  • Step 2. Based on the value of t * , compute the optimal value of y * using Equation (11).
  • Step 3. Based on the value of t * and y * , determine the optimal value of T C U * using Equation (5).

3.2. General EGQ

In this model, in addition to considering the slaughter age and the order size as decision variables, the purchasing cost is considered as not constant and as a dependent variable of the order size. This model is an extension of the EGQ model without discount. Like in the previous solution procedure, the partial derivatives of objective function (7) are set equal to zero as follows (see Appendix B to review the procedure to calculate the number of orders):
T C U i t = D S k b e k t A n y α 1 1 + b e k t 1 n + 1                                             D p i m 0 k b e k t A n α 1 1 + b e k t 1 n + 1                                                                               + A k p i r y b e k t α 1 1 + b e k t 1 n 1 2 n                                                                                                     + D z θ 0 e θ 1 t A 1 + b e k t 1 n + D z b k θ 0 e k t e θ 1 t 1 A n θ 1 1 + b e k t 1 n + 1                                         + D x a 0 + a 1 t + a 2 t 2 + a 3 t 3 A 1 + b e k t 1 n                                                                                     + D x b k e k t a 0 t + a 1 t 2 2 + a 2 t 3 3 + a 3 t 4 4 A n 1 + b e k t 1 n + 1 = 0
T C U i y = 0 y = 2 D S r p i 1 α 2 A 2 1 + b e k t 2 n
By substituting the order size from Equation (13) into Equation (12), the optimal slaughter age can be obtained by the Newton–Raphson method. It is necessary to remark that t i , y i t , and T C U i t are the slaughter age, the order size, and the total cost at the breakpoint i ( u i y i t < u i + 1 ), respectively. Therefore, the steps of the proposed solution procedure are as follows:
  • Step 1. For i = 1,2 J
    (a)
    Substituting the value of Equation (13) into Equation (12), find the value of t i from Equation (12) using the Newton–Raphson method.
    (b)
    Based on the value of t i * , compute the value of y i using Equation (13).
    (c)
    If u i > y i , then y i * = u i , else if u i + 1 < y i , then y i * = u i + 1 1 , and hence go to Step 1 (d). Otherwise, y i * = y i ( u i y i < u i + 1 ), then go to Step 1 (e).
    (d)
    Based on the value of y i , find the value of t i from Equation (12) using the Newton–Raphson method
    (e)
    Having t i * and y i * , obtain T C U i * from Equation (7)
  • Step 2. Compare the total cost ( T C U i * ) of the solutions ( t i * , y i * ) for i = 1,2 J . The pair ( t i * , y i * ) that gives the lowest cost is the optimal solution for the problem.

3.3. EGQ with Known Slaughter Age

As stated before, the slaughter age is defined by market behavior and customer purchase data. Therefore, this model has only one decision variable, which is the order size ( y ). Now, the optimal order size is obtained by setting the derivation of objective function (8) regarding y i equal to zero, as follows:
T C U i y = D S 1 α y 2 m t + r p i 1 α m t 2 = 0
Therefore, the solution is:
y = 2 D S r p i 1 α 2 m t 2
According to Equation (15), the order size value is like the classic EOQ model with all-units discount. Based on the classic EOQ, the solution procedure is as follows:
  • Step 1. Define the slaughter age and compute the values of m t , f ( t ) , and v ( t ) by using Equations (1)–(3), respectively.
  • Step 2. Set i = J
  • Step 3. Determine the optimal value of the breakpoint i ( y i ) using Equation (15)
  • Step 4. If u j > y i , then y i * = u j and go to Step 5. Otherwise, y i * = y i ; then, go to Step 6.
  • Step 5. Set i = i 1 and go to Step 3
  • Step 6. Obtain the optimal value of T C U i * form Equation (8) according to all the y i * s obtained.
  • Step 7. Compare the total cost ( T C U i * ) of the optimal solution ( y i * ) for every breakpoint ( J , J 1 , . i ). The y i * that provides the lowest cost is the optimal solution for the problem.

4. Numerical Examples

As is clear from the previous descriptions, the model with the all-units discount is a more comprehensive and complete model than the model without the discount, and its solutions will always be better or equal to the case without the discount. Additionally, the last model, which has a general discount and a specific slaughter age, is a special case of the model with a general discount, in which the slaughter time is known in advance based on market needs and customer requests. To demonstrate the application of the solution procedures for the three proposed models, a general numerical instance for a particular type of broiler inventory system is developed in this section. The general instance has been considered for all models; however, some parameters may change in each model. Based on the report of the European Food Safety Authority [33], farms can slaughter broilers from 21 days up to 150 days old. The usual slaughter age is 47 days in the US, while the slaughter age is 42 days in European countries. Therefore, the values of the minimum ( O ) and maximum ( W ) allowable length of breeding period are set equal to 21 and 55 respectively. The numerical instance uses the parameters of Richard’s growth curve and the feed consumption curve that are estimated by Goliomytis et al. [30], as follows:
k = 0.036 ; n = 0.0087 ; A = 6870.2 ; b = 0.043 , and
a 0 = 532.2 ; a 1 = 67.15 ; a 2 = 0.651 ; a 3 = 0.0018 .
Accordingly, the growth curve and the feed consumption curve applied for this instance are as follows (see Figure 2 and Figure 3):
m t = 6870.2 1 + 0.043 e 0.036 t 1 0.0087 ; and
f t = 532.2 + 67.15 t 0.651 t 2 + 0.0018 t 3 .
Moreover, the exponential function v t for NH3 production (Microliter/day·m2·bird) is based on the work of Leonard et al. [32] with b 0 = 19.44 and b 1 = 0.078 . Thus, the ammonia production curve is v t = 19.44 e 0.078 t (see Figure 4). The rest of the parameters in the general instance are D = 1,000,000,000   g , m 0 = 45   g , r = 0.02 / year , x = 0.0001   $ / g / day , S = 5.000   $ / cycle , and z = 0.00001   $ / μ L / d · m 2 · bird .

4.1. EGQ with No Discount

Based on the general model parameters given above, the purchasing cost per gram is p = 0.01   $ / g . The optimal solutions for the EGQ model with no discount are obtained with the Section 3.1 solution procedure. Applying the Newton–Raphson method, the optimal slaughter age is t * = 53.26 53 . Next, substituting t * into Equations (11) and (5), the optimal order size and total cost are determined as follows: y * 75,472 chicks, and T C U * 3,091,349 $.

4.2. General EGQ

The same conditions and parameters from the EGQ model with no discount are used for the general EGQ model, except for one difference: the supplier sells newborn chicks at different prices based on the quantity ordered. Thus, p = ( p 1 , p 2 , p 3 , p 4 ) = ( 0.010,0.008,0.006,0.004 )   $ / g and u = ( u 1 , u 2 , u 3 , u 4 ) = ( 0,50000,80000,100000 )   g . Applying the proposed solution procedure explained in Section 3.2, the optimal solution for each breakpoint is as follows:
  • First Point ( 0 y 50,000 ): t 1 * 53 , y 1 * 50,000 ; and T C U 1 * 3,095,060 $.
  • Second Point ( 50,000 y 80,000 ): t 2 * 53 , y 2 * 80,000 ; and T C U 2 * 3,056,307 $.
  • Third Point ( 80,000 y 100,000 ): t 3 * 52 , y 3 * 100,000 ; and T C U 3 * 3,020,000 $.
  • Fourth Point ( y 100,000 ): t 4 * 52 , y 4 * 122,605 ; and T C U 4 * 2,982,433 $.
  • Breakpoint 4 produces the lowest total cost ( T C U 4 * 2,982,433 ), so the optimal solutions are t * = 52 days, y * = 122,605 chicks, and T C U * = 2,982,433 $.

4.3. EGQ with Known Slaughter Age

The EGQ model conditions and parameters (Section 4.2) are applied for this last numerical example, just with one difference: the slaughter age is considered equal to 62 days. Now, the optimal solution is easily obtained as follows by applying the solution procedure detailed in Section 3.3:
  • Step 1. Using Equations (1)–(3)
m t = 3959 g, f t = 2614.5 , and v t = 2265 .
  • Step 2 to 3. Using Equation (15)
u 4 = 0.004 y 4 = 97,303 ;
u 3 = 0.006 y 3 = 79,448 ;
u 2 = 0.008 y 2 = 68,804 ;
  • Step 4 to 5.
Since y 4 * < u 4 , then y 4 * = u 4 = 100,000 .
Since y 3 * < u 3 , then y 3 * = u 3 = 80,000 .
Since u 2 y 3 * < u 3 , so the value of y 2 * is equal to y 2 = 68,804 .
It is not necessary to check the lower breakpoint ( i = 1 ) as the order size value is between its range.
  • Step 6. Using Equation (8)
T C U 2 * = 3,114,003 , T C U 3 * = 3,083,876 ; and T C U 4 * = 3,052,760 .
  • Step 7.
Since T C U 4 * < T C U 3 * < T C U 2 * , so the optimal total cost is equal to T C U * = T C U 4 * = 3,052,760 $, and thus, the optimal order size is y * = y 4 * = 100,000 chicks.

5. Conclusions

NH3 in poultry is a crucial biological process effect of the protein, litter, and manure composition generated on farms. When ammonia gas is exposed to moisture, it reacts and forms a basic, corrosive solution named ammonium. This aqueous solution causes harm to broilers. NH3 production increases as poultry ages, so managers must use methods to regulate the temperature, adjust the ventilation rate, and regulate diet and feed supplements to counteract it. These methods impose additional costs on poultry farm operations, which must be included in the inventory models. In addition, the process of shipping newborn chicks from suppliers to farms causes a percentage of deaths due to stressors and incidents such as catching, loading, transporting, and unloading. Therefore, a general EGQ model with mortality in newborn chicks and ammonia production inhibition cost under an all-units discount policy is presented in this study. In this general model, the poultry farm acquires newborn items at a discount price, which is lower than the wholesale price if the quantity ordered is above a breakpoint. The optimal solutions are obtained by using an analytical approach. Additionally, two special cases of the general model are formulated and explained based on their widespread use in the real world. The first model, EGQ with no discount, assumes that the supplier does not offer any discount policy. This condition is common in real life as some items are special or in limited supply, so the supplier does not offer any discount for them. This model is also solved using an analytical approach. In the last model, the consumption weight of poultry items is known in advance in most countries due to customer interest and pre-existing purchasing data. Once the consumption weight is set, the slaughter age is calculated according to a biological weight function. Therefore, the second case is an EGQ with known slaughter age. It is proved that this model has a solution like the one of a classic EOQ model with an all-units discount policy.
The proposed models can be extended (a) by using amelioration and deterioration rates, (b) by considering an expiration date in the consumption period, (c) by introducing a mortality rate during the breeding period, among others. Additionally, this research can be extended to other farm operations like livestock and fishery. Finally, it is possible to investigate the environmental effects of other harmful gases produced in livestock and poultry farms.

Author Contributions

Conceptualization, A.H.N., E.N., L.E.C.-B., D.G.-N., G.T.-G., A.C.-M., I.d.J.L.-H. and N.R.S.; Methodology, A.H.N., E.N., L.E.C.-B., D.G.-N., G.T.-G., A.C.-M., I.d.J.L.-H. and N.R.S.; Software, A.H.N., E.N. and L.E.C.-B.; Validation, A.H.N., E.N., L.E.C.-B., D.G.-N., G.T.-G., A.C.-M. and I.d.J.L.-H.; Formal analysis, A.H.N., L.E.C.-B., D.G.-N., G.T.-G., A.C.-M. and N.R.S.; Investigation, A.H.N., E.N., L.E.C.-B., D.G.-N., A.C.-M., I.d.J.L.-H. and N.R.S.; Data curation, A.H.N., E.N., L.E.C.-B., G.T.-G., I.d.J.L.-H. and N.R.S.; Writing—original draft, A.H.N. and E.N.; Writing—review & editing, L.E.C.-B.; Supervision, L.E.C.-B. 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

Data are in the paper.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. The Proof of Convexity of the EGQ Model Without Discount

Proving the convexity of this problem mathematically is a difficult task, but by considering a hypothetical situation, it can be proved that this problem is convex. Consider the example given in Section 4, if the number ordered is equal to 75,000, then the objective function for different values of bird ages is shown in Figure A1, so the function is convex.
On the other hand, if the age of the bird is given, the second derivative of the objective function (5) with respect to the number of orders becomes Equation (A1). As a result, the objective function is always positive with respect to y , and subsequently, the function is convex. In addition, the shape of the objective function with respect to both decision variables is shown in Figure A2. Therefore, it can be said that the EGQ problem is convex.
T C U y = 2 D S 1 α y 3 A 1 + b e k t 1 n 0
Figure A1. The graph of objective function (5) with respect to the bird age.
Figure A1. The graph of objective function (5) with respect to the bird age.
Sustainability 15 08086 g0a1
Figure A2. The graph of objective function (5) for the numerical example.
Figure A2. The graph of objective function (5) for the numerical example.
Sustainability 15 08086 g0a2aSustainability 15 08086 g0a2b

Appendix B. The Calculation of the Number of Orders of the General EGQ

The partial derivative of objective function (7) with respect to the order size becomes Equation (B1).
T C U i y = 2 D S 1 α y 2 A 1 + b e k t 1 n                                                         + r p i 1 α A 1 + b e k t 1 n 2 = 0
According to Equation (B2), the value of the order size is defined as:
y = 4 D S r p i 1 α 2 A 2 1 + b e k t 2 n

References

  1. Rezaei, J. Economic order quantity for growing items. Int. J. Prod. Econ. 2014, 155, 109–113. [Google Scholar] [CrossRef]
  2. Zhang, Y.; Li, L.Y.; Tian, X.Q.; Feng, C. Inventory management research for growing items with carbon-constrained. In 2016 35th Chinese Control Conference (CCC), Chengdu, China, 27–29 July 2016; IEEE: Piscataway, NJ, USA, 2016. [Google Scholar]
  3. De-la-Cruz-Márquez, C.G.; Cárdenas-Barrón, L.E.; Mandal, B. An inventory model for growing items with imperfect quality when the demand is price sensitive under carbon emissions and shortages. Math. Probl. Eng. 2021, 2021, 6649048. [Google Scholar] [CrossRef]
  4. Rana, K.; Singh, S.R.; Saxena, N.; Sana, S.S. Growing items inventory model for carbon emission under the permissible delay in payment with partially backlogging. Green Financ. 2021, 3, 153–174. [Google Scholar] [CrossRef]
  5. Gharaei, A.; Almehdawe, E. Optimal sustainable order quantities for growing items. J. Clean. Prod. 2021, 307, 127216. [Google Scholar] [CrossRef]
  6. Nobil, A.H.; Sedigh, A.H.A.; Cárdenas-Barrón, L.E. A generalized economic order quantity inventory model with shortage: Case study of a poultry farmer. Arab. J. Sci. Eng. 2019, 44, 2653–2663. [Google Scholar] [CrossRef]
  7. Nobil, A.; Taleizadeh, A.A. Economic Order Quantity for Growing Items with Discrete Orders. J. Model. Eng. 2019, 17, 123–129. [Google Scholar]
  8. Khalilpourazari, S.; Pasandideh, S.H.R. Modeling and optimization of multi-item multi-constrained EOQ model for growing items. Knowl.-Based Syst. 2019, 164, 150–162. [Google Scholar] [CrossRef]
  9. Malekitabar, M.; Yaghoubi, S.; Gholamian, M.R. A novel mathematical inventory model for growing-mortal items (case study: Rainbow trout). Appl. Math. Model. 2019, 71, 96–117. [Google Scholar] [CrossRef]
  10. Gharaei, A.; Almehdawe, E. Economic growing quantity. Int. J. Prod. Econ. 2020, 223, 107517. [Google Scholar] [CrossRef]
  11. Sebatjane, M.; Adetunji, O. Optimal inventory replenishment and shipment policies in a three-echelon supply chain for growing items with expiration dates. Opsearch 2022, 59, 809–838. [Google Scholar] [CrossRef]
  12. Sebatjane, M.; Adetunji, O. Economic order quantity model for growing items with imperfect quality. Oper. Res. Perspect. 2019, 6, 100088. [Google Scholar] [CrossRef]
  13. Nobil, A.H.; Nobil, E.; Cárdenas-Barrón, L.E.; Garza-Núñez, D.; Treviño-Garza, G.; Céspedes-Mota, A.; Loera-Hernández, I.D.J.; Smith, N.R. Economic Order Quantity for Growing Items with Mortality Function under Sustainable Green Breeding Policy. Mathematics 2023, 11, 1039. [Google Scholar] [CrossRef]
  14. Khan, M.A.A.; Cárdenas-Barrón, L.E.; Treviño-Garza, G.; Céspedes-Mota, A.; Loera-Hernández, I.D.J. Integrating prepayment installment, pricing and replenishment decisions for growing items with power demand pattern and non-linear holding cost under carbon regulations. Comput. Oper. Res. 2023, 156, 106225. [Google Scholar] [CrossRef]
  15. Sebatjane, M.; Adetunji, O. Economic order quantity model for growing items with incremental quantity discounts. J. Ind. Eng. Int. 2019, 15, 545–556. [Google Scholar] [CrossRef]
  16. Sebatjane, M.; Adetunji, O. Three-echelon supply chain inventory model for growing items. J. Model. Manag. 2019, 15, 567–587. [Google Scholar] [CrossRef]
  17. Sebatjane, M.; Adetunji, O. Optimal inventory replenishment and shipment policies in a four-echelon supply chain for growing items with imperfect quality. Prod. Manuf. Res. 2020, 8, 130–157. [Google Scholar] [CrossRef]
  18. Sebatjane, M.; Adetunji, O. Optimal lot-sizing and shipment decisions in a three-echelon supply chain for growing items with inventory level-and expiration date-dependent demand. Appl. Math. Model. 2021, 90, 1204–1225. [Google Scholar] [CrossRef]
  19. Pourmohammad-Zia, N.; Karimi, B.; Rezaei, J. Dynamic pricing and inventory control policies in a food supply chain of growing and deteriorating items. Ann. Oper. Res. 2021, 304, 1–40. [Google Scholar] [CrossRef]
  20. Pourmohammad-Zia, N.; Karimi, B.; Rezaei, J. Food supply chain coordination for growing items: A trade-off between market coverage and cost-efficiency. Int. J. Prod. Econ. 2021, 242, 108289. [Google Scholar] [CrossRef]
  21. Sebatjane, M. The impact of preservation technology investments on lot-sizing and shipment strategies in a three-echelon food supply chain involving growing and deteriorating items. Oper. Res. Perspect. 2022, 9, 100241. [Google Scholar] [CrossRef]
  22. Wang, Y.; Wang, S.; Jing, H.; Zhang, T.; Song, N.; Xu, S. CircRNA-IGLL1/miR-15a/RNF43 axis mediates ammonia-induced autophagy in broilers jejunum via Wnt/β-catenin pathway. Environ. Pollut. 2022, 292, 118332. [Google Scholar] [CrossRef] [PubMed]
  23. Koerkamp, P.G. Review on emissions of ammonia from housing systems for laying hens in relation to sources, processes, building design and manure handling. J. Agric. Eng. Res. 1994, 59, 73–87. [Google Scholar] [CrossRef]
  24. Miles, D.M.; Branton, S.L.; Lott, B.D. Atmospheric ammonia is detrimental to the performance of modern commercial broilers. Poult. Sci. 2004, 83, 1650–1654. [Google Scholar] [CrossRef] [PubMed]
  25. Webb, J.; Menzi, H.; Pain, B.F.; Misselbrook, T.H.; Dämmgen, U.; Hendriks, H.; Döhler, H. Managing ammonia emissions from livestock production in Europe. Environ. Pollut. 2005, 135, 399–406. [Google Scholar] [CrossRef]
  26. Wang, J.; Zhang, P.; Wang, S.; Yang, L.; Luo, J.; He, W.; Du, G.; Wang, X.; Zhang, Z.; Yang, X. Evaluation of a green-sustainable industrialized cleaner production for FeV50 and FeV80 alloys from vanadium slag by calcification roasting–ammonia on-line cycle. J. Clean. Prod. 2021, 320, 128896. [Google Scholar] [CrossRef]
  27. Brink, M.; Janssens, G.P.; Demeyer, P.; Bağci, Ö.; Delezie, E. Reduction of dietary crude protein and feed form: Impact on broiler litter quality, ammonia concentrations, excreta composition, performance, welfare, and meat quality. Anim. Nutr. 2022, 9, 291–303. [Google Scholar] [CrossRef] [PubMed]
  28. Yılmaz, A.R.A.L.; Arikan, M.S.; Akin, A.C.; Kuyululu, Ç.Y.K.; Güloğlu, S.C.; Sakarya, E. Economic losses due to live weight shrinkage and mortality during the broiler transport. Ank. Üniversitesi Vet. Fakültesi Derg. 2014, 61, 205–210. [Google Scholar] [CrossRef]
  29. De-la-Cruz-Márquez, C.G.; Cárdenas-Barrón, L.E.; Mandal, B.; Smith, N.R.; Bourguet-Díaz, R.E.; Loera-Hernández, I.D.J.; Céspedes-Mota, A.; Treviño-Garza, G. An inventory model in a three-echelon supply chain for growing items with imperfect quality, mortality, and shortages under carbon emissions when the demand is price sensitive. Mathematics 2022, 10, 4684. [Google Scholar] [CrossRef]
  30. Goliomytis, M.; Panopoulou, E.; Rogdakis, E. Growth curves for body weight and major component parts, feed consumption, and mortality of male broiler chickens raised to maturity. Poult. Sci. 2003, 82, 1061–1068. [Google Scholar] [CrossRef]
  31. Shah, S.W.A.; Chen, D.; Zhang, J.; Liu, Y.; Ishfaq, M.; Tang, Y.; Teng, X. The effect of ammonia exposure on energy metabolism and mitochondrial dynamic proteins in chicken thymus: Through oxidative stress, apoptosis, and autophagy. Ecotoxicol. Environ. Saf. 2020, 206, 111413. [Google Scholar] [CrossRef]
  32. Leonard, J.J.; Feddes, J.J.R.; McQuitty, J.B. Air quality in commercial broiler housing. Can. Agric. Eng. 1984, 26, 65–72. [Google Scholar]
  33. European Food Safety Authority. Technical assistance to the Commission (Article 31 of Regulation (EC) No 178/2002) for the preparation of a data collection system of welfare indicators in EU broilers’ slaughterhouses. EFSA J. 2013, 11, 3299. [Google Scholar]
Figure 1. The weight inventory level of broilers.
Figure 1. The weight inventory level of broilers.
Sustainability 15 08086 g001
Figure 2. Richard’s growth curve.
Figure 2. Richard’s growth curve.
Sustainability 15 08086 g002
Figure 3. Feed consumption curve.
Figure 3. Feed consumption curve.
Sustainability 15 08086 g003
Figure 4. Ammonia curve.
Figure 4. Ammonia curve.
Sustainability 15 08086 g004
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

Nobil, A.H.; Nobil, E.; Cárdenas-Barrón, L.E.; Garza-Núñez, D.; Treviño-Garza, G.; Céspedes-Mota, A.; Loera-Hernández, I.d.J.; Smith, N.R. Economic Growing Quantity Model with Mortality in Newborn Items and Inhibition Cost of Ammonia Production under All-Units Discount Policy. Sustainability 2023, 15, 8086. https://doi.org/10.3390/su15108086

AMA Style

Nobil AH, Nobil E, Cárdenas-Barrón LE, Garza-Núñez D, Treviño-Garza G, Céspedes-Mota A, Loera-Hernández IdJ, Smith NR. Economic Growing Quantity Model with Mortality in Newborn Items and Inhibition Cost of Ammonia Production under All-Units Discount Policy. Sustainability. 2023; 15(10):8086. https://doi.org/10.3390/su15108086

Chicago/Turabian Style

Nobil, Amir Hossein, Erfan Nobil, Leopoldo Eduardo Cárdenas-Barrón, Dagoberto Garza-Núñez, Gerardo Treviño-Garza, Armando Céspedes-Mota, Imelda de Jesús Loera-Hernández, and Neale R. Smith. 2023. "Economic Growing Quantity Model with Mortality in Newborn Items and Inhibition Cost of Ammonia Production under All-Units Discount Policy" Sustainability 15, no. 10: 8086. https://doi.org/10.3390/su15108086

APA Style

Nobil, A. H., Nobil, E., Cárdenas-Barrón, L. E., Garza-Núñez, D., Treviño-Garza, G., Céspedes-Mota, A., Loera-Hernández, I. d. J., & Smith, N. R. (2023). Economic Growing Quantity Model with Mortality in Newborn Items and Inhibition Cost of Ammonia Production under All-Units Discount Policy. Sustainability, 15(10), 8086. https://doi.org/10.3390/su15108086

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