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.     

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 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.     

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.     

Pore‐scale modeling of fluid‐particles interaction and emerging poromechanical effects

A micro‐hydromechanical model for granular materials is presented. It combines the discrete element method for the modeling of the solid phase and a pore‐scale finite volume formulation for the flow of an incompressible pore fluid. The coupling equations are derived and contrasted against the equations of conventional poroelasticity. An analogy is found between the discrete element method pore‐scale finite volume coupling and Biot's theory in the limit case of incompressible phases. The simulation of an oedometer test validates the coupling scheme and demonstrates the ability of the model to capture strong poromechanical effects. A detailed analysis of microscale strain and stress confirms the analogy with poroelasticity. An immersed deposition problem is finally simulated and shows the potential of the method to handle phase transitions. Copyright © 2013 John Wiley & Sons, Ltd.     

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.     

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.     

Two-dimensional Dynamics Simulation of Two-phase Debris Flow

To investigate the movement mechanism of debris flow, a two‐dimensional, two‐phase, depth‐integrated model is introduced. The model uses Mohr‐Coulomb plasticity for the solid rheology, and the fluid stress is modeled as a Newtonian fluid. The interaction between solid and liquid phases, which plays a major role in debris flow movement, is assumed to consist of drag and buoyancy forces. The applicability of drag force formulas is discussed. Considering the complex interaction between debris flow and the bed surface, a combined friction boundary condition is imposed on the bottom, and this is also discussed. To solve the complex model equations, a numerical method with second‐order accuracy based on the finite volume method is proposed. Several numerical experiments are performed to verify the feasibilities of model and numerical schemes. Numerical results demonstrate that different solid volume fractions substantially affect debris flow movement.     

Modeling of fluid–solid interaction in granular media with coupled lattice Boltzmann/discrete element methods: application to piping erosion

In this article, we present a numerical method to deal with fluid–solid interactions and simulate particle–fluid systems as encountered in soils. This method is based on a coupling between two methods, now widely used in mechanics of granular media and fluid dynamics respectively: the discrete element (DE) method and the lattice Boltzmann (LB) method. The DE method is employed to model interactions between particles, whereas the LB method is used to describe an interstitial Newtonian fluid flow. The coupling presented here is a full one in the sense that particle motions act on fluid flow and reciprocally. This article presents in details each of the two methods and the principle of the coupling scheme. Determination of hydrodynamic forces and torques is also detailed, and the treatment of boundaries is explained. The coupled method is finally illustrated on a simple example of piping erosion, which puts in evidence that the combined LB–DE scheme constitutes a promising tool to study coupled problems in geomechanics. Copyright © 2011 John Wiley & Sons, Ltd.     

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.     

Frictionless contact formulation for dynamic analysis of nonlinear saturated porous media based on the mortar method

A finite element algorithm for frictionless contact problems in a two‐phase saturated porous medium, considering finite deformation and inertia effects, has been formulated and implemented in a finite element programme. The mechanical behaviour of the saturated porous medium is predicted using mixture theory, which models the dynamic advection of fluids through a fully saturated porous solid matrix. The resulting mixed formulation predicts all field variables including the solid displacement, pore fluid pressure and Darcy velocity of the pore fluid. The contact constraints arising from the requirement for continuity of the contact traction, as well as the fluid flow across the contact interface, are enforced using a penalty approach that is regularised with an augmented Lagrangian method. The contact formulation is based on a mortar segment‐to‐segment scheme that allows the interpolation functions of the contact elements to be of order N. The main thrust of this paper is therefore how to deal with contact interfaces in problems that involve both dynamics and consolidation and possibly large deformations of porous media. The numerical algorithm is first verified using several illustrative examples. This algorithm is then employed to solve a pipe‐seabed interaction problem, involving large deformations and dynamic effects, and the results of the analysis are also compared with those obtained using a node‐to‐segment contact algorithm. The results of this study indicate that the proposed method is able to solve the highly nonlinear problem of dynamic soil–structure interaction when coupled with pore water pressures and Darcy velocity. Copyright © 2015 John Wiley & Sons, Ltd.     

Coupled Model of Two-phase Debris Flow, Sediment Transport and Morphological Evolution

