Polynomial chaos representation of a stochastic preconditioner

Christophe Desceliers (1) , R. Ghanem (2) , Christian Soize (1)

1 LaM - Laboratoire de Mécanique
2 USC - Department of aerospace and mechanical engineering

A method is developed in this paper to accelerate the convergence in computing the solution of stochastic algebraic systems of equations. The method is based on computing, via statistical sampling, a polynomial chaos decomposition of a stochastic preconditioner to the system of equations. This preconditioner can subsequently be used in conjunction with either chaos representations of the solution or with approaches based on Monte Carlo sampling. In addition to presenting the supporting theory, the paper also presents a convergence analysis and an example to demonstrate the significance of the proposed algorithm.

Mots clés: uncertainty quantification; polynomial chaos; stochastic systems; preconditioners; linear solvers; incomplete factorization

Domaines: Mécanique [physics.med-ph]; Probabilités [math.PR]; Statistiques [math.ST]

Christophe Desceliers, R. Ghanem, Christian Soize. Polynomial chaos representation of a stochastic preconditioner. International Journal for Numerical Methods in Engineering, 2005, 64 (5), pp.618-634. ⟨10.1002/nme.1382⟩. ⟨hal-00686183⟩