Constructing adaptive generalized polynomial chaos method to measure the uncertainty in continuous models: A computational approach

Due to errors in measurements and inherent variability in the quantities of interest, models based on random differential equations give more realistic results than their deterministic counterpart. The generalized polynomial chaos (gPC) is a powerful technique used to approximate the solution of these equations when the random inputs follow standard probability distributions. But in many cases these random inputs do not have a standard probability distribution. In this paper, we present a step-by-step constructive methodology to implement directly a useful version of adaptive gPC for arbitrary distributions, extending the applicability of the gPC. The paper mainly focuses on the computational aspects, on the implementation of the method and on the creation of a useful software tool. This tool allows the user to easily change the types of distributions and the order of the expansions, and to study their effects on the convergence and on the results. Several examples illustrating the usefulness of the method are included.

[1]  R. H. Myers,et al.  STAT 319 : Probability & Statistics for Engineers & Scientists Term 152 ( 1 ) Final Exam Wednesday 11 / 05 / 2016 8 : 00 – 10 : 30 AM , 2016 .

[2]  George Adomian,et al.  Solving Frontier Problems of Physics: The Decomposition Method , 1993 .

[3]  P. Kloeden,et al.  Numerical Solution of Stochastic Differential Equations , 1992 .

[4]  I. Gihman,et al.  Stochastic Differential Equations , 1975 .

[5]  Michel Loève,et al.  Probability Theory I , 1977 .

[6]  D. Williams.  STOCHASTIC DIFFERENTIAL EQUATIONS: THEORY AND APPLICATIONS , 1976 .

[7]  Matthew P. Wand,et al.  Kernel Smoothing , 1995 .

[8]  J. Simonoff.  Smoothing Methods in Statistics , 1998 .

[9]  Juan Carlos Cortés,et al.  Mean square power series solution of random linear differential equations , 2010, Comput. Math. Appl..

[10]  Ali H. Nayfeh,et al.  Problems in Perturbation , 1985 .

[11]  A. Weber,et al.  Modeling epidemics caused by respiratory syncytial virus (RSV). , 2001, Mathematical biosciences.

[12]  W. T. Martin,et al.  The Orthogonal Development of Non-Linear Functionals in Series of Fourier-Hermite Functionals , 1947 .

[13]  E. Renshaw,et al.  STOCHASTIC DIFFERENTIAL EQUATIONS , 1974 .

[14]  N. Bershad,et al.  Random differential equations in science and engineering , 1975, Proceedings of the IEEE.

[15]  Nicola Bellomo,et al.  Nonlinear Stochastic Systems in Physics and Mechanics , 1987 .

[16]  Wilfred D. Iwan,et al.  On the existence and uniqueness of solutions generated by equivalent linearization , 1978 .

[17]  O. Ernst,et al.  ON THE CONVERGENCE OF GENERALIZED POLYNOMIAL CHAOS EXPANSIONS , 2011 .

[18]  Habib N. Najm,et al.  Numerical Challenges in the Use of Polynomial Chaos Representations for Stochastic Processes , 2005, SIAM J. Sci. Comput..

[19]  R. Ghanem,et al.  Stochastic Finite Elements: A Spectral Approach , 1990 .

[20]  Dirk P. Kroese,et al.  Handbook of Monte Carlo Methods , 2011 .

[21]  Juan Carlos Cortés,et al.  Numerical solution of random differential initial value problems: Multistep methods , 2011 .

[22]  N. Wiener.  The Homogeneous Chaos , 1938 .

[23]  Goodarz Ahmadi,et al.  Application of Wiener-Hermite Expansion to Nonstationary Random Vibration of a Duffing Oscillator , 1983 .

[24]  Dongbin Xiu,et al.  The Wiener-Askey Polynomial Chaos for Stochastic Differential Equations , 2002, SIAM J. Sci. Comput..

[25]  George E. Karniadakis,et al.  Time-dependent generalized polynomial chaos , 2010, J. Comput. Phys..