Matrix exponential based electromagnetic transients simulation algorithm with Krylov subspace approximation and accurate dense output

The ever-increasing scale of system interconnection and diverse dynamic time constants of power equipment pose challenges to power system electromagnetic transients simulation algorithm. In this paper, an advanced numerical integration method based on matrix exponential operator is proposed. Krylov subspace approximation is utilized to effectively reduce the problem dimension, and an accurate dense output formula is established to enable the method adopting larger simulation step size and revealing higher frequency dynamic details simultaneously. Accuracy, numerical stability and computation efficiency is demonstrated through numerical case studies.

[1]  Awad H. Al-Mohy,et al.  A New Scaling and Squaring Algorithm for the Matrix Exponential , 2009, SIAM J. Matrix Anal. Appl..

[2]  Awad H. Al-Mohy,et al.  Computing the Action of the Matrix Exponential, with an Application to Exponential Integrators , 2011, SIAM J. Sci. Comput..

[3]  C. Loan,et al.  Nineteen Dubious Ways to Compute the Exponential of a Matrix , 1978 .

[4]  Marlis Hochbruck,et al.  Exponential Rosenbrock-Type Methods , 2008, SIAM J. Numer. Anal..

[5]  N. Higham.  The Scaling and Squaring Method for the Matrix Exponential Revisited , 2005, SIAM J. Matrix Anal. Appl..

[6]  F. Aluffi-Pentini,et al.  A Novel Algorithm for the Numerical Integration of Systems of Ordinary Differential Equations Arising in Chemical Problems , 2003 .

[7]  Oleg Wasynczuk,et al.  Automated state model generation algorithm for power circuits and systems , 1996 .

[8]  Yoshihiro Yamane,et al.  Explicit Time Integration Scheme Using Krylov Subspace Method for Reactor Kinetics Equation , 2011 .

[9]  C. Lubich,et al.  On Krylov Subspace Approximations to the Matrix Exponential Operator , 1997 .

[10]  Quan Chen,et al.  Time-Domain Analysis of Large-Scale Circuits by Matrix Exponential Method With Adaptive Control , 2012, IEEE Transactions on Computer-Aided Design of Integrated Circuits and Systems.

[11]  Jitse Niesen,et al.  Algorithm 919: A Krylov Subspace Algorithm for Evaluating the ϕ-Functions Appearing in Exponential Integrators , 2009, TOMS.

[12]  Y. Saad.  Analysis of some Krylov subspace approximations to the matrix exponential operator , 1992 .

[13]  Roger B. Sidje,et al.  Expokit: a software package for computing matrix exponentials , 1998, TOMS.

[14]  A. Ostermann,et al.  Implementation of exponential Rosenbrock-type integrators , 2009 .

[15]  M. Hochbruck,et al.  Exponential integrators , 2010, Acta Numerica.

[16]  W. H. Kersting,et al.  Radial distribution test feeders , 1991, 2001 IEEE Power Engineering Society Winter Meeting. Conference Proceedings (Cat. No.01CH37194).