Abstract: Modeling sediment transport processes by means of numerical approach often poses severe challenges. In this way, a number of techniques have been suggested to solve flow and sediment equations in decoupled, semi-coupled or fully coupled forms. Furthermore, in order to capture flow discontinuities, a number of techniques, like artificial viscosity and shock fitting, have been proposed for solving these equations which are mostly required careful calibration processes. In this research, a numerical scheme for solving shallow water and Exner equations in fully coupled form is presented. First-Order Centered scheme is applied for producing required numerical fluxes and the reconstruction process is carried out toward using Monotonic Upstream Scheme for Conservation Laws to achieve a high order scheme. In order to satisfy C-property of the scheme in presence of bed topography, Surface Gradient Method is proposed. Combining the presented scheme with fourth order Runge-Kutta algorithm for time integration yields a competent numerical scheme. In addition, to handle non-prismatic channels problems, Cartesian Cut Cell Method is employed. A trained Multi-Layer Perceptron Artificial Neural Network which is of Feed Forward Back Propagation (FFBP) type estimates sediment flow discharge in the model rather than usual empirical formulas. Hydrodynamic part of the model is tested for showing its capability in simulation of flow discontinuities, transcritical flows, wetting/drying conditions and non-prismatic channel flows. In this end, dam-break flow onto a locally non-prismatic converging-diverging channel with initially dry bed conditions is modeled. The morphodynamic part of the model is verified simulating dam break on a dry movable bed and bed level variations in an alluvial junction. The results show that the model is capable in capturing the flow discontinuities, solving wetting/drying problems even in non-prismatic channels and presenting proper results for movable bed situations. It can also be deducted that applying Artificial Neural Network, instead of common empirical formulas for estimating sediment flow discharge, leads to more accurate results.
Abstract: A fast finite volume solver for multi-layered shallow
water flows with mass exchange and an erodible bed is developed.
This enables the user to solve a number of complex sediment-based
problems including (but not limited to), dam-break over an erodible
bed, recirculation currents and bed evolution as well as levy and
dyke failure. This research develops methodologies crucial to the
under-standing of multi-sediment fluvial mechanics and waterway
design. In this model mass exchange between the layers is allowed
and, in contrast to previous models, sediment and fluid are able
to transfer between layers. In the current study we use a two-step
finite volume method to avoid the solution of the Riemann problem.
Entrainment and deposition rates are calculated for the first time in
a model of this nature. In the first step the governing equations are
rewritten in a non-conservative form and the intermediate solutions
are calculated using the method of characteristics. In the second stage,
the numerical fluxes are reconstructed in conservative form and are
used to calculate a solution that satisfies the conservation property.
This method is found to be considerably faster than other comparative
finite volume methods, it also exhibits good shock capturing. For most
entrainment and deposition equations a bed level concentration factor
is used. This leads to inaccuracies in both near bed level concentration
and total scour. To account for diffusion, as no vertical velocities
are calculated, a capacity limited diffusion coefficient is used. The
additional advantage of this multilayer approach is that there is a
variation (from single layer models) in bottom layer fluid velocity:
this dramatically reduces erosion, which is often overestimated in
simulations of this nature using single layer flows. The model is
used to simulate a standard dam break. In the dam break simulation,
as expected, the number of fluid layers utilised creates variation in
the resultant bed profile, with more layers offering a higher deviation
in fluid velocity . These results showed a marked variation in erosion
profiles from standard models. The overall the model provides new
insight into the problems presented at minimal computational cost.
Abstract: The one of most important objects in implementation
of damage analysis observations is manner of dam break wave
propagation. In this paper velocity and wave height due dam break in
with and without tailwater states for appointment hazardous lands
and flood radius are investigate. In order to modeling above
phenomenon finite volume method of Roe type for solving shallow
water equations is used. Results indicated that in the dry bed state
risk radius due to dam break is too high. While in the wet bed risk
radius has a less wide. Therefore in the first state constructions and
storage facilities are encountered with destruction risk. Further
velocity due to dam break in the second state is more comparing to
the first state. Hence erosion and scour the river bed in the dry bed is
too more compare to the wet bed.
Abstract: In this paper, a two-dimensional (2D) numerical
model for the tidal currents simulation in Persian Gulf is presented.
The model is based on the depth averaged equations of shallow water
which consider hydrostatic pressure distribution. The continuity
equation and two momentum equations including the effects of bed
friction, the Coriolis effects and wind stress have been solved. To
integrate the 2D equations, the Alternative Direction Implicit (ADI)
technique has been used. The base of equations discritization was
finite volume method applied on rectangular mesh. To evaluate the
model validation, a dam break case study including analytical
solution is selected and the comparison is done. After that, the
capability of the model in simulation of tidal current in a real field is
represented by modeling the current behavior in Persian Gulf. The
tidal fluctuations in Hormuz Strait have caused the tidal currents in
the area of study. Therefore, the water surface oscillations data at
Hengam Island on Hormoz Strait are used as the model input data.
The check point of the model is measured water surface elevations at
Assaluye port. The comparison between the results and the
acceptable agreement of them showed the model ability for modeling
marine hydrodynamic.
Abstract: Today, numerical simulation is a powerful tool to
solve various hydraulic engineering problems. The aim of this
research is numerical solutions of shallow water equations using
finite volume method for Simulations of dam break over wet and dry
bed. In order to solve Riemann problem, Roe-s approximate solver is
used. To evaluate numerical model, simulation was done in 1D and
2D states. In 1D state, two dam break test over dry bed (with and
without friction) were studied. The results showed that Structural
failure around the dam and damage to the downstream constructions
in bed without friction is more than friction bed. In 2D state, two
tests for wet and dry beds were done. Generally in wet bed case,
waves are propagated to canal sides but in dry bed it is not
significant. Therefore, damage to the storage facilities and
agricultural lands in wet bed case is more than in dry bed.
Abstract: To estimate the risks of dam failure phenomenon, it is necessary to understand this phenomenon and the involved governing factors. Overtopping and piping are the two main reasons of earthdam failures. In the piping context, the piping is determined as a phenomenon which is occurred between two phases, the water liquid and the solid soil. In this investigation, the onset of piping and its development, as well as the movement of water in soil, are numerically approached. In this regard, a one-dimensional numerical model based on the mass-conserving finite-volume method is developed and applied in order to simulate the piping phenomenon in a continuous circular tunnel of given initial length and radius, located between upstream and downstream. The simulation result includes the time-variations of radius along the tunnel until the radius value reaches its critical and the piping phenomenon converts to overtopping.