scholarly journals Simulation of Dam Breaks on Dry Bed Using Finite Volume Roe-TVD Method

Hydrology ◽  
2021 ◽  
Vol 8 (2) ◽  
pp. 88
Author(s):  
Ebrahim Alamatian ◽  
Sara Dadar ◽  
Bojan Đurin

Dams are one of the most important hydraulic structures. In view of unrecoverable damages occurring after a dam failure, analyzing a dams’ break is necessary. In this study, a dam located in Iran is considered. According to adjacent tourist and entertainment zones, the breaking of the dam could lead to severe problems for the area and bridges downstream of the river. To investigate the issue, a numerical FORTRAN code based on the 2D finite volume Roe-TVD method on a fixed bed is provided to assess the effects of the dam break. Turbulence terms and dry bed conditions were considered in the code. A numerical wave tank (NWT) with a triangular barrier in the bed was numerically modeled and compared with analytical models to verify the capability of the code. Comparing numerical, experimental and analytical results showed that estimated water level and mass conservation in the numerical model is in good agreement with the experimental data and analytical solutions. The 2D approach used has reduced the cost of computing compared to a 3D approach while obtaining accurate results. The code is finally applied to a full-scale dam-break flood. Six KM of the natural river downstream of the dam, including two bridges, B1 and B2, is considered. Flood flow hydrographs and water level variations at bridges B1 and B2 are presented. The results denoted that bridges B1 and B2 will be flooded after 12 and 21 min, respectively, and are at risk of the potential break. Thus, it is necessary to announce and possibly evacuate the resort area alongside the dam in order to decrease losses.

2019 ◽  
Vol 27 (1) ◽  
pp. 344-353
Author(s):  
Abdul-Hassan K. Al-Shukur ◽  
Ranya Badea’ Mahmoud

One of the most common type of embankment dam failure is the dam-break due to overtopping. In this study, the finite elements method has been used to analyze seepage and limit equilibrium method to study stability of the body of an earthfill dam during the flood condition. For this purpose, the software Geostudio 2012 is used through its subprograms SEEP/W and SLOPE/W. Al-Adhaim dam in Iraq has been chosen to analysis the 5 days of flood. It was found that the water flux of seepage during the flood reaches about 8.772*10-5. m3/sec when the water level 146.5 m at 2nd day. Seepage through the embankment at maximum water level increased by 55.1 % from maximum water level. It was concluded that the factor of safety against sliding in downstream side decrease with increasing water level and vice versa. It was also concluded that the deposits are getting more critical stability during the conditions of flood when the factor of safety value reaches 1.219 at 2nd day.


2018 ◽  
Author(s):  
Alfredo L. Aretxabaleta ◽  
Neil K. Ganju ◽  
Zafer Defne ◽  
Richard P. Signell

Abstract. Water level in semi-enclosed bays, landward of barrier islands, is mainly driven by offshore sea level fluctuations that are modulated by bay geometry and bathymetry, causing spatial variability in the ensuing response (transfer). Local wind setup can have a secondary role that depends on wind speed, fetch, and relative orientation of the wind direction and the bay. Inlet geometry and bathymetry primarily regulate the magnitude of the transfer between open ocean and bay. Tides and short-period offshore oscillations are more damped in the bays than longer-lasting offshore fluctuations, such as storm surge and sea level rise. We compare observed and modeled water levels at stations in a mid-Atlantic bay (Barnegat Bay) with offshore water level proxies. Observed water levels in Barnegat Bay are compared and combined with model results from the Coupled Ocean–Atmosphere–Wave–Sediment Transport (COAWST) modeling system to evaluate the spatial structure of the water level transfer. Analytical models based on the dimensional characteristics of the bay are used to combine the observed data and the numerical model results in a physically consistent approach. Model water level transfers match observed values at locations inside the Bay in the storm frequency band (transfers ranging from 70–100 %) and tidal frequencies (10–55 %). The contribution of frequency-dependent local setup caused by wind acting along the bay is also considered. The approach provides transfer estimates for locations inside the Bay where observations were not available resulting in a complete spatial characterization. The approach allows for the study of the Bay response to alternative forcing scenarios (landscape changes, future storms, and rising sea level). Detailed spatial estimates of water level transfer can inform decisions on inlet management and contribute to the assessment of current and future flooding hazard in back-barrier bays and along mainland shorelines.


Water ◽  
2020 ◽  
Vol 12 (6) ◽  
pp. 1639
Author(s):  
Abdelkrim Aharmouch ◽  
Brahim Amaziane ◽  
Mustapha El Ossmani ◽  
Khadija Talali

