A continuum‐discrete model using Darcy's law: formulation and verification

This paper presents a numerical scheme for fluid‐particle coupled discrete element method (DEM), which is based on poro‐elasticity. The motion of the particles is resolved by means of DEM. While within the proposition of Darcian regime, the fluid is assumed as a continuum phase on a Eulerian mesh, and the continuity equation on the fluid mesh for a compressible fluid is solved using the FEM. Analytical solutions of traditional soil mechanics examples, such as the isotropic compression and one‐dimensional upward seepage flow, were used to validate the proposed algorithm quantitatively. The numerical results showed very good agreement with the analytical solutions, which show the correctness of this algorithm. Sensitivity studies on the effect of some influential factors of the coupling scheme such as pore fluid bulk modulus, volumetric strain calculation, and fluid mesh size were performed to display the accuracy, efficiency, and robustness of the numerical algorithm. It is revealed that the pore fluid bulk modulus is a critical parameter that can affect the accuracy of the results. Because of the iterative coupling scheme of these algorithms, high value of fluid bulk modulus can result in instability and consequently reduction in the maximum possible time‐step. Furthermore, the increase of the fluid mesh size reduces the accuracy of the calculated pore pressure. This study enhances our current understanding of the capacity of fluid‐particle coupled DEM to simulate the mechanical behavior of saturated granular materials. Copyright © 2014 John Wiley & Sons, Ltd.     

Coupled bonded particle and lattice Boltzmann method for modelling fluid–solid interaction

This paper presents a two‐dimensional coupled bonded particle and lattice Boltzmann method (BPLBM) developed to simulate the fluid–solid interactions in geomechanics. In this new technique, the bonded particle model is employed to describe the inter‐particle movement and forces, and the bond between a pair of contacting particles is assumed to be broken when the tensile force or tangential force reaches a certain critical value. As a result the fracture process can be delineated based on the present model for the solid phase comprising particles, such as rocks and cohesive soils. In the meantime, the fluid phase is modelled by using the LBM, and the immersed moving boundary scheme is utilized to characterize the fluid–solid interactions. Based on the novel technique case studies have been conducted, which show that the coupled BPLBM enjoys substantially improved accuracy and enlarged range of applicability in characterizing the mechanics responses of the fluid–solid systems. Indeed such a new technique is promising for a wide range of application in soil erosion in Geotechnical Engineering, sand production phenomenon in Petroleum Engineering, fracture flow in Mining Engineering and fracture process in a variety of engineering disciplines. Copyright © 2016 John Wiley & Sons, Ltd.     

A coupled 3‐dimensional bonded discrete element and lattice Boltzmann method for fluid‐solid coupling in cohesive geomaterials

This paper presents a 3D bonded discrete element and lattice Boltzmann method for resolving the fluid‐solid interaction involving complicated fluid‐particle coupling in geomaterials. In the coupled technique, the solid material is treated as an assembly of bonded and/or granular particles. A bond model accounting for strain softening in normal contact is incorporated into the discrete element method to simulate the mechanical behaviour of geomaterials, whilst the fluid flow is solved by the lattice Boltzmann method based on kinetic theory and statistical mechanics. To provide a bridge between theory and application, a 3D algorithm of immersed moving boundary scheme was proposed for resolving fluid‐particle interaction. To demonstrate the applicability and accuracy of this coupled method, a benchmark called quicksand, in which particles become fluidised under the driving of upward fluid flow, is first carried out. The critical hydraulic gradient obtained from the numerical results matches the theoretical value. Then, numerical investigation of the performance of granular filters generated according to the well‐acknowledged design criteria is given. It is found that the proposed 3D technique is promising, and the instantaneous migration of the protected soils can be readily observed. Numerical results prove that the filters which comply with the design criteria can effectively alleviate or eliminate the appearance of particle erosion in dams.     

Pore‐scale modeling of surface erosion in a particle bed

