Mixed finite element and TPSA finite volume methods for linearized elasticity and Cosserat materials

Jan Martin Nordbotten,Wietse M. Boon,Omar Duran,Eirik Keilegavlen

2024-09-20

Abstract:Cosserat theory of elasticity is a generalization of classical elasticity that allows for asymmetry in the stress tensor by taking into account micropolar rotations in the medium. The equations involve a rotation field and associated "couple stress" as variables, in addition to the conventional displacement and Cauchy stress fields. In recent work, we derived a mixed finite element method (MFEM) for the linear Cosserat equations that converges optimally in these four variables. The drawback of this method is that it retains the stresses as unknowns, and therefore leads to relatively large saddle point system that are computationally demanding to solve. As an alternative, we developed a finite volume method in which the stress variables are approximated using a minimal, two-point stencil (TPSA). The system consists of the displacement and rotation variables, with an additional "solid pressure" unknown. Both the MFEM and TPSA methods are robust in the incompressible limit and in the Cauchy limit, for which the Cosserat equations degenerate to classical linearized elasticity. We report on the construction of the methods, their a priori properties, and compare their numerical performance against an MPSA finite volume method.

Numerical Analysis

What problem does this paper attempt to address?

The main problems that this paper attempts to solve are how to effectively discretize and solve the linear elastic equations and Cosserat material equations in geological reservoir simulation. Specifically, the paper explores the following two key issues:

1. **Improving the Robustness and Efficiency of Numerical Methods**:
- When dealing with linear elastic problems, the Finite Volume Methods (FVM) have advantages over the traditional Finite Element Methods (FEM) such as being compatible with the mesh structure and being able to explicitly and robustly represent the traction on internal surfaces (such as fractures).
- In particular, ensure the robustness of the numerical method in the incompressible limit (\(\lambda \to \infty\)) and the Cauchy limit (\(\ell \to 0\)).

2. **Introducing New Discretization Methods**:
- A Finite Volume Method based on the Two - Point Stress Approximation (TPSA) for linear elastic equations and Cosserat material equations is proposed. This method reduces the computational complexity by simplifying the stress approximation while maintaining robustness in limit cases.
- The performance of the Mixed Finite Element Method (MFEM) and the Multi - Point Stress Approximation (MPSA) methods are compared and analyzed.

### Specific Problem Description

#### Linear Elastic Equations
Linear elastic equations are usually discretized by the Finite Element Method, but in some cases, the Finite Volume Method has unique advantages. For example, when dealing with fluid flow equations in underground applications, the Finite Volume Method can better adapt to the mesh structure and provide a more explicit representation of the internal surface traction.

#### Cosserat Material Equations
Cosserat elastic theory is an extension of the classical elastic theory, which allows the asymmetry of the stress tensor and takes into account the micropolar rotation in the medium. This is especially important when dealing with problems such as granular media. The paper proposes a Mixed Finite Element Method and a Finite Volume Method based on TPSA to solve the Cosserat material equations and compares their performance.

### Numerical Experiments
To verify the effectiveness and robustness of these methods, the paper conducts three numerical experiments:
1. **Linear Elastic Medium with Uniform Parameters**: Evaluate the robustness of the method in the incompressible limit.
2. **Linear Elastic Medium with Parameter Heterogeneity**: Study the performance of the method in the case of discontinuous material coefficients.
3. **Composite Cosserat/Elastic Material**: Consider the case where some areas are linear elastic and some areas are Cosserat materials, and evaluate the convergence of the method under different material parameters.

Through these experiments, the paper shows the convergence and robustness of different methods under different conditions, especially in limit cases.

In summary, this paper aims to develop and compare efficient and robust numerical methods suitable for linear elastic equations and Cosserat material equations, with particular attention to the application of the Finite Volume Method in these fields.