The hp-mortar finite-element method for the mixed elasticity and Stokes problems

被引:15
作者
Ben Belgacem, F
Chilton, LK
Seshaiyer, P
机构
[1] Univ Toulouse 3, CNRS, INSAT, Math Ind & Phys UMR 5640,UPS,UT1, F-31062 Toulouse 04, France
[2] USAF, Inst Technol, Dept Math & Stat, Wright Patterson AFB, OH 45433 USA
[3] Texas Tech Univ, Dept Math & Stat, Lubbock, TX 79409 USA
关键词
D O I
10.1016/S0898-1221(03)90079-4
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
The motivation of this work is to apply the hp-version of the mortar finite-element method to the nearly incompressible elasticity model formulated as a mixed displacement-pressure problem as well as to Stokes equations in primal velocity-pressure variables. Within each subdomain, the local approximation is designed using div-stable hp-mixed finite elements. The displacement is computed in a mortared space, while the pressure is not subjected to any constraints across the interfaces. By a Boland-Nicolaides argument, we prove that the discrete saddle-point problem satisfies a Babuska-Brezzi inf-sup condition. The inf-sup constant is optimal in the sense that it depends only on the local (to the subdomains) characteristics of the mixed finite elements and, in particular, it does not increase with the total number of the subdomains. The consequences, that we are aware of, such an important result are twofold. The numerical analysis of the approximability properties of the hp-mortar discretization for the mixed elasticity problem allows us to derive an asymptotic rate of convergence that is optimal up to rootlogp in the displacement; this is addressed in the present contribution. When the mortar discrete problem is inverted by substructured iterative methods based on Krylov subspaces with block preconditioners, in view of the results for conforming finite elements [1], the condition number of the solver should grow logarithmically on (p, h) and not depend on the total number of the subdomains. (C) 2003 Elsevier Science Ltd. All rights reserved.
引用
收藏
页码:35 / 55
页数:21
相关论文
共 39 条
[1]   Iterative substructuring preconditioners for mortar element methods in two dimensions [J].
Achdou, Y ;
Maday, Y ;
Widlund, OB .
SIAM JOURNAL ON NUMERICAL ANALYSIS, 1999, 36 (02) :551-580
[2]  
[Anonymous], 1975, SOBOLEV SPACES
[3]  
[Anonymous], RAIRO MATH MODEL NUM
[4]  
Arrow KJ., 1958, Studies in Nonlinear Programming
[5]   THE P AND H-P VERSIONS OF THE FINITE-ELEMENT METHOD, BASIC PRINCIPLES AND PROPERTIES [J].
BABUSKA, I ;
SURI, M .
SIAM REVIEW, 1994, 36 (04) :578-632
[6]   THE OPTIMAL CONVERGENCE RATE OF THE P-VERSION OF THE FINITE-ELEMENT METHOD [J].
BABUSKA, I ;
SURI, M .
SIAM JOURNAL ON NUMERICAL ANALYSIS, 1987, 24 (04) :750-776
[7]   LOCKING EFFECTS IN THE FINITE-ELEMENT APPROXIMATION OF ELASTICITY PROBLEMS [J].
BABUSKA, I ;
SURI, M .
NUMERISCHE MATHEMATIK, 1992, 62 (04) :439-463
[8]  
Belgacem FB, 2000, NUMER MATH, V85, P257, DOI 10.1007/s002110000136
[9]   Optimal convergence rates of hp mortar finite element methods for second-order elliptic problems [J].
Ben Belgacem, F ;
Seshaiyer, P ;
Suri, M .
ESAIM-MATHEMATICAL MODELLING AND NUMERICAL ANALYSIS-MODELISATION MATHEMATIQUE ET ANALYSE NUMERIQUE, 2000, 34 (03) :591-608
[10]   Mixed mortar finite element method for the incompressible stokes problem: Convergence analysis [J].
Ben Belgacem, F .
SIAM JOURNAL ON NUMERICAL ANALYSIS, 2000, 37 (04) :1085-1100