1,280 research outputs found
A Bramble-Pasciak conjugate gradient method for discrete Stokes equations with random viscosity
We study the iterative solution of linear systems of equations arising from
stochastic Galerkin finite element discretizations of saddle point problems. We
focus on the Stokes model with random data parametrized by uniformly
distributed random variables and discuss well-posedness of the variational
formulations. We introduce a Bramble-Pasciak conjugate gradient method as a
linear solver. It builds on a non-standard inner product associated with a
block triangular preconditioner. The block triangular structure enables more
sophisticated preconditioners than the block diagonal structure usually applied
in MINRES methods. We show how the existence requirements of a conjugate
gradient method can be met in our setting. We analyze the performance of the
solvers depending on relevant physical and numerical parameters by means of
eigenvalue estimates. For this purpose, we derive bounds for the eigenvalues of
the relevant preconditioned sub-matrices. We illustrate our findings using the
flow in a driven cavity as a numerical test case, where the viscosity is given
by a truncated Karhunen-Lo\`eve expansion of a random field. In this example, a
Bramble-Pasciak conjugate gradient method with block triangular preconditioner
outperforms a MINRES method with block diagonal preconditioner in terms of
iteration numbers.Comment: 19 pages, 1 figure, submitted to SIAM JU
Factorizing the Stochastic Galerkin System
Recent work has explored solver strategies for the linear system of equations
arising from a spectral Galerkin approximation of the solution of PDEs with
parameterized (or stochastic) inputs. We consider the related problem of a
matrix equation whose matrix and right hand side depend on a set of parameters
(e.g. a PDE with stochastic inputs semidiscretized in space) and examine the
linear system arising from a similar Galerkin approximation of the solution. We
derive a useful factorization of this system of equations, which yields bounds
on the eigenvalues, clues to preconditioning, and a flexible implementation
method for a wide array of problems. We complement this analysis with (i) a
numerical study of preconditioners on a standard elliptic PDE test problem and
(ii) a fluids application using existing CFD codes; the MATLAB codes used in
the numerical studies are available online.Comment: 13 pages, 4 figures, 2 table
IGA-based Multi-Index Stochastic Collocation for random PDEs on arbitrary domains
This paper proposes an extension of the Multi-Index Stochastic Collocation
(MISC) method for forward uncertainty quantification (UQ) problems in
computational domains of shape other than a square or cube, by exploiting
isogeometric analysis (IGA) techniques. Introducing IGA solvers to the MISC
algorithm is very natural since they are tensor-based PDE solvers, which are
precisely what is required by the MISC machinery. Moreover, the
combination-technique formulation of MISC allows the straight-forward reuse of
existing implementations of IGA solvers. We present numerical results to
showcase the effectiveness of the proposed approach.Comment: version 3, version after revisio
Efficient stochastic finite element methods for flow in heterogeneous porous media. Part 2: random lognormal permeability
Efficient and robust iterative methods are developed for solving the linear systems of equations arising from stochastic finite element methods for single phase fluid flow in porous media. Permeability is assumed to vary randomly in space according to some given correlation function. In the companion paper, herein referred to as Part 1, permeability was approximated using a truncated Karhunen‐Loève expansion (KLE). The stochastic variability of permeability is modeled using lognormal random fields and the truncated KLE is projected onto a polynomial chaos basis. This results in a stochastic nonlinear problem since the random fields are represented using polynomial chaos containing terms that are generally nonlinear in the random variables. Symmetric block Gauss‐Seidel used as a preconditioner for CG is shown to be efficient and robust for stochastic finite element method
Hierarchical Schur complement preconditioner for the stochastic Galerkin finite element methods
Use of the stochastic Galerkin finite element methods leads to large systems
of linear equations obtained by the discretization of tensor product solution
spaces along their spatial and stochastic dimensions. These systems are
typically solved iteratively by a Krylov subspace method. We propose a
preconditioner which takes an advantage of the recursive hierarchy in the
structure of the global matrices. In particular, the matrices posses a
recursive hierarchical two-by-two structure, with one of the submatrices block
diagonal. Each one of the diagonal blocks in this submatrix is closely related
to the deterministic mean-value problem, and the action of its inverse is in
the implementation approximated by inner loops of Krylov iterations. Thus our
hierarchical Schur complement preconditioner combines, on each level in the
approximation of the hierarchical structure of the global matrix, the idea of
Schur complement with loops for a number of mutually independent inner Krylov
iterations, and several matrix-vector multiplications for the off-diagonal
blocks. Neither the global matrix, nor the matrix of the preconditioner need to
be formed explicitly. The ingredients include only the number of stiffness
matrices from the truncated Karhunen-Lo\`{e}ve expansion and a good
preconditioned for the mean-value deterministic problem. We provide a condition
number bound for a model elliptic problem and the performance of the method is
illustrated by numerical experiments.Comment: 15 pages, 2 figures, 9 tables, (updated numerical experiments
- …