PRECONDITIONED NULLSPACE METHOD FOR THE TWO-DIMENSIONAL OSEEN PROBLEM

被引:11
|
作者
Le Borne, Sabine [1 ]
机构
[1] Tennessee Technol Univ, Dept Math, Cookeville, TN 38505 USA
来源
SIAM JOURNAL ON SCIENTIFIC COMPUTING | 2009年 / 31卷 / 04期
基金
美国国家科学基金会;
关键词
nullspace method; saddle point problem; QR factorization; preconditioning; DIVERGENCE-FREE BASES; CONSTRUCTION; COMPUTATION; ELEMENT;
D O I
10.1137/070691577
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
The Oseen problem, which arises in the simulation of the time-dependent Navier-Stokes equations for incompressible fluid flow, leads to indefinite, nonsymmetric and possibly ill-conditioned linear systems of equations. This paper presents a method for obtaining a reduced linear system from the original system, which is then solved by the preconditioned BiCGStab method. The system reduction is obtained through an efficient implicit representation of a basis of discretely divergence-free functions, also known as the nullspace method. We will show a close relationship between this preconditioned nullspace method and the BFBt-preconditioned pressure Schur complement problem arising in typical block preconditioners. We will use this relationship to analyze the spectrum of the preconditioned reduced system. We will also present numerical tests for the two-dimensional Oseen problem to illustrate the performance of this method, which is similar to the performance of the BFBt-preconditioned pressure Schur complement solver. In particular, it is robust w.r.t. the strength of the viscosity for simple constant wind but not for more complicated flows, and it shows only a slight dependence on the mesh size.
引用
收藏
页码:2494 / 2509
页数:16
相关论文
共 50 条