Stability and Convergency Exploration of Matrix Exponential Integration on Power Delivery Network Transient Simulation

We propose a stability preserved Arnoldi algorithm for matrix exponential in the time domain simulation of large-scale power delivery networks (PDNs), which are formulated as semi-explicit differential-algebraic equations (DAEs). The matrix exponential and vector products (MEVPs) compose the solution of DAEs in multistep integration methods and can be efficiently approximated with the rational Krylov subspace. To produce stable simulation results for the ill-conditioned system from semi-explicit DAEs, the revised Arnoldi algorithm introduces a new structured orthogonalization process to construct the Krylov subspace. We demonstrate the performance of the new algorithm with theoretical proof and experiments. In the computation of MEVPs, we utilize the exponential related $\varphi $ functions to improve the numerical accuracy. We further explore the optimal ratio to confine the spectrum in the rational Krylov subspace. Finally, the transient framework is tested on a group of system-level PDNs, showing that matrix exponential-based algorithms could achieve high efficiency and accuracy.

[1]  Karl Meerbergen,et al.  Implicitly restarted Arnoldi with purification for the shift-invert transformation , 1997, Math. Comput..

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

[3]  Xiang Zhang,et al.  Simulation Algorithms With Exponential Integration for Time-Domain Analysis of Large-Scale Power Delivery Networks , 2015, IEEE Transactions on Computer-Aided Design of Integrated Circuits and Systems.

[4]  Marlis Hochbruck,et al.  Preconditioning Lanczos Approximations to the Matrix Exponential , 2005, SIAM J. Sci. Comput..

[5]  Timo Reis,et al.  Surveys in Differential-Algebraic Equations III , 2015 .

[6]  W. Marsden.  I and J , 2012 .

[7]  Mattan Kamon,et al.  A coordinate-transformed Arnoldi algorithm for generating guaranteed stable reduced-order models of RLC circuits , 1996, Proceedings of International Conference on Computer Aided Design.

[8]  Gerhard Wanner,et al.  Dahlquist’s classical papers on stability theory , 2006 .

[9]  Chung-Kuan Cheng,et al.  Transient Circuit Simulation for Differential Algebraic Systems using Matrix Exponential , 2018, 2018 IEEE/ACM International Conference on Computer-Aided Design (ICCAD).

[10]  Farid N. Najm,et al.  Circuit Simulation , 2010 .

[11]  M.L. Liou,et al.  Computer-aided analysis of electronic circuits: Algorithms and computational techniques , 1977, Proceedings of the IEEE.

[12]  Axel Ruhe,et al.  The spectral transformation Lánczos method for the numerical solution of large sparse generalized symmetric eigenvalue problems , 1980 .

[13]  Cleve B. Moler,et al.  Nineteen Dubious Ways to Compute the Exponential of a Matrix, Twenty-Five Years Later , 1978, SIAM Rev..

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

[15]  Chung-Kuan Cheng,et al.  Exploring the exponential integrators with Krylov subspace algorithms for nonlinear circuit simulation , 2017, 2017 IEEE/ACM International Conference on Computer-Aided Design (ICCAD).

[16]  Satoru Iwata,et al.  Index characterization of differential-algebraic equations in hybrid analysis for circuit simulation , 2010, Int. J. Circuit Theory Appl..

[17]  B. Parlett,et al.  How to implement the spectral transformation , 1987 .

[18]  Yiyu Shi,et al.  Modeling and design for beyond-the-die power integrity , 2010, 2010 IEEE/ACM International Conference on Computer-Aided Design (ICCAD).

[19]  Lawrence T. Pileggi,et al.  PRIMA: passive reduced-order interconnect macromodeling algorithm , 1998, 1997 Proceedings of IEEE International Conference on Computer Aided Design (ICCAD).

[20]  R. Freund.  Krylov-subspace methods for reduced-order modeling in circuit simulation , 2000 .

[21]  Sani R. Nassif,et al.  Power grid analysis benchmarks , 2008, 2008 Asia and South Pacific Design Automation Conference.

[22]  Joost Rommes,et al.  Exploiting structure in large-scale electrical circuit and power system problems , 2009 .

[23]  R. Winkler.  Stochastic differential algebraic equations of index 1 and applications in circuit simulation , 2003 .

[24]  P. Rentrop,et al.  The Drazin inverse in multibody system dynamics , 1993 .

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

[26]  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.

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