[Paper Review] Investigating the Maximum Number of Real Solutions to the Power Flow Equations: Analysis of Lossless Four-Bus Systems

Daniel K. Molzahn, Matthew Niemerg+2 more|arXiv (Cornell University)|Mar 18, 2016

Polynomial and algebraic computation37 references5 citations

TL;DR

This paper investigates the maximum number of real solutions to the power flow equations in lossless four-bus systems composed of PV buses with unity voltage magnitudes. Using numerical algebraic geometry and Galois group analysis, it shows that the system can have up to 16 real solutions—strictly fewer than the 20 complex solutions—providing strong evidence that a gap exists between real and complex solution counts in larger systems.

ABSTRACT

The power flow equations model the steady-state relationship between the power injections and voltage phasors in an electric power system. By separating the real and imaginary components of the voltage phasors, the power flow equations can be formulated as a system of quadratic polynomials. Only the real solutions to these polynomial equations are physically meaningful. This paper focuses on the maximum number of real solutions to the power flow equations. An upper bound on the number of real power flow solutions commonly used in the literature is the maximum number of complex solutions. There exist two- and three-bus systems for which all complex solutions are real. It is an open question whether this is also the case for larger systems. This paper investigates four-bus systems using techniques from numerical algebraic geometry and conjectures a negative answer to this question. In particular, this paper studies lossless, four-bus systems composed of PV buses connected by lines with arbitrary susceptances. Computing the Galois group, which is degenerate, enables conversion of the problem of counting the number of real solutions to the power flow equations into counting the number of positive roots of a univariate sextic polynomial. From this analysis, it is conjectured that the system has at most 16 real solutions, which is strictly less than the maximum number of complex solutions, namely 20. We also provide explicit parameter values where this system has 16 real solutions so that the conjectured upper bound is achievable.

Motivation & Objective

• To determine the maximum number of real solutions to the power flow equations in lossless four-bus systems.
• To investigate whether the number of real solutions can reach the upper bound of complex solutions (20) in four-bus systems.
• To analyze the structure of power flow solutions using algebraic geometry techniques, particularly Galois group theory.
• To provide explicit parameter configurations where the maximum number of real solutions is achieved.
• To conjecture that 16 is the tightest upper bound on real solutions for such systems, based on numerical and theoretical evidence.

Proposed method

• Formulates the power flow equations as a system of quadratic polynomials by separating real and imaginary parts of voltage phasors.
• Applies numerical polynomial homotopy continuation (NPHC) to compute all complex solutions and track solution paths.
• Uses Galois group theory to analyze the symmetry of the system and reduce the problem to counting positive roots of a univariate sextic polynomial.
• Employs the implicit function theorem to extend results from special cases (e.g., zero line susceptance) to nearby parameter regimes.
• Constructs explicit parameter sets (e.g., b₁₂ = 0) where all 16 complex solutions are real, confirming the bound is achievable.
• Analyzes solutions at infinity via homogenization to show that four solutions diverge and remain nonreal when a line susceptance approaches zero.

Experimental results

Research questions

• RQ1 Can all 20 complex solutions of the four-bus power flow equations be real in a lossless system with PV buses?
• RQ2 What is the maximum number of real solutions possible in a lossless four-bus system with unity voltage magnitudes and zero power injections?
• RQ3 Does the number of real solutions remain bounded when line susceptances are varied, particularly when one is near zero?
• RQ4 Can algebraic geometry techniques such as Galois group analysis be used to reduce the solution counting problem to univariate polynomial root counting?
• RQ5 Is there a strict gap between the number of complex solutions (20) and the maximum number of real solutions in four-bus systems?

Key findings

• The system achieves 16 real solutions when one line susceptance is zero, demonstrating that the upper bound of 16 is achievable.
• All 16 complex solutions can be real when b₁₂ = 0, as confirmed by the positive roots of a quartic polynomial derived from the sextic univariate reduction.
• When a line susceptance is zero, the number of complex solutions drops from 20 to 16, with the remaining four solutions diverging to infinity and remaining nonreal.
• The four solutions at infinity are nonreal for all values of the parameter near zero, as they satisfy V_di² + V_qi² = 0, which only has the trivial real solution.
• The Galois group analysis confirms that the system’s solution structure allows the number of real solutions to be determined by counting positive roots of a univariate sextic polynomial.
• Based on numerical experiments and theoretical analysis, the paper conjectures that 16 is the maximum number of real solutions for general lossless four-bus PV systems, even when no parameter is near zero.