A probabilistic approach for computing water retention of particulate systems from statistics of grain size and tessellated pore network

The paper offers an analytical determination of the hydraulic properties of an unsaturated soil with reference to its retention curve, which describes the relationship between the volumetric water content and capillarity through matric suction. The analysis combines a particulate approach focused on the physics at the pore scale, including microstructural aspects, with a probabilistic approach where the void space and grain size are considered as random variables. In the end, the soil water characteristic curve of an unsaturated granular medium along a drying path can be derived analytically based on the sole information of particle size distribution. The analysis hinges on the tessellation of a wet granular system into an assemblage of tetrahedral unit cells revealing a pore network upon which capillary physics are computed with respect to pore throat invasion by a non-wetting fluid with evolving pendular capillary bridges. The crux of the paper is to pass from particle size probability distribution to a matching void space distribution to eventually reveal key information such as void cell and solid volume statistics. Making reasonable statistically based assumptions to render calculations tractable, the water retention curve can be readily constructed. Model predictions compare quite favourably with experimental data available for actual soils, especially in the high saturation range. Having a sound scientific basis, the model can be made amenable to address a variety of soils with a wider range of particle sizes.     

Improved algorithms for generalized finite element simulations of three-dimensional hydraulic fracture propagation

This paper reports improvements to algorithms for the simulation of 3-D hydraulic fracturing with the Generalized Finite Element Method (GFEM). Three optimizations are presented and analyzed. First, an improved initial guess based on solving a 3-D elastic problem with the pressure from the previous step is shown to decrease the number of Newton iterations and increase robustness. Second, an improved methodology to find the time step that leads to fracture propagation is proposed and shown to decrease significantly the number of iterations. Third, reduced computational cost is observed by properly recycling the linear part of the coupled stiffness matrix. Two representative examples are used to analyze these improvements. Additionally, a methodology to include the leak-off term is presented and verified against asymptotic analytical solutions. Conservation of mass is shown to be well satisfied in all examples.     

A new numerical approach to the hyperstatic reaction method for segmental tunnel linings

This paper proposes a numerical approach to the hyperstatic reaction method (HRM) for the analysis of segmental tunnel linings. The influence of segmental joints has been considered directly using a fixity ratio that is determined on the basis of the rotational stiffness. The parameters necessary for the calculation are presented. A specific implementation has been developed using a FEM framework. This code is able to consider the three‐dimensional (3D) effect of segment joints in successive rings on the tunnel lining behaviour. The present HRM allows one to take an arbitrary distribution of segment joints along the tunnel boundary into consideration. In addition, the rotational stiffness of segment joints has been simulated using nonlinear behaviour, as it is closer to the true behaviour of a joint than linear or bilinear behaviour. The numerical results of three hypotheses on ring interaction, which allow the 3D effect of a segmental tunnel lining to be taken into account, have been compared with data obtained from the shield‐driven tunnel of the Bologna–Florence high‐speed railway line project. The numerical results presented in the paper show that the proposed HRM can be used to effectively estimate the behaviour of a segmental tunnel lining. Copyright © 2014 John Wiley & Sons, Ltd.     

Granular element method for three‐dimensional discrete element calculations

This paper endows the recently‐proposed granular element method (GEM) with the ability to perform 3D discrete element calculations. By using non‐uniform rational B‐Splines to accurately represent complex grain geometries, we proposed an alternative approach to clustering‐based and polyhedra‐based discrete element methods whereby the need for complicated and ad hoc approaches to construct 3D grain geometries is entirely bypassed. We demonstrate the ability of GEM in capturing arbitrary‐shaped 3D grains with great ease, flexibility, and without excessive geometric information. Furthermore, the applicability of GEM is enhanced by its tight integration with existing non‐uniform rational B‐Splines modeling tools and ability to provide a seamless transition from binary images of real grain shapes (e.g., from 3D X‐ray CT) to modeling and discrete mechanics computations.Copyright © 2013 John Wiley & Sons, Ltd.     

New graphical methods for estimating aquifer hydraulic parameters using pumping tests with exponentially decreasing rates

