Numerical solution for consolidation and desiccation of soft soils

The consolidation and desiccation behaviour of soft soils can be described by two time‐dependent non‐linear partial differential equations using the finite strain theory. Analytical solutions do not exist for these governing equations. In this paper, we develop efficient numerical methods and software for finding the numerical solutions. We introduce a semi‐implicit time integration scheme, and show numerically that our method converges. In addition, the numerical solution matches well with the experimental result. A boundary refinement method is also developed to improve the convergence and stability for the case of Neumann type boundary conditions. Interface governing equations are derived to maintain the continuity of consolidation and desiccation processes. This is useful because the soil column can undergo desiccation on top and consolidation on the bottom simultaneously. The numerical algorithms has been implemented into a computer program and the results have been verified with centrifuge test results conducted in our laboratory. Copyright © 2001 John Wiley & Sons, Ltd.     

Numerical modelling of fluid flow in microscopic images of granular materials

A program for the simulation of two‐dimensional (2‐D) fluid flow at the microstructural level of a saturated anisotropic granular medium is presented. The program provides a numerical solution to the complete set of Navier–Stokes equations without a priori assumptions on the viscous or convection components. This is especially suited for the simulation of the flow of fluids with different density and viscosity values and for a wide range of granular material porosity. The analytical solution for fluid flow in a simple microstructure of porous medium is used to verify the computer program. Subsequently, the flow field is computed within microscopic images of granular material that differ in porosity, particle size and particle shape. The computed flow fields are shown to follow certain paths depending on air void size and connectivity. The permeability tensor coefficients are derived from the flow fields, and their values are shown to compare well with laboratory experimental data on glass beads, Ottawa sand and silica sands. The directional distribution of permeability is expressed in a functional form and its anisotropy is quantified. Permeability anisotropy is found to be more pronounced in the silica sand medium that consists of elongated particles. Copyright © 2001 John Wiley & Sons, Ltd.     

Three‐dimensional transient thermo‐hydro‐mechanical fundamental solutions of unsaturated soils

The governing differential equations of unsaturated soils considering the thermo‐poro‐mechanical behaviour consist of equilibrium, moisture air and heat transfer equations. In this paper at first, following some necessary simplifications, the thermal three‐dimensional fundamental solution for an unsaturated deformable porous medium with linear elastic behaviour in Laplace transform domain is presented. Subsequently, the closed‐form time domain fundamental solutions are derived by analytical inversion of the Laplace transform domain solutions. Then a set of numerical results are presented, which demonstrate the accuracies and some salient features of the derived analytical transient fundamental solutions. Finally, the closed‐form time domain fundamental solution will be verified mathematically by comparison with the previously introduced corresponding fundamental solution. Copyright © 2009 John Wiley & Sons, Ltd.     

Steady‐state analytical solutions of flow and deformation coupling due to a point sink within a finite fluid‐saturated poroelastic layer

An exact steady‐state closed‐form solution is presented for coupled flow and deformation of an axisymmetric isotropic homogeneous fluid‐saturated poroelastic layer with a finite radius due to a point sink. The hydromechanical behavior of the poroelastic layer is governed by Biot's consolidation theory. Boundary conditions on the lateral surface are specifically chosen to match the appropriate finite Hankel transforms and simplify the transforms of the governing equations. Ordinary differential equations in the transformed domain are solved, and then the analytical solutions in the physical space for the pore pressure and the displacements are finally obtained by using finite Hankel inversions. The analytical solutions at some special locations such as the top and bottom surfaces, lateral surface, and the symmetrical axis are given and analyzed. And a case study for the consolidation of a water‐saturated soft clay layer due to pumping is conducted. The analytical solution is verified against the finite element solution. Meanwhile, an analysis of coupled hydromechanical behavior is carried out herein. The presented analytical solution is an exact solution to the practical poroelastic problem within an axisymmetric finite layer. It can provide us a better understanding of the poroelastic behavior of the finite layer due to fluid extraction. Besides, it can be applied to calibrate numerical schemes of axisymmetric poroelasticity within finite domains. Copyright © 2017 John Wiley & Sons, Ltd.     

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.     

A constitutive model for mechanical and chemo‐mechanical compaction in sedimentary basins and finite element analysis

