*Article* **Effects of Carbon-Emission and Setup Cost Reduction in a Sustainable Electrical Energy Supply Chain Inventory System**

#### **Umakanta Mishra 1, Jei-Zheng Wu 1,\* and Anthony Shun Fung Chiu <sup>2</sup>**


Received: 15 February 2019; Accepted: 28 March 2019; Published: 29 March 2019

**Abstract:** This article develops a sustainable electricity supply chain mathematical model that assumes linear price-dependent customer demands where the price is a decision variable under setup cost and carbon emission. The sustainable electrical supply chain system contained: (a) power generation; (b) transmission substations; (c) distribution substations; and (d) customer. The production rates depend on the demand rate, and demand for electricity by the customers is dependent on the price of electricity where the electrical energy was generated and transmitted through multiple substations to customers. Moreover, we considered that the capacities of transmission rates, power generation, and distances in between two stations are associated with the distribution costs and transmission cost. Here, we used the theory of inventory to develop a new model and suggested a procedure to deduce an optimal solution for this model. Finally, a numerical example and sensitivity analysis are employed to illustrate the present study and with managerial insights.

**Keywords:** sustainable electrical energy supply chain; inventory; price-dependent demand; transmission and distribution costs; carbon-emission

#### **1. Introduction**

#### *1.1. Background of the Research*

Generally, a supply chain inventory model defines the retailer-customer relationship. Setup cost plays a crucial role in today's supply chain inventory management system for shipment of items on time. The setup procedure is not evaluated as a fixed/known constraint, but needs to be considered at a time of minimizing waste, satisfying deadlines, productivity, and elaborating resource utilization. To minimize the total cost investment, the manufacturer needs to reduce setup costs. In a supply chain inventory system, the capacity of electrical energy has the same process for determining the order quantity and total profit. In an electrical supply chain inventory system, electricity generated in a power generation plant will be transmitted through a transmission line and distribution substation to the customers whose demand is influenced by the price of electricity in order to maximize the profit.

The consumption of electricity energy has rapidly increased. The total electricity consumption worldwide in 2015 was greater than in 1980. Electricity consumption for industrial, commercial, residential and transportation sectors has been enhanced since 2015. This is due to increasing climate changes upon home electricity use at a rate of more than 100 million kWh/day, due to the need for air-conditioning (mainly in summer time) in homes. The electricity price may have influenced

the demand for electricity. The price of electricity is set to reduce the demand. Generally, the price is usually the maximum for commercial customers and the residential sector because of the cost to distribute electricity supply to them. Furthermore, an increase of global greenhouse gas carbon emissions (GGCE) in all countries of world is farther than ever from reaching the goals of the Paris Agreement, according to the new United Nations report, released in 2018, prior to a meeting of officials and environmental experts from all around of the world in Katowice, Poland, for climate negotiations. A few years ago at the environmental negotiations in Paris all countries agreed to bring down GGCE sufficient to maintain global warming below 2 degrees Celsius, and under 1.5 degrees. A report from the UN's top climate panel published earlier this decrease found that an additional 0.5 degrees of warming would have forceful and unfortunate effects on the environment, allowing for more empirical support for those countries pressing for more challenging targets during the negotiations at COP24. But to maintain temperature increases under 1.5 degrees, global carbon emissions are required to peak by 2020. In between 2014 to 2016, global carbon emissions stayed comparatively flat, and negotiators hope that the decade long trend of increasing carbon emissions is about to reverse. But after few years of stagnation, global carbon emissions grew again in 2017, reaching a record 53.4 gigatons of carbon emission equivalent. To meet the 1.5-degree mark, the world requires to jointly bring emissions down in the next 12 years. Therefore, governments all over the world should implement cap-and trade regulations policies to keep the emissions low. For reducing carbon emissions, production companies can monitor and enhance the emission performance of their products during their life-cycle stages. The carbon-emission assessment provides a possible mechanism to serve companies with some emission reduction. During the production process, the manufacturer should formulate a low carbon system.

#### *1.2. Research Questions, Motivation and Contribution of the Model*

In the literature review section (Section 1.3) it can be clearly seen that some works seek to determine sustainable electrical supply chain inventory models and no research tries to determine a sustainable electrical supply chain inventory model with setup cost reduction and CO2 emissions. Our research questions this model: (1) what is the electrical power distribution factor's effect on the distribution substation, and the electrical power transmission factor's impact on the transmission substation and electrical power generation factor; (2) how much is the ordering quantity?; (3) what is customer's average electricity consumption time?; and (3) what is the retailer's selling price? To answer these questions in this study, we developed and solved a sustainable electricity supply chain inventory model with setup cost reduction and CO2 emissions while considering environmental parameters. The aim of this study is to examine the effect of price-dependent demand on the sustainable electrical supply chain inventory system under setup cost reduction and carbon emissions. The energy is transmitted via distribution networks to customers whose demand is determined by the price of electricity in order to find an optimal solution. The contributions of this paper are presented in Table 1.

#### *1.3. Literature Review*

Several research papers highlighted various integrated inventory models with various key parameters. Yang [1] deduced an inventory model considering lead time and crashing cost. Hoque [2] studied an integrated inventory model considering a normal distribution lead time, setup time, batch time and the cost of transportation. Sarkar et al. [3] deduced an inventory model by adding an imperfect production concept. A non-defective product adopts a binomial distribution function and demand adopts a mixture of the normal distribution function in their model. Mishra [4] formulated a production-inventory model with price dependent demand where the production depends on the rate of demand. Multiple buyers and a single-vendor model apply in a food inventory system studied by Fauza et al. [5]. Denizel et al. [6] formulated a lot size inventory model with setup costs decreased by various amounts depending upon the raw-materials and discussed the shortest path problem. Diaby [7] demonstrated a complete model to reduce setup cost and time with determining cut setup time and

minimizing the total cost. Nyea et al. [8] studied several inventory models for optimal investment of setup cost reduction or optimal setup times; and also discussed the queuing model to estimate work in the process level in their models. Later, Freimer et al. [9] demonstrated improvement in quality and setup cost reduction of the process in their model. Huang et al. [10] assumed setup cost policy reduction by an added investment cost. Sarkar and Moon [11] deduced a model by putting the policy of reorder point, quality improvement and lead time by considering that backorder rate has a great impact under a production process and is imperfect. Sarkar et al. [12] studied the concept of setup cost-reduction policy under quality improvement. Sarkar et al. [13] studied an effect of setup cost-reduction process in a two-echelon supply chain inventory model under deterioration and using the technique of quality improvement. Sarkar et al. [14] developed an integrated inventory model for a setup cost-reduction policy, carbon-emission policy and used the technique of the Stackelberg game approach to find the total cost. Little research has been done on a supply chain electrical energy inventory model. Banbury [15] was the first researcher to developed an electricity supply chain model. Thereafter, (Schneider et al. [16], Schneider et al. [17]) first developed this using a very simple inventory policy to find the electrical supply chain policy. Later, Taylor et al. [18] studied capacity and price competition in an electricity market by using a two-stage game theory model. Wu et al. [19] studied a new model to examine the formal generation with sporadic supply. Ouedraogo [20] formulated an electricity supply chain with demand in an African power system. Wangsa and Wee [21] studied an electrical supply chain inventory policy assuming the blackout cost. Recently, interesting research by Wangsa et al. [22] assumed a sustainable supply chain inventory model and the effect of price-dependent demand.


