Simulation of two-phase debris flow scouring bridge pier

Debris flows are typical two-phase flows, which commonly accompany erosion in mountainous areas, and may destroy bridge engineering by scouring. In this study, a physically-based two-phase model is applied for the simulation of debris flow scouring of bridge pier. In this model, the shear stress of debris flow on an erodible bed is considered to be a function of the solid shear stress, fluid shear stress, and volume fraction; accordingly, the erosion is incorporated into the two-phase model. Using a highaccuracy computational scheme based on the finite volume method, the model is employed for simulating a dynamic debris flow over an erodible bed. The numerical results are consistent with the experimental data, and verify the feasibility of the two-phase model. Moreover, a simple numerical test is performed to exhibit the fundamental behaviour of debris flow scouring of bridge pier, which shows that the degree of erosion on each side of the pier is higher compared to other areas. The scouring depth is influenced by the variations of solid volume fraction and velocity of debris flow and pier width.     

Velocity profile of debris flow based on quadratic rheology model

The quadratic rheology model considers the yield stress, viscous stress, turbulent stress and disperse stress, so it is used in this study to derive the velocity profile of debris flows. The quadratic model with the parabolic eddy viscosity was numerically solved, and an analytical solution was derived for the quadratic model with a constant eddy viscosity. These two solutions were compared with the Arai-Takahashi model that excluded the viscous stress and the yield stress. The three models were tested by using 17 experiment cases of debris flows over rigid beds. The results prove that the quadratic model with parabolic and constant eddy viscosities is applicable to muddy and granular flows, whereas the Arai-Takahashi model tends to overestimate the flow velocity near the water surface if a plug-like layer exists. In addition, the von Karman constant and the zero-velocity elevation in the three models are related to sediment concentration. The von Karman constant decreases first and then increases as the sediment concentration increases. The zero-velocity elevation is below the bed surface, likely due to the invalidity of the non-slip boundary condition for the debris flows over fixed beds.     

Two-dimensional numerical model for debris flows in the Jiangjia Gully, Yunnan Province

Debris flows are recurrent natural hazards in many mountainous regions.This paper presents a numerical study on the propagation of debris flows in natural erodible open channels,in which the bed erosion and sedimentation processes are important.Based on the Bingham fluid theory,a mathematical model of the two-dimensional non-constant debris flow is developed.The governing equations include the continuity and momentum conservation equations of debris flow,the sediment convection-diffusion equation,the bed erosion-deposition equation and the bed-sediment size gradation adjustment equation.The yield stress and shear stress components are included to describe the dynamic rheological properties.The upwind control-volume Finite Volume Method (FVM) is applied to discretize the convection terms.The improved SIMPLE algorithm with velocity-free-surface coupled correction is developed to solve the equations on non-orthogonal,quadrilateral grids.The model is applied to simulate a debris flow event in Jiangjia Gully,Yunnan Province and to predict the flow pattern and bed erosion-deposition processes.The results show the effectiveness of the proposed numercial model in debris flow simulation and potential hazard analysis.     

Deposition Morphology of Non-homogeneous Debris Flow and Its Energy Characteristics

Non-homogeneous two-phase debris flows are widely found in the western mountainous regions of China. To investigate the characteristics of the debris flow deposition process related to the morphology and extent of the debris fan, a series of physical experiments were carried out using an experimental flume. Some useful relationships were obtained to link the flow velocity with the geometric characteristics of deposition morphology and the corresponding area or volume. Based on these, some expressions about energy dissipation process in both the transport-deposition zone and deposition zone are presented, and improved equations describing solidliquid two-phase energy transformations in the specific deposition zone are also established. These results provide a basis for further investigating the underlying mechanisms of non-homogeneous debris flows, based upon which effective disaster control measures can be undertaken.     

Impacts of future climate change (2030-2059) on debris flow hazard: A case study in the Upper Minjiang River basin,China

