Convergence of a cell-centered finite volume discretization for linear elasticity

We show convergence of a cell-centered finite volume discretization for linear elasticity. The discretization, termed the MPSA method, was recently proposed in the context of geological applications, where cell-centered variables are often preferred. Our analysis utilizes a hybrid variational formulation, which has previously been used to analyze finite volume discretizations for the scalar diffusion equation. The current analysis deviates significantly from previous in three respects. First, additional stabilization leads to a more complex saddle-point problem. Secondly, a discrete Korn's inequality has to be established for the global discretization. Finally, robustness with respect to the Poisson ratio is analyzed. The stability and convergence results presented herein provide the first rigorous justification of the applicability of cell-centered finite volume methods to problems in linear elasticity

Stable cell-centered finite volume discretization for Biot equations

In this paper we discuss a new discretization for the Biot equations. The discretization treats the coupled system of deformation and flow directly, as opposed to combining discretizations for the two separate sub-problems. The coupled discretization has the following key properties, the combination of which is novel: 1) The variables for the pressure and displacement are co-located, and are as sparse as possible (e.g. one displacement vector and one scalar pressure per cell center). 2) With locally computable restrictions on grid types, the discretization is stable with respect to the limits of incompressible fluid and small time-steps. 3) No artificial stabilization term has been introduced. Furthermore, due to the finite volume structure embedded in the discretization, explicit local expressions for both momentum-balancing forces as well as mass-conservative fluid fluxes are available. We prove stability of the proposed method with respect to all relevant limits. Together with consistency, this proves convergence of the method. Finally, we give numerical examples verifying both the analysis and convergence of the method

Mixed-dimensional models for real-world applications

We explore mathematical models for physical problems in which it is necessary to simultaneously consider equations in different dimensions; these are called mixed-dimensional models. We first give several examples, and then an overview of recent progress made towards finding a general method of solution of such problems

Adaptive asynchronous time-stepping, stopping criteria, and a posteriori error estimates for fixed-stress iterative schemes for coupled poromechanics problems

In this paper we develop adaptive iterative coupling schemes for the Biot system modeling coupled poromechanics problems. We particularly consider the space-time formulation of the fixed-stress iterative scheme, in which we first solve the problem of flow over the whole space-time interval, then exploiting the space-time information for solving the mechanics. Two common discretizations of this algorithm are then introduced based on two coupled mixed finite element methods in-space and the backward Euler scheme in-time. Therefrom, adaptive fixed-stress algorithms are build on conforming reconstructions of the pressure and displacement together with equilibrated flux and stresses reconstructions. These ingredients are used to derive a posteriori error estimates for the fixed-stress algorithms, distinguishing the different error components, namely the spatial discretization, the temporal discretization, and the fixed-stress iteration components. Precisely, at the iteration $k\geq 1$ of the adaptive algorithm, we prove that our estimate gives a guaranteed and fully computable upper bound on the energy-type error measuring the difference between the exact and approximate pressure and displacement. These error components are efficiently used to design adaptive asynchronous time-stepping and adaptive stopping criteria for the fixed-stress algorithms. Numerical experiments illustrate the efficiency of our estimates and the performance of the adaptive iterative coupling algorithms

Stable mixed finite elements for linear elasticity with thin inclusions

We consider mechanics of composite materials in which thin inclusions are modeled by lower-dimensional manifolds. By successively applying the dimensional reduction to junctions and intersections within the material, a geometry of hierarchically connected manifolds is formed which we refer to as mixed-dimensional. The governing equations with respect to linear elasticity are then defined on this mixed-dimensional geometry. The resulting system of partial differential equations is also referred to as mixed-dimensional, since functions defined on domains of multiple dimensionalities are considered in a fully coupled manner. With the use of a semi-discrete differential operator, we obtain the variational formulation of this system in terms of both displacements and stresses. The system is then analyzed and shown to be well-posed with respect to appropriately weighted norms. Numerical discretization schemes are proposed using well-known mixed finite elements in all dimensions. The schemes conserve linear momentum locally while relaxing the symmetry condition on the stress tensor. Stability and convergence are shown using a priori error estimates and confirmed numerically.publishedVersio

High-accuracy phase-field models for brittle fracture based on a new family of degradation functions

Phase-field approaches to fracture based on energy minimization principles have been rapidly gaining popularity in recent years, and are particularly well-suited for simulating crack initiation and growth in complex fracture networks. In the phase-field framework, the surface energy associated with crack formation is calculated by evaluating a functional defined in terms of a scalar order parameter and its gradients, which in turn describe the fractures in a diffuse sense following a prescribed regularization length scale. Imposing stationarity of the total energy leads to a coupled system of partial differential equations, one enforcing stress equilibrium and another governing phase-field evolution. The two equations are coupled through an energy degradation function that models the loss of stiffness in the bulk material as it undergoes damage. In the present work, we introduce a new parametric family of degradation functions aimed at increasing the accuracy of phase-field models in predicting critical loads associated with crack nucleation as well as the propagation of existing fractures. An additional goal is the preservation of linear elastic response in the bulk material prior to fracture. Through the analysis of several numerical examples, we demonstrate the superiority of the proposed family of functions to the classical quadratic degradation function that is used most often in the literature.Comment: 33 pages, 30 figure

Domain decomposition preconditioning for non-linear elasticity problems

We consider domain decomposition techniques for a non-linear elasticity problem. Our main focus is on non-linear preconditioning, realized in the framework of additive Schwarz preconditioned inexact Newton (ASPIN) methods. The standard 1-level ASPIN method is extended to a 2-level method by adding a non-linear coarse solver. Numerical experiments show that the coarse component is necessary for scalability in terms of linear iterations inside the Newton loop. Moreover, for problems that are dominated by nonlinearities that are not localized in space the non-linear coarse iterations are crucial for achieving computational efficiency.publishedVersio