A fully coupled transient two‐dimensional model was employed to study fundamentals of flood‐induced surface erosion in a particle bed. The interaction of the liquid and solid phases is the key mechanism related to surface erosion. The solid phase was idealized at a particle scale by using the discrete element method. The fluid phase was modeled at a mesoscale level and solved using the lattice Boltzmann method. The fluid forces applied on the particles were calculated on the basis of the momentum the fluid exchanges with the particle. The proposed approach was used to model both single particles and particle beds subjected to Couette flow conditions. The behavior of both the single particle and the particle bed depended on particle diameter and surface shear fluid velocity. The conducted simulations show that the fluid flow profile penetrates the bed for a small distance. This penetration initiates sheet‐flow and surface erosion as the fluid interacts with particles. The effect of suppressing particle rotation on the fluid‐induced forces on the particle was also examined. Suppressing particle spinning may lead to underestimated erosion rate. Results of fluid and particle velocities were compared against experimental results and appeared to agree with the observed trends.Copyright © 2013 John Wiley & Sons, Ltd.     

On the implicit immersed boundary method in coupled discrete element and lattice Boltzmann method

The coupled discrete element method and lattice Boltzmann method (DEMLBM) has increasingly drawn attention of researchers in geomechanics due to its mesoscopic nature since 2000. Immersed boundary method (IBM) and immersed moving boundary (IMB) are two popular schemes for coupling fluid particle in DEMLBM. This work aims at coupling DEM and LBM using the latest IBM algorithm and investigating its accuracy, computational efficiency, and applicability. Two benchmark tests, interstitial fluid flow in an ideal packing and single particle sedimentation in viscous fluid, are carried out to demonstrate the accuracy of IBM through semi-empirical Ergun equation, finite element method (FEM), and IMB. Then, simulations of particle migration with relatively large velocity in Poiseuille flow are utilized to address limitations of IBM in DEMLBM modeling. In addition, advantages and deficiencies of IBM are discussed and compared with IMB. It is found that the accuracy of IBM can be only guaranteed when sufficient boundary points are used and it is not suitable for geomechanical problems involving large fluid or particle velocity.     

LBM–DEM modeling of fluid–solid interaction in porous media

Three porous media flow problems, in which the fluid mechanical interactions are critical, are studied in a mesoscopic–microscopic coupling system. In this system, fluid flow in the pore space is explicitly modeled at mesoscopic level by the lattice Boltzmann method, the geometrical representation and the mechanical behavior of the solid skeleton are modeled at microscopic level by the particulate distinct element method (DEM), and the interfacial interaction between the fluid and the solids is resolved by an immersed boundary scheme. In the first benchmark problem, the well‐known and frequently utilized Ergun equation is validated in periodic particle and periodic pore models. In the second problem, the upward seepage problem is simulated over three stages: The settlement of the column of sphere under gravity loading is measured to illustrate the accuracy of the DEM scheme; the system is solved to hydrostatic state with pore space filled with fluid, showing that the buoyancy effect is captured correctly in the mesoscopic–microscopic coupling system; then, the flow with constant rate is supplied at the bottom of the column; the swelling of the ground surface and pore pressure development from the numerical simulation are compared with the predictions of the macroscopic consolidation theory. In the third problem, the fluid‐flow‐induced collapse of a sand arch inside a perforation cavity is tested to illustrate a more practical application of the developed system. Through comparing simulation results with analytical solutions, empirical law and physical laboratory observations, it is demonstrated that the developed lattice Boltzmann–distinct element coupling system is a powerful fundamental research tool for investigating hydromechanical physics in porous media flow. Copyright © 2012 John Wiley & Sons, Ltd.     

A coupled discrete element model for the simulation of soil and water flow through an orifice

Soil erosion around defective underground pipes can cause ground collapses and sinkholes in urban areas. Most of these soil erosion events are caused by fluidization of the surrounding soil with subsequent washing into defective sewer pipes. In this study, this soil erosion process is simplified as the gradual washout of sand particles mixed with water through an orifice. The discrete element method is used to simulate the large deformation behavior of the sand particles, and the Darcy fluid model is coupled with this approach to simulate fluid flow through porous sand media. A coupled 3D discrete element model is developed and implemented based on this scheme. To simulate previous experiments using this coupled model considering the current computing capacity, we incorporated a ‘supply layer’ to study the continuous erosion process. The coupled model can predict the erosion flow rates of sand and water and the shape of erosion void. Thus, the model can be used as an effective and efficient tool to investigate the soil erosion process around defective pipes. Copyright © 2017 John Wiley & Sons, Ltd.     

A unified CFD‐DEM approach for modeling of debris flow impacts on flexible barriers

