Solving a significant number of differential-algebraic equations is essential for conducting transient stability simulations on a complex and interconnected electrical power system. Because of the increasing size and intricacy of modern grids, traditional sequential simulation methods are becoming more computationally challenging and time-consuming to use for dynamic simulations. This research aims to examine parallelism techniques from an innovative algorithm viewpoint to the hardware level by utilizing multiple graphics processing units as the main computational tool designed to simulate the dynamic behaviour of power systems on a large scale. Our approach features a novel hybrid algorithm to improve the dynamic simulations of grids with significant renewable energy integration. The algorithm employs a Schwarz-based approach in its initial decomposition phase to isolate synchronous machines and renewable sources from the grid. This separation enhances the algorithm's ability to manage these components effectively. Subsequently, each subsystem undergoes further division into subdomains in the second stage, with separate computations carried out using the Schur-complement technique. The simulations were conducted on various test systems, with a maximum of 25,000 buses, 8,000 synchronous generators, and 256 PV farms all modelled in detail. We introduce a GPU-oriented preprocessing and vectorization parallelization method for the implementation of our two-stage algorithm. This approach provides the flexibility to tailor parallelization techniques, including nested for loops and kernel pragmas, to craft a hardware solution that enhances the performance of our dynamic domain decomposition algorithm. Simulation outcomes have verified that this approach can yield a notable 7.8x acceleration in speed when executed on an NVIDIA GeForce RTX 2070 SUPER GPU.
Read full abstract