Computational & Technology Resources
an online resource for computational,
engineering & technology publications
PROCEEDINGS OF THE FIFTH INTERNATIONAL CONFERENCE ON ENGINEERING COMPUTATIONAL TECHNOLOGY
Edited by: B.H.V. Topping, G. Montero and R. Montenegro
A Parallel Model Reduction Algorithm for Efficient Power Flow Prediction in Vibrating Systems
Department of Naval Architecture and Ocean Engineering, Dalian University of Technology, China
X. Li, "A Parallel Model Reduction Algorithm for Efficient Power Flow Prediction in Vibrating Systems", in B.H.V. Topping, G. Montero, R. Montenegro, (Editors), "Proceedings of the Fifth International Conference on Engineering Computational Technology", Civil-Comp Press, Stirlingshire, UK, Paper 116, 2006. doi:10.4203/ccp.84.116
Keywords: power flow, Krylov subspace, model reduction, parallel implementation, adaptive interpolation.
This paper focuses on power flow prediction in vibrating systems over a frequency band of interest. An efficient method is proposed to predict power flow from finite element models based on the rational Krylov subspace projection . It belongs to a family of well-established moment-matching techniques, which have been successfully used in circuit simulation , electromagnetic problems  and acoustic field computation [4,5]. In this paper, general formulation of power flow is derived as multiple-input multiple-output transfer functions of linear time-invariant systems. Instead of constructing the reduced-order models (ROMs) by directly projecting the original systems onto the rational Krylov subspaces, a matrix-free algorithm is developed to construct the ROMs from the interpolation values of the transfer functions and their higher-order derivatives . The ROMs match the moments of the transfer functions of the original systems at interpolation frequencies to a specified order, which guarantees the preservation of power flow and its higher-order derivatives up to the same order. The matrix-free algorithm reduces both the computational cost and the memory requirement in constructing the ROMs.
The accuracy of the ROMs is mainly determined by the choice of interpolation frequencies. This paper investigates how to improve the accuracy of the ROMs by adaptive interpolation. After an initial ROM is constructed, its accuracy is improved by adding new interpolation points at the frequencies corresponding to large approximation error. Unfortunately, the exact error is not known without solving the original system. The results presented in this paper show that the residual norms provide a good way to identify the frequencies with large approximation error. Since the residual norms are obtained directly from the ROM, this approach minimizes the approximation error at a reasonable cost. The performance of the matrix-free algorithm can be enhanced by adapting it for simultaneous execution in multiple processes. The main process adaptively chooses interpolation frequencies according to the current error estimate, evenly distributes them among all processes to achieve load balancing, gathers computation results from the rest processes, and constructs the ROM according to the matrix-free algorithm. Inter-process communication is implemented by using the message passing interface . Significant parallel speedup is achieved because the major computational task in each process is independent.
Case studies are conducted on the testbeds from Harwell-Boeing collection . Four numerical examples are presented to illustrate the efficiency of adaptive interpolation, the effect of matching-moment order, and the effect of different combinations of initial and updating sizes. For the ROMs of the same size, the one using more interpolation frequencies is usually more accurate. Initial and updating sizes also affect the accuracy of the resulting ROMs. It is observed that starting from a small-sized ROM and iterating many times with small updating size does not necessarily guarantee an accurate ROM. It is preferable to use a large initial size and update only one or two times in order to construct an accurate ROM efficiently.
purchase the full-text of this paper (price £20)