The Distribution of the Number of Real Solutions to the Power Flow Equations

In this paper we study the distributions of the number of real solutions to the power flow equations over varying electrical parameters. We introduce a new monodromy and parameter homotopy continuation method for quickly finding all solutions to the power flow equations. We apply this method to find distributions of the number of real solutions to the power flow equations and compare these distributions to those of random polynomials. It is observed that while the power flow equations tend to admit many fewer real-valued solutions than a bound on the total number of complex solutions, for low levels of load they tend to admit many more than a corresponding random polynomial. We show that for cycle graphs the number of real solutions can achieve the maximum bound for specific parameter values and for complete graphs with four or more vertices there are susceptance values that give infinitely many real solutions.

[1]  Yong Min,et al.  A Continuation-Based Method to Compute the Relevant Unstable Equilibrium Points for Power System Transient Stability Analysis , 2009 .

[2]  J. Kiefer,et al.  Asymptotic Minimax Character of the Sample Distribution Function and of the Classical Multinomial Estimator , 1956 .

[3]  Abraham Martín del Campo,et al.  Critical points via monodromy and local methods , 2015, 1503.01662.

[4]  Alisha Zachariah,et al.  Distributions of the Number of Solutions to the Network Power Flow Equations , 2018, 2018 IEEE International Symposium on Circuits and Systems (ISCAS).

[5]  Julia Lindberg,et al.  The Geometry of Real Solutions to the Power Flow Equations , 2018, 2018 56th Annual Allerton Conference on Communication, Control, and Computing (Allerton).

[6]  Chen-Sung Chang,et al.  Toward a CPFLOW-based algorithm to compute all the type-1 load-flow solutions in electric power systems , 2005, IEEE Transactions on Circuits and Systems I: Regular Papers.

[7]  Dan Wu,et al.  An efficient method to locate all the load flow solutions - revisited , 2015, 2015 53rd Annual Allerton Conference on Communication, Control, and Computing (Allerton).

[8]  Ronald L. Graham,et al.  Concrete mathematics - a foundation for computer science , 1991 .

[9]  James S. Thorp,et al.  An efficient algorithm to locate all the load flow solutions , 1993 .

[10]  Andrew J. Sommese,et al.  The numerical solution of systems of polynomials - arising in engineering and science , 2005 .

[11]  Konstantin Turitsyn,et al.  Numerical polynomial homotopy continuation method to locate all the power flow solutions , 2014, 1408.2732.

[12]  Dhagash Mehta,et al.  On the Network Topology Dependent Solution Count of the Algebraic Load Flow Equations , 2015, IEEE Transactions on Power Systems.

[13]  L. Ni,et al.  Parallel processing for the load flow of power systems: the approach and applications , 1989, Proceedings of the 28th IEEE Conference on Decision and Control,.

[14]  Nigel Boston,et al.  Exploiting Symmetry in the Power Flow Equations Using Monodromy , 2020, ArXiv.

[15]  Timothy Duff,et al.  Solving polynomial systems via homotopy continuation and monodromy , 2016, ArXiv.

[16]  D. N. Bernshtein.  The number of roots of a system of equations , 1975 .

[17]  Dinesh Manocha,et al.  SOLVING SYSTEMS OF POLYNOMIAL EQUATIONS , 2002 .

[18]  Enrique Mallada,et al.  Distributed Synchronization of Heterogeneous Oscillators on Networks With Arbitrary Topology , 2014, IEEE Transactions on Control of Network Systems.

[19]  Bernard Lesieutre,et al.  On the Distribution of Real-Valued Solutions to the Power Flow Equations , 2019, 2019 57th Annual Allerton Conference on Communication, Control, and Computing (Allerton).

[20]  K. Farahmand.  On the Average Number of Real Roots of a Random Algebraic Equation , 1986 .

[21]  Fabrice Rouillier,et al.  Solving Zero-Dimensional Systems Through the Rational Univariate Representation , 1999, Applicable Algebra in Engineering, Communication and Computing.

[22]  B. C. Lesieutre,et al.  Counterexample to a Continuation-Based Algorithm for Finding All Power Flow Solutions , 2013, IEEE Transactions on Power Systems.

[23]  Lei Chen,et al.  A Continuation-Based Method to Compute the Relevant Unstable Equilibrium Points for Power System Transient Stability Analysis , 2009, IEEE Transactions on Power Systems.

[24]  Thomas J. Overbye,et al.  A new method for finding low-voltage power flow solutions , 2000, 2000 Power Engineering Society Summer Meeting (Cat. No.00CH37134).

[25]  A. G. Kouchnirenko.  Polyèdres de Newton et nombres de Milnor , 1976 .

[26]  Bernd Sturmfels,et al.  A polyhedral method for solving sparse polynomial systems , 1995 .

[27]  B. F. Caviness,et al.  Quantifier Elimination and Cylindrical Algebraic Decomposition , 2004, Texts and Monographs in Symbolic Computation.

[28]  T. Y. Li.  Numerical solution of multivariate polynomial systems by homotopy continuation methods , 2008 .

[29]  J. Lambert.  Numerical Methods for Ordinary Differential Equations , 1991 .

[30]  Jonathan D. Hauenstein,et al.  Multiprojective witness sets and a trace test , 2015, Advances in Geometry.

[31]  Paul Breiding,et al.  HomotopyContinuation.jl: A Package for Homotopy Continuation in Julia , 2017, ICMS.

[32]  Dhagash Mehta,et al.  Counting Equilibria of the Kuramoto Model Using Birationally Invariant Intersection Index , 2017, SIAM J. Appl. Algebra Geom..

[33]  Jean‐Claude Sabonnadière,et al.  Stability Analysis of Power Systems , 2013 .

[34]  A. Khovanskii.  Newton polyhedra and the genus of complete intersections , 1978 .

[35]  J. Baillieul,et al.  Geometric critical point analysis of lossless power system models , 1982 .

[36]  J. Rodriguez,et al.  Solving Parameterized Polynomial Systems with Decomposable Projections , 2016, 1612.08807.