Jan Martin Nordbotten

Jan Martin Nordbotten

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.

Eirik Keilegavlen, Jan Martin Nordbotten

We introduce a new cell-centered finite volume discretization for elasticity with weakly enforced symmetry of the stress tensor. The method is motivated by the need for robust discretization methods for deformation and flow in porous media, and falls in the category of multi-point stress approximations (MPSA). By enforcing symmetry weakly, the resulting method has additional flexibility beyond previous MPSA methods. This allows for a construction of a method which is applicable to all grid types, and in particular the method amends a crucial shortcoming in previous MPSA methods for simplex grids. By formulating the method as a discrete variational problem, we prove convergence of the new method for a wide range of problems, with conditions that can be verified at the time of discretization. We present the first set of comprehensive numerical tests for the MPSA methods in three dimensions, covering Cartesian and simplex grids, with both heterogeneous and nearly incompressible media. The tests show that the new method consistently is second order convergent in displacement, despite being lowest order, with a rate that mostly is between 1 and 2 for stresses. The results further show that the new method is more robust and computationally cheaper than previous MPSA methods.

Erlend Storvik, Jakub Wiktor Both, Kundan Kumar et al.

In this work we are interested in effectively solving the quasi-static, linear Biot model for poromechanics. We consider the fixed-stress splitting scheme, which is a popular method for iteratively solving Biot's equations. It is well-known that the convergence of the method is strongly dependent on the applied stabilization/tuning parameter. In this work, we propose a new approach to optimize this parameter. We show theoretically that it depends also on the fluid flow properties and not only on the mechanics properties and the coupling coefficient. The type of analysis presented in this paper is not restricted to a particular spatial discretization. We only require it to be inf-sup stable. The convergence proof applies also to low-compressible or incompressible fluids and low-permeable porous media. Illustrative numerical examples, including random initial data, random boundary conditions or random source terms and a well-known benchmark problem, i.e. Mandel's problem are performed. The results are in good agreement with the theoretical findings. Furthermore, we show numerically that there is a connection between the inf-sup stability of discretizations and the performance of the fixed-stress splitting scheme.

Jakub Wiktor Both, Kundan Kumar, Jan Martin Nordbotten et al.

In this paper, we study the robust linearization of nonlinear poromechanics of unsaturated materials. The model of interest couples the Richards equation with linear elasticity equations, employing the equivalent pore pressure. In practice a monolithic solver is not always available, defining the requirement for a linearization scheme to allow the use of separate simulators, which is not met by the classical Newton method. We propose three different linearization schemes incorporating the fixed-stress splitting scheme, coupled with an L-scheme, Modified Picard and Newton linearization of the flow. All schemes allow the efficient and robust decoupling of mechanics and flow equations. In particular, the simplest scheme, the Fixed-Stress-L-scheme, employs solely constant diagonal stabilization, has low cost per iteration, and is very robust. Under mild, physical assumptions, it is theoretically shown to be a contraction. Due to possible break-down or slow convergence of all considered splitting schemes, Anderson acceleration is applied as post-processing. Based on a special case, we justify theoretically the general ability of the Anderson acceleration to effectively accelerate convergence and stabilize the underlying scheme, allowing even non-contractive fixed-point iterations to converge. To our knowledge, this is the first theoretical indication of this kind. Theoretical findings are confirmed by numerical results. In particular, Anderson acceleration has been demonstrated to be very effective for the considered Picard-type methods. Finally, the Fixed-Stress-Newton scheme combined with Anderson acceleration provides a robust linearization scheme, meeting the above criteria.

Mats Kirkesæther Brun, Thomas Wick, Inga Berre et al.

This paper concerns the analysis and implementation of a novel iterative staggered scheme for quasi-static brittle fracture propagation models, where the fracture evolution is tracked by a phase field variable. The model we consider is a two-field variational inequality system, with the phase field function and the elastic displacements of the solid material as independent variables. Using a penalization strategy, this variational inequality system is transformed into a variational equality system, which is the formulation we take as the starting point for our algorithmic developments. The proposed scheme involves a partitioning of this model into two subproblems; phase field and mechanics, with added stabilization terms to both subproblems for improved efficiency and robustness. We analyze the convergence of the proposed scheme using a fixed point argument, and find that under a natural condition, the elastic mechanical energy remains bounded, and, if the diffusive zone around crack surfaces is sufficiently thick, monotonic convergence is achieved. Finally, the proposed scheme is validated numerically with several bench-mark problems.

Mats Kirkesæther Brun, Elyes Ahmed, Inga Berre et al.

This paper concerns splitting-based iterative procedures for the coupled nonlinear thermo-poroelasticity model problem. The thermo-poroelastic model problem we consider is formulated as a three-field system of PDE's, consisting of an energy balance equation, a mass balance equation and a momentum balance equation, where the primary variables are temperature, fluid pressure, and elastic displacement. Due to the presence of a nonlinear convective transport term in the energy balance equation, it is convenient to have access to both the pressure and temperature gradients. Hence, we introduce these as two additional variables and extend the original three-field model to a five-field model. For the numerical solution of this five-field formulation, we compare three approaches that differ by how we treat the coupling/decoupling between the flow and/from heat and/from mechanics; these approaches have in common a simultaneous application of the fixed-stress splitting scheme on both the non-linearity and the coupling structure of the problem. More precisely, the derived procedures transform a nonlinear and fully coupled problem into a set of simpler subproblems to be solved sequentially in an iterative fashion. We provide a convergence proof for the derived algorithms, and validate our results through several numerical examples.