An increase in extreme precipitation events due to future climate change will have a decisive influence on the formation of debris flows in earthquake-stricken areas. This paper aimed to describe the possible impacts of future climate change on debris flow hazards in the Upper Minjiang River basin in Northwest Sichuan of China, which was severely affected by the 2008 Wenchuan earthquake. The study area was divided into 1285 catchments, which were used as the basic assessment units for debris flow hazards. Based on the current understanding of the causes of debris flows, a binary logistic regression model was used to screen key factors based on local geologic, geomorphologic, soil, vegetation, and meteorological and climatic conditions. We used the weighted summation method to obtain a composite index for debris flow hazards, based on two weight allocation methods: Relative Degree Analysis and rough set theory. Our results showed that the assessment model using the rough set theory resulted in better accuracy. According to the bias corrected and downscaled daily climate model data, future annual precipitation (2030-2059) in the study area are expected to decrease, with an increasing number of heavy rainfall events. Under future climate change, areas with a high-level of debris flow hazard will be even more dangerous, and 5.9% more of the study area was categorized as having a high-level hazard. Future climate change will cause an increase in debris flow hazard levels for 128 catchments, accounting for 10.5% of the total area. In the coming few decades, attention should be paid not only to traditional areas with high-level of debris flow hazards, but also to those areas with an increased hazard level to improve their resilience to debris flow disasters.     

Changes in runout distances of debris flows over time in the Wenchuan earthquake zone

A large number of debris flows occurred in the Wenchuan earthquake zone after the 12 May 2008 earthquake.The risks posed by these debris flows were rather high.An appropriate model is required to predict the possible runout distance and impacted area.This paper describes a study on the runout characteristics of the debris flows that occurred in the Wenchuan earthquake zone over the past four years.A total of 120 debris flows are analyzed.Separate multivariate regression models are established for the runout distances of hill-slope debris flows and channelized debris flows.The control variables include type of debris flow,debris flow volume,and elevation difference.Comparison of the debris flows occurring before and after the earthquake shows that the runout distance increased after the earthquake due to sufficient material supply and increased mobility of the source materials.In addition,the runout distances of annual debris flow events in 2008,2010 and 2011 are analyzed and compared.There is a tendency that the runout distance decreases over time due to the decreasing source material volume and possible changes of debris flow type.Comparison between the debris flows in the earthquake zone and the debris flows in Swiss Alps,Canada,Austria,and Japan shows that the former have a smaller mobility.     

Volume estimation of small scale debris flows based on observations of topographic changes using airborne LiDAR DEMs

This paper describes a geographic information system(GIS)-based method for observing changes in topography caused by the initiation, transport, and deposition of debris flows using highresolution light detection and ranging(LiDAR) digital elevation models(DEMs) obtained before and after the debris flow events. The paper also describes a method for estimating the volume of debris flows using the differences between the LiDAR DEMs. The relative and absolute positioning accuracies of the LiDAR DEMs were evaluated using a real-time precise global navigation satellite system(GNSS) positioning method. In addition, longitudinal and cross-sectional profiles of the study area were constructed to determine the topographic changes caused by the debris flows. The volume of the debris flows was estimated based on the difference between the LiDAR DEMs. The accuracies of the relative and absolute positioning of the two LiDAR DEMs were determined to be ±10 cm and ±11 cm RMSE, respectively, which demonstrates the efficiency of the method for determining topographic changes at an scale equivalent to that of field investigations. Based on the topographic changes, the volume of the debris flows in the study area was estimated to be 3747 m3, which is comparable with the volume estimated based on the data from field investigations.     

Slope debris flows in the Wenchuan Earthquake area

Avalanches and landslides, induced by the Wenchuan Earthquake on May 12, 2008, resulted in a lot of disaggregated, solid material on slopes that could be readily mobilized as source material for debris flows. Rainstorms triggered numerous slope debris flows with great damage to highways and rivers over the subsequent two years. Slope debris flows (as opposed to channelized debris flows) are defined as phenomena in which high-concentration mixtures of debris and water flow down slopes for short distances to highways and river banks. Based on field investigations and measurements of 19 slope debris flows, their main characteristics and potential mitigation strategies were studied. High rainfall intensity is the main triggering factor. Critical rainfall intensities for simultaneous occurrence of single, several and numerous slope debris flow events were 20 mm/day, 30mm/day, and 90 mm/day, respectively. Field investigations also revealed that slope debris flows consist of high concentrations of cobbles, boulders and gravel. They are two-phase debris flows. The liquid phase plays the role of lubrication instead of transporting medium. Solid particles collide with each other and consume a lot of energy. The velocities of slope debris flows are very low, and their transport distances are only several tens of meters. Slope debris flows may be controlled by construction of drainage systems and by reforestation.     

Numerical modeling of bed deposition in rapid flow-like landslides: a case study of the Gaolou landslide in Shaanxi Province,China

