A locally conservative variational multiscale method for the simulation of porous media flow with multiscale source terms

We present a variational multiscale mixed finite element method for the solution of Darcy flow in porous media, in which both the permeability field and the source term display a multiscale character. The formulation is based on a multiscale split of the solution into coarse and subgrid scales. This decomposition is invoked in a variational setting that leads to a rigorous definition of a (global) coarse problem and a set of (local) subgrid problems. One of the key issues for the success of the method is the proper definition of the boundary conditions for the localization of the subgrid problems. We identify a weak compatibility condition that allows for subgrid communication across element interfaces, a feature that turns out to be essential for obtaining high-quality solutions. We also remove the singularities due to concentrated sources from the coarse-scale problem by introducing additional multiscale basis functions, based on a decomposition of fine-scale source terms into coarse and deviatoric components. The method is locally conservative and employs a low-order approximation of pressure and velocity at both scales. We illustrate the performance of the method on several synthetic cases and conclude that the method is able to capture the global and local flow patterns accurately.     

A coupling of mixed and continuous Galerkin finite element methods for poroelasticity II: the discrete-in-time case

In this paper, we formulate a finite element procedure for approximating the coupled fluid and mechanics in Biot’s consolidation model of poroelasticity. Here, we approximate the pressure by a mixed finite element method and the displacements by a Galerkin method. Theoretical convergence error estimates are derived in a discrete-in-time setting. Of particular interest is the case when the lowest-order Raviart–Thomas approximating space or cell-centered finite differences are used in the mixed formulation and continuous piecewise linear approximations are used for displacements. This approach appears to be the one most frequently applied to existing reservoir engineering simulators.     

An unconditionally stable explicit and precise multiple timescale finite element modeling scheme for the fully coupled hydro-mechanical analysis of saturated poroelastic media

An unconditionally stable, fully explicit and highly precise multiple timescale finite element modeling scheme is described for a fully coupled hydro-mechanical (FCHM) analysis of saturated poroelastic media. The finite element method (FEM) is used for the discretization of the FCHM differential equation in the space domain. Direct integration is performed based on the precise time step integration method (PTSIM) for the time derivatives. Two configurations for the proposed scheme are constructed (abbreviated as PTSIM-f1 and -f2, respectively). The stability and convergence of the PTSIM-f1 and -f2 are proved using a matrix-based spectral analysis in the time domain. It is demonstrated that the explicit scheme proposed in this paper is unconditionally stable and independent of the time-step size. The algorithmic error estimation results indicate that the numerical modeling performed using PTSIM-f1 and -f2 in the time domain match the computer precision. Theoretically, the algorithmic error is caused by only the mesh discretization. Therefore, the proposed modeling scheme is a semi-analytical scheme. The applicability and accuracy of the proposed scheme are examined using sample calculations. By comparing with the analytical solutions, it is indicated that the modeling results have significant advantages over the standard FEM in terms of precision and computational efficiency for large timescales.     

A coupling of mixed and continuous Galerkin finite element methods for poroelasticity I: the continuous in time case

In this paper, we formulate a finite element procedure for approximating the coupled fluid and mechanics in Biot’s consolidation model of poroelasticity. Here, we approximate the pressure by a mixed finite element method and the displacements by a Galerkin method. Theoretical convergence error estimates are derived in a continuous in-time setting for a strictly positive constrained specific storage coefficient. Of particular interest is the case when the lowest-order Raviart–Thomas approximating space or cell-centered finite differences are used in the mixed formulation, and continuous piecewise linear approximations are used for displacements. This approach appears to be the one most frequently applied to existing reservoir engineering simulators.     

A fracture mapping and extended finite element scheme for coupled deformation and fluid flow in fractured porous media

This paper presents a fracture mapping (FM) approach combined with the extended finite element method (XFEM) to simulate coupled deformation and fluid flow in fractured porous media. Specifically, the method accurately represents the impact of discrete fractures on flow and deformation, although the individual fractures are not part of the finite element mesh. A key feature of FM‐XFEM is its ability to model discontinuities in the domain independently of the computational mesh. The proposed FM approach is a continuum‐based approach that is used to model the flow interaction between the porous matrix and existing fractures via a transfer function. Fracture geometry is defined using the level set method. Therefore, in contrast to the discrete fracture flow model, the fracture representation is not meshed along with the computational domain. Consequently, the method is able to determine the influence of fractures on fluid flow within a fractured domain without the complexity of meshing the fractures within the domain. The XFEM component of the scheme addresses the discontinuous displacement field within elements that are intersected by existing fractures. In XFEM, enrichment functions are added to the standard finite element approximation to adequately resolve discontinuous fields within the simulation domain. Numerical tests illustrate the ability of the method to adequately describe the displacement and fluid pressure fields within a fractured domain at significantly less computational expense than explicitly resolving the fracture within the finite element mesh. Copyright © 2013 John Wiley & Sons, Ltd.     

Hybrid time integration and coupled solution methods for nonlinear finite element analysis of partially saturated deformable porous media at small strain

