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.     

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

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

颗粒离散单元法动力人工边界设置方法

采用颗粒离散单元法进行动力计算时,人工截断边界上需设置吸收边界条件,以防止波的反射。鉴于颗粒离散单元数值计算模型的人工边界上颗粒单元半径大小不一、边界面凸凹不平,在连续介质的黏性、黏弹性、自由场边界条件方程基础之上,推导出适用于离散介质的等效方程。在离散介质的黏性边界条件等效方程中引入微调系数,提出比值迭代法以快速确定其最优值,以实现对波的最佳吸收。采用二维颗粒离散单元计算软件PFC2D,分别建立黏性、黏弹性、自由场边界条件相关数值分析模型,探讨颗粒分布模式对黏性边界上颗粒单元半径、速度分布及比值迭代过程的影响;采用外源波动算例及经典Lamb问题算例验证黏弹性边界设置方法的正确性;通过隧洞算例检验提出的自由场边界条件设置方法的正确性。     

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 numerical study of particle motion and two‐phase interaction in aeolian sand transport using a coupled large eddy simulation – discrete element method

Wind‐blown sand movement, considered as a particle‐laden two‐phase flow, was simulated by a new numerical code developed in the present study. The discrete element method was employed to model the contact force between sand particles. Large eddy simulation was used to solve the turbulent atmospheric boundary layer. Motions of sand particles were traced in the Lagrangian frame. Within the near‐surface region of the atmospheric boundary layer, interparticle collisions will significantly alter the velocity of sand. The sand phase is quite dense in this region, and its feedback force on fluid motion cannot be ignored. By considering the interparticle collision and two‐phase interaction, four‐way coupling was achieved in the numerical code. Profiles of sand velocity from the simulations were in good agreement with experimental measurements. The mass flux shows an exponential decay and is comparable to reported experimental and field measurements. The turbulence intensities and shear stress of sand particles were estimated from particle root‐mean‐square velocities. Distributions of slip velocity and feedback force were analysed to reveal the interactions between sand particles and the continuous fluid phase.     

Transient seepage analysis in zoned anisotropic soils based on the scaled boundary finite‐element method

The scaled boundary finite‐element method, a semi‐analytical computational scheme primarily developed for dynamic stiffness of unbounded domains, is applied to the analysis of unsteady seepage flow problems. This method is based on the finite‐element technology and gains the advantages of the boundary element method as well. Only boundary of the domain is discretized, no fundamental solution is required and singularity problems can be modeled rigorously. Anisotropic and non‐homogeneous materials satisfying similarity are modeled with no additional efforts. In this study, firstly, formulation of the method for the transient seepage flow problems is derived followed by its solution procedures. The accuracy, simplicity and applicability of the method are demonstrated via four numerical examples of transient seepage flow – three of them are available in the literature. Homogenous, non‐homogenous, isotropic and anisotropic material properties are considered to show the versatility of the technique. Excellent agreement with the finite‐element method is observed. The method out‐performs the finite‐element method in modeling singularity points. Copyright © 2014 John Wiley & Sons, Ltd.     

基于区域分解法的地下水有限元与边界元耦合模型—淄博市王旺庄水源地地下水数值模拟

区域分解法（DDM）是20世纪90年代兴起的一种求解偏微分方程的新方法，方法本身独到的耦合思想和高效的并行计算机理，对于求解复杂的、大型的地下水问题具有相当的优势和广阔的应用前景。本文以淄博市王旺庄水源地地下水流模型为例，应用重叠型区域分解法（DDM）构造了边界单元法（BEM）与有限单元法（FEM）耦合模型，在两种数值方法各自优点的基础上，更形象地再现了实际水文地质原型，有效地消除了人为边界造成的流场失真。     

Sensitivity of incipient particle motion to fluid flow penetration depth within a packed bed

