A Low-rank solver for the Navier--Stokes equations with uncertain viscosity

Kookjin Lee,Howard C. Elman,Bedřich Sousedík
DOI: https://doi.org/10.1137/17M1151912
2017-10-17
Abstract:We study an iterative low-rank approximation method for the solution of the steady-state stochastic Navier--Stokes equations with uncertain viscosity. The method is based on linearization schemes using Picard and Newton iterations and stochastic finite element discretizations of the linearized problems. For computing the low-rank approximate solution, we adapt the nonlinear iterations to an inexact and low-rank variant, where the solution of the linear system at each nonlinear step is approximated by a quantity of low rank. This is achieved by using a tensor variant of the GMRES method as a solver for the linear systems. We explore the inexact low-rank nonlinear iteration with a set of benchmark problems, using a model of flow over an obstacle, under various configurations characterizing the statistical features of the uncertain viscosity, and we demonstrate its effectiveness by extensive numerical experiments.
Numerical Analysis
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: **How to efficiently solve the steady - state Navier - Stokes equations under random viscosity conditions**. Specifically, the author focuses on the efficient computation of the steady - state Navier - Stokes equations in the case of uncertain viscosity. This uncertainty may stem from measurement errors or the uncertainty of the proportion of each phase in a multiphase medium. The paper proposes an iterative low - rank approximation method to solve this problem, which is based on the linearization schemes of Picard and Newton iterations and combines the stochastic finite - element discretization technique. Through this method, the author can significantly reduce the computational cost while ensuring accuracy. ### Problem Background 1. **Navier - Stokes Equations**: The fundamental equations that describe fluid motion. 2. **Random Viscosity**: Viscosity, as a positive random field, depends on a set of random variables, which makes the problem random. 3. **High - Dimensional Problem**: As the number of random variables increases, the dimension of the problem will increase rapidly, leading to a substantial increase in computational complexity. ### Solutions To address these problems, the author proposes the following methods: - **Low - Rank Approximation**: Utilize the low - rank tensor format to represent and solve the linearized system, thereby reducing the amount of computation. - **GMRES Method**: Use the Generalized Minimum Residual Method (GMRES) in tensor format to solve the linear system. - **Preconditioning**: Introduce preconditioners to accelerate convergence. - **Inexact Nonlinear Iteration**: During the nonlinear iteration process, dynamically adjust the solution accuracy of the linear system according to the accuracy of the current solution to further improve efficiency. ### Main Contributions 1. **Low - Rank Approximation Framework**: Propose a low - rank approximation framework suitable for stochastic Navier - Stokes equations, which can effectively handle high - dimensional stochastic problems. 2. **Efficient Solver**: Develop an efficient solver based on the tensor Krylov subspace method, which significantly reduces the computational cost. 3. **Numerical Verification**: Verify the effectiveness and superiority of the proposed method through a series of benchmark tests. In summary, this paper aims to solve the problem of efficiently solving the Navier - Stokes equations under random viscosity conditions and provides an innovative low - rank approximation method for this purpose.