The volume fraction of the solid and liquid phase of debris flows,which evolves simultaneously across terrains,largely determines the dynamic property of debris flows. The entrainment process significantly influences the amplitude of the volume fraction. In this paper,we present a depth-averaged two-phase debris-flow model describing the simultaneous evolution of the phase velocity and depth,the solid and fluid volume fractions and the bed morphological evolution. The model employs the Mohr–Coulomb plasticity for the solid stress,and the fluid stress is modeled as a Newtonian viscous stress. The interfacial momentum transfer includes viscous drag and buoyancy. A new extended entrainment rate formula that satisfies the boundary momentum jump condition(Iverson and Ouyang,2015) is presented. In this formula,the basal traction stress is a function of the solid volume fraction and can take advantage of both the Coulomb and velocity-dependent friction models. A finite volume method using Roe's Riemann approximation is suggested to solve the equations. Three computational cases are conducted and compared with experiments or previous results. The results show that the current computational model and framework are robust and suitable for capturing the characteristics of debris flows.     

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.     

颗粒沉降问题的高解析度CFD-DEM-IBM流固耦合数值模拟方法

为实现颗粒沉降问题的精确数值模拟,探索稠密颗粒两相流中流体与颗粒间的影响机制,基于浸入边界法（IBM）,建立了高解析度CFD-DEM-IBM流固耦合数值模拟方法。分别通过流体动力学方法（CFD）及离散单元法（DEM）描述连续流体及非连续颗粒,引入浸入边界法处理颗粒移动边界,并在Navier-Stokes动量方程中附加体力项以体现流固相互作用,采用交错迭代算法在同一时间步内多次迭代求解直至收敛以实现流固间的强耦合。通过单、双、群颗粒的沉降行为模拟,结果表明:与传统CFD-DEM方法相比,该方法能够准确考虑颗粒间及流固两相间的相互作用,获得高解析度的流场信息。模拟结果与前人数值计算结果吻合,验证了该方法的准确性与有效性及其对稠密颗粒两相流问题的适用性与优越性。     

Coupled three-dimensional discrete element–finite difference simulation of dynamic compaction

Discrete element method has been widely adopted to simulate processes that are challenging to continuum-based approaches. However, its computational efficiency can be greatly compromised when large number of particles are required to model regions of less interest to researchers. Due to this, the application of DEM to boundary value problems has been limited. This paper introduces a three-dimensional discrete element–finite difference coupling method, in which the discrete–continuum interactions are modeled in local coordinate systems where the force and displacement compatibilities between the coupled subdomains are considered. The method is validated using a model dynamic compaction test on sand. The comparison between the numerical and physical test results shows that the coupling method can effectively simulate the dynamic compaction process. The responses of the DEM model show that dynamic stress propagation (compaction mechanism) and tamper penetration (bearing capacity mechanism) play very different roles in soil deformations. Under impact loading, the soil undergoes a transient weakening process induced by dynamic stress propagation, which makes the soil easier to densify under bearing capacity mechanism. The distribution of tamping energy between the two mechanisms can influence the compaction efficiency, and allocating higher compaction energy to bearing capacity mechanism could improve the efficiency of dynamic compaction.

     

A practical and efficient numerical scheme for the analysis of steady state unconfined seepage flows

The scaled boundary finite‐element method (SBFEM), a novel semi‐analytical technique, is applied to the analysis of the confined and unconfined seepage flow. This method combines the advantages of the finite‐element method and the boundary element method. In this method, only the boundary of the domain is discretized; no fundamental solution is required, and singularity problems can be modeled rigorously. Anisotropic and nonhomogeneous materials satisfying similarity are modeled without additional efforts. In this paper, SBFE equations and solution procedures for the analysis of seepage flow are outlined. The accuracy of the proposed method in modeling singularity problems is demonstrated by analyzing seepage flow under a concrete dam with a cutoff at heel. As only the boundary is discretized, the variable mesh technique is advisable for modeling unconfined seepage analyses. The accuracy, effectiveness, and efficiency of the method are demonstrated by modeling several unconfined seepage flow problems. Copyright © 2011 John Wiley & Sons, Ltd.     

水气二相流方程的一种离散数值解法

提出了水气二相流方程的一种数值解法.在利用有限元方法求解水气二相流方程时,引入了离散Newton迭代方法,用于非线性有限元方程组的线性化处理,将这一步计算的收敛阶由原有研究的线性收敛提高到平方收敛,并避免了直接应用Newton迭代方法给编程带来的不便.同时在求解两相的有限元方程组时,采用两相方程组并行迭代的方法,与联立计算相比节省了大量的内存空间.     

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.