Assessment of acceleration modelling for fluid‐filled porous media subjected to dynamic loading

The purpose of this paper is to examine the importance of different possible simplifying approximations when performing numerical simulations of fluid‐filled porous media subjected to dynamic loading. In particular, the relative importance of the various acceleration terms for both the solid and the fluid, especially the convective contribution, is assessed. The porous medium is modelled as a binary mixture of a solid phase, in the sense of a porous skeleton, and a fluid phase that represents both liquid and air in the pores. The solid particles are assumed to be intrinsically incompressible, whereas the fluid is assigned a finite intrinsic compressibility. Finite element (FE) simulations are carried out while assuming material properties and loading conditions representative for a road structure. The results show that, for the range of the material data used in the simulations, omitting the relative acceleration gives differences in the solution of the seepage velocity field, whereas omitting only the convective term does not lead to significant differences. Copyright © 2007 John Wiley & Sons, Ltd.     

A model for moisture and heat flow in fractured unsaturated porous media

The present study investigates propagation of a cohesive crack in non‐isothermal unsaturated porous medium under mode I conditions. Basic points of skeleton deformation, moisture, and heat transfer for unsaturated porous medium are presented. Boundary conditions on the crack surface that consist of mechanical interaction of the crack and the porous medium, water, and heat flows through the crack are taken into consideration. For spatial discretization, the extended finite element method is used. This method uses enriched shape functions in addition to ordinary shape functions for approximation of displacement, pressure, and temperature fields. The Heaviside step function and the distance function are exploited as enrichment functions for representing the crack surfaces displacement and the discontinuous vertical gradients of the pressure and temperature fields along the crack, respectively. For temporal discretization, backward finite difference scheme is applied. Problems solved from the literature show the validity of the model as well as the dependency of structural response on the material properties and loading. Copyright © 2016 John Wiley & Sons, Ltd.     

Fluctuation responses of saturated porous media subjected to cyclic thermal loading

Normalized, coupled governing equations for one-dimensional thermal consolidation problems are established. The non-dimensional coefficients of thermal consolidation and thermal diffusivity are defined accordingly. An analytical solution is deduced by using the Laplace transform and the Gauss–Legendre method of Laplace transform inversion. The responses of saturated porous media subjected to cyclic thermal loading are studied. The evolution of temperature, pore pressure and displacement from instantaneous state to quasi-steady state, with elapsed time, are analysed. The characteristics of cyclic fluctuation and the attenuation of the field variables with increased depth are also analysed. The influences of the permeability of media on thermal responses are discussed.     

Numerical study of two‐phase fluid distributions in fractured porous media

Two‐phase fluid distributions in fractured porous media were studied using a single‐component multiphase (SCMP) lattice Boltzmann method (LBM), which was selected among three commonly used numerical approaches through a comparison against the available results of micro X‐ray computed tomography. The influence of the initial configuration and the periodic boundary conditions in the SCMP LBM for the fluid distribution analysis were investigated as well. It was revealed that regular porous media are sensitive to the initial distribution, whereas irregular porous media are insensitive. Moreover, to eliminate the influence of boundaries, the model's buffer size of an SCMP LBM simulation was suggested to be taken as approximately 12.5 times the average particle size. Then, the two‐phase fluid distribution of a porous medium was numerically studied using the SCMP LBM. Both detailed distribution patterns and macroscopic morphology parameters were reasonably well captured. Finally, the two‐phase fluid distributions in a fractured porous media were investigated. The influence of the degree of saturation, fracture length, and fracture width on the fluid distributions and migration was explored. Copyright © 2015 John Wiley & Sons, Ltd.     

A fully coupled element‐free Galerkin model for hydro‐mechanical analysis of advancement of fluid‐driven fractures in porous media

Hydraulic fracturing (HF) of underground formations has widely been used in different fields of engineering. Despite the technological advances in techniques of in situ HF, the industry uses semi‐analytical tools to design HF treatment. This is due to the complex interaction among various mechanisms involved in this process, so that for thorough simulations of HF operations a fully coupled numerical model is required. In this study, using element‐free Galerkin (EFG) mesh‐less method, a new formulation for numerical modeling of hydraulic fracture propagation in porous media is developed. This numerical approach, which is based on the simultaneous solution of equilibrium and continuity equations, considers the hydro‐mechanical coupling between the crack and its surrounding porous medium. Therefore, the developed EFG model is capable of simulating fluid leak‐off and fluid lag phenomena. To create the discrete equation system, the Galerkin technique is applied, and the essential boundary conditions are imposed via penalty method. Then, the resultant constrained integral equations are discretized in space using EFG shape functions. For temporal discretization, a fully implicit scheme is employed. The final set of algebraic equations that forms a non‐linear equation system is solved using the direct iterative procedure. Modeling of cracks is performed on the basis of linear elastic fracture mechanics, and for this purpose, the so‐called diffraction method is employed. For verification of the model, a number of problems are solved. According to the obtained results, the developed EFG computer program can successfully be applied for simulating the complex process of hydraulic fracture propagation in porous media. Copyright © 2016 John Wiley & Sons, Ltd.     

