ABSTRACT
Background and Objective: By increasing the penetration level of renewable energies on the generation side and the emergence of new variable load on the demand side, stochastic analysis of the conventional security constrained unit commitment problem has become more important for the secure optimal operation of the electricity market. Today, the increasing utilization of plug-in electric vehicles, which consume electricity rather than fossil fuel for driving, offers new opportunities and challenges to the operation of electric power system. By appropriate managing and day-ahead scheduling of these types of vehicles, challenges can be replaced by opportunities for the power system operation and planning. Methodology: In this study, a new method is proposed for stochastic security-constrained unit commitment problem in the presence of wind power generations and plug-in electric vehicles. The method enjoys the advantages of conventional scenario-based approaches and mitigates their barriers by using interior point optimization techniques. The proposed algorithm is implemented on two standard networks: A 6-bus test system and a large-scale 118-bus system. Results: This study demonstrate the accuracy and efficiency of the proposed method, especially in large-scale power systems with different types of uncertainties. Conclusion: By increasing the speed of simulation, more uncertainties can be modeled and therefore, more realistic and accurate results can be obtained.
PDF Abstract XML References Citation
How to cite this article
DOI: 10.3923/jas.2016.189.200
URL: https://scialert.net/abstract/?doi=jas.2016.189.200
INTRODUCTION
Independent System Operators (ISO) accomplish Security-Constrained Unit Commitment (SCUC) in order to design economic and secure generation schedules for the daily electricity market. They apply the comprehensive market information provided by participants, like performance of generating units, generation offer and demand order, capability of transmission lines and so on1. The SCUC provides an economically workable unit commitment which is practically acceptable. The corresponding market participants are provided by SCUC-based generation dispatch. Commonly, a satisfactory SCUC solution could be achieved as long as daily market is proper and robust1.
A large volume of literature is devoted to SCUS solving methods. In Wang et al.2 and Fu et al.3, Augmented Lagrangian Relaxation (ALR) method and an algorithm according to Benders Decomposition (BD) technique have been proposed, respectively. Constraint reduction and computation speed increase were investigated4. Non-convex SCUC and a new method based on quadratic programming were proposed5. Effectiveness of two algorithms extensively applied in SCUC, i.e., Mixed-Integer Programing (MIP) and Lagrangian Relaxation (LR), were evaluated6.
In recent years, uncertainties of the components of power systems have increased drastically. Two most important sources of power system uncertainties are wind generations and price-based loads. Wind energy usage has grown increasingly in power systems in recent years. However, inherent probabilistic and non-dispatchable characteristics of wind energy could lead to problem in power system features such as frequency, voltage and generation sufficiency. In addition, price-based load is considered another uncertain parameter in power system operation and planning. Therefore, regarding Stochastic SCUC as a tool for optimizing power system operation is completely inevitable7.
Stochastic SCUC problem which is closer to the real performance of large power systems than deterministic SCUC has attracted great interest of researchers. Distribution generations widely utilized in power systems were taken into consideration for SCUC8. Scenario-based methods are widely applied to solve SCUC with uncertainties; however, they result in escalating the computational burden. In Mehrtash et al.9, stochastic SCUC was solved according to point estimation method and BD, including significant reduction in the No. of scenarios and calculation time.
Plug-in Electric Vehicles (PEVs) are considered hourly spread and mobile demand in power systems. The accumulated storage ability of PEVs could alter hourly generation portfolio and decrease the performance costs of the grid10.
Due to the importance of PEV concept and application, a large number of recent studies have concentrated on their related subjects. In Tomic and Kempton11, the economic capability of PEVs for contributing to regulation services was evaluated. The integration of PEV in power systems was considered12 and their electricity market problem was investigated13.
Similar to wind generation, PEV is an uncertain source for grid. Therefore, its random behavior is investigated in different power system problems, especially Stochastic SCUC. The PEVs random behavior in conjunction with stochastic characteristic of wind energy was evaluated in stochastic SCUC14. A cooperation model and constraints of PEVs and wind generation in the presence of battery storage were proposed15. In the last reference, the model included the distribution pattern of user trips and evaluated the dynamic process of stored energy. In Krad and Gao16, the capability of PEVs in providing contingency reserves was investigated. Power system operation and control from the aspect of being influenced by contributing battery-based energy storage transportation by railway transportation network was studied17.
New Interior Point Optimization Techniques (IPOPT) are capable to solve large-scale Mixed-Integer Non-Linear Programing (MINLP) problems18. In this study, a new method is proposed for Stochastic SCUC problem in the presence of wind power generations and PEVs. This method enjoys the advantages of conventional scenario-based approaches and mitigates their barriers using IPOPT techniques. Instead of linearizing, in the proposed method, non-linear problem is solved via IPOPT techniques. Moreover, since the number of Benderss cuts is reduced, the number of master problem solving iterations in BD is lessened and the simulation time is decreased considerably.
The rest of this study is organized as follows: Interior point method is presented in Section 2. Problem formulation and solution methodology are presented in Sections 3 and 4, respectively. The result of the case studies is given in Section 5. Finally, the conclusion is presented in Section 6.
MATERIALS AND METHODS
Interior point method: Advantages of IPOPT for solving large-scale problems will be more understandable when it is compared with Simplex method. Simplex, introduced by Dantzig19, is a method for solving Linear Programing (LP). Algorithms based on Simplex method search vertices of the feasible region; therefore, for large-scale problems, in which the number of vertices is high, the solution time is considerable. As a result, Simplex is categorized among exponential time algorithms. In other words, time needed by Simplex for solving an LP problem is an exponential function of the problem scale20.
Interior point approach was proposed by Karmarkar21. Its algorithms search path is inside the feasible region; in fact, central points should be investigated (Fig. 1). Thus, the problem scale is not too critical and it can be considered among polynomial time algorithms20. A complete comparison between Simplex and IPOPT was proposed22.
At first, interior point approach was proposed only for LP problems. However, the idea of using this method for Non-Linear Programing (NLP) became popular immediately. Recently, COUENNE, which is a new algorithm based on IPOPT, has been proposed by Belotti et al.23. The COUENNE is an open-source solver for non-convex MINLPs18. In order to overcome MINLP challenges, COUENNE uses disjunctive cuts and a tight branch-and-bound approach24. This solver was implemented in General Algebraic Modeling System (GAMS) as a non-commercial solver. In this study, GAMS/COUENNE is used to solve Stochastic SCUC as an MINLP problem.
Stochastic SCUC problem formulation: Concepts of corrective and preventive actions and definition of Stochastic SCUC problem have been mentioned in authors previous study9. In this section, only the formulation of the problem is discussed.
Objective function: Objective function, Eq. 1 consists of generation cost, startup and shutdown costs of thermal and hydro units, operation cost of PEV fleets, availability cost for providing spinning reserve in scenarios and expected cost of corrective actions in scenarios for mitigating uncertainties.
Availability cost is the payment to generators that provides reserves as a corrective action in response to uncertainties. The providing reserve of generators is limited by their ramp up/down capability25. The availability cost is considered one-third of the marginal cost of a generator. Thermal units are formulated as non-quick start units; therefore, their scenario commitment status is similar to the base case. Thus, the startup/shutdown costs should not be introduced in the scenarios14.
![]() | |
Fig. 1: | Search path of IPOPT method |
![]() | (1) |
Constraints for the base case: In this stage, the values related to the wind, load, PEV, etc. are fixed to their forecasted values. Thus, these formulas are deterministic.
The system power balance constraint is presented in Eq. 2. A detailed formulation of thermal and hydro units constraints is available25,26. The base case PEV fleet constraints are mentioned in Eq. 3-10. The net hourly absorbed/delivered energy from/to grid by PEV battery is given in Eq. 3. The charging efficiency, ηv, is defined as the ratio of the energy stored in the PEV battery to the energy drawn from the grid. Equation 4 represents the hourly charge/discharge/idle status of fleets. Obviously, these status are mutually exclusive. Charge/discharge power limitations are given in Eq. 5 and 6. Equation 7 represents the hourly energy balance in PEV batteries. Parameter Nv,t denotes connectivity status of PEV fleet to the grid. When PEV fleet is connected to the grid , Nv,t = 1, it is in the either charge or discharge or idle mode. Once PEV fleet is decoupled from the grid, Nv,t = 0, the charge/discharge power will be zero according to Eq. 4-614. The capacity limitation of the PEV batteries is mentioned in Eq. 8 and 9. The cost function of the PEV batteries, which is considered as a non-linear function, is given in Eq. 10. The cost function parameters depend on the depth of discharge and cycles to failure of the battery. A complete discussion about PEV battery technologies can be find27. Network constraint based on DC power flow is shown in Eq. 11 and 12:
![]() | (2) |
![]() | (3) |
![]() | (4) |
![]() | (5) |
![]() | (6) |
![]() | (7) |
![]() | (8) |
![]() | (9) |
![]() | (10) |
![]() | (11) |
![]() | (12) |
Constraints for each scenario: At this stage, which is the stochastic part of the formulation, the values related to wind, load, PEV, etc. are derived from Monte Carlo scenarios. The initial scenarios are generated and then reduced to final scenarios via scenario reduction techniques. Afterwards, the final scenarios are used in the following formulation.
The systems power balance constraint for each scenario is given in Eq. 13. Scenario constraint related to the thermal and hydro units is mentioned28. The PEV scenario constraint is formulated in Eq. 14-21. The equations are similar to those described in the base case section, with the only difference that the deterministic forecasted values for wind, load, PEV, etc. are replaced by their value obtained from the final Monte Carlo scenarios. Parameter , one of the uncertain inputs, denotes ratio of the number of PEVs in scenarios to the number of base case PEVs. The corrective action for the scenarios is enforced by Eq. 22, while the hourly cost of corrective action,
, is included in the objective function14. Equation 22 shows the correlation between generation of a unit in the base case, Pi,t and its generation in scenarios,
. Network constraint based on DC power flow for each scenario is presented in Eq. 23 and 24:
![]() | (13) |
![]() | (14) |
![]() | (15) |
![]() | (16) |
![]() | (17) |
![]() | (18) |
![]() | (19) |
![]() | (20) |
![]() | (21) |
![]() | (22) |
![]() | (23) |
![]() | (24) |
Solution methodology: The scenario-based Stochastic SCUC problem in the presence of PEVs, modeled in Eq. 1-24, is a non-linear, non-convex, large-scale and polynomial-time hard (NP-hard) problem.
![]() | |
Fig. 2: | Flowchart of the proposed method |
Solving such a problem for large-scale power systems would be intractable without decomposition9.
The flowchart of the proposed method is shown in Fig. 2. According to the flowchart, Stochastic SCUC problem is solved at three stages.
At the first stage, Monte Carlo scenarios are generated and reduced to the final scenarios. Due to the high simulation time of Stochastic SCUC problem, scenario reduction must be done before solving the main problem26.
Afterwards, BD is used for solving the problem14. Although, some new IPOPT solvers, such as COUENNE, are able to solve this problem without BD, it can be concluded from the simulation of the sample case studies that the simulation time decreases significantly using BD.
At the second stage, a deterministic SCUC problem, master problem, is solved using IPOPT.
At the third stage, the results of the master problem should be evaluated in the scenario security check subproblem for each final scenario. Every violation will be referred to the master problem, second stage, for adding new constraints based on the generated Benders cuts. Finally, all the constraints are met in all final scenarios.
This algorithm has two main improvement compared to the previous algorithms proposed for Stochastic SCUC problem14. First, the master problem is solved as an MINLP by IPOPT, instead of a, MIP, by Simplex based methods. As mentioned in Section 2, interior point solvers, such as COUENNE, include better features to handle large-scale problems. The ability of these solvers for solving MINLPs is illustrated in the next section by a case study simulation. Second, due to the ability of the IPOPT solvers, UC and network evaluation check are solved at a single stage. However, in previous algorithms, UC is solved at stage one and the network evaluation check is solved as a BD subproblem at stage two. Combination of the UC and network evaluation check at a single stage can improve the simulation time, due to a considerable amount of the total simulation time belongs to master problem. Reducing the subproblem stages leads to fewer benders cuts; therefore, fewer iterations of master problem solving are needed. This improvement can be concluded from the simulation results of the case studies in the next section.
RESULTS AND DISCUSSION
In this section, two case studies, a 6-bus test system and a 118-bus power system, are studied to demonstrate the validity of the proposed method for Stochastic SCUC problem in the presence of PEVs. The proficiency of the proposed algorithm with IPOPT is illustrated in comparison to the conventional Simplex-based methods.
![]() | |
Fig. 3: | One-line diagram of 6-bus test system |
Table 1: | Thermal unit characteristic |
![]() | |
Table 2: | Transmission line characteristics |
![]() | |
6-bus system: First, the proposed method is implemented on the 6-bus test system shown in Fig. 3. The system includes three thermal generators, one wind farm, seven transmission lines, three load points and five PEV fleets. The loads located on buses 3, 4 and 5 consume 20, 40 and 40% of the total load, respectively.
The parameters of the thermal units and transmission lines are shown in Table 1 and 2, respectively. Table 3 shows the hourly forecasted load and wind power generation of the system. In order to show the capability of the proposed method at high volatility of power generation, the penetration level of the wind power (ratio of the wind power generation to the total generation) is increased to 40% at some hours. As a result, a large amount of generation is probabilistic at those hours. The hourly penetration level of the wind power is given in Fig. 4.
The PEV fleet characteristics and their travel pattern are presented in Table 4 and 5, respectively. Charging efficiency of each fleet, which is the ratio of energy stored in the battery to the energy drawn from the grid, is assumed 85%14. The annular driving distance by a PEV fleet is 12000 mile and its average per day is 32.88 mile29,30. The energy required by a single PEV is about 9 kWh/day with the average 3.65 mile/kWh11.
Table 3: | Hourly forecasted load and wind power for 6-bus system |
![]() | |
Therefore, the hourly energy required by fleet one to five is 7.65, 9.00, 2.25, 7.20 and 4.50 MWh, respectively14.
Four cases are studied on the 6-bus test system:
• | Deterministic SCUC with Simplex (conventional method) |
• | Deterministic SCUC with IPOPT (proposed method) |
• | Stochastic SCUC with Simplex (conventional method) |
• | Stochastic SCUC with IPOPT (proposed method) |
Case 1: In this case, the uncertainties are modeled with their forecasted values. Due to the ignorance of the uncertainties, there is no need to generate either Monte Carlo scenarios or scenario security check subproblem.
![]() | |
Fig. 4: | Penetration level of wind power in 6-bus test system |
![]() | |
Fig. 5: | Generation scheduling of thermal units in case 1 |
Table 4: | PEV fleet characteristics for 6-bus system |
![]() | |
Table 5: | PEV fleet travel characteristics for 6-bus system |
![]() | |
The deterministic SCUC problem is solved as an MIP problem using CPLEX solver of GAMS18. Quadratic cost functions are linearized in order to solve the problem using a Simplex solver, i.e., CPLEX. The BD is used to decompose the problem into a master UC problem and a network security check subproblem as explained14. The operation cost is 97091.741 $. Generation scheduling of thermal units and power dispatch of PEVs is demonstrated in Fig. 5 and Table 6, respectively.
![]() | |
Fig. 6: | Generation scheduling of thermal units in case 2 |
Table 6: | Power dispatch of PEVs in case 1 |
![]() | |
Table 7: | Power dispatch of PEVs in case 2 |
![]() | |
Table 8: | Probability of final scenarios |
![]() | |
Case 2: In this case, the deterministic SCUC problem is solved as an MINLP problem using COUENNE solver of GAMS18. As mentioned in Section 2, due to the capability of this solver, there is no need to either linearize or decompose the problem. The operation cost is 97067.693 $. Power dispatch of PEVs and generation scheduling of thermal units are given in Table 7 and Fig. 6, respectively.
As given in Fig. 5 and 6, using the proposed method, the third thermal unit schedule at h 6 is on, while it is off in case 1. Moreover, the total operation cost is decreased by 24.048 $, which is an improvement in the schedule that leads to less operation cost. This decrease in the operation cost is the result of solving the problem as an MINLP without accepting the errors of linearizing the cost functions.
Case 3: For considering uncertainties, wind power generation, load points and No. of PEVs in fleets are modeled by normal distribution functions. Their mean is set to their forecasted values and their standard deviations are assumed as 10, 5 and 10% of the forecasted values, respectively. Stochastic SCUC problem is solved as an MIP problem using CPLEX solver of GAMS. Quadratic cost functions are linearized in order to solve the problem using Simplex solver, i.e., CPLEX. The BD is used to decompose the problem into a master UC problem and two subproblems, a network security check subproblem and a scenario security check subproblem14. To run the Monte Carlo simulation, initially, 10000 random scenarios are generated and then, they are reduced to 15 final scenarios using the fast backward/forward technique of SCENRED library of GAMS software GAMS/SCENRED (http://www.gams.com/). The probability of each final scenario is given in Table 8. The total operation cost is 93311.122 $. Generation scheduling of thermal units and power dispatch of PEVs are shown in Fig. 7 and Table 9, respectively.
Case 4: In this case, Stochastic SCUC problem is solved as an MINLP problem using COUENNE solver of GAMS. As mentioned in Section 2, due to the capability of this solver, there is no need to linearize the problem.
![]() | |
Fig. 7: | Generation scheduling of thermal units in case 3 |
![]() | |
Fig. 8: | Generation scheduling of thermal units in case 4 |
Table 9: | Power dispatch of PEVs in case 3 |
![]() | |
Table 10: | Power dispatch of PEVs in case 4 |
![]() | |
Table 11: | Total operation cost and simulation time for all cases |
![]() | |
The problem is solved using BD based on the proposed algorithm as shown in Fig. 2. The operation cost is 93282.094 $. Power dispatch of PEVs and generation scheduling of thermal units are given in Table 10 and Fig. 8, respectively. The operation cost and simulation time of all the four cases are presented in Table 11.
By comparing Fig. 7 and 8, we understand that, in contrast to case 3, in case 4, the third thermal unit is off at 3 h.
Table 12: | Wind farms hourly generation power for 118-bus system |
![]() | |
Table 13: | PEV fleet characteristics for 118-bus system |
![]() | |
Table 14: | Total operation cost and simulation time for all cases of 118-bus system |
![]() | |
Moreover, the total operation cost is decreased by 29.028 $. Therefore, as shown in Table 11, it can be concluded that solving Stochastic SCUC problem as an MINLP with IPOPT leads to less final operation cost due to the elimination of linearizing errors. Besides, the simulation time decreases due to mitigating the iterations of master UC problem by decreasing the number of the required Benders cuts.
118-bus system: For scenario-based Stochastic SCUC problem in large-scale power systems, the computational burden is drastically intensified and efficient reduction of CPU time is necessary. Moreover, the advantage of the proposed method in a large-scale power system with numerous uncertainties can be better comprehended.
In the second part of the case studies, the proposed method is implemented on a 118-bus test system, which includes 54 thermal generators, 9 transformers, 186 transmission lines and 91 load points. The system is modified by adding three wind farms and five PEV fleets. The location and the hourly generation of the three wind farms are given in Table 12. Totally, 100000 PEVs are divided into 5 fleets. The number of PEVs in fleet 1 through 5 is 34000, 20000, 10000, 16000 and 20000, respectively14. The fleet characteristics are shown in Table 13. The PEVs travel characteristics are similar to those in the previous case study (Table 5).
Four cases are implemented on this system and the results are given in Table 14. As shown in this table, the proposed method based on IPOPT reduces both operation cost and simulation time for the Stochastic SCUC of such large-scale power systems.
CONCLUSION
In this study, a new method based on interior point optimization has been proposed for scenario-based Stochastic SCUC in the presence of PEVs. Whereas, the proposed method enjoys the advantages of the scenario-based method, it also mitigates the obstacles of this method, e.g., time consuming simulation due to computational burden and linearizing errors. By increasing the speed of simulation, more uncertainties can be modeled and, therefore, more realistic and accurate results can be obtained. The proposed algorithm has been implemented on two standard networks: A 6-bus test system and a large-scale 118-bus system. These case studies demonstrate the accuracy and efficiency of the proposed method especially in the large-scale power systems with different type of uncertainties.
ACKNOWLEDGMENT
This research was supported by Young Researchers and Elite Club, Dariun Branch, Islamic Azad University, Dariun, Fars, Iran. We thank Hamid Dehghanfar, Dean of Young Researchers and Elite Club, for assistance with supporting our study.
APENDIX NOMENCLATURE![]() |
REFERENCES
- Wang, S.J., S.M. Shahidehpour, D.S. Kirschen, S. Mokhtari and G.D. Irisarri, 1995. Short-term generation scheduling with transmission and environmental constraints using an augmented Lagrangian relaxation. IEEE Trans. Power Syst., 10: 1294-1301.
CrossRefDirect Link - Fu, Y., M. Shahidehpour and Z. Li, 2006. AC contingency dispatch based on security-constrained unit commitment. IEEE Trans. Power Syst., 21: 897-908.
CrossRefDirect Link - Zhai, Q., X. Guan, J. Cheng and H. Wu, 2010. Fast identification of inactive security constraints in SCUC problems. IEEE Trans. Power Syst., 25: 1946-1954.
CrossRefDirect Link - Shafie-Khah, M., M.P. Moghaddam and M.K. Sheikh-El-Eslami, 2011. Unified solution of a non-convex SCUC problem using combination of modified Branch-and-Bound method with Quadratic Programming. Energy Conv. Manage., 52: 3425-3432.
CrossRefDirect Link - Fu, Y., Z. Li and L. Wu, 2013. Modeling and solution of the large-scale security-constrained unit commitment. IEEE Trans. Power Syst., 28: 3524-3533.
CrossRefDirect Link - Khodaei, A., M. Shahidehpour and S. Bahramirad, 2011. SCUC with hourly demand response considering intertemporal load characteristics. IEEE Trans. Smart Grid, 2: 564-571.
CrossRefDirect Link - Kargarian, A. and Y. Fu, 2014. System of systems based security-constrained unit commitment incorporating active distribution grids. IEEE Trans. Power Syst., 29: 2489-2498.
CrossRefDirect Link - Mehrtash, M., M. Raoofat, M. Mohammadi and H. Zareipour, 2016. Fast stochastic security-constrained unit commitment using point estimation method. Int. Trans. Electr. Energy Syst., 26: 671-688.
CrossRefDirect Link - Lund, H. and W. Kempton, 2008. Integration of renewable energy into the transport and electricity sectors through V2G. Energy Policy, 36: 3578-3587.
CrossRefDirect Link - Tomic, J. and W. Kempton, 2007. Using fleets of electric-drive vehicles for grid support. J. Power Sources, 168: 459-468.
CrossRefDirect Link - Kempton, W. and J. Tomic, 2005. Vehicle-to-grid power fundamentals: Calculating capacity and net revenue. J. Power Sources, 144: 268-279.
CrossRefDirect Link - Kempton, W. and J. Tomic, 2005. Vehicle-to-grid power implementation: From stabilizing the grid to supporting large-scale renewable energy. J. Power Sources, 144: 280-294.
CrossRefDirect Link - Khodayar, M.E., L. Wu and M. Shahidehpour, 2012. Hourly coordination of electric vehicle operation and volatile wind power generation in SCUC. IEEE Trans. Smart Grid, 3: 1271-1279.
CrossRefDirect Link - Shao, C., X. Wang, X. Wang, C. Du, C. Dang and S. Liu, 2014. Cooperative dispatch of wind generation and electric vehicles with battery storage capacity constraints in SCUC. IEEE Trans. Smart Grid, 5: 2219-2226.
CrossRefDirect Link - Krad, I. and D.W. Gao, 2013. Impact of phev on reserve scheduling: A milp-scuc method. Proceedings of the North American Power Symposium (NAPS), September 22-24, 2013, Manhattan, KS., pp: 1-6.
CrossRef - Sun, Y., Z. Li, M. Shahidehpour and B. Ai, 2015. Battery-based energy storage transportation for enhancing power system economics and security. IEEE Trans. Smart Grid, 6: 2395-2402.
CrossRefDirect Link - Dantzig, G.B., 1951. Maximization of a Linear Function of Variables Subject to Linear Inequalities. In: Activity Analysis of Production and Allocation, Koopmans, T.C. (Ed.). Chapter 21, John Wiley and Sons, New York, USA., pp: 339-347.
Direct Link - Karmarkar, N., 1984. A new polynomial-time algorithm for linear programming. Combinatorica, 4: 373-395.
CrossRefDirect Link - Szabo, Z. and M. Kovacs, 2003. On interior-point methods and simplex method in linear programming. Analele Stiintifice ale Universitatii Ovidius Constanta, 11: 155-162.
Direct Link - Belotti, P., J. Lee, L. Liberti, F. Margot and A. Wachter, 2009. Branching and bounds tighteningtechniques for non-convex MINLP. Optim. Methods Software, 24: 597-634.
CrossRefDirect Link - Belotti, P., 2012. Disjunctive Cuts for Nonconvex MINLP. In: Mixed Integer Nonlinear Programming, Lee, J. and S. Leyffer (Eds.). Springer, New York, USA., ISBN: 9781461419273, pp: 117-144.
Direct Link - Fu, Y., M. Shahidehpour and Z. Li, 2005. Security-constrained unit commitment with AC constraints. IEEE Trans. Power Syst., 20: 1001-1013.
CrossRefDirect Link - Wang, J., M. Shahidehpour and Z. Li, 2008. Security-constrained unit commitment with volatile wind power generation. IEEE Trans. Power Syst., 23: 1319-1327.
CrossRefDirect Link - Li, T. and M. Shahidehpour, 2005. Price-based unit commitment: A case of lagrangian relaxation versus mixed integer programming. IEEE Trans. Power Syst., 20: 2015-2025.
CrossRefDirect Link - Saber, A.Y. and G.K. Venayagamoorthy, 2011. Plug-in vehicles and renewable energy sources for cost and emission reductions. IEEE Trans. Ind. Electron., 58: 1229-1238.
CrossRefDirect Link - Roe, C., A. Meliopoulos, J. Meisel and T. Overbye, 2008. Power system level impacts of plug-in hybrid electric vehicles using simulation data. Proceedings of the IEEE Energy 2030 Conference, 2008, ENERGY, November 17-18, 2008, Atlanta, GA., pp: 1-6.
CrossRef