FETI-DP algorithms for 2D Biot model with discontinuous Galerkin discretization

Pilhwa Lee

2023-06-22

Abstract:Dual-Primal Finite Element Tearing and Interconnecting (FETI-DP) algorithms are developed for a 2D Biot model. The model is formulated with mixed-finite elements as a saddle-point problem. The displacement $\mathbf{u}$ and the Darcy flux flow $\mathbf{z}$ are represented with $P_1$ piecewise continuous elements and pore-pressure $p$ with $P_0$ piecewise constant elements, {\it i.e.}, overall three fields with a stabilizing term. We have tested the functionality of FETI-DP with Dirichlet preconditioners. Numerical experiments show a signature of scalability of the resulting parallel algorithm in the compressible elasticity with permeable Darcy flow as well as almost incompressible elasticity.

Numerical Analysis,Analysis of PDEs

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to develop and validate the Dual - Primal Finite Element Tearing and Interconnecting (FETI - DP) algorithm for solving the two - dimensional Biot model. Specifically, the main research problems include:

1. **Discretization of the linear poroelastic model**: Represent the displacement \( \mathbf{u} \) and Darcy flow \( \mathbf{z} \) as piecewise continuous elements and the pore pressure \( p \) as a piecewise constant element through the mixed finite element method, thereby constructing a saddle - point problem.

2. **Numerical stability and convergence**: Test the functionality and scalability of the FETI - DP algorithm for compressible and nearly incompressible elastic problems and permeable Darcy flow. In particular, study how to accelerate the iterative solution process through the Dirichlet pre - conditioner.

3. **Elastic locking and non - physical oscillatory pressure distribution**: Deal with the problems of elastic locking and non - physical oscillatory pressure distribution common in the traditional continuous Galerkin method to ensure the stability and accuracy of the numerical solution.

4. **Scalability of parallel computing**: Verify the parallel performance of the proposed FETI - DP algorithm in a distributed memory system, especially its scalability when the number of sub - domains is increased.

The paper demonstrates the effectiveness and scalability of the FETI - DP algorithm in solving the two - dimensional Biot model through a series of numerical experiments. In particular, it performs well in dealing with compressible and nearly incompressible elastic problems and solves some key problems in traditional methods.