Numerical three-dimensional modeling of earthen dam piping failure

A physically-based numerical three-dimensional earthen dam piping failure model is developed for homogeneous and zoned soil dams.This model is an erosion model,coupled with force/moment equilibrium analyses.Orifice flow and two-dimensional(2D)shallow water equations(SWE)are solved to simulate dam break flows at different breaching stages.Erosion rates of different soils with different construction compaction efforts are calculated using corresponding erosion formulae.The dam's real shape,soil properties,and surrounding area are programmed.Large outer 2D-SWE grids are used to control upstream and downstream hydraulic conditions and control the boundary conditions of orifice flow,and inner 2D-SWE flow is used to scour soil and perform force/moment equilibrium analyses.This model is validated using the European Commission IMPACT(Investigation of Extreme Flood Processes and Uncertainty)Test#5 in Norway,Teton Dam failure in Idaho,USA,and Quail Creek Dike failure in Utah,USA.All calculated peak outflows are within 10%errors of observed values.Simulation results show that,for a V-shaped dam like Teton Dam,a piping breach location at the abutment tends to result in a smaller peak breach outflow than the piping breach location at the dam's center;and if Teton Dam had broken from its center for internal erosion,a peak outflow of 117851 m'/s,which is 81%larger than the peak outflow of 65120 m3/s released from its right abutment,would have been released from Teton Dam.A lower piping inlet elevation tends to cause a faster/earlier piping breach than a higher piping inlet elevation.

Analysis and Modeling of Slope Stability in the Three-Gorges Dam Reservoir (China)——The Case of Huangtupo Landslide

The water level in the Three Gorges Dam reservoir is expected to change between the elevations of 145 m and 175 m, as a function of the flood control implementation and the intensity of the annual flood. As a matter of fact, the hydraulical and mechanical loadings, related to the water level modifications, will result in alterations in the slope stability conditions. The town of Badong (Hubei), of 20 000 inhabitants, is one of the towns which was submerged by the impoundment of the reservoir. As a consequence, the new town of Badong was constructed on a nearby site which appeared to be partly an unstable site. A part of this site corresponds to an old landslide, the Huangtupo landslide, the base of which had to be submerged by the water of the reservoir. The analysis of the Huangtupo landslide, taking into account various events scenarios, drainage and reinforcement measures and monitoring devices, allows to illustrate the general process implemented all along the reservoir in order to mitigate the landslide hazard.

Monitoring models for base flow effect and daily variation of dam seepage elements considering time lag effect

Affected by external environmental factors and evolution of dam performance, dam seepage behavior shows nonlinear time-varying characteristics. In this study, to predict and evaluate the long-term development trend and short-term fluctuation of the dam seepage behavior, two monitoring models were developed, one for the base flow effect and one for daily variation of dam seepage elements. In the first model, to avoid the influence of the time lag effect on the evaluation of seepage variation with the time effect component of seepage elements, the base values of the seepage element and the reservoir water level were extracted using the wavelet multi-resolution analysis method, and the time effect component was separated by the established base flow effect monitoring model. For the development of the daily variation monitoring model for dam seepage elements, all the previous factors, of which the measured time series prior to the dam seepage element monitoring time may have certain influence on the monitored results, were considered. Those factors that were positively correlated with the analyzed seepage element were initially considered to be the support vector machine(SVM) model input factors, and then the SVM kernel function-based sensitivity analysis was performed to optimize the input factor set and establish the optimized daily variation SVM model. The efficiency and rationality of the two models were verified by case studies of the water level of two piezometric tubes buried under the slope of a concrete gravity dam.Sensitivity analysis of the optimized SVM model shows that the influences of the daily variation of the upstream reservoir water level and rainfall on the daily variation of piezometric tube water level are processes subject to normal distribution.

Dam-break Modeling in Alpine Valleys

