Optimization Theory, Methods, and Applications in Engineering 2014
View this Special IssueResearch Article  Open Access
Zhan’an Zhang, Xingguo Cai, "EnergySaving Generation Dispatch Using Minimum Cost Flow", Mathematical Problems in Engineering, vol. 2015, Article ID 562462, 9 pages, 2015. https://doi.org/10.1155/2015/562462
EnergySaving Generation Dispatch Using Minimum Cost Flow
Abstract
This study uses a minimum cost flow method to solve a dispatch problem in order to minimize the consumption of coal in the dispatching of a thermal power system. Lowcarbon generation dispatching is also considered here since the scheduling results are consistent with energysaving generation dispatch. Additionally, this study employs minimum coal consumption as an objective function in considering the output constraints, load balance constraints, line loss, ramp rate limits, spinning reverse needs, prohibited operating zone requirements, security constraints, and other common constraints. The minimum cost flow problem, considering the loss of network flow, is known as a generalized network flow problem, which can be expressed as a quadratic programming problem in mathematics. Accordingly, the objective function was solved by LINGO11, which was used to calculate a network in a single time; a continuous period dispatch plan was obtained by accumulating each period network flow together. This analysis proves it feasible to solve a minimal cost flow problem with LINGO11. Theoretical analysis and numerical results prove the correctness and effectiveness of the proposed method.
1. Introduction
There seems to be rather compelling evidence that global warming is an issue that we seriously need to be concerned about today [1, 2]. Carbon dioxide accounts for 80% of the greenhouse effect, and rising carbon dioxide levels are the main cause of global warming [3, 4]. China pledged to reduce its carbon intensity by 40–45% by 2020 based on 2005 levels. By the end of 2012, 71.5% of China’s generation capacity was from thermal power, of which 92.5% was the product of coalfired generation [5]. Coalfired electricity consumes about 50% of China’s coal production, and the CO_{2} emissions from power generation account for 40% of the total CO_{2} emissions in China. Accordingly, China’s power industry has implemented energysaving generation dispatching (ESGD) and lowcarbon generation dispatching (LCGD).
ESGD is one of the most important problems in power system operations requiring load demand at minimum total fuel cost while accounting for various unit and system constraints. The ESGD model is an optimization problem that considers linear and nonlinear characteristics, including power balance constraints, generation limit constraints, node voltage constraints, ramp rate limits, spinning reverse needs, prohibited operating zone requirements, and security constraints, among others [6]. In this paper, lowcarbon generation dispatching is also considered, since the scheduling results are consistent with energysaving generation dispatching. In practical scheduling applications, a daily scheduling period is generally divided into 24 or more intervals; therefore, dispatching of each period can be solved as a static optimization problem.
Many mathematical techniques have been developed and applied to dispatch problem such as linear programming [7], interiorpoint method [8], Lagrangian relaxation algorithm [9], quadratic programming [10] and other traditional algorithms. These algorithms essentially need some problem simplification such that the problem is linear or convex. Thus, a true global minimum cannot be guaranteed [11]. The dynamic programming method [12] has also been successfully used in solving the dispatch problems; however, this method may result in “curse of dimensionality.” More recently, the metaheuristic algorithms such as particle swarm optimization (PSO) [13], genetic algorithm (GA) [14], simulated annealing (SA) [15], evolutionary programming (EP) [16], and ant colony optimization (ACO) [17], have also been considered in the context of economic dispatch. However, the biggest problem that metaheuristic algorithms faced is that the computational efficiency is rather low. The related parameters are not easy to set up and the computational time is long. It is not easy to implement when applied to large electrical power systems. Their highly heuristic nature usually leads to suboptimal solutions. Differential evolution (DE) is a stochastic search based method [18], which can present a simple structure, fast convergence speed, and robustness. However, DE fast convergence might lead the direction of the search toward a local optimal and premature solution. Essentially, the economic dispatch problem is a large scale nonlinear programming problem. In pursuit of the optimal solution for economic dispatch, various hybrid methods have been investigated and implemented [5, 19–22]. Unfortunately, these hybrid algorithms normally take lengthy calculation time when compared with the mathematical optimization methods. Moreover, previous algorithms rarely considered network structures, we use a network flow method here to make up for the deficiency, because the network flow method can well retain the topology of the network.
Networks provide a useful way for modeling power system problems and are used extensively in power system dispatching [23, 24]. As an important network problem, ESGD problems can be formulated and solved as minimum cost flow problems when the cost is a quadratic function of the power, which is a nonlinear, minimum cost flow problem.
The main objective of this study is to introduce convex quadratic programming to solve the ESGD problem, since coal consumption and network losses are all convex functions of the power flow through a network. In order to do that this study employed LINGO11 to solve quadratic programming problems accounting for linear and nonlinear equality constraints and inequality constraints. Accordingly, a minimum cost flow algorithm was used to solve the ESGD model and calculate a network within a single moment. From that process, a continuous period of ESGD planning was obtained by accumulating period network flow results. This process confirmed that the minimal cost flow method was successful for solving the ESGD problem and, therefore, has value for these types of applications.
2. Problem Formulation
2.1. The Mathematical Model of ESGD Problem
The ESGD problem determines the optimal schedule of the available generating units to simultaneously minimize the generation cost rate and meet the load demand of a power system while meeting various equality and inequality constraints. This mathematical model can be described as follows:where is the cost function of the th generator; , , and are the cost coefficients of the th generator; is the power of the th generator at time; and represents the number of generators committed to the operating system.
2.2. The Mathematical Model of LowCarbon Generation Dispatching
Achieving the lowest carbon emission is the target of lowcarbon generation dispatching. This mathematical model can be described as follows:where is the electricalcarbon characteristic function. This formula represents the CO_{2} emissions when the output of unit is at time (t/h), which can be expressed as where is the CO_{2} emission coefficient of the fuel used in a power source; the standard coal emission factor is 2.77, which means that 2.77 kg CO_{2} can be discharged for every 1 kg of standard coal burnt; is the calorific value of unit fuel, which is 8.14 kWh/kg of standard coal; and is power generation efficiency, which can be expressed aswhere 3600 is the electric heating value (kJ/kWh) and 29308 is the calorific value of standard coal (kJ/kg).
If the coal consumption function and electricalcarbon characteristic satisfies the following relationshipwhere is a constant, then we can get then the 2 kinds of scheduling results are consistent [25].
2.3. Constraints
The objective function needs to satisfy the following constraints.
(1) Power Constraints. The power constraints include generator output, transformer capacity, and line transmission limits. Considerwhere is the lower and is the upper output limit of unit , respectively.
(2) Node Voltage Constraints. Considerwhere is the voltage of node at time and is the lower and is the upper limit of node , respectively.
(3) Power Balance Constraints. Considerwhere is the total load demand and is the transmission network losses, which is a function of unit power outputs that can be represented using the coefficients:where is the loss coefficient square matrix; is the loss coefficient vector; and is the loss coefficient constant [26].
(4) Operation Ramp Rate Limits. The power output of a practical generator cannot be adjusted instantaneously without limits. The operating range for all online units is restricted by their ramp rate limits during each dispatch period. Therefore, the dispatch output of a generator should be limited by the constraints of up and down ramp rates [27], which are given as follows:where and are the rampdown and rampup rate limits of the th thermal unit, respectively [28]. If the unit ramp rate limits are considered, the real power operating limits are modified as follows:
(5) Spinning Reserve. The added spinning reserve factor must be considered to prevent a sudden large load to the system or a failure in a certain large unit requirement. This condition can be explained as follows:where is the spinning reserve in the th hour and is the ON and OFF status of the th conventional unit at the period ( represents OFF status and represents ON status) [6].
(6) Prohibited Operating Zone. The prohibited operating zone is the range of prohibited output power resulting from the physical limitations of machine components, steam valves, vibration in the shaft bearing, and other conditions that can cause discontinuity in the electrical energy cost curve. Therefore, some units must be considered as prohibited zones in practical operation. The feasible operating zones of thermal units can be described as follows [29]: where is the number of prohibited zones in the th generator curve; is the index of prohibited zone of the th generator; and and are the upper and lower limits of the th prohibited zone of unit , respectively.
(7) Security Constraints. For secure operation, the transmission line loading () is restricted by its upper limit as follows:where is the total number of lines [30].
3. The Minimal Cost Flow Method
3.1. Network Flow Theory Introduction
Network flow problems can be described with graph theory, where a graph represents a network of nodes and connecting arcs. If a path exists between any two pairs of vertices in a graph, then that graph is a connected graph where each arc has a specific direction.
As depicted in Figure 1, the common abstraction that models a flow network is a directed graph , where is the set of vertices and is the set of edges over these vertices.
is the edge connecting nodes and , which has a predetermined direction, flow , the cost of the unit flow , and flow rate limit (as shown in Figure 1). Each edge has a flow that defines the number of units of the commodity that flows from to . An edge also has a capacity limit that constrains the number of units that can flow over that edge. Here, denotes the cost per unit flow transported directly from vertex to vertex [31]. When no units are flowing over an edge, then the arc is outlined in a dotted line.
A special source vertex produces units of a commodity that flow through the edges of the graph to be consumed by a sink vertex known as the receiving node. The rest nodes, for which there are both outflow and inflow arcs, are termed intermediate points.
The following criteria must be satisfied for any feasible flow through a network.
(1) Capacity Constraint. The flow through an edge cannot exceed the capacity of the edge , . If an edge does not exist in the network, then .
(2) Flow Conservation. Aside from the source vertex and sink vertex , each vertex must satisfy the property for which the sum of for all edges in (the flow into ) must equal the sum of for all edges (the flow out of ). This property ensures that flow is neither produced nor consumed in the network, except at and .
(3) Skew Symmetry. For consistency, the quantity represents the net flow from vertex to . This means that it must be the case that , which holds, even if both edges and exist in a directed graph (Figure 1) [32].
3.2. The Minimal Cost Flow Method
The minimum cost flow problem determines the minimum total cost under the specified flow while taking into account the arc cost. The objective function can be expressed as follows:
If is a constant, the minimum cost flow problem can be solved efficiently since it can be formulated as a linear programming problem.
Both the coal consumption and the network losses are convex functions of power flow within a power network. In minimal cost flow problems, the cost of the arc is a convex function of the flow on the arc and that is a convex cost flow problem. When is a linear function of flow , then the objective function changes into a quadratic programming problem. Not only the objective function but also the constraints include quadratic terms of . Accordingly, the cost function of unit flow can be expressed aswhere stands for or . Hence, the objective function model can be expressed as follows:
In this model, the value of depends on the nature of node , where if node is a supply node, if node is a demand node, and when node is a transshipment node [33].
In a general minimum cost flow problem, the arc is the conservation of the flow, and the flow entering an arc equals the flow leaving the arc. This assumption is reasonable in many practical application scenarios; however, the power flow will diminish when flowing through the grid due to resistance, which, in this study, is considered an issue associated with generalized flow problems. In generalized flow problems, arcs might consume or generate flow. If units of flow enter an , then units arrive at node , where is a positive multiplier associated with the arc. If , the arc is lossy, and, if , the arc is gainy. The problem becomes a general minimal cost flow problem when and . In this study, all the arcs are loss arcs due to the resistance, which means that [34].
There is no single standard algorithm that can always be used to solve convex programming problems. There have been many algorithms for solving convex quadratic programming problems, such as the Lemke method, the interiorpoint method [35, 36], the effective set method [37, 38], and an ellipsoid algorithm, for example, each having its own advantages and disadvantages.
In this study, LINGO11, very mature software widely used in various mathematical optimization problems, is used to solve the quadratic programming algorithm problems.
3.3. Multiperiod Network Flow Model
A continuous time network flow is constituted by combining each network flow chart for a single period of time. A multiperiod model is a threedimensional model with each layer corresponding to a singleperiod model (as shown in Figure 2). Each is united to form new , and, in the same way, to form a new . represents the output and represents the load. Considering each period together, a continuoustime power scheduling plan can be obtained. It is merely a schematic diagram, and the practical applications in power system may be more complex than that in Figure 2. For instance, the direction of the arcs in each layer may vary with different operating conditions. In order to simplify the analysis, we assume that the lines and transformers meet their capacity constraints.
The above structure model (Figure 2) represents a typical generalized network flow programming problem. It can be expressed assatisfied withwhere represents the dimensional arc cost vector; represents the dimensional arc flow; represents an dimensional vector injected nodes; and represents dimensional nodearc incidence matrix [39].
4. The Solution Process Diagram
The steps in the procedure of ESGD (the lowcarbon dispatching calculation process is the same) can be described as follows (as shown in Figure 3).
Step 1. Determine the load at time according to the load forecasting curve, where represents the total of time periods.
Step 2. Put as the flow at sink vertex.
Step 3. Solve the solution according to (1), satisfying the constraints (see (7) to (15)).
Step 4. If , then output ESGD scheduling results; otherwise, if , then let and go to Step 2.
5. Example Analysis
For this study, we calculated a modified IEEE 6bus power system, which contains 2 generators, 2 transformers and 5 transmission lines (as shown in Figure 4).
The parameters in the modified system are shown in Tables 1, 2, and 3. In the tables, is the resistivity of a conductor per unit length, is the length of lines, and are the resistance of lines and transformers, respectively, represents short circuit loss of a transformer, and , , and refer to the cost coefficients of the th generator (in Table 3).