Actual pumping tests may involve continuously decreasing rates over a certain period of time, and the hydraulic conductivity (K) and specific storage (S_s) of the tested confined aquifer cannot be interpreted from the classical constant‐rate test model. In this study, we revisit the aquifer drawdown characteristics of a pumping test with an exponentially decreasing rate using the dimensionless analytical solution for such a variable‐rate model. The drawdown may decrease with time for a short period of time at intermediate pumping times for such pumping tests. A larger ratio of initial to final pumping rate and a smaller radial distance of the observation well will enhance the decreasing feature. A larger decay constant results in an earlier decrease, but it weakens the extent of such a decrease. Based on the proposed dimensionless transformation, we have proposed two graphical methods for estimating K and S_s of the tested aquifer. The first is a new type curve method that does not employ the well function as commonly done in standard type curve analysis. Another is a new analytic method that takes advantage of the decreasing features of aquifer drawdown during the intermediate pumping stage. We have demonstrated the applicability and robustness of the two new graphical methods for aquifer characterization through a synthetic pumping test.     

Modelling landslides in unsaturated slopes subjected to rainfall infiltration using material point method

This paper presents a dynamic fully coupled formulation for saturated and unsaturated soils that undergo large deformations based on material point method. Governing equations are applied to porous material while considering it as a continuum in which the pores of the solid skeleton are filled with water and air. The accuracy of the developed method is tested with available experimental and numerical results. The developed method has been applied to investigate the failure and post‐failure behaviour of rapid landslides in unsaturated slopes subjected to rainfall infiltration using two different bedrock geometries that lie below the top soil. The models show different failure and post‐failure mechanisms depending on the bedrock geometry and highlight the negative effects of continuous rain infiltrations. Copyright © 2016 John Wiley & Sons, Ltd.     

Large‐scale poroelastic fractured reservoirs modeling using the fast multipole displacement discontinuity method

An effective approach to modeling the geomechanical behavior of the network and its permeability variation is to use a poroelastic displacement discontinuity method (DDM). However, the approach becomes rather computationally intensive for an extensive system of cracks, particularly when considering coupled diffusion/deformation processes. This is because of additional unknowns and the need for time‐marching schemes for the numerical integration. The Fast Multipole Method (FMM) is a technique that can accelerate the solution of large fracture problems with linear complexity with the number of unknowns both in memory and CPU time. Previous works combining DDM and FMM for large‐scale problems have accounted only for elastic rocks, neglecting the fluid leak‐off from the fractures into the matrix and its influence on pore pressure and stress field. In this work we develop an efficient geomechanical model for large‐scale natural fracture networks in poroelastic reservoirs with fracture flow in response to injection and production operations. Accuracy and computational performance of the proposed method with those of conventional poroelastic DDM are compared through several case studies involving up to several tens of thousands of boundary elements. The results show the effectiveness of the FMM approach to successfully evaluate field‐scale problems for the design of exploitation strategies in unconventional geothermal and petroleum reservoirs. An example considering faults reveals the impact of reservoir compartmentalization because of sealing faults for both geomechanical and flow variables under elastic and poroelastic rocks. Copyright © 2015 John Wiley & Sons, Ltd.     

Parallelized combined finite-discrete element (FDEM) procedure using multi-GPU with CUDA

This paper focuses on the efficiency of finite discrete element method (FDEM) algorithmic procedures in massive computers and analyzes the time-consuming part of contact detection and interaction computations in the numerical solution. A detailed operable GPU parallel procedure was designed for the element node force calculation, contact detection, and contact interaction with thread allocation and data access based on the CUDA computing. The emphasis is on the parallel optimization of time-consuming contact detection based on load balance and GPU architecture. A CUDA FDEM parallel program was developed with the overall speedup ratio over 53 times after the fracture from the efficiency and fidelity performance test of models of in situ stress, UCS, and BD simulations in Intel i7-7700K CPU and the NVIDIA TITAN Z GPU. The CUDA FDEM parallel computing improves the computational efficiency significantly compared with the CPU-based ones with the same reliability, providing conditions for achieving larger-scale simulations of fracture.     

Numerical and analytical observations on long and infinite slopes

Interest in the mechanics of landslides has led to renewed evaluation of the infinite slope equations, and the need for a more general framework for estimating the factor of safety of long and infinite slopes involving non‐homogeneous soil profiles. The paper describes finite element methods that demonstrate the potential for predicting failure in long slope profiles where the critical mechanism is not necessarily at the base of the soil layer. The influence of slope angle is also examined in long slopes, leading to some counter‐intuitive conclusions about the impact of slope steepness on the factor of safety. Copyright © 2010 John Wiley & Sons, Ltd.