This paper presents a unified modeling framework to investigate the impacts of debris flow on flexible barriers, based on coupled computational fluid dynamics and discrete element method (CFD‐DEM). We consider a debris flow as a mixture of fluid and particles where the fluid and particle phases are modeled by the CFD and the DEM, respectively. The fluid‐particle coupling is considered by the exchange of interaction forces between CFD and DEM calculations. The flexible barrier is simulated by the DEM as a network of bonded particles with remote interactions. The proposed coupled CFD‐DEM approach enables us to conveniently handle the complicated three‐way interactions among the fluid, the particles, and the flexible barrier structure for debris flow impact simulations. The proposed approach is first used to investigate the influences of channel inclination and the volumetric solid fraction in a debris mixture on the impact force, the resultant deformation, and the retained mass in a flexible barrier. The predictions agree well with existing experimental and numerical studies. We further examine the possible failure modes of a flexible barrier under debris flow impact and their underlying mechanisms. The performance of different components in a flexible barrier system, including single wires, double twists and cables, and their load sharing mechanisms, are carefully evaluated. The proposed unified framework offers a novel, promising pathway towards physically based, quantitative analysis and design of flexible barriers for debris flow mitigation.     

Meshless numerical modeling of brittle–viscous deformation: first results on boudinage and hydrofracturing using a coupling of discrete element method (DEM) and smoothed particle hydrodynamics (SPH)

We have developed a new approach for the numerical modeling of deformation processes combining brittle fracture and viscous flow. The new approach is based on the combination of two meshless particle-based methods: the discrete element method (DEM) for the brittle part of the model and smooth particle hydrodynamics (SPH) for the viscous part. Both methods are well established in their respective application domains. The two methods are coupled at the particle scale, with two different coupling mechanisms explored: one is where DEM particles act as virtual SPH particles and one where SPH particles are treated like DEM particles when interacting with other DEM particles. The suitability of the combined approach is demonstrated by applying it to two geological processes, boudinage, and hydrofracturing, which involve the coupled deformation of a brittle solid and a viscous fluid. Initial results for those applications show that the new approach has strong potential for the numerical modeling of coupled brittle–viscous deformation processes.     

The influence of the void fraction on the particle migration: A coupled computational fluid dynamics–discrete element method study about drag force correlations

Granular soils subjected to flow through their soil skeleton can show a behaviour in which fine particles migrate through the pore space between coarser particles. This process is called internal instability or suffusion. This contribution deals with the numerical analysis of the migration of fine particles in a soil column subjected to fluid flow with unresolved coupled computational fluid dynamics–discrete element method (CFD–DEM) with special regards to the used drag force correlation. The contribution investigates the influence of the Schiller–Naumann model and its extension with a voidage term on the migration behaviour of fine particles. The voidage term is further varied with a parameter, which controls the impact of the change of the void fraction on the drag force. It could be observed that the Schiller–Naumann model does not yield in a suffusive behaviour while the extended models show significant particle migration. Thereby, increasing the impact of the void fraction on the drag force results in stronger particle migration. These results reveal the need for good validation techniques. They indicate how the drag force correlation can be adapted to depict the correct particle migration behaviour.     

Numerical convergence study of iterative coupling for coupled flow and geomechanics

In this paper, we consider algorithms for modeling complex processes in porous media that include fluid and structure interactions. Numerous field applications would benefit from a better understanding and integration of porous flow and solid deformation. Important applications in environmental and petroleum engineering include carbon sequestration, surface subsidence, pore collapse, cavity generation, hydraulic fracturing, thermal fracturing, wellbore collapse, sand production, fault activation, and waste disposal, while similar issues arise in biosciences and chemical sciences as well. Here, we consider solving iteratively the coupling of flow and mechanics. We employ mixed finite element method for flow and a continuous Galerkin method for elasticity. For single-phase flow, we demonstrate the convergence and convergence rates for two widely used schemes, the undrained split and the fixed stress split. We discuss the extension of the fixed stress iterative coupling scheme to an equation of state compositional flow model coupled with elasticity and a single-phase poroelasticity model on general hexahedral grids. Computational results are presented.     

Extended CFD–DEM for free‐surface flow with multi‐size granules