The corresponding cost flow network diagram is shown in Figure 5, in which the directed diagram of the system includes 8 vertexes and 12 arcs.
In order to consider network losses, the net loss unit is turned into a standard coal consumption unit; that is, 1 MW = 0.123 t/h. In this model, the equivalent calorific value of 1 kg standard coal is 7000 kcal. The totalcost optimization model can be expressed as follows:where refers to the flow in the network; then the objective function is obviously a convex cost flow problem and a quadratic programming method can be used to solve it. When the load changes, as shown in Figure 8, the differing distribution of the minimum cost flow can be concluded and, therefore, a continuous dispatch result can be obtained.
In order to simplify the analysis, we assumed that all units run within a safe operation area. The ramp rate limit of 100 MW unit is 2 MW/min, so the ramp rate limits can meet the requirements created when the time interval is an hour.
We also assumed that the reactive power can be compensated locally and the system does not transmit reactive power. The voltage of nodes was kept within their specified ratings, and the reserve capacity of the grid was assumed to meet operational requirements. Additionally, the test was performed on an Intel (R) Pentium (R) CPU 2.13 GHz 2.0 GB RAM with LINGO11, and the average iteration times were 58. The maximum flow that can pass the network and the corresponding minimal cost distribution is shown in Figure 6. The minimum flow distribution of the minimum cost flow is shown in Figure 7.
The load and corresponding minimal coal consumption are listed in Table 4 and shown in Figure 8, respectively. Where is the time period (h), is the power of load (MW) and is the coal consumption (t/h).