Compaction and associated fluid flow are fundamental processes in sedimentary basin deformation. Purely mechanical compaction originates mainly from pore fluid expulsion and rearrangement of solid particles during burial, while chemo‐mechanical compaction results from Intergranular Pressure‐Solution (IPS) and represents a major mechanism of deformation in sedimentary basins during diagenesis. The aim of the present contribution is to provide a comprehensive 3D framework for constitutive and numerical modeling of purely mechanical and chemo‐mechanical compaction in sedimentary basins. Extending the concepts that have been previously proposed for the modeling of purely mechanical compaction in finite poroplasticity, deformation by IPS is addressed herein by means of additional viscoplastic terms in the state equations of the porous material. The finite element model integrates the poroplastic and poroviscoplastic components of deformation at large strains. The corresponding implementation allows for numerical simulation of sediments accretion/erosion periods by progressive activation/deactivation of the gravity forces within a fictitious closed material system. Validation of the numerical approach is assessed by means of comparison with closed‐form solutions derived in the context of a simplified compaction model. The last part of the paper presents the results of numerical basin simulation performed in one dimensional setting, demonstrating the ability of the modeling to capture the main features in elastoplastic and viscoplastic compaction. Copyright © 2016 John Wiley & Sons, Ltd.     

Coupled formulation and algorithms for the simulation of non‐planar three‐dimensional hydraulic fractures using the generalized finite element method

This paper presents a coupled hydro‐mechanical formulation for the simulation of non‐planar three‐dimensional hydraulic fractures. Deformation in the rock is modeled using linear elasticity, and the lubrication theory is adopted for the fluid flow in the fracture. The governing equations of the fluid flow and elasticity and the subsequent discretization are fully coupled. A Generalized/eXtended Finite Element Method (G/XFEM) is adopted for the discretization of the coupled system of equations. A Newton–Raphson method is used to solve the resulting system of nonlinear equations. A discretization strategy for the fluid flow problem on non‐planar three‐dimensional surfaces and a computationally efficient strategy for handling time integration combined with mesh adaptivity are also presented. Several three‐dimensional numerical verification examples are solved. The examples illustrate the generality and accuracy of the proposed coupled formulation and discretization strategies. Copyright © 2015 John Wiley & Sons, Ltd.     

Some numerical issues using element‐free Galerkin mesh‐less method for coupled hydro‐mechanical problems

A new formulation of the element‐free Galerkin (EFG) method is developed for solving coupled hydro‐mechanical problems. The numerical approach is based on solving the two governing partial differential equations of equilibrium and continuity of pore water simultaneously. Spatial variables in the weak form, i.e. displacement increment and pore water pressure increment, are discretized using the same EFG shape functions. An incremental constrained Galerkin weak form is used to create the discrete system equations and a fully implicit scheme is used for discretization in the time domain. Implementation of essential boundary conditions is based on a penalty method. Numerical stability of the developed formulation is examined in order to achieve appropriate accuracy of the EFG solution for coupled hydro‐mechanical problems. Examples are studied and compared with closed‐form or finite element method solutions to demonstrate the validity of the developed model and its capabilities. The results indicate that the EFG method is capable of handling coupled problems in saturated porous media and can predict well both the soil deformation and variation of pore water pressure over time. Some guidelines are proposed to guarantee the accuracy of the EFG solution for coupled hydro‐mechanical problems. Copyright © 2008 John Wiley & Sons, Ltd.     

Numerical modelling of multiphase immiscible flow in double‐porosity featured groundwater systems

A numerical model describing the flow of multiphase, immiscible fluids in a deformable, double‐porosity featured soil has been developed. The model is focused on the modelling of the secondary porosity features in soil, which is more relevant to groundwater contamination problems. The non‐linear saturation and relative permeabilities were expressed as functions of the capillary pressures. The governing partial differential equations in terms of soil displacement and fluid pressures were solved numerically. Galerkin's weighted‐residual finite element method was employed to obtain the spatial discretization whereas temporal discretization was achieved using a fully implicit scheme. The model was verified against established, peer‐reviewed works, and the assumption that the immiscible fluids (non‐aqueous phase liquids) will flow preferentially through the secondary porosity features in soil was validated. Copyright © 2011 John Wiley & Sons, Ltd.     

Coupled hydromechanical‐fracture simulations of nonplanar three‐dimensional hydraulic fracture propagation