Dam-break analysis is of great importance in mountain environment,especially where reservoirs are located upstream of densely populated areas and hydraulic hazard should be assessed for land planning purposes.Accordingly,there is a need to identify suitable operative tools which may differ from the ones used in flat flood-prone areas.This paper shows the results provided by a 1D and a 2D model based on the Shallow Water Equations(SWE) for dam-break wave propagation in alpine regions.The 1D model takes advantage of a topographic toolkit that includes an algorithm for pre-processing the Digital Elevation Model(DEM) and of a novel criterion for the automatic cross-section space refinement.The 2D model is FLO-2D,a commercial software widely used for flood routing in mountain areas.In order to verify the predictive effectiveness of these numerical models,the test case of the Cancano dam-break has been recovered from the historical study of De Marchi(1945),which provides a unique laboratory data set concerning the consequences of the potential collapse of the former Cancano dam(Northern Italy).The measured discharge hydrograph at the dam also provides the data to test a simplified method recently proposed for the characterization of the hydrograph following a sudden dam-break.

Cloud-Verhulst hybrid prediction model for dam deformation under uncertain conditions

Uncertainties existing in the process of dam deformation negatively influence deformation prediction. However, existing deformation pre- diction models seldom consider uncertainties. In this study, a cloud-Verhulst hybrid prediction model was established by combing a cloud model with the Verhulst model. The expectation, one of the cloud characteristic parameters, was obtained using the Verhulst model, and the other two cloud characteristic parameters, entropy and hyper-entropy, were calculated by introducing inertia weight. The hybrid prediction model was used to predict the dam deformation in a hydroelectric project. Comparison of the prediction results of the hybrid prediction model with those of a traditional statistical model and the monitoring values shows that the proposed model has higher prediction accuracy than the traditional sta- tistical model. It provides a new approach to predicting dam deformation under uncertain conditions.

Modeling of Wind Erosion of the Aitik Tailings Dam Using SWEEP Model

Waste impoundments are usually not protected by any growing crops, soil roughness, or wind barriers, and dust due to wind erosion of the surface waste material can be a major environmental problem. The tailings dam in Aitik is highly susceptible to wind erosion when it is dry and windy. Strong dust storms may occur during dry seasons. These events are short-lived but may transport a huge amount of particulate matter and adversely impact air quality in downwind areas. This study estimated material loss from the Aitik tailings dam using SWEEP model. The simulated total material loss, saltation and creep loss, suspension loss and PM10 loss were 4.4941 kg/m2, 0.042 kg/m2, 4.4559 kg/m2, 0.264057 kg/m2 respectively. To control the dust generation, wind breaks could be suggested to be installed at 800 m downwind the west side edge of the tailings dam to reduce the saltation and creep and it would in turn reduce suspended emission. The measured total suspended material for five days by exposure profiling method was 0.475 kg/m2, which had a big difference with a simulate value of 4.4559 kg/m2. Many reasons contributed to the disagreement of simulated value and measured value in this study. In order to use the model for mining industry, more tests are needed to validate the modeling result, and calibration methods could be useful in adjusting the internal coefficients and empirical equations.

Centrifugal model tests and numerical simulations for barrier dam break due to overtopping

The present study focuses on the breaching process and failure of barrier dams due to overtopping. In this work, a series of centrifugal model tests is presented to examine the failure mechanisms of landslide dams. Based on the experimental results, failure process and mechanism of barrier dam due to overtopping are analyzed and further verified by simulating the experimental overtopping failure process. The results indicate that the barrier dam will develop during the entire process of overtopping in the width direction, whereas the breach will cease to develop at an early stage in the depth direction because of the large particles that accumulate on the downstream slope. Moreover, headcut erosion can be clearly observed in the first two stages of overtopping, and coarsening on the downstream slope occurs in the last stage of overtopping. Thus, the bottom part of the barrier dam can survive after dam breaching and full dam failure becomes relatively rare for a barrier dam. Furthermore, the remaining breach would be smaller than that of a homogeneous cohesive dam under the same conditions.

A Study on Physical Model Test for Cheongpyeong Dam Discharge Recalculation