Thermo‐hydro‐mechanical modeling of damage in unsaturated porous media: Theoretical framework and numerical study of the EDZ

The damage model presented in this article (named ‘THHMD’ model) is dedicated to non‐isothermal unsaturated porous media. It is formulated by means of three independent strain state variables, which are the thermodynamic conjugates of net stress, suction and thermal stress. The damage variable is a second‐order tensor. Stress/strain relationships are derived from Helmholtz free energy, which is assumed to be the sum of damaged elastic potentials and ‘crack‐closure energies’. Damage is assumed to grow with tensile strains due to net stress, with pore shrinkage due to suction and with thermal dilatation. Specific conductivities are introduced to account for the effects of cracking on the intensification and on the orientation of liquid water and vapor flows. These conductivities depend on damage and internal length parameters. The mechanical aspects of the THHMD model are validated by comparing the results of a triaxial compression test with experimental measurements found in the literature. Parametric studies of damage are performed on three different heating problems related to nuclear waste disposals. Several types of loading and boundary conditions are investigated. The thermal damage potential is thoroughly studied. The THHMD model is expected to be a useful tool in the assessment of the Excavation Damaged Zone, especially in the vicinity of nuclear waste repositories. Copyright © 2011 John Wiley & Sons, Ltd.     

A numerical model for nonlinear large deformation dynamic analysis of unsaturated porous media including hydraulic hysteresis

A numerical model based on the theory of mixtures is proposed for the nonlinear dynamic analysis of flow and deformation in unsaturated porous media. Starting from the conservation laws, the governing differential equations and the finite element incremental approximations suitable for nonlinear large deformation static and dynamic analyses are derived within the updated Lagrangian framework. The coupling between solid and fluid phases is enforced according to the effective stress principle taking suction dependency of the effective stress parameter into account. The effect of hydraulic hysteresis on the effective stress parameter and soil water characteristic curve is also taken into account. The application of the approach is demonstrated through numerical analyses of several fundamental nonlinear problems and the results are compared to the relevant analytical solutions. The effects of suction, large deformations and hydraulic hysteresis on static and dynamic response of unsaturated soils are particularly emphasized.     

Hydro‐mechanical modeling of cohesive crack propagation in multiphase porous media using the extended finite element method

In this paper, a numerical model is developed for the fully coupled hydro‐mechanical analysis of deformable, progressively fracturing porous media interacting with the flow of two immiscible, compressible wetting and non‐wetting pore fluids, in which the coupling between various processes is taken into account. The governing equations involving the coupled solid skeleton deformation and two‐phase fluid flow in partially saturated porous media including cohesive cracks are derived within the framework of the generalized Biot theory. The fluid flow within the crack is simulated using the Darcy law in which the permeability variation with porosity because of the cracking of the solid skeleton is accounted. The cohesive crack model is integrated into the numerical modeling by means of which the nonlinear fracture processes occurring along the fracture process zone are simulated. The solid phase displacement, the wetting phase pressure and the capillary pressure are taken as the primary variables of the three‐phase formulation. The other variables are incorporated into the model via the experimentally determined functions, which specify the relationship between the hydraulic properties of the fracturing porous medium, that is saturation, permeability and capillary pressure. The spatial discretization is implemented by employing the extended finite element method, and the time domain discretization is performed using the generalized Newmark scheme to derive the final system of fully coupled nonlinear equations of the hydro‐mechanical problem. It is illustrated that by allowing for the interaction between various processes, that is the solid skeleton deformation, the wetting and the non‐wetting pore fluid flow and the cohesive crack propagation, the effect of the presence of the geomechanical discontinuity can be completely captured. Copyright © 2012 John Wiley & Sons, Ltd.     

Extension of Biot theory to the problem of saturated microporous elastic media with isolated cracks or/and vugs

The purpose of this paper is to develop the macroscopic model of hydro‐mechanical coupling for the case of a porous medium containing isolated cracks or/and vugs. In the development, we apply the asymptotic expansion homogenization method. It is shown that the general structure of Biot's model is the same as in the case of homogeneous medium, but the poro‐elastic parameters are modified. Two numerical examples are presented. They concern the computations of Biot's parameters in isotropic and anisotropic cases. It can also be seen how the presence of near‐zero‐volume cracks influences Biot's parameters of the porous matrix. It can significantly affect the coupled hydro‐mechanical behaviour of damaged porous medium. Copyright © 2012 John Wiley & Sons, Ltd.     

Modelling of steady‐state fluid flow in 3D fractured isotropic porous media: application to effective permeability calculation

In this paper, 3D steady‐state fluid flow in a porous medium with a large number of intersecting fractures is derived numerically by using collocation method. Fluid flow in the matrix and fractures is described by Darcy's law and Poiseuille's law, respectively. The recent theoretical development presented a general potential solution to model the steady‐state flow in fractured porous media under a far‐field condition. This solution is a hypersingular integral equation with pressure field in the fracture surfaces as the main unknown. The numerical procedure can resolve the problem for any form of fractures and also takes into account the interactions and the intersection between fractures. Once the pressure field and then the flux field in fractures have been determined, the pressure field in the porous matrix is computed completely. The basic problem of a single fracture is investigated, and a semi‐analytical solution is presented. Using the solution obtained for a single fracture, Mori‐Tanaka and self‐consistent schemes are employed for upscaling the effective permeability of 3D fractured porous media. Copyright © 2012 John Wiley & Sons, Ltd.     

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