This paper presents an algorithm and a fully coupled hydromechanical‐fracture formulation for the simulation of three‐dimensional nonplanar hydraulic fracture propagation. The propagation algorithm automatically estimates the magnitude of time steps such that a regularized form of Irwin's criterion is satisfied along the predicted 3‐D fracture front at every fracture propagation step. A generalized finite element method is used for the discretization of elasticity equations governing the deformation of the rock, and a finite element method is adopted for the solution of the fluid flow equation on the basis of Poiseuille's cubic law. Adaptive mesh refinement is used for discretization error control, leading to significantly fewer degrees of freedom than available nonadaptive methods. An efficient computational scheme to handle nonlinear time‐dependent problems with adaptive mesh refinement is presented. Explicit fracture surface representations are used to avoid mapping of 3‐D solutions between generalized finite element method meshes. Examples demonstrating the accuracy, robustness, and computational efficiency of the proposed formulation, regularized Irwin's criterion, and propagation algorithm are presented.     

Analytical 3D transient elastodynamic fundamental solution of unsaturated soils

Unsaturated soils are considered as porous continua, composed of porous skeleton with its pores filled by water and air. The governing partial differential equations (PDE) are derived based on the mechanics for isothermal and infinitesimal evolution of unsaturated porous media in terms of skeleton displacement vector, liquid, and gas scalar pressures. Meanwhile, isotropic linear elastic behavior and liquid retention curve are presented in terms of net stress and capillary pressure as constitutive relations. Later, an explicit 3D Laplace transform domain fundamental solution is obtained for governing PDE and then closed‐form analytical transient 3D fundamental solution is presented by means of analytical inverse Laplace transform technique. Finally, a numerical example is presented to validate the assumptions used to derive the analytical solution by comparing them with the numerically inverted ones. The transient fundamental solutions represent important features of the elastic wave propagation theory in the unsaturated soils. Copyright © 2010 John Wiley & Sons, Ltd.     

FFT‐based bounds on the permeability of complex microstructures

A numerical scheme for the computation of the permeability of complex microstructures is presented. As a darcean counterpart of the FFT‐based scheme in elasticity, the method is designed to be directly coupled with 3D imaging techniques of porous samples, without meshing or definition of an equivalent pore network. The method relies on the variational principle of Hashin and Shtrikman, which ensures a rigorous upper bound status to the estimated permeabilities and provides an energetically consistent rule for heterogeneous voxels comprising both the solid and fluid phases. Copyright © 2014 John Wiley & Sons, Ltd.     

Study on low‐strain integrity testing of pipe‐pile using the elastodynamic finite integration technique

This study focuses on the three‐dimensional (3‐D) characteristics of wave propagation in pipe‐pile using elastodynamic finite integration technique. First, a real 3‐D pile‐soil model in cylindrical coordinate system is presented. Then, the governing equations are established. With the boundary and initial conditions, the numerical solution is obtained. The accuracy and feasibility of the self‐written code are further verified via comparing with the measured data. Velocity histories at different angles of pile top and pile tip are illustrated, and the snapshots reflecting the 3‐D characteristics of wave propagation are given. It shows that the interferences of Rayleigh waves can confuse the result interpretation for pile integrity testing. The increase of hammer contact time can effectively mitigate the interferences, and the interferences of Rayleigh waves are weakest at an angle of 90° from where hammer hits. Besides, surrounding soil can partly mitigate the wave interferences. Copyright © 2012 John Wiley & Sons, Ltd.     

Responses of chemically active and naturally fractured shale under time‐dependent mechanical loading and ionic solution exposure

In this paper, the analytical dual‐porosity dual‐permeability poromechanics solution for saturated cylinders is extended to account for electrokinetic effects and material transverse isotropy, which simulate the responses of chemically active naturally fractured shale under time‐dependent mechanical loading and ionic solution exposure. The solution addresses the stresses, fracture pore pressure, matrix pore pressure, fluid fluxes, ion concentration evolution, and displacements due to the applied stress, pore pressure, and solute concentration difference between the sample and the circulation fluid. The presented solution will not only help validate numerical simulations but also assist in calibrating and interpreting laboratory results on dual‐porosity dual‐permeability shale. It is recommended that the analytical solutions of radial and axial displacements be used to match the corresponding laboratory‐recorded data to determine shale dual permeability and chemo‐electrical parameters including membrane coefficient, ions diffusion coefficients, and electro‐osmotic permeability.     

Coupled consolidation of a porous medium with a cylindrical or a spherical cavity

