Mitigation of the amplitude and celerity of dam break shock wave using combination of groyne and vertical pier

This paper aims to numerically investigate the impact of the combination of groynes and vertical piers on mitigating the amplitude and celerity of dam break wave-front.First,a robust numerical model based on the Finite Volume Method(FVM)with the Harten,Lax and van Leer(HLL)approximate Reimann’s solver’s technique is developed and validated with other research studies.Then many different scenarios including group of piers,group of groynes with different shapes and combination of groynes and piers are considered for dam break flow modeling.The final results showed that group of piers alone or group of groynes alone have no significant effects on the wave-front depth and celerity but,combination of L-head groyne with piers can decrease the amplitude and celerity of the dam break wavefront or flood wave.In this case,a blockage zone between the groynes and intensive backwater effects behind them were observed.

Modelling Dam Break Evolution over a Wet Bed with Smoothed Particle Hydrodynamics: A Parameter Study

When investigating water flow in spillways and energy dissipation, it is important to know the behavior of the free surfaces. To capture the real dynamic behavior of the free surfaces is therefore crucial when performing simulations. Today, there is a lack in the possibility to model such phenomenon with traditional methods. Hence, this work focuses on a parameter study for one alternative simulation tool available, namely the meshfree, Lagrangian particle method Smoothed Particle Hydrodynamics (SPH). The parameter study includes the choice of equation-of-state (EOS), the artificial viscosity constants, using a dynamic versus a static smoothing length, SPH particle spatial resolution and the finite element method (FEM) mesh scaling of the boundaries. The two dimensional SPHERIC Benchmark test case of dam break evolution over a wet bed was used for comparison and validation. The numerical results generally showed a tendency of the wave front to be ahead of the experimental results, i.e. to have a greater wave front velocity. The choice of EOS, FEM mesh scaling as well as using a dynamic or a static smoothing length showed little or no significant effect on the outcome, though the SPH particle resolution and the choice of artificial viscosity constants had a major impact. A high particle resolution increased the number of flow features resolved for both choices of artificial viscosity constants, but at the expense of increasing the mean error. Furthermore, setting the artificial viscosity constants equal to unity for the coarser cases resulted in a highly viscous and unphysical solution, and thus the relation between the artificial viscosity constants and the particle resolution and its impact on the behavior of the fluid needed to be further investigated.

A Dam Break Analysis Using HEC-RAS

This paper firstly described the dam break in the aspects of theories and models. Break parameters prediction, the understanding of dam break mechanics, peak outflow prediction were shown as the essential for the dam break analysis, and eventually determined the loss of the damages. Secondly, as an application example, Foster Joseph Sayers Dam break was further modeled and analyzed using USACE Hydrologic Engineering Center’s River Analysis System (HEC-RAS) model based on available geometry data. The results show that dam break is a complicated and comprehensive process involving lots of principles. Combination of mechanics and case studies, reflection of predominant mechanisms of headcut erosion, more specific categorization of dam, prudent investigation and inference of dam break process are needed in developing a satisfactory dam break simulation model. Foster Joseph Sayers Dam break due to piping elongates the time period of high water surface level, which increases the duration of risk. However, the dam break does not increase the downstream maximum water surface elevation (Max. W.S. Elev) significantly at previous design Probable Maximum Flood (PMF). Dam break has a greater impact on the downstream location where is closer to the dam in accordance with the comparison of the hydrographs at different downstream locations. Sensitivity analysis demonstrates that the changes of dam break parameters had no much influence on the downstream Max. W.S. Elev.

Numerical Study of Dam-Break Induced Tsunami-Like Bore with A Hump of Different Slopes

Numerical simulation of dam-break wave, as an imitation of tsunami hydraulic bore, with a hump of different slopes is performed in this paper using an in-house code, named a Constrained Interpolation Profile （CIP）-based model. The model is built on a Cartesian grid system with the Navier Stokes equations using a CIP method for the flow solver, and employs an immersed boundary method （IBM） for the treatment of solid body boundary. A more accurate interface capturing scheme, the Tangent of hyperbola for interface capturing/Slope weighting （THINC/SW） scheme, is adopted as the interface capturing method. Then, the CIP-based model is applied to simulate the dam break flow problem in a bumpy channel. Considerable attention is paid to the spilling type reflected bore, the following spilling type wave breaking, free surface profiles and water level variations over time. Computations are compared with available experimental data and other numerical results quantitatively and qualitatively. Further investigation is conducted to analyze the influence of variable slopes on the flow features of the tsunami-like bore.

