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 CO
2 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 (NH
3) as a by-product. Additionally, excreted nitrogen (N) is associated with deteriorating litter quality, which exists in NH
3 content and contributes to the spread of hock and foot lesions in broilers. High NH
3 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], NH
3 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 NH
3 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 NH
3 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.
Paper | Year | Objective Function | Solution Method | Growth Curve | Dead Newborn Animals | NH3 Emission | Discount | Types of Items |
---|
Rezaei [1] | 2014 | PF | Analytical | BW | | | | P |
Zhang et al. [2] | 2016 | CF | Analytical | BW | | | | P |
Sebatjane and Adetunji [12] | 2019 | PF | Analytical | BW/LR | | | | P |
Sebatjane and Adetunji [15] | 2019 | CF | Analytical | BW | | | ID | L |
Sebatjane and Adetunji [16] | 2019 | CF | Analytical | BW | | | | L |
Khalilpourazari and Pasandideh [8] | 2019 | PF | Meta-heuristic | BW | | | | P |
Malekitabar et al. [9] | 2019 | PF | Game theory | BW | * | | | F |
Nobil et al. [6] | 2019 | CF | Analytical | LR | | | | P |
Nobil and Taleizadeh [7] | 2019 | CF | Closed form | LR | | | | P |
Sebatjane and Adetunji [17] | 2020 | PF | Analytical | BW | | | | L |
Gharaei and Almehdawe [10] | 2020 | CF | Closed form | LR | | | | P |
De-la-Cruz-Márquez et al. [29] | 2021 | PF | Analytical | BW | | | | P |
Gharaei and Almehdawe [5] | 2021 | CF | Meta-heuristic | BW | | | | P |
Pourmohammad-Zia et al. [19] | 2021 | PF | Analytical | BW | | | | P |
Pourmohammad-Zia et al. [20] | 2021 | PF | Game theory | BW | | | | P |
Rana et al. [4] | 2021 | CF | Analytical | BW | | | | P |
Sebatjane and Adetunji [19] | 2021 | PF | Analytical | BW | | | | P |
De-la-Cruz-Márquez et al. [3] | 2022 | PF | Analytical | BW | | | | P |
Sebatjane and Adetunji [11] | 2022 | CF | Analytical | BW | | | | P |
Sebatjane [21] | 2022 | CF | Analytical | BW | | | | P |
Nobil et al. [13] | 2023 | CF | Analytical | BW | | | | P |
Khan et al. [14] | 2023 | PF | Analytical | BW | | | | P |
This paper | | CF | Analytical | BW | * | * | AD | P |
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 and , respectively, and the rest of notation of the model is as follows:
: the number of discount categories based on a price;
: the weight of a broiler at time (weight unit)
: the highest length of breeding time (day)
: the lowest length of the breeding time (day)
: 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.
: the holding cost percentage rate per year.
: the weight volume breakpoint in units for a unit purchasing cost (weight unit)
: the purchasing cost per weight unit for price range (monetary unit/weight unit) if
: the setup cost per growing cycle (monetary unit/setup)
: the cost of production (feeding) during growing period per unit item (monetary unit/unit item)
: the inhibition cost of ammonia production (monetary unit/microliter/day·square meter·weight unit)
: the polynomial function of feeding consumption (unit items)
: the exponential function of ammonia production (microliter/day·square meter·bird).
: the growth rate.
: the asymptotic weight.
: the integration constant of the growth function.
: the shape parameter of the growth function.
: the demand periods.
: the total weight at the slaughter time (weight unit).
: the total cost per year (monetary unit).
: the total cost of breakpoint per year (monetary unit).
: the total number of growth items ordered at the beginning of each cycle (unit items)
: the slaughter age (day)
Figure 1 illustrates the weight stock level of broilers in the proposed model. In this inventory system,
newborn chicks with an initial weight
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
. The poultry farm raises them for a period until they reach the ideal slaughter weight (
). At the time of slaughter (
), all broilers begin processed and shipped to the market for consumption, during the consumption period (
).
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:
where
is the body weight of a live bird at age
.
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:
where
,
,
, and
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, NH
3, emissions in poultry farms. NH
3 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:
where
and
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 (
) is received with initial weight (
) from the supplier, so the purchasing cost per growth cycle is equal to
, where
is the purchasing cost per weight unit without discount. Additionally, the length of the consumption period is
, while the average weight of the inventory during the consumption period is
. The holding cost per cycle is
, where
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
and
are considered as the feeding cost per unit and the length of the breeding period, respectively. The total feeding cost per cycle is
=
.
Finally, there is a cost of dealing with ammonia production defined by
. Therefore, the total cost per cycle with no discount is as follows:
To meet the total cost function per year,
should be divided by
. This means that the total cost function per year without discount is:
where
and
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:
where
is the order quantity,
and
.
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 (
) is defined as follows:
where
and
are the decision variables and the price (
) depends on the order size of
.
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 (
) 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
, and
are obtained by defining the slaughter age from Equations (2) and (3), respectively. Based on the value of
, the total cost in Equation (7) is adapted as follows:
where
is the only decision variable, and
is a dependent variable.
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 (
) and maximum (
) 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:
, and
.
Accordingly, the growth curve and the feed consumption curve applied for this instance are as follows (see
Figure 2 and
Figure 3):
; and
.
Moreover, the exponential function
for NH
3 production (Microliter/day·m
2·bird) is based on the work of Leonard et al. [
32] with
and
. Thus, the ammonia production curve is
(see
Figure 4). The rest of the parameters in the general instance are
,
,
,
,
, and
.
4.1. EGQ with No Discount
Based on the general model parameters given above, the purchasing cost per gram is
. 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
. Next, substituting
into Equations (11) and (5), the optimal order size and total cost are determined as follows:
chicks, and
$.
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,
and
. Applying the proposed solution procedure explained in
Section 3.2, the optimal solution for each breakpoint is as follows:
First Point (): , ; and $.
Second Point (): , ; and $.
Third Point (): , ; and $.
Fourth Point (): , ; and $.
Breakpoint 4 produces the lowest total cost (), so the optimal solutions are 52 days, chicks, and $.
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:
g, , and .
;
;
;
Since , then .
Since , then .
Since , so the value of is equal to .
It is not necessary to check the lower breakpoint () as the order size value is between its range.
, ; and .
Since , so the optimal total cost is equal to $, and thus, the optimal order size is 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.