Some flow-like landslides tend to lose materials while moving on a relatively dry sliding surface.This phenomenon is called bed deposition.In contrast to the bed entrainment phenomenon,bed deposition is relatively poorly understood.Therefore,an improved depth-averaged model is proposed to quantify this phenomenon.The deposition depth is calculated according to the momentum conservation of the deposited mass,and the rheological property of the sliding mass on the bottom is modified considering an abrupt increase in the depth of the sliding surface after deposition.Utilizing the proposed model,the Gaolou landslide,a typical flowlike landside occurring on October 6 th,2006 with an obvious bed deposition phenomenon in Shaanxi Province of China,is simulated to investigate the influence of bed deposition on its propagation process.The results indicate that the proposed model can effectively depict the bed deposition phenomenon in the Gaolou landslide.Bed deposition dissipates part of the kinetic energy of this landslide;thus,the simulated debris inundation area would likely be overestimated when this effect is neglected.On the other hand,the thin liquefied layer formed by the bed deposition process reduces the friction energy dissipation,contributing to the high mobility of this landslide.     

Effects of a flexible net barrier on the dynamic behaviours and interception of debris flows in mountainous areas

Flexible net barriers are a new type of effective mitigation measure against debris flows in valleys and can affect the kinematic energy and mass of debris flows. Here, ten flume tests were performed to study the dynamic behaviours of debris flows with differences in volumes, concentrations (solid volume fraction), and travel distances after interception by a uniform flexible net barrier. A high-speed camera was used to monitor the whole test process, and their dynamic behaviours were recorded. A preliminary computational framework on energy conversion is proposed according to the deposition mechanisms and outflow of debris flow under the effects of the flexible net barrier. The experimental results show that the dynamic interaction process between a debris flow and the flexible net barrier can be divided into two stages: (a) the two-phase impact of the leading edge of the debris flow with the net and (b) collision and friction between the body of the debris flow and intercepted debris material. The approach velocity of a debris flow decreases sharply (a maximum of 63%) after the interception by the net barrier, and the mass ratio of the debris material being intercepted and the kinetic energy ratio of the debris material being absorbed by the net barrier are close due to the limited interception efficiency of the flexible net barrier, which is believed to be related to the flexibility. The energy ratio of outflow is relative small despite the large permeability of the flexible net barrier.     

Strongly nonlinear long interfacial waves in uniform basic flows and their dispersion relationship

In this paper, long interfacial waves of finite amplitude in uniform basic flows are considered with the assumption that the aspect ratio between wavelength and water depth is small. A new model is derived using the velocities at arbitrary distances from the still water level as the velocity variables instead of the commonly used depth-averaged velocities. This significantly improves the dispersion properties and makes them applicable to a wider range of water depths. Since its derivation requires no assumption on wave amplitude, the model thus can be used to describe waves with arbitrary amplitude.     

GIS tools for preliminary debris-flow assessment at regional scale