Experimental Investigation on the Interactions Between Dam-Break Flow and a Floating Box

An experiment that investigates the interactions between dam-break flow and a floating box is reported in this paper.The interaction processes are described in detail,and the effects of reservoir depth,clearance between the structure and vertical wall,and structure motion on the observed interactions are discussed.Results indicate that interactions are mild if the reservoir depth is relatively small and most of the water flows through under the structure.However,as the reservoir depth increases,the interactions intensify,and overtopping and wave splashing occur.The incident wave is divided into two parts by the structure:1)a jet generated by wave impact and 2)the water flowing through under the structure.The water flowing through under the structure is the main factor influencing the motion of the structure.The phenomenon of secondary slamming is observed in all of the test cases,and the relative motion of the structure and free surface explains the occurrence of secondary slamming.

Draining Tangjiashan Barrier Lake after Wenchuan Earthquake and the flood propagation after the dam break

Tangjiashan Barrier Lake is one of the largest barrier lakes caused by the Wenchuan Earthquake. Its risk analysis, emergency plan and effect of the emergency plan are introduced in this paper. The dam height of Tangjiashan Barrier Dam is about 105 m, and the reservoir storage capacity is 3.2×108 m3. When the dam broke the flood peak were estimated to be larger than 48000 m3/s, which might cause a enormous disaster to the downstream cities and residents. A discharge channel with 13 m deep and 8 m wide was drug, so that the water may flow out of the lake before the dam breaks. As a result, the drainage and risk mitigation project are successful. During the drainage process, the flood peak was about 6500 m3/s, and about 1.6×108 m3 of water was drained off and the residual reservoir capacity was only 8.97×107 m3. A new channel with average width 100 m was formed, which can bear floods of 200 years frequency. The successful experience and the collected data can be used to deal with the similar natural disasters in future.

ONE-AND TWO-DIMENSIONAL COUPLED HYDRODYNAMICS MODEL FOR DAM BREAK FLOW

1-D and 2-D mathematical models for dam break flow were established and verified with the measured data in laboratory. The 1-D and 2-D models were then coupled, and used to simulate the dam break flow from the reservoir tail to the dam site, the propagation of dam break waves in the downstream channel, and the submergence of dam break flow in the downstream town with the hydrodynamics method. As a numerical example, the presented model was employed to simulate dam break flow of a hydropower station under construction. In simulation, different dam-break durations, upstream flows and water levels in front of dam were considered, and these influencing factors of dam break flow were analyzed, which could be referenced in planning and designing hydropower stations.

Analysis on Dam-Breaking Mode of Tangjiashan Barrier Dam in Beichuan County

Tangjiashan landslide is a typical high-speed consequent landslide of medium-steep dip angle. This landslide triggered by earthquake took place in about semi-minute. The relative sliding displacement is 900 meters, so average sliding speed is about 30 meters per second. The longitudinal length of barrier dam which is formed by high-speed landslide along river is 803.4 meters; and maximum width crossing river is 611.8 meters. And its volume is estimated about 20.37 million steres. Through detailed geological investigation of the barrier dam, together with early geological information before earthquake, geological structures of the barrier dam and its stability of upstream and downstream slopes are studied when water level reaches different elevations in condition of continual after shocks with seismic intensity of 7 or 8 Richter scale. On this basis, dam-breaking mode of barrier dam is discussed deeply. Thereby, analytic results provide significant guidance and advices to front headquarters of Tangjiashan barrier dam, so that some proper engineering measures can be implemented and flood discharge can be carried out well.

Dam-break Flood Simulation under Various Likely Scenarios and Mapping Using GIS:Case of a Proposed Dam on River Yamuna,India

The precision modeling of dam break floods can lead to formulation of proper emergency action plan to minimize flood impacts within the economic lifetime of the assets.Application of GIS techniques in integration with hydrological modeling for mapping of the flood inundated areas can play a momentous role in further minimizing the risk and likely damages.In the present study,dam break analysis using DAMBRK model was performed under various likely scenarios.Probable Maximum Flood (PMF)calculated for a return period of 1000 years using deterministic approach was adopted for dam break analysis of the proposed dam under various combinations of breach dimensions.The available downstream river cross-sections data sets were used as input in the model to generate the downstream flood profile.Dam break flow depths generated by the DAMBRK model under various combinations of structural failure are subsequently plotted on Digital Elevation Model(DEM)of the downstream of dam site to map the likely affected area.The simulation results reveals that in one particular case the flood without dam may be more intense if a rainfall of significant intensity takes place.

