﻿ Quantification of Node Wise Commodity in Supply Chain and Its Cost Analysis

American Journal of Operations Research
Vol.07 No.01(2017), Article ID:73808,19 pages
10.4236/ajor.2017.71005

Quantification of Node Wise Commodity in Supply Chain and Its Cost Analysis

Prem Prakash Mishra1, S. S. Mishra2*, S. K. Yadav2, R. S. Singh3, Ravendra Kumar4

1Department of Mathematics, National Institute of Technology, Chumukedima, Nagaland, India

3Kamla Nehru Institute of Management, Dr. APJ Abdul Kalam Technical University, Lucknow, India

4Department of Mathematics, Dr. Ram Manohar Lohia Government Degree College, Aonla, India    Received: November 11, 2016; Accepted: January 21, 2017; Published: January 24, 2017

ABSTRACT

In this paper, we focus on the intermediate nodes of network and quantification of level of commodity and its cost on each node because intermediate nodes have stocking capacities which we generally see in the supply chain network. The commodity is supplied from a node to node in response to the power form of demand at a particular time. Since the traffic intensity of the demand of commodity also affects the flow of the commodity in the network, hence study of flow of commodity in the network is believed to be a significant contribution in this area. Several cases of quantifying the level of commodity in different situations as well as the cost analysis of incoming and outgoing commodity at a particular node have been thoroughly discussed in the paper. The present problem, presumably seeks to contribute to managerial decision making in supply chain network.

Keywords:

Supply Chain, Traffic Intensity, Node Wise Amount of Flow and Cost Analysis 1. Introduction

A fresh attempt needs to be made to quantify the level of commodity on the nodes of supply chain network at any instant of time with a particular amount of flow. Here, we seek to characterize the nodal properties of the network which in turn helps in making the strategies for future plans about a particular node of network and gauge the degree of importance of a particular node in whole network. Node wise quantification commodity paves the way for obtaining the economic inflow quantity at the node as well as the total optimal cost of the network system. In this sequence of knowledge, we wish to report some important developments which are relevant to it. Ahuja et al.  have combined the several techniques to produce an algorithm for non-sparse graphs with large maximum arc capacity and also discussed about the bound for the transportation problems. Ahuja et al.  defined the flow in a network N (V, E) as a real function $x:V×V\to R$ with the three properties

1) Capacity constraints: ${x}_{ij}\le C\left(i,j\right)$ . The flow along an edge cannot exceed its capacity, $\forall i,j\in V$ .

2) Skew symmetric: ${x}_{ij}=-{x}_{ji}.$ The flow from $i$ to $j$ is ${x}_{ij}$ and the flow from $j$ to $i$ is ${x}_{ij}$ .

3) Flow conservation: $\sum _{j\in \text{\hspace{0.17em}}\text{\hspace{0.17em}}V}{x}_{ij}=0$ unless $i$ = source or $j$ = sink.

The net flow to a node is zero except for the source that “produces” flow and the sink which “consumes” flow.

This paper focuses on the third condition at the intermediate nodes of the network where summation of all incoming flow and outgoing flow is zero, but in some case of the supply chain network this conservation law at instantaneous incoming and outgoing of flow does not hold good. In this situation, value of the flow is given as below at the intermediate node which appears like source or sink for short length of time that is given below, where ${x}_{i}$ is called the amount of flow. The amount of flow at ${i}^{\text{th}}$ node is treated as variable in case of the random network of the supply chain. During the random alteration in incoming flow and outgoing flow of commodity at the ${i}^{\text{th}}$ node, the amount of flow changes with time.

To manage the fluctuation in the amount of flow at intermediate nodes require storage capacity. This phenomenon provides me with the insight to measure the level of commodity which will be stored during this random change in the network. In case of the information flow, we will be able to know about the level of information at particular time according to that we may solve the routing problem of information flow with performance objectives.

Consider a node in computer communication network that receives the packet of information from various sources and must forward them toward their unlimited destination. Generally, the node has a fixed or at least a maximum rate at which it can transmit data. Since the arrival of packet to a node will be quite random, the node will usually have some buffering capacity, allowing the node to temporarily store packets which it can forward immediately.

Robert  discussed about the buffer of data packet at the node during the flow of the packets of information in network. They have mentioned the rate control schemes when the queue size grows infinitely then the buffer space at some nodes may be exhausted. For this reason, some of the packets at these nodes will have to be discarded and later retransmitted, thereby wasting communication resources. Proper buffer management can also help avoid deadlocks due to buffer overflow. Some of them are given as call blocking, packet discarding, packet blocking and packet scheduling. Mishra and Mishra  attempted an economic order quantity model for queued customers with normal life time of inventory flowing in supply chain. Mishra and Singh  worked on an economic order quantity (EOQ) model with quadratic deterioration items allowing shortages for computing the total optimal cost. There are other several researchers Mishra and Mishra  investigated an inventory model for deteriorating items under cobweb phenomenon and permissible delay in payment for fuzzy environment. And Mishra et al.  illustrated the inventory flow in supply chain with deteriorating items for customers in queue and optimized the profit under in fuzzy environment. Jose et al.  have worked on the online measurement of large traffic aggregates on commodity switches such as anomaly detection and traffic engineering and also proposed a practical measurement model based on emerging open flow switches. Gill et al.  have compiled models of information flow through queuing network and the multi-hop flow routing dimension. They also proposed it as convex optimization problem. Thomas et al.  analyzed passive measurement system for network test beds. The ability to capture and process packet-level data is of intrinsic importance in network test beds that offer broad experimental capabilities to researchers. Since we know that commodity may be the energy or oil, gas, fluid etc. Narasimhan  and Farlow  worked on the Fourier’s heat conduction equation wrote that Fourier’s heat

conduction equation i.e. $\frac{\partial }{\partial t}u\left(x,t\right)=k\frac{{\partial }^{2}}{\partial {x}^{2}}u\left(x,t\right)$ continues to constitute the