In this contribution an algorithm for parameter identification of geometrically linear Terzaghi–Biot‐type fluid‐saturated porous media is proposed, in which non‐uniform distributions of the state variables such as stresses, strains and fluid pore pressure are taken into account. To this end a least‐squares functional consisting of experimental data and simulated data is minimized, whereby the latter are obtained with the finite element method. This strategy allows parameter identification based on in situ experiments. In order to improve the efficiency of the minimization process, a gradient‐based optimization algorithm is applied, and therefore the corresponding sensitivity analysis for the coupled two‐phase problem is described in a systematic manner. For illustrative purpose, the performance of the algorithm is demonstrated for a slope stability problem, in which a quadratic Drucker–Prager plasticity model for the solid and a linear Darcy law for the fluid are combined. Copyright © 2001 John Wiley & Sons, Ltd.     

A three‐dimensional staggered finite element approach for random parametric modeling of thermo‐hygral coupled phenomena in porous media

The aim of this paper is to present a three‐dimensional (3D) finite element modeling of heat and mass transfer phenomena in partially saturated open porous media with random fields of material properties. Randomness leads to transfer processes within the porous medium that naturally need a full 3D modeling for any quantitative assessment of these processes. Nevertheless, the counterpart of 3D modeling is a significant increase in computations cost. Therefore, a staggered solution strategy is adopted which permits to solve the equations sequentially. This appropriate partitioning reduces the size of the discretized problem to be solved at each time step. It is based on a specific iterative algorithm to account for the interaction between all the transfer processes. Accordingly, a suitable linearization of mass convective boundary conditions, consistent with the staggered algorithm, is also derived. After some validation tests, the 3D numerical model is used for studying the drying process of a cementitious material with regard to its intrinsic permeability randomness. Copyright © 2011 John Wiley & Sons, Ltd.     

Numerical modeling of multiphase fluid flow in deforming porous media: A comparison between two- and three-phase models for seismic analysis of earth and rockfill dams

In this paper, a fully coupled numerical model is presented for the finite element analysis of the deforming porous medium interacting with the flow of two immiscible compressible wetting and non-wetting pore fluids. The governing equations involving coupled fluid flow and deformation processes in unsaturated soils are derived within the framework of the generalized Biot theory. The displacements of the solid phase, the pressure of the wetting phase and the capillary pressure are taken as the primary unknowns of the present formulation. The other variables are incorporated into the model using the experimentally determined functions that define the relationship between the hydraulic properties of the porous medium, i.e. saturation, relative permeability and capillary pressure. It is worth mentioning that the imposition of various boundary conditions is feasible notwithstanding the choice of the primary variables. The modified Pastor–Zienkiewicz generalized constitutive model is introduced into the mathematical formulation to simulate the mechanical behavior of the unsaturated soil. The accuracy of the proposed mathematical model for analyzing coupled fluid flows in porous media is verified by the resolution of several numerical examples for which previous solutions are known. Finally, the performance of the computational algorithm in modeling of large-scale porous media problems including the large elasto-plastic deformations is demonstrated through the fully coupled analysis of the failure of two earth and rockfill dams. Furthermore, the three-phase model is compared to its simplified one which simulates the unsaturated porous medium as a two-phase one with static air phase. The paper illustrates the shortcomings of the commonly used simplified approach in the context of seismic analysis of two earth and rockfill dams. It is shown that accounting the pore air as an independent phase significantly influences the unsaturated soil behavior.     

A new computational strategy for solving two‐phase flow in strongly heterogeneous poroelastic media of evolving scales

We develop a new computational methodology for solving two‐phase flow in highly heterogeneous porous media incorporating geomechanical coupling subject to uncertainty in the poromechanical parameters. Within the framework of a staggered‐in‐time coupling algorithm, the numerical method proposed herein relies on a Petrov–Galerkin postprocessing approach projected on the Raviart–Thomas space to compute the Darcy velocity of the mixture in conjunction with a locally conservative higher order finite volume discretization of the nonlinear transport equation for the saturation and an operator splitting procedure based on the difference in the time‐scales of transport and geomechanics to compute the effects of transient porosity upon saturation. Notable features of the numerical modeling proposed herein are the local conservation properties inherited by the discrete fluxes that are crucial to correctly capture the fingering patterns arising from the interaction between heterogeneity and nonlinear viscous coupling. Water flooding in a poroelastic formation subject to an overburden is simulated with the geology characterized by multiscale self‐similar permeability and Young modulus random fields with power‐law covariance structure. Statistical moments of the poromechanical unknowns are computed within the framework of a high‐resolution Monte Carlo method. Numerical results illustrate the necessity of adopting locally conservative schemes to obtain reliable predictions of secondary recovery and finger growth in strongly heterogeneous deformable reservoirs. Copyright © 2011 John Wiley & Sons, Ltd.