Numerical Simulation of Transient Stability of Power Systems by an Adaptive Runge‐Kutta Method

A new algorithm for the transient stability simulation of power systems, which are represented by a set of nonlinear different ml and algebraic equations, is proposed. This algorithm is based on the adaptive continuous Runge‐Kutta method of order three with a step changing strategy based on estimation of the local discretization error by an embedded two‐step Runge‐Kutta method of order four. A shift time threshold and a dead band are used to improve both the accuracy and efficiency of the algorithm. This method was tested on the 118 bus power system with 19 generators and was found faster than the usually preferred method based on a simultaneous approach. A bigger speedup is expected for larger systems.