Abstract

A two-dimensional numerical simulation of the water flow in a rectangular channel with submerged obstacles distributed alternately along its banks is presented. The governing equations of flow are the shallow water equations, which will be solved by the Boltzmann
lattice method (LBM) with multiple relaxation times (MRT). The non-slip bounce-back scheme was used on walls and obstacles, constant discharge at the inlet and fixed depth at the outlet of the channel. Due to the characteristics of the problem to be simulated, a large eddy simulation (LES) technique was incorporated into the computational code, which allows to obtain results
that are closer to the actual behavior of the flow. In addition, the stability of the simulation at all points of the mesh is evaluated for each step of time and, together with the property of the consistency of the LBM, the convergence of the solution is obtained. The simulation
provides the depth, velocities in the x and y directions, and the magnitude of water vorticity.

Abstract

This paper analyzes the ability of the Lattice Boltzmann method (LBM) with multiple relaxation times (MRT) in the simulation of flow in practical engineering problems. The case study covered refers to the first section of the initiation channel, which is part of the piracema channel, located in the Itaipu Hydroelectric Power Plant. The initiation channel has submerged obstacles distributed from one margin to the other, in order to reduce water velocity and allow the piracema cycle to occur. The governing equations of flow are the shallow water equations, which will be solved through the LBM-MRT. The non-slip bounce-back scheme was used on walls and obstacles, constant discharge at the inlet and fixed depth at the outlet of the channel. Due to the characteristics of the problem to be simulated, a large eddy simulation (LES) technique was incorporated into the computational code, which allows to obtain results that are closer to the actual behavior of the flow. In addition, the stability of the simulation at all points of the mesh is evaluated for each step of time and, together with the property of the consistency of the LBM, the convergence of the solution is obtained. The simulation provides the depth, velocities in the x and y directions, and the magnitude of water vorticity.

Abstract

An adaptive Finite Point Method (FPM) for solving shallow water problems is presented. The numerical methodology we propose, which is based on weighted‐least squares approximations on clouds of points, adopts an upwind‐biased discretization for dealing with the convective terms in the governing equations. The viscous and source terms are discretized in a pointwise manner and the semi‐discrete equations are integrated explicitly in time by means of a multi‐stage scheme. Moreover, with the aim of exploiting meshless capabilities, an adaptive h‐refinement technique is coupled to the described flow solver. The success of this approach in solving typical shallow water flows is illustrated by means of several numerical examples and special emphasis is placed on the adaptive technique performance. This has been assessed by carrying out a numerical simulation of the 26th December 2004 Indian Ocean tsunami with highly encouraging results. Overall, the adaptive FPM is presented as an accurate enough, cost‐effective tool for solving practical shallow water problems. Copyright © 2011 John Wiley & Sons, Ltd.

Abstract

In this paper, "finite point method" (FPM) is presented for modeling 2D shallow water flow problem. The method is based on the use of a weighted least-square approximation procedure, incorporating QR factorization and an iterative adjustment of local approximation parameters. The stabilization of the convective term in this present work is derived from the approximate Riemann solver proposed by Roe. The present method is shown to produce competitive accuracy in the comparisons with the analytical solutions and the well-known Galerkin characteristic-based split (CBS) algorithm.

Abstract

Due to the importance of the shallow-water equations in models of real-life phenomena, in recent years the study and model of problems that involve them have been the object of interest of many people. By reason of this, it is imperative to have efficient numerical methods to obtain an approximation of the solutions of the shallow-water equations.

Several authors have worked in approximations using the well-known finite volume and finite element methods, nevertheless, even when these methods compute good approximations to real-life behavior, the computational cost is usually high, which could be a limitation to the application of these methods.

This paper presents an explicit Generalized Finite Difference-Volume Hybrid approximation to the solution of the shallow-water equations, solved on irregular regions meshed with logically rectangular grids; the numerical results show the accuracy obtained with a low-cost implementation. The proposed scheme is a hybridization of a generalized finite difference scheme with the finite volume method.

Abstract

This work presents the numeric simulation of flow in a reservoir, using the finite volume method for solver the system of equations that model the two-dimensional flow in shallow water, neglecting the tangential tensions. The solution of the system of linearized equations was obtained by computational implementation of the Gauss-Seidel iterative method. To illustrate the application of the numerical scheme in hydraulic engineering, it was considered the study of flow in an underground reservoir with internal pillars, whose flow is generated by the opening of two outlet gates. The employed code allowed the temporal analysis of the volume and the flow in the reservoir, the graphical interpretation of the investigated physical phenomenon as well as, the calculation of the precision of the model. The results obtained show the expected physical interpretation, good agreement with literature data, good precision,stability and low computational cost.

Abstract

We present a stable finite element formulation for the shallow water equations using the finite increment calculus (FIC) procedure. This research is focused on the stability properties of the FIC technique and uses linear triangles for the spatial discretization with an equal order of interpolation for all the variables. The extension to higher order polynomial interpolation functions and different geometries is straightforward. The present FIC-FEM procedure is also able to introduce artificial viscosity for an adequate shock capturing. An implicit time integration has been used. Special attention has been payed to the dry domain in order to solve the moving boundaries with a fixed mesh Eulerian approach. Three academic examples are included in order to test the capabilities of the FIC-FEM procedure: the global stabilization, the shock capturing technique and the dry-wet interface. An experimental benchmark tests the overall accuracy of the present formulation.