A discrete element method is applied to a three‐dimensional analysis related to sediment entrainment on a micro‐scale. Sediment entrainment is the process by which a fluid medium accelerates particles from rest and advects them upward until they are either transported as bedload or suspended by the flow. Modelling of the entrainment process is a critically important aspect for studies of erosion, pollutant resuspension and transport, and formation of bedforms in environmental flows. Previous discrete element method studies of sediment entrainment have assumed the flow within the particle bed to be negligible and have only allowed for the motion of the topmost particles. At the same time, micro‐scale experimental studies indicate that there is a small slip of the fluid flow at the top of the bed, indicating the presence of non‐vanishing fluid velocity within the topmost bed layers. The current study demonstrates that the onset of particle incipient motion, which immediately precedes particle entrainment, is highly sensitive to this small fluid flow within the topmost bed layers. Using an exponential decay profile for the inner‐bed fluid flow, the discrete element method calculations are repeated with different fluid penetration depths within the bed for several small particle Reynolds numbers. For cases with slip velocity corresponding to that observed in previous experiments with natural sediment, the predicted particle velocity is found to be a few percent of the fluid velocity at the top of the viscous wall layer, which is a reasonable range of velocities for observation of incipient particle motion. This method for prescribing the fluid flow within the particle bed allows for the current discrete element method to be extended in future studies to the analysis of sediment entrainment under the influence of events such as turbulent bursting. Additionally, predictions for the slip velocities and fluid flow profile within the bed suggest the need for further experimental studies to provide the data necessary for additional improvement of the discrete element method models.     

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.     

Development of a discontinuous approach for modeling fluid flow in heterogeneous media using the numerical manifold method

In the numerical modeling of fluid flow in heterogeneous geological media, large material contrasts associated with complexly intersected material interfaces are challenging, not only related to mesh discretization but also for the accurate realization of the corresponding boundary constraints. To address these challenges, we developed a discontinuous approach for modeling fluid flow in heterogeneous media using the numerical manifold method (NMM) and the Lagrange multiplier method (LMM) for modeling boundary constraints. The advantages of NMM include meshing efficiency with fixed mathematical grids (covers), the convenience of increasing the approximation precision, and the high integration precision provided by simplex integration. In this discontinuous approach, the elements intersected by material interfaces are divided into different elements and linked together using the LMM. We derive and compare different forms of LMMs and arrive at a new LMM that is efficient in terms of not requiring additional Lagrange multiplier topology, yet stringently derived by physical principles, and accurate in numerical performance. To demonstrate the accuracy and efficiency of the NMM with the developed LMM for boundary constraints, we simulate a number of verification and demonstration examples, involving a Dirichlet boundary condition and dense and intersected material interfaces. Last, we applied the developed model for modeling fluid flow in heterogeneous media with several material zones containing a fault and an opening. We show that the developed discontinuous approach is very suitable for modeling fluid flow in strongly heterogeneous media with good accuracy for large material contrasts, complex Dirichlet boundary conditions, or complexly intersected material interfaces. Copyright © 2015 John Wiley & Sons, Ltd.     

SPH-DEM coupling method based on GPU and its application to the landslide tsunami. Part I: method and validation

Landslide-induced tsunami is a complex fluid–solid coupling process that plays a crucial role in the study of a disaster chain. To simulate the coupling behaviors between the fluid and solid, a graphics processing unit-based coupled smoothed particle hydrodynamics (SPH)-discrete element method (DEM) code is developed. A series of numerical tests, which are based on the laboratory test by Koshizuka et al. (Particle method for calculating splashing of incompressible viscous fluid, 1995) and Kleefsman et al. (J Comput Phys 206:363–393, 2005), are carried out to study the influence of the parameters, and to verify the accuracy of the developed SPH code. To ensure accurate results of the SPH simulation, the values for the diffusion term, particle resolution (1/25 characteristic length), and smoothing length (1.2 times of particle interval) are suggested. The ratio of the SPH particle size and the DEM particle’s diameter influences the accuracy of the coupling simulation between solid particles and water. For the coupling simulation of a single particle or a loose particle assembly (not contact each other) with fluid, this ratio should be smaller than 1/20; for a dense particle assembly, a ratio of smaller than 1/6 will be good.

     

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.     

基于B-DEM的颗粒聚合体应力计算和破碎路径模拟

