Making use of BDF-GMRES methods for solving short and long-term dynamics in power systems

Abstract This paper applies the BDF-GMRES methods for solving the Differential Algebraic Equations (DAEs) associated to the simulation of short and long-term dynamics in power systems. The investigations are concentrated on the construction of a fine ILU-GMRES preconditioner for solving efficiently not only the well-conditioned coefficient matrices but specially the ill-ones. It is shown that, if the image matrix (preconditioner origin) is firstly preprocessed (scaled, normalized and reordered), a high quality ILU preconditioner is achieved. Numerical experiments considering different test-systems and different operation conditions illustrate how tricky can be the simulation of power system dynamics if the Jacobian matrix (coefficient matrix) is ill-conditioned, normally associated to an adverse operation condition. It is shown that a traditional implicit integration method may fail in this case, whereas the combination BDF-GMRES presents an outstanding performance.

[1]  S. McFee,et al.  On the application of a pre-conditioned conjugate gradient algorithm to power network analysis , 1993 .

[2]  Ying Chen,et al.  Distributed Transient Stability Simulation of Power Systems Based on a Jacobian-Free Newton-GMRES Method , 2009 .

[3]  R. Bacher,et al.  Application of non-stationary iterative methods to an exact Newton-Raphson solution process for power flow equations , 1996 .

[4]  I. C. Decker,et al.  Parallel implementation of a power system dynamic simulation methodology using the conjugate gradient method , 1991 .

[5]  Peter W. Sauer,et al.  A preconditioned iterative solver for dynamic simulation of power systems , 1995, Proceedings of ISCAS'95 - International Symposium on Circuits and Systems.

[6]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[7]  Ahmad Sadeghi Yazdankhah,et al.  Impacts of distributed generations on power system transient and voltage stability , 2012 .

[8]  R. Prada,et al.  Critical investigation of preconditioned GMRES via incomplete LU factorization applied to power flow simulation , 2011 .

[9]  Richard Barrett,et al.  Templates for the Solution of Linear Systems: Building Blocks for Iterative Methods , 1994, Other Titles in Applied Mathematics.

[10]  Sudipta Ghosh,et al.  The localness of electromechanical oscillations in power systems , 2012 .

[11]  M. A. Pai,et al.  A new preconditioning technique for the GMRES algorithm in power flow and P-V curve calculations , 2003 .

[12]  Peter W. Sauer,et al.  Power System Dynamics and Stability , 1997 .

[13]  Linda R. Petzold,et al.  Using Krylov Methods in the Solution of Large-Scale Differential-Algebraic Systems , 1994, SIAM J. Sci. Comput..

[14]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[15]  J. Pasciak,et al.  Computer solution of large sparse positive definite systems , 1982 .

[16]  A. Bihain,et al.  The mixed Adams-BDF variable step size algorithm to simulate transient and long term phenomena in power systems , 1994 .

[17]  Iain S. Duff,et al.  On Algorithms For Permuting Large Entries to the Diagonal of a Sparse Matrix , 2000, SIAM J. Matrix Anal. Appl..

[18]  Iain S. Duff,et al.  The Design and Use of Algorithms for Permuting Large Entries to the Diagonal of Sparse Matrices , 1999, SIAM J. Matrix Anal. Appl..

[19]  Ying Chen,et al.  A Jacobian-free Newton-GMRES(m) method with adaptive preconditioner and its application for power flow calculations , 2006, IEEE Transactions on Power Systems.

[20]  Siddhartha Kumar Khaitan,et al.  A Class of New Preconditioners for Linear Solvers Used in Power System Time-Domain Simulation , 2010, IEEE Transactions on Power Systems.

[21]  Michael A. Saunders,et al.  SNOPT: An SQP Algorithm for Large-Scale Constrained Optimization , 2005, SIAM Rev..

[22]  A. Bose,et al.  A highly parallel method for transient stability analysis , 1989, Conference Papers Power Industry Computer Application Conference.

[23]  Elizabeth Cuthill,et al.  Several Strategies for Reducing the Bandwidth of Matrices , 1972 .

[24]  H.R. Fankhauser,et al.  Advanced simulation techniques for the analysis of power system dynamics , 1990, IEEE Computer Applications in Power.

[25]  M. A. Pai,et al.  Iterative solver techniques in large scale power system computation , 1997, Proceedings of the 36th IEEE Conference on Decision and Control.

[26]  C. W. Gear,et al.  Numerical initial value problem~ in ordinary differential eqttations , 1971 .

[27]  Yousef Saad,et al.  ILUT: A dual threshold incomplete LU factorization , 1994, Numer. Linear Algebra Appl..

[28]  Jawad Talaq,et al.  Optimal power system stabilizers for multi machine systems , 2012 .

[29]  J. Cash.  Modified extended backward differentiation formulae for the numerical solution of stiff initial value problems in ODEs and DAEs , 2000 .

[30]  Adam Semlyen.  Fundamental concepts of a Krylov subspace power flow methodology , 1996 .

[31]  C. Oldenburg,et al.  Numerical Simulation of Ferrofluid Flow for Subsurface Environmental Engineering Applications , 2000 .

[32]  Linda R. Petzold,et al.  Numerical solution of initial-value problems in differential-algebraic equations , 1996, Classics in applied mathematics.

[33]  H. Mori,et al.  An ILU(p)-Preconditoner Bi-CGStab Method for Power Flow Calculation , 2007, 2007 IEEE Lausanne Power Tech.

[34]  I. Duff,et al.  Direct Methods for Sparse Matrices , 1987 .

[35]  Carl D. Meyer,et al.  Matrix Analysis and Applied Linear Algebra , 2000 .

[36]  M. Benzi.  Preconditioning techniques for large linear systems: a survey , 2002 .

[37]  Solving the nonlinear power flow equations with a Newton process and GMRES , 1996, 1996 IEEE International Symposium on Circuits and Systems. Circuits and Systems Connecting the World. ISCAS 96.