solidot新版网站常见问题,请点击这里查看。

A Bramble-Pasciak conjugate gradient method for discrete Stokes problems with lognormal random viscosity. (arXiv:1712.06472v1 [math.NA])

来源于:arXiv
We study linear systems of equations arising from a stochastic Galerkin finite element discretization of saddle point problems with random data and its iterative solution. We consider the Stokes flow model with random viscosity described by the exponential of a correlated Gaussian process and shortly discuss the discretization framework and the representation of the emerging matrix equation. Due to the high dimensionality and the coupling of the associated symmetric, indefinite, linear system, we resort to iterative solvers and problem-specific preconditioners. As a standard iterative solver for this problem class, we consider the block diagonal preconditioned MINRES method and further introduce the Bramble-Pasciak conjugate gradient method as a promising alternative. This special conjugate gradient method is formulated in a non-standard inner product with a block triangular preconditioner. From a structural point of view, such a block triangular preconditioner enables a better approxi 查看全文>>