Efficient stochastic Galerkin methods for random diffusion equations
- Department of Mathematics, Purdue University, 150 N. University Street, West Lafayette, IN 47907 (United States), E-mail: dxiu@math.purdue.edu
- Department of Mathematics, Purdue University, 150 N. University Street, West Lafayette, IN 47907 (United States), E-mail: shen@math.purdue.edu
We discuss in this paper efficient solvers for stochastic diffusion equations in random media. We employ generalized polynomial chaos (gPC) expansion to express the solution in a convergent series and obtain a set of deterministic equations for the expansion coefficients by Galerkin projection. Although the resulting system of diffusion equations are coupled, we show that one can construct fast numerical methods to solve them in a decoupled fashion. The methods are based on separation of the diagonal terms and off-diagonal terms in the matrix of the Galerkin system. We examine properties of this matrix and show that the proposed method is unconditionally stable for unsteady problems and convergent for steady problems with a convergent rate independent of discretization parameters. Numerical examples are provided, for both steady and unsteady random diffusions, to support the analysis.
- OSTI ID:
- 21167746
- Journal Information:
- Journal of Computational Physics, Vol. 228, Issue 2; Other Information: DOI: 10.1016/j.jcp.2008.09.008; PII: S0021-9991(08)00477-4; Copyright (c) 2008 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved; Country of input: International Atomic Energy Agency (IAEA); ISSN 0021-9991
- Country of Publication:
- United States
- Language:
- English
Similar Records
A well-posed and stable stochastic Galerkin formulation of the incompressible Navier–Stokes equations with random data
DG-IMEX Stochastic Galerkin Schemes for Linear Transport Equation with Random Inputs and Diffusive Scalings