Cost reduction of stochastic Galerkin method by adaptive identification of significant polynomial chaos bases for elliptic equations

Pranesh, Srikara and Ghosh, Debraj (2018) Cost reduction of stochastic Galerkin method by adaptive identification of significant polynomial chaos bases for elliptic equations. In: COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING, 340 . pp. 54-69.

PDF Com_Met_App_Cha_Eng_340_54_2018.pdf - Published Version Restricted to Registered users only Download (545kB) | Request a copy

Abstract

One widely used and computationally efficient method for uncertainty quantification using spectral stochastic finite element is the stochastic Galerkin method. Here the solution is represented in polynomial chaos expansion, and the residual of the discretized governing equation is projected on the polynomial chaos bases. This results in a system of deterministic algebraic equations with the polynomials chaos coefficients as unknown. However, one impediment for its large scale applications is the curse of dimensionality, that is, the exponential growth of the number of polynomial chaos bases with the stochastic dimensionality and degree of expansion. Here, for a stochastic elliptic problem, an adaptive selection of polynomial chaos bases is proposed. Accordingly, during the first few iterations in the preconditioned conjugate gradient method for solving the system of linear algebraic equations, the chaos bases with maximal contribution -in an appropriately defined metric - to the solution are first identified. Subsequently, only these bases are retained for further iterations until convergence is achieved. Using numerical studies a three times cost saving over the existing method is observed. Furthermore, for enhancing the computational cost gain, the stochastic Galerkin method is reformulated as a generalized Sylvester equation. This step allowed efficient usage of the sparsity of moments of product of polynomial chaos bases. Through numerical studies on problems with large stochastic dimensionality, an additional cost saving of up to one order of magnitude -twenty times -is observed. This amounts to sixty times speedup over the existing method, when adaptive selection and generalized Sylvester equation formulation are used together. The proposed methodology can be easily incorporated in an existing standard stochastic Galerkin method solver for elliptic problems. (C) 2018 Elsevier B.V. All rights reserved.

Item Type:	Journal Article
Publication:	COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING
Publisher:	ELSEVIER SCIENCE SA
Additional Information:	Copy right for this article belong to ELSEVIER SCIENCE SA, PO BOX 564, 1001 LAUSANNE, SWITZERLAND
Keywords:	Uncertainty quantification; Stochastic finite element; Stochastic Galerkin method; Conjugate gradient; Generalized Sylvester equation; Polynomial chaos expansion
Department/Centre:	Division of Mechanical Sciences > Civil Engineering
Date Deposited:	11 Sep 2018 14:57
Last Modified:	06 Nov 2018 07:13
URI:	http://eprints.iisc.ac.in/id/eprint/60641

Actions (login required)

View Item