scholarly journals 2D Shallow Water Model for Dam Break and Column Interactions

2020 ◽  
Vol 6 (3) ◽  
pp. 237
Author(s):  
Putu Indah Dianti Putri ◽  
Rifqi Fauzan Iskandar ◽  
Mohammad Bagus Adityawan ◽  
Hadi Kardhana ◽  
Dian Indrawati

Dam break causes disastrous effects on the surrounding area, especially at the downstream, therefore, there is a need for accurate and timely predictions of dam break propagation to prevent both property damage and loss of life. This study aimed to determine the movement of dam-break flow in the downstream area by solving the Shallow Water Equations (SWE) or Saint Venant Equations which are based on the conservation of mass and momentum derived from Navier Stokes equation. The model was generated using a finite difference scheme which is the most common and simplest method for dam-break modeling while Forward Time Central Space (FTCS) numerical scheme was applied to simulate two-dimensional SWE. Moreover, the accuracy of the numerical model was checked by comparing its results with the analytic results of one-dimensional cases and a relatively small value of error was found in comparison to the analytic models as indicated with the RMSE values close to 0. The numerical to the two-dimensional models were also compared to a simple dam break in a flume and dam break with column interactions and the wave propagation in both cases was observed to become very close at a certain time.  The model, however, used numerical filter (Hansen) to reduce the oscillations or numerical instability. The simulation and analysis, therefore, showed the ability of the numerical scheme of FTCS to resolve both cases of the simple dam break and dam break with column interactions in the Two-dimensional Shallow Water. 

2021 ◽  
Vol 2021 ◽  
pp. 1-18
Author(s):  
Nugool Sataporn ◽  
Worasait Suwannik ◽  
Montri Maleewong

Compute Unified Device Architecture (CUDA) implementations are presented of a well-balanced finite volume method for solving a shallow water model. The CUDA platform allows programs to run parallel on GPU. Four versions of the CUDA algorithm are presented in addition to a CPU implementation. Each version is improved from the previous one. We present the following techniques for optimizing a CUDA program: limiting register usage, changing the global memory access pattern, and using loop unroll. The accuracy of all programs is investigated in 3 test cases: a circular dam break on a dry bed, a circular dam break on a wet bed, and a dam break flow over three humps. The last parallel version shows 3.84x speedup over the first CUDA implementation. We use our program to simulate a real-world problem based on an assumed partial breakage of the Srinakarin Dam located in Kanchanaburi province, Thailand. The simulation shows that the strong interaction between massive water flows and bottom elevations under wet and dry conditions is well captured by the well-balanced scheme, while the optimized parallel program produces a 57.32x speedup over the serial version.


2009 ◽  
Vol 6 (6) ◽  
pp. 6759-6793 ◽  
Author(s):  
C. Biscarini ◽  
S. Di Francesco ◽  
P. Manciola

Abstract. This paper presents numerical simulations of free surface flows induced by a dam break comparing the shallow water approach to fully three-dimensional simulations. The latter are based on the solution of the complete set of Reynolds-Averaged Navier-Stokes (RANS) equations coupled to the Volume of Fluid (VOF) method. The methods assessment and comparison are carried out on a dam break over a flat bed without friction and a dam break over a triangular bottom sill. Experimental and numerical literature data are compared to present results. The results demonstrate that the shallow water approach loses some three-dimensional phenomena, which may have a great impact when evaluating the downstream wave propagation. In particular, water wave celerity and water depth profiles could be underestimated due to the incorrect shallow water idealization that neglects the three-dimensional aspects due to the gravity force, especially during the first time steps of the motion.


2010 ◽  
Vol 663 ◽  
pp. 456-477 ◽  
Author(s):  
A. FERRARI ◽  
L. FRACCAROLLO ◽  
M. DUMBSER ◽  
E. F. TORO ◽  
A. ARMANINI

