NUMERICAL SOLUTION OF THE TWO LAYER SHALLOW WATER EQUATION USING FINITE VOLUME METHOD

2020 ◽  
Vol 19 (1) ◽  
pp. 57-71
Author(s):  
Arnasyita Yulianti Soelistya ◽  
Sumardi
2012 ◽  
Vol 2012 ◽  
pp. 1-14 ◽  
Author(s):  
Szu-Hsien Peng

The purpose of this study is to model the flow movement in an idealized dam-break configuration. One-dimensional and two-dimensional motion of a shallow flow over a rigid inclined bed is considered. The resulting shallow water equations are solved by finite volumes using the Roe and HLL schemes. At first, the one-dimensional model is considered in the development process. With conservative finite volume method, splitting is applied to manage the combination of hyperbolic term and source term of the shallow water equation and then to promote 1D to 2D. The simulations are validated by the comparison with flume experiments. Unsteady dam-break flow movement is found to be reasonably well captured by the model. The proposed concept could be further developed to the numerical calculation of non-Newtonian fluid or multilayers fluid flow.


2018 ◽  
Vol 7 (2) ◽  
pp. 94
Author(s):  
EKO MEIDIANTO N. R. ◽  
P. H. GUNAWAN ◽  
A. ATIQI ROHMAWATI

This one-dimensional simulation is performed to find the convergence of different fluxes on the water wave using shallow water equation. There are two cases where the topography is flat and not flat. The water level and grid of each simulation are made differently for each case, so that the water waves that occur can be analyzed. Many methods can be used to approximate the shallow water equation, one of the most used is the finite volume method. The finite volume method offers several numerical solutions for approximate shallow water equation, including Rusanov and HLLE. The derivation result of the numerical solution is used to approximate the shallow water equation. Differences in numerical and topographic solutions produce different waves. On flat topography, the rusanov flux has an average error of 0.06403 and HLLE flux with an average error of 0.06163. While the topography is not flat, the rusanov flux has a 1.63250 error and the HLLE flux has an error of 1.56960.


2020 ◽  
Vol 50 (3) ◽  
pp. 287-302
Author(s):  
Róbert ČUNDERLÍK ◽  
Matej MEDĽA ◽  
Karol MIKULA

The paper presents local quasigeoid modelling in Slovakia using the finite volume method (FVM). FVM is used to solve numerically the fixed gravimetric boundary value problem (FGBVP) on a 3D unstructured mesh created above the real Earth's surface. Terrestrial gravimetric measurements as input data represent the oblique derivative boundary conditions on the Earth's topography. To handle such oblique derivative problem, its tangential components are considered as surface advection terms regularized by a surface diffusion. The FVM numerical solution is fixed to the GOCE-based satellite-only geopotential model on the upper boundary at the altitude of 230 km. The horizontal resolution of the 3D computational domain is 0.002 × 0.002 deg and its discretization in the radial direction is changing with altitude. The created unstructured 3D mesh of finite volumes consists of 454,577,577 unknowns. The FVM numerical solution of FGBVP on such a detailed mesh leads to large-scale parallel computations requiring 245 GB of internal memory. It results in the disturbing potential obtained in the whole 3D computational domain. Its values on the discretized Earth's surface are transformed into the local quasigeoid model that is tested at 404 GNSS/levelling benchmarks. The standard deviation of residuals is 2.8 cm and decreases to 2.6 cm after removing 9 identified outliers. It indicates high accuracy of the obtained FVM-based local quasigeoid model in Slovakia.


2016 ◽  
Vol 2016 ◽  
pp. 1-7 ◽  
Author(s):  
Sudi Mungkasi

This paper presents a numerical entropy production (NEP) scheme for two-dimensional shallow water equations on unstructured triangular grids. We implement NEP as the error indicator for adaptive mesh refinement or coarsening in solving the shallow water equations using a finite volume method. Numerical simulations show that NEP is successful to be a refinement/coarsening indicator in the adaptive mesh finite volume method, as the method refines the mesh or grids around nonsmooth regions and coarsens them around smooth regions.


2014 ◽  
Vol 2014 ◽  
pp. 1-20 ◽  
Author(s):  
Sheng Bi ◽  
Jianzhong Zhou ◽  
Yi Liu ◽  
Lixiang Song

A second-order accurate, Godunov-type upwind finite volume method on dynamic refinement grids is developed in this paper for solving shallow-water equations. The advantage of this grid system is that no data structure is needed to store the neighbor information, since neighbors are directly specified by simple algebraic relationships. The key ingredient of the scheme is the use of the prebalanced shallow-water equations together with a simple but effective method to track the wet/dry fronts. In addition, a second-order spatial accuracy in space and time is achieved using a two-step unsplit MUSCL-Hancock method and a weighted surface-depth gradient method (WSDM) which considers the local Froude number is proposed for water depths reconstruction. The friction terms are solved by a semi-implicit scheme that can effectively prevent computational instability from small depths and does not invert the direction of velocity components. Several benchmark tests and a dam-break flooding simulation over real topography cases are used for model testing and validation. Results show that the proposed model is accurate and robust and has advantages when it is applied to simulate flow with local complex topographic features or flow conditions and thus has bright prospects of field-scale application.


2013 ◽  
Vol 392 ◽  
pp. 100-104 ◽  
Author(s):  
Fareed Ahmed ◽  
Faheem Ahmed ◽  
Yong Yang

In this paper we present a robust, high order method for numerical solution of multidimensional compressible inviscid flow equations. Our scheme is based on Nodal Discontinuous Galerkin Finite Element Method (NDG-FEM). This method utilizes the favorable features of Finite Volume Method (FVM) and Finite Element Method (FEM). In this method, space discretization is carried out by finite element discontinuous approximations. The resulting semi discrete differential equations were solved using explicit Runge-Kutta (ERK) method. In order to compute fluxes at element interfaces, we have used Roe Approximate scheme. In this article, we demonstrate the use of exponential filter to remove Gibbs oscillations near the shock waves. Numerical predictions for two dimensional compressible fluid flows are presented here. The solution was obtained with overall order of accuracy of 3. The numerical results obtained are compared with experimental and finite volume method results.


2006 ◽  
Vol 134 (9) ◽  
pp. 2515-2526 ◽  
Author(s):  
F. Xiao ◽  
X. D. Peng ◽  
X. S. Shen

Abstract This paper presents a novel finite-volume grid that uses not only the volume-integrated average (VIA) like the traditional finite-volume method, but also the surface-integrated average (SIA) as the model variables. The VIA and SIA are generically called “moments” in the context used here and are carried forward in time separately as the prognostic quantities. With the VIA defined in the control volume while the SIA is on the surface of the control volume, the discretization based on VIA and SIA leads to some new features in the numerical dispersions. A simple formulation using both VIA and SIA for shallow-water equations is presented. The numerical dispersion of the resulting grid, which is denoted as the “M grid,” is discussed with comparisons to the existing ones.


Sign in / Sign up

Export Citation Format

Share Document