A Krylov enhanced proper orthogonal decomposition method for efficient nonlinear model reduction

A model order reduction method is proposed for approximating nonlinear partial differential equations (PDEs). The method attempts to combine the desirable attributes of Krylov reduction and proper orthogonal decomposition (POD) and is entitled Krylov enhanced POD (KPOD). The method approximates nonlinear input/output (I/O) behavior using a sequence of state dependent Krylov subspaces which are obtained via simulation of a nonlinear finite element discretized system. POD is then implemented to extract, from the sequence of subspaces, a projection basis that describes the state dependent variance of the Krylov subspaces. Due to Krylov's moment matching property, the variance describes how the I/O of a sequence of Krylov reduced models would vary as the original model's state varies. Galerkin projection is performed using the extracted basis to yield a low dimensional I/O approximation of the original nonlinear model. Reduced order models of an electro-thermal switch generated using the conventional POD and KPOD are presented and used to compare the I/O approximation capabilities of each method. Our results indicate that KPOD reduced models are capable of accurately modeling I/O using a sequence of Krylov subspaces generated from a single nonlinear discretized model simulation while POD can require multiple nonlinear discretized model simulations to generate a reduced model with similar I/O accuracy. The proposed method can be a potential asset for practical applications in the area of I/O behavior prediction and design optimization of nonlinear systems.

[1]  J. Korvink,et al.  MST MEMS model order reduction: Requirements and benchmarks , 2006 .

[2]  Zu-Qing Qu,et al.  Model Order Reduction Techniques with Applications in Finite Element Analysis , 2004 .

[3]  Siak Piang Lim,et al.  PROPER ORTHOGONAL DECOMPOSITION AND ITS APPLICATIONS – PART II: MODEL REDUCTION FOR MEMS DYNAMICAL ANALYSIS , 2002 .

[4]  C. Rowley,et al.  Modeling of transitional channel flow using balanced proper orthogonal decomposition , 2007, 0707.4112.

[5]  Adam Fic,et al.  Solving Transient Nonlinear Heat Conduction Problems by Proper Orthogonal Decomposition and the Finite-Element Method , 2005 .

[6]  G. Kerschen,et al.  The Method of Proper Orthogonal Decomposition for Dynamical Characterization and Order Reduction of Mechanical Systems: An Overview , 2005 .

[7]  J. Marsden,et al.  Structure-preserving Model Reduction of Mechanical Systems , 2000 .

[8]  Siak Piang Lim,et al.  Proper orthogonal decomposition and component mode synthesis in macromodel generation for the dynamic simulation of a complex MEMS device , 2003 .

[9]  J. Korvink,et al.  Dynamic electro-thermal simulation of microsystems—a review , 2005 .

[10]  Serkan Gugercin,et al.  Smith-Type Methods for Balanced Truncation of Large Sparse Systems , 2005 .

[11]  Nisar Ahmed,et al.  Dimensionally reduced Krylov subspace model reduction for large scale systems , 2007, Appl. Math. Comput..

[12]  Wilkins Aquino,et al.  An object-oriented framework for reduced-order models using proper orthogonal decomposition (POD) , 2007 .

[13]  Jan G. Korvink,et al.  Efficient optimization of transient dynamic problems in MEMS devices using model order reduction , 2005 .

[14]  D. V. Lyubimov,et al.  Thermosolutal convection in a horizontal porous layer heated from below in the presence of a horizontal through flow , 2008 .

[15]  M. Parameswaran,et al.  Modelling surface-micromachined electrothermal actuators , 2004, Canadian Journal of Electrical and Computer Engineering.

[16]  Khalide Jbilou,et al.  Matrix Krylov subspace methods for large scale model reduction problems , 2006, Appl. Math. Comput..

[17]  B. Moore Principal component analysis in linear systems: Controllability, observability, and model reduction , 1981 .

[18]  Z. Bai Krylov subspace techniques for reduced-order modeling of large-scale dynamical systems , 2002 .

[19]  R. Cook,et al.  Concepts and Applications of Finite Element Analysis , 1974 .

[20]  Mats G. Larson,et al.  Adaptive finite element approximation of multiphysics problems , 2007 .

[21]  Jan G. Korvink,et al.  Review: Automatic Model Reduction for Transient Simulation of MEMS‐based Devices , 2002 .

[22]  J. Korvink,et al.  Model order reduction of 3D electro-thermal model for a novel micromachined hotplate gas sensor , 2004, 5th International Conference on Thermal and Mechanical Simulation and Experiments in Microelectronics and Microsystems, 2004. EuroSimE 2004. Proceedings of the.

[23]  Gene H. Golub,et al.  Matrix computations , 1983 .

[24]  M. Gunzburger,et al.  Reduced-order modeling of time-dependent PDEs with multiple parameters in the boundary data , 2007 .

[25]  A. Antoulas,et al.  A Survey of Model Reduction by Balanced Truncation and Some New Results , 2004 .

[26]  Xiaolin Chen,et al.  Coupled electrothermal-mechanical analysis for MEMS via model order reduction , 2010 .

[27]  Kurt Busch,et al.  A Krylov‐subspace based solver for the linear and nonlinear Maxwell equations , 2007 .

[28]  Hamid Zahrouni,et al.  A model reduction method for the post-buckling analysis of cellular microstructures , 2007 .

[29]  Steven C. Chapra,et al.  Numerical Methods for Engineers , 1986 .

[30]  S. D. Senturia,et al.  Generating efficient dynamical models for microelectromechanical systems from a few finite-element simulation runs , 1999 .

[31]  R. Guyan Reduction of stiffness and mass matrices , 1965 .

[32]  Clarence W. Rowley,et al.  Model Reduction for fluids, Using Balanced Proper Orthogonal Decomposition , 2005, Int. J. Bifurc. Chaos.