Computational fluid dynamics and discrete element method (CFD–DEM) is extended with the volume of fluid (VOF) method to model free‐surface flows. The fluid is described on coarse CFD grids by solving locally averaged Navier–Stokes equations, and particles are modelled individually in DEM. Fluid–particle interactions are achieved by exchanging information between DEM and CFD. An advection equation is applied to solve the phase fraction of liquid, in the spirit of VOF, to capture the dynamics of free fluid surface. It also allows inter‐phase volume replacements between the fluid and solid particles. Further, as the size ratio (SR) of fluid cell to particle diameter is limited (i.e. no less than 4) in coarse‐grid CFD–DEM, a porous sphere method is adopted to permit a wider range of particle size without sacrificing the resolution of fluid grids. It makes use of more fluid cells to calculate local porosities. The developed solver (cfdemSolverVOF) is validated in different cases. A dam break case validates the CFD‐component and VOF‐component. Particle sedimentation tests validate the CFD–DEM interaction at various Reynolds numbers. Water‐level rising tests validate the volume exchange among phases. The porous sphere model is validated in both static and dynamic situations. Sensitivity analyses show that the SR can be reduced to 1 using the porous sphere approach, with the accuracy of analyses maintained. This allows more details of the fluid phase to be revealed in the analyses and enhances the applicability of the proposed model to geotechnical problems, where a highly dynamic fluid velocity and a wide range of particle sizes are encountered. Copyright © 2015 John Wiley & Sons, Ltd.     

A meshless method for axisymmetric problems in continuously nonhomogeneous saturated porous media

A meshless method based on the local Petrov–Galerkin approach is proposed to analyze 3-d axisymmetric problems in porous functionally graded materials. Constitutive equations for porous materials possess a coupling between mechanical displacements for solid and fluid phases. The work is based on the u–u formulation and the incognita fields of the coupled analysis in focus are the solid skeleton displacements and the fluid displacements. Independent spatial discretization is considered for each phase of the model, rendering a more flexible and efficient methodology. Both displacements are approximated by the moving least-squares (MLS) scheme. The paper presents in the first time a general meshless method for the numerical analysis of axisymmetric problems in continuously nonhomogeneous saturated porous media. Numerical results are given for boreholes in continuously nonhomogeneous porous medium with prescribed misfit and exponential variation of material parameters in the excavation zone.     

A coupled two‐phase fluid flow and elastoplastic deformation model for unsaturated soils: theory,implementation, and application

Although numerous numerical models have been proposed for simulating the coupled hydromechanical behaviors in unsaturated soils, few studies satisfactorily reproduced the soil–water–air three‐phase coupling processes. Particularly, the impacts of deformation dependence of water retention curve, bonding stress, and gas flow on the coupled processes were less examined within a coupled soil–water–air model. Based on our newly developed constitutive models (Hu et al., 2013, 2014, 2015) in which the soil–water–air couplings have been appropriately captured, this study develops a computer code named F²Mus3D to investigate the coupled processes with a focus on the above impacts. In the numerical implementation, the generalized‐α time integration scheme was adopted to solve the equations, and a return‐mapping implicit stress integration scheme was used to update the state variables. The numerical model was verified by two well‐designed laboratory tests and was applied for modeling the coupled elastoplastic deformation and two‐phase fluid flow processes in a homogenous soil slope induced by rainfall infiltration. The simulation results demonstrated that the numerical model well reproduces the initiation of a sheared zone at the toe of the slope and its propagation toward the crest as the rain infiltration proceeds, which manifests a typical mechanism for rainfall‐induced shallow landslides. The simulated plastic strain and deformation would be remarkably underestimated when the bonding stress and/or the deformation‐dependent nature of hydraulic properties are ignored in the coupled model. But on the contrary, the negligence of gas flow in the slope soil results in an overestimation of the rainfall‐induced deformation. Copyright © 2015 John Wiley & Sons, Ltd.     

Coupled discrete element and finite volume solution of two classical soil mechanics problems

One dimensional solutions for the classic critical upward seepage gradient/quick condition and the time rate of consolidation problems are obtained using coupled routines for the finite volume method (FVM) and discrete element method (DEM), and the results compared with the analytical solutions. The two phase flow in a system composed of fluid and solid is simulated with the fluid phase modeled by solving the averaged Navier–Stokes equation using the FVM and the solid phase is modeled using the DEM. A framework is described for the coupling of two open source computer codes: YADE-OpenDEM for the discrete element method and OpenFOAM for the computational fluid dynamics. The particle–fluid interaction is quantified using a semi-empirical relationship proposed by Ergun [12]. The two classical verification problems are used to explore issues encountered when using coupled flow DEM codes, namely, the appropriate time step size for both the fluid and mechanical solution processes, the choice of the viscous damping coefficient, and the number of solid particles per finite fluid volume.     