In this paper, the wave propagation on a plane dry bottom after a dam break is analysed. Two mathematical models have been used and compared with each other for simulating such a dam-break scenario. First, the fully three-dimensional Navier–Stokes equations for a weakly compressible fluid have been solved using the new smooth particle hydrodynamics formulation, recently proposed by Ferrari et al. (Comput. Fluids, vol. 38, 2009, p. 1203). Second, the two-dimensional shallow water equations (SWEs) are solved using a third-order weighted essentially non-oscillatory finite-volume scheme. The numerical results are critically compared against the laboratory measurements provided by Fraccarollo & Toro (J. Hydraul. Res., vol. 33, 1995, p. 843). The experimental data provide the temporal evolution of the pressure field, the water depth and the vertical velocity profile at 40 gauges, located in the reservoir and in front of the gate. Our analysis reveals the shortcomings of SWEs in the initial stages of the dam-break phenomenon in reproducing many important flow features of the unsteady free-surface flow: the shallow water model predicts a complex wave structure and a wavy evolution of local free-surface elevations in the reservoir that can be clearly identified to be only model artefacts. However, the quasi-incompressible Navier–Stokes model reproduces well the high gradients in the flow field and predicts the cycles of simultaneous rapid decreasing and frozen stages of the free surface in the tank along with the velocity oscillations. Asymptotically, i.e. for ‘large times’, the shallow water model and the weakly compressible Navier–Stokes model agree well with the experimental data, since the classical SWE assumptions are satisfied only at large times.


2010 ◽  
Vol 14 (4) ◽  
pp. 705-718 ◽  
Author(s):  
C. Biscarini ◽  
S. Di Francesco ◽  
P. Manciola

Abstract. This paper presents numerical simulations of free surface flows induced by a dam break comparing the shallow water approach to fully three-dimensional simulations. The latter are based on the solution of the complete set of Reynolds-Averaged Navier-Stokes (RANS) equations coupled to the Volume of Fluid (VOF) method. The methods assessment and comparison are carried out on a dam break over a flat bed without friction, a dam break over a triangular bottom sill and a dam break flow over a 90° bend. Experimental and numerical literature data are compared to present results. The results demonstrate that the shallow water approach, even if able to sufficiently reproduce the main aspects of the fluid flows, loses some three-dimensional phenomena, due to the incorrect shallow water idealization that neglects the three-dimensional aspects related to the gravity force.


Water ◽  
2021 ◽  
Vol 13 (16) ◽  
pp. 2152
Author(s):  
Gonzalo García-Alén ◽  
Olalla García-Fonte ◽  
Luis Cea ◽  
Luís Pena ◽  
Jerónimo Puertas

2D models based on the shallow water equations are widely used in river hydraulics. However, these models can present deficiencies in those cases in which their intrinsic hypotheses are not fulfilled. One of these cases is in the presence of weirs. In this work we present an experimental dataset including 194 experiments in nine different weirs. The experimental data are compared to the numerical results obtained with a 2D shallow water model in order to quantify the discrepancies that exist due to the non-fulfillment of the hydrostatic pressure hypotheses. The experimental dataset presented can be used for the validation of other modelling approaches.


Water ◽  
2021 ◽  
Vol 13 (3) ◽  
pp. 344
Author(s):  
Le Thi Thu Hien ◽  
Nguyen Van Chien

The aim of this paper was to investigate the ability of some 2D and 3D numerical models to simulate flood waves in the presence of an isolated building or building array in an inundated area. Firstly, the proposed 2D numerical model was based on the finite-volume method (FVM) to solve 2D shallow-water equations (2D-SWEs) on structured mesh. The flux-difference splitting method (FDS) was utilized to obtain an exact mass balance while the Roe scheme was invoked to approximate Riemann problems. Secondly, the 3D commercially available CFD software package was selected, which contained a Flow 3D model with two turbulent models: Reynolds-averaged Navier-Stokes (RANs) with a renormalized group (RNG) and a large-eddy simulation (LES). The numerical results of an impact force on an obstruction due to a dam-break flow showed that a 3D solution was much better than a 2D one. By comparing the 3D numerical force results of an impact force acting on building arrays with the existence experimental data, the influence of velocity-induced force on a dynamic force was quantified by a function of the Froude number and the water depth of the incident wave. Furthermore, we investigated the effect of the initial water stage and dam-break width on the 3D-computed results of the peak value of force intensity.


2011 ◽  
Vol 21 (03) ◽  
pp. 421-457 ◽  
Author(s):  
RAPHAËL DANCHIN ◽  
MARIUS PAICU

Models with a vanishing anisotropic viscosity in the vertical direction are of relevance for the study of turbulent flows in geophysics. This motivates us to study the two-dimensional Boussinesq system with horizontal viscosity in only one equation. In this paper, we focus on the global existence issue for possibly large initial data. We first examine the case where the Navier–Stokes equation with no vertical viscosity is coupled with a transport equation. Second, we consider a coupling between the classical two-dimensional incompressible Euler equation and a transport–diffusion equation with diffusion in the horizontal direction only. For both systems, we construct global weak solutions à la Leray and strong unique solutions for more regular data. Our results rest on the fact that the diffusion acts perpendicularly to the buoyancy force.


Sign in / Sign up

Export Citation Format

Share Document