20,800 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
Towards tensor-based methods for the numerical approximation of the Perron-Frobenius and Koopman operator
The global behavior of dynamical systems can be studied by analyzing the
eigenvalues and corresponding eigenfunctions of linear operators associated
with the system. Two important operators which are frequently used to gain
insight into the system's behavior are the Perron-Frobenius operator and the
Koopman operator. Due to the curse of dimensionality, computing the
eigenfunctions of high-dimensional systems is in general infeasible. We will
propose a tensor-based reformulation of two numerical methods for computing
finite-dimensional approximations of the aforementioned infinite-dimensional
operators, namely Ulam's method and Extended Dynamic Mode Decomposition (EDMD).
The aim of the tensor formulation is to approximate the eigenfunctions by
low-rank tensors, potentially resulting in a significant reduction of the time
and memory required to solve the resulting eigenvalue problems, provided that
such a low-rank tensor decomposition exists. Typically, not all variables of a
high-dimensional dynamical system contribute equally to the system's behavior,
often the dynamics can be decomposed into slow and fast processes, which is
also reflected in the eigenfunctions. Thus, the weak coupling between different
variables might be approximated by low-rank tensor cores. We will illustrate
the efficiency of the tensor-based formulation of Ulam's method and EDMD using
simple stochastic differential equations
- …