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.     

On the pore‐scale mechanisms leading to brittle and ductile deformation behavior of crystalline rocks

Deformation mechanisms at the pore scale are responsible for producing large strains in porous rocks. They include cataclastic flow, dislocation creep, dynamic recrystallization, diffusive mass transfer, and grain boundary sliding, among others. In this paper, we focus on two dominant pore‐scale mechanisms resulting from purely mechanical, isothermal loading: crystal plasticity and crofracturing. We examine the contributions of each mechanism to the overall behavior at a scale larger than the grains but smaller than the specimen, which is commonly referred to as the mesoscale. Crystal plasticity is assumed to occur as dislocations along the many crystallographic slip planes, whereas microfracturing entails slip and frictional sliding on microcracks. It is observed that under combined shear and tensile loading, microfracturing generates a softer response compared with crystal plasticity alone, which is attributed to slip weakening where the shear stress drops to a residual level determined by the frictional strength. For compressive loading, however, microfracturing produces a stiffer response than crystal plasticity because of the presence of frictional resistance on the slip surface. Behaviors under tensile, compressive, and shear loading invariably show that porosity plays a critical role in the initiation of the deformation mechanisms. Both crystal plasticity and microfracturing are observed to initiate at the peripheries of the pores, consistent with results of experimental studies. Copyright © 2015 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.     

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.     

A discontinuous finite element suspended sediment transport model for water quality assessments in river networks

Suspended sediment plays an important role in the distribution and transport of many pollutants (such as radionuclides) in rivers. Pollutants may adsorb on fine suspended particles (e.g. clay) and spread according to the suspended sediment movement. Hence, the simulation of the suspended sediment mechanism is indispensable for realistic transport modelling. This paper presents and tests a simple mathematical model for predicting the suspended sediment transport in river networks. The model is based on the van Rijn suspended load formula and the advection–diffusion equation with a source or sink term that represents the erosion or deposition fluxes. The transport equation is solved numerically with the discontinuous finite element method. The model evaluation was performed in two steps, first by comparing model simulations with the measured suspended sediment concentrations in the Grote Nete–Molse Nete River in Belgium, and second by a model intercomparison with the sediment transport model NST MIKE 11. The simulations reflect the measurements with a Nash‐Sutcliffe model efficiency of 0.6, while the efficiency between the proposed model and the NST MIKE 11 simulations is 0.96. Both evaluations indicate that the proposed sediment transport model, that is sufficiently simple to be practical, is providing realistic results.     

Analysis of shaft resistance of jacked piles in sands

In recent years, pile jacking has become a viable alternative installation method for displacement piles. Pile jacking produces minimal noise, vibration and air pollution during installation. In addition, it is possible, at the end of jacking, to have a good estimate of the ultimate static capacity of the pile. In this paper, the shaft resistance of piles jacked into sand is studied using one‐dimensional finite element analysis. The finite element simulations, using a two‐surface plasticity model, demonstrate the effects of relative density and confinement on the unit shaft resistance of piles jacked in sand. The impact of the number of jacking strokes on the unit shaft capacity is also assessed. Based on the numerical results, we developed equations for shaft resistance quantifying the effects of relative density, initial confinement and number of jacking strokes. Predictions using these equations are compared with data obtained from centrifuge tests and field tests. Copyright © 2010 John Wiley & Sons, Ltd.     

Review and enhancement of 3D concrete models for large‐scale numerical simulations of concrete structures

The present paper focuses on selected plasticity and damage‐plasticity models for describing the 3D material behavior of concrete. In particular, a plasticity model and a damage‐plasticity model are reviewed and evaluated. Based on the results of the evaluation, enhancements are proposed, aiming at improving the correspondence between predicted and observed material behavior and aiming at implementing a robust and efficient stress update algorithm in a finite element program for performing large‐scale 3D numerical simulations of concrete structures. The capabilities of the concrete models are demonstrated by 3D numerical simulations of benchmark tests with combined bending and torsional loading and combined compression and shear loading and by a large‐scale 3D finite element analysis of a model test of a concrete arch dam. Copyright © 2011 John Wiley & Sons, Ltd.