High‐order ADE scheme for solving the fluid diffusion equation in non‐uniform grids and its application in coupled hydro‐mechanical simulation

To improve the stability and efficiency of explicit technique, one proposed method is to use an unconditionally stable alternating direction explicit (ADE) scheme. However, the standard ADE scheme is only moderately accurate and restricted to uniform grids. This paper derives a novel high‐order ADE scheme capable of solving the fluid diffusion equation in non‐uniform grids. The new scheme is derived by performing a fourth‐order finite difference approximation to the spatial derivatives of the diffusion equation in non‐uniform grid. The implicit Crank‐Nicolson technique is then applied to the resulting approximation, and the subsequent equation is split into two alternating direction sweeps, giving rise to a new high‐order ADE scheme. Because the new scheme can be potentially applied in coupled hydro‐mechanical (H‐M) simulation, the pore pressure solutions from the new scheme are then sequentially coupled with an existing geomechanical simulator in the computer program Fast Lagrangian Analysis of Continua. This coupling procedure is called the sequentially explicit coupling technique based on the fourth‐order ADE scheme (SEA‐4). Verifications of well‐known consolidation problems showed that the new ADE scheme and SEA‐4 can reduce computer runtime by 46% to 75% to that of Fast Lagrangian Analysis of Continua's basic scheme. At the same time, the techniques still maintained average percentage error of 1.6% to 3.5% for pore pressure and 0.2% to 1.5% for displacement solutions and were still accurate under typical grid non‐uniformities. This result suggests that the new high‐order ADE scheme can provide an efficient explicit technique for solving the flow equation of a coupled H‐M problem, which will be beneficial for large‐scale and long‐term H‐M problems in geoengineering.     

Modelling hydraulic and capillary-driven two-phase fluid flow in unsaturated concretes at the meso-scale with a unique coupled DEM-CFD technique

The goal of the research was to demonstrate the impact of thin porous interfacial transition zones (ITZs) between aggregates and cement matrix on fluid flow in unsaturated concrete caused by hydraulic/capillary pressure. To demonstrate this impact, a novel coupled approach to simulate the two-phase (water and moist air) flow of hydraulically and capillary-driven fluid in unsaturated concrete was developed. By merging the discrete element method (DEM) with computational fluid dynamics (CFD) under isothermal settings, the process was numerically studied at the meso-scale in two-dimensional conditions. A flow network was used to describe fluid behaviour in a continuous domain between particles. Small concrete specimens of a simplified particle mesostructure were subjected to fully coupled hydro-mechanical simulation tests. A simple uniaxial compression test was used to calibrate the pure DEM represented by bonded spheres, while a permeability and sorptivity test for an assembly of spheres was used to calibrate the pure CFD. For simplified specimens of the pure cement matrix, cement matrix with aggregate, and cement matrix with aggregate and ITZ of a given thickness, DEM/CFD simulations were performed sequentially. The numerical results of permeability and sorptivity were directly compared to the data found in the literature. A satisfactory agreement was achieved. Porous ITZs in concrete were found to reduce sorption by slowing the capillary-driven fluid flow, and to speed the full saturation of pores when sufficiently high hydraulic water pressures were dominant.     

A continuum‐discrete hydromechanical analysis of granular deposit liquefaction

A coupled continuum‐discrete hydromechanical model was employed to analyse the liquefaction of a saturated loose deposit of cohesionless particles when subjected to a dynamic base excitation. The pore fluid flow was idealized using averaged Navier–Stokes equations and the discrete element method was employed to model the solid phase particles. A well established semi‐empirical relationship was utilized to quantify the fluid–particle interactions. The conducted simulations revealed a number of salient micro‐mechanical mechanisms and response patterns associated with the deposit liquefaction. Space and time variation of porosity was a major factor which affected the coupled response of the solid and fluid phases. Pore fluid flow was within Darcy's regime. The predicted response exhibited macroscopic patterns consistent with experimental results and case histories of the liquefaction of granular soil deposits. Copyright © 2004 John Wiley & Sons, Ltd.