Research Article  Open Access
Sujit Kumar De, Shib Sankar Sana, Adrijit Goswami, "An EOQ Model for Phase Inventory with Induced Demand and Periodic Cycle Time", Journal of Industrial Engineering, vol. 2014, Article ID 605178, 14 pages, 2014. https://doi.org/10.1155/2014/605178
An EOQ Model for Phase Inventory with Induced Demand and Periodic Cycle Time
Abstract
This paper deals with a stock flow of an inventory problem over induced demand. The inventory is consumed through “core customer” or chain marketing system in an induced environment (inductance) to exhaust all the items of the stock inventory in an indefinite time. The demand rate is depicted due to induced factor which is generated from the same inventory presented nearby. The inventory cycle time is split into several periodic times due to oscillatory feature of the inventory which is called phase inventory. Considering uniform demand, this cycle time splits into two basic parts, namely, “first shift” (phase) and “second shift” (phase). Since the process dampens over time, so the whole inventory will exhaust after few periods. A cost function consisted of inventory cost, setup cost, and loss for induced items is minimized to obtain optimal order quantity and replenishment time. The multivariate lagrange interpolation (MLI) over the average values of the postsensitivity analysis is developed here. Finally, graphical illustrations are made to justify the model.
1. Introduction
Recently, a great evolution has been made in the traditional EOQ model. Through its process, a survey work of the industrial system over Harris’s work [1] in the last century has been made by Andriolo et al. [2]. A capacitated lot sizing problem for coordination of transportation services is developed by Bruno et al. [3]. A review of consignment stocking policy under advertising, pricing, and bargaining strategies is also discussed by Sarker [4]. Moreover, several articles of the inventory in supply chain [5–11] have suggested that discounting on sales price is the emergent way to exhaust all items from the inventory in due course of time. Bozorgi et al. [12] develop a new model where the inventory moves with environmental temperature, controlled by emission function. Sodhi et al. [13] state the bullwhip effect on maintenance repair and overhaul (MRO) customers. Analysis of batch ordering inventory for infinite horizon with capacity constraints is studied by Yang et al. [14]. Coelho and Laporte [15] derive an inventory routing problem (IRP) for the inequalities, based on demand and available capacities. Also, Avinadav et al. [16] are able to develop a multiplicative demand function relating to age and price effect, whereas Chand and Sethi [17] made an analysis for dynamic and stationary demand in multiperiod lot sizing problem for both finite and infinite time horizons in the same time. Keskin and Capar [18] describe a decision support system model on the basis of the utility of three pillars of supply chain management, namely, locationtransportationinventory (LTI). Gambini et al. [19] investigate an EOQ model where stocking charge is independent of time and the stock flow down depends on its level only. Dolan [20] observes the effect of quantity discount on the joint EOQ model and suggests an optimal strategy. Monahan [21] develops a quantity discount model from the perspective of the seller with the constraints which are imposed to ensure sufficient benefit to the buyer. A continuous discount pricing model is explained by Dave et al. [22]. Tripathy and Pattanaik [23] analyze a fuzzy entropic model with stock dependent demand for perishable items. Syed and Aziz [24] develop a fuzzy price break model also. A model with perishable asset revenue management is studied by Weatherford and Bodily [25]. All of the above authors assure that the discounting on sales prices of the items in stock increases selling and it would be continued until the whole inventory is exhausted.
However, the present world economy is very much crucial and the uncertainty in marketing system is increasing day by day. In this situation, nobody can assure that salesman’s commodities would be purchased just after the opening of his/her shop. In our observations, such facts occur mainly due to two basic reasons.(a)Firstly, the purchasing capacity of the common people is decreasing day by day due to profit making business of multinational companies (MNC) throughout the world, especially in third world countries. As a result, a huge number of products are being piled as time progresses. This phenomenon is commonly known as “crisis of capitalism.”(b)The second one is the high competition among various companies and few reputations of them. A reputed company attracts customers easily rather than nonreputed ones as they have their own “core customer circle” (see Definition 2). These core customers may be influenced by other MNC at the time of purchasing the goods just to have separate tastes and feelings of new emerging items of the market. On the other hand, for a newly setup shop/company, naturally they have neither “reputation” nor “core customer circle” but they have to struggle through this system. Such types of companies/salesmen begin to start their business with few peer (core, friendrelative circle) customers. The customers (core/noncore) who are influenced or motivated by the advertisement through the printed/digital medias/sales teams are called “floating” or “surface” (see Definition 1) customers.
Behavioral study of the customers indicates that the floating customers would not always purchase goods at discounted price but to oppose many cases although they were the surface customers in few days ago. The general motives of the customers are “purchase quality goods at fair price”; their views on discounted items are usually as follows: “unfair, bad quality or very near to outdated time bar.” So, a negative motivation arises through the announcement “discounted items.” Hence, the unsold items are rotated periodically for an indefinite time.
The consumption rate of common assets like mobile phone is increasing day by day. At the same time, critical risks and challenges are observed among the various mobile service providers. The main challenges are coming from the frequent change of SIM (subscriber identity module) cards among the customers. For this reason, several companies are offering price discount to capture more customers in favor of them. In this way, for instance, Vodafone Company has captured India’s largest telephone company Bharat Sanchar Nigam Limited, BSNL’s market. Likewise, Airtel, MTS, DOCOMO, UNINOR, and so forth have captured partially Vodafone’s market. Under the circumstances, a SIM card seller of any company, after purchasing (manufacturing) a certain amount of items, will have to wait a long time until all the items are exhausted from its own inventory. On the other hand, the demand of assets/goods like car, motor bike, watch, ceiling fan, almirah, utensils, soaps, and so forth in a particular enterprise may fall under this consideration. Thus, the lingering time to exhaust all items is basically due to the presence of another shop (inventory setup that offers price discount or not) in which every customer (core/noncore) has multiple options to choose their better service provider. Now, we define floating customers as induced customers in which the products are not outflowed from the inventory system for opposed demand due to “induced factor,” and the inventory relating to this system is termed as phase or periodic inventory. Therefore, the behavior of the inventory function is quite similar to the discharging process in a capacitor of an electrical circuit through inductance. In this literature, several attempts in trade economics have been made by several researchers. Bryant’s [26] thermodynamic approach to the financial economics is quite familiar. Similarly, a classical thermodynamic and general equilibrium theory has been developed by Smith and Foley [27] and the role of money is developed by Soddy [28].
In our present study, we made an attempt for analyzing an inventory model under electrical circuit theory (ECT). Here, the order quantity is independent of time, so we have studied the MLI function over the postsensitivity analysis of the model to make a final conclusion.
2. Assumptions and Notations
The following notations and assumptions are used to develop the model.
Assumptions(1)The inventory starts at the beginning (time at the time phase difference) with maximum onhand inventory and depletes through natural system. After adjusting the demand of the customers at two phases of each cycle time, it reaches zero level at the end of the cycle time.(2)The time horizon is infinite (days).(3)Shortages are not allowed.(4)Demand rate is a time dependent phase inventory function.(5)Holding cost is uniform over the cycle time.(6)Inductance starts at the linear time phase difference and ends at the time of the twice multiple of linear time phase difference in the first shift of the last period.(7)Conducting factor is low.(8)Items are similar which are stored in the shops nearby.(9)Total number of customers for particular types of product is fixed.(10)Expiry date of the items is indefinite.
Notation : The maximum order quantity per cycle : Instantaneous demand at time in : The instantaneous inventory at time in : The damping factor : Linear time phase difference : Inventory holding cost per unit item per time unit ($) : Setup cost per cycle ($) : Cost of induced factor per unit induced demand item ($) : The number of days per cycle, a positive integer : Replenishment time in days : Capacity factor : Induced factor : Resisting index : Average inventory cost per period (day) ($).
3. Formulation of the Model
3.1. Stock Flow Inventory Model under Inductance
Let us consider an inventory chain having no external supply of goods. The circuit has some resisting factor due to the weak presence of core customers and induced factor due to the presence of other shops nearby. In this system, a seller after purchasing a maximum order quantity began to sale through the core customers with the presence of floating customers. The inventory is not exhausted after a long elapsed time whenever the core customers are absent in the system (see (3)–(6)), but the system dampens slowly the process and all the items would be exhausted few days later (see Figure 1).
Definition 1 (floating or surface customers (induced factor)). Such customers do not have any faith or belief in any kind of branded company. They can influence core customers and are influenced by the upcoming advertisement of several companies. Therefore, their presence in the inventory system can create a harmonious induced factor or threat in the system itself. It is defined as a percent of customers who do not purchase items due to the influences of other customers/shops within a certain interval of time.
Definition 2 (core customers (conducting factor)). Every inventory circuit has specific challenges or limitations (resisting factor ). Conversely, every inventory circuit (seller’s network) has its own peer/friend circle or chain marketing system (however small it may be) who will purchase goods/commodities from his (her) own system. They are called core customers. The resisting factor can be defined as a percent of customers who do not like to purchase the items from the particular inventory system due to dissatisfaction of service providers, lack of goodwill and misconduct from the selling personnel, and so forth. Therefore, the % of satisfied customers acts as a conducting factor of the inventory system.
Definition 3 (induced items). It is defined as the amount of goods that are not sold due to induced factor. If is the instantaneous demand rate at time , then the total amount of inventory that is unsold due to induced factor is given by where is the instantaneous inventory at time .
Definition 4 (capacity factor). It is used to define the degree in which the items are accommodating in the shop for sale in the indefinite time horizon. It is denoted by and is defined by resisting index , induced factor , and the damping factor of the system (see Table 1).