This paper presents a theoretical approach to analyse coupled, linear thermoporoelastic fields in a saturated porous medium under radial and spherical symmetry. The governing equations account for compressibility and thermal expansion of constituents, heat sink due to thermal dilatation of water and thermal expansion of the medium, and thermodynamically coupled heat–water flow. It has been reported in the literature that thermodynamically coupled heat–water flows known as thermo-osmosis and thermal filtration have the potential to significantly alter the flow fields in clay-rich barriers in the near field of a underground waste containment scheme. This study presents a mathematical model and examines the effects of thermo-osmosis and thermal-filtration on coupled consolidation fields in a porous medium with a cavity. Analytical solutions of the governing equations are presented in the Laplace transform space. A numerical inversion scheme is used to obtain the time-domain solutions for a cylindrical cavity in a homogeneous or a non-homogeneous medium. A closed form time-domain solution is presented for a spherical cavity in a homogeneous medium. Selected numerical solutions for homogeneous and non-homogeneous media show a significant increase in pore pressure and displacements due to the presence of thermodynamically coupled flows and a negligible influence on temperature. © 1998 John Wiley & Sons, Ltd.     

Finite element modelling of three-phase flow in deforming saturated oil reservoirs

A numerical simulation is presented for three-dimensional three-phase fluid flow in a deforming saturated oil reservoir. The mathematical formulation describes a fully coupled governing equation systen which consists of the equilibrium and continuity equations for three immiscible fluids flowing in a porous medium. An elastoplastic soil model, based on a Mohr–Coulomb yield surface, is used. The finite element method is applied to obtain simultaneous solutions to the governing equations where displacement and fluid pressures are the primary unknowns. The final discretized equations are solved by a direct solver using fully implicit procedures. The developed model is used to investigate the problems of three-phase fluid flow in a deforming saturated oil reservoir.     

Theoretical and numerical study of the steady‐state flow through finite fractured porous media

This paper investigates the two‐dimensional flow problem through an anisotropic porous medium containing several intersecting curved fractures. First, the governing equations of steady‐state fluid flow in a fractured porous body are summarized. The flow follows Darcy's law in matrix and Poiseuille's law in fractures. An infinite transversal permeability is considered for the fractures. A multi‐region boundary element method is used to derive a general pressure solution as a function of discharge through the fractures and the pressure and the normal flux on the domain boundary. The obtained solution fully accounts for the interaction and the intersection between fractures. A numerical procedure based on collocation method is presented to compute the unknowns on the boundaries and on the fractures. The numerical solution is validated by comparing with finite element solution or the results obtained for an infinite matrix. Pressure fields in the matrix are illustrated for domains containing several interconnected fractures, and mass balance at the intersection points is also checked. Copyright © 2013 John Wiley & Sons, Ltd.     

Multi-scale modelling of the permeability evolution of fine sands during cement suspension grouting with filtration

Fluid flow during permeation grouting of fine sands with a microcement-based grout is studied by assuming that the heterogeneous medium composed of the initial granular skeleton, filtered cement and the interstitial fluid phase can be replaced by a continuous equivalent medium at the macroscopic level. Consequently, the method of Homogenization of Periodic Structures (HPS) is used to identify the effective permeability tensor evolution under the effect of cement filtration. The expression of the macroscopic permeability tensor derived through the HPS procedure is shown to depend on the permeating fluid viscosity and the geometrical arrangement of the sand grains and cement deposit within the microstructure. Numerical computations are made using various two-dimensional and three-dimensional microstructures, and the model results are confronted with grouting experiments performed on small scale columns in the laboratory.     

Numerical solutions for flow in porous media

A numerical approach is proposed to model the flow in porous media using homogenization theory. The proposed concept involves the analyses of micro‐true flow at pore‐level and macro‐seepage flow at macro‐level. Macro‐seepage and microscopic characteristic flow equations are first derived from the Navier–Stokes equation at low Reynolds number through a two‐scale homogenization method. This homogenization method adopts an asymptotic expansion of velocity and pressure through the micro‐structures of porous media. A slightly compressible condition is introduced to express the characteristic flow through only characteristic velocity. This characteristic flow is then numerically solved using a penalty FEM scheme. Reduced integration technique is introduced for the volumetric term to avoid mesh locking. Finally, the numerical model is examined using two sets of permeability test data on clay and one set of permeability test data on sand. The numerical predictions agree well with the experimental data if constraint water film is considered for clay and two‐dimensional cross‐connection effect is included for sand. Copyright © 2003 John Wiley & Sons, Ltd.