《求解二阶混合有限体元离散系统高效预条件子》由会员分享，可在线阅读，更多相关《求解二阶混合有限体元离散系统高效预条件子（35页珍藏版）》请在金锄头文库上搜索。
1、?2011-4-20Efficient Preconditioners for SolvingSecond Order Mixed-type Finite VolumeElement Discretization SystemsCandidateXu SenlinSupervisorProfessor Shu ShiCollegeMathematics and Computational ScienceProgramThe Computational MathematicsSpecializationMultigrid Method and Domain Decomposition Metho
2、dDegreeMaster of ScienceUniversityXiangtan UniversityDate2011-4-20? ? ? ? ?(FVE)?.?,?,?.?,?,?- ILU?AMG?PGMRES(m)?.?,?,?.?,?.?,?:?.? ?,?,?.?,?,?PGMRES(m)?,?,?m?.? ?;?; PGMRES(m)?; ILU?; AMG?.IAbstractFinite Volume Element (FVE) method is one of the discretization methods for partialdifferential equat
3、ions. Second order mixed-type finite volume element method which wasdeveloped recently has many advantages, such as being able to preserve local conservationof certain physical quantities. Fast algorithms for the corresponding discretization systemswill be discussed in this paper.Firstly, ILU- and A
4、MG-GMRES(m) are employed to solve the linear systems arisingfrom second order mixed-type finite volume element scheme for an elliptic problem withjump coefficients. Numerical results show that the numbers of iteration of the two PGM-RES (m) are unstable which strongly depend on not only the mesh siz
5、e but also the jumpcoefficient. Therefore, it is quite necessary to develop newly efficient preconditioners.Secondly,twonewpreconditionersaredesignedfortheabovelinearsystems,i.e., blockdiagonal and two-level preconditioner. For consistent mesh, theoretical analysis of the for-mer is given, and it is
6、 proved that the spectral condition number of the preconditioned sys-tem is uniformly bounded. Numerical results show that our preconditioners are stable andeffi cient. The numbers of iteration of the corresponding PGMRES (m) are significantly re-duced, furthermore, they are independent of mesh size
7、 and insensitiveto jump coefficient orrestart parameter.Key Words: finite volume element; preconditioner; PGMRES(m) method; ILU factoriza-tion; AMG method.II?.1?.32.1?.32.2?.4?PGMRES(m)?.83.1 PGMRES(m)?.83.2?.11?.154.1?.154.2?.21?.25?.26?.30?.31III?7, 11, 19, 24, 2729, 35,?(?)?,?.?(FVE)?.? ,?FVE?25,
8、 26, 31, 40, 45, 48?,?FVE?.?,?FVE?8, 41, 49?.?8?,?,?,?,?.?,?,?,?.?,?,?GMRES(m)?4, 5, 39?,?,?,?.?,?.? ?:?LU (Incomplete LU Factorization, ILU)?1, 12, 15, 18, 22, 30, 33,?(Algebraic Multigrid, AMG)?2, 3, 6, 9, 10, 13, 14, 16, 20, 21, 32, 36, 37, 43, 44, 46, 47, 50,?,?.?, ILU?,?ILU(k)?ILUT?,?ILU(0)?,?. AMG?,?,?,?.?ILU(0)?AMG?PGMRES(m)?,?,?PGMRES(m)?,?,?,?.?,?,?m.?,?,?.1?,?A?,?.?,?,?A?Ad= diag(A11,A22)?,?A11?A22?.?A11,?,?AMG?,?AMG?A11?;?A22,?diag(A22)?,?,?