Definition 5 (shift frequency/time phase). Here, the frequency is the number of cycles (periods) that occurred in a particular time. Let, at least, one period occur in a day; more than one day may correspond to a single period also. It is defined as the total number of turnovers per day. We may assume that the inventory curve completes at least one cycle (reach a similar phase point) in a day and gives two different phase intervals in which the curve gets positive (first shift) and negative values (second shift) successively (in a day, repeated phase time is not possible in general).
So for a complete oscillation, the inventory curve will take just one day; hence, its frequency is one; that is, . Thus to have circular frequency we take
Definition 6 (periodic time ). It is defined as the time (day) required to complete one closed path (period) of the inventory curve and is denoted by If the inventory circuit has no core customers, then the periodic time is ; and . And
3.2. Multivariate Lagrange Interpolation (MLI)
Let be an variable multinomial function of degree . Since there are terms in , it is a necessary condition that we have distinct points , , , for to be uniquely defined. In other words, , where the are the coefficients in . is the tuple of independent variables of . ) is an exponent vector with nonnegative integer entries, consisting of an ordered partition of an integer between 0 and inclusively. . is the usual vector dot product; . Now we wish to write in the form , where is a multinomial function in the independent variables with the property that when is equal to th data value or , then and . Let us consider the system of linear equations From this system, we construct the sample matrix Assuming , we could solve for the coefficients in by inverting .
If is singular, then the coefficients of are not uniquely determined, characterizing the geometric configuration of the points so that appears to be an intricate problem.
Let . Substituting in , we have Next, substituting in , , we have the following matrix : Here, the appears twice in that means . In other words, at . By construction, that provides . Therefore,
Applying the above definitions and interpolation function, we will discuss our proposed model as follows.
Let the inventory start at a time phase point and exhaust at the first shift of the last day after covering full days of the time horizon . The induction starts after time later and ends at the time of the first shift of the last day. Then, the governing differential equation is with and at time .
Now, substituting in the above, we have Now, if no external supply is supplied to the chain, then (11) reduces to where is called the damping factor, and the initial conditions are given by The solution of (12) for is given by where Here, is the phase difference, is the linear distance, is the web length of the oscillatory curve, and the period of oscillation is given by Also, Again from (19),
Therefore, the oninventory is as follows:
We see (Figure 2) that the time domain of the first shift [ is +ve] and the second shift [ is −ve] is as follows: Now, the inventory holding cost for the first shift (phase) per cycle is The inventory holding cost for the second shift (phase) per cycle is given by Now using (14), the loss per cycle due to the induced factor is Now, using (23), (24), and (25), the average profit per day is given by Now using Table 1, our problem reduces to where
Particular Cases: Since the model is developed with the specific cases.
(1) Cases of Optimality
Case I: when induced factor is fixed and the resisting index and conducting factor vary.
Case II: when resisting index is fixed and the induced factor and conducting factor vary.
Case III: when conducting factor is fixed and the resisting index and induced factor vary.
Case IV: when conducting factor is fixed and the resisting index and induced factor vary.
All measurements with different parameters are summarized in Table 1.
4. Numerical Example
We have considered the values of the parameters in appropriate units as follows.
Example 1. Let the values of parameters for Case I be , , , and . Then, we have the result (see Table 2).

