Numerical Simulation of Dam Breaking Using Smoothed Particle Hydrodynamics and Viscosity Behavior

Smoothed particle hydrodynamics (SPH) is a Lagrangian meshless particle method. It is one of the best method for simulating violent free surface flows in fluids and solving large fluid deformations. Dam breaking is a typical example of these problems. The basis of SPH was reviewed, including some techniques for governing equation resolution, such as the stepping method and the boundary handling method. Then numerical results of a dam breaking simulation were discussed, and the benefits of concepts like artificial viscosity and position correction were analyzed in detail. When compared with dam breaking simulated by the volume of fluid (VOF) method, the wave profile generated by SPH had good agreement, but the pressure had only reasonable agreement. Improving pressure results is clearly an important next step for research.

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.

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.

Analysis of the wave propagation and the flow char-acteristics around cylinder due to dam break water

Dam breaks are easily triggered by heavy rains due to extreme weather such as typhoons,causing serious economic losses and casualties.Through the investigation of Chaoshan coastal zone,it is found that there have been dam breaks caused by geological disasters.In the design and management of water conservancy project,it is very important to analyze the effect of disastrous flow caused by dam break on the building.In this paper,the effect of the dam break flow on the cylinder is simulated numerically by taking the water body with initial velocity as the dam break flow,and the flow characteristics around the cylinder and the water body are analyzed.Numerical model adopted the Renault Average Navier-Stokes(RANS)model and volume of fluid(VOF)method to analyze the evolution of free water surface.It is found that there are different patterns of water movement in the process of dam break resulting in the creation of several isolated convex hull forms of dam-break waves on the stationary water surface,which causes longer disturbances in the water near the cylinder and makes the cylinder more vulnerable to fatigue damage.The increase of the height of the dam breaking water will lead to the increase of the hydrodynamic force on the pipeline.This study has guiding significance for the study of dam break and dam body design in water conservancy projects.

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.

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.

Safety of Dams： A Pathological Approach of Qualitative and Quantitative Risks

Dams are critical and essential elements in any infrastructure and, in front of accidents occurred in many countries, it is extremely important to know the risk of these structures. Inserted in this context, it was found in the technical literature, methods and tools capable of measuring the exposure value by means of indicators. In the study, the highlights were 12 methods of qualitative, semiquantitative and quantitative risk analysis, representing an overview of risk analysis methods available in the literature with potential use in dams, that it has been done into electronic spreadsheets. The case study is performed on a sample of concrete dam and earth/rockfill built and operated by Eletrobr^s Furnas Company, supported by documentary research, projects, field inspections and interviews with experts. After applying the methods and the analysis thereof, has been prepared the Eletrobras Fumas dam risk analysis method which is characterized by adapting the criteria analyzed to the reality of the company＇s dams and it was also performed the portfolio risk analysis of 18 dams. In spite of the variety and subjectivity of qualitative and semiquantitative methods, the results show that they tend to converge on the analysis of dam based on risk. The application methodology demonstrates the feasibility assessment stage, covering the preliminary analysis for portfolio dams, followed by formal and individual risk analyzes for the most critical structures. These results confirm the applicability of risk analysis techniques, contributing to the consolidation of this toot as fundamental in the dam safety.

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 catastrophic natural disaster chain of typhoon-rainstorm-landslide-barrier lake-flooding in Zhejiang Province, China

On August 10,2019,due to the effect of a rainstorm caused by Super Typhoon Lekima,a landslide occurred in Shanzao Village,China.It blocked the Shanzao stream,forming a barrier lake,and then the barrier lake burst.This is a rare natural disaster chain of typhoon-rainstorm-landslide-barrier lake-flooding.This study was built on field surveys,satellite image interpretation,the digital elevation model(DEM),engineering geological analysis and empirical regression.The purpose was to reveal the characteristics and causes of the landslide,the features and formation process of the barrier lake and the dam break flooding discharge.The results show that the volume of the landslide deposit is approximately 2.4×105 m3.The burst mode of the landslide dam is overtopping,which took only 22 minutes from the formation of the landslide dam to its overtopping.The dam-break peak flow was 1353 m3/s,and the average velocity was 2.8–3.0 m/s.This study shows that the strongly weathered rock and soil slope has low strength and high permeability under the condition of heavy rainfall,which reminds us the high risk of landslides and the importance of accurate early warning of landslides under heavy rainfalls in densely populated areas of Southeast China,as well as the severity of the disaster chain of typhoon-rainstorm-landslide-barrier lake-flooding.

