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.     

A comprehensive study on the parameters setting in smoothed particle hydrodynamics (SPH) method applied to hydrodynamics problems

Smoothed particle hydrodynamics (SPH) is a meshfree, Lagrangian particle method which has advantages in handling solids with extremely large deformation. Like any other numerical methods, cares must be taken to ensure its desirable accuracy and stability through considering several correction techniques in calculation. The selection of values for parameters in those correction approaches is a key step in SPH simulation, which is always difficult for new beginners to deal well with effectively. This paper examines the common inconsistency and instability problems in SPH method and studies its computational efficiency when applied to hydrodynamics problems with material strength like soil column collapse. We analyzed in detail how the correction techniques mitigate these inconsistency and instability problems. Also, the numerical testing results associate with different values for the parameters used in the correction techniques are provided for better understanding the influence of these parameters and for finding out the desirable values. It is found that (1) the SPH method is easily subjected to an inconsistency problem in the boundary area due to the boundary deficiency, and it can be treated well by adopting “virtual particles” contributing to the particle summations. (2) The numerical oscillation in SPH simulation can be mitigated effectively by artificial viscosity with the suggested parameter values. (3) The tension cracking treatment, artificial viscosity and artificial stress work well in removing the tensile instability problem in SPH method. In addition, the nearest neighboring particle searching (NNPS) algorithm, spacing ratio, smoothing length and time step influence the efficiency and accuracy of SPH method significantly. It is shown that SPH method with suggested parameters values can produce a very good result compared with the experimental result.     

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.     

Hierarchical multiscale numerical modelling of internal erosion with discrete and finite elements

This paper presents a coupled finite and discrete-element model (FEM and DEM) to simulate internal erosion. The model is based on ICY, an interface between COMSOL, an FEM engine, and YADE, a DEM code. With this model, smaller DEM subdomains are generated to simulate particle displacements at the grain scale. Particles in these small subdomains are subjected to buoyancy, gravity, drag and contact forces for short time steps (0.1 s). The DEM subdomains provide the macroscale (continuum) model with a particle flux distribution. Through a mass conservation equation, the flux distribution allows changes in porosity, hydraulic conductivity and hydraulic gradient to be evaluated for the same time steps at a larger, continuum scale. The updated hydraulic gradients from the continuum model provide the DEM subdomains with updated hydrodynamic forces based on a coarse-grid method. The number of particles in the DEM subdomains is also updated based on the new porosity distribution. The hierarchical multiscale model (HMM) was validated with the simulation of suffusion. Results for the proposed HMM algorithm are consistent with results based on a DEM model incorporating the full sample and simulation duration. The proposed HMM algorithm could enable the modelling of internal erosion for soil volumes that are too large to be modelled with a single DEM subdomain.

     

SPH–DEM coupling method based on GPU and its application to the landslide tsunami. Part II: reproduction of the Vajont landslide tsunami

Acta Geotechnica - Landslide tsunamis are complex fluid–solid coupling processes that often cause enormous catastrophes. In this study, the smooth particle hydrodynamics (SPH) and discrete...     

Simulation of debris flow on an instrumented test slope using an updated Lagrangian continuum particle method

We present an updated Lagrangian continuum particle method based on smoothed particle hydrodynamics (SPH) for simulating debris flow on an instrumented test slope. The site is a deforested area near the village of Ruedlingen, a community in the canton of Schaffhausen in Switzerland. Artificial rainfall experiments were conducted on the slope that led to failure of the sediment in the form of a debris flow. We develop a 3D mechanistic model for this test slope and conduct numerical simulations of the flow kinematics using an SPH formulation that captures large deformation, material nonlinearity, and the complex post-failure movement of the sediment. Two main simulations explore the impact of changes in the mechanical properties of the sediment on the ensuing kinematics of the flow. The first simulation models the sediment as a granular homogeneous material, while the second simulation models the sediment as a heterogeneous material with spatially varying cohesion. The variable cohesion is meant to represent the effects of root reinforcement from vegetation. By comparing the numerical solutions with the observed failure surfaces and final free-surface geometries of the debris deposit, as well as with the observed flow velocity, flow duration, and hot spots of strain concentration, we provide insights into the accuracy and robustness of the SPH framework for modeling debris flows.

     

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.     

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.     

颗粒接触无网格法及其在边坡成灾范围模拟中的应用