Eren Ucar, Eirik Keilegavlen, Inga Berre et al.

Simulating the deformation of fractured media requires the coupling of different models for the deformation of fractures and the formation surrounding them. We consider a cell-centered finite-volume approach, termed the multipoint stress approximation (MPSA) method, which is developed in order to discretize coupled flow and mechanical deformation in the subsurface. Within the MPSA framework, we consider fractures as co-dimension one inclusions in the domain, with the fracture surfaces represented as line pairs in 2D (faces in 3D) that displace relative to each other. Fracture deformation is coupled to that of the surrounding domain through internal boundary conditions. This approach is natural within the finite-volume framework, where tractions are defined on surfaces of the grid. The MPSA method is capable of modeling deformation considering open and closed fractures with complex and nonlinear relationships governing the displacements and tractions at the fracture surfaces. We validate our proposed approach using both problems for which analytical solutions are available and more complex benchmark problems, including comparison with a finite-element discretization.

Manuel Antonio Borregales, Kundan Kumar, Jan Martin Nordbotten et al.

We consider L-scheme and Newton based solvers for Biot model under small or large deformation. The mechanical deformation follows the Saint Venant-Kirchoff constitutive law. Further, the fluid compressibility is assumed to be nonlinear. A Lagrangian frame of reference is used to keep track of the deformation. We perform an implicit discretization in time (backward Euler) and propose two linearization schemes for solving the nonlinear problems appearing within each time step: Newton's method and L-scheme. The linearizations are used monolithically or in combination with a splitting algorithm. The resulting schemes can be applied for any spatial discretization. The convergences of all schemes are shown analytically for cases under small deformation. Illustrative numerical examples are presented to confirm the applicability of the schemes, in particular, for large deformation.

Florin Adrian Radu, Kundan Kumar, Jan Martin Nordbotten et al.

In this work we present a mass conservative numerical scheme for two-phase flow in porous media. The model for flow consists on two fully coupled, non-linear equations: a degenerate parabolic equation and an elliptic equation. The proposed numerical scheme is based on backward Euler for the temporal discretization and mixed finite element method (MFEM) for the discretization in space. Continuous, semi-discrete (continuous in space) and fully discrete variational formulations are set up and the existence and uniqueness of solutions is discussed. Error estimates are presented to prove the convergence of the scheme. The non-linear systems within each time step are solved by a robust linearization method. This iterative method does not involve any regularization step. The convergence of the linearization scheme is rigorously proved under the assumption of a Lipschitz continuous saturation. The case of a Hölder continuous saturation is also discussed, a rigorous convergence proof being given for Richards' equation. Numerical results are presented to sustain the theoretical findings.

Daniel Førland Holmen, Jan Martin Nordbotten, Jon Eivind Vatne

We consider the simplicial de Rham complex and the Čech-de Rham complex, two bigraded Hilbert complexes whose Hodge-Laplace problems govern spatially coupled problems in mixed dimension and homogeneous dimension, respectively. The former complex can be realized as a subcomplex of the latter. In this paper, we quantify how close these complexes are to each other by constructing bounded cochain complexes between them, and thus we quantify how close a mixed-dimensional formulation of a problem is to an equidimensionally coupled formulation of the same problem. From this construction, we derive a priori- and a posteriori error estimates between the associated Hodge-Laplace problems on the two complexes. These estimates represent the error which is introduced by treating a spatially coupled problem as mixed-dimensional, rather than an equidimensional problem with thin overlaps.

Juan Michael Sargado, Eirik Keilegavlen, Inga Berre et al.

Numerical simulations of brittle fracture using phase-field approaches often employ a discrete approximation framework that applies the same order of interpolation for the displacement and phase-field variables. Most common is to use linear finite elements to discretize the linear momentum and phase-field equations. However the use of $P_1$ Lagrange shape functions to model the phase-field is not optimal, since the latter develops cusps for fully developed cracks that in turn occur at locations correspoding to Gauss points of the associated FE model for the mechanics. Such feature is challenging to reproduce accurately with low order elements, and consequently element sizes must be made very small relative to the phase-field regularization parameter in order to achieve convergence of results with respect to the mesh. In this paper, we combine the standard $P_1$ FE discretization of stress equilibrium with a cell-centered finite volume approximation of the phase-field evolution equation based on the two-point flux approximation that is constructed on the same simplex mesh. Compared to a pure FE formulation utilizing linear elements, the proposed framework results in looser restrictions on mesh refinement with respect to the phase-field length scale. Furthermore, initialization of the history field is straightforward and accomplished through a local procedure. The ability to employ a coarser mesh relative to the traditional implementation is shown for several numerical examples, demonstrating savings in computational cost on the order of 50 to 80 percent for the studied cases.

Jan Martin Nordbotten

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.