TY - JOUR
T1 - New hybridized mixed methods for linear elasticity and optimal multilevel solvers
AU - Gong, Shihua
AU - Wu, Shuonan
AU - Xu, Jinchao
N1 - Funding Information:
Acknowledgements The work of the first and third authors was supported in part by National Natural Science Foundation of China (NSFC) (Grant Nos. 91430215, 41390452) and by Beijing International Center for Mathematical Research of Peking University, China. The work of the second and third authors was supported in part by the DOE Grant DE-SC0009249 as part of the Collaboratory on Mathematics for Mesoscopic Modeling of Materials and by DOE Grant DE-SC0014400 and NSF Grant DMS-1522615.
Publisher Copyright:
© 2018, Springer-Verlag GmbH Germany, part of Springer Nature.
PY - 2019/2/13
Y1 - 2019/2/13
N2 - In this paper, we present a family of new mixed finite element methods for linear elasticity for both spatial dimensions n= 2 , 3 , which yields a conforming and strongly symmetric approximation for stress. Applying Pk + 1- Pk as the local approximation for the stress and displacement, the mixed methods achieve the optimal order of convergence for both the stress and displacement when k≥ n. For the lower order case (n- 2 ≤ k< n) , the stability and convergence still hold on some special grids. The proposed mixed methods are efficiently implemented by hybridization, which imposes the inter-element normal continuity of the stress by a Lagrange multiplier. Then, we develop and analyze multilevel solvers for the Schur complement of the hybridized system in the two dimensional case. Provided that no nearly singular vertex on the grids, the proposed solvers are proved to be uniformly convergent with respect to both the grid size and Poisson’s ratio. Numerical experiments are provided to validate our theoretical results.
AB - In this paper, we present a family of new mixed finite element methods for linear elasticity for both spatial dimensions n= 2 , 3 , which yields a conforming and strongly symmetric approximation for stress. Applying Pk + 1- Pk as the local approximation for the stress and displacement, the mixed methods achieve the optimal order of convergence for both the stress and displacement when k≥ n. For the lower order case (n- 2 ≤ k< n) , the stability and convergence still hold on some special grids. The proposed mixed methods are efficiently implemented by hybridization, which imposes the inter-element normal continuity of the stress by a Lagrange multiplier. Then, we develop and analyze multilevel solvers for the Schur complement of the hybridized system in the two dimensional case. Provided that no nearly singular vertex on the grids, the proposed solvers are proved to be uniformly convergent with respect to both the grid size and Poisson’s ratio. Numerical experiments are provided to validate our theoretical results.
UR - http://www.scopus.com/inward/record.url?scp=85055722816&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85055722816&partnerID=8YFLogxK
U2 - 10.1007/s00211-018-1001-3
DO - 10.1007/s00211-018-1001-3
M3 - Article
AN - SCOPUS:85055722816
SN - 0029-599X
VL - 141
SP - 569
EP - 604
JO - Numerische Mathematik
JF - Numerische Mathematik
IS - 2
ER -