conceptual foundation on which rests the analysis of many physical, biological, and social systems. This equation of heat flow may be treated as the partial differential equation of the commodity flow which has most important role in this paper. Sifaleras  presented a wide range of problems concerning minimum cost network flow and gave an overview of the classic linear single commodity minimum cost network flow problems and some other closely related cases. Fonoberova and Rebennack  discussed an approach for solving some power system problems by using optimal dynamic flow problems. All parameters of networks are assumed to be dependent of time. They used the algorithm of special dynamic programming techniques based on the time-expanded network method together with classical optimization method. M. Fonoberova and Dmitrii  attempted the dynamic version of the minimum cost multi commodity flow problem. Sisman  illustrated about the determination of the location of the facilities, size of the warehouses and its numbers. Shah and Ward  also defined and developed the measures of lean production. Mishra et al.  have analyzed the traits of investment packages and obtained the equation of demand in terms of safety, security and growth. Mishra and Mishra  have reported on the demand as a power of level of inventory with shape, size parameters and determined the price of an economic order quantity for deteriorating items under perfect competition.

The present paper focuses on the alteration of level of node wise flowing commodity due to random change in the incoming and outgoing flow of commodity at particular time as novel and useful contribution to the existing knowledge. The demand rate of the commodity is in power form of available commodity with shape and size parameter. Here, we seek to measure and compute the level of commodity at any node at a particular time in terms of shape, size parameters and amount of flow of commodity. The cost analysis of incoming and outgoing node wise flowing commodities has been discussed under the different cases of size parameter.

2. Notations

There are following notations and its abbreviation for whole paper

${x}_{ij}$ = Amount of outgoing commodity from ${i}^{\text{th}}$ node to ${j}^{\text{th}}$ , $\forall j\in V$ .

${x}_{ji}$ = Amount of incoming commodity to ${i}^{\text{th}}$ node from ${j}^{\text{th}}$ , $\forall j\in V$ .

${x}_{i}=\sum _{j|\left(i,j\right)\in E}{x}_{ij}-\sum _{j|\left(i,j\right)\in E}{x}_{ji}$ = Amount of flow at ${i}^{\text{th}}$ node.

${I}_{0,t}$ = Level of commodity at time t with zero amount of flow.

${I}_{\epsilon ,t}$ = Level of commodity at time t with amount of flow ( $\in >\text{Capacity}\text{\hspace{0.17em}}$ ).

${I}_{0,t}$ = Level of commodity at initial time with amount of flow x.

${I}_{x,t}\text{\hspace{0.17em}}$ = Level of commodity, at time t and amount of flow at x.

$T$ = Tension at any time on the node (to manage the level of commodity).

$\delta x$ = Small change in amount of flow.

$\alpha$ = Shape parameter which measures the shape of commodity packets i.e. $\alpha >0.$

$\beta$ = Size parameter which measures the size of commodity packets i.e. $\beta >0.$

$\theta$ = Rate of deterioration or evaporation or obsolescence of the commodity.

$\frac{\partial \theta }{\partial t}=\varphi$ = Acceleration of deterioration or evaporation or obsolescence of the commodity.

${D}_{x,\text{\hspace{0.17em}}t}$ = Demand at time t and x amount of flow.

$\rho$ = Traffic intensity of commodity at node i.e. ratio of rate of incoming commodity with rate of outgoing commodity.

$\xi$ = Component of tension in horizontal direction i.e. parallel to the amount of flow axis.

${T}_{1}$ = Tension at any point P at $\left(x,t\right)$ co-ordinate.

${T}_{2}$ = Tension at any point Q at $\left(x+\delta x,t\right)$ co-ordinate.

${C}_{{h}_{ji}}$ = Carrying cost of commodity which is coming from ${j}^{\text{th}}$ nodes to ${i}^{\text{th}}$ node.

${C}_{{O}_{ij}}$ = Ordering cost per unit of commodity which is outgoing from i node to jth nodes.

${C}_{{p}_{ji}}$ = Purchasing cost per unit of commodity by node ${i}^{\text{th}}$ from node ${j}^{\text{th}}$ .

${C}_{{h}_{i}}$ = Holding cost per unit of commodity at ${i}^{\text{th}}$ node.

TVC = Total variance cost of the network system.

${q}_{i}^{*}$ = The optimal inflow quantity at the ${i}^{\text{th}}$ node.

3. Assumptions

The following assumptions have been thoroughly used in the paper.

1) Let the behavior of the intermediate nodes (supplier nodes) are appearing to be change according to the change of the amount of flow. Either it may play the role of source or sink.

2) The number of incoming arc flow and outgoing arc flow are randomly changing with respect to time at an intermediate node.

3) Let the initial amount of flow on a node is zero at any time t. Hence, it implies that level of commodity is zero on it.

4) The level of a commodity at a node increases with the increase of the amount of flow later on it starts decrease with an increase in the amount of flow of commodity.

5) When the amount of flow become more than capacity, the availability level of a commodity, becomes zero. It may be understood as when the demand of the commodity becomes more than the available commodity, in this situation amount of flow becomes more than the capacity i.e. node is now at the level of saturation.

6) The commodity is coming gradually from source and supply to the other nodes. We find that the finite level of commodity is maintained at the supplier node for fulfill the demand.

7) The elastic behavior of level of commodity with respect to the amount of flow. Due to this elasticity, the change in the amount of flow $x$ $x+\delta x$ at time t, node has sophisticated tension to manage the level of the commodity that changes ${I}_{x,t}$ to ${I}_{x+\delta x,t}$ . This tension can be realized on the basis of mainly two forces which are working on level of commodity.

8) The intermediate node manages the flow of commodities to fulfill the demand with its available stock of commodity at an instant of time. This creates a managerial tension which may be realized by two forces that are working on level of commodity. One of them is due to deterioration rate or the rate of evaporation (or obsolescence) that is working as the rate of retardation on level of

commodity i.e. $\frac{\partial {D}_{x,t}}{\partial t}=\alpha {I}_{x,t}^{\beta }$ and another is the cause of the rate of demand of commodity, i.e. $\frac{\partial {D}_{x,t}}{\partial t}=\alpha {I}_{x,t}^{\beta }$ (power form of level of commodity).

9) These forces create a tension on the level of commodity. Let T1 and T2 be tensions, making $\psi \text{and}\vartheta$ angles from the horizontal axis parallel to the amount of flow axis.