We present a numerical framework for efficiently simulating seawater flow in coastal aquifers using a finite volume method. The mathematical model consists of coupled and nonlinear partial differential equations. Difficulties arise from the nonlinear structure of the system and the complexity of natural fields, which results in complex aquifer geometries and heterogeneity in the hydraulic parameters. When numerically solving such a model, due to the mentioned feature, attempts to explicitly perform the time integration result in an excessively restricted stability condition on time step. An implicit method, which calculates the flow dynamics at each time step, is needed to overcome the stability problem of the time integration and mass conservation. A fully implicit finite volume scheme is developed to discretize the coupled system that allows the use of much longer time steps than explicit schemes. We have developed and implemented this scheme in a new module in the context of the open source platform DuMu X . The accuracy and effectiveness of this new module are demonstrated through numerical investigation for simulating the displacement of the sharp interface between saltwater and freshwater in groundwater flow. Lastly, numerical results of a realistic test case are presented to prove the efficiency and the performance of the method.


2021 ◽  
Vol 11 (12) ◽  
pp. 5638
Author(s):  
Selahattin Kocaman ◽  
Stefania Evangelista ◽  
Hasan Guzel ◽  
Kaan Dal ◽  
Ada Yilmaz ◽  
...  

Dam-break flood waves represent a severe threat to people and properties located in downstream regions. Although dam failure has been among the main subjects investigated in academia, little effort has been made toward investigating wave propagation under the influence of tailwater depth. This work presents three-dimensional (3D) numerical simulations of laboratory experiments of dam-breaks with tailwater performed at the Laboratory of Hydraulics of Iskenderun Technical University, Turkey. The dam-break wave was generated by the instantaneous removal of a sluice gate positioned at the center of a transversal wall forming the reservoir. Specifically, in order to understand the influence of tailwater level on wave propagation, three tests were conducted under the conditions of dry and wet downstream bottom with two different tailwater depths, respectively. The present research analyzes the propagation of the positive and negative wave originated by the dam-break, as well as the wave reflection against the channel’s downstream closed boundary. Digital image processing was used to track water surface patterns, and ultrasonic sensors were positioned at five different locations along the channel in order to obtain water stage hydrographs. Laboratory measurements were compared against the numerical results obtained through FLOW-3D commercial software, solving the 3D Reynolds-Averaged Navier–Stokes (RANS) with the k-ε turbulence model for closure, and Shallow Water Equations (SWEs). The comparison achieved a reasonable agreement with both numerical models, although the RANS showed in general, as expected, a better performance.


2019 ◽  
Vol 17 (08) ◽  
pp. 1950053
Author(s):  
Yanyun Wu ◽  
Liquan Mei ◽  
Meilan Qiu ◽  
Yuchuan Chu

We present a stabilized finite volume element method for the coupled Stokes–Darcy problem with the lowest order [Formula: see text] element for the Stokes region and [Formula: see text] element for the Darcy region. Based on adding a jump term of discrete pressure to the approximation equation, a discrete inf-sup condition is established for the proposed method. The optimal error estimates in the [Formula: see text]-norm for the velocity and piezometric head and in the [Formula: see text]-norm for the pressure are proved. And they are also verified through some numerical experiments. Two figures are given to show the full comparison for the local mass conservation between the proposed method and the stabilized finite element method. And this method can also be computed directly in the irregular domain according to the last experiment.


2014 ◽  
Vol 7 (3) ◽  
pp. 909-929 ◽  
Author(s):  
J. Thuburn ◽  
C. J. Cotter ◽  
T. Dubos

Abstract. A new algorithm is presented for the solution of the shallow water equations on quasi-uniform spherical grids. It combines a mimetic finite volume spatial discretization with a Crank–Nicolson time discretization of fast waves and an accurate and conservative forward-in-time advection scheme for mass and potential vorticity (PV). The algorithm is implemented and tested on two families of grids: hexagonal–icosahedral Voronoi grids, and modified equiangular cubed-sphere grids. Results of a variety of tests are presented, including convergence of the discrete scalar Laplacian and Coriolis operators, advection, solid body rotation, flow over an isolated mountain, and a barotropically unstable jet. The results confirm a number of desirable properties for which the scheme was designed: exact mass conservation, very good available energy and potential enstrophy conservation, consistent mass, PV and tracer transport, and good preservation of balance including vanishing ∇ × ∇, steady geostrophic modes, and accurate PV advection. The scheme is stable for large wave Courant numbers and advective Courant numbers up to about 1. In the most idealized tests the overall accuracy of the scheme appears to be limited by the accuracy of the Coriolis and other mimetic spatial operators, particularly on the cubed-sphere grid. On the hexagonal grid there is no evidence for damaging effects of computational Rossby modes, despite attempts to force them explicitly.


