Research Article  Open Access
QingYun Zeng, JianPing Pan, HanZheng Sun, "SPH Simulation of Structures Impacted by Tailing Debris Flow and Its Application to the Buffering Effect Analysis of Debris Checking Dams", Mathematical Problems in Engineering, vol. 2020, Article ID 9304921, 17 pages, 2020. https://doi.org/10.1155/2020/9304921
SPH Simulation of Structures Impacted by Tailing Debris Flow and Its Application to the Buffering Effect Analysis of Debris Checking Dams
Abstract
Since the tailing dam fails catastrophically with substantial instantaneous deformation, it is difficult to measure the migration of debris flow caused by the failure of the tailings dam. A simulation model of tailing debris flow based on Smoothed Particle Hydrodynamics (SPH) theory of elasticplastic constitutive equation has been established by considering the viscoplasticity of mud and the elasticplastic characteristics of tailing sand to investigate the impact effect of tailing flow on the downstream structures. By comparing the experimental and two different simulation results obtained, it can be concluded that SPH elasticplastic constitutive model can effectively simulate the accumulation and migration processes of the tailing debris flow, which indicates that the SPH model has good applicability to solve geotechnical large deformation problems of similar tailings flow slide. Then, the verified simulation model developed based on a series of simulations of tailing debris flow propagations was used to determine the momentum reduction on the downstream structure resulting from the presence of a simple checking dam perpendicular to the direction of propagation and to determine the characteristics of stresses applied to this structure in terms of peak impact force and evolution over time to the main flow direction.
1. Introduction
Tailings dam is one of the most complex geotechnical structures. As the rainfall increases, the tailing infiltration line in the reservoir increases, and once the saturation is reached, dam break occurs in the form of highspeed moving sandbearing mudflow [1, 2]. These flows commonly occur without a prior warning in mountainous regions around the world and therefore cause massive loss of both human lives and properties. Debris flow caused due to the collapse of saturated tailing soil is usually difficult to solve in fluid dynamics. The relevant parameters are not only difficult to be measured in the field, but the calculation is also pretty complicated. Therefore, the numerical simulation calculation has become the first choice to study the tailing debris flow. The literature shows that most of the existing simulation methods applied to viscous debris flow are primarily gridbased finite element method (FEM) [3] and Boltzmann algorithm [4]. FEM is bounded by the element mesh to deal with large soil deformation. FEM has cell grid constraints that either interrupt the calculation or provide the inaccurate result of the calculation [5]. Similarly, the Boltzmann algorithm often requires a perfect grid to solve the higherorder differential equations, which greatly increases the computational cost [6].
Nowadays, the meshless method has been developed to overcome the deficiencies of the abovestated calculation method applied in the mechanical behavior of geomaterial analysis. The whole solution region in the meshless method is discretized into independent nodes and has no dependence on the grid, and the calculation methods under the different framework of the meshless theory are established, such as DEM (discrete element method), EFG (elementfree galerkin method, improved DEM), and SPH (smoothed particle hydrodynamics method), especially the DEM and EFG methods show a highly computational accuracy in the application of rock displacement and deformation in rock mechanics or geotechnical mechanics [7–10]. However, the applicability of methods such as DEM and EFG in large deformation of soil flow problems of similar debris flow or landslide have some defects, because these methods are not easy to determine some nonphysical parameters such as the friction and elastic parameters of the contact zone which are to be calibrated as they are not directly obtained from experiments and they are difficult to directly describe the stressstrain relationship of soil, which limited their application in this research direction [11–13]. The saturated tailing soil mainly experiences the stage of soil stiffness degradation and the stage of tailing debris flow movement in the process of large deformation of soil flow. When the saturated soil has been transformed from solid to liquid, it will be difficult reasonably to calculate and predict the complex dynamic behavior of soil by using traditional solid mechanics [13]. Therefore, it is necessary to introduce elasticplastic mechanics and fluid dynamics to analyze the different phase states in the process of saturated tailing soil deformation.
Recently, a meshless particle method known as SPH has attracted more attention. The SPH method can get rid of the finite element method's dependence on the grid and shows obvious advantages in the problems related to grid distortion and grid movement [14, 15]. On the other hand, in dealing with the problems related to dynamic response which occurs during large deformation of debris flow, SPH can easily introduce constitutive equations describing the properties of various materials including solid mechanics and hydrodynamics, and the advantages of Lagrangian particularization will describe the moving state of a solidfluid mixture in the form of particles, which is beneficial to accurately grasp the deformation characteristics of different positions in the process of large deformation of debris flow [16, 17]. Then, several researchers had applied SPH to the migration impact simulation of viscous debris flow. Antonio et al. [18] have successfully established the simulation calculation model of viscous debris flow by referring to the test parameters of indoor debris flow sliding experiments conducted by Laigle and Coussot [19]. Huang et al. [20] utilized SPH to carry out the migration rule analysis of landslide and debris flow disaster under the impact of the Wenchuan earthquake. Dai et al. [21] established the coupled SPH model to simulate the mountain debris flowstructure interaction. However, the SPH method is not widely used in the case of nonNewtonian flow such as tailing debris flow. It can be concluded from the related literature about the SPH control equation that the characteristics of the stressstrain relation of tailings debris flow are rarely reported. As a result, many researchers pay attention to the research and application of debris flow constitutive model based on SPH theory.
The present paper focused on the application of the SPH method to the simulations of structures impacted by tailing debris flow that involves large deformations of soil flow. Firstly, the main features of the SPH model were presented with specific attention considered to the treatment of boundary conditions and rheological behavior and SPH fluidsolid coupled model based on the constitution of elasticplastic constitutive mechanics for tailings debris flow in Section 2 was established. Then, in Section 3, the model is validated by simulating the impact effect of tailing flow on the downstream structures with the indoor model test and two different numerical simulations (SPH and FEM) to check the applicability of SPH in simulating the tailing debris flow. In Section 4, the verified SPH simulation model has been applied to derive the respective relationship between the height of checking dam and height of mud surface, impact force, and momentum reduction of the tailing debris flow that occurs at the downstream structure. The results obtained from the calculations were then introduced and discussed in detail. Finally, several conclusions were summarized in Section 5.
2. SPH Formulations of Structures Impacted by Tailing Debris Flow
To get rid of the grid constraint caused by the mesh division in the FEM model, the SPH method uses the free distribution of discrete particles that carry field variables instead of its problem domain. The aforementioned method also introduces the kernel function to simplify the complex interactive relationship between particles [22]. Therefore, tailing debris flow problems can be solved by tracking the movements of the particles and the changes of the carried field variables.
2.1. SPH Integral Interpolation Function Approximation
For continuous fluid motion, the SPH integral interpolation function f (x) can be accurately simulated by a large number of particles containing mass and momentum [23]. The function f (x) at position x (x, y, z) is approximated by the following expression:where Ω represents the computational domain at position x, is the element of volume over the domain, is the weight equation or smooth kernel function, and h is smooth length, which determines the size of the computational domain of the kernel function .
Since the distribution of the initial interpolation points is not random, it changes constantly with the motion of the fluid after the initial value is taken. To describe the fluid motion more accurately, (1) is discretized into particle form, which is given by the following equation:where i represents the concerning particle, j is a neighboring particle in Ω, and and represent the mass and density of particle j, respectively. Smooth length h is 1.05–1.5 times of particle spacing , which determines the number of influenced particles in the computational domain (Figure 1).
The threedimensional form of cubic spline kernel function is adopted in the following form [16], where relative displacement vector , and :
The smooth length h should satisfy two requirements [24]. First, the value of h should be smaller than the influence radius of the individual particles, so that less number of iterations is required for interpolation calculation. Second, h must be larger than the particle spacing to ensure the existence of correlation particles in the range of influence; otherwise, the lack of correlation between particles might lead to the particle penetration of the finite element mesh in the coupling calculation of SPH and FEM. The purpose of applying the smooth length is to ensure that the total mass of the nearest neighbor particles is constant or maintained within an allowable range. According to the discussion on the relationship between smooth length and particle spacing, the rational range of smooth length h should be [25]. Based on these requirements and the cost of calculation, the economical and reasonable value of has been selected for the present research work.
2.2. Governing Equations of Viscous Debris Flow
When the tailing dam runs at a high water level, the tailings are usually close to saturation. Since the shear strength of the tailings dam is very low, the sliding surface appears easily, which in turn forms the dam breaking debris flow. The tailings debris flow satisfies the characteristic description of high sediment flow in terms of mechanical properties and movement mechanism, showing obvious incompressible viscous nonNewtonian flow characteristics [26]. The most widely used nonNewtonian models for debris flow simulations are the Bingham [27, 28] and HerschelBulkley [18, 29] models.
In the present work, referring to the study of rheological control of mountain debris flow [20, 30, 31], the viscosity between tailing particles in highspeed motion is very small in comparison to the gravity of the tailing fluid and the viscosity resistance between slurry; therefore, it can be neglected. At the same time, the mixture of tailing particles and rainwater has been considered as a uniform continuous incompressible Bingham fluid. The governing equations of viscous debris flow for the Bingham model are the continuity equation and momentum equation. The following equation expresses the application of the aforementioned equations (continuity and momentum equation) in the present research work:where is the velocity vector, represents a divergence symbol, is the density of debris flow, is hydrostatic pressure, is the total stress tensor, is the interface curvature, represents the viscosity coefficient, is the gravitational acceleration vector, and is the distance function of the contact surface between air and mud.
To maintain a good correlation between the governing equations and the SPH kernel function, governing equations are written in the following Lagrange form [23, 32]:where is the gravitational acceleration, is the velocity of debris flow, is the principal stress tensor, is the shear stress tensor, and represent the pressure on particles i and j, respectively, and and denote the Cartesian components x, y, and z with the Einstein convention applied to repeated indices.
To improve the stability of numerical calculations and avoid the shock of solidliquid coupling contact of dam break fluid, an artificial viscosity is added to the momentum equation which denotes dissipations of the kinetic energy of the shocked region into heat. Monaghan and Gingold [33] adopted this approach and introduced it in the following form:where is the relative distance, is the relative velocity, and and represent the mean sound velocity and average density, respectively. To avoid the overflow of the calculated results, the parameters are satisfied under the condition that the smooth length is dynamic. For the treatment of dam, break debris flows under viscoelastic characteristics, and are considered to be constant. In the current work, we have assumed q_{2} = 0.01–0.1 and q_{2} = 0 [34].
2.3. Constitutive Equations of Tailing Debris Flow
Debris flow, although possesses a complex composition, can be broadly considered as a mixture of water, clay, and granular matter. To simplify the calculation, the stressstrain relation of viscoplastic fluids such as the tailing debris flow is assumed to follow the Bingham model [21, 27, 28].
According to the characteristics of Bingham flow, the liquidphase stressstrain relation of tailing debris flow is described by the following expression:where is the shear stress tensor, is the viscosity coefficient, represents the yield strength, is the shear strain tensor, is the second invariant of the tensor of shear strain, and and denote the Cartesian components x, y, and z with the Einstein convention applied to repeated indices.
For ideal viscoplastic materials, the shear strain rate is given by the following equation:where of the particle i in the SPH form can be expressed as the following equation [35]:
Considering the solidliquid cohesive effect of saturated tailings, the molar Coulomb yield criterion is introduced to solve the shear stress:where is the hydrostatic pressure, is the internal friction angle, and is the cohesive strength.
For the solidphase part of debris flow, the stressstrain relation of elasticplastic material in solid mechanics principle is given by the following equation:where the rate of partial shear stress tensor can be modified according to Jaumann criterion [12] which is as follows:where is the shear modulus of tailings and represents the spin rate tensor. For ideal elastoplastic materials, the spin rate tensor must satisfy the following condition:where and of the particle i in SPH form can be written as follows [35]:
According to the variation of particle unit velocity in fluid mechanics [36], the relative velocity between two particles can be expressed as a matrix:where the square matrix in the upper formula can be decomposed into
The first matrix in the upper formula (17) is (antisymmetric), and and matrices are symmetric. These three matrices are also called secondorder tensors in hydrodynamics [36].
Among the nine components in the antisymmetric matrix , there are only three independent components which are as follows:
These three components are exactly the three components of the rotational angular velocity vector of the fluid particle, so is called the rotational velocity tensor. At the same time, , that is, half the curl of the velocity vector is satisfied (equation (14)).
Among the nine components of symmetric matrices and , there are only six independent components as follows:then and :where is the linear strain rate tensor of fluid particle in fluid mechanics in the direction of three axes and is the fluid particle’s angular deformation rate tensor. Therefore, the tensor is called the strain rate tensor in fluid mechanics, and is called the velocity gradient tensor.
Moreover, to ensure the weak compressibility and continuity of tailing debris flow, changes in hydrostatic pressure within the particle motion domain can be described directly by the equation of state (EOS) in the form of particle density. The SPH equation of state in the form of weak compression was proposed by Monaghan [37], which is as follows:where constant , is the initial average density of tailing debris flow, and is the sound velocity in the corresponding fluid medium. In the current work, c_{0} = 340 m/s. In the numerical model, it is assumed that the tailing is elastoplastic materials and the mud is weakly compressible viscoplastic fluid. The elasticplastic constitutive model of the tailings flow is finally established by using the equation of state (21) coupled with (8) and (12) to determine the contact correlation between the mud and the tailings (Figure 2).
2.4. Boundary Condition Used in Current SPH Model
While applying SPH, to save computation time, the boundary was assumed to be pressurefree or constituted of weakpressure interfaces. However, the idea of “nodetosurface” unidirectional contact control based on compensation by introducing the gradient of momentum and mass into the interpolation kernel function was chosen as it becomes the core idea of dealing with boundary problems [38]. Meanwhile, to avoid particle penetration at the boundary, an improved virtual particle compensation method was introduced. The simulation accuracy was improved substantially by setting virtual particles with no mass and no volume at the boundary to avoid the interference of real boundary particles and the moving fluid particles [39, 40].
While applying the virtual particle compensation method, a discrete layer of particles is set at the boundary where the displacement and mass of the boundary particles are assumed to be zero. According to Newton's third law, when the boundary particles come in contact with the fluid particles, repulsive forces are produced in the vertical direction. Consequently, the velocity vector and the displacement vector of the fluid particles change, and the fluid particles are separated and they do not penetrate, as presented in Figure 3. This force contribution from the boundary particles is added to the righthand side of the momentum equation ((6)).
2.5. Operational Framework of SPH Model
The SPH governing equations of the final elastoplastic constitutive model are obtained by substituting the SPH form variable ((8) and ((12)) from Section 2.3 into the momentum (6) and combining with continuity (5). According to the SPH governing equations, the flow chart of the SPH program is shown in Figure 4.
3. Verification of the Developed SPH Model
3.1. Experimental Setup for Tailing Debris Flow Model Test
Since the failure of tailing dams take place catastrophically, the relevant parameters are difficult to measure in the field, and therefore, the laboratory model test is one of the most widespread methods to research tailing debris flow. The tailing debris flow model test has been simulated using the fluidstructure coupled SPH model to verify its numerical accuracy. The schematic of the tailing flow modeling test is presented in Figure 5. The experiments were performed in a 4.3 m long × 0.35 m wide × 0.8 m highdischarge flume, in which the horizontal section length was 2.3 m and the slope horizontal projection length was 2.0 m. The stainless steel bottom plate of the flume was 0.35 m wide and two sidewalls were made of 0.25 m high clear acrylic plate to facilitate camera recording. Rigid checking dam with a height of 3.5 cm at 65 cm downstream of the liftgate and concrete structure of length, width, and height of 0.15 m, 0.15 m, and 0.55 m, respectively, were placed at Z = 2.2 m (note that for all the figures in which zaxis is displayed, the origin has been considered corresponding to the liftgate of tailing debris flow reservoir).
The weight of tailings and mud mixture of 30 kg with a volume concentration ratio of 1 : 3 were placed at the tailing debris flow reservoir with an initial slope towards the outflow boundary on the righthand side. The average particle size of tailing, d_{50} = 0.23 m, was obtained from a mine tailing pond in Jiangxi Province, China, and the physical and mechanical properties were estimated in the laboratory. Dynamic load cells and displacement sensors were set up at the dam and downstream structures, and the data were collected simultaneously under multichannels. The camera and slurry height measuring the scale were set up at the relevant places to record the slide and impact behavior of tailing flow.
3.2. SPH Simulation Model for Tailing Debris Flow Impact
In this section, the SPH model has been used to simulate the tailing flow impact downstream structure under the action of check dam blocking. The SPH model is consistent with the size of the indoor experimental model and the initial calculation domain of SPH shown in Figure 6. The bottom plate of the flume, checking dam, and structure are divided into FEM grids, and SPH particles assignment is used to describe the tailing flow. The debris flow around the tailing pond in the field consists of silt and heavy sand clay, and the behaviors of tailing flow are viscous matrix flow. According to the indoor rheological experiments and realtime debris flow observation in the field, the coefficient of viscosity varies from 0.5 to 2.0 Pa·s [41, 42]. So, the viscosity coefficient of tailing flow in the current work was chosen, μ = 1.0 Pa·s. Internal friction angle between tailing particles ϕ = 26.1 and cohesive strength C = 0.8 kPa. However, the internal friction angle between the mud and the tailing particles was ignored, and only the cohesive strength was considered. Accordingly, the spin parameters, displacement parameters, and other relevant initial parameters were set ((8) and (12)) to obtain the solution of momentum equation (6) and realize the position tracking of the tailing flow particles. To treat the coupling contact between FEM mesh and SPH particle, the former was set as a master segment and the latter was assigned to be the slave segment. Also, a dynamic friction coefficient of 0.12 between each other was assumed, and once the calculation begins, the tailing flow under the action of gravity will start to move from rest.
To carry out the FEM analysis, “solid” rigid material with density, Young’s modulus, and Poisson’s ratio of 7860 kg/m^{3}, 100 GPa, and 0.25, respectively, was selected. For the SPH validation model, the “viscoplastic” fluid materials were assumed to be comprised of about 9000 discretized SPH particles, and the particles were assumed to be uniformly distributed in the whole reservoir. There were 3000 particles, which represent tailing and are made up of “elasticplastic” material. After defining the initial density and particle arrangement of the material, the continuity equation (5) was solved to realize the mass tracking of tailing flow particles. It was observed from the results of SPH simulation of debris flow [18] that smaller the “MACH number” of controlling fluid compressibility, more accurate are the results of numerical calculation; therefore, MACH number less than 0.3 was recommended. To satisfy the weak compressibility of tailing flow, a pressure reduction coefficient (PC) = 0 was defined in the SPH model. The boundary constraint setting was carried out by applying a layer of virtual particles (about 2000) at X = 0.175 m plane and X = −0.175 m plane, respectively. The distance between the SPH particles was assumed to be 0.01 m. The parameters of the relevant SPH model are calculated and presented in Table 1.

3.3. Numerical Comparison with Physical Results
The comparison of the flow patterns of the dam impacted by tailings flow at 0.65 m downstream of the liftgate is presented in Figure 7. It was observed that the tested and the simulated results are in good agreement with each other. Further, it was observed from SPH simulation that the tailings with elasticplastic characteristics carry ascertain the degree of deposition, while the mud particles with viscoplastic characteristics carry some suspended tailings across the blocking device. Array flow characteristics such as “debris flow head” and “debris flow body” of mountain debris flow appeared in the model experiment as well as in the SPH and FEM simulations (Figures 7(b) and 7(c)). As the liftgate was opened, the tailing flow was rapidly flushed forward in the form of a lump of the watersand mixture, which accumulated strong kinetic energy of impact and formed “debris flow head” (Time = 2.5 s). When time = 2.8 s, the subsequent tailing flow movement became smooth, which is known as “debris flow body.” Besides, the above characteristics also confirmed that it is easy for the saturated tailing reservoir under longterm high water level operation to develop into viscous debris flow after unstable failure [43].
(a)
(b)
(c)
Meanwhile, a comparison between the CFX simulation results based on the FEM grid division and the SPH simulation results is shown in Figure 8, and the accumulation and migration processes of the tailing debris flow were consistent. It could be seen that the SPH method was feasible for simulating large deformation of fluidsolid coupling. Besides, in the view of the apparent particle spatter and the discriminating degree of fluid accumulation at the boundary of the solid wall, the SPH method showed better adaptability (Figures 7 and 8). It was noticed that, in the next frame where the time, t = 1.5 s, the “debris flow head” touches the checking dam, which results in a splash of the tailing flow particles and when the time, t = 2.5 s, the peak slurry height reaches 12 cm. The damage caused by the instantaneous impact of the “debris flow head” on the structure is significant, as evident in Figures 8 and 9. Also, the peak slurry height reaches 20 cm at a time, t = 4.5 s. However, increasing the time further to 5.25 s results in a reduction in the height of the slurry present in front of the structure and the particles gather and subsequently flatten.
(a)
(b)
(c)
(d)
(e)
(f)
The graph presented in Figure 10 depicts a comparison of the numerical and experimental time history of slurry surface elevation acting on checking dam and structure. It can be seen that the SPH and FEM numerical calculation results are consistent with the experimental data, especially at the time range of start to peak appearance. Although the liquid level entered the gentle stage, at time t = 4.5 s and time t = 7 s, the liquid level fluctuation occurred at the checking dam and structure, respectively. This behavior is also evident in the SPH simulation presented in Figures 8(d) and 8(f), respectively. Nevertheless, after the peak value appearance, FEM simulation results were larger than those of SPH and experimental results, and the liquid level fluctuation in FEM simulation was found to be insignificant. The SPH and experimental results have shown better consistency after the peak value appearance. Dai et al. [21] while studying the behavior of debris flow observed that the sudden liquid level fluctuations in the stationary phase are caused by the solidliquid coupling cohesive action between tailings and mud, which eventually results in large specific gravity tailing transport lag. So, it becomes imperative to complete the material composite setting of viscoplastic mud and elasticplastic tailings instead of considering the debris flow to be a homogeneous liquid substance.
By controlling the ratio of particle density to initial density, (21) was utilized to plot the pressure distribution curve. Meanwhile, the experimental pressure was measured by the load cells (electric resistance strain gauge) at the checking dam and the structure. It was noticed from the presented Figure 11 that the numerically calculated pressure was in good agreement with the experimental data, especially the SPH results showed better consistency of the peak values and more obvious characteristic of pressure fluctuation caused by liquid level fluctuation. Also, the instantaneous increase of the impact pressure in a short time is due to the strong destructive energy released at the “debris flow head.” Therefore, the tailing dam and the downstream structure should be separated by a safe distance so that the kinetic energy of the “debris flow head” does not exert a detrimental impact. From this section, it can be concluded that the solidliquid coupling SPH model can simulate the impact process of tailings flow slide which is similar to mudflow mountain migration characteristics having intermittent continuous flow characteristics and estimate the impact force on checking dam and downstream structure.
4. Numerical Simulation of Checking Dam Height Effect
When the tailing flow slide across the checking dam is stopped, then the intercepted sediment approximates a triangular wedge as presented in Figure 8(f). The tailing flow interception volume under the action of checking dam which is 3.5 cm in height accounts for about 30% of the tailing debris flow reservoir volume. Therefore, to reduce the impact of tailing debris flow on downstream structures, it is necessary to study the effect of checking dam height on debris flow disaster. In this context, the verified SPH solidliquid coupling simulation model can be used as presented in Figure 6 to modify the height of checking dam h from 0 cm, 2 cm, 3.5 cm, and 5 cm to 6.5 cm at 65 cm downstream of the liftgate. The calculation parameters are consistent with Table 1.
4.1. Effect of Checking Dam Height on Debris Flow
The simulation results of the impact pattern of tailing flow under different checking dams heights are presented in Figure 12. At the time, t = 1.75 s, four groups of numerical simulation barring simulation 1 (h = 0) produce particles spatter of varying degrees under the action of checking the dam. With the increase in checking dam height, the impact of “debris flow head” is enhanced, and it also takes more time for tailing flow to pass through the dam. At the time, t = 4.9 s, under the action of checking dam and tailingmud cohesion, some of the tailings were accumulated near the dam.
4.2. RunUp and Impact Force
The influence of checking dam height exerted on both slurry runup and the impact force acting on the downstream structure have been presented in Figures 13 and 14. It was noted that the density field variable of tailing flow was obtained by solving the continuity equation (5). To assess the pressure distribution, the equation of state (21) was solved. Meanwhile, the area integral of pressure was performed to estimate the magnitude of impact force acting on the structure. Due to the buffer effect, while checking the dam, a slightly different spatial distribution of these particles could perhaps result in a different pressure determination at the front end of the downstream structure. The slurry runup peak on the structure, shown in Figure 13, is greater for h = 0 cm than it is for the other values of h. Consequently, the associated impact force is also greater for the condition of checking dam existent. With the increase in the checking dam height, the rate of the slurry runup in front of the downstream structure decreased, and a lag in the peak value was observed. Nevertheless, further investigation is required to properly elucidate this variability.
So, the influence of checking dam height on both slurry runup peak and the impact force peak acting on the checking dam and downstream structure have been investigated and presented in Figures 15 and 16. The blocking height has a significant buffering effect on the breaking tailing flow, and it varies linearly with dam height. Increasing the dam height by 1.5 cm, the peak value of slurry runup at the downstream structure is reduced by 15.95%, and the associated average reduction of the peak value of the impact force is 30.75%. On the other hand, the effective blocking efficiency at checking the dam also shows a linear distribution and the absolute slope value is consistent with the associated linear distribution at the downstream structure.
4.3. Velocity Distribution around the Structure
As the impact forces, presented in Figure 14 at the time, t = 4.55 s, for different checking dam heights correspond to key moments in the tailing flowstructure interaction, the velocity vectors in the xdirection are presented sequentially in Figure 17 to investigate the velocity distribution around the structure. When tailing debris flow impinges on the structure, the forward flow particles rapidly extrude and get divided into either side and the velocity of flow in the xdirection increases rapidly, which results in forming a nearly symmetrical “velocitydense zone” (shown in yellow and blue color) at ±45° corners in front of the structure. However, at the rear end of the structure, there is a slight presence of the “velocitydense zone” in the opposite direction, and an intermediate confluence is formed as the wake downstream from the structure develops. This occurs because the reverse flow rate component is produced by the reaction force of the boundary particles according to Newton’s third law [44]. The divided tailing flow is then reoriented towards the centerline of the flume, and since the X component of velocity is small as compared with the Z component, the divided tailing flow eventually combines further while moving downstream.
(a)
(b)
(c)
(d)
Figure 18 shows the comparison of numerically calculated pressure crosssectional view of tailing flowstructure interaction 0.01 m upstream of structure for time t = 4.55 s under the action of different checking dam heights. As the checking dam height increases, at the time, t = 4.55 s, the range of “velocitydense zone” decreases substantially. The divided tailing flow flowing around the structure reconnect later in the simulations is lagged (Figures 17(b)–17(d)) and less impact via kinetic energy is exerted on the side pillars of the structure for dam height h = 2 cm (Figures 17(a) and 18(a)). Although the pressure distribution before the structure is weakened by the checking dam, the concentrated distribution of pressure at side pillars is obvious (Figure 18)). From the point view of disaster prevention and control, the appearance of “velocitydense zone” and crosssectional pressure provides a certain reference basis for increasing the stiffness of both sides of the pillar in the structural design. Combining with the blocking dam set up upstream, the system can slow down the erosion of the side pillar foundation in the front of the structure by tailing flow and therefore avoid cracking and failure of the walls near the pillars.
(a)
(b)
(c)
(d)
5. Conclusions
In the present research work, the Bingham flow equation of motion and mudtailing constitutive model was modified, and smooth particle hydrodynamics (SPH) discrete formula of NS equation is established to evaluate the impact simulation of tailing flow on the downstream structure under the action of checking dam. The following conclusions can be drawn from the aforementioned study:(1)The impact process of the tailing flow showed the characteristics of “debris flows head” and “debris flows body” which are similar to the migration characteristics of mountain debris flow. The mud properties showed the formation of continuous flow, which has typical Bingham flow characteristics. The “debris flows head” exerted a substantial instantaneous impact on the structure; consequently, impact damage was easily caused to the downstream structure. To avoid the aforementioned problem, the tailing dam and the downstream structure should maintain a safe distance and meet the requirements of safety regulations.(2)The numerical calculation results, based on the SPH method, of impact process of the tailing flow, were in good agreement with the experimental results, especially the depth of mud inundation and the impact strength before the structure. The simulation program by using an SPH method based on the elasticplastic constitutive model had good applicability in solving geotechnical large deformation problems of similar tailing flow slides. The SPH algorithm can be used to dynamically reflect the disasterforming process, predict and evaluate the impact damage of the affected bodies under different disaster intensities, and lay a foundation for early warning and scientific guidance for disaster prevention and mitigation.(3)The SPH simulation model was used to analyze the buffer effect of checking dam height on the dam break tailing flow. The results showed that the blocking height had a linear distribution on the buffer effect of the tailing flow, and changes in the “velocitydense zone” in the range of ±45° on the front part of the structure directly affect the disasterbearing capacity of the side pillar foundation of the structure.
However, the SPH numerical model can further be developed better to employ the pore water pressure in tailing flowlike regime since the present SPH model is still controlled by hydrostatic pressure which ignores the influence of pore water pressure in tailingmud interaction. Besides, the detailed analysis of the viscosity coefficient and its effects on the tailings flow process is to be conducted in the further works.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Disclosure
JianPing Pan and HanZheng Sun are coauthors.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This research was supported by the National Natural Science Foundation of China (No. 51675059), Key Research and Development Project of Science and Technology Plan of Hunan Province, China (No. 2019SK2174), and Science and Technology Project of Education Department, Jiangxi Province, China (No. GJJ181337).
References
 G. Yin, G. Li, Z. Wei, L. Wan, G. Shui, and X. Jing, “Stability analysis of a copper tailings dam via laboratory model tests: a Chinese case study,” Minerals Engineering, vol. 24, no. 2, pp. 122–130, 2011. View at: Publisher Site  Google Scholar
 R. Aleixo, Y. Ozeren, and M. Altinakar, “PIVPTV measurements of a tailings dambreak flow,” Hydrodynamic and Mass Transport at Freshwater Aquatic Interfaces, Springer, Cham, Switzerland, 2016. View at: Publisher Site  Google Scholar
 M. Quecedo, M. Pastor, M. I. Herreros, and J. A. Fernández Merodo, “Numerical modelling of the propagation of fast landslides using the finite element method,” International Journal for Numerical Methods in Engineering, vol. 59, no. 6, pp. 755–794, 2004. View at: Publisher Site  Google Scholar
 L. Alessandro, K. W. Falk, M. Miller, and H. J. Herrmann, “LatticeBoltzmann method for geophysical plastic flows,” in Recent Advances in Modeling Landslides and Debris Flows, pp. 131–140, Springer, Cham, Switzerland, 2015, Springer Series in Geomechanics and Geoengineering. View at: Publisher Site  Google Scholar
 G. H. Shi, “Discontinuous deformation analysis: a new numerical model for the statics and dynamics of deformable block structures,” Engineering Computations, vol. 9, no. 2, pp. 157–168, 1992. View at: Publisher Site  Google Scholar
 F. J. Alexander, A. L. Garcia, and B. J. Alder, “A consistent Boltzmann algorithm,” Physical Review Letters, vol. 74, no. 26, pp. 5212–5215, 1995. View at: Publisher Site  Google Scholar
 M. Tapias, E. E. Alonso, and J. Gili, “A particle model for rockfill behaviour,” Géotechnique, vol. 65, no. 12, pp. 975–994, 2015. View at: Publisher Site  Google Scholar
 S. Arimoto and A. Murakami, “Characteristics of localized behavior of saturated soil by EFG,” in Proceedings of 1st International Workshop on New Frontiers in Computational. Geomechanics, pp. 169–172, World Scientific, Banff, Canada, 2002. View at: Google Scholar
 S. Soodeh and P Ali, “Threedimensional simulation of fully coupled hydromechanical behavior of saturated porous media using the element free galerkin (EFG) method,” Computers and Geotechnics, vol. 46, pp. 75–83, 2012. View at: Publisher Site  Google Scholar
 H. H. Bui and R. Fukagawa, “An improved SPH method for saturated soils and its application to investigate the mechanisms of embankment failure: case of hydrostatic porewater pressure,” International Journal for Numerical and Analytical Methods in Geomechanics, vol. 37, no. 1, pp. 31–50, 2013. View at: Publisher Site  Google Scholar
 Y. An, Q. Wu, C. Shi, and Q. Liu, “Threedimensional smoothedparticle hydrodynamics simulation of deformation characteristics in slope failure,” Géotechnique, vol. 66, no. 8, pp. 670–680, 2016. View at: Publisher Site  Google Scholar
 M. Hu, M. B. Liu, M. W. Xie, and G. R. Liu, “Threedimensional runout analysis and prediction of flowlike landslides using smoothed particle hydrodynamics,” Environmental Earth Sciences, vol. 73, no. 4, pp. 1629–1640, 2015. View at: Publisher Site  Google Scholar
 Y. Huang, L. Had, X. Pan, and Q Xu, “Numerical simulation of large deformation of soil flow based on the SPH method,” Chinese Journal of Geotechnical Engineering, vol. 31, no. 10, pp. 1520–1524, 2009, in Chinese. View at: Google Scholar
 N. Prapanch and T. Gaurav, “A study of energy transfer during water entry of solids using incompressible SPH simulations,” SadhanaAcademy Proceedings in Engineering Sciences, vol. 42, no. 4, pp. 517–531, 2017. View at: Publisher Site  Google Scholar
 C. Hermange, G. Oger, and T. D. Le, “Energy considerations in the SPH method with deformable boundaries and application to FSI problems,” Journal of Computational Physics, vol. 1, Article ID 100008, 2019. View at: Publisher Site  Google Scholar
 M. B. Liu and G. R. Liu, “Smoothed particle hydrodynamics (SPH): an overview and recent developments,” Archives of Computational Methods in Engineering, vol. 17, no. 1, pp. 25–76, 2010. View at: Publisher Site  Google Scholar
 M. A. Rahman and K. Konagai, “A handson approach to estimate debris flow velocity for rational mitigation of debris hazard,” Canadian Geotechnical Journal, vol. 55, no. 1, pp. 1–52, 2017. View at: Publisher Site  Google Scholar
 P. Antonio, M. Lorenzo, S. Nicola, and E. Paris, “SPH modeling of fast muddy debris flow: numerical and experimental comparison of certain commonly utilized approaches,” Italian Journal of Geosciences, vol. 132, no. 3, pp. 350–365, 2013. View at: Publisher Site  Google Scholar
 D. Laigle and P. Coussot, “Numerical modeling of mudflows,” Journal of Hydraulic Engineering, vol. 123, no. 7, pp. 617–623, 1997. View at: Publisher Site  Google Scholar
 Y. Huang, H. Cheng, Z. Dai et al., “SPHbased numerical simulation of catastrophic debris flows after the 2008 Wenchuan earthquake,” Bulletin of Engineering Geology and the Environment, vol. 74, no. 4, pp. 1137–1151, 2015. View at: Publisher Site  Google Scholar
 Z. Dai, Y. Huang, H. Cheng, and Q. Xu, “SPH model for fluidstructure interaction and its application to debris flow impact estimation,” Landslides, vol. 14, no. 3, pp. 917–928, 2017. View at: Publisher Site  Google Scholar
 R. L. Yaidel, R. Dirk, and R. M. Carlos, “Dynamic particle refinement in SPH: application to freesurface flow and noncohesive soil simulations,” Computational Mechanics, vol. 51, pp. 731–741, 2013. View at: Publisher Site  Google Scholar
 G. R. Liu and M. B. Liu, Smoothed Particle Hydrodynamics: A MeshFree Particle Method, World Scientific, Singapore, 2003.
 W. Shi, J. Chen, Y. Shen, and Y. Jiang, “Sph method with spacebased variable smoothing length and its applications to freesurface flow,” International Journal of Computational Methods, vol. 16, no. 2, pp. 640–657, 2019. View at: Publisher Site  Google Scholar
 J. Z. Xu and W. H. Tang, “Initial smoothing length and space between particles in SPH method for numerical simulation of highspeed impacts,” Chinese Journal of Computational Physics, vol. 26, no. 4, pp. 548–552, 2009, in Chinese. View at: Google Scholar
 L. T. Zhang, Q. Qi, Q. Li, and S. Zhang, “3D numerical simulation on the evolution of debris flow from dam break of tailing pond,” International Journal of Earth Sciences and Engineering, vol. 7, no. 1, pp. 111–115, 2014. View at: Google Scholar
 Y. Huang, W. Zhang, Q. Xu, P. Xie, and L. Xie, “Runout analysis of flowlike landslides triggered by the Ms 8.0 2008 Wenchuan earthquake using smoothed particle hydrodynamics,” Landslides, vol. 9, no. 2, pp. 275–283, 2012. View at: Publisher Site  Google Scholar
 S. Hao, T. W. J. Van Asch, J.P. Malet, and S. Gröndahl, “A GISbased numerical model for simulating the kinematics of mud and debris flows over complex terrain,” Natural Hazards and Earth System Science, vol. 9, no. 6, pp. 1897–1909, 2009. View at: Publisher Site  Google Scholar
 D. Laigle, P. Lachamp, and M. Naaim, “SPHbased numerical investigation of mudflow and other complex fluid flow interactions with structures,” Computational Geosciences, vol. 11, no. 4, pp. 297–306, 2007. View at: Publisher Site  Google Scholar
 D. Rickenmann, D. Laigle, B. W. McArdell, and J. Hübl, “Comparison of 2D debrisflow simulation models with field events,” Computational Geosciences, vol. 10, no. 2, pp. 241–264, 2006. View at: Publisher Site  Google Scholar
 C.L. Shieh, C.D. Jan, and Y.F. Tsai, “A numerical simulation of debris flow and its application,” Natural Hazards, vol. 13, no. 1, pp. 39–54, 1996. View at: Publisher Site  Google Scholar
 M. A. Ahmed, K. Ichiro, and S. Yasuyuki, “Simulation of threedimensional rapid freesurface granular flow past different types of obstructions using the SPH method,” Journal of Glaciology, vol. 62, no. 232, pp. 1–13, 2016. View at: Publisher Site  Google Scholar
 J. J. Monaghan and R. A. Gingold, “Shock simulation by the particle method SPH,” Journal of Computational Physics, vol. 52, no. 2, pp. 374–389, 1983. View at: Publisher Site  Google Scholar
 J. J. Monaghan, “On the integration of the SPH equations for a dusty fluid with high drag,” European Journal of Mechanics—B/Fluids, vol. 79, pp. 454–462, 2020. View at: Publisher Site  Google Scholar
 W. T. Zhang, C. Q. Shi, Y. An, S. H. Yang, and Q. Q. Liu, “Viscous elastoplastic SPH model for the longdistance highspeed landslide,” International Journal of Computational Methods, vol. 16, no. 2, Article ID 1846011, 2018. View at: Publisher Site  Google Scholar
 K. K. Pijush, M. C. Ira, and R. D. David, Fluid Mechanics, Academic Press, Cambridge, MA, USA, 2012.
 J. J. Monaghan, “Simulating free surface flows with SPH,” Journal of Computational Physics, vol. 110, no. 2, pp. 399–406, 1994. View at: Publisher Site  Google Scholar
 P. M. Campbell, “Some new algorithms for boundary value problems in smooth particle hydrodynamic,” Tech. Rep., Mission Research Corporation, Albuquerque, New Mexico, 1989, Technical Report DNATR88286. View at: Google Scholar
 J. C. Marongiu, F. Leboeuf, and E. Parkinson, “Numerical simulation of the flow in a Pelton turbine using the meshless method smoothed particle hydrodynamics: a new simple solid boundary treatment,” Proceedings of the Institution of Mechanical Engineers, Part A: Journal of Power and Energy, vol. 221, no. 6, pp. 849–856, 2007. View at: Publisher Site  Google Scholar
 C. X. Liang, C. Y. Zhang, H. Y. Liu, Z. R. Qin, and Q. Wang, “Modeling low Reynolds number incompressible flows with curved boundaries using SPH,” International Journal for Numerical Methods in Fluids, vol. 68, no. 9, pp. 1173–1188, 2012. View at: Publisher Site  Google Scholar
 K. Hu, F. Wei, and Y. Li, “Realtime measurement and preliminary analysis of debrisflow impact force at Jiangjia Ravine, China,” Earth Surface Processes and Landforms, vol. 36, no. 9, pp. 1268–1278, 2011. View at: Publisher Site  Google Scholar
 S. Zhang, L. Zhang, Q. Qi, Q. Li, and P. Shi, “Unsteady hydromagnetic forced convective heat transfer flow of a micropolar fluid along a porous wedge with convective surface boundary condition,” International Journal of Heat and Technology, vol. 33, no. 3, pp. 127–132, 2015. View at: Publisher Site  Google Scholar
 G. Wang, S. Tian, B. Hu, Z. Xu, J. Chen, and X. Kong, “Evolution pattern of tailings flow from dam failure and the buffering effect of debris blocking dams,” Water, vol. 11, no. 11, pp. 2388–2401, 2019. View at: Publisher Site  Google Scholar
 S. G. Philippe, N. Loan, T. Ronald, and S. Tomoya, “Smoothedparticle hydrodynamics numerical modeling of structures impacted by tsunami bores,” Journal of Waterway, Port, Coastal and Ocean Engineering, vol. 140, no. 1, pp. 68–81, 2014. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2020 QingYun Zeng et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.