GPU-based power flow analysis with Chebyshev preconditioner and conjugate gradient method

a b s t r a c t Traditionally, linear equations in power system applications are solved by direct methods based on LU decomposition. With the development of advanced power system controls, the industrial and research community is more interested in simulating larger, interconnected power grids. Iterative methods such as the conjugate gradient method have been applied to power system applications in the literature for its parallelism potential with larger systems. Preconditioner, used for preconditioning the linear system for a better convergence rate in iterative computations, is an indispensable part of iterative solving process. This work implemented a polynomial preconditioner Chebyshev preconditioner with graphic processing unit (GPU), and integrated a GPU-based conjugate gradient solver. Results show that GPU-based Chebyshev preconditioner can reach around 46× speedup for the largest test system, and conjugate gradient can gain more than 4× speedup. This demonstrates great potentials for GPU application in power system

[1]  Jonas Koko,et al.  Parallel preconditioned conjugate gradient algorithm on GPU , 2012, J. Comput. Appl. Math..

[2]  Janusz Bialek,et al.  Approximate model of European interconnected system as a benchmark system to study effects of cross-border trades , 2005 .

[3]  Robert C. Green,et al.  High performance computing for electric power systems: Applications and trends , 2011, 2011 IEEE Power and Energy Society General Meeting.

[4]  L. van der Sluis,et al.  Scalable Newton-Krylov Solver for Very Large Power Flow Problems , 2010, IEEE Transactions on Power Systems.

[5]  Amirhassan Asgari Kamiabad.  Implementing a Preconditioned Iterative Linear Solver Using Massively Parallel Graphics Processing Units , 2011 .

[6]  Tom R. Halfhill.  NVIDIA's Next-Generation CUDA Compute and Graphics Architecture, Code-Named Fermi, Adds Muscle for Parallel Processing , 2009 .

[7]  M. A. Pai,et al.  Iterative solver techniques in large scale power system computation , 1997, Proceedings of the 36th IEEE Conference on Decision and Control.

[8]  Norberto Garcia.  Parallel power flow solutions using a biconjugate gradient algorithm and a Newton method: A GPU-based approach , 2010, IEEE PES General Meeting.

[9]  Richard F. Barrett,et al.  Matrix Market: a web resource for test matrix collections , 1996, Quality of Numerical Software.

[10]  D. Niebur,et al.  DC Power Flow Based Contingency Analysis Using Graphics Processing Units , 2007, 2007 IEEE Lausanne Power Tech.

[11]  P. W. Sauer,et al.  Conjugate Gradient Approach to Parallel Processing in Dynamic Simulation of Power Systems , 1992, 1992 American Control Conference.

[12]  Fangxing Li,et al.  Exploration of multifrontal method with GPU in power flow computation , 2013, 2013 IEEE Power & Energy Society General Meeting.

[13]  R D Zimmerman,et al.  MATPOWER: Steady-State Operations, Planning, and Analysis Tools for Power Systems Research and Education , 2011, IEEE Transactions on Power Systems.

[14]  Zhao Li,et al.  On limitations of traditional multi-core and potential of many-core processing architectures for sparse linear solvers used in large-scale power system applications , 2011, 2011 IEEE Power and Energy Society General Meeting.

[15]  F. de Leon,et al.  Iterative solvers in the Newton power flow problem: preconditioners, inexact solutions and partial Jacobian updates , 2002 .

[16]  Lei Zhang,et al.  Efficient CUDA Polynomial Preconditioned Conjugate Gradient Solver for Finite Element Computation of Elasticity Problems , 2013 .

[17]  A. Semlyen,et al.  A New Preconditioned Conjugate Gradient Power Flow , 2002, IEEE Power Engineering Review.