Domain decomposition preconditioners for non-selfconjugate second order elliptic problems

被引:0
作者
Zhang, HY [1 ]
Sun, JC [1 ]
机构
[1] Chinese Acad Sci, Inst Software, Beijing, Peoples R China
来源
NUMERICAL TREATMENT OF MULTIPHASE FLOWS IN POROUS MEDIA | 2000年 / 552卷
关键词
non-selfconjugate; elliptic equation; domain decomposition; Schur complement; preconditioner;
D O I
暂无
中图分类号
O59 [应用物理学];
学科分类号
摘要
A non-symmetric interface Schur complement arises from non-selfconjugate second order elliptic problems with domain decomposition methods. The usual numerical methods for solving it axe GMRES, ORTHOMIN, and BICGSTAB, but they take a large amount of computer time and memory. The authors find in this paper that the nonsymmetric Schur complement can in fact be changed into a symmetric one by scaling. Then an efficient preconditioner can be provided by which the preconditioned system can be solved iteratively by a modified PCG method. When the problem is imposed on a rectangular region, the condition number is estimated and is nearly one. Numerical experiments axe also presented. Non-selfconjugate problems arise in mathematical modeling and numerical simulation of fluid flows and transport in porous media.
引用
收藏
页码:409 / 418
页数:10
相关论文
共 6 条
[1]  
BRAMBLE JH, 1992, CORNELL MATH DEP LEC
[2]  
CHAN TF, 1990, THIRD INTERNATIONAL SYMPOSIUM ON DOMAIN DECOMPOSITION METHODS FOR PARTIAL DIFFERENTIAL EQUATIONS, P245
[3]  
CHAN TF, 1992, 6 INT C DOM DEC SIAM, P157
[4]   A CAPACITANCE MATRIX-METHOD FOR DIRICHLET PROBLEM ON POLYGON REGION [J].
DRYJA, M .
NUMERISCHE MATHEMATIK, 1982, 39 (01) :51-64
[5]  
SAAD Y, 1982, ITERATIVE METHODS SP
[6]  
Smith B., 1996, Domain decomposition