2018 ◽  
Vol 147 ◽  
pp. 03009 ◽  
Author(s):  
Bagus Pramono Yakti ◽  
Mohammad Bagus Adityawan ◽  
Mohammad Farid ◽  
Yadi Suryadi ◽  
Joko Nugroho ◽  
...  

A dam break induced-flood propagation modeling is needed to reduce the losses of any potential dam failure. On the 25 July 2013, there was a dam break generated flood due to the failure of Way Ela Natural Dam that severely damaged houses and various public facilities. This study simulated the flooding induced by the failure of Way Ela Natural Dam. A two-dimensional (2D) numerical model, HEC-RAS v.5, is used to simulate the overland flow. The dam failure itself is simulated using HECHMSv.4. The results of this study, the flood inundation, flood depth, and flood arrival time are verified by using available secondary data. These informations are very important to propose mitigation plans with respect to possible dam break in the future.


2019 ◽  
Vol 8 (6) ◽  
pp. 250 ◽  
Author(s):  
Raffaele Albano ◽  
Leonardo Mancusi ◽  
Jan Adamowski ◽  
Andrea Cantisani ◽  
Aurelia Sole

Mapping the delineation of areas that are flooded due to water control infrastructure failure is a critical issue. Practical difficulties often present challenges to the accurate and effective analysis of dam-break hazard areas. Such studies are expensive, lengthy, and require large volumes of incoming data and refined technical skills. The creation of cost-efficient geospatial tools provides rapid and inexpensive estimates of instantaneous dam-break (due to structural failure) flooded areas that complement, but do not replace, the results of hydrodynamic simulations. The current study implements a Geographic Information System (GIS) based method that can provide useful information regarding the delineation of dam-break flood-prone areas in both data-scarce environments and transboundary regions, in the absence of detailed studies. Moreover, the proposed tool enables, without advanced technical skills, the analysis of a wide number of case studies that support the prioritization of interventions, or, in emergency situations, the simulation of numerous initial hypotheses (e.g., the modification of initial water level/volume in the case of limited dam functionality), without incurring high computational time. The proposed model is based on the commonly available data for masonry dams, i.e., dam geometry (e.g., reservoir capacity, dam height, and crest length), and a Digital Elevation Model. The model allows for rapid and cost-effective dam-break hazard mapping by evaluating three components: (i) the dam-failure discharge hydrograph, (ii) the propagation of the flood, and (iii) the delineation of flood-prone areas. The tool exhibited high accuracy and reliability in the identification of hypothetical dam-break flood-prone areas when compared to the results of traditional hydrodynamic approaches, as applied to a dam in Basilicata (Southern Italy). In particular, the over- and under-estimation rates of the proposed tool, for the San Giuliano dam, Basilicata, were evaluated by comparing its outputs with flood inundation maps that were obtained by two traditional methods whil using a one-dimensional and a two-dimensional propagation model, resulting in a specificity value of roughly 90%. These results confirm that most parts of the flood map were correctly classified as flooded by the proposed GIS model. A sensitivity value of over 75% confirms that several zones were also correctly identified as non-flooded. Moreover, the overall effectiveness and reliability of the proposed model were evaluated, for the Gleno Dam (located in the Central Italian Alps), by comparing the results of literature studies concerning the application of monodimensional numerical models and the extent of the flooded area reconstructed by the available historical information, obtaining an accuracy of around 94%. Finally, the computational efficiency of the proposed tool was tested on a demonstrative application of 250 Italian arch and gravity dams. The results, when carried out using a PC, Pentium Intel Core i5 Processor CPU 3.2 GHz, 8 GB RAM, required about 73 min, showing the potential of such a tool applied to dam-break flood mapping for a large number of dams.


2011 ◽  
Vol 243-249 ◽  
pp. 3877-3882
Author(s):  
Rong Yong Ma ◽  
Xiang Chen ◽  
Lei Lei Yang ◽  
Xiao Qing Zhang

In this paper,the basic situation of Guangxi luocheng Kama reservoir and the potential dangers are introduced firstly, and then,different models often used to estimate the loss of life in dam failure at home and abroad are introduced and analysed.Based on the above, a suitable model for national actual situation is chosen to predict the loss due to this reservoir dam breach,and the severity according to the calculation result of assuming Kama reservoir dam-break is assessed. Finally, the relevant departments’s measures of rescuing and evacuating downstream masses etc. at danger in 2009 are proved to be correct,necessary and timely.


Sign in / Sign up

Export Citation Format

Share Document