Cheongpyeong Dam was built in 1943 for the purpose of power generation. Since its construction, discharge flow data based on a theoretical formula have been used to the present times and this leads to a problem of accurate discharge flow information not being available. In particular, Cheongpyeong Dam has been partially repaired and modified to maintain the dam structure over a long period of time and is not being properly reflected with changes at the downstream of the river caused by river improvement projects and sedimentation in reservoir. With a goal to improve Cheongpyeong Dam discharge flow calculation, this study aimed at verifying discharge capability and discharge flow by damper opening in relation to the previously suggested discharge flow through a hydraulic model test based on an accurate reproduction of the dam structure and surrounding topographies as in the present conditions. In this study, a hydraulic model test was conducted to examine the discharge flow of Cheongpyeong Dam. In addition, a comparative examination was carried out against the existing discharge flow proposed using theoretical equations. As a verification of the discharge flow of Cheongpyeong Dam, discharge flows in all sluices and a single sluice were examined. Then, the impact of sluice interference caused by the dam structure consisting with 24 sluices was investigated. As a result of the examination, it was found that the difference between discharge flow calculated using the existing theoretical equations and discharge flow derived from the hydraulic model test was insignificant. Based on the results of hydraulic model test, a formula to estimate stage-discharge flow at a sluice was derived and suggested.

Evaluation of behaviors of earth and rockfill dams during construction and initial impounding using instrumentation data and numerical modeling

In this study,the behavior of Gavoshan dam was evaluated during construction and the first impounding.A two-dimensional（2D） numerical analysis was conducted based on a finite difference method on the largest cross-section of the dam using the results of instrument measurements and back analysis.These evaluations will be completed in the case that back analysis is carried out in order to control the degree of the accuracy and the level of confidence of the measured behavior since each of the measurements could be controlled by comparing it to the result obtained from the numerical model.Following that,by comparing the results of the numerical analysis with the measured values,it is indicated that there is a proper consistency between these two values.Moreover,it was observed that the dam performance was suitable regarding the induced pore water pressure,the pore water pressure ratio r;,settlement,induced stresses,arching degree,and hydraulic fracturing probability during the construction and initial impounding periods.The results demonstrated that the maximum settlement of the core was 238 cm at the end of construction.In the following 6 years after construction（initial impounding and exploitation period）,the accumulative settlement of the dam was 270 cm.It is clear that 88% of the total settlement of the dam took place during dam construction.The reason is that the clay core was smashed in the wet side,i.e.the optimum moisture content.Whereas the average curving ratio was 0.64 during dam construction; at the end of the initial impounding,the maximum amount of curving ratio in the upstream was 0.81,and the minimum（critical） amount in the downstream was 0.52.It was also concluded that this dam is safe in comparison with the behaviors of other similar dams in the world.

A simplified physically-based breach model for a high concrete-faced rockfill dam:A case study

A simplified physically-based model was developed to simulate the breaching process of the Gouhou concrete-faced rockfill dam （CFRD）, which is the only breach case of a high CFRD in the world. Considering the dam height, a hydraulic method was chosen to simulate the initial scour position on the downstream slope, with the steepening of the downstream slope taken into account; a headcut erosion formula was adopted to simulate the backward erosion as well. The moment equilibrium method was utilized to calculate the ultimate length of a concrete slab under its self-weight and water loads. The calculated results of the Gouhou CFRD breach case show that the proposed model provides reasonable peak breach flow, final breach width, and failure time, with relative errors less than 15% as compared with the measured data. Sensitivity studies show that the outputs of the proposed model are more or less sensitive to different parameters. Three typical parametric models were compared with the proposed model, and the comparison demonstrates that the proposed physically-based breach model performs better and provides more detailed results than the parametric models.

SEFIDROOD RIVER SUB-WATERSHED-DAM-ESTUARY AND DEGRADATION MODEL:A HOLISTIC APPROACH IN IRAN

The major concern of this article is to address the shortcoming and outgoing effects of the human activities on the landscape patterns and their consequences in the Sefidrood River watershed in Iran. A flow of data includes three inputs; each of them belongs to one part of three zones of a fluvial system. The three parts of the Sefidrood River fluvial system include Zone 1,a sub-watershed as degradation modeling site,Zone 2,Sefidrood Dam as dam site,and Zone 3,17km away from the Sefidrood River path to the Caspian Sea as ending point site. The degradation model in the Zone 1 provides a suitable mean for decision support system to decrease the human impacts on each small district. The maximum number for degradation coefficient belongs to the small district with the highest physiographic density,relatively cumulative activities,and a lower figure for the habitat vulnerability. The human degradation impact were not limited to the upstream. The investigation to the Sefidrood Dam and ending point of the Sefidrood River depicts that sedimentation continues as a significant visual impact in the Sefidrood Dam reservoir and the estuary.

Application of thermodynamics-based rate-dependent constitutive models of concrete in the seismic analysis of concrete dams