The goal of the paper is to determine the most efficient, yet accurate and stable, finite element nonlinear solution method for analysis of partially saturated deformable porous media at small strain. This involves a comparison between fully implicit, semi‐implicit, and explicit time integration schemes, with monolithically coupled and staggered‐coupled nonlinear solution methods and the hybrid combination thereof. The pore air pressure p_a is assumed atmospheric, that is, p_a=0 at reference pressure. The solid skeleton is assumed to be pressure‐sensitive nonlinear isotropic elastic. Coupled partially saturated ‘consolidation’ in the presence of surface infiltration and traction is simulated for a simple one‐dimensional uniaxial strain example and a more complicated plane strain slope example with gravity loading. Three mixed plane strain quadrilateral elements are considered: (i) Q4P4; (ii) stabilized Q4P4S; and (iii) Q9P4; “Q” refers to the number of solid skeleton displacement nodes, and “P” refers to the number of pore fluid pressure nodes. The verification of the implementation against an analytical solution for partially saturated pore water flow (no solid skeleton deformation) and comparison between the three time integration schemes (fully implicit, semi‐implicit, and explicit) are presented. It is observed that one of the staggered‐coupled semi‐implicit schemes (SIS(b)), combined with the fully implicit monolithically coupled scheme to resolve sharp transients, is the most efficient computationally. Copyright © 2015 John Wiley & Sons, Ltd.     

A micromechanical finite element model for linear and damage‐coupled viscoelastic behaviour of asphalt mixture

This study presents a finite element (FE) micromechanical modelling approach for the simulation of linear and damage‐coupled viscoelastic behaviour of asphalt mixture. Asphalt mixture is a composite material of graded aggregates bound with mastic (asphalt and fine aggregates). The microstructural model of asphalt mixture incorporates an equivalent lattice network structure whereby intergranular load transfer is simulated through an effective asphalt mastic zone. The finite element model integrates the ABAQUS user material subroutine with continuum elements for the effective asphalt mastic and rigid body elements for each aggregate. A unified approach is proposed using Schapery non‐linear viscoelastic model for the rate‐independent and rate‐dependent damage behaviour. A finite element incremental algorithm with a recursive relationship for three‐dimensional (3D) linear and damage‐coupled viscoelastic behaviour is developed. This algorithm is used in a 3D user‐defined material model for the asphalt mastic to predict global linear and damage‐coupled viscoelastic behaviour of asphalt mixture. For linear viscoelastic study, the creep stiffnesses of mastic and asphalt mixture at different temperatures are measured in laboratory. A regression‐fitting method is employed to calibrate generalized Maxwell models with Prony series and generate master stiffness curves for mastic and asphalt mixture. A computational model is developed with image analysis of sectioned surface of a test specimen. The viscoelastic prediction of mixture creep stiffness with the calibrated mastic material parameters is compared with mixture master stiffness curve over a reduced time period. In regard to damage‐coupled viscoelastic behaviour, cyclic loading responses of linear and rate‐independent damage‐coupled viscoelastic materials are compared. Effects of particular microstructure parameters on the rate‐independent damage‐coupled viscoelastic behaviour are also investigated with finite element simulations of asphalt numerical samples. Further study describes loading rate effects on the asphalt viscoelastic properties and rate‐dependent damage behaviour. Copyright © 2006 John Wiley & Sons, Ltd.     

A coupled discrete element model for the simulation of soil and water flow through an orifice

Soil erosion around defective underground pipes can cause ground collapses and sinkholes in urban areas. Most of these soil erosion events are caused by fluidization of the surrounding soil with subsequent washing into defective sewer pipes. In this study, this soil erosion process is simplified as the gradual washout of sand particles mixed with water through an orifice. The discrete element method is used to simulate the large deformation behavior of the sand particles, and the Darcy fluid model is coupled with this approach to simulate fluid flow through porous sand media. A coupled 3D discrete element model is developed and implemented based on this scheme. To simulate previous experiments using this coupled model considering the current computing capacity, we incorporated a ‘supply layer’ to study the continuous erosion process. The coupled model can predict the erosion flow rates of sand and water and the shape of erosion void. Thus, the model can be used as an effective and efficient tool to investigate the soil erosion process around defective pipes. Copyright © 2017 John Wiley & Sons, Ltd.     

A non‐linear coupled finite element–boundary element model for the prediction of vibrations due to vibratory and impact pile driving

This paper presents a non‐linear coupled finite element–boundary element approach for the prediction of free field vibrations due to vibratory and impact pile driving. Both the non‐linear constitutive behavior of the soil in the vicinity of the pile and the dynamic interaction between the pile and the soil are accounted for. A subdomain approach is used, defining a generalized structure consisting of the pile and a bounded region of soil around the pile, and an unbounded exterior linear soil domain. The soil around the pile may exhibit non‐linear constitutive behavior and is modelled with a time‐domain finite element method. The dynamic stiffness matrix of the exterior unbounded soil domain is calculated using a boundary element formulation in the frequency domain based on a limited number of modes defined on the interface between the generalized structure and the unbounded soil. The soil–structure interaction forces are evaluated as a convolution of the displacement history and the soil flexibility matrices, which are obtained by an inverse Fourier transformation from the frequency to the time domain. This results in a hybrid frequency–time domain formulation of the non‐linear dynamic soil–structure interaction problem, which is solved in the time domain using Newmark's time integration method; the interaction force time history is evaluated using the θ‐scheme in order to obtain stable solutions. The proposed hybrid formulation is validated for linear problems of vibratory and impact pile driving, showing very good agreement with the results obtained with a frequency‐domain solution. Linear predictions, however, overestimate the free field peak particle velocities as observed in reported field experiments during vibratory and impact pile driving at comparable levels of the transferred energy. This is mainly due to energy dissipation related to plastic deformations in the soil around the pile. Ground vibrations due to vibratory and impact pile driving are, therefore, also computed with a non‐linear model where the soil is modelled as an isotropic elastic, perfectly plastic solid, which yields according to the Drucker–Prager failure criterion. This results in lower predicted free field vibrations with respect to linear predictions, which are also in much better agreement with experimental results recorded during vibratory and impact pile driving. Copyright © 2008 John Wiley & Sons, Ltd.