Electricity and minimal coal consumption accumulation are listed in Table 5 and shown in Figure 10, respectively. refers to the accumulation of coal consumption (t) and refers to accumulation of electricity (MWh).

It can be seen that the accumulated electricity will gradually increase with the increase of time. The same trend occurs with the accumulated minimal coal consumption. Unit scheduling results are shown in Figure 9.
According to the electricalcarbon conversion relations, 2.77 tons of CO_{2} can be emitted for 1 ton of standard coal consumed. In addition, the minimum CO_{2} emissions and the accumulation curve can also be obtained (as shown in Figures 11 and 12).
The value of CO_{2} emissions and its accumulation are listed in Table 6, where and refer to CO_{2} emissions and its accumulation value, respectively.

We can see that the emissions of CO_{2} change with the power load by the same regularity, which reflects that the objectives of lowcarbon generation dispatching and energysaving generation dispatching are consistent.
In order to support our findings, we also solved this same problem by invoking the CPLEX solver in the MATLAB toolbox YALMIP [40]. The extended solver of LINGO11 includes BranchandBound solver, Global solver, and Multistart solver, while the solver in CPLEX is mainly based on interiorpoint method. If these two different methods come to a same conclusion, then the result concluded is credible. The result by CPLEX produced the same solution garnered using LINGO11, which proves the correctness of the method used in this study.
6. Conclusions
The minimal cost flow algorithm can be used to dispatch the power system and can reflect the system network topology, making it fairly easy to consider system constraints. The method is simple, rapid, and clear in terms basic concepts. The example analysis proves that the method is feasible and practical.
In order to simplify the analysis, this study only considers the energysaving scheduling of thermal units, without considering the effects of hydropower, wind power, and other renewable energy sources: although they do not consume coal and discharge CO_{2}, they do have an impact on loadflow distribution and network losses in power systems. Accordingly, the next step will be to initiate research to study the influences of renewable energy sources on ESGD in power systems.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The authors would like to acknowledge the guidance and support from their colleagues and supervisor, along with great thanks to Prof. Galen Leonhardy for editorial assistance.
References
 W. Yu and B. G. Xin, “Governance mechanism for global greenhouse gas emissions: a stochastic differential game approach,” Mathematical Problems in Engineering, vol. 2013, Article ID 312585, 13 pages, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 H. M. Ndiritu, K. Kibicho, and B. B. Gathitu, “Influence of flow parameters on capture of carbon dioxide gas by a wet scrubber,” Journal of Power Technologies, vol. 93, no. 1, pp. 9–15, 2013. View at: Google Scholar
 Y.C. Chang, T.S. Chan, and W.S. Lee, “Economic dispatch of chiller plant by gradient method for saving energy,” Applied Energy, vol. 87, no. 4, pp. 1096–1101, 2010. View at: Publisher Site  Google Scholar
 J. Kotowicz and P. Lukowicz, “Influence of chosen parameters on economic effectiveness of a supercritical combined heat and power plant,” Journal of Power Technologies, vol. 93, no. 5, pp. 323–329, 2013. View at: Google Scholar
 C.T. Cheng, S.S. Li, and G. Li, “A hybrid method of incorporating extended priority list into equal incremental principle for energysaving generation dispatch of thermal power systems,” Energy, vol. 64, no. 1, pp. 688–696, 2014. View at: Publisher Site  Google Scholar
 G.C. Liao, “A novel evolutionary algorithm for dynamic economic dispatch with energy saving and emission reduction in power system integrated wind power,” Energy, vol. 36, no. 2, pp. 1018–1029, 2011. View at: Publisher Site  Google Scholar
 R. A. Jabr, A. H. Coonick, and B. J. Cory, “A homogeneous linear programming algorithm for the security constrained economic dispatch problem,” IEEE Transactions on Power Systems, vol. 15, no. 3, pp. 930–936, 2000. View at: Publisher Site  Google Scholar
 R. A. Jabr, A. H. Coonick, and B. J. Cory, “A primaldual interior point method for optimal power flow dispatching,” IEEE Transactions on Power Systems, vol. 17, no. 3, pp. 654–662, 2002. View at: Publisher Site  Google Scholar
 W. Ongsakul and N. Petcharaks, “Unit commitment by enhanced adaptive Lagrangian relaxation,” IEEE Transactions on Power Systems, vol. 19, no. 1, pp. 620–628, 2004. View at: Publisher Site  Google Scholar
 H. Zhong, Q. Xia, Y. Wang, and C. Kang, “Dynamic economic dispatch considering transmission losses using quadratically constrained quadratic program method,” IEEE Transactions on Power Systems, vol. 28, no. 3, pp. 2232–2241, 2013. View at: Publisher Site  Google Scholar
 N. Sinsuphan, U. Leeton, and T. Kulworawanichpong, “Optimal power flow solution using improved harmony search method,” Applied Soft Computing Journal, vol. 13, no. 5, pp. 2364–2374, 2013. View at: Publisher Site  Google Scholar
 J. J. Hargreaves and B. F. Hobbs, “Commitment and dispatch with uncertain wind generation by dynamic programming,” IEEE Transactions on Sustainable Energy, vol. 3, no. 4, pp. 724–734, 2012. View at: Publisher Site  Google Scholar
 J. Soares, M. Silva, T. Sousa, Z. Vale, and H. Morais, “Distributed energy resource shortterm scheduling using Signaled Particle Swarm Optimization,” Energy, vol. 42, no. 1, pp. 466–476, 2012. View at: Publisher Site  Google Scholar
 C. Yaşar and S. Özyön, “Solution to scalarized environmental economic power dispatch problem by using genetic algorithm,” International Journal of Electrical Power and Energy Systems, vol. 38, no. 1, pp. 54–62, 2012. View at: Publisher Site  Google Scholar
 M. H. Gomes and J. T. Saraiva, “A market based active/reactive dispatch including transformer taps and reactor and capacitor banks using simulated Annealing,” Electric Power Systems Research, vol. 79, no. 6, pp. 959–972, 2009. View at: Publisher Site  Google Scholar
 P. Somasundaram and K. Kuppusamy, “Application of evolutionary programming to security constrained economic dispatch,” International Journal of Electrical Power and Energy Systems, vol. 27, no. 56, pp. 343–351, 2005. View at: Publisher Site  Google Scholar
 S. Pothiya, I. Ngamroo, and W. Kongprawechnon, “Ant colony optimisation for economic dispatch problem with nonsmooth cost functions,” International Journal of Electrical Power and Energy Systems, vol. 32, no. 5, pp. 478–487, 2010. View at: Publisher Site  Google Scholar
 M. Sharma, M. Pandit, and L. Srivastava, “Reserve constrained multiarea economic dispatch employing differential evolution with timevarying mutation,” International Journal of Electrical Power and Energy Systems, vol. 33, no. 3, pp. 753–766, 2011. View at: Publisher Site  Google Scholar
 J. S. Alsumait, J. K. Sykulski, and A. K. AlOthman, “A hybrid GAPSSQP method to solve power system valvepoint economic dispatch problems,” Applied Energy, vol. 87, no. 5, pp. 1773–1781, 2010. View at: Publisher Site  Google Scholar
 J. Cai, Q. Li, L. Li, H. Peng, and Y. Yang, “A hybrid CPSOSQP method for economic dispatch considering the valvepoint effects,” Energy Conversion and Management, vol. 53, no. 1, pp. 175–181, 2012. View at: Publisher Site  Google Scholar
 C. C. A. Rajan, “A solution to the economic dispatch using EP based SA algorithm on large scale power system,” International Journal of Electrical Power & Energy Systems, vol. 32, no. 6, pp. 583–591, 2010. View at: Publisher Site  Google Scholar
 M. Basu, “Hybridization of bee colony optimization and sequential quadratic programming for dynamic economic dispatch,” International Journal of Electrical Power and Energy Systems, vol. 44, no. 1, pp. 591–596, 2013. View at: Publisher Site  Google Scholar
 A. R. L. Oliveira, S. Soares, and L. Nepomuceno, “Optimal active power dispatch combining network flow and interior point approaches,” IEEE Transactions on Power Systems, vol. 18, no. 4, pp. 1235–1240, 2003. View at: Publisher Site  Google Scholar
 J. Zhu and J. A. Momoh, “Multiarea power systems economic dispatch using nonlinear convex network flow programming,” Electric Power Systems Research, vol. 59, no. 1, pp. 13–20, 2001. View at: Publisher Site  Google Scholar
 C. Li, Y. Liu, Y. Cao, Y. Tan, C. Xue, and S. Tang, “Consistency evaluation of lowcarbon generation dispatching and energysaving generation dispatching,” Proceedings of the Chinese Society of Electrical Engineering, vol. 31, no. 31, pp. 94–101, 2011. View at: Google Scholar
 S. S. Reddy, B. K. Panigrahi, R. Kundu, R. Mukherjee, and S. Debchoudhury, “Energy and spinning reserve scheduling for a windthermal power system using CMAES with mean learning technique,” International Journal of Electrical Power and Energy Systems, vol. 53, no. 1, pp. 113–122, 2013. View at: Publisher Site  Google Scholar
 C.C. Kuo, “Generation dispatch under large penetration of wind energy considering emission and economy,” Energy Conversion and Management, vol. 51, no. 1, pp. 89–97, 2010. View at: Publisher Site  Google Scholar
 A. Immanuel Selvakumar, “Enhanced crossentropy method for dynamic economic dispatch with valvepoint effects,” International Journal of Electrical Power and Energy Systems, vol. 33, no. 3, pp. 783–790, 2011. View at: Publisher Site  Google Scholar
 R. AzizipanahAbarghooee, T. Niknam, A. Roosta, A. R. Malekpour, and M. Zare, “Probabilistic multiobjective windthermal economic emission dispatch based on point estimated method,” Energy, vol. 37, no. 1, pp. 322–335, 2012. View at: Publisher Site  Google Scholar
 M. A. Abido, “Multiobjective particle swarm optimization for environmental/economic dispatch problem,” Electric Power Systems Research, vol. 79, no. 7, pp. 1105–1113, 2009. View at: Publisher Site  Google Scholar
 F.R. Xie and R.A. Jia, “Nonlinear fixed charge transportation problem by minimum cost flowbased genetic algorithm,” Computers & Industrial Engineering, vol. 63, no. 4, pp. 763–778, 2012. View at: Publisher Site  Google Scholar
 G. T. Heineman, G. Pollice, and S. Selkow, Algorithms in a Nutshell, O’Reilly Media, Sebastopol, Calif, USA, 2008.
 F. S. Hillier and G. J. Lieberman, Introduction to Operations Research, McGrawHill, 9th edition, 2010.
 R. K. Ahuja, T. L. Magnanti, and J. Orlin, Network Flows: Theory, Algorithms and Applications, Prentice Hall, Englewood Cliffs, NJ, USA, 1993. View at: MathSciNet
 F. Curtis and J. Nocedal, “Steplength selection in interiorpoint methods for quadratic programming,” Applied Mathematics Letters, vol. 20, no. 5, pp. 516–523, 2007. View at: Publisher Site  Google Scholar  MathSciNet
 M. D’Apuzzo and M. Marino, “Parallel computational issues of an interior point method for solving large boundconstrained quadratic programming problems,” Parallel Computing: Theory and Applications, vol. 29, no. 4, pp. 467–483, 2003. View at: Publisher Site  Google Scholar  MathSciNet
 W. Cheng, Z. Chen, and D.H. Li, “An active set truncated Newton method for largescale bound constrained optimization,” Computers & Mathematics with Applications, vol. 67, no. 5, pp. 1016–1023, 2014. View at: Publisher Site  Google Scholar
 M.T. Yu, T.Y. Lin, and C. Hung, “Activeset sequential quadratic programming method with compact neighbourhood algorithm for the multipolygon mass production cuttingstock problem with rotatable polygons,” International Journal of Production Economics, vol. 121, no. 1, pp. 148–161, 2009. View at: Publisher Site  Google Scholar
 Z.A. Zhang and X.G. Cai, “Power purchase plan using minimal cost flow,” Computer Modelling and New Technologies, vol. 18, no. 2, pp. 281–285, 2014. View at: Google Scholar
 J. Löfberg, “YALMIP: a toolbox for modeling and optimization in MATLAB,” in Proceedings of the IEEE International Symposium on Computer Aided Control System Design, pp. 284–289, Taipei, Taiwan, September 2004. View at: Google Scholar
Copyright
Copyright © 2015 Zhan’an Zhang and Xingguo Cai. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.