Example 2. We consider the values of parameters for Case II in appropriate units as follows: , , , and . Then, we have the result (see Table 3).

Example 3. For Case III, let , , , and . Then, the result (see Table 4) is obtained.

Example 4. In Case IV, we consider , , , and . Then, we have the result (see Table 5).

4.1. Discussion of Tables 2–5
From Table 2, we see that if we take 20% induced factor over the whole time period, then the minimum average cost is $1207.94 at total conducting factor 8.74 and resisting index 30%, exhausting all the items within 20.03 days. In Table 3, at 20% resisting index, the average minimum cost decreases to $747.76 for total low induced factor 3% and the total conducting factor 3.16, exhausting all the items within 5.02 days. Table 4 shows that if the total conducting factor is 15, then the average minimum cost comes down to $355.43 for the total induced factor 0.12% and the resisting index 0.08%, exhausting all the items within 6.38 days. However, Table 5 shows that the minimum average cost is $303.74 at very high resisting factor (69.05) and average induced factor (11.80) for negligible conducting factor (0.009), exhausting all the items within 5.06 days.
5. Sensitivity Analysis
Considering , , , , and with appropriate units in the model and taking the range of variations of the parameters from −50% to +50% we have obtained optimal Table 6.

5.1. Implication to Introduce Multivariate Lagrange Interpolation (MLI)
In the sensitivity analysis table, we usually see the maximum profit/minimum cost of the model for specific change of a particular parameter. This kind of result is quite limited to its local limitations. However, to have a global impact of the model we require a general objective function for the sensitivities of the whole parameters imparting in the system. But to increase our level of confidence, we need to know the range of the average cost function, the essential parts of any kind of inventory management.
From the solution Tables 2–5, it is observed that the average cost attains its minimum value at the fixed values of conducting factor . As a result, we have considered the sensitivity analysis of and we have constructed the MLI polynomial function in this case also.
First we construct the MLI for the six datasets that are shown in Table 7.