10) Actually, this tension is a hypothetical tension and it may be measured as

$T=\rho \delta x\left(\frac{{\partial }^{2}D}{\partial {t}^{2}}±\frac{{\partial }^{2}\theta }{\partial {t}^{2}}\right)$

11) When level of commodity increases, then $T=\rho \delta x\left(\frac{{\partial }^{2}D}{\partial {t}^{2}}-\frac{{\partial }^{2}\theta }{\partial {t}^{2}}\right)$ and when level of commodity decreases then $T=\rho \delta x\left(\frac{{\partial }^{2}D}{\partial {t}^{2}}+\frac{{\partial }^{2}\theta }{\partial {t}^{2}}\right)\text{\hspace{0.17em}}$ .

If we resolve $T$ into two components, we find that the component along “amount of flow” axis direction as below

$-{T}_{1}\mathrm{cos}\psi +{T}_{2}\mathrm{cos}\vartheta =0⇒\xi ={T}_{1}\mathrm{cos}\psi ={T}_{2}\mathrm{cos}\vartheta$

And another component ${T}_{2}\mathrm{sin}\vartheta -{T}_{1}\mathrm{sin}\psi$ = Force for uplift the level of commodity.

The boundary conditions for the level of commodity at a supplier node is assumed as below

1) ${I}_{0,t}=0$ 2) ${I}_{\epsilon ,t}=0,\text{}\epsilon >C=\text{Capacity}$ 3) ${I}_{x,0}=0\text{or}f\left(x\right)4\right)$

${\left(\frac{\partial {I}_{x,t}}{\partial t}\right)}_{t=0}=0\text{or}g\left(x\right)$

4. Development of Mathematical Model and Its Solution

We have obtained the following partial differential equation of the level of commodity at a particular node

$\partial x\frac{\alpha \beta \rho }{\xi }{I}_{x,t}^{\beta -1}\frac{\partial {i}_{x,t}}{\partial t}±\partial x\frac{\rho }{\xi }\frac{\partial \varphi }{\partial t}=\frac{\left({T}_{2}\mathrm{sin}\vartheta -{T}_{1}\mathrm{sin}\psi \right)}{\xi }$

${I}_{x,t}^{\beta -1}\frac{\partial {I}_{x,t}}{\partial t}±\frac{1}{\alpha \beta }\frac{\partial \varphi }{\partial t}=\left(\frac{\xi }{\alpha \beta \rho }\right)\frac{\left[{\left(\frac{\partial {I}_{x,t}}{\partial x}\right)}_{x+\delta x}-{\left(\frac{\partial {I}_{x,t}}{\partial x}\right)}_{x}\right]}{\partial x}$ (1)

Let $\frac{\xi }{\alpha \beta \rho }={C}^{2}$ and assume that deterioration rate $\varphi =\frac{\partial \theta }{\partial t}$ = Constant i.e. $\frac{\partial \varphi }{\partial t}=0$

$\frac{\partial {I}_{x,t}}{\partial t}=\frac{{C}^{2}}{{I}_{x,t}^{\beta -1}}\frac{{\partial }^{2}{I}_{x,t}}{\partial {x}^{2}}$ (2)

4.1. Case I When β = 1 Then, from Equation (2) We Get

$\frac{\partial {I}_{x,t}}{\partial t}={C}^{2}\frac{{\partial }^{2}{I}_{x,t}}{\partial {x}^{2}}$ , Let the solution of the differential equation ${I}_{x,t}\text{\hspace{0.17em}}=X\left(x\right)\cdot T\left(t\right)$

and applying method of separation of variable to solve the partial differential equation, we get

$\frac{1}{X}\frac{{\text{d}}^{2}X}{\text{d}{x}^{2}}=\frac{1}{{}^{{C}^{2}T}}\frac{\text{d}T}{\text{d}t}=k$

There are three possible conditions for $k.$

1) If $k=0$ we get, ${I}_{x,t}\text{\hspace{0.17em}}=Ax+B$ where $A$ and $B$ are constants. It implies that level of commodity is independent of time only linear function of amount of flow. This is suitable for the steady state flow but not for the transient flow model.

2) For $k<0$ , let $k=-{p}^{2}$ , we get ${I}_{x,t}\text{\hspace{0.17em}}=\left(a\mathrm{cos}px+b\mathrm{sin}px\right){\text{e}}^{-{p}^{2}{c}^{2}t}$ , where $a$ , $b$ and $p$ are constants. This implies that the level of commodity is periodic with amount of flow as we as converges at zero when $t\to \infty$ . This model is generally suitable for measurement of amount of energy level at supplier node.

3) For $k>0=-{p}^{2}$ , we get ${I}_{x,t}\text{\hspace{0.17em}}=\left(d{\text{e}}^{px}+f{\text{e}}^{-px}\right){\text{e}}^{{p}^{2}{c}^{2}t}$ , where $d$ and $f$ are constants. It shows that when $t\to \infty$ , the level of commodity also tend to infinite. It is not valid for physical model.

4.2. Case II When $\beta \ne 1$ and Let $\varphi =\frac{\partial \theta }{\partial t}$ = Constant

$\frac{\partial {I}_{x,t}}{\partial t}=\frac{{C}^{2}}{{I}_{x,t}^{\beta -1}}\frac{{\partial }^{2}{I}_{X,T}}{\partial {x}^{2}}$ , ${I}_{x,t}={\left\{{C}_{2}x\right\}}^{\left(\frac{2}{1-\beta }\right)}{\left\{\left(\beta -1\right)\left(\frac{2\xi \left(1+\beta \right)}{\alpha \beta \rho {\left(1-\beta \right)}^{2}}{\left({C}_{2}\right)}^{2}t\right)\right\}}^{\frac{1}{\beta -1}}$

${I}_{x,t}={\left(\frac{2\xi \left(1+\beta \right)}{\alpha \beta \rho \left(1-\beta \right)}\right)}^{\frac{1}{\beta -1}}{\left(\frac{{x}^{2}}{t}\right)}^{\frac{1}{1-\beta }}$ (3)

