Abstract:In this paper, we propose a numerical method for computing solutions to Biot's fully dynamic model of incompressible saturated porous media [Biot;1956]. Our spatial discretization scheme is based on the three-field formulation (u-w-p) and the coupling of a lowest order Raviart-Thomas mixed element [Raviart,Thomas;1977] for fluid variable fields (w, p ) and a nodal Galerkin finite element for skeleton variable field (u). These mixed spaces are constructed based on the natural topology of the variables; hence, are physically compatible and able to exactly model the kind of continuity which is expected. The method automatically satisfies the well known LBB (inf-sup) stability condition and avoids locking that usually occurs in the numerical computations in the incompressible limit and very low hydraulic conductivity. In contrast to the majority of approaches, our three-field formulation can fully capture dynamic behavior of porous media even in high frequency loading phenomena with considerable fluid acceleration such as liquefaction and biomechanics of porous tissues under rapid external loading. Moreover, we address the importance of consistent initial conditions for poroelasticity equations with the incompressibility constraint, which represent a system of differential algebraic equations. The energy balance equation is derived for the full porous medium and used to assess the stability and accuracy of our time integration. To highlight the capabilities of our method, a variety of numerical studies are provided including verification with analytical and boundary element solutions, wave propagation analyses, hydraulic conductivity effects on damping and frequency content, energy balance analyses, mass lumping considerations, effects of mesh pattern and size, and stability analyses. We also explain some discrepancies commonly found in dynamic poroelasticity results in the literature.

A finite element algorithm for parameter identification of material models for fluid saturated porous media

A methodology for parameter identification of poro-viscoelastic models of soft soil

Mixed finite element method for saturated poroelastoplastic media at large strains

Parameter Identification for Two-Dimensional Orthotropic Material Bodies by the Boundary Element Method

A topology-motivated mixed finite element method for dynamic response of porous media

A finite element analysis of the influence of effective stress principle on coupled problems of flow and deformation in porous media

Parameter identification for coupled finite element–boundary element models

Parameter identification by optimization method for a pollution problem in porous media

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

Identifiability of material parameters in solid mechanics

Exact Solutions for One‐dimensional Transient Response of Fluid‐saturated Porous Media

Identification of elastic orthotropic material parameters by the scaled boundary finite element method

Optimization of the porous material described by the Biot model

Parameter Identification of the Fractional-Order Mathematical Model for Convective Mass Transfer in a Porous Medium

Nonlinear material identification of heterogeneous isogeometric Kirchhoff-Love shells

A hybrid algorithm of particle swarm optimization and finite element method to identify local mesoscopic damage of concrete-like materials

Identification of material parameters in low-data limit: application to gradient-enhanced continua

Dynamic Analysis of Fluid Saturated Porous Viscoelastic Media by Explicit Finite Element Method

Computational identification of adsorption and desorption parameters for pore scale transport in periodic porous media

A gradient based numerical algorithm to solve inverse dynamic viscoelastic problems of multi-variable identification

A mixed nonlocal finite element model for thermo‐poro‐elasto‐plastic simulation of porous media with multiphase fluid flow