Smoothed aggregation algebraic multigrid for stochastic PDE problems with layered materials

SUMMARY We consider two algebraic multilevel solvers for the solution of discrete problems arising from PDEs with random inputs. Our focus is on problems with large jumps in material coefficients. The model problem considered is that of a diffusion problem with uncertainties in the diffusion coefficients and realization values differing dramatically between different layers in the spatial domain, with the location of the interfaces between the layers assumed to be known. The stochastic discretization is based on the generalized polynomial chaos, and the spatial problem is discretized using conforming finite elements. A multigrid solver based on smoothed aggregation is presented, and numerical experiments are provided demonstrating convergence properties of the multigrid solver. The observed convergence is shown to depend only weakly on the stochastic discretization and, for the discretized stochastic PDE problem, generally mimics that of the corresponding deterministic problem. Copyright © 2014 John Wiley & Sons, Ltd.

[1]  J. Mandel,et al.  An Iterative Method with Convergence Rate Chosen a priori , 1999 .

[2]  Gianluca Iaccarino,et al.  A least-squares approximation of partial differential equations with high-dimensional random inputs , 2009, J. Comput. Phys..

[3]  Roger Ghanem,et al.  Numerical solution of spectral stochastic finite element systems , 1996 .

[4]  Marian Brezina,et al.  Convergence of algebraic multigrid based on smoothed aggregation , 1998, Numerische Mathematik.

[5]  C. Schwab,et al.  Sparse high order FEM for elliptic sPDEs , 2009 .

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

[7]  Habib N. Najm,et al.  A multigrid solver for two-dimensional stochastic diffusion equations , 2003 .

[8]  H. Matthies,et al.  Hierarchical parallelisation for the solution of stochastic finite element equations , 2005 .

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

[10]  D. Xiu,et al.  Modeling Uncertainty in Steady State Diffusion Problems via Generalized Polynomial Chaos , 2002 .

[11]  Roger G. Ghanem,et al.  Hierarchical Schur complement preconditioner for the stochastic Galerkin finite element methods , 2012, Numer. Linear Algebra Appl..

[12]  Stefan Vandewalle,et al.  Iterative Solvers for the Stochastic Finite Element Method , 2008, SIAM J. Sci. Comput..

[13]  Houman Owhadi,et al.  A non-adapted sparse approximation of PDEs with stochastic inputs , 2010, J. Comput. Phys..

[14]  Elisabeth Ullmann,et al.  A Kronecker Product Preconditioner for Stochastic Galerkin Finite Element Discretizations , 2010, SIAM J. Sci. Comput..

[15]  Howard C. Elman,et al.  Solving the Stochastic Steady-State Diffusion Problem using , 2006 .

[16]  Marian Brezina,et al.  Algebraic multigrid by smoothed aggregation for second and fourth order elliptic problems , 2005, Computing.

[17]  Stefan Vandewalle,et al.  Algebraic multigrid for stationary and time‐dependent partial differential equations with stochastic coefficients , 2008, Numer. Linear Algebra Appl..

[18]  Howard C. Elman,et al.  Block-diagonal preconditioning for spectral stochastic finite-element systems , 2008 .