Since, ${I}_{x,t}={\left(\frac{2\xi \left(1+\beta \right)}{\alpha \beta \rho \left(1-\beta \right)}\right)}^{\frac{1}{\beta -1}}{\left(\frac{{x}^{2}}{t}\right)}^{\frac{1}{1-\beta }}$ , so demand rate can be expressed as

below.

$\frac{\partial {D}_{x,t}}{\partial t}=\alpha {\left\{\frac{2\xi \left(1+\beta \right)}{\alpha \beta \rho \left(1-\beta \right)}\right\}}^{\frac{\beta }{\beta -1}}{\left(\frac{{x}^{2}}{t}\right)}^{\frac{\beta }{1-\beta }}$

It implies that

$\frac{\partial {D}_{x,t}}{\partial t}=\alpha {\left\{\frac{2\xi \left(1+\beta \right)}{\alpha \beta \rho \left(1-\beta \right)}\right\}}^{\frac{\beta }{\beta -1}}{\left(\frac{{x}^{2}}{t}\right)}^{\frac{\beta }{1-\beta }}{x}^{\frac{2\beta }{1-\beta }}\cdot {t}^{\frac{1-2\beta }{1-\beta }}+C$

Let the initial condition be ${D}_{0,0}=0$ then C = 0,

Hence, ${D}_{x,t}=\alpha {\left\{\frac{2\xi \left(1+\beta \right)}{\alpha \beta \rho \left(1-\beta \right)}\right\}}^{\frac{\beta }{\beta -1}}\frac{\left(1-\beta \right)}{\left(1-2\beta \right)}{x}^{\frac{2\beta }{1-\beta }}\cdot {t}^{\frac{1-2\beta }{1-\beta }}$ (4)

The above expression shows the relation between demand and amount of flow.

The level of commodity, for particular amount of flow at sink as well as for same amount of flow at the source, may be same.

4.2.1. Case-A (If β < 1)

If $\beta <1$ then ${I}_{x,t}={\left(\frac{\alpha \beta \rho \left(1-\beta \right)}{2\xi \left(1+\beta \right)}\right)}^{\frac{1}{1-\beta }}{\left(\frac{{x}^{2}}{t}\right)}^{\frac{1}{1-\beta }}$

Economic Inflow quantity (Optimal Inflow quantity)

Since total incoming amount of flow is ${q}_{i}=\sum _{j|\left(i,j\right)\in E}{x}_{ji}$ , it takes $t$ time in

supply to other nodes. Here, we consider the instantaneous fulfillment of commodity at a node

Hence, $\frac{{C}_{{h}_{ji}}}{2}×{q}_{i}×t={C}_{{o}_{ij}}×\frac{t×\frac{\partial {D}_{x,t}}{\partial t}}{{q}_{i}}$

If $\beta <1$ then ${I}_{x,t}={\left(\frac{\alpha \beta \rho \left(1-\beta \right)}{2\xi \left(1+\beta \right)}\right)}^{\frac{1}{1-\beta }}{\left(\frac{{x}^{2}}{t}\right)}^{\frac{1}{1-\beta }}$

${q}_{i}^{*}={q}_{i}=\sqrt{\frac{{C}_{{O}_{ij}}}{{C}_{{h}_{ji}}}\frac{\partial {D}_{x,t}}{\partial t}}={x}^{\frac{\beta }{1-\beta }}\sqrt{\frac{\alpha {C}_{{O}_{ij}}}{{C}_{{h}_{ji}}}{\left\{\frac{\left(1-\beta \right)\alpha \beta \rho }{2\left(\beta +1\right)\xi t}\right\}}^{\frac{\beta }{1-\beta }}}$

${q}_{i}^{*}={x}^{\frac{\beta }{1-\beta }}\sqrt{\frac{\alpha {C}_{{O}_{ij}}}{{C}_{{h}_{ji}}}{\left\{\frac{\left(1-\beta \right)\alpha \beta \rho }{2\left(\beta +1\right)\xi t}\right\}}^{\frac{\beta }{1-\beta }}}$ (5)

Cost analysis at a particular node of a network

We here minimize TVC as

Minimize $TVC=\sum _{i\in V}\sum _{j|\left(i,j\right)\in E}C{o}_{ij}{x}_{ij}+\sum _{i}{C}_{{h}_{i}}{I}_{{x}_{i},t}+\sum _{j|\left(i,j\right)\in E}\sum _{i\in V}{C}_{{p}_{ji}}{x}_{ji}$

Subject to constraints

$\sum _{j|\left(i,j\right)\in E}{x}_{ij}-\sum _{j|\left(i,j\right)\in E}{x}_{ji}={x}_{i}$

${I}_{{x}_{i},t}={\left\{\frac{\left(1-\beta \right)\alpha \beta \rho }{2\left(\beta +1\right)\xi }\right\}}^{\frac{1}{1-\beta }}{\left(\frac{{x}^{2}}{t}\right)}^{\frac{1}{1-\beta }}$

$\sum _{i\in V}{x}_{i}=0,\text{\hspace{0.17em}}{x}_{ij}\ge 0{\text{andx}}_{ji}\ge 0,\text{\hspace{0.17em}}0\le {x}_{i}\le \epsilon$ .

4.2.2. Case-IIB (If β > 1)

If $\beta >1$ then ${I}_{x,t}={\left\{-\frac{\left(\beta -1\right)\alpha \beta \rho }{2\left(\beta +1\right)\xi }\right\}}^{\frac{1}{\beta -1}}{\left(\frac{{x}^{2}}{t}\right)}^{\frac{1}{1-\beta }}$

${I}_{x,t}={\left\{\frac{\left(\beta -1\right)\alpha \beta \rho }{2\left(\beta +1\right)\xi }\right\}}^{\frac{1}{\beta -1}}{\left(\frac{1}{x}\right)}^{\frac{2}{\beta -1}}{\text{e}}^{\left[\frac{2\mathrm{log}\left(i{t}^{\frac{1}{2}}\right)}{\beta -1}\right]}$ (6)

Under the above condition, from the expression (6) it is clear that we get complex values which are not possible in the practical situations. But in the definite range of size parameter $\beta$ i.e. $1<\beta <1.5$ we get real values of Ix,t as