This paper discusses the seismic analysis of concrete dams with consideration of material nonlinearity. Based on a consistent rate-dependent model and two thermodynamics-based models, two thermodynamics-based rate-dependent constitutive models were developed with consideration of the influence of the strain rate. They can describe the dynamic behavior of concrete and be applied to nonlinear seismic analysis of concrete dams taking into account the rate sensitivity of concrete. With the two models, a nonlinear analysis of the seismic response of the Koyna Gravity Dam and the Dagangshan Arch Dam was conducted. The results were compared with those of a linear elastic model and two rate-independent thermodynamics-based constitutive models, and the influences of constitutive models and strain rate on the seismic response of concrete dams were discussed. It can be concluded from the analysis that, during seismic response, the tensile stress is the control stress in the design and seismic safety evaluation of concrete dams. In different models, the plastic strain and plastic strain rate of concrete dams show a similar distribution. When the influence of the strain rate is considered, the maximum plastic strain and plastic strain rate decrease.

Fluid solid coupling model based on endochronic damage for roller compacted concrete dam

According to the characteristics of thin-layer rolling and pouring construction technology and the complicated mechanical behavior of the roller compacted concrete dam(RCCD)construction interface,a constitutive model of endochronic damage was established based on the endochronic theory and damage mechanics.The proposed model abandons the traditional concept of elastic-plastic yield surface and can better reflect the real behavior of rolled control concrete.Basic equations were proposed for the fluid solid coupling analysis,and the relationships among the corresponding key physical parameters were also put forward.One three-dimensional finite element method(FEM)program was obtained by studying the FEM type of the seepage-stress coupling intersection of the RCCD.The method was applied to an actual project,and the results show that the fluid solid interaction influences dam deformation and dam abutment stability,which is in accordance with practice.Therefore,this model provides a new method for revealing the mechanical behavior of RCCD under the coupling field.

Coupling model for assessing anti-seepage behavior of curtain under dam foundation

The seepage under a dam foundation is mainly controlled by the performance of the curtain.Its anti-seepage behavior may be weakened by the long term physic-chemical actions from groundwater.According to seepage hydraulics and geochemistry theory,a coupling model for assessing the behavior of the curtain under a dam foundation is set up,which consists of seepage module,solute transport module,geochemistry module and curtain erosion module,solved by FEM.A case study was carried out.The result shows that the curtain efficiency is weakened all the time.Aqueous calcium from the curtain is always in dissolution during the stress period for simulation,which leads to the increasing amount in groundwater reaching 846.35-865.312 g/m3.Within the domain,reaction extent differs in different parts of the curtain.The dissolution of Ca(OH)2 accounts to 877.884 g/m3 near the bottom and is much higher than that of the other parts.The erosion is much more serious near the bottom of the curtain than the other parts,which is the same to the upstream and downstream.Calcium dissolution is mainly controlled by hydraulic condition and dispersion,and it varies in a non-linear way within the domain.

Use River Pollutant Modeling to Simulate and Predict the Change in the Damietta Branch Water Quality before and after Construction of the Ethiopian Dam

This research was conducted on the Damietta branch of the Nile River, Egypt. The Damietta branch receives pollution loadings from the Omar-Bek drain and two power stations located along the path of the branch. The main objective of this research consisted of comparing between the Damietta branch water quality before and after the Ethiopian Dam is built. This comparison was conducted by using the river pollutant (RP) modeling. First, the actual data and the modeling results were compared in order to prove the efficiency and validity of the RP modeling. Findings from regression analysis yielded a strong positive linear relationship (r = 0.987) between the two results. The modeling results showed that Omar-Bek drain had less impact on the Damietta branch water quality. The results also showed that the effluent discharge from the two power stations affected water quality and aquatic life because large quantities of warm and polluted water discharged back into the Damietta branch. The results also showed that constructing the Ethiopian Renaissance Dam would slightly increase pollutants concentrations in the Damietta branch and that this increase would cause a slight deterioration in water quality.

Seepage simulation of high concrete-faced rockfill dams based on generalized equivalent continuum model

