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

Multiscale Method for Elastic Wave Propagation in the Heterogeneous, Anisotropic Media

Seismic wave simulation in realistic Earth media with full wavefield methods is a fundamental task in geophysical studies. Conventional approaches such as the finite-difference method and the finite-element method solve the wave equation in geological models represented with discrete grids and elements. When the Earth model includes complex heterogeneities at multiple spatial scales, the simulation requires fine discretization and therefore a system with many degrees of freedom, which often exceeds current computational abilities. In this dissertation, I address this problem by proposing new multiscale methods for simulating elastic wave propagation based on previously developed algorithms for solving the elliptic partial differential equations and the acoustic wave equation. The fundamental motivation for developing the multiscale method is that it can solve the wave equation on a coarsely discretized mesh by incorporating the effects of fine-scale medium properties using so-called multiscale basis functions. This can greatly reduce computation time and degrees of freedom compared with conventional methods. I first derive a numerical homogenization method for arbitrarily heterogeneous, anisotropic media that utilizes the multiscale basis functions determined from a local linear elasticity equation to compute effective, anisotropic properties, and these equivalent elastic medium parameters can be used directly in existing elastic modeling algorithms. Then I extend the approach by constructing multiple basis functions using two types of appropriately defined local spectral linear elasticity problems. Given the eigenfunctions determined from local spectral problems, I develop a generalized multiscale finite-element method (GMsFEM) for elastic wave propagation in heterogeneous, anisotropic media in both continuous Galerkin (CG) and discontinuous Galerkin (DG) formulations. The advantage of the multiscale basis functions is they are model-dependent, unlike the predefined polynomial basis functions applied in conventional finite-element methods. For this reason, the GMsFEM can effectively capture the influence of fine-scale variation of the media. I present results for several numerical experiments to verify the effectiveness of both the numerical homogenization method and GMsFEM. These tests show that the effectiveness of the multiscale method relies on the appropriate choice of boundary conditions that are applied for the local problem in numerical homogenization method and on the selection of basis functions from a large set of eigenfunctions contained in local spectral problems in GMsFEM. I develop methods for solving both these problems, and the results confirm that the multiscale method can be powerful tool for providing accurate full wavefield solutions in heterogeneous, anisotropic media, yet with reduced computation time and degrees of freedom compared with conventional full wavefield modeling methods. Specially, I applied the DG-GMsFEM to the Marmousi-2 elastic model, and find that DG-GMsFEM can greatly reduce the computation time compared with continuous Galerkin (CG) FEM

A first order hyperbolic framework for large strain computational solid dynamics: An upwind cell centred Total Lagrangian scheme

This paper builds on recent work developed by the authors for the numerical analysis of large strain solid dynamics, by introducing an upwind cell centred hexahedral Finite Volume framework implemented within the open source code OpenFOAM [http://www.openfoam.com/http://www.openfoam.com/]. In Lee, Gil and Bonet [1], a first order hyperbolic system of conservation laws was introduced in terms of the linear momentum and the deformation gradient tensor of the system, leading to excellent behaviour in two dimensional bending dominated nearly incompressible scenarios. The main aim of this paper is the extension of this algorithm into three dimensions, its tailor-made implementation into OpenFOAM and the enhancement of the formulation with three key novelties. First, the introduction of two different strategies in order to ensure the satisfaction of the underlying involutions of the system, that is, that the deformation gradient tensor must be curl-free throughout the deformation process. Second, the use of a discrete angular momentum projection algorithm and a monolithic Total Variation Diminishing Runge-Kutta time integrator combined in order to guarantee the conservation of angular momentum. Third, and for comparison purposes, an adapted Total Lagrangian version of the Hyperelastic-GLACE nodal scheme of Kluth and Despr´es [2] is presented. A series of challenging numerical examples are examined in order to assess the robustness and accuracy of the proposed algorithm, benchmarking it against an ample spectrum of alternative numerical strategies developed by the authors in recent publications

Computational Engineering

The focus of this Computational Engineering Workshop was on the mathematical foundation of state-of-the-art and emerging finite element methods in engineering analysis. The 52 participants included mathematicians and engineers with shared interest on discontinuous Galerkin or Petrov-Galerkin methods and other generalized nonconforming or mixed finite element methods

Poroelasticity problem: numerical difficulties and efficient multigrid solution

This work contains some of the more relevant results obtained by the author regarding the numerical solution of the Biot’s consolidation problem. The emphasis here is on the stable discretization and the highly efficient solution of the resulting algebraic system of equations, which is of saddle point type. On the one hand, a stabilized linear finite element scheme providing oscillation-free solutions for this model is proposed and theoretically analyzed. On the other hand, a monolithic multigrid method is considered for the solution of the resulting system of equations after discretization by using the stabilized scheme. Since this system is of saddle point type, special smoothers of “Vanka”-type have to be considered. This multigrid method is designed with the help of an special local Fourier analysis that takes into account the specific characteristics of the considered block-relaxations. Results from this analysis are presented and compared with those experimentally computed