DeepAI AI Chat
Log In Sign Up

A Low-rank Solver for the Stochastic Unsteady Navier-Stokes Problem

by   Howard C. Elman, et al.
University of Maryland

We study a low-rank iterative solver for the unsteady Navier-Stokes equations for incompressible flows with a stochastic viscosity. The equations are discretized using the stochastic Galerkin method, and we consider an all-at-once formulation where the algebraic systems at all the time steps are collected and solved simultaneously. The problem is linearized with Picard's method. To efficiently solve the linear systems at each step, we use low-rank tensor representations within the Krylov subspace method, which leads to significant reductions in storage requirements and computational costs. Combined with effective mean-based preconditioners and the idea of inexact solve, we show that only a small number of linear iterations are needed at each Picard step. The proposed algorithm is tested with a model of flow in a two-dimensional symmetric step domain with different settings to demonstrate the computational efficiency.


page 1

page 2

page 3

page 4


Krylov techniques for low-rank ADI

One of the most computationally expensive steps of the low-rank ADI meth...

Stochastic Discontinuous Galerkin Methods with Low–Rank Solvers for Convection Diffusion Equations

We investigate numerical behaviour of a convection diffusion equation wi...

Low-rank solutions to the stochastic Helmholtz equation

In this paper, we consider low-rank approximations for the solutions to ...

Stochastic dynamical modeling of turbulent flows

Advanced measurement techniques and high performance computing have made...

A neural network multigrid solver for the Navier-Stokes equations

We present a deep neural network multigrid solver (DNN-MG) that we devel...

Bayesian inversion for electromyography using low-rank tensor formats

The reconstruction of the structure of biological tissue using electromy...