提出了一种基于颗粒接触的无网格方法（PCMM）,并编制了相应的C++程序,解决了有限元等网格类方法在模拟边坡失稳滑动过程中的网格畸变问题。该方法利用颗粒离散元中的接触拓扑创建连续介质单元,通过颗粒的运动演化实现连续介质单元的自动删除及重建,通过在连续介质单元中引入考虑应变软化效应的Mohr-Coulomb模型及最大拉应力模型,实现边坡的弹塑性分析及失稳滑移过程的模拟。利用PCMM分析了均质边坡的弹塑性场、沙堆的形成过程及土质边坡的失稳过程。计算结果表明,PCMM在小变形下具有足够的精度,且在模拟材料大变形方面具有明显优势,是一种模拟边坡成灾范围的有效方法。     

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 three-dimensional discrete element-lattice Boltzmann methods for fluid-solid interaction with polyhedral particles

Interaction between solid particles and fluid is of fundamental interest to scientists and engineers in many different applications—cardiopulmonary flows, aircraft and automobile aerodynamics, and wind loading on buildings to name a few. In geomechanics, particle shape significantly affects both particle-particle and particle-fluid interaction. Herein, we present a generalized method for modeling the interaction of arbitrarily shaped polyhedral particles and particle assemblages with fluid using a coupled discrete element method (DEM) and lattice Boltzmann method (LBM) formulation. The coupling between DEM and LBM is achieved through a new algorithm based on a volume-fraction approach to consider three-dimensional convex polyhedral particles moving through fluid. The algorithm establishes the interaction using linear programming and simplex integration and is validated against experimental data. This approach to modeling the interaction between complex polyhedral particles and fluid is shown to be accurate for directly simulating hydrodynamic forces on the particles.     

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.     

颗粒材料数值样本的坐标排序生成技术

颗粒材料离散颗粒模型的数值模拟结果与颗粒材料的数值样本密切相关,随着离散单元在颗粒材料数值模拟领域的广泛应用,颗粒材料的数值样本生成技术日益受到重视。基于RSA模型研究如何使随机生成的颗粒材料更密实,对均匀颗粒而言亦即如何在指定区域内生成更多的颗粒,讨论了4类修正方案,并建议了一种基于坐标排序的样本生成技术。研究表明,在传统的颗粒体随机生成技术基础上,通过对随机生成的x坐标序列或y坐标序列进行排序,可使生成的颗粒材料数值样本更密实。     

Frictional contact algorithms in SPH for the simulation of soil–structure interaction

Simulation of frictional contact between soils and rigid or deformable structure in the framework of smoothed particle hydrodynamics (SPH) is presented in this study. Two algorithms are implemented into the SPH code to describe contact behavior, where the contact forces are calculated using the law of conservation of momentum based on ideal plastic collision or using the criteria of partial penetrating. In both algorithms, the problem of boundary deficiency inherited from SPH is properly handled so that the particles located at contact boundary can have precise acceleration, which is critical for contact detection. And the movement and rotation of the rigid structure are taken into account so that it is easy to simulate the process of pile driving or movement of a retaining wall in geotechnical engineering analysis. Furthermore, the capability of modeling deformability of a structure during frictional contact simulations broadens the fields of SPH application. In contrast to previous work dealing with contact in SPH, which usually use particle‐to‐particle contact or ignoring sliding between particles and solid structure, the method proposed here is more efficient and accurate, and it is suitable to simulate interaction between soft materials and rigid or deformable structures, which are very common in geotechnical engineering. A number of numerical tests are carried out to verify the accuracy and stability of the proposed algorithms, and their results are compared with analytical solutions or results from finite element method analysis. Good agreement obtained from these comparisons suggests that the proposed algorithms are robust and can be applied to extend the capability of SPH in solving geotechnical problems. Copyright © 2013 John Wiley & Sons, Ltd.     

基于GIS空间数据的滑坡SPH粒子模型研究

针对光滑粒子流体动力学方法（SPH）在滑坡模拟中建立粒子模型的难题,提出了基于地理信息系统（GIS）栅格数据的粒子排列与插入方法。根据该方法,建立了滑坡SPH粒子模型及相关粒子生成程序,进一步以结合摩尔-库仑破坏准则的SPH宾汉流体模型为核心,实现了运用SPH方法模拟滑坡破坏后三维运动的过程。该SPH模型在对唐家山滑坡的模拟中得到了验证,并预测了金坪子滑坡破坏后的影响范围。结果表明：基于GIS空间数据的滑坡SPH粒子模型具有可行性与良好的适用性。以GIS数据库为基础,开展滑坡灾害的模拟研究,将大大提高对滑坡等地质灾害的仿真分析,为滑坡灾害的预测与防治提供参考。     

Study on the particle breakage of ballast based on a GPU accelerated discrete element method