结合边界元法和离散元法,提出一种可以进行计算颗粒内部应力和破碎路径的方法。该方法利用离散元法求解颗粒的相互作用和每个颗粒上的荷载。然后利用边界元法计算颗粒的应力分布,为了实现动态平衡,将颗粒的加速度视为恒定大小的体力。但体力导致边界积分方程中出现域积分,故采用直线积分法将域积分转化为边界积分,以保证边界元法降维的优势。为了提高边界元的计算效率,对于几何形状相似的颗粒,以其中一个颗粒作为模板颗粒,只需要计算模板颗粒在局部坐标系中的系数矩阵,其他相似颗粒可以通过局部和全局坐标系之间的映射获得。在得到应力后,基于Hoek-Brown准则来判断颗粒是否破碎。此外,将破坏路径简化为直线,并采用最小二乘法拟合得到破坏路径。     

基于边界元法的边坡矢量和稳定分析

矢量和法物理力学意义明确,计算简单,且能根据边坡当前的应力分布状态合理地评价其整体稳定性状态。其中边坡的应力状态通常是采用有限元法来求解。由于边界元法具有研究问题降阶、离散化带来的误差值仅产生在边界以及计算量小等优点,在工程中得到了广泛应用;对于平面问题,以源点作为原点,以所积分单元的切向和法向为坐标轴建立局部坐标系,对于线性单元可以得到所有积分的解析解。因此,可以得到计算区域内部任意点的场变量的解析解,这就保证了位于边界附近区域场变量的精度。利用边界元法得到二维边坡体内连续的应力分布状态,使用矢量和法对该边坡进行稳定性分析,并且与基于有限元的矢量和法、极限平衡法进行对比分析。边坡圆弧滑面和折线滑面的计算结果表明,基于边界元法得到的矢量和安全系数和基于有限元的矢量和法、极限平衡法基本一致;边界元法对应的矢量和安全系数对边界单元尺寸不敏感。     

Numerical simulation of particle plugging in hydraulic fracture by element partition method

Hydraulic fracturing (HF) treatment often involves particle migration and is applied for propping or plugging fractures. Particle migration behaviors, e.g., bridging, packing, and plugging, significantly affect the HF process. Hence, it is crucial to effectively simulate particle migration. In this study, a new numerical approach is developed based on a coupled element partition method (EPM). The EPM is used to model natural and hydraulic fractures, in which a fracture is allowed to propagate across an element, thereby avoiding remeshing in fracture simulations. To characterize the water flow process in a fracture, a fully hydromechanical coupled equation is adopted in the EPM. To model particle transportation in fractures with water flow, each particle is treated as a discrete element. The particles move in the fracture as a result of being dragged by fluid. Their movement, contact, and packing behaviors are simulated using the discrete element method. To reflect the plugging effect, an equivalent aperture approach is proposed. Using this method, the particle migration and its effect on water flow are well simulated. The simulation results show that this method can effectively reproduce particle bridging, plugging, and unblocking in a hydraulic fracture. Furthermore, it is demonstrated that particle plugging significantly affects water flow in a fracture and hence the propagation of hydraulic fracture. This method provides a simple and feasible approach for the simulation of particle migration in a hydraulic fracture.     

Time-domain dynamic analysis of dam–reservoir–foundation interaction including the reservoir bottom absorption

In this paper, time-domain dynamic analysis of dam–reservoir–foundation interaction is presented by coupling the dual reciprocity boundary element method (DRBEM) in the infinite reservoir and foundation domain and the finite element method in the finite dam domain. An efficient coupling procedure is formulated by using the sub-structuring method. The effects of the reservoir bottom absorption are included in the formulations. Sharan's boundary condition for the far-end of the infinite fluid domain is implemented. To verify the proposed scheme, numerical examples are carried out by comparing with the available exact solutions and finite–finite element coupling results for the dam–reservoir interaction. A complete dam–reservoir–foundation interaction is also studied by including the bottom absorption effects. Copyright © 2004 John Wiley & Sons, Ltd.     

用空气单元法求解渗流场的逸出边界问题

