Analysis of a 2-field finite element solver for poroelasticity on quadrilateral meshes

JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS

卷:393

Article

Wang, Zhuoran¹  Tavener, Simon²  Liu, Jiangguo² 

[1] Sun Yat Sen Univ, Sch Math Zhuhai, Zhuhai 519082, Guangdong, Peoples R China

[2] Colorado State Univ, Dept Math, Ft Collins, CO 80523 USA

关键词: Arbogast-Correa spaces;    Enriched Lagrangian finite elements;    Poroelasticity;    Quadrilateral meshes;    Weak Galerkin;   

DOI  :  10.1016/j.cam.2021.113539

来源: Elsevier

【 摘 要 】

This paper presents a novel 2-field finite element solver for linear poroelasticity on convex quadrilateral meshes. The Darcy flow is discretized for fluid pressure by a lowest-order weak Galerkin (WG) finite element method, which establishes the discrete weak gradient and numerical velocity in the lowest-order Arbogast-Correa space. The linear elasticity is discretized for solid displacement by the enriched Lagrangian finite elements with a special treatment for the volumetric dilation. These two types of finite elements are coupled through the implicit Euler temporal discretization to solve poroelasticity problems. A rigorous error analysis is presented along with numerical tests to demonstrate the accuracy and locking-free property of this new solver. (c) 2021 Elsevier B.V. All rights reserved.

【 授权许可】

Free