Computational & Technology Resources
an online resource for computational,
engineering & technology publications
PROCEEDINGS OF THE SIXTH INTERNATIONAL CONFERENCE ON ENGINEERING COMPUTATIONAL TECHNOLOGY
Edited by: M. Papadrakakis and B.H.V. Topping
Order Reduction of Linear Time-Invariant Finite-Integration-Technique Discretized Structures with Length and Frequency Parameterization
K.S. Stavrakakis1, T. Wittig2, W. Ackermann1 and T. Weiland1
1Institute for Electromagnetic Field Theory, TU Darmstadt, Germany
K.S. Stavrakakis, T. Wittig, W. Ackermann, T. Weiland, "Order Reduction of Linear Time-Invariant Finite-Integration-Technique Discretized Structures with Length and Frequency Parameterization", in M. Papadrakakis, B.H.V. Topping, (Editors), "Proceedings of the Sixth International Conference on Engineering Computational Technology", Civil-Comp Press, Stirlingshire, UK, Paper 125, 2008. doi:10.4203/ccp.89.125
Keywords: model order reduction, multivariate Krylov subspace methods, parameterized systems, finite integration technique.
In this paper we consider large linear time invariant systems which depend on several parameters s1,s2,...,sn. These may be the frequency, material or geometric parameters. It is presented how they arise from the FIT-discretized [1,2] Maxwell equations. In order to solve such systems in a reasonable time, model order reduction is used which aims to reduce the n dimensional system to an m<<n dimensional system which preserves some qualities of the original system. Especially model order reduction via projection on Krylov-subspaces has proven to be very effective and its characteristics for one-parameter linear systems are briefly highlighted here.
In the single parameter case [3,4] the order reduction is carried out via projection on a subspace with dimension m<<n. The solution vector x of the original system is restricted on a space V spanned by orthonormal trial vectors v1...v2. By projection with a test matrix W in Cl*m the dimension of the system matrix is reduced from n*n to m*m. The matrices V and W are chosen so that the resulting reduced order model is a partial realization, Padé approximation or rational interpolant of the transfer function of the original system.
Details concerning the relation of the approximations and the projection matrices are found in  as well as in [3,4]. It is a common characteristic of all these cases, that colsp(V) and colsp(W) must contain unions of Krylov subspaces of the system matrices and vectors . An analogous method exists for multivariate linear systems (see also [6,7,8]) and is presented also in this paper. Unfortunately it cannot be used for FIT-systems directly as they are nonlinear with respect to line length. This work shows a linearization method. We consider the mesh being stretched according to the length variation, so that the mesh topology is maintained and expresses all matrices according to this mesh variation. The procedure presented here is demonstrated on a numerical example.
purchase the full-text of this paper (price £20)