A CIP-Based Numerical Simulation of Free Surface Flow Related to Freak Waves

An enhanced numerical model for simulating two-dimensional incompressible viscous flow with distorted free surface is reported. The numerical simulation is carried out through the CIP （Constrained Interpolation Profile）-based method, which is described in the paper. A more accurate interface capturing scheme, the VOF/WLIC scheme （VOF：Volume-of-Fluid;WLIC：weighed line interface calculation）, is adopted as the interface capturing method. To assess the developed algorithm and its versatility, a selection of test problems are examined, i.e. the square wave propagation, the Zalesak’s rigid body rotation, dam breaking problem with and without obstacles, wave sloshing in an excited wave tank and interaction between extreme waves and a floating body. Excellent agreements are obtained when numerical results are compared with available analytical, experimental, and other numerical results. These examples demonstrate that the use of the VOF/WLIC scheme in the free surface capturing makes better results and also the proposed CIP-based model is capable of predicting the freak wave-related phenomena.

Viscous flow with free surface motion by least square finite element method

A method for simulation of free surface problems is presented. Based on the viscous incompressible Navier-Stokes equations, space discretization of the flow is obtained by the least square finite element method. The time evolution is obtained by the finite difference method. Lagrangian description is used to track the free surface. The results are compared with the experimental dam break results, including water collapse in a 2D rectangular section and in a 3D cylinder section. A good agreement is achieved for the distance of surge front as well as the height of the residual column.

EXPERIMENTAL STUDY OF DAM-BREAK FLOW IN CASCADE RESERVOIRS WITH STEEP BOTTOM SLOPE

Dam break can cause a significant disaster in the downstream, especially, in a valley with cascade reservoirs, which would aggravate the disaster extent. The experimental studies of the dam-break flow of cascade reservoirs are few and far between at the present, Most of related studies concern the failure of a single dam.. This article presents an experimental study of the characteristics of an instantly filled dam-break flow of cascade reservoirs in a rectangular glass flume with a steep bottom slope. A new method was used to simulate the sudden collapse of the dam. A series of sensors for automatic water-levels were deployed to record the rapid water depth fluctuation. The experimental results show that, the ratio of the initial water depth of the downstream reservoir to that of the upstream reservoir would greatly affect the flood peak water depth in the downstream reservoir area and in the stream channel behind the downstream dam, while the influence of the dam spacing is insignificant. In addition, the comparison between the single reservoir and the cascade reservoirs shows some difference in the dam-break flow pattern and the stage hydrograph at the corresponding gauging points.

Numerical and experimental investigation of three-dimensionality in the dam-break flow against a vertical wall

The three-dimensionality extent of the dam break flow over a vertical wall is investigated numerically and experimentally in this paper. The numerical method is based on Reynolds averaged Navier-Stokes（RANS） equation that describes the three-dimensional incompressible turbulent flow. The free surface is captured by using the unstructured multi-dimensional interface capturing（UMTHINC） scheme. The equations are discretized on 2-D and 3-D unstructured grids using finite volume method. The numerical simulations are compared with newly conducted experiment with emphasis on the effect of three-dimensionality on both free surface and impact pressure. The comparison between the numerical and experimental results shows good agreement. Furthermore, the results also show that 3-D motion of the flow originates at the moment of impact at the lower corners of the impact wall and propagates to the inner region as time advances. The origin of the three-dimensionality is found to be the turbulence development as well as the relative velocity between the side wall region and the inner region of the wave front at the moment of impact.

A virtual geographic environment for dynamic simulation and analysis of tailings dam failure

It is of great significance for disaster prevention and mitigation to carry out disaster simulations for dam failure accidents in advance,but at present,there are few professional systems for disaster simulations of tailings dams.In this paper,we focused on the construction of a virtual geographic environment(VGE)system that provides an effective tool for visualizing the dam-break process of a tailings pond.The dam-break numerical model of the tailings dam based on computational fluid dynamics(CFD)was integrated into the VGE system.The infrastructure of the VGE was supported by a 3-D geographic information system(GIS)with a user-friendly interface for the initiation,visualization,and analysis of the dynamic process of tailings dam failure.Key technologies,including the integration of numerical models,rendering of large-scale scenes,and optimizations of disaster simulation and visualization,were discussed in detail.In the prototype system,information on the run-out path,travel distance,etc.can be obtained to visually describe the flow motion released by two dam failure cases.The simulation results showed that the VGE can be used for the multidimensional,dynamic and interactive visualization of dam-break disasters,and can also be useful for assessing the risk associated with tailings dams.