${I}_{x,t}={\left\{\frac{\left(\beta -1\right)\alpha \beta \rho }{2\left(\beta +1\right)\xi }\right\}}^{\frac{1}{\beta -1}}{\left(\frac{t}{{x}^{2}}\right)}^{\frac{1}{\beta -1}}$

${q}^{*}={q}_{i}=\sqrt{\frac{{C}_{{O}_{ij}}}{{C}_{{h}_{ji}}}\frac{\partial {D}_{x,t}}{\partial t}}=\frac{1}{{x}^{\frac{\beta }{\beta -1}}}\sqrt{\alpha {\left\{\frac{\left(\beta -1\right)\alpha \beta \rho }{2\left(\beta +1\right)\xi }\right\}}^{\frac{\beta }{\beta -1}}\frac{{C}_{{O}_{ij}}}{{C}_{{h}_{ji}}}}$

Let $\beta =1+\frac{1}{2n},\text{where}n=1,2,3\cdots .$

Consequently, we get the series solution

${I}_{x,t}=\sum _{n=1}^{\infty }{\left\{\frac{\alpha \left(2n+1\right)\rho }{4n\left(1+4n\right)\xi }\right\}}^{2n}{\left(\frac{t}{{x}^{2}}\right)}^{2n}\text{,where}\text{}n=1,2,3\cdots$

Economic inflow quantity (optimal inflow quantity)

${q}^{*}={q}_{i}=\sqrt{\frac{{C}_{{O}_{ij}}}{{C}_{{h}_{ji}}}\frac{\partial {D}_{x,t}}{\partial t}}=\frac{1}{{x}^{\frac{\beta }{\beta -1}}}\sqrt{\alpha {\left\{\frac{\left(\beta -1\right)\alpha \beta \rho t}{2\left(\beta +1\right)\xi }\right\}}^{\frac{\beta }{\beta -1}}\frac{{C}_{{O}_{ij}}}{{C}_{{h}_{ji}}}}$

Cost analysis at a particular node of a network

In this section, we minimize the total variance cost as

Minimize

$TVC=\sum _{i}\sum _{j|\left(i,j\right)\in E}{C}_{{O}_{ij}}\text{ }{x}_{ij}+\sum _{i}{C}_{{h}_{i}}{I}_{{x}_{i},t}+\sum _{j|\left(i,j\right)\in E}\sum _{i}{C}_{{p}_{ji}}\text{ }{x}_{ji}$

Subject to the constraints

$\sum _{j|\left(i,j\right)\in E}{x}_{ij}-\sum _{j|\left(i,j\right)\in E}{x}_{ji}={x}_{i}$

${I}_{x,t}={\left\{\frac{\left(\beta -1\right)\alpha \beta \rho }{2\left(\beta +1\right)\xi }\right\}}^{\frac{1}{\beta -1}}{\left(\frac{t}{{x}^{2}}\right)}^{\frac{1}{\beta -1}}$

$\sum _{i=1}{x}_{i}=0,\text{\hspace{0.17em}}{x}_{ij}\ge 0\text{and}{x}_{ji}\ge 0,\text{\hspace{0.17em}}0\le {x}_{i}\le \epsilon$

5. Sensitivity Analysis

The sensitivity analysis seeks to study, the variational approach of various parameters involved in the model which is used as an important instrument to judge the validity of the model for its future application. The following observations are worth mentioning to discuss the sensitivity analysis of the model. Matlab has been used for computing the values and drawing the graphs in this section. Tables and graphs are arranged under Appendix at last.

5.1. Case-I

5.1.1. Observations from the Tables

From Table A1, an observation is drawn that if we increase the value of shape parameter of commodity then the level of commodity decreases but if we simultaneously decrease the size parameter the level of commodity gradually increases. Secondly, we may observe that if only stress parameter $\xi$ increases and others remain constant then level of commodity decreases.

From the Table A2, we conclude that when size parameter $\beta$ and stress parameter $\xi$ time $t$ remain constant then increase in the amount of flow becomes a cause of increase in the level of commodity at the node. And if all $\beta$ , $\xi$ , $t$ are increased simultaneously, then it shows the decrease in level of commodity at node. It is also observed that the increment in the size parameter $\beta$ of commodity decreases the level of commodity but increase in amount of flow, increases the level of commodity.

The observations are drawn from Table A3 to Table A6(c) showing that the optimal inflow quantity increases when 1) the amount of flow increases as well as the size parameter decreases; 2) the increment occurs in the amount of flow as well as in time; 3) the increment occurs in the value of stress parameter as well as in the amount of flow; 4) the increment occurs in shape parameter as well as in stress parameter; 5) the increment occurs in the traffic intensity as well as the value of stress parameter.

But the optimal inflow quantity decreases when 1) the value of stress parameter increases for fixed amount of flow; 2) the increment occurs in the size parameter as well as amount of flow.

5.1.2. Observations from Graphs

From the Figure A1, we may report that the increase in the amount of flow of commodity at node increases the level of commodity but after some time it decreases very rapidly. It is also observed that when the amount of flow reaches its maximum value then the level of commodity gets small which is true for ideal intermediate nodes of the network.

The Figure A2 illustrates about the effect of traffic intensity of commodity. Initially, when it increases then corresponding level of commodity also increases but after some time duration it may be cause of decrease of level of commodity at the node.

In the Figure A3, when the size parameter remains constant then we observe that with the increase in the traffic intensity of the commodity, the level of commodity also increases. It is seen that traffic intensity increases till certain level and after that it starts decreasing but level of commodity does not decrease, it gradually increases.

From Figure A4, it is clear that at a particular time when the node has constant value of stress parameter then the level of commodity previously increases with increase in the size parameter but after certain limit it begins decreasing in the level of commodity with increase in the size parameter.

5.2. Case-II

5.2.1. Observations from the Tables

From Table A7(a), it shows that if $\beta$ lies from 1 to 1.25 and traffic intensity remains constant then the level of commodity becomes zero. After this range, increase in the amount of flow is reported but stress parameter remains constant then the level of commodity decreases.