**Table 1.** Other authors contribution to this theme.

Finally, in this study we investigate many research articles involving supply chain inventory model-related carbon emissions, setup cost reduction and electricity energy. Research paper related to the above are the following: Hammami et al. [23] developed a multi-echelon supply chain model with reducing carbon emission, several manufacturing facilities, different outside suppliers, and distinct

distribution centers. Tang et al. [24] studied a carbon-emission policy with minimal frequency of shipments for a periodic inventory review system. Thereafter, Tang et al. [25] developed a sustainable supply chain network for consumers, and environmental manners are added by inventory, routing, and location.

#### *1.4. Methodology*

In this section, we consider a sustainable electricity supply chain power system with price-dependent demand electricity demand. In reality, it is shown that determining the capacity of a sustainable electrical supply chain system has the same methodology as finding the order quantity *q* in the supply chain inventory system. The supply chain inventory system involves a vendor–buyer coordination and freight forwarding. The buyer sells items to the customers and orders items from the vendor. The vendor produces the items and sends in batch to the buyer. The buyer will then sell the items to the customers. But in sustainable electrical supply chain system case, the electricity generated from a power generation will be transmitted through a transmission line and distribution substation to the customers whose demands are influenced by the price of electricity where the electricity is continuously supplied to consumers without any interruptions. The electricity demand is considered as *D*(*p*) kWh/year. The power generation produces the electricity in a batch size of *qTζηρ* kWh where *ζ* (positive integer) is the distribution factor's effect on the distribution substation, *η* (positive integer) is the transmission factor's impact on the transmission substation and *ρ* (positive integer) is a power generation factor. The finite power supply rate is *P* = *λD*(*p*) kWh/year, [*λ* > 1] and a setup cost. The electricity energy of *qTζη* kWh is supplied by the power generator to the transmission substation, then *qTζ* kWh of electricity is supplied to the distribution substation and *E* = *qT* kWh of electricity is consumed by the customers. Hence, to maximize the profit of a sustainable electrical supply chain system, we consider the sales revenue, production cost, setup cost reduction of the power generation, ordering cost of customers and transmission/distribution costs of substations. The transmission and distribution costs are functions of the power plant, the transmission substation and the distribution substation with maximum capacities of in *z<sup>x</sup> <sup>p</sup>* kWA, *z<sup>x</sup> <sup>t</sup>* kWA, *z<sup>x</sup> <sup>d</sup>* kWA, respectively. The comparison in supply chain inventory system and sustainable electrical supply chain system are shown in Figure 1.

**Figure 1.** Supply chain inventory and sustainable electrical supply chain.

The rest of the paper is prepared as follows: in Section 2, the assumption and notation for model formulation of the electrical supply chain inventory system are given; in Section 3, we describe the model formulation of the electrical supply chain inventory system. In Section 4, a numerical example is presented for validation of this model. A discussion of the managerial implications is presented in Section 5. Section 6 concludes the work by offering directions for future research work.

#### **2. Assumptions and Notations**

In this section, some notation can be used for development of the mathematical model (see Section 3) using the following assumptions:

#### *2.1. Assumptions*


#### *2.2. Notation*

Notations are used in the model are shown in Table 2 as follows:


#### **Table 2.** Notations.

#### **3. Model Formulation**

In this section, we explain how to find the total cost function with regard to customers, the distribution and transmission substations, carbon emission cost, and the total profit function for power generation. The total cost functions are given by the following components:

#### *3.1. Ordering Cost*

Ordering cost is defined as the customer's total cost per unit time:

$$TC\_O = \frac{AD(p)}{qT} \tag{1}$$

#### *3.2. Distribution Cost*

The same methodology as (Wangsa et al. [22]) to determine the distribution cost is obtained from the distribution substation. The cost of distribution for partial load *<sup>G</sup>* can be written as *<sup>G</sup>* <sup>=</sup> *<sup>c</sup>*1*z<sup>x</sup> d z y d* . The transmission function, increase in rate per kVA/mile when *z y <sup>d</sup>* increases. The distribution cost per kVA/mile is *GZ* = *βG* + (1 − *β*)*c*1, where 0 ≤ *β* ≤ 1 is the coefficient of the adjusted inverse function. Therefore *GZ* = *βc*<sup>1</sup> - *zx <sup>d</sup>*−*z y d z y d* + *c*1. The estimated total cost for the distribution substation as a function of demand, corrections factor and distance yields: *GD* = - *βc*<sup>1</sup> - *zx <sup>d</sup>*−*z y d z y d* + *c*<sup>1</sup> *D*(*p*)*m*3Ω. The actual power supply capacity is given by *z y <sup>d</sup>* = *qTζ* Ω; therefore, the cost for distribution can be written as in the following equation:

$$TC\_D = \frac{D(p)\beta c\_1 z\_d^x m\_3}{qT\zeta} + D(p)m\_3\Omega(1-\beta)c\_1 \tag{2}$$

#### *3.3. Carbon Emission Distribution Cost*

The carbon emission distribution cost has one component. This carbon emission component is related to total distribution rates. The component is summated and briefly represented by:

$$CE\_D = \stackrel{\circ}{c}\_1^{\circ} \left[ \frac{D(p) \mathfrak{F} \mathfrak{F} z\_d^{\chi} m\_3}{q T \zeta} + D(p) \mathfrak{F} m\_3 \Omega(1 - \beta) \right] \tag{3}$$

#### *3.4. Transmission Substation Cost*

This total cost comprises of the energy holding cost and transmission cost. The cost of energy holding at the transmission substation is presented by *<sup>r</sup>*<sup>1</sup> *<sup>p</sup>ζη<sup>E</sup>* <sup>2</sup> <sup>=</sup> *<sup>r</sup>*<sup>1</sup> *<sup>p</sup>ζηqT* <sup>2</sup> . The transmission cost can be calculated as: *Gt* <sup>=</sup> *<sup>D</sup>*(*p*)*βc*1*z<sup>x</sup> <sup>t</sup> m*<sup>2</sup> *qTζη* + *D*(*p*)*m*2Ω(1 − *β*)*c*1. Therefore, the total cost obtained at the transmission substation is the sum of energy holding cost and transmission cost. Therefore, the total transmission substations cost;

