Abstract: Modeling dam-break flows over non-flat beds requires
an accurate representation of the topography which is the main
source of uncertainty in the model. Therefore, developing robust
and accurate techniques for reconstructing topography in this class
of problems would reduce the uncertainty in the flow system. In
many hydraulic applications, experimental techniques have been
widely used to measure the bed topography. In practice, experimental
work in hydraulics may be very demanding in both time and cost.
Meanwhile, computational hydraulics have served as an alternative
for laboratory and field experiments. Unlike the forward problem,
the inverse problem is used to identify the bed parameters from the
given experimental data. In this case, the shallow water equations
used for modeling the hydraulics need to be rearranged in a way
that the model parameters can be evaluated from measured data.
However, this approach is not always possible and it suffers from
stability restrictions. In the present work, we propose an adaptive
optimal control technique to numerically identify the underlying bed
topography from a given set of free-surface observation data. In this
approach, a minimization function is defined to iteratively determine
the model parameters. The proposed technique can be interpreted
as a fractional-stage scheme. In the first stage, the forward problem
is solved to determine the measurable parameters from known data.
In the second stage, the adaptive control Ensemble Kalman Filter is
implemented to combine the optimality of observation data in order to
obtain the accurate estimation of the topography. The main features
of this method are on one hand, the ability to solve for different
complex geometries with no need for any rearrangements in the
original model to rewrite it in an explicit form. On the other hand, its
achievement of strong stability for simulations of flows in different
regimes containing shocks or discontinuities over any geometry.
Numerical results are presented for a dam-break flow problem over
non-flat bed using different solvers for the shallow water equations.
The robustness of the proposed method is investigated using different
numbers of loops, sensitivity parameters, initial samples and location
of observations. The obtained results demonstrate high reliability and
accuracy of the proposed techniques.
Abstract: Urban flooding resulting from a sudden release of
water due to dam-break or excessive rainfall is a serious threatening
environment hazard, which causes loss of human life and large
economic losses. Anticipating floods before they occur could
minimize human and economic losses through the implementation
of appropriate protection, provision, and rescue plans. This work
reports on the numerical modelling of flash flood propagation
in urban areas after an excessive rainfall event or dam-break.
A two-dimensional (2D) depth-averaged shallow water model is
used with a refined unstructured grid of triangles for representing
the urban area topography. The 2D shallow water equations are
solved using a second-order well-balanced discontinuous Galerkin
scheme. Theoretical test case and three flood events are described
to demonstrate the potential benefits of the scheme: (i) wetting and
drying in a parabolic basin (ii) flash flood over a physical model of
the urbanized Toce River valley in Italy; (iii) wave propagation on
the Reyran river valley in consequence of the Malpasset dam-break
in 1959 (France); and (iv) dam-break flood in October 1982 at the
town of Sumacarcel (Spain). The capability of the scheme is also
verified against alternative models. Computational results compare
well with recorded data and show that the scheme is at least as
efficient as comparable second-order finite volume schemes, with
notable efficiency speedup due to parallelization.
Abstract: In this paper, we show shallow water in a tin box as an analogous simulation tool for high-speed aerodynamics education and research. It is customary that we use a water tank to create shallow water flow. While a flow in a water tank is not necessarily uniform and is sometimes wavy, we can visualize a clear supercritical flow even when we move a body manually in stationary water in a simple shallow tin box. We can visualize a blunt shock wave around a moving circular cylinder together with a shock pattern around a diamond airfoil. Another interesting analogous experiment is a hydrodynamic shock tube with water and tea. We observe the contact surface clearly due to color difference of the two liquids those are invisible in the real gas dynamics experiment. We first revisit the similarities between high-speed aerodynamics and shallow water hydraulics. Several educational and research experiments are then introduced for engineering students. Shallow water experiments in a tin box simulate properly the high-speed flows.
Abstract: The study mapped selected wells in Inisa town, Osun state, in the guinea savanna region of southwest Nigeria, and determined the water quality considering certain elements. It also assessed the variation in the elevation of the water table surface to depth of the wells in the months of August and November. This is with a view to determine the level of contamination of the water with respect to land use and anthropogenic activities, and also to determine the variation that occurs in the quantity of well water in the rainy season and the start of the dry season. Results show a random pattern of the distribution of the mapped wells and shows that there is a shallow water table in the study area. The temporal changes in the elevation show that there are no significant variations in the depth of the water table surface over the period of study implying that there is a sufficient amount of water available to the town all year round. It also shows a high concentration of sodium in the water sample analyzed compared to other elements that were considered, which include iron, copper, calcium, and lead. This is attributed majorly to anthropogenic activities through the disposal of waste in landfill sites. There is a low concentration of lead which is a good indication of a reduced level of pollution.
Abstract: A coupled two-layer finite volume/finite element
method was proposed for solving dam-break flow problem
over deformable beds. The governing equations consist of the
well-balanced two-layer shallow water equations for the water flow
and a linear elastic model for the bed deformations. Deformations
in the topography can be caused by a brutal localized force or
simply by a class of sliding displacements on the bathymetry.
This deformation in the bed is a source of perturbations, on
the water surface generating water waves which propagate with
different amplitudes and frequencies. Coupling conditions at the
interface are also investigated in the current study and two mesh
procedure is proposed for the transfer of information through the
interface. In the present work a new procedure is implemented at
the soil-water interface using the finite element and two-layer finite
volume meshes with a conservative distribution of the forces at
their intersections. The finite element method employs quadratic
elements in an unstructured triangular mesh and the finite volume
method uses the Rusanove to reconstruct the numerical fluxes. The
numerical coupled method is highly efficient, accurate, well balanced,
and it can handle complex geometries as well as rapidly varying
flows. Numerical results are presented for several test examples of
dam-break flows over deformable beds. Mesh convergence study is
performed for both methods, the overall model provides new insight
into the problems at minimal computational cost.
Abstract: We present a new class of numerical techniques to
solve shallow water flows over dry areas including run-up. Many
recent investigations on wave run-up in coastal areas are based on
the well-known shallow water equations. Numerical simulations have
also performed to understand the effects of several factors on tsunami
wave impact and run-up in the presence of coastal areas. In all these
simulations the shallow water equations are solved in entire domain
including dry areas and special treatments are used for numerical
solution of singularities at these dry regions. In the present study we
propose a new method to deal with these difficulties by reformulating
the shallow water equations into a new system to be solved only in the
wetted domain. The system is obtained by a change in the coordinates
leading to a set of equations in a moving domain for which the
wet/dry interface is the reconstructed using the wave speed. To solve
the new system we present a finite volume method of Lax-Friedrich
type along with a modified method of characteristics. The method is
well-balanced and accurately resolves dam-break problems over dry
areas.
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: A numerical technique in a boundary-fitted curvilinear grid model is developed to simulate the extent of inland inundation along the coastal belts of Peninsular Malaysia and Southern Thailand due to 2004 Indian ocean tsunami. Tsunami propagation and run-up are also studied in this paper. The vertically integrated shallow water equations are solved by using the method of lines (MOL). For this purpose the boundary-fitted grids are generated along the coastal and island boundaries and the other open boundaries of the model domain. A transformation is used to the governing equations so that the transformed physical domain is converted into a rectangular one. The MOL technique is applied to the transformed shallow water equations and the boundary conditions so that the equations are converted into ordinary differential equations initial value problem. Finally the 4th order Runge-Kutta method is used to solve these ordinary differential equations. The moving boundary technique is applied instead of fixed sea side wall or fixed coastal boundary to ensure the movement of the coastal boundary. The extent of intrusion of water and associated tsunami propagation are simulated for the 2004 Indian Ocean tsunami along the west coast of Peninsular Malaysia and southern Thailand. The simulated results are compared with the results obtained from a finite difference model and the data available in the USGS website. All simulations show better approximation than earlier research and also show excellent agreement with the observed data.
Abstract: The aerial photogrammetry of shallow water bottoms has the potential to be an efficient high-resolution survey technique for shallow water topography, thanks to the advent of convenient UAV and automatic image processing techniques Structure-from-Motion (SfM) and Multi-View Stereo (MVS)). However, it suffers from the systematic overestimation of the bottom elevation, due to the light refraction at the air-water interface. In this study, we present an empirical method to correct for the effect of refraction after the usual SfM-MVS processing, using common software. The presented method utilizes the empirical relation between the measured true depth and the estimated apparent depth to generate an empirical correction factor. Furthermore, this correction factor was utilized to convert the apparent water depth into a refraction-corrected (real-scale) water depth. To examine its effectiveness, we applied the method to two river sites, and compared the RMS errors in the corrected bottom elevations with those obtained by three existing methods. The result shows that the presented method is more effective than the two existing methods: The method without applying correction factor and the method utilizes the refractive index of water (1.34) as correction factor. In comparison with the remaining existing method, which used the additive terms (offset) after calculating correction factor, the presented method performs well in Site 2 and worse in Site 1. However, we found this linear regression method to be unstable when the training data used for calibration are limited. It also suffers from a large negative bias in the correction factor when the apparent water depth estimated is affected by noise, according to our numerical experiment. Overall, the good accuracy of refraction correction method depends on various factors such as the locations, image acquisition, and GPS measurement conditions. The most effective method can be selected by using statistical selection (e.g. leave-one-out cross validation).
Abstract: The exhaustion of oil and the environmental pollution from the use of fossil fuel are increasing. Tidal current power (TCP) has been proposed as an alternative energy source because of its predictability and reliability. By applying a duct and single point mooring (SPM) system, a TCP device can amplify the generating power and keep its position properly. Because the generating power is proportional to cube of the current stream velocity, amplifying the current speed by applying a duct to a TCP system is an effective way to improve the efficiency of the power device. An SPM system can be applied at any water depth and is highly cost effective. Simple installation and maintenance procedures are also merits of an SPM system. In this study, we designed an SPM system for a duct-type TCP device for use in shallow water. Motions of the duct are investigated to obtain the response amplitude operator (RAO) as the magnitude of the transfer function. Parameters affecting the stability of the SPM system such as the fairlead departure angle, current velocity, and the number of clamp weights are analyzed and/or optimized. Wadam and OrcaFlex commercial software is used to design the mooring line.
Abstract: Tsunami and inundation modelling due to far field tsunami propagation in a limited area is a very challenging numerical task because it involves many aspects such as the formation of various types of waves and the irregularities of coastal boundaries. To compute the effect of far field tsunami and extent of inland inundation due to far field tsunami along the coastal belts of west coast of Malaysia and Southern Thailand, a formulated boundary condition and a moving boundary condition are simultaneously used. In this study, a boundary fitted curvilinear grid system is used in order to incorporate the coastal and island boundaries accurately as the boundaries of the model domain are curvilinear in nature and the bending is high. The tsunami response of the event 26 December 2004 along the west open boundary of the model domain is computed to simulate the effect of far field tsunami. Based on the data of the tsunami source at the west open boundary of the model domain, a boundary condition is formulated and applied to simulate the tsunami response along the coastal and island boundaries. During the simulation process, a moving boundary condition is initiated instead of fixed vertical seaside wall. The extent of inland inundation and tsunami propagation pattern are computed. Some comparisons are carried out to test the validation of the simultaneous use of the two boundary conditions. All simulations show excellent agreement with the data of observation.
Abstract: This paper describes the development of a boundary
fitted nested grid (BFNG) model to compute tsunami propagation of
2004 Indonesian tsunami in Southern Thailand coastal waters. We
develop a numerical model employing the shallow water nested
model and an orthogonal boundary fitted grid to investigate the
tsunami impact on the Southern Thailand due to the Indonesian
tsunami of 2004. Comparisons of water surface elevation obtained
from numerical simulations and field measurements are made.
Abstract: Oases are complex and fragile agro-ecosystems. They
have always existed in environments characterized by an arid climate,
scarcity of rainfall, high temperatures and high evaporation. These
palms have grown up despite the severity of the physical
characteristics thanks to the water's existence and irrigation practice.
The oases are generally spread along non-perennial rivers (wadis),
shallow water table or deep artesian groundwater. However, the
sustainability of oasis system is threatened by water scarcity and
declining of water table levels particularly in arid areas. Located in
the southern east area of Morocco, Tafilalet plain encompasses one of
the largest palm groves in the kingdom. In recent years, this area has
become increasingly threatened by water shortage and has seen a
sharp deterioration under the effect of several combined
anthropogenic and climatic factors. The Bayoud disease, successive
years of drought, Hassan Addakhil dam construction etc are all
factors that have affected both water and phoenicicole heritage of the
area. The objective of this study is to understand the interaction
between qualitative and quantitative degradation of groundwater
resources, and the palm grove dynamics, while reviewing the
assumption that groundwater resources contribute in a direct way to
the conservation of this oasis agroecosystem. A historical analysis
tracing both the oasis dynamics and the groundwater evolution has
been established. Data were collected from satellite images, surveys
with different actors (farmers, Regional Office for Agricultural
Development, Basin agency...). They were complemented by a
synthesis of numerous technical reports in the area. The results
showed that within 40 years, the thickness of the groundwater table
has dropped in 50 %. Along with this, there has been a downsizing of
date palm by 50 %. Areas with higher groundwater level were the
least affected by the downsizing. So we can say that the shallow
groundwater contribute significantly and directly to the water supply
of date palm through its root system, and largely ensures the oasis
ecosystem sustainability.
Abstract: This paper focuses on the development of a 2-D boundary fitted and nested grid (BFNG) model to compute the tsunami propagation of Indonesian tsunami 2004 along the coastal region of Penang in Peninsular Malaysia.
In the presence of a curvilinear coastline, boundary fitted grids are suitable to represent the model boundaries accurately. On the other hand, when large gradient of velocity within a confined area is expected, the use of a nested grid system is appropriate to improve the numerical accuracy with the least grid numbers.
This paper constructs a shallow water nested and orthogonal boundary fitted grid model and presents computational results of the tsunami impact on the Penang coast due to the Indonesian tsunami of 2004. The results of the numerical simulations are compared with available data.
Abstract: The effect of small non-parallelism of the base flow
on the stability of slightly curved mixing layers is analyzed in the
present paper. Assuming that the instability wavelength is much
smaller than the length scale of the variation of the base flow we
derive an amplitude evolution equation using the method of multiple
scales. The proposed asymptotic model provides connection between
parallel flow approximations and takes into account slow
longitudinal variation of the base flow.
Abstract: Method of multiple scales is used in the paper in order
to derive an amplitude evolution equation for the most unstable mode
from two-dimensional shallow water equations under the rigid-lid
assumption. It is assumed that shallow mixing layer is slightly curved
in the longitudinal direction and contains small particles. Dynamic
interaction between carrier fluid and particles is neglected. It is
shown that the evolution equation is the complex Ginzburg-Landau
equation. Explicit formulas for the computation of the coefficients of
the equation are obtained.
Abstract: This paper presents a new function expansion method for finding traveling wave solution of a non-linear equation and calls it the (G'/G)-expansion method. The shallow water wave equation is reduced to a non linear ordinary differential equation by using a simple transformation. As a result the traveling wave solutions of shallow water wave equation are expressed in three forms: hyperbolic solutions, trigonometric solutions and rational solutions.
Abstract: Linear stability analysis of wake-shear layers in twophase
shallow flows is performed in the present paper. Twodimensional
shallow water equations are used in the analysis. It is
assumed that the fluid contains uniformly distributed solid particles.
No dynamic interaction between the carrier fluid and particles is
expected in the initial moment. The stability calculations are
performed for different values of the particle loading parameter and
two other parameters which characterize the velocity ratio and the
velocity deficit. The results show that the particle loading parameter
has a stabilizing effect on the flow while the increase in the velocity
ratio or in the velocity deficit destabilizes the flow.
Abstract: The spectral action balance equation is an equation that
used to simulate short-crested wind-generated waves in shallow water
areas such as coastal regions and inland waters. This equation consists
of two spatial dimensions, wave direction, and wave frequency which
can be solved by finite difference method. When this equation with
dominating convection term are discretized using central differences,
stability problems occur when the grid spacing is chosen too coarse.
In this paper, we introduce the splitting upwind schemes for avoiding
stability problems and prove that it is consistent to the upwind scheme
with same accuracy. The splitting upwind schemes was adopted
to split the wave spectral action balance equation into four onedimensional
problems, which for each small problem obtains the
independently tridiagonal linear systems. For each smaller system
can be solved by direct or iterative methods at the same time which
is very fast when performed by a multi-processor computer.