From Table A8(b) and Table A9 we conclude that if we increase in the traffic intensity for the size parameter lies between 1 and 1.5, then the level of commodity increases. But if we increase in the amount of flow then it has decreasing pattern of level of commodity.

From the Table A9, Table A10(a) and Table A10(b), we observe that the optimal inflow quantity increases if 1) the traffic intensity as well as size parameter, stress parameter increases; 2) the value of size parameter and the stress parameter decrease as well as the amount of flow and traffic intensity increase; 3) the amount of flow increases but the optimal inflow quantity decreases when the amount of flow, size parameter and stress parameter increase.

5.2.2. Observations from Figures

In the Figure A5, it shows that for beta (size parameter) lies from 1 to1.25, there is no change in the level of commodity but when size parameter increases from 1.25 to 1.5, the level of commodity also increases.

Figure A6 shows that change in the amount of flow from 0 to 1.25, there occurs no change in the level of commodity. But from change in amount of flow 1.25 to 1.3, a very small change in the level of commodity takes place and after change 1.3 to 1.4, in the amount of flow the level of commodity also linearly increases but change in the amount of flow from 1.4 to 1.5 the level of commodity very rapidly increases. We may conclude from the observation that the level of commodity increases linearly with amount of flow. Moreover, if certain amount of flow linearity increases then this level of commodity increases rapidly.

6. Conclusion

Upon deeply analyzing the node wise level of commodity in the supply chain and total minimum variance cost of the network, we have been able to draw close observations under the sensitivity analysis of the model. By which we conclude that this model is very sensitive about the size parameter of the system which proves the rationale of the value of the flow at any node as an important performance measure of the system. On the basis of level of commodity, we have computed the economic inflow quantities (optimal inflow quantity), and obtained the total minimum variance cost of the network under the subject to the conditions with non-negative constraints. These results have a vast spectrum of application in various fields of network such as network of communication, supply chain network and commodity (gas, oil and energy etc.) flow network. In order to control the flow of commodity (overflow on buffer stock), this model is very useful. This model paves the path to measure the level of commodity during flow in the random network. To measure and quantify the level of commodity at any node under uncertain environment of flow of commodity and to find the effect of the price on the flow of commodity would be our future research.

Cite this paper

Mishra, P.P., Mishra, S.S., Yadav, S.K., Singh, R.S. and Kumar, R. (2017) Quantification of Node Wise Commodity in Supply Chain and Its Cost Analysis. American Journal of Operations Research, 7, 64-82. http://dx.doi.org/10.4236/ajor.2017.71005

References

1. 1. Ahuja, R.K., Goldberg, A.V. Orlin, J.B. and Tarjan, R.E. (1992) Finding Minimum-Cost Flows by Double Scaling, Mathematical Programming, 53, 243-266. https://doi.org/10.1007/BF01585705

2. 2. Ahuja, R.K., Magnanti, T.L. and Orlin, J.B. (1993) Network Flows: Theory, Algorithms and Applications. Prentice Hall, Englewood Cliffss, NJ.

3. 3. Bertsekas, D. and Gallager, R. (1992) Data Networks. 2nd Edition, Prentice Hall, Englewood Cliffs, New Jersey.

4. 4. Mishra, S.S. and Mishra, P.P. (2012) Phase Wise Supply Chain Model of EOQ with Normal Life Time for Queued Customers: A Computational Approach. American Journal of Operation Research, 2, 296-307. https://doi.org/10.4236/ajor.2012.23036

5. 5. Mishra, S.S. and Singh, P.K. (2011) Computational Approach to an Economic Order Quantity Model with Power form Stock Dependent Demand and Cubic Deterioration. American Journal of Operation Research, 1, 5-13.

6. 6. Mishra, S.S. and Mishra, P.P. (2011) A (Q, R) Model for Fuzzyfied Deterioration under Cobweb Phenomenon and Permissible Delay in Payment. Computers and Mathematics with Applications, 61, 921-932. https://doi.org/10.1016/j.camwa.2010.12.040

7. 7. Mishra, S.S., Yadav, S.K. and Rawat, S. (2015) Inventory Flow in Supply Chain with Deteriorating Items for Customer in Queue: Computation of Profit Optimization in Fuzzy Environment. American Journal of Operations Research, 5, 1-8.

8. 8. Jose, L., Yu, M. and Rexford, J. (2011) Online Measurement of Large Traffic Aggregates on Commodity Switches. Princeton University, Princeton.

9. 9. Gill, P., et al. (2008) Network Information Flow in Network of Queues. IEEE International Symposium on Modeling, Analysis and Simulation of Computers and Telecommunication Systems, MASCOTS 2008, Baltimore, 8-10 September 2008, 1-8.

10. 10. Thomas, C., et al. (2012) A Passive Measurement System for Network Test Beds. Book on Lecture Notes of the Institute for Computer Sciences, Social Informatics a Telecommunications Engineering, 44, 130-145. https://doi.org/10.1007/978-3-642-35576-9_14

11. 11. Narasimhan, T.N. (1999) Fourier’s Heat Conduction Equation: History, Influence, and Connections. Reviews of Geophysics, 37, 151-172. https://doi.org/10.1029/1998RG900006

12. 12. Farlow, S.J. (1993) Partial Differential Equations for Scientists and Engineers. Dover, Mineola.

13. 13. Sifaleras, A. (2013) Minimum Cost Network Flows: Problems, Algorithms and Software. Yugoslav Journal of Operation Research, 23, 3-17. https://doi.org/10.2298/YJOR121120001S

14. 14. Fonoberova, M. and Rebennack, F.S. (2010) Alogorithms for Finding Optimal Flows in Dynamic Networks. In: Rebennack, S., Pardalos, P.M., Pereira, M.V.F. and Iliadis, N.A., Eds., Handbook of Power System II, Springer, Berlin, 31-54. https://doi.org/10.1007/978-3-642-12686-4_2

15. 15. Fonoberova, M. and Lozovanu, D.D. (2005) The Minimum Cost Multicommodity Flow Problem in Dynamic Networks and an Algorithm for Its Solving. Computer Science Journal of Moldova, 13, 29-36.

16. 16. Sisman, B. (2012) Supply Chain Network Design Considering Customer Service Level. Proceeding of the International Conference on Industrial Engineering and Operations Management, Istanbul, 3-6 July 2012, 2250-2258.