$$TC\_t = \frac{r\_1 p\_\zeta^\gamma \eta qT}{2} + \frac{D(p)\beta c\_1 z\_t^\chi m\_2}{qT\tilde{\chi}\eta} + D(p)m\_2\Omega(1-\beta)c\_1\tag{4}$$

#### *3.5. Carbon Emission Transmission Substation Cost*

The carbon emission transmission substation cost has three components. The first component is related to total energy holding rates at the transmission substation, the second and third components are related to the total transmission rates. These components are summated and briefly represented by

$$CE\_t = \frac{\stackrel{\wedge}{r\_1}\!\!\!\!p\!\!\!\!\!p\!\!\!\!\!\/\eta \!\!T}{2} + \frac{D(p)\beta\!\!\!\!\!\!\!\!\!\!\/\eta \!\!\!\!\!\!\/)z\_1^x m\_2 + D(p)m\_2\Omega(1-\beta)\mathbb{J}\!\!\!\!\!\/\eta \,\tag{5}$$

#### *3.6. Profit of Power Generation*

The total profit of power generation can be presented by:

$$
\Pi\_{PG} = \mathcal{S}R - PC - \mathcal{S}C - INVC - HC - TRC - CE\_{TR} - CE\_{PG} \tag{6}
$$

Sales revenue

$$SR = D(p)p\tag{7}$$

Production cost

$$PC = D(p)k \tag{8}$$

Power generation creates electrical energy in (*qTζηρ*) kWh in one production run. Therefore, the setup cost for power generation can be found by using the equation below. Setup cost is:

$$SC = \frac{D(p)v\_0e^{-\tau I}}{qT\zeta\eta\rho} \tag{9}$$

Total setup investment cost is:

$$INVC = \frac{D(p)I}{qT\zeta\eta\rho} \tag{10}$$

Energy holding cost is:

$$\mathcal{H}\mathbb{C} = r2k \frac{\left[E\_r^\prime \eta \rho \left(\frac{E\_r^\prime \eta}{P} + (\rho - 1) - \frac{E}{D(p)}\right) - \frac{\rho^2 \left[E\_r^\prime \eta\right]^2}{2P}\right] - \left[\frac{E\_r^\prime \eta}{D(p)} \left(1 + 2 + 3 \cdot \dots \dots + (\rho - 1)\right)E\right]}{\frac{E\_r^\prime \eta \rho}{D(p)}}$$

*HC* can be written as rewritten as *HC* <sup>=</sup> *<sup>r</sup>*2*<sup>k</sup> <sup>E</sup>ζη* 2 *ρ* <sup>1</sup> <sup>−</sup> *<sup>D</sup>*(*p*) *P* <sup>−</sup> <sup>1</sup> <sup>+</sup> <sup>2</sup>*D*(*p*) *P* . Putting *P* = *λD*(*p*) in *HC*. Therefore, the holding cost for is

$$\text{HCC} = r\_2 k \frac{qT\zeta\eta}{2\lambda} [\rho(\lambda - 1) - \lambda + 2] \tag{11}$$

Transmission cost can be calculated as:

$$G\_p = \frac{D(p)\beta c\_1 z\_p^\times m\_1}{qT\zeta\eta} + D(p)m\_1\Omega(1-\beta)c\_1\tag{12}$$

#### *3.7. Total Carbon Emission of Transmission Cost*

The carbon emission transmission cost has one component. This transmission emission component is related to total transmission rates. The component is summated and briefly represented by:

$$CE\_{TR} = \stackrel{\wedge}{c\_1} \left[ \frac{D(p)\xi\beta z\_p^{\times}m\_1}{qT\zeta\eta} + D(p)\xi m\_1\Omega(1-\beta) \right] \tag{13}$$

#### *3.8. Total Carbon Emission of Energy Holding Cost for Power Generation*

The carbon emission energy holding cost has one component. The holding emission component is related to total energy hold in a power generation plant. The component is summated and briefly represented by:

$$CE\_{PG} = r\_2^{\wedge} \left[ k \frac{qT \xi \zeta \eta}{2\lambda} [\rho(\lambda - 1) - \lambda + 2] \right] \tag{14}$$

Therefore, the total profit for power generation Π*PG* can be defined as:

$$\begin{split} \Pi\_{PG} = D(p)(p-k) - \frac{D(p)v\eta e^{-\tau l}}{qT\zeta\eta\rho} - \frac{D(p)l}{qT\zeta\eta\rho} - (r\_2 + \stackrel{\wedge}{r}\_2^{\gamma}\xi)k\frac{qT\zeta\eta}{2\lambda}[\rho(\lambda-1) - \lambda + 2] \\ - \left[\frac{D(p)\ell(c\_1+c\_1^{\gamma}\xi)z\_p^rx\_1}{qT\zeta\eta} + D(p)m\_1\Omega(1-\beta)(c\_1+\stackrel{\wedge}{c}\_1^{\gamma}\xi)\right] \end{split} \tag{15}$$

Therefore, the total profit is:

$$\begin{split} \Pi &= \Pi\_{PG} - T\mathbb{C}\_{O} - T\mathbb{C}\_{D} - C\mathbb{E}\_{t} - T\mathbb{C}\_{t} - C\mathbb{E}\_{t} \\ &= D(p) \left[ p - k - (m\_{1} + m\_{2} + m\_{3})\Omega(1 - \beta)(c\_{1} + \stackrel{\wedge}{c}\_{1}\mathbb{I}\_{\xi}) \right] \\ &- \frac{D(p)}{qT\zeta\eta\rho} \left[ A\xi\eta\rho + \beta(c\_{1} + \stackrel{\wedge}{c}\_{1}\mathbb{I}\_{\xi})z\_{d}^{x}m\_{3}\eta\rho + \beta(c\_{1} + \stackrel{\wedge}{c}\_{1}\mathbb{I}\_{\xi})z\_{m}^{x}m\_{2}\rho + v\_{0}e^{-\tau I} + I + \beta(c\_{1} + \stackrel{\wedge}{c}\_{1}\mathbb{I}\_{\xi}^{\times})z\_{p}^{x}m\_{1} \right] \\ &\qquad - \frac{qT\zeta\eta}{2\lambda} \left[ \Lambda(r\_{1} + \stackrel{\wedge}{r}\_{1}\mathbb{I}\_{\xi})p + (r\_{2} + \stackrel{\wedge}{r}\_{2}\mathbb{I}\_{\xi})k(\rho(\lambda - 1) - \lambda + 2) \right] \end{split} \tag{16}$$

Next, we analyze the consequence of *I*, *ζ*, *η*, *ρ*, *T* and *p* on ∏ fixed *q* by using the second order partial derivatives of Equation (16) with respect to *I*, *ζ*, *η*, *ρ*, *T* and *p*:

