scholarly journals Simulation of Mosul Dam Break Using Finite Volume Method

2020 ◽  
Vol 10 (2) ◽  
pp. 10-20
Author(s):  
Ayhan H. Saleem ◽  
Jowhar R. Mohammad

Mosul dam is an earth-fill embankment located north of Iraq on the Tigris River forming a reservoir with 11.11 km3 water storage capacity which is the largest dam in the country. The dam is built on a rock bed foundation, in which the dissolution process is dynamic in the zone where gypsum and anhydrite layers present. During the construction development seepage locations were found in the dam foundation and the grouting process is in progress until now to control this problem. Therefore, the possibility of the Mosul dam break is highlighted by previous studies. In this research, a FORTRAN code based on the finite volume method is modified to solve the two-dimensional shallow water equations and simulating the Mosul dam break. The computational domain discretized using unstructured triangular mesh. The solver applied Harten lax van leer with contact (HLLC) wave approximate Riemann solver to calculate the cell interface fluxes, and the semi-implicit scheme employed to solve the friction source term. The numerical scheme applied to two benchmark test cases, and the results showed that the presented model was robust and accurate especially in handling wet/dry beds, mixed flow regimes, discontinuities, negative water depths, and complex topography. The results of this study demonstrate that flood waves may reach the center of Mosul city in < 6 h and water depth may rise to 34 m after 7 h of Mosul dam breaking. Finally, the simulation results of the Mosul dam break were used to prepare an emergency action plan.

2020 ◽  
Vol 50 (3) ◽  
pp. 287-302
Author(s):  
Róbert ČUNDERLÍK ◽  
Matej MEDĽA ◽  
Karol MIKULA

The paper presents local quasigeoid modelling in Slovakia using the finite volume method (FVM). FVM is used to solve numerically the fixed gravimetric boundary value problem (FGBVP) on a 3D unstructured mesh created above the real Earth's surface. Terrestrial gravimetric measurements as input data represent the oblique derivative boundary conditions on the Earth's topography. To handle such oblique derivative problem, its tangential components are considered as surface advection terms regularized by a surface diffusion. The FVM numerical solution is fixed to the GOCE-based satellite-only geopotential model on the upper boundary at the altitude of 230 km. The horizontal resolution of the 3D computational domain is 0.002 × 0.002 deg and its discretization in the radial direction is changing with altitude. The created unstructured 3D mesh of finite volumes consists of 454,577,577 unknowns. The FVM numerical solution of FGBVP on such a detailed mesh leads to large-scale parallel computations requiring 245 GB of internal memory. It results in the disturbing potential obtained in the whole 3D computational domain. Its values on the discretized Earth's surface are transformed into the local quasigeoid model that is tested at 404 GNSS/levelling benchmarks. The standard deviation of residuals is 2.8 cm and decreases to 2.6 cm after removing 9 identified outliers. It indicates high accuracy of the obtained FVM-based local quasigeoid model in Slovakia.


2012 ◽  
Vol 2012 ◽  
pp. 1-14 ◽  
Author(s):  
Szu-Hsien Peng

The purpose of this study is to model the flow movement in an idealized dam-break configuration. One-dimensional and two-dimensional motion of a shallow flow over a rigid inclined bed is considered. The resulting shallow water equations are solved by finite volumes using the Roe and HLL schemes. At first, the one-dimensional model is considered in the development process. With conservative finite volume method, splitting is applied to manage the combination of hyperbolic term and source term of the shallow water equation and then to promote 1D to 2D. The simulations are validated by the comparison with flume experiments. Unsteady dam-break flow movement is found to be reasonably well captured by the model. The proposed concept could be further developed to the numerical calculation of non-Newtonian fluid or multilayers fluid flow.


Author(s):  
Yanbo Huang ◽  
Shanshan Li ◽  
Zhenhai Pan

Abstract Fluid-structure interaction (FSI) is an important fundamental problem with wide scientific and engineering applications. The immersed boundary method has proved to be an effective way to model the interaction between a moving solid and its surrounding fluid. In this study, a novel modeling approach based on the coupled immersed-boundary and finite-volume method is proposed to simulate fluid-structure interaction problems. With this approach, the whole computational domain is treated as fluid and discretized by only one set of Eulerian grids. The computational domain is divided into solid parts and fluid parts. A goal velocity is locally determined in each cell inside the solid part. At the same time, the hydrodynamic force exerted on the solid structure is calculated by integrating along the faces between the solid cells and fluid cells. In this way, the interaction between the solid and fluid is solved explicitly and the costly information transfer between Lagranian grids and Eulerian grids is avoided. The interface is sharply restricted into one single grid width throughout the iterations. The proposed modeling approach is validated by conducting several classic numerical experiments, including flow past static and freely rotatable square cylinders, and sedimentation of an ellipsoid in finite space. Throughout the three numerical experiments, satisfying agreements with literatures have been obtained, which demonstrate that the proposed modeling approach is accurate and robust for simulating FSI problems.


Author(s):  
Valery Ponyavin ◽  
Roald Akberov ◽  
Yitung Chen ◽  
Hsuan-Tsung Hsieh ◽  
Darrell W. Pepper

The calculation of gas flow during the motion of a projectile in the gun barrel is a complicated computational task due of the presence of numerous factors, such as nonisothermicity, turbulence, changes in the shape of the computational domain with time, etc. In this study, an attempt to calculate the characteristics of gas flow around a projectile during the motion of the projectile in the gun barrel is undertaken. The flow is considered axisymmetrical, nonstationary, nonisothermal, compressible, and turbulent. For calculating the flow around the projectile, the finite volume method was employed. During the motion of the projectile, the flow pattern behind it changed from subsonic to supersonic. The results of the calculations are represented in figures depicting the flow at different moments of time. The figures show the fields of velocity, pressure and density, as well as the appearance of shock waves inside the gun barrel at subsonic and supersonic speeds.


2021 ◽  
Author(s):  
Farid Boushaba ◽  
Salah Daoudi ◽  
Ahmed Yachouti ◽  
Youssef Regad

Abstract This paper presents numerical solvers, based on the finite volume method. This scheme solves dam break problems on the dry bottom in 2D configuration. The difficulty of the simulation of this type of problem lies in the propagation of shocks on the dry bottom. The equation model used is the shallow water equations written in conservative form. The scheme used is second order in space and time. The method is modified to treat dry bottoms. The validity of the method is demonstrated over the dam break example. A comparison with finite elements shows the weakness and robustness of each method.


2013 ◽  
Vol 2013 ◽  
pp. 1-12 ◽  
Author(s):  
Hamid Reza Vosoughifar ◽  
Azam Dolatshah ◽  
Seyed Kazem Sadat Shokouhi

This paper was concerned to simulate both wet and dry bed dam break problems. A high-resolution finite volume method (FVM) was employed to solve the one-dimensional (1D) and two-dimensional (2D) shallow water equations (SWEs) using an unstructured Voronoi mesh grid. In this attempt, the robust local Lax-Friedrichs (LLxF) scheme was used for the calculating of the numerical flux at cells interfaces. The model named V-Break was run under the asymmetry partial and circular dam break conditions and then verified by comparing the model outputs with the documented results. Due to a precise agreement between those output and documented results, the V-Break could be considered as a reliable method for dealing with shallow water (SW) and shock problems, especially those having discontinuities. In addition, statistical observations indicated a good conformity between the V-Break and analytical results clearly.


Sign in / Sign up

Export Citation Format

Share Document