In the optimal flow problem, the generation and transmission of electric energy is optimised, taking into account the active and reactive power generation limits, demand requirements, bus voltage limits, and network flow limits. Due to the non-linear power flow constraints, OPF is formulated as a non-convex optimisation problem that is generally difficult to solve.

The problem was first formulated in 1962 and a large number of optimization algorithms and relaxations have been proposed since then. A popular approximation, the directed-current optimal power flow (DCOPF), is based on the linear programming problem that is obtained through the linearisation of the power flow equations. While DCOPF is useful in a wide variety of applications, a solution of DCOPF may not satisfy the non-linear power flow equations and hence the resulting solution may be infeasible and may be of limited utility. Other types of algorithms that were studied for the OPF include Newton Raphson, quadratic programming, nonlinear programming, Lagrangian relaxation, interior point methods, genetic algorithms, and other heuristics. Although some of these algorithms can handle large-scale networks most them can only compute a local optimal usually without assurance on the quality of the solution.  That is because most of the algorithms rely on the Karush-Kuhn-Tucker necessary conditions of optimality, which cannot even guarantee a locally optimal solution, in the non-convex problem. Alternatively, the OPF can be formulated as a quadratically constrained quadratic program. There, convex relaxations within second-order cone (SOCP) programming and semidefinite programming (SDP) can be applied. In contrast to the other proposed approaches, convex relaxations make it possible to check if a solution is globally optimal. If the solution is not optimal, the relaxations provide a lower bound and hence a bound on how far any feasible solution is from optimality. In particular, Bai et al. proposed a semidefinite programming relaxation for the ACOPF for general networks, which makes it possible to find globally optimal solutions for several well-known instances. More recently, the moments and sum-of-square decomposition have been used to build hierarchies of improving SDP relaxations for a polynomial programming formulation of ACOPF. To overcome the computational complexity of using SDP and polynomial programming, graph sparsity has been exploited to simplify the SDP relaxation of the OPF.

To further improve the scalability of SDP relaxations, ADMM-based computation can be used to solve sparse, large-scale SDPs. Alternatively, cheaper hierarchies are being investigated based on LP and SOCP relaxations. In the future, a combination of hierarchies mixing constraints from different cones may be envisioned.

Another issue to address is development of techniques to certify infeasibility of optimal power flow instances.

From an industrial point of view, dealing with incomplete data is one of the issues models and tools have to address. Aggregations of industrial data may lead to physically non-meaningful models, since some section of the power network are not represented in the data.



J. Carpentier, Contribution to the economic dispatch problem, Bulletin Society Francaise Electriciens, vol. 3, no. 8, pp. 431-447, 1962.

X. Bai, H. Wei, K. Fujisawa and Y. Wang, Semidefinite programming for optimal power flow problems, International Journal of Electric Power & Energy Systems, vol. 30, no. 6-7, pp. 383-392, 2008.

J. Lavaei and S. H. Low. "Zero duality gap in optimal power flow problem." IEEE Transactions on Power Systems, vol. 27, no 1, pp. 92-107, 2012.

S.H. Low . "Convex relaxation of optimal power flow, part I: Formulations and equivalence." IEEE Transactions on Control of Network Systems, vol. 1, no. 1, pp. 15-27, 2014.

X. Kuang, approximating the ACOPF problem with hierarchy of SOCP problems, 2015.

Madani, Kalbat, Lavaei, ADMM for Sparse Semidefinite Programming with Applications to Optimal Power Flow problem.

Molzahn (2013), Sufficient conditions for power flow in-solvability with applications to voltage stability margins.

Coffrin, Hijazi (2015), The QC relaxation: Theoreticl and Computational Results on Optimal Power Flow, IEEE Trans. Power Systems, Vol. PP, Issue 99



OPF specific software

1) PowerTools

PowerTools ( is a mathematical programming library for Power Systems optimization. PowerTools is an open source c++ mathematical modeling platform linking to state-of-the-art optimization solvers. It implements problems such as Optimal Power Flow, and Optimal Transmission Switching, while benefiting from cutting-edge research results. These include smart on/off constraints modeling, efficient bound propagation, domain-specific cut generation techniques and tight convexification procedures providing global optimality guarantees. PowerTools is still under active development, future versions will include Gas expansion and Unit Commitment models, decomposition methods, multi-threading, and support to additional solvers such as Cplex, Bonmin and others. Note that PowerTools is already used by another open-source project for Smart Grid simulation:  

2) Matpower

- Zimmerman, R., Murillo-Sanchez, C., Thomas, R.: Matpower: Steady-state operations, planning, and analysis tools for power systems research and education. Power Systems, IEEE Transactions on 26(1), 2011.

- D. K. Molzahn, J. T. Holzer, B. C. Lesieutre and C. L. DeMarco, "Implementation of a Large-Scale Optimal Power Flow Solver Based on Semidefinite Programming," in IEEE Transactions on Power Systems, vol. 28, no. 4, pp. 3987-3998, Nov. 2013.

- Kristiansen (2003), Utilising Matpower in Optimal Power Flow, Model Identification and Control, Vol 24, No 1, pp. 49-59.

- Lavei, Rantzer, Low (2011), Power Flow Optimisation using positive quadratic programming

3) OPFsolver

- Madani, R., Ashraphijuo, M., Lavaei, J. (2014)  SDP solver for Optimal Power Flow User'92s Manual.

- Madani, Ashraphijuo, Lavaei (2014), Promises of Conic Relaxations for Contingency constrained OPF, IEEE Trans. Power Systems.


G. P. Harrison and A. R. Wallace, "Optimal power flow evaluation of distribution network capacity for the connection of distributed generation," in IEEE Proceedings - Generation, Transmission and Distribution, vol. 152, no. 1, pp. 115-122, 10 Jan. 2005.

5) Plexos

- PLEXOS for power systems (2013)

- Foley, A. (2010), A strategic review of electricity system models, Energy, Vol 35 (12), No. 12, pp. 4522-4530.

General purpose solver applied to OPF


Alguacil, N., & Conejo, A. J. (2000). Multiperiod optimal power flow using Benders decomposition. Power Systems, IEEE Transactions on, 15(1), 196-201.

Jabr, R. (2008). Optimal power flow using an extended conic quadratic formulation. Power Systems, IEEE Transactions on, 23(3), 1000-1008.

Alsac, O., Bright, J., Prais, M., & Stott, B. (1990). Further developments in LP-based optimal power flow. Power Systems, IEEE Transactions on, 5(3), 697-711.


 R. A. Jabr, "Radial distribution load flow using conic programming," in IEEE Transactions on Power Systems, vol. 21, no. 3, pp. 1458-1459, Aug. 2006.

B. Kocuk, S. S. Dey and X. A. Sun, "Inexactness of SDP Relaxation and Valid Inequalities for Optimal Power Flow," in IEEE Transactions on Power Systems, vol. 31, no. 1, pp. 642-651, Jan. 2016.

Molzahn (2015), Solution of optimal power flow problems using moment relaxations augmented with objective function penalisation

3) Xpress



Dr Martin Mevissen, IBM

Dr Cedric Josz, LAAS CNRS, Toulouse

Dr Bissan Ghaddar, University of Waterloo

Dr Fabrizio Lacalandra, QuanTek