$$\begin{split} \frac{d\Pi}{d\xi} &= \frac{D(p)}{qT\xi^{2}\eta\rho} \left[ \beta(c\_{1} + \stackrel{\wedge}{c}\_{1}\xi)z\_{d}^{\scriptscriptstyle{x}}m\_{3}\eta\rho + \beta(c\_{1} + \stackrel{\wedge}{c}\_{1}\xi)z\_{d}^{\scriptscriptstyle{x}}m\_{2}\rho + v\_{0}e^{-\tau I} + I + \beta(c\_{1} + \stackrel{\wedge}{c}\_{1}\xi)z\_{p}^{\scriptscriptstyle{x}}m\_{1} \right] \\ & \qquad - \frac{qT\eta}{2\lambda} \left[ \lambda(r\_{1} + \stackrel{\wedge}{r}\_{1}\xi)p + (r\_{2} + \stackrel{\wedge}{r}\_{2}\xi)k(\rho(\lambda - 1) - \lambda + 2) \right] \end{split} \tag{17}$$

$$\frac{\partial^2 \Pi}{\partial \xi^2} = -\frac{2D(p)}{qT\xi^3 \eta \rho} \left[ \beta(\mathbf{c}\_1 + \mathbf{\hat{c}}\_1 \underline{\mathbf{y}}) z\_d^x m\_3 \eta \rho + \beta(\mathbf{c}\_1 + \mathbf{\hat{c}}\_1 \underline{\mathbf{y}}) z\_d^x m\_2 \rho + \upsilon\_0 e^{-\tau I} + I + \beta(\mathbf{c}\_1 + \mathbf{\hat{c}}\_1 \underline{\mathbf{y}}) z\_p^x m\_1 \right] < 0 \tag{18}$$

$$\begin{array}{l} \frac{\partial \Pi}{\partial \eta} = \frac{D(p)}{qT\xi\eta^{2}\rho} \left[ \beta(c\_{1} + \stackrel{\wedge}{c}\_{1}^{\wedge}\underline{\xi})z\_{1}^{x}m\_{2}\rho + v\_{0}e^{-\pi I} + I + \beta(c\_{1} + \stackrel{\wedge}{c}\_{1}^{\wedge}\underline{\xi})z\_{p}^{x}m\_{1} \right] \\\ -\frac{qT\zeta}{2\lambda} \left[ \lambda(r\_{1} + \stackrel{\wedge}{r}\_{1}^{\wedge}\underline{\xi})p + (r\_{2} + \stackrel{\wedge}{r}\_{2}^{\wedge}\underline{\xi})k(\rho(\lambda - 1) - \lambda + 2) \right] \end{array} \tag{19}$$

$$\frac{\partial^2 \Pi}{\partial \eta^2} = -\frac{2D(p)}{q T \zeta \eta^3 \rho} \left[ \beta (\mathbf{c}\_1 + \stackrel{\wedge}{c}\_1 \stackrel{\wedge}{\xi}) z\_t^x m\_2 \rho + \nu\_0 e^{-\tau I} + I + \beta (\mathbf{c}\_1 + \stackrel{\wedge}{c}\_1 \stackrel{\wedge}{\xi}) z\_p^x m\_1 \right] < 0 \tag{20}$$

$$\frac{\partial \Pi}{\partial \rho} = \frac{D(p)}{q T \zeta \eta \rho^2} \left[ v\_0 e^{-\tau I} + I + \beta (c\_1 + \stackrel{\wedge}{c\_1} \zeta) z\_p^x m\_1 \right] - \frac{q T \zeta \eta (r\_2 + \stackrel{\wedge}{r\_2} \xi) k (\lambda - 1)}{2 \lambda} \tag{21}$$

$$\frac{\partial^2 \Pi}{\partial \rho^2} = -\frac{2D(p)}{qT\tilde{\zeta}\eta\rho^3} \left[ v\_0 e^{-\tau I} + I + \beta(c\_1 + \hat{c}\_1^\times \underline{\chi}) z\_p^\times m\_1 \right] < 0 \tag{22}$$

$$\begin{aligned} \frac{d\Pi}{d\Gamma} &= \frac{D(p)}{q^{12}\zeta\eta\rho} \left[ A\zeta\eta\rho + \beta(c\_1 + \stackrel{\wedge}{c}\_1\xi)z\_d^x m\_3\eta\rho + \beta(c\_1 + \stackrel{\wedge}{c}\_1\xi)z\_f^x m\_2\rho + v\_0 e^{-\pi I} + I + \beta(c\_1 + \stackrel{\wedge}{c}\_1\xi)z\_p^x m\_1 \right] \\ &- \frac{q\zeta\eta}{2\lambda} \left[ \lambda(r\_1 + \stackrel{\wedge}{r}\_1\xi)p + (r\_2 + r\_2\xi)k(\rho(\lambda - 1) - \lambda + 2) \right] \end{aligned} \tag{23}$$

$$\frac{\partial^2 \Pi}{\partial T^2} = -\frac{2D(p)}{qT^3 \zeta \eta \rho} \left[ \begin{array}{c} A\zeta \eta \rho + \beta(c\_1 + \stackrel{\wedge}{c}\_1 \xi) z\_d^x m\_3 \eta \rho + \beta(c\_1 + \stackrel{\wedge}{c}\_1 \xi) \\ \times z\_t^x m\_2 \rho + v\_0 e^{-\tau I} + I + \beta(c\_1 + \stackrel{\wedge}{c}\_1 \xi) z\_p^x m\_1 \end{array} \right] < 0 \tag{24}$$

$$\begin{aligned} \overset{\text{dII}}{\frac{\partial \Pi}{\partial p}} &= a - bp - b \left[ p - k + (\beta - 1) \Omega (c\_1 + \stackrel{\wedge}{c\_1} \xi) (m\_1 + m\_2 + m\_3) \right] - \frac{q \Gamma \underline{\eta} (r\_1 + \stackrel{\wedge}{c\_1} \xi)}{2} \\ &+ \frac{b}{q \Gamma \underline{\eta} \rho \rho} \left[ A \zeta \eta \rho + \beta (c\_1 + \stackrel{\wedge}{c\_1} \xi) z\_d^x m\_3 \eta \rho + \beta (c\_1 + \stackrel{\wedge}{c\_1} \xi) z\_d^x m\_2 \rho + \nu \alpha e^{-\pi I} + I + \beta (c\_1 + \stackrel{\wedge}{c\_1} \xi) z\_p^x m\_1 \right] \end{aligned} \tag{25}$$

and

$$\frac{\partial \Pi}{\partial p} = -2b < 0\tag{26}$$

