Search

In this paper, we consider the optimal power flow (OPF) problem which consists in determining the power production at each bus of an electric network by minimizing the production cost. Our contribution is an exact solution algorithm for the OPF problem. It consits in a spatial branch-and-bound algorithm based on a compact quadratically constrained convex relaxation. This compact relaxation is computed by solving the rank relaxation once at the beginning of the algorithm. The key point of this approach is that the lower bound at the root node of the branch-and-bound tree is equal to the rank relaxation value, but is obtained by solving a quadratic convex problem which is much faster than solving a SDP. To construct this compact relaxation, we add only O(n) variables that model the squares of the initial variables, where $n$ is the number of buses in the power system. The relations between the initial and auxiliary variables are therefore non-convex. By relaxing them in our relaxation, we have only O(n) equalities to force by the branch-and-bound algorithm to prove global optimality. Our first experiments show that our new algorithm Compact OPF (COPF) performs better than the methods of the literature we compare it with

We consider the exact solution of problem $(QP)$ that consists in minimizing a quadratic function subject to quadratic constraints. Starting from the classical convex relaxation that uses the McCormick's envelopes, we introduce 12 inequalities that are derived from the ranges of the variables of $(QP)$. We prove that these general Triangle inequalities cut feasible solutions of the McCormick's envelopes. We then show how we can compute a convex relaxation $(P^*)$ which optimal value equals to the "Shor's plus RLT plus Triangle" semi-definite relaxation of $(QP)$ that includes the new inequalities. We also propose a heuristic for solving this huge semi-definite program that serves as a separation algorithm. We then solve $(QP)$ to global optimality with a branch-and-bound based on $(P^*)$. Moreover, as the new inequalities involved the lower and upper bounds on the original variables of $(QP)$, their use in a branch-and-bound framework accelerates the whole process. We show on the unitbox instances that our method outperforms the compared solvers., preprint