Index-aware model order reduction for differential-algebraic equations

We introduce a model order reduction (MOR) procedure for differential-algebraic equations, which is based on the intrinsic differential equation contained in the starting system and on the remaining algebraic constraints. The decoupling procedure in differential and algebraic part is based on the projector and matrix chain which leads to the definition of tractability index. The differential part can be reduced by using any MOR method, we use Krylov-based projection methods to illustrate our approach. The reduction on the differential part induces a reduction on the algebraic part. In this paper, we present the method for index-1 differential-algebraic equations. We implement numerically this procedure and show numerical evidence of its validity.

[1]  Lawrence T. Pileggi,et al.  Asymptotic waveform evaluation for timing analysis , 1990, IEEE Trans. Comput. Aided Des. Integr. Circuits Syst..

[2]  Nonlinear functional analysis and its applications, part I: Fixed-point theorems , 1991 .

[3]  C. W. Gear,et al.  Simultaneous Numerical Solution of Differential-Algebraic Equations , 1971 .

[4]  S. Selberherr Analysis and simulation of semiconductor devices , 1984 .

[5]  Roland W. Freund,et al.  SPRIM: structure-preserving reduced-order interconnect macromodeling , 2004, IEEE/ACM International Conference on Computer Aided Design, 2004. ICCAD-2004..

[6]  U. Ascher,et al.  Numerical Methods for Differential-Algebraic Equations , 1998 .

[7]  Michael Günther,et al.  A PDAE Model for Interconnected Linear RLC Networks , 2001 .

[8]  Joseph W. Jerome Analysis of Charge Transport: A Mathematical Study of Semiconductor Devices , 1995 .

[9]  Roswitha März,et al.  Basic Properties of Some Differential-Algebraic Equations , 1989 .

[10]  Peter A. Markowich,et al.  The Stationary Semiconductor Device Equations. , 1987 .

[11]  E. Hairer,et al.  Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems , 2010 .

[12]  R. März,et al.  Projector based treatment of linear constant coefficient DAEs , 2011 .

[13]  Athanasios C. Antoulas,et al.  Approximation of Large-Scale Dynamical Systems , 2005, Advances in Design and Control.

[14]  H. V. D. Vorst,et al.  Model Order Reduction: Theory, Research Aspects and Applications , 2008 .

[15]  L. Sirovich Turbulence and the dynamics of coherent structures. I. Coherent structures , 1987 .

[16]  D. Stewart,et al.  Error estimates for Krylov subspace approximations of matrix exponentials , 1996 .

[17]  Roswitha März,et al.  Canonical projectors for linear differential algebraic equations , 1996 .

[18]  H. Tal-Ezer Spectral methods in time for parabolic problems , 1989 .

[19]  Joost Rommes,et al.  Computing Rightmost Eigenvalues for Small-Signal Stability Assessment of Large-Scale Power Systems , 2010, IEEE Transactions on Power Systems.

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

[21]  William J. McCalla Fundamentals of Computer-Aided Circuit Simulation , 1987 .

[22]  Albert E. Ruehli,et al.  The modified nodal approach to network analysis , 1975 .

[23]  Ngai Wong,et al.  An Efficient Projector-Based Passivity Test for Descriptor Systems , 2010, IEEE Transactions on Computer-Aided Design of Integrated Circuits and Systems.

[24]  B. Moore Principal component analysis in linear systems: Controllability, observability, and model reduction , 1981 .

[25]  Karen Willcox,et al.  Model Reduction of an Actively Controlled Supersonic Diffuser , 2005 .

[26]  Tatjana Stykel,et al.  Gramian-Based Model Reduction for Descriptor Systems , 2004, Math. Control. Signals Syst..

[27]  Roland W. Freund,et al.  Efficient linear circuit analysis by Pade´ approximation via the Lanczos process , 1994, EURO-DAC '94.

[28]  K. Glover All optimal Hankel-norm approximations of linear multivariable systems and their L, ∞ -error bounds† , 1984 .

[29]  N. Martins,et al.  Reduced-order transfer matrices from RLC network descriptor models of electric power grids , 2012, 2012 IEEE Power and Energy Society General Meeting.

[30]  N. Martins,et al.  Gramian-Based Reduction Method Applied to Large Sparse Power System Descriptor Models , 2008, IEEE Transactions on Power Systems.

[31]  L. Sirovich Turbulence and the dynamics of coherent structures. II. Symmetries and transformations , 1987 .

[32]  Joost Rommes,et al.  Efficient Methods for Large Resistor Networks , 2010, IEEE Transactions on Computer-Aided Design of Integrated Circuits and Systems.