Using Section 3.2, we have and using (6), the sample matrix is and the submatrices are as follows: Then, the corresponding determinants are Therefore, we have Similarly, we have for thedata sets that Then, the sample matrix is and the submatrices are as follows: Now, the corresponding determinants are Therefore, we have Here, considering average conducting factor , our problem reduces to To solve the above, we will maximize and separately subject to the given condition and obtain the result in Table 8.

5.2. Interpretation on Optimal Values in Table 8
We have constructed an MLI function to make a significant decision over the model, so we need not discuss the other tables. From Table 8, we see that the range of the average cost is [$2878.99, $3551.95], for the range of induced factor [0.86%, 2.12%], the range of resisting index is [4.80%, 7.52%], and all quantity exhausted time range is [8.71 days, 12.00 days]. This information can be helpful enough for a decision maker of a global trading system.
5.3. Graphical Interpretation of the Model
From Figures 3 and 4, it is observed that the average minimum cost is nominal (Figure 3) whenever we take the induced factor around 2.5% and the exhaust time is nearer to 5 days for the constant values of the resisting index and conducting factor . Moreover, the low induced factor (around 0.5%) and long exhausting time cause a high cost when the and values are kept constant.
Figure 4 shows that the average minimum cost occurs at around 4% resisting index and the inventory level is empty within 14 days for constant values of the conducting factor and induced factor . The average cost is around $6000, whenever the resisting index lies around 8% near 5day exhausting time for the fixed values of and . Analyzing both the figures, we may conclude that Figure 3 can give more economic result than the other.
6. Conclusion
In oligopolistic marketing system, competitions among the companies/industries are very much tuff. In these situations, it is hardly possible to survive and sustain a company under globalized trading system. The traditional concept was that the company will build inventory setup to provide customers’ needs under smooth and nice path without getting any challenge from other setups. Breaking this myth, we have observed in a bird’s view that the present time is actually of “sleep with challenge and awake to challenge others” or “hire and fire.” The proposed model helps the industrial engineers in knowing how to control the inventory model under several limitations which are newly introduced in the inventory literature.
The basic features of our model are as follows.(1) Peoples’ motivations result in induced factor form generation, means demand rate becomes oscillatory in nature that resulting the inventory circuit as discharging of a condenser in an electrical circuit.(2)Key of the inventory chain indicates whether the shop is open or closed. No induced factor will come if the shop is closed, but it will take little time (phase time difference) to generate a considerable induced factor.(3)As the inventory is depicted with time, the induced factor began to die out and its presence will disappear in the first shift at the end of cycle time which is a common phenomenon of any kind of physical system.(4)Frequency of the oscillatory inventory curve generates phases, namely, first shift for a day time and the second shift for the night or other similar means.(5)The definitions like core customers, surface customers, and induced and conducting factors have been discussed under competitive sphere of market and investments.(6)Numerical multivariate interpolation has been utilized to overcome any kind of discrepancies over the decision variables of the model.The major limitations of the model are as follows.(1)The model is of multivalued decision variables, so no finite search direction will be observed. For this reason we have to rely on some line search technique also.(2)In many cases the relational of equations among the variables , , and may not effect remarkably on the cost optimization model.(3)Several definitions may not support the socalled “inventory management problems.” Therefore, our model of study represents a special category of EOQ model, and we rename it as phase inventory model.
7. Scope of Future Work
Anybody may link up the inventory process with the trade economics, thermodynamic process, and the electrical circuit theory in a new outlook also. A concept of steady and equilibrium production and supply system can be chalked out from this paper.
Appendix
(1) The time distance for a complete oscillation (wave length) is calculated as follows:
Again, we have , where and are measured in time units.
(2) To discuss the dimension of both sides of the following equation: we take So,
(3) Using (12), the simplified form of (17) is obtained as follows: Therefore,
(4) Using (12), the simplified form of (19) is obtained as follows: Thus, Now,