The parallel implementation of the waveform relaxation method for transient stability simulations

The WR (waveform relaxation) algorithm is extended to a structure-preserving power system model in which the loads are retained. This results in a system of differential/algebraic equations (DAEs). Power system exhibit several unique dynamic properties which may be exploited in an advantageous manner by the WR algorithm. These physical properties include the coherency properties of the power system which lead to the partitions for the textured model approach, the near diagonal dominance which leads to longer windows for uniform convergence, and the localized response from which the multirate capabilities of the WR method can be used. These characteristics enable power systems to obtain more favorable results than were obtained in VLSI simulations. The authors present several theoretical results as well as computational results on parallel implementation.