From Equations (18), (20), (22), (24) and (26) it can be concluded that, for any feasible solution of *q* the total profit function ∏ (Equation (16) is a concave function of *ζ*, *η*, *ρ*, *T* and *p*.

**Theorem 1.** *For the any positive integer (ζ*, *η*, *ρ*, *T*, *p), the required objective function* ∏ *(Equation (16) is concave function of I and q. Therefore, the maximum value* ∏ *(Equation (16) is settled at the point I and q which satisfies <sup>∂</sup>* <sup>∏</sup> *<sup>∂</sup><sup>I</sup>* <sup>=</sup> <sup>0</sup> *and <sup>∂</sup>* <sup>∏</sup> *<sup>∂</sup><sup>q</sup>* = 0*.*

**Proof.** To calculate the optimal solution for fixed-integers (*ζ*, *η*, *ρ*, *T*, *p*), used the partial derivatives with respect to *I* and *q*, as shown in the following equations:

$$\frac{\partial \Pi}{\partial I} = -\frac{D(p)\left(1 - \tau v\_0 \varepsilon^{-\tau I}\right)}{qT\zeta\eta\rho} \tag{27}$$

$$\begin{aligned} \frac{d\Pi}{d\eta} &= \frac{D(Ap)}{q^T T\_\xi \eta \rho} \left[ A\xi\eta\rho + \beta(c\_1 + \stackrel{\wedge}{c\_1}\xi\rangle z\_d^x m\_3 \eta\rho + \beta(c\_1 + c\_2) z\_t^x m\_2 \rho + v\_0 e^{-\tau I} + I + \beta(c\_1 + \stackrel{\wedge}{c\_1}\xi\rangle z\_p^x m\_1) \right] \\ &- \frac{T\_\xi \eta}{2\lambda} \left[ \lambda(r\_1 + \stackrel{\wedge}{r}\_1\xi)p + (r\_2 + \stackrel{\wedge}{r}\_2\xi)k(\rho(\lambda - 1) - \lambda + 2) \right] \end{aligned} \tag{28}$$

Now, set Equations (27) and (28) equal to zero and solve for *I* and *q*:

$$-\frac{D\left(p\right)\left(1-\tau v\_0 e^{-\tau I}\right)}{qT\zeta\eta\rho} = 0 \Rightarrow I = \frac{1}{\pi}\log[\tau v\_0] \tag{29}$$

Then:

*D*(*Ap*) *<sup>q</sup>*2*Tζηρ <sup>A</sup>ζηρ* + *<sup>β</sup>*(*c*<sup>1</sup> + <sup>∧</sup> *c*1*ξ*)*z<sup>x</sup> <sup>d</sup>m*3*ηρ* + *<sup>β</sup>*(*c*<sup>1</sup> + *<sup>c</sup>*2*ξ*)*z<sup>x</sup> <sup>t</sup> <sup>m</sup>*2*<sup>ρ</sup>* <sup>+</sup> *<sup>v</sup>*0*e*−*τ<sup>I</sup>* <sup>+</sup> *<sup>I</sup>* <sup>+</sup> *<sup>β</sup>*(*c*<sup>1</sup> <sup>+</sup> <sup>∧</sup> *c*1*ξ*)*z<sup>x</sup> pm*<sup>1</sup> −*Tζη* 2*λ <sup>λ</sup>*(*r*<sup>1</sup> + <sup>∧</sup> *<sup>r</sup>*1*ξ*)*<sup>p</sup>* + (*r*<sup>2</sup> + <sup>∧</sup> *r*2*ξ*)*k*(*ρ*(*λ* − 1) − *λ* + 2) = 0 <sup>⇒</sup> *<sup>q</sup>* <sup>=</sup> <sup>1</sup> *Tζη* 2*λ* ⎡ ⎢ ⎣*D*(*p*) ⎡ ⎢ ⎣ *<sup>A</sup>ζηρ* + *<sup>β</sup>*(*c*<sup>1</sup> + <sup>∧</sup> *c*1*ξ*)*z<sup>x</sup> <sup>d</sup>m*3*ηρ* + *β*(*c*<sup>1</sup> + *c*2) ×*z<sup>x</sup> <sup>t</sup> <sup>m</sup>*2*<sup>ρ</sup>* <sup>+</sup> *<sup>v</sup>*0*e*−*τ<sup>I</sup>* <sup>+</sup> *<sup>I</sup>* <sup>+</sup> *<sup>β</sup>*(*c*<sup>1</sup> <sup>+</sup> <sup>∧</sup> *c*1*ξ*)*z<sup>x</sup> pm*<sup>1</sup> ⎤ ⎥ ⎦ ⎤ ⎥ ⎦ *ρ <sup>λ</sup>*(*r*1+ <sup>∧</sup> *<sup>r</sup>*1*ξ*)*p*+(*r*2+ <sup>∧</sup> *r*2*ξ*)*k*(*ρ*(*λ*−1)−*λ*+2) (30)

In order to prove the concavity of the required objective profit function ∏; can show that the following conditions:

$$\frac{\partial^2 \Pi}{\partial I^2} = -\frac{D(p)\left(1 + \tau^2 v\_0 e^{-\tau I}\right)}{q T \zeta \eta \rho} < 0\tag{31}$$

$$\frac{\partial^2 \Pi}{\partial \eta^2} = -\frac{2D(p)}{q^3 T \zeta \eta \rho} \left[ \begin{array}{c} A\zeta \eta \rho + \beta(c\_1 + \stackrel{\wedge}{c}\_1 \stackrel{\chi}{\varphi}) z\_d^x m\_3 \eta \rho + \beta(c\_1 + c\_2) z\_t^x \\ \times m\_2 \rho + \upsilon\_0 e^{-\tau I} + I + \beta(c\_1 + \stackrel{\wedge}{c}\_1 \stackrel{\chi}{\varphi}) z\_p^x m\_1 \end{array} \right] < 0 \tag{32}$$

$$\frac{\partial \Pi}{\partial I \partial q} = \frac{\partial \Pi}{\partial q \partial I} = \frac{D(p)\left(1 - e^{-\tau I} \tau v\_0\right)}{q^2 T \zeta \eta \rho} \tag{33}$$

The Hessian matrix can be found as follows:

$$H = \begin{bmatrix} \frac{\partial^2 \Pi}{\partial \xi^2} & \frac{\partial^2 \Pi}{\partial \xi \partial \eta} \\ \frac{\partial^2 \Pi}{\partial \eta \partial I} & \frac{\partial^2 \Pi}{\partial \eta^2} \end{bmatrix} = \begin{bmatrix} -\frac{D(p)\left(1 + \tau^2 \tau\_0 \varepsilon^{-1}\right)}{q^T \xi \eta \rho} & \frac{D(p)\left(1 - \varepsilon^{-1} \tau\_0 \eta\right)}{q^2 T \xi \eta \rho} \\\\ \frac{D(p)\left(1 - \varepsilon^{-1} \tau\_0 \eta\right)}{q^T T \xi \eta \rho} & -\frac{2D(p)}{q^T T \xi \eta \rho} \\\\ +\beta \left(c\_1 + c\_2\right) \varpi\_1^x m\_2 \rho + v\_0 e^{-\pi} I \\ + I + \beta (c\_1 + c\_2^\dagger \xi) \varpi\_p^x m\_1 \end{bmatrix}.$$

The determinant of |*H*| is followed:

<sup>|</sup>*H*<sup>|</sup> <sup>=</sup> <sup>−</sup> [*D*(*p*)] 2 *<sup>q</sup>*4*T*2*ζ*2*η*2*ρ*<sup>2</sup> <sup>+</sup> <sup>2</sup>[*D*(*p*)] 2 *e*−*τIτv*<sup>0</sup> *<sup>q</sup>*4*T*2*ζ*2*η*2*ρ*<sup>2</sup> <sup>+</sup> <sup>2</sup>[*D*(*p*)] 2 *e*−*τ<sup>I</sup> Iτ*2*v*<sup>0</sup> *<sup>q</sup>*4*T*2*ζ*2*η*2*ρ*<sup>2</sup> <sup>+</sup> <sup>2</sup>*A*[*D*(*p*)] 2 *e*−*τIτ*2*v*<sup>0</sup> *q*4*T*2*ζηρ* +[*D*(*p*)] 2 *e*−2*τIτ*2*v*<sup>2</sup> 0 *<sup>q</sup>*4*T*2*ζ*2*η*2*ρ*<sup>2</sup> <sup>+</sup> <sup>2</sup>[*D*(*p*)] 2 *e*−*τ<sup>I</sup> βτ*2*c*1*m*3*v*0*z<sup>x</sup> d <sup>q</sup>*4*T*2*ζ*2*ηρ* <sup>+</sup> <sup>2</sup>[*D*(*p*)] 2 *e*−*τ<sup>I</sup> βτ*<sup>2</sup> <sup>∧</sup> *c*1*ξm*3*v*0*z<sup>x</sup> d q*4*T*2*ζ*2*ηρ* +2[*D*(*p*)] 2 *e*−*τ<sup>I</sup> βτ*2*c*1*m*2*v*0*z<sup>x</sup> t <sup>q</sup>*4*T*2*ζ*2*η*2*<sup>ρ</sup>* <sup>+</sup> <sup>2</sup>[*D*(*p*)] 2 *e*−*τ<sup>I</sup> βτ*<sup>2</sup> <sup>∧</sup> *c*1*ξm*2*v*0*z<sup>x</sup> t <sup>q</sup>*4*T*2*ζ*2*η*2*<sup>ρ</sup>* <sup>+</sup> <sup>2</sup>[*D*(*p*)] 2 *e*−*τ<sup>I</sup> βτ*2*c*1*ξm*1*v*0*z<sup>x</sup> p q*4*T*2*ζ*2*η*2*ρ*<sup>2</sup> <sup>+</sup>2[*D*(*p*)] 2 *e*−*τ<sup>I</sup> βτ*<sup>2</sup> <sup>∧</sup> *c*1*ξm*1*v*0*z<sup>x</sup> p q*4*T*2*ζ*2*η*2*ρ*<sup>2</sup> = [*D*(*p*)] 2 ⎡ ⎢ ⎢ ⎣*τ*2*v*<sup>2</sup> <sup>0</sup>*e*−2*τI*+2*e*−*τIτv*<sup>0</sup> ⎡ ⎢ ⎢ ⎣ 1 + *τI* + *Aτζηρ* + *τβ* <sup>∧</sup> *c*1*ξ ηρm*3*z<sup>x</sup> <sup>d</sup>* + *<sup>ρ</sup>m*2*z<sup>x</sup> <sup>t</sup>* + *m*1*z<sup>x</sup> p* +*τβ* <sup>∧</sup> *c*1*ξ ηρm*3*z<sup>x</sup> <sup>d</sup>* + *<sup>ρ</sup>m*2*z<sup>x</sup> <sup>t</sup>* + *m*1*z<sup>x</sup> p* ⎤ ⎥ ⎥ ⎦−1 ⎤ ⎥ ⎥ ⎦ *q*4*T*2*ζ*2*η*2*ρ*<sup>2</sup> <sup>=</sup> [*D*(*p*)] 2 [*τ*2*v*<sup>2</sup> <sup>0</sup>*e*−2*τI*+2*B*−1] *<sup>q</sup>*4*T*2*ζ*2*η*2*ρ*<sup>2</sup> > <sup>0</sup>

$$\text{where } B = \varepsilon^{-\tau I} \tau v\_0 \left[ \begin{array}{c} 1 + \tau I + A\tau \zeta \eta \rho + \tau \beta \mathbf{c}\_1^\top \xi \left( \eta \rho m\_3 z\_d^\mathbf{x} + \rho m\_2 z\_t^\mathbf{x} + m\_1 z\_p^\mathbf{x} \right) \\ + \tau \beta \mathbf{c}\_1^\top \xi \left( \eta \rho m\_3 z\_d^\mathbf{x} + \rho m\_2 z\_t^\mathbf{x} + m\_1 z\_p^\mathbf{x} \right) \end{array} \right] > 1.$$

The behavior of the concavity for objective function ∏ with respect to our decision variables *ζ*, *η*, *ρ*, *T*, *p*, *I* and *q*, the algorithm is similar to Wangsa et al. [22], and was developed to draw the global maximum feasible solutions for *I*∗, *ζ*∗, *η*∗, *ρ*∗, *T*∗, *p*∗, *q*∗ and ∏∗.

*3.9. Algorithm*

Step 1. First compute *I* from Equation (29). Step 2. a. set *ζ* = 1. b. set *η* = 1. c. set *ρ* = 1. d. set *T* = 1. e. set *p* = 20. Step 3. Compute the optimal *q*∗ by using Equation (30). Step 4. After completed Step 3, compute the actual electrical power capacities. a. Power generation The actual capacity of the power generation *z y <sup>p</sup>* = *qTζηρ* Ω. If *z y <sup>p</sup>* ≤ *<sup>z</sup><sup>x</sup> <sup>p</sup>* then find *<sup>q</sup>* <sup>=</sup> *<sup>z</sup> y p <sup>T</sup>ζηρ* <sup>Ω</sup> and go to Step 5. b. Transmission substation The actual capacity of the transmission substation *z y <sup>t</sup>* = *qTζη* Ω. If *z y <sup>t</sup>* ≤ *<sup>z</sup><sup>x</sup> <sup>t</sup>* then find *<sup>q</sup>* <sup>=</sup> *<sup>z</sup> y t <sup>T</sup>ζη* <sup>Ω</sup> and go to Step 5. c. Distribution substation The actual capacity of the distribution substation *z y <sup>d</sup>* = *qTζ* Ω. If *z y <sup>t</sup>* ≤ *<sup>z</sup><sup>x</sup> <sup>t</sup>* then find *<sup>q</sup>* <sup>=</sup> *<sup>z</sup> y d <sup>T</sup><sup>ζ</sup>* <sup>Ω</sup> and go to Step 5. Step 5. Computed ∏ from Equation (16). Step 6. Set *ζ* = 1 + 1 and repeated Step 3 to Step 5. Step 7. If ∏(*q*∗ *<sup>ζ</sup>* , *ζ*, *ηζ* , *ρζ* , *T<sup>ζ</sup>* , *p<sup>ζ</sup>* ) ≥ ∏(*q*<sup>∗</sup> *<sup>ζ</sup>*−1, *<sup>ζ</sup>* <sup>−</sup> 1, *ηζ*−1, *ρζ*−1, *<sup>T</sup>ζ*−1, *<sup>p</sup>ζ*−1) then go to Step 8. Otherwise go to Step 6. Step 8. Set *η* = 1 + 1 and repeated Step 2b to Step 7. Step 9. If ∏(*q*∗ *<sup>η</sup>*, *ζη* <sup>∗</sup>, *η*, *ρη*, *Tη*, *pη*) ≥ ∏(*q*<sup>∗</sup> *<sup>η</sup>*−1, *ζη*−<sup>1</sup> <sup>∗</sup>, *η* − 1, *ρη*−1, *Tη*−1, *pη*−1) then go to Step 9. Otherwise go to Step 8. Step 10. Set *ρ* = 1 + 1 and repeated Step 2c to Step 9. Step 11. If ∏(*q*∗ *<sup>ρ</sup>*, *ζρ* ∗, *ηρ* <sup>∗</sup>, *ρ*, *Tρ*, *pρ*) ≥ ∏(*q*<sup>∗</sup> *<sup>ρ</sup>*−1, *ζρ*−<sup>1</sup> <sup>∗</sup>, *ηρ*−1, *ρ* − 1, *Tη*−1, *pη*−1) then go to Step 12. Otherwise go to Step 10.

Step 12. Set *T* = 1 + 1 and repeated Step 2d to Step 11. Step 13. If ∏(*q*∗ *<sup>T</sup>*, *ζT*∗, *ηT*∗, *ρ*<sup>∗</sup> *<sup>T</sup>*, *T*, *pT*) ≥ ∏(*q*<sup>∗</sup> *<sup>T</sup>*−1, *<sup>ζ</sup>T*−<sup>1</sup> <sup>∗</sup>, *ηT*−<sup>1</sup> ∗, *ρ*∗ *<sup>T</sup>*−1, *<sup>T</sup>* <sup>−</sup> 1, *pT*−1) then go to Step 14. Otherwise go to Step 12. Step 14. Set *p* = 1 + 1 and repeated Step 2e to Step 13. Step 15. If ∏(*q*∗ *<sup>p</sup>*, *ζ <sup>p</sup>* <sup>∗</sup>, *η<sup>p</sup>* ∗, *ρ*∗ *<sup>p</sup>*, *T*<sup>∗</sup> *<sup>p</sup>* , *p*) ≥ ∏(*q*<sup>∗</sup> *<sup>p</sup>*−1, *<sup>ζ</sup> <sup>p</sup>*−<sup>1</sup> <sup>∗</sup>, *ηp*−<sup>1</sup> ∗, *ρ*∗ *<sup>p</sup>*−1, *<sup>T</sup>*<sup>∗</sup> *<sup>p</sup>*−1, *<sup>p</sup>* <sup>−</sup> <sup>1</sup>) then go to Step 16. Otherwise go to Step 14. Step 15. If ∏(*q*∗ *<sup>p</sup>*, *ζ <sup>p</sup>* <sup>∗</sup>, *η<sup>p</sup>* ∗, *ρ*∗ *<sup>p</sup>*, *T*<sup>∗</sup> *<sup>p</sup>* , *p*∗) ≥ ∏(*q*<sup>∗</sup> *<sup>p</sup>*−1, *<sup>ζ</sup> <sup>p</sup>*−<sup>1</sup> <sup>∗</sup>, *ηp*−<sup>1</sup> ∗, *ρ*∗ *<sup>p</sup>*−1, *<sup>T</sup>*<sup>∗</sup> *<sup>p</sup>*−1, *<sup>p</sup>* <sup>−</sup> <sup>1</sup>) then find *<sup>I</sup>*∗, *ζ*∗, *η*∗, *ρ*∗, *T*∗, *p*∗, *q*∗ and go to Step 16. Step 16. Stop.

#### **4. Numerical Example**

In this section, used data to demonstrate the application of the model. This study considers the sustainable electrical energy supply chain inventory system in the Taiwanese electrical production industry to determine the optimal ordering quantity and total profit. The parameters in this section are assumed from a previous published paper. The data values are:

#### *4.1. Customer Data*

Let *a* = 80 units, *b* = 2 units, *A* = \$100/*unit*/*order*, *r*<sup>1</sup> = \$ 0.004/*unit*/*year*, ∧ *r*<sup>1</sup> = 0.002 *ton*/*year*, *r*<sup>2</sup> = \$ 0.002/*unit*/*year*, <sup>∧</sup> *r*<sup>2</sup> = 0.003 *ton*/*year*, *β* = 0.01 *unit*, Ω = 1.2 *kVA*/*kWh*, *v*<sup>0</sup> = \$7/*setup*, *ξ* = \$ 1/*ton*/*year* and *τ* = 0.2 *unit*.

#### *4.2. Transmission Substation, Distribution Substation and Power Generation Data*

Let *k* = \$4/*kWh*, *λ* = 2 *units*, *c*<sup>1</sup> = \$0.00011/*kVA*/*mile*, *z<sup>x</sup> <sup>p</sup>* = 6 *kVA*, *z<sup>x</sup> <sup>t</sup>* = 5 *kVA*, *z<sup>x</sup> <sup>d</sup>* = 2 *kVA*, *m*<sup>1</sup> = 0.2 *mile*, *m*<sup>2</sup> = 0.15 *mile* and *m*<sup>3</sup> = 0.1 *mile*.

Based on input data, by using Algorithm and Mathematica 9.0, the optimal solutions (see Tables 3 and 4); *ζ*∗ = 2, *η*∗ = 1, *ρ*∗ = 1, *T*∗ = 1, *p*∗ = 22, *D*(*p*) <sup>∗</sup> = 36, *I*<sup>∗</sup> = 1.68236, *q*<sup>∗</sup> = 183.875, *z y*∗ *<sup>p</sup>* = 441.30, *z y*∗ *<sup>t</sup>* = 441.30, *z y*∗ *<sup>d</sup>* = 441.30 and Π<sup>∗</sup> = 601.653. Furthermore, Figures 2 and 3 for graphical representations of the total profit function is a concave with respect to feasible optimal values *p*∗ and *q*∗.


**Table 3.** Details of the procedures for the solution.

\* The local maximum solution; ← the optimal solution.


**Table 4.** The results.

**Figure 2.** The contour of the objective function ∏∗ with respect to *ζ* = 2, *η* = 1, *ρ* = 1, *I* = 1.68236 and *T* = 1.

From Table 3, we obtain the *ρ* = 1 batch. It is intended that the electricity produced by a power generator is 8.55471 kWh. But not all the electricity energy (8.55471 kWh) induced is transmitted at once, but in periods with 8.55471 kWh each. As the generator has a device to minimize the total energy holding cost of the electricity, for each batch, the transmission substation obtains 4.277 kWh of electricity; it then transmits in 2 batches to the distribution station at 6.41593 kWh each. This is done to minimize the transmission cost and distribution cost. Based on those results, the total electrical power consumption of customer 183.875 kW. The demand of customer is 36 kWh/year.

**Figure 3.** Concavity of ∏∗ with respect to *p* and *q* at *ζ* = 2, *η* = 1, *ρ* = 1, *I* = 1.68236 and *T* = 1.

#### **5. Discussion and Managerial Implication**

In this section, we consider the effect of changes in the main parameters as well as summarize the results of the sensitivity analysis, which shows (see Table 5) the impact of each of the parameters which are *a*, *b*, *k*, *λ*, *v*0, *τ*, *A*, ∧ *c*1, ∧ *r*1, ∧ *r*<sup>2</sup> and *β*, respectively, on the total profit.


**Table 5.** Sensitivity analysis of key parameters of the model.

**Table 5.** *Cont*.


Based on the computational results (Table 5), the following managerial insights can be obtained:

#### *5.1. Impact on Demand Parameters*

An increase in the value of *a* results in an increase in demand, which forces the retailer to increase the selling price *p*, fixed *ζ*, *η*, *ρ*, *T*, *I*, increases the actual capacity of the distribution, transmission, power generation substation and the customer's power consumption in order to increase the electrical supply chain profit. On the other hand, an increase in the value of *b* could reduce the demand. So, the retailer then reduces the selling price *p*, fixed *η*, *ρ*, *T*, *I*, reduces the actual capacity of the distribution, transmission, power generation substation and the customer's power consumption. In this case, although the demand rate could be maintained on the higher side with the high value of *b*. Therefore, the profit of the electrical supply chain system would be decreased.

#### *5.2. Impact on Production Parameters*

An increase in the value of *λ* results in an increase in production with fix at demand rates, which forces the retailer to fix the selling price *p*, *ζ*, *η*, *ρ*, *T*, *I*, decreases the actual capacity of the distribution, transmission, power generation substation and the customer's power consumption in order to increase the electrical supply chain profit. An increase in the value of *k* could reduce the demand. So, the retailer then increases the selling price *p* with fixed *ζ*, *η*, *ρ*, *T*, *I*, and increases an actual capacity of

the distribution, transmission, power generation substation and the customer's power consumption. In this case, although the demand rate could be maintained in the less with the high value of *k*, customer's power consumption is higher, and actual capacity higher. Therefore, the profit of the electrical supply chain system is decreased.

#### *5.3. Impact on Setup Cost Reduction Parameters*

An increase in the value of *v*<sup>0</sup> results in an increase in setup cost with fixed at demand rates, which forces the retailer to fix the *p*, *ζ*, *η*, *ρ*, *T*, increases the actual capacity of the distribution, transmission, power generation substation and the customer's power consumption, increases the setup investment cost, and decreases the electrical supply chain profit. In this case, the profit of the electrical supply chain system decreased because all actual capacity and customer's power consumption with fixed selling price increase. On the other hand, an increase in the value of *τ* results in a decrease in setup cost with a fix at demand rates, which forces the retailer to fix the *p*, *ζ*, *η*, *ρ*, *T*, decreases the actual capacity of the distribution, transmission, power generation substation and the customer's power consumption, increases the setup investment cost and increases the electrical supply chain profit. In this case, the profit of the electrical supply chain system is increased because the demand rate maintained in the fixed with the high value of *τ*, customer's power consumption is lower, and all actual capacities are lower.

#### *5.4. Impact on Ordering Cost*

An increase in the value of *A* results in an increases the customer's power consumption rate with fixed demand rates; selling price increases the actual capacity of the distribution, transmission, and power generation substation. This result indicates that profit of the electrical supply chain system decreased because of higher of all actual capacity and customer's power consumption rates with fixed selling price.

#### *5.5. Impact on Loss Factor*

An increase in the value of all carbon emission parameters *β*, results in an no change of the customer's power consumption rate; demand rates; selling price, increases the actual capacity of the distribution, transmission, and power generation substation. This result indicates that, total profit can be unchanged because of all actual capacity and customer power consumption rates with fixed selling price are unchanged when small changes of *β*. Therefore, small changes of loss factor result in unchanged profit.

#### *5.6. Impact on Carbon Emission Parameters*

An increase in the value of all carbon emission parameters <sup>∧</sup> *<sup>c</sup>*1, <sup>∧</sup> *<sup>r</sup>*<sup>1</sup> and <sup>∧</sup> *r*<sup>2</sup> results in an increase the customer's power consumption rate with fixed demand rates; selling price increases the actual capacity of the distribution, transmission, and power generation substation. This result indicates that profit of the electrical supply chain system decreased because of higher actual capacity and customer power consumption rates at a fixed selling price.

#### **6. Conclusions and Future Research**

In this paper, a sustainable electricity supply chain mathematical model that assumes linear price-dependent customer demands where the price is a decision variable with reduction of setup cost under carbon emission, is considered. This model has been developed based on the inventory management theory, and examined how the all optimal decision variables and the total profits for sustainable electrical supply chain are affected by key parameters. Based on our computational results, it supplies managerial insights to the production system and marketing managers to help in planning a successful and sustainable electrical energy supply chain. For a future study, researchers can

extend the present model to include green technology investment under carbon emissions regulation. Researchers can also study incorporating price discount strategies as well as the effect of green technology investment under carbon emissions with a multi-transmission and distribution substation.

**Author Contributions:** Data Curation, Software U.M. and J.-Z.W.; Formal Analysis, Investigation, Writing Original Draft Preparation, Writing-Review, Visualization and Editing U.M.; Conceptualization, Methodology, Validation A.C and J.-Z.W.; Funding Acquisition, Supervision, Project Administration J.-Z.W.

**Funding:** This research was funded by Ministry of Science and Technology, Taiwan, grant number MOST106-2218-E-007-024 & MOST108-2911-I-031-502; and The APC was funded by MOST106-2218-E-007-024 & MOST108-2911-I-031-502.

**Acknowledgments:** The authors would like to thank the editor and anonymous reviewers for their valuable and constructive comments, which have led to a significant improvement in the manuscript.

**Conflicts of Interest:** The authors declare no conflict of interest.

#### **References**


© 2019 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).