This research focused on the three-dimensional(3 D) seepage field simulation of a high concrete-faced rockfill dam(CFRD) under complex hydraulic conditions. A generalized equivalent continuum model of fractured rock mass was used for equivalent continuous seepage field analysis based on the improved node virtual flow method. Using a high CFRD as an example, the generalized equivalent continuum range was determined, and a finite element model was established based on the terrain and geological conditions, as well as structural face characteristics of the dam area. The equivalent seepage coefficients of different material zones or positions in the dam foundation were calculated with the Snow model or inverse analysis. Then, the 3 D seepage field in the dam area was calculated under the normal water storage conditions, and the corresponding water head distribution, seepage flow, seepage gradient, and seepage characteristics in the dam area were analyzed. The results show that the generalized equivalent continuum model can effectively simulate overall seepage patterns of the CFRD under complex hydraulic conditions and provide a reference for seepage analysis of similar CFRDs.

Numerical modeling of earthen dam breach due to piping failure

Based on model tests of earthen dam breach due to piping failure, a numerical model was developed.A key difference from previous research is the assumption that the cross-section of the pipe channel is an arch, with a rectangle at the bottom and a semicircle at the top before the collapse of the pipe roof, rather than a rectangular or circular cross-section.A shear stress-based erosion rate formula was utilized, and the arched pipe tunnel was assumed to enlarge along its length and width until the overlying soil could no longer maintain stability.Orifice flow and open channel flow were adopted to calculate the breach flow discharge for pressure and free surface flows, respectively.The collapse of the pipe roof was determined by comparing the weight of the overlying soil and the cohesion of the soil on the two sidewalls of the pipe.After the collapse, overtopping failure dominated, and the limit equilibrium method was adopted to estimate the stability of the breach slope when the water flow overtopped.In addition, incomplete and base erosion, as well as one-and two-sided breaches were taken into account.The USDAARS-HERU model test P1, with detailed measured data, was used as a case study, and two artificially filled earthen dam failure cases were studied to verify the model.Feedback analysis demonstrates that the proposed model can provide satisfactory results for modeling the breach flow discharge and breach development process.Sensitivity analysis shows that the soil erodibility and initial piping position significantly affect the prediction of the breach flow discharge.Furthermore, a comparison with a well-known numerical model shows that the proposed model performs better than the NWS BREACH model.

DAM大气模型构建与验证

根据热层物理、经验模型原理和代码分析,研究模型构建方法。进而剖析国内热层大气探测和热层模型构建现状,提出存在的困难和未来的发展思路。以GOST模型为基础,分析模型的工作机制、地磁扰动期大气密度预报误差来源和密切相关的模型系数,推导模型密度对相关模型系数的偏导数矩阵。利用天基实测密度,有针对性地构建磁暴期大气模型(DAM)。并通过独立于建模的实测密度数据,验证DAM模型性能。统计发现,地磁活动指数Ap介于100~132时GOST,MSIS00和DAM模型的相对误差均值依次为64.32%,-176.72%,-14.83%。Ap指数80~132时,相对误差均值对应为77.44%,-136.74%,-14.14%,DAM模型性能较GOST和MSIS00均有明显提升。证明通过搭建大气模型框架和实测密度数据估计模型参数的建模方法是可行和有效的。

Comparison of depth-averaged concentration and bed load flux sediment transport models of dam-break flow

This paper presents numerical simulations of dam-break flow over a movable bed. Two different mathematical models were compared: a fully coupled formulation of shallow water equations with erosion and deposition terms(a depth-averaged concentration flux model), and shallow water equations with a fully coupled Exner equation(a bed load flux model). Both models were discretized using the cell-centered finite volume method, and a second-order Godunov-type scheme was used to solve the equations. The numerical flux was calculated using a Harten, Lax, and van Leer approximate Riemann solver with the contact wave restored(HLLC). A novel slope source term treatment that considers the density change was introduced to the depth-averaged concentration flux model to obtain higher-order accuracy. A source term that accounts for the sediment flux was added to the bed load flux model to reflect the influence of sediment movement on the momentum of the water. In a onedimensional test case, a sensitivity study on different model parameters was carried out. For the depth-averaged concentration flux model,Manning's coefficient and sediment porosity values showed an almost linear relationship with the bottom change, and for the bed load flux model, the sediment porosity was identified as the most sensitive parameter. The capabilities and limitations of both model concepts are demonstrated in a benchmark experimental test case dealing with dam-break flow over variable bed topography.