A Block-Parallel Newton Method via Overlapping Epsilon Decompositions

The purpose of this paper is to present a block-parallel Newton method for solving large nonlinear systems. A graph-theoretic decomposition algorithm is first used to partition the Jacobian into weakly coupled, possibly overlapping blocks. It is then shown that it suffices to invert only the diagonal blocks to carry out the Newton iterates. A rigorous justification of this practice is provided by using a convergence result of Kantorovich in the expanded space of the iterates, where overlapping blocks appear as disjoint. The individual blocks, or a group of blocks, can be inverted by a dedicated processor, making the new block-diagonal Newton method ideally suited for parallel processing. Applications to the power flow problems are presented, and parallelization issues are discussed briefly.

[1]  H. H. Happ,et al.  Multicomputer Configurations and Diakoptics: Dispatch of Real Power in Power Pools , 1969 .

[2]  Fred C. Schweppe,et al.  Distance Measures and Coherency Recognition for Transient Stability Equivalents , 1973 .

[3]  H. H. Happ,et al.  Diakoptics—The solution of system problems by tearing , 1974 .

[4]  B. Stott,et al.  Review of load-flow calculation methods , 1974 .

[5]  J. Medanic,et al.  Solution of load-flow problems in power systems by ϵ-coupling method , 1975 .

[6]  R. Kasturi,et al.  Piecewise Newton-Raphson load flow ߞAn exact method using ordered elimination , 1976, IEEE Transactions on Power Apparatus and Systems.

[7]  Ikeda,et al.  Overlapping decompositions, expansions, and contractions of dynamic systems , 1979 .

[8]  M. Pai,et al.  Electromechanical distance measure for decomposition of power systems , 1984 .

[9]  J. Zaborszky,et al.  A Textured Model for Computationally Efficient Reactive Power Control and Management , 1985, IEEE Transactions on Power Apparatus and Systems.

[10]  Dragoslav D. Siljak,et al.  Nested ε-decompositions and clustering of complex systems , 1986, Autom..

[11]  M. Pai,et al.  Two time scale decomposition and stability analysis of power systems , 1988 .

[12]  K. A. Gallivan,et al.  Parallel Algorithms for Dense Linear Algebra Computations , 1990, SIAM Rev..

[13]  D. Siljak,et al.  On the convergence of parallel asynchronous block-iterative computations , 1990 .

[14]  Mariesa L. Crow,et al.  The parallel implementation of the waveform relaxation method for transient stability simulations , 1990 .

[15]  A. Bhaya,et al.  Asynchronous block-iterative methods for almost linear equations , 1991 .

[16]  D. Siljak,et al.  Nested epsilon decompositions of linear systems: weakly coupled and overlapping blocks , 1991 .