Coupled 2D Hydrodynamic and Sediment Transport Modeling of Megaflood due to Glacier Dam-break in Altai Mountains,Southern Siberia

One of the largest known megafloods on earth resulted from a glacier dam-break,which occurred during the Late Quaternary in the Altai Mountains in Southern Siberia.Computational modeling is one of the viable approaches to enhancing the understanding of the flood events.The computational domain of this flood is over 9460 km2 and about 3.784 × 106 cells are involved as a 50 m × 50 m mesh is used,which necessitates a computationally efficient model.Here the Open MP(Open Multiprocessing) technique is adopted to parallelize the code of a coupled 2D hydrodynamic and sediment transport model.It is shown that the computational efficiency is enhanced by over 80% due to the parallelization.The floods over both fixed and mobile beds are well reproduced with specified discharge hydrographs at the dam site.Qualitatively,backwater effects during the flood are resolved at the bifurcation between the Chuja and Katun rivers.Quantitatively,the computed maximum stage and thalweg are physically consistent with the field data of the bars and deposits.The effects of sediment transport and morphological evolution on the flood are considerable.Sensitivity analyses indicate that the impact of the peak discharge is significant,whilst those of the Manningroughness,medium sediment size and shape of the inlet discharge hydrograph are marginal.

Analysis of seismic disaster failure mechanism and dam-break simulation of high arch dam

Based on a Chinese national high arch dam located in a meizoseismal region, a nonlinear numerical analysis model of the damage and failure process of a dam-foundation system is established by employing a 3-D deformable distinct element code(3DEC) and its re-development functions. The proposed analysis model considers the dam-foundation-reservoir coupling effect, infl uence of nonlinear contact in the opening and closing of the dam seam surface and abutment rock joints during strong earthquakes, and radiation damping of far fi eld energy dissipation according to the actual workability state of an arch dam. A safety assessment method and safety evaluation criteria is developed to better understand the arch dam system disaster process from local damage to ultimate failure. The dynamic characteristics, disaster mechanism, limit bearing capacity and the entire failure process of a high arch dam under a strong earthquake are then analyzed. Further, the seismic safety of the arch dam is evaluated according to the proposed evaluation criteria and safety assessment method. As a result, some useful conclusions are obtained for some aspects of the disaster mechanism and failure process of an arch dam. The analysis method and conclusions may be useful in engineering practice.

Optimal Evacuation Scheme Based on Dam-Break Flood Numerical Simulation

The optimal evacuation scheme is studied based on the dam-break flood numerical simulation. A threedimensional dam-break mathematical model combined with the volume of fluid (VOF) method is adopted. According to the hydraulic information obtained from numerical simulation and selecting principles of evacuation emergency scheme, evacuation route analysis model is proposed, which consists of the road right model and random degree model. The road right model is used to calculate the consumption time in roads, and the random degree model is used to judge whether the roads are blocked. Then the shortest evacuation route is obtained based on Dijstra algorithm. Gongming Reservoir located in Shenzhen is taken as a case to study. The results show that industrial area Ⅰ is flooded at 2 500 s, and after 5 500 s, most of industrial area Ⅱ is submerged. The Hushan, Loucun Forest and Chaishan are not flooded around industrial area Ⅰand Ⅱ. Based on the above analysis, the optimal evacuation scheme is determined.

Experimental Study on the Dam-Break Hydrographs at the Gate Location

When studying the dam-break flow phenomenon,the basic hydrodynamic features of the dam-break flow at the gate location should be verified primarily.In this study,laboratory experiments were performed in a rectangular and horizontal flume with the same initial water head setting on the dry and wet downstream bed conditions.Water surface elevation was extracted through image analysis and validated by comparing with the data measured using a wave gauge.Temporal variation of the water surface elevation at the gate location,quantified in terms of high-speed video recorded images,can be divided into three stages,the sharp decreasing stage,the relatively steady stage,and the gradually decreasing stage.Applicability of several classic analytical solutions of the dam-break problem at the gate location was validated using present experimental data.Ritter's solution is effective for the dry bed condition while Stoker's solution could be applied to the wet bed case,and both are only applicable during the steady stage.Lin' solution reproduces the gate-site hydrographs well during both the relatively steady and the gradually decreasing stages,especially for the condition under which the down-upstream water depth ratio is smaller than 0.138.

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.

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.

Study of total variation diminishing (TVD) slope limiters in dam-break flow simulation