The assessment of the areas endangered by debris flows is a major issue in the context of mountain watershed management. Depending on the scale of analysis, different methods are required for the assessment of the areas exposed to debris flows. While 2-D numerical models are advised for detailed mapping of inundation areas on individual alluvial fans, preliminary recognition of hazard areas at the regional scale can be adequately performed by less data-demanding methods, which enable priority ranking of channels and alluvial fans at risk by debris flows. This contribution focuses on a simple and fast procedure that has been implemented for regional-scale identification of debris-flow prone channels and prioritization of the related alluvial fans. The methodology is based on the analysis of morphometric parameters derived from Digital Elevation Models (DEMs). Potential initiation sites of debris flows are identified as the DEM cells that exceed a threshold of slope-dependent contributing area. Channel reaches corresponding to debris flows propagation, deceleration and stopping conditions are derived from thresholds of local slope. An analysis of longitudinal profiles is used for the computation of the runout distance of debris flows. Information on erosion-resistant bedrock channels and sediment availability surveyed in the field are taken into account in the applications. A set of software tools was developed and made available (https://github.com/HydrogeomorphologyTools) to facilitate the application of the procedure. This approach, which has been extensively validated by means of field checks, has been extensively applied in the eastern Italian Alps. This contribution discusses potential and limitations of the method in the frame of the management of small mountain watersheds.     

Artificial Neural Network-based prediction of glacial debris flows in the ParlungZangbo Basin,southeastern Tibetan Plateau,China

Accurate prediction on geological hazards can prevent disaster events in advance and greatly reduce property losses and life casualties.Glacial debris flows are the most serious hazards in southeastern Tibet in China due to their complexity in formation mechanism and the difficulty in prediction.Data collected from 102 glacier debris flow events from 31 gullies since 1970 and regional meteorological data from 1970 to 2019 in ParlungZangbo River Basin in southeastern Tibet were used for Artificial Neural Network(ANN)-based prediction of glacial debris flows.The formation mechanism of glacial debris flows in the ParlungZangbo Basin was systematically analyzed,and the calculations involving the meteorological data and disaster events were conducted by using the statistical methods and two layers fully connected neural networks.The occurrence probabilities and scales of glacial debris flows(small,medium,and large)were predicted,and promising results have been achieved.Through the proposed model calculations,a prediction accuracy of 78.33%was achieved for the scale of glacial debris flows in the study area.The prediction accuracy for both large-and medium-scale debris flows are higher than that for small-scale debris flows.The debris flow scale and the probability of occurrence increase with increasing rainfall and temperature.In addition,the K-fold cross-validation method was used to verify the reliability of the model.The average accuracy of the model calculated under this method is about 93.3%,which validates the proposed model.Practices have proved that the combination of ANN and disaster events can provide sound prediction on geological hazards under complex conditions.     

Debris flows introduced in landslide deposits under rainfall conditions: The case of Wenjiagou gully

Debris flows often occur in landslide deposits during heavy rainstorms. Debris flows are initiated by surface water runoff and unsaturated seepage under rainfall conditions. A physical model based on an infinitely long, uniform and void-rich sediment layer was applied to analyze the triggering of debris-flow introduced in landslide deposits. To determine the initiation condition for rainfall-induced debris flows, we conducted a surface water runoff and saturated-unsaturated seepage numerical program to model rainfall infiltration and runoff on a slope. This program was combined with physical modeling and stability analysis to make certain the initiation condition for rainfall-introduced debris flows. Taking the landslide deposits at Wenjiagou gully as an example, the initiation conditions for debris flow were computed. The results show that increase height of surface-water runoff and the decrease of saturated sediment shear strength of are the main reasons for triggering debris-flows under heavy rainfall conditions. The debris-flow triggering is affected by the depth of surface-water runoff, the slope saturation and shear strength of the sediment.     

Impact failure models and application condition of trees in debris-flow hazard mitigation

Forestry has played an important role in hazard mitigation associated with debris flows. Most forest mitigation measures refer to the experience of soil and water conservation, which disregard the destructive effect of debris flows, causing potentially serious consequences. Determination of the effect of a forest on reducing debris-flow velocity and even stopping debris flows requires distinguishing between when the debris flow will destroy the forest and when the trees will withstand the debris-flow impact force. In this paper, we summarized two impact failure models of a single tree: stem breakage and overturning. The influences of different tree sizes characteristics(stem base diameter, tree weight, and root failure radius) and debris-flow characteristics(density, velocity, flow depth, and boulder diameter) on tree failure were analyzed. The observations obtained from the model adopted in this study show that trees are more prone to stem breakage than overturning. With an increase in tree size, the ability to resist stem breakage and overturning increases. Debris-flow density influences the critical failure conditions of trees substantially less than the debrisflow velocity, depth, and boulder diameter. The application conditions of forests in debris-flow hazard mitigation were proposed based on the analysis of the model results. The proposed models were applied in the Xiajijiehaizi Gully as a case study, and the results explain the destruction of trees in the forest dispersing zone. This work provides references for implementing forest measures for debris-flow hazard mitigation.     

Debris flow susceptibility assessment using a probabilistic approach: A case study in the Longchi area,Sichuan province,China

The Longchi area with the city of Dujiangyan, in the Sichuan province of China, is composed of Permian stone and diorites and Triassic sandstones and mudstones intercalated with slates. An abundance of loose co-seismic materials were present on the slopes after the May 12, 2008 Wenchuan earthquake, which in later years served as source material for rainfall-induced debris flows or shallow landslides. A total of 48 debris flows, all triggered by heavy rainfall on 13th August 20l0, are described in this paper. Field investigation, supported by remote sensing image interpretation, was conducted to interpret the co-seismic landslides in the debris flow gullies. Specific characteristics of the study area such as slope, aspect, elevation, channel gradient, lithology, and gully density were selected for the evaluation of debris flow susceptibility. A score was given to all the debris flow gullies based on the probability of debris flow occurrence for the selected factors. In order to get the contribution of the different factors, principal component analyses were applied. A comprehensive score was obtained for the 48 debris flow gullies which enabled us to make a susceptibility map for debris flows with three classes. Twenty-two gullies have a high susceptibility, twenty gullies show a moderate susceptibility and six gullies have a low susceptibility for debris flows.     

Curvature derived from LiDAR digital elevation models as simple indicators of debris-flow susceptibility

To mitigate the damage caused by debris flows resulting from heavy precipitation and to aid in evacuation plan preparation, areas at risk should be mapped on a scale appropriate for affected individuals and communities. We tested the effectiveness of simply identifying debris-flow hazards through automated derivation of surface curvatures using LiDAR digital elevation models. We achieved useful correspondence between plan curvatures and areas of existing debris-flow damage in two localities in Japan using the analysis of digital elevation models(DEMs). We found that plan curvatures derived from 10 m DEMs may be useful to indicate areas that are susceptible to debris flow in mountainous areas. In residential areas located on gentle sloping debris flow fans, the greatest damage to houses was found to be located in the elongated depressions that are connected to mountain stream valleys. Plan curvaturederived from 5 m DEM was the most sensitive indicators for susceptibility to debris flows.     

应用Level Set方法模拟三维气泡与液滴的运动

将Level Set方法与ICE数值格式相结合,建立了求解含有运动界面的不可压缩气液两相流动问题的三维模型及数值解法。应用所建立的模型和算法三维模拟了4种典型气泡和液滴的运动过程,包括气泡从单个喷嘴生成和上升的过程,两个大小不同的空气气泡在水中自由上升时的相互作用过程,水滴下落冲击水面的过程和水滴与颗粒正向等温碰撞过程。模拟结果表明,所建立的模型与算法能准确方便地捕捉到三维气泡和液滴界面的破碎与融合等拓扑结构的变化,说明Level Set方法是一种行之有效且大有前途的处理气液运动界面问题的方法。     

Early warning model for slope debris flow initiation

Early warning model of debris flow is important for providing local residents with reliable and accurate warning information to escape from debris flow hazards. This research studied the debris flow initiation in the Yindongzi gully in Dujiangyan City, Sichuan province, China with scaled-down model experiments. We set rainfall intensity and slope angle as dominating parameters and carried out 20 scaled-down model tests under artificial rainfall conditions. The experiments set four slope angles(32°, 34°, 37°, 42°) and five rainfall intensities(60 mm/h, 90 mm/h, 120 mm/h, 150 mm/h, and 180 mm/h) treatments. The characteristic variables in the experiments, such as, rainfall duration, pore water pressure, moisture content, surface inclination, and volume were monitored. The experimental results revealed the failure mode of loose slope material and the process of slope debris flow initiation, as well as the relationship between the surface deformation and the physical parameters of experimental model. A traditional rainfall intensity-duration early warning model(I-D model) was firstly established by using a mathematical regression analysis, and it was then improved into ISD model and ISM model(Here, I is rainfall Intensity, S is Slope angle, D is rainfall Duration, and M is Moisture content). The warning model can provide reliable early warning of slope debris flow initiation.     

On the criteria to create a susceptibility map to debris flow at a regional scale using Flow-R

Studies on susceptibility to debris flows at regional scale(100-1000 km~2) are important for the protection and management of mountain areas. To reach this objective, routing models, mainly based on land topography, can be used to predict susceptible areas rapidly while necessitating few input data. In this research, Flow-R model is implemented to create the susceptibility map for the debris flow of the Vizze Valley(BZ, North-Eastern Italy; 134 km~2). The analysis considers the model application at local scale for three sub-catchments and then it explores the model upscaling at the regional scale by verifying two methods to generate the source areas of debris-flow initiation. Using data of an extreme event occurred in the Vizze Valley(4 August 2012) and historical information, the modeling verification highlights that the propagation parameters are relatively simple to set in order to obtain correct runout distances. A double DTM filtering-using a threshold for the upslope contributing area(0.1 km~2) and a threshold for the terrain-slope angle(15°)-provides a satisfactory prediction of source areas and susceptibility map within the geological conditions of the Vizze Valley.