A decoupled time-domain Simulation method via Invariant subspace partition for power system analysis

A decoupled method is proposed to deal with time-domain simulation for power system dynamic analysis. Traditionally, there are two main categories of numerical integration methods: explicit methods and implicit methods. The implicit methods are numerically stable but require more computational time to solve the nonlinear equations, while explicit methods are relatively efficient but may cause a numerical stability problem. This paper proposes a new hybrid method to take advantage of both explicit and implicit methods based on the invariant subspace partition. The original power system equations are decoupled into two parts that correspond to the stiff and nonstiff subspaces. For the stiff invariant subspace, the implicit method is applied to achieve numerical stability, and the explicit method is employed to handle nonstiff invariant subspace for the computational efficiency. As a result, the new hybrid method is both numerically stable and efficient. The approach is demonstrated through New England 39-bus and IEEE 118-bus systems.

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

[2]  E. Hofer,et al.  A Partially Implicit Method for Large Stiff Systems of ODEs with Only Few Equations Introducing Small Time-Constants , 1976 .

[3]  Mohamed S. Kamel,et al.  Automatic Partitioning of Stiff Systems and Exploiting the Resulting Structure , 1979, TOMS.

[4]  Gautam M. Shroff,et al.  Stabilization of unstable procedures: the recursive projection method , 1993 .

[5]  J. Lambert Numerical Methods for Ordinary Differential Systems: The Initial Value Problem , 1991 .

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

[7]  Yuan Zhou,et al.  Local parameterization approach for unified time domain simulation of power system dynamics , 2002, 2002 IEEE Power Engineering Society Winter Meeting. Conference Proceedings (Cat. No.02CH37309).

[8]  A projection method for the numerical solution of linear systems in separable stiff differential equations , 1989 .

[9]  Uri M. Ascher,et al.  Computer methods for ordinary differential equations and differential-algebraic equations , 1998 .

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

[11]  J. Butcher Numerical methods for ordinary differential equations in the 20th century , 2000 .

[12]  G. Akrivis A First Course In The Numerical Analysis Of Differential Equations [Book News & Reviews] , 1998, IEEE Computational Science and Engineering.

[13]  Peter Deuflhard,et al.  Scientific Computing with Ordinary Differential Equations , 2002 .

[14]  Vladimír Janovský,et al.  Continuation of Invariant Subspaces via the Recursive Projection Method , 2003 .

[15]  C. W. Gear,et al.  Iterative Solution of Linear Equations in ODE Codes , 1983 .

[16]  Y. Saad Krylov subspace methods for solving large unsymmetric linear systems , 1981 .

[17]  M. Stubbe,et al.  STAG-A New Unified Software Program for the Study of the Dynamic Behaviour of Electrical Power Systems , 1989, IEEE Power Engineering Review.

[18]  Venkataramana Ajjarapu,et al.  Identification of voltage collapse through direct equilibrium tracing , 2000 .

[19]  F. Iavernaro,et al.  Boundary values methods for time-domain simulation of power system dynamic behavior , 1998 .

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

[21]  David S. Watkins,et al.  The Numerical Solution of Separably Stiff Systems by Precise Partitioning , 1983, TOMS.

[22]  Åke Björck,et al.  A blockQR algorithm for partitioning stiff differential systems , 1983 .

[23]  Chih-Wen Liu,et al.  New methods for computing power system dynamic response for real-time transient stability prediction , 2000 .

[24]  H. Okubo,et al.  Multiple time-scale power system dynamic simulation , 1993 .