Search

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.

Derived from simplifications of the Saint-Venant equations, the kinematic wave model has the ability to describe the behavior of surface runoff in watersheds. This paper aims to obtain the numerical simulation of the flow routing in a natural watershed, by using lattice Boltzmann method. In the computational model, the surface of the basin will be represented by a V-shaped segmented in two lateral planes and one main channel. The simulation considers the effective precipitation flowing on the watershed per unit of width at the exit of each of the planes that represent the surface of the basin. The water flowing from the planes enters the main channel in the form of lateral contribution. Hydrograms of two rain events are obtained, which present the volume drained in the outlet corresponding to the whole basin in each event. Two equilibrium distribution functions were developed by Chapmann-Enskog expansion at time scales and model D1Q3, one suitable for flow on the basin surface and another for the main channel, in order to obtain the variables of interest in each case. The numerical results obtained were compared with the KINEROS2 hydrological model.

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.

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.