依据空气单元模拟排水孔的思路,在渗流场任意可能逸出边界外布置一层渗透性较强的单元,与其他实体单元一起进行渗流计算,可以实现该边界的自动求解。与现行方法相比,不仅省去了逸出边界迭代计算的麻烦,而且避免了由于逸出点定位不准确而可能引起的渗流计算不合理或收敛困难等缺陷。算例分析表明,空气单元法的计算精度与空气单元的相对渗透系数R有关,与厚度L无关。当选取了合适的R值后,可以较好地逼近现行方法。对于因引入空气单元而造成的强弱透水介质界面处数值拟真性差的现象,还提出了加入接触面薄层单元的改进方法。方法改进后,计算精度可进一步提高。最后,以小湾水电站22#坝段渗流场计算为例,说明该法的实用性和可行性。     

A new development of the scaled boundary finite element method for wave motion in layered half-space

The scaled boundary finite element method (SBFEM) developed by Wolf and Song has shown certain parallels to the finite element method (FEM) and boundary element method (BEM). Because of its semi-analytical nature, SBFEM is particularly suitable for the analysis of wave propagation in unbounded domains. This paper makes a certain modification of the standard SBFEM. A new idea of scaling surface instead of a scaling center is introduced to formulate the governing SBFE equations for the analysis of wave propagation in multilayered half-space, which leads to simplifying the modeling and saving considerably the computational effort. In addition, by employing the proposed approach, some problems encountered in engineering practice, which are difficult to deal with by the conventional SBFEM, for example, 3D foundation impedance on half-space with irregular geographical features, can be effectively solved. The proposed approach also helps to simplify the solution of shell structures. Numerical examples are provided to validate the accuracy and efficiency of the proposed approach.     

A novel coupled FDM-DEM modelling method for flexible membrane boundary in laboratory tests

Difficulties are involved in discrete element method (DEM) modelling of the flexible boundary, that is, the membranes covering the soil sample, which can be commonly found in contemporary laboratory soil tests. In this paper, a novel method is proposed wherein the finite difference method (FDM) and DEM are coupled to simulate the rubber membrane and soil body, respectively. Numerical plane strain and triaxial tests, served by the flexible membrane, are implemented and analysed later. The effect of the membrane modulus on the measurement accuracy is considered, with analytical formulae derived to judge the significance of this effect. Based on an analysis of stress-strain responses and the grain rotation field, the mechanical performances produced by the flexible and rigid lateral boundaries are compared for the plane strain test. The results show that (1) the effect of the membrane on the test result becomes more significant at larger strain level because the membrane applies additional lateral confining pressure to the soil body; (2) the tested models reproduce typical stress and volumetric paths for specimens with shear bands; (3) for the plane strain test, the rigid lateral boundary derives a much higher peak strength and larger bulk dilatation, but a similar residual strength, compared with the flexible boundary. The latter produces a more uniform (or ‘diffuse') rotation field and more mobilised local kinematics than does the former. All simulations show that the proposed FDM-DEM coupling method is able to simulate laboratory tests with a flexible boundary membrane.     

Hydro-micromechanical modeling of wave propagation in saturated granular crystals

Biot theory predicts wave velocities in a saturated granular medium using the pore geometry, viscosity, densities, and elastic moduli of the solid skeleton and pore fluid, neglecting the interaction between constituent particles and local flow, which becomes essential as the wavelength decreases. Here, a hydro-micromechanical model, for direct numerical simulations of wave propagation in saturated granular media, is implemented by two-way coupling the lattice Boltzmann method (LBM) and the discrete element method (DEM), which resolve the pore-scale hydrodynamics and intergranular behavior, respectively. The coupling scheme is benchmarked with the terminal velocity of a single sphere settling in a fluid. In order to mimic a small amplitude pressure wave entering a saturated granular medium, an oscillating pressure boundary on the fluid is implemented and benchmarked with the one-dimensional wave equation. The effects of input waveforms and frequencies on the dispersion relations in 3D saturated poroelastic media are investigated with granular face-centered-cubic crystals. Finally, the pressure and shear wave velocities predicted by the numerical model at various effective confining pressures are found to be in excellent agreement with Biot analytical solutions, including his prediction for slow compressional waves.