17. 17. Shah, R. and Ward, P.T. (2009) Defining and Developing Measures of Lean Production. The TQM Journal, 21, 127-142.

18. 18. Mishra, S.S., Mishra, P.P. and Sharma, S.K. (2012) Trait Analysis of Investment Packages as Using Computational Technique: A Case Study of Insurance Companies. ESMSJ, 2, 22-28.

19. 19. Mishra, S.S. and Mishra, P.P. (2008) Price Determination for an Economic Order Quantity (EOQ) Model for Deteriorating Items under Perfect Competition. International Journal of Computer and Mathematics with Application, 56, 1082-1101. https://doi.org/10.1016/j.camwa.2008.01.034

Appendix 1

Table A1. $\left({I}_{x,t}\text{\hspace{0.17em}}\text{verses}\text{\hspace{0.17em}}\beta ,\alpha ,\rho ,\xi \right)$ .

Table A2. $\left({I}_{x,t}\text{\hspace{0.17em}}\text{verses}\text{\hspace{0.17em}}\beta ,x,t,\text{and}\xi \right)$ .

Table A3. $\left({q}^{*}\text{\hspace{0.17em}}\text{verses}\text{\hspace{0.17em}}\beta ,\alpha ,\rho ,x,t\text{and}\xi \right)$ Ch = 50 cost unit, Co = 100 cost unit).

Table A4. $\left({q}^{*}\text{\hspace{0.17em}}\text{verses}\text{\hspace{0.17em}}\beta ,\alpha ,\rho ,x,t\text{and}\xi \right)$ .

Table A5. $\left({q}^{*}\text{\hspace{0.17em}}\text{verses}\text{\hspace{0.17em}}\beta ,\alpha ,\rho ,x,t\text{and}\xi \right)$ .

(a) (b) (c)

Table A6. (a)-(c) $\left({q}^{*}\text{\hspace{0.17em}}\text{verses}\text{\hspace{0.17em}}\beta ,\alpha ,\rho ,x,t\text{and}\xi \right)$ .

Table A7. $\left({I}_{x,t}\text{ }\text{verses}\text{\hspace{0.17em}}\beta ,x,t,\xi \text{}\text{ }\text{and}\alpha =100\right)$ .

(a) (b)

Table A8. (a)-(b) $\left(\text{\hspace{0.17em}}{I}_{x,t}\text{ }\text{verses}\text{\hspace{0.17em}}\beta ,x,t,\rho ,\xi \text{ }\right)$ .

Table A9. $\left({q}^{*}\text{ }\text{verses}\text{\hspace{0.17em}}\beta ,\alpha ,\rho ,x,t,\xi \right)$ .

Table A10. (a) $\left({q}^{*}\text{ }\text{verses}\text{\hspace{0.17em}}\text{}\beta ,\alpha ,\rho ,x,t,\xi \text{ }\right)$ . (b) $\left({q}^{*}\text{ }\text{verses}\text{\hspace{0.17em}}\beta ,\alpha ,\rho ,x,t,\xi \right)$ .

Case I:

Figure A1. Level of commodity, traffic intensity and time.

Figure A2. Rho (traffic intensity), time and level of commodity when size parameter $\beta$ varies.

Figure A3. Traffic intensity (rho), level of commodity and time when size parameter $\beta$ is constant.

Figure A4. Size parameter (beta) and level of commodity.

Case II:

Figure A5. Level of commodity and size parameter (beta = b).

Figure A6. Amount of flow and level of commodity.

Appendix 2

Some proofs are given below

Case II When $\beta \ne 1$ and Let $\varnothing =\frac{\partial \theta }{\partial t}=\text{Constant}$

$\frac{\partial {I}_{x,t}}{\partial t}=\frac{{C}^{2}}{{I}_{x,t}^{\beta -1}}\frac{{\partial }^{2}{I}_{x,t}}{\partial {x}^{2}}$ , ${I}_{x,t}={\left\{{C}_{2}x\right\}}^{\left(\frac{2}{1-\beta }\right)}{\left\{\left(\beta -1\right)\left(\frac{2\xi \left(1+\beta \right)}{\alpha \beta \rho {\left(1-\beta \right)}^{2}}{\left({C}_{2}\right)}^{2}t\right)\right\}}^{\frac{1}{\beta -1}}$

${I}_{x,t}={\left(\frac{2\xi \left(1+\beta \right)}{\alpha \beta \rho \left(1-\beta \right)}\right)}^{\frac{1}{\beta -1}}{\left(\frac{{x}^{2}}{t}\right)}^{\frac{1}{1-\beta }}$

Since, ${I}_{x,t}={\left(\frac{2\xi \left(1+\beta \right)}{\alpha \beta \rho \left(1-\beta \right)}\right)}^{\frac{1}{\beta -1}}{\left(\frac{{x}^{2}}{t}\right)}^{\frac{1}{1-\beta }}$ , so demand rate can be expressed as

below.

$\frac{\partial {D}_{x,t}}{\partial t}=\alpha {\left\{\frac{2\xi \left(1+\beta \right)}{\alpha \beta \rho \left(1-\beta \right)}\right\}}^{\frac{\beta }{\beta -1}}{\left(\frac{{x}^{2}}{t}\right)}^{\frac{\beta }{1-\beta }}$

It implies that

$\begin{array}{l}{D}_{x,t}=\alpha {\left\{\frac{2\xi \left(1+\beta \right)}{\alpha \beta \rho \left(1-\beta \right)}\right\}}^{\frac{\beta }{\beta -1}}{\left(\frac{{x}^{2}}{t}\right)}^{\frac{\beta }{1-\beta }}{x}^{\frac{2\beta }{1-\beta }}\cdot {t}^{\frac{1-2\beta }{1-\beta }}+C\\ ⇒\end{array}$

Let the initial condition be ${D}_{0,0}=0$ then $C=0$ ,

