An efficient iterative method for solving parameter-dependent and random diffusion problems

by   Xiaobing Feng, et al.

This paper develops and analyzes a general iterative framework for solving parameter-dependent and random diffusion problems. It is inspired by the multi-modes method of [7,8] and the ensemble method of [19] and extends those methods into a more general and unified framework. The main idea of the framework is to reformulate the underlying problem into another problem with a parameter-independent diffusion coefficient and a parameter-dependent (and solution-dependent) right-hand side, a fixed-point iteration is then employed to compute the solution of the reformulated problem. The main benefit of the proposed approach is that an efficient direct solver and a block Krylov subspace iterative solver can be used at each iteration, allowing to reuse the LU matrix factorization or to do an efficient matrix-matrix multiplication for all parameters, which in turn results in significant computation saving. Convergence and rates of convergence are established for the iterative method both at the variational continuous level and at the finite element discrete level under some structure conditions. Several strategies for establishing reformulations of parameter-dependent and random diffusion problems are proposed and their computational complexity is analyzed. Several 1-D and 2-D numerical experiments are also provided to demonstrate the efficiency of the proposed iterative method and to validate the theoretical convergence results.


page 1

page 2

page 3

page 4


Truncation Preconditioners for Stochastic Galerkin Finite Element Discretizations

Stochastic Galerkin finite element method (SGFEM) provides an efficient ...

A class of GADI methods for time-dependent linear systems with multitask kernel-learning parameter prediction

This paper develops a class of general alternating-direction implicit (G...

A modified Kačanov iteration scheme with application to quasilinear diffusion models

The classical Kačanov scheme for the solution of nonlinear variational p...

An efficient nonlinear solver and convergence analysis for a viscoplastic flow model

This paper studies a finite element discretization of the regularized Bi...

A Unified Approach for Solving Sequential Selection Problems

In this paper we develop a unified approach for solving a wide class of ...

An efficient multi-modes Monte Carlo homogenization method for random materials

In this paper, we propose and analyze a new stochastic homogenization me...

Computational issues by interpolating with inverse multiquadrics: a solution

We consider the interpolation problem with the inverse multiquadric radi...