Convergence acceleration of polynomial chaos solutions via sequence transformation

Abstract We investigate convergence acceleration of the solution of stochastic differential equations characterized by their polynomial chaos expansions. Specifically, nonlinear sequence transformations are adapted to these expansions, viewed as a one-parameter family of functions with the parameter being the polynomial degree of the expansion. These transformations can be generally viewed as nonlinear generalizations of Richardson Extrapolation and permit the estimation of coefficients in higher order expansions having knowledge of the coefficients in lower order ones. Stochastic Galerkin closure that typically characterizes the solution of such equations yields polynomial chaos representations that have the requisite analytical properties to ensure suitable convergence of these nonlinear sequence transformations. We investigate specifically Shanks and Levin transformations, and explore their properties in the context of a stochastic initial value problem and a stochastic elliptic problem.

[1]  Claude Brezinski,et al.  Extrapolation algorithms and Pade´ approximations: a historical survey , 1996 .

[2]  Raúl Tempone,et al.  Galerkin Finite Element Approximations of Stochastic Elliptic Partial Differential Equations , 2004, SIAM J. Numer. Anal..

[3]  D. Shanks Non‐linear Transformations of Divergent and Slowly Convergent Sequences , 1955 .

[4]  William B. Gragg,et al.  On Extrapolation Algorithms for Ordinary Initial Value Problems , 1965 .

[5]  Baskar Ganapathysubramanian,et al.  Sparse grid collocation schemes for stochastic natural convection problems , 2007, J. Comput. Phys..

[6]  I. Babuska,et al.  Solving elliptic boundary value problems with uncertain coefficients by the finite element method: the stochastic formulation , 2005 .

[7]  Claude Brezinski,et al.  Convergence acceleration during the 20th century , 2000 .

[8]  Christian Soize,et al.  Computational Aspects for Constructing Realizations of Polynomial Chaos in High Dimension , 2010, SIAM J. Sci. Comput..

[9]  M. Lemaire,et al.  Stochastic finite element: a non intrusive approach by regression , 2006 .

[10]  John Red-Horse,et al.  Propagation of probabilistic uncertainty in complex physical systems using a stochastic finite element approach , 1999 .

[11]  Roger Ghanem,et al.  Stochastic model reduction for chaos representations , 2007 .

[12]  O. L. Maître,et al.  Spectral Methods for Uncertainty Quantification: With Applications to Computational Fluid Dynamics , 2010 .

[13]  L. Richardson The Approximate Arithmetical Solution by Finite Differences of Physical Problems Involving Differential Equations, with an Application to the Stresses in a Masonry Dam , 1911 .

[14]  R. Ghanem,et al.  Stochastic Finite-Element Analysis of Seismic Soil-Structure Interaction , 2002 .

[15]  Roger Ghanem,et al.  Characterization of reservoir simulation models using a polynomial chaos‐based ensemble Kalman filter , 2009 .

[16]  Roger G. Ghanem,et al.  Physical Systems with Random Uncertainties: Chaos Representations with Arbitrary Probability Measure , 2005, SIAM J. Sci. Comput..

[17]  B. Burrows,et al.  Lower bounds for quartic anharmonic and double‐well potentials , 1993 .

[18]  A. C. Aitken XXV.—On Bernoulli's Numerical Solution of Algebraic Equations , 1927 .

[19]  Christian Soize,et al.  Random field representations for stochastic elliptic boundary value problems and statistical inverse problems , 2013, European Journal of Applied Mathematics.

[20]  J. Stoer,et al.  Numerical treatment of ordinary differential equations by extrapolation methods , 1966 .

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

[22]  Christian Soize,et al.  A computational inverse method for identification of non-Gaussian random fields using the Bayesian approach in very high dimension , 2011, Computer Methods in Applied Mechanics and Engineering.

[23]  Roger G. Ghanem,et al.  Asymptotic Sampling Distribution for Polynomial Chaos Representation of Data: A Maximum Entropy and Fisher information approach , 2006, Proceedings of the 45th IEEE Conference on Decision and Control.

[24]  R. Ghanem,et al.  Multi-resolution analysis of wiener-type uncertainty propagation schemes , 2004 .

[25]  R. Ghanem,et al.  Iterative solution of systems of linear equations arising in the context of stochastic finite elements , 2000 .

[26]  J. Nash,et al.  PARABOLIC EQUATIONS. , 1957, Proceedings of the National Academy of Sciences of the United States of America.

[27]  Roger G. Ghanem,et al.  On the construction and analysis of stochastic models: Characterization and propagation of the errors associated with limited data , 2006, J. Comput. Phys..

[28]  Jacques Dutka,et al.  Richardson extrapolation and Romberg integration , 1984 .

[29]  G. Karniadakis,et al.  Multi-Element Generalized Polynomial Chaos for Arbitrary Probability Measures , 2006, SIAM J. Sci. Comput..

[30]  David Levin,et al.  Development of non-linear transformations for improving convergence of sequences , 1972 .

[31]  Ernst Joachim Weniger,et al.  Nonlinear sequence transformations for the acceleration of convergence and the summation of divergent series , 1989 .

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

[33]  Roger Ghanem,et al.  A general purpose library for stochastic finite element computations , 2003 .

[34]  Christian Soize,et al.  Maximum likelihood estimation of stochastic chaos representations from experimental data , 2006 .

[35]  F. Benth,et al.  Convergence Rates for Finite Element Approximations of Stochastic Partial Differential Equations , 1998 .