LATTICE BOLTZMANN SIMULATIONS OF TRIAGULAR CAVITY FLOW AND FREE-SURFACE PROBLEMS

The Lattice Boltzmann Method （LBM） was investigated to solve triangular cavity flow and free-surface problems in hydraulic dynamics. Some cases of triangular cavity flow and backward step flow were simulated to show the efficiency and stability of this method. Two-dimensional partial dam breaking problem and the propagation and diffraction of dam-break wave around rectangular and circular cylinder were numerically studied successfully. Excellent agreement was obtained between numerical predictions and available results.

Numerical simulation of 3-D free surface flows by overlapping MPS

An overlapping moving particle semi-implicit （MPS） method is applied for 3-D free surface flows based on our in-house particle solver MLParticle-SJTU. In this method, the coarse particles are distributed in the whole domain and the fine particles are distributed in the local region of interest at the same time. With the fine particles being generated and removed dynamically, an algorithm of generating particles based on the 3-D overlapping volume is developed. Then, a 3-D dam break flow with an obstacle is simulated to validate the overlapping MPS. The qualitative comparison among experimental data and the results obtained by the VOF and the MPS shows that the shape of the free surface obtained by the overlapping MPS is more accurate than that obtained by the UNI-coarse and close to that obtained by the UNI-fine in the overlapping domain. In addition, the water height and the impact pressure at Pi are also in an overall agreement with experimental data. Finally, the CPU time required by the overlapping MPS is about half of that required by the UNl-fine.

AN UNSTRUCTURED FINITE-VOLUME ALGORITHM FOR NONLINEAR TWO-DIMENSIOAL SHALLOW WATER EQUATION

An unstructured finite-volume numerical algorithm was presented for solution of the two-dimensional shallow water equations, based on triangular or arbitrary quadrilateral meshes. The Roe type approximate Riemann solver was used to the system. A second-order TVD scheme with the van Leer limiter was used in the space discretization and a two-step Runge-Kutta approach was used in the time discretization. An upwind, as opposed to a pointwise, treatment of the slope source terms was adopted and the semi-implicit treatment was used for the friction source terms. Verification for two-dimension dam-break problems are carried out by comparing the present results with others and very good agreement is shown.

A NEW HIGH ORDER SCHEME FOR CONVECTION EQUATION AND ITS APPLICATION

With the undetermined coefficient method, a new high order scheme for convection equation named as HAUC was worke d out. The accuracy of the new scheme was analyzed by comparing the computed res ults with the exact solutions of 1-D pure convection equation and nonlinear con vection equation. The effectiveness of the HAUC was also examined by comparing w ith the results obtained by other schemes. Finally, it was applied to simulate 1 -D dam break flow with differenct ratios of initial upstream water depth to do wnstream one. The results show that the scheme has the ability to simulate both undular bores and moving hydraulic jumps.

A HIGH RESOLUTION FINITE VOLUME METHOD FOR SOLVING SHALLOW WATER EQUATIONS

A high resolution finite volume numerical method for solving the shallow water equations is developed in this paper. In order to extend finite difference TVD scheme to finite volume method, a new geometry and topology of control bodies is defined by considering the corresponding relationships between nodes and elements. This solver is implemented on arbitrary quadrilateral meshes and their satellite elements, and based on a second order hybrid type of TVD scheme in space discretization and a two step Runge Kutta method in time discretization. Then it is used to deal with two typical dam break problems and very satisfactory results are obtained comparied with other numerical solutions. It can be considered as an efficient implement for the computation of shallow water problems, especially concerning those having discontinuities, subcritical and supercritical flows and complex geometries.

Numerical modeling of sediment transport based on unsteady and steady flows by incompressible smoothed particle hydrodynamics method

The purpose of the present paper is to introduce a simple two-part multi-phase model for the sediment transport problems based on the incompressible smoothed particle hydrodynamics（ISPH） method. The proposed model simulates the movement of sediment particles in two parts. The sediment particles are classified into three categories, including the motionless particles, moving particles behave like a rigid body, and moving particles with a pseudo fluid behavior. The criterion for the classification of sediment particles is the Bingham rheological model. Verification of the present model is performed by simulation of the dam break waves on movable beds with different conditions and the bed scouring under steady flow condition. Comparison of the present model results, the experimental data and available numerical results show that it has good ability to simulate flow pattern and sediment transport.