Hence, ${D}_{x,t}=\alpha {\left\{\frac{2\xi \left(1+\beta \right)}{\alpha \beta \rho \left(1-\beta \right)}\right\}}^{\frac{\beta }{\beta -1}}\frac{\left(1-\beta \right)}{\left(1-2\beta \right)}{x}^{\frac{2\beta }{1-\beta }}\cdot {t}^{\frac{1-2\beta }{1-\beta }}$

$\frac{\partial {I}_{x,t}}{\partial t}=\frac{{C}^{2}}{{I}_{x,t}^{\beta -1}}\frac{{\partial }^{2}{I}_{x,t}}{\partial {x}^{2}}$ , let ${I}_{x,t}=X\left(x\right)\cdot T\left(t\right)$ be a solution and use the method of

separation of variables to solve the partial differential equation. We get

${X}^{\beta }\frac{\text{d}T}{\text{d}t}={C}^{2}{T}^{2-\beta }\frac{{\text{d}}^{2}X}{\text{d}{x}^{2}}$ and $\frac{1}{{C}^{2}{T}^{2-\beta }}\frac{\text{d}T}{\text{d}t}=\frac{1}{{X}^{\beta }}\frac{{\text{d}}^{2}X}{\text{d}{x}^{2}}=l$ which in turn gives

us

$T={\left[\left(\beta -1\right)\left({C}^{2}lt+{c}_{1}\right)\right]}^{\frac{1}{\beta -1}}$

$\frac{1}{{X}^{\beta }}\frac{{\text{d}}^{2}X}{\text{d}{x}^{2}}=l$

If $l=\frac{\left(\beta +1\right)}{{\left(1-\beta \right)}^{2}}{\left({C}_{2}\right)}^{2}$ then $X={\left\{{C}_{1}+{C}_{2}x\right\}}^{\left(\frac{2}{1-\beta }\right)}$ is a solution to differential

equation.

Further,

$\frac{1}{{C}^{2}{T}^{2-\beta }}\frac{\text{d}T}{\text{d}t}=l⇒T={\left[\left(\beta -1\right)\left({C}^{2}lt+{C}_{3}\right)\right]}^{\frac{1}{\beta -1}}$

${I}_{x,t}={\left\{{C}_{1}+{C}_{2}x\right\}}^{\left(\frac{2}{1-\beta }\right)}{\left[\left(\beta -1\right)\left({C}^{2}lt+{C}_{3}\right)\right]}^{\frac{1}{\beta -1}}$

${I}_{0,t}={\left\{{C}_{1}\right\}}^{\left(\frac{2}{1-\beta }\right)}{\left[\left(\beta -1\right)\left({C}^{2}lt+{C}_{3}\right)\right]}^{\frac{1}{\beta -1}}=0⇒{C}_{1}=0$

${I}_{x,0}={\left\{{C}_{2}x\right\}}^{\left(\frac{2}{1-\beta }\right)}{\left[\left(\beta -1\right)\left({C}_{3}\right)\right]}^{\frac{1}{\beta -1}}=0⇒{C}_{3}=0$

${I}_{x,t}={\left\{{C}_{2}x\right\}}^{\left(\frac{2}{1-\beta }\right)}{\left[\left(\beta -1\right)\frac{2\xi \left(\beta +1\right)}{\alpha \beta \rho {\left(\beta -1\right)}^{2}}{\left({C}_{2}\right)}^{2}t\right]}^{\frac{1}{\beta -1}}$

${I}_{x,t}={\left[\frac{2\xi \left(\beta +1\right)}{\alpha \beta \rho \left(1-\beta \right)}\right]}^{\frac{1}{\beta -1}}{\left(\frac{{x}^{2}}{t}\right)}^{\frac{1}{1-\beta }}$

Since ${I}_{x,t}={\left[\frac{2\xi \left(\beta +1\right)}{\alpha \beta \rho \left(1-\beta \right)}\right]}^{\frac{1}{\beta -1}}{\left(\frac{{x}^{2}}{t}\right)}^{\frac{1}{1-\beta }}$ so demand rate can be expressed as

below.

$\frac{\partial {D}_{x,t}}{\partial t}=\alpha {\left\{\frac{2\xi \left(1+\beta \right)}{\alpha \beta \rho \left(1-\beta \right)}\right\}}^{\frac{\beta }{\beta -1}}{\left(\frac{{x}^{2}}{t}\right)}^{\frac{\beta }{1-\beta }}$

${D}_{x,t}=\alpha {\left\{\frac{2\xi \left(1+\beta \right)}{\alpha \beta \rho \left(1-\beta \right)}\right\}}^{\frac{\beta }{\beta -1}}{\left(\frac{{x}^{2}}{t}\right)}^{\frac{\beta }{1-\beta }}{x}^{\frac{2\beta }{1-\beta }}\cdot {t}^{\frac{1-2\beta }{1-\beta }}+C$

Let the initial condition be ${D}_{0,0}=0$ then C = 0,

${D}_{x,t}=\alpha {\left\{\frac{2\xi \left(1+\beta \right)}{\alpha \beta \rho \left(1-\beta \right)}\right\}}^{\frac{\beta }{\beta -1}}{\left(\frac{{x}^{2}}{t}\right)}^{\frac{\beta }{1-\beta }}{x}^{\frac{2\beta }{1-\beta }}\cdot {t}^{\frac{1-2\beta }{1-\beta }}$

When $\beta \ne 0$ then there are following two cases arise as

Case-IIB $\left(\text{If}\text{\hspace{0.17em}}\beta >1\right)$ , If $\beta >1$ then

${I}_{x,t}={\left\{\frac{\left(\beta -1\right)\alpha \beta \rho }{2\left(\beta +1\right)\xi }\right\}}^{\frac{1}{\beta -1}}{\left(\frac{{x}^{2}}{t}\right)}^{\frac{1}{1-\beta }}{\left(-1\right)}^{\frac{1}{\beta -1}}⇒{I}_{x,t}={\left\{\frac{\left(\beta -1\right)\alpha \beta \rho }{2\left(\beta +1\right)\xi }\right\}}^{\frac{1}{\beta -1}}{\left(\frac{1}{x}\right)}^{\frac{2}{\beta -1}}{\text{e}}^{\left[\frac{2\mathrm{log}\left(i{t}^{2}\right)}{\beta -1}\right]}$