Breakage of particles will have greatly influence on mechanical behavior of granular material(GM)under external loads,such as ballast,rockfill and sand.The discrete element method(DEM)is one of the most popular methods for simulating GM as each particle is represented on its own.To study breakage mechanism of particle breakage,a cohesive contact mode is developed based on the GPU accelerated DEM code-Blaze-DEM.A database of the 3D geometry model of rock blocks is established based on the 3D scanning method.And an agglomerate describing the rock block with a series of non-overlapping spherical particles is used to build the DEM numerical model of a railway ballast sample,which is used to the DEM oedometric test to study the particles’breakage characteristics of the sample under external load.Furthermore,to obtain the meso-mechanical parameters used in DEM,a black-analysis method is used based on the laboratory tests of the rock sample.Based on the DEM numerical tests,the particle breakage process and mechanisms of the railway ballast are studied.All results show that the developed code can better used for large scale simulation of the particle breakage analysis of granular material.     

Lagrangian meshfree particles method (SPH) for large deformation and failure flows of geomaterial using elastic–plastic soil constitutive model

Simulation of large deformation and post‐failure of geomaterial in the framework of smoothed particle hydrodynamics (SPH) are presented in this study. The Drucker–Prager model with associated and non‐associated plastic flow rules is implemented into the SPH code to describe elastic–plastic soil behavior. In contrast to previous work on SPH for solids, where the hydrostatic pressure is often estimated from density by an equation of state, this study proposes to calculate the hydrostatic pressure of soil directly from constitutive models. Results obtained in this paper show that the original SPH method, which has been successfully applied to a vast range of problems, is unable to directly solve elastic–plastic flows of soil because of the so‐called SPH tensile instability. This numerical instability may result in unrealistic fracture and particles clustering in SPH simulation. For non‐cohesive soil, the instability is not serious and can be completely removed by using a tension cracking treatment from soil constitutive model and thereby give realistic soil behavior. However, the serious tensile instability that is found in SPH application for cohesive soil requires a special treatment to overcome this problem. In this paper, an artificial stress method is applied to remove the SPH numerical instability in cohesive soil. A number of numerical tests are carried out to check the capability of SPH in the current application. Numerical results are then compared with experimental and finite element method solutions. The good agreement obtained from these comparisons suggests that SPH can be extended to general geotechnical problems. Copyright © 2008 John Wiley & Sons, Ltd.     

Study of mechanical behavior and strain localization of methane hydrate bearing sediments with different saturations by a new DEM model

This paper presents a numerical investigation into mechanical behavior and strain localization in methane hydrate (MH) bearing sediments using the distinct element method (DEM). Based on the results of a series of laboratory tests on the bonded granules idealized by two glued aluminum rods and the available experimental data of methane hydrate samples, a pressure and temperature dependent bond contact model was proposed and implemented into a two-dimensional (2D) DEM code. This 2D DEM code was then used to numerically carry out a series of biaxial compression tests on the MH samples with different methane hydrate saturations, whose results were then compared with the experimental data obtained by Masui et al. [9]. In addition, stress, strain, void ratio and velocity fields, the distributions of bond breakage and averaged pure rotation rate (APR) as well as the evolution of strain localization were examined to investigate the relationships between micromechanical variables and macromechanical responses in the DEM MH samples. The numerical results show that: (1) the shear strength increases as methane hydrate saturation S_MH increases, which is in good agreement with the experimental observation; (2) the strain localization in all the DEM MH samples develops with onset of inhomogeneity of void ratio, velocity, strain, APR, and distortion of stress fields and contact force chains; and (3) the methane hydrate saturation affects the type of strain localization, with one shear band developed in the case of 40.9% and 67.8% methane saturation samples, and two shear bands formed for 50.1% methane saturation sample.     

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.     

Multiphase SPH modeling of free surface flow in porous media with variable porosity

This paper presents a numerical model for simulating free surface flow in porous media with spatially varying porosity. The governing equations are based on the mixture theory. The resistance forces between solid and fluid is assumed to be nonlinear. A multiphase SPH approach is presented to solve the governing equations. In the multiphase SPH, water is modeled as a weakly compressible fluid, and solid phase is discretized by fixed solid particles carrying information of porosity. The model is validated by several numerical examples including seepage through specimen, fast flow through rockfill dam and wave interaction with porous structure. Good agreements between numerical results and experimental data are obtained in terms of flow rate and evolution of free surface. Parameter study shows that (1) the nonlinear resistance law provides more accurate results; (2) particle size and porosity have significant influence on the porous flow.