Stability-preserving model order reduction for linear stochastic Galerkin systems

被引:0
作者
Roland Pulch
机构
[1] University of Greifswald,Institute of Mathematics and Computer Science
来源
Journal of Mathematics in Industry | / 9卷
关键词
Linear dynamical system; Polynomial chaos; Stochastic Galerkin method; Model order reduction; Asymptotic stability; Lyapunov equation; 65L05; 65L20; 65L80; 34C20; 34D20;
D O I
暂无
中图分类号
学科分类号
摘要
Mathematical modeling often yields linear dynamical systems in science and engineering. We change physical parameters of the system into random variables to perform an uncertainty quantification. The stochastic Galerkin method yields a larger linear dynamical system, whose solution represents an approximation of random processes. A model order reduction (MOR) of the Galerkin system is advantageous due to the high dimensionality. However, asymptotic stability may be lost in some MOR techniques. In Galerkin-type MOR methods, the stability can be guaranteed by a transformation to a dissipative form. Either the original dynamical system or the stochastic Galerkin system can be transformed. We investigate the two variants of this stability-preserving approach. Both techniques are feasible, while featuring different properties in numerical methods. Results of numerical computations are demonstrated for two test examples modeling a mechanical application and an electric circuit, respectively.
引用
收藏
相关论文
共 50 条
[1]  
Benner P(2015)A survey of projection-based model order reduction methods for parametric dynamical systems SIAM Rev 57 483-531
[2]  
Gugercin S(2012)Stability preservation in projection-based model order reduction of large scale systems Eur J Control 18 122-132
[3]  
Willcox K(2012)On the convergence of generalized polynomial chaos expansions ESAIM: M2AN 46 317-339
[4]  
Castañé Selga R(2016)Fast and accurate model reduction for spectral methods in uncertainty quantification Int J Uncertain Quantificat 6 271-286
[5]  
Lohmann B(2003)Model reduction methods based on Krylov subspaces Acta Numer 12 267-319
[6]  
Eid R(1982)Numerical solution of stable non-negative definite Lyapunov equation IMA J Numer Anal 2 303-323
[7]  
Ernst OG(1975)The modified nodal approach to network analysis IEEE Trans Circuits Syst 22 504-509
[8]  
Mugler A(2016)A POD projection method for large-scale algebraic Riccati equations Numer Algebra Control Optim 6 413-435
[9]  
Starkloff HJ(1991)Solution of Lyapunov equations by alternating direction implicit iteration Comput Math Appl 21 43-58
[10]  
Ullmann E(2013)On the passivity of polynomial chaos-based augmented models for stochastic circuits IEEE Trans Circuits Syst I, Regul Pap 60 2998-3007