A two-dimensional （2D） dam-break flow numerical model was developed based on the finite-volume total variation diminishing （TVD） and monotone upstream-centered scheme for conservation laws （MUSCL）-Hancock scheme, which has second-order accuracy in both time and space. A Harten-Lax-van Leer-contact （HLLC） approximate Riemann solver was used to evaluate fluxes. The TVD MUSCL-Hancock numerical scheme utilizes slope limiters, such as the minmod, double minmod, superbee, van Albada, and van Leer limiters, to prevent spurious oscillations and maintain monotonicity near discontinuities. A comparative study of the impact of various slope limiters on the accuracy of the numerical flow model was conducted with several dam-break examples including wet and dry bed cases. The numerical results of the superbee and double minmod limiters agree better with the theoretical solution and have higher accuracy than other limiters in one-dimensional （1D） space. The ratio of the downstream water depth to the upstream water depth was used to select the proper slope limiter. For the 2D numerical model, the superbee limiter should not be used, owing to significant numerical dispersion.

城市水库溃坝耦合极端降雨前后模拟对比分析

以探究城市上游水库溃坝耦合极端降雨及仅溃坝两种工况下洪水灾害程度的变化为目标,确定郑州市二七区、中原区、金水区及惠济区为研究区域,结合芝加哥雨型公式及长历时暴雨强度公式推求设计暴雨,依托“7·20郑州特大暴雨”实测数据推求研究区域产汇流特征,建立极端降雨与洪水过程的响应关系。在此基础上,以郑州市二七区尖岗水库为研究对象,构建两种工况下二维溃坝洪水演进模型,对水情信息进行对比分析。结果表明:相比于仅溃坝,降雨耦合溃坝工况下淹没范围更大,其中增加的淹没范围内各点的最大水深、流速、洪水严重性程度均处于较低等级。在河道两侧,两种工况下各等级水深、流速及洪水严重性淹没范围大致相同;地势较低的区域降雨耦合溃坝工况下淹没水深和流速增大较多,易涝点水深增加约0.6 m,选取断面水面线升高约1 m。

基于CEL方法的带结构物溃坝流数值模拟

随着全球气温的不断升高,极端天气现象增多,溃坝现象成为水工结构方向需要重点防范的自然灾害。溃坝流的发展会对下游结构物产生巨大的破坏作用,因此深入探究溃坝流现象的水动力学特性以及溃坝流下游结构物所受到的冲击力影响势在必行。本文针对带有结构物的溃坝流动问题现象,采用CEL方法(Coupled Eulerian-Lagrangian analysis,耦合欧拉-拉格朗日方法),对于溃坝流的发展、演化进行了数值模拟。研究结果表明,模拟结果中溃坝流自由表面及流体尖端时程变化与文献试验结果拟合较好,CEL方法可应用在溃坝流现象的模拟,并且具有较好的准确性和精度。结构物材料性能、高度及与结构物与溃坝流间的位置,均会影响溃坝流的运动状态以及结构物的受力变形状态,在不考虑结构物破坏的情况下,提高结构物弹性模量,增加结构物高度,以及减小结构物与溃坝流间的距离均可提高结构物对于溃坝流的阻挡作用。

极端工况下梯级水电枢纽群大坝连溃情景分析

为分析极端工况下梯级水电枢纽群可能出现的最大连溃灾难,提出了两种溃坝工况:一种针对土石坝或受损及病险的重力坝与拱坝,认为大坝漫顶即溃,且全溃;另一种针对正常运行的重力坝和拱坝情形,认为漫顶不溃,采用堰流公式计算下泄水量,并采用洪峰展平公式计算两种工况下溃坝洪水向下游的演进过程。以大渡河流域13座梯级水电枢纽群为例展开分析。结果表明:下游泸定、大岗山、瀑布沟大坝能够抵挡和削弱一部分上游溃坝洪水,是流域的控制性大坝;工况一(漫顶即溃)发生时的溃坝影响比工况二(漫顶不溃)发生时的影响更为恶劣,水量更大,持续时间更长,产生的峰顶流量更大;漫顶即溃时的下泄流量过程与库容、距离以及起溃水位有关,而漫顶不溃坝时,下泄流量过程与库容、起溃水位、坝型选择密切相关;“上小下大型”梯级电站布局和以重力坝、拱坝为主的坝型选择将有效降低梯级水电枢纽群发生大坝连溃的风险。研究成果可为梯级水电枢纽群大坝连溃风险分析提供理论支撑及技术参考。