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


Kaczmarz-Type Methods for Solving Matrix Equations

In this paper, several Kaczmarz-type numerical methods for solving the m...

Truncation Preconditioners for Stochastic Galerkin Finite Element Discretizations

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

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

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

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 Unified Approach for Solving Sequential Selection Problems

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

Novel ensemble algorithms for random two-domain parabolic problems

In this paper, three efficient ensemble algorithms are proposed for fast...

Symbol based convergence analysis in multigrid methods for saddle point problems

Saddle point problems arise in a variety of applications, e.g., when sol...

Please sign up or login with your details

Forgot password? Click here to reset