G. Tinoco-Guerrero*, F. J. Domínguez-Mota, J. G. Tinoco-Ruiz, J. S. Lucas-Martínez, N. S. Tinoco-Guerrero
Facultad de Ciencias Físico-Matemáticas |
Universidad Michoacana de San Nicolás de Hidalgo |
Morelia, Mexico |
keywords
hybrid method; finite difference; finite volume; shallow-water equations; irregular regions
In nature, there exist many types of flow that can be characterized as shallow-water flows. The main characteristic of these kinds of flows is that the vertical scales are much smaller than the horizontal ones. This happens in many regions all over the world, such as lakes, some rivers and, in some special cases, in parts of the oceans.
Since the flow in these cases is almost horizontal, a great number of simplifications can be done in the physical and mathematical formulation taking into account that the value of the pressure is essentially the hydrostatic one. It is important to remark that, even with these simplifications, the formulations are not two-dimensional, since the bottom friction must be taken into account in the boundary layers. In the literature, these effects are often not essential for the models, and it is only necessary to consider the two-dimensional depth-averaged form [1].
Due to the simplicity of the finite difference methods, some important advances have been done in the approximation to the solution of these equations; Po-Wei and Chia-Ming present in [2] a generalized finite difference method that can produce very good results, with the limitation that the average flow direction has to be known a priori in order to use the generalized finite difference method.
Young discusses a meshless method in [3,4] that can produce very good results in very irregular domains by using a local radial-basis-functions differential-quadrature approach; the results presented in his papers have very good quality and can be applied to real-life scenarios. One more time, even when this method produces very good results, even for inflow problems, the computational cost can be very high.
One of the most common problems that are presented in many works is the treatment of discontinuities and singularities; it is well known that one way to overcome these problems is using a conservative form of the equations. Nevertheless, Ulrik Skre Fjordhol [5], Bruno Gabutti [6], and Carlos Parés [7] have proposed accurate numerical schemes to approximate the solution of non-conservative Hyperbolic Equations, discussing upwind like and splitting schemes for solving numerically this version of the equations. Following the previous ideas, in this paper we present numerical schemes for the conservative form and non-conservative form of the shallow-water equations; the main aim of this work is to obtain an explicit method that can be applied to irregular domains.
In order to do that, let us first consider the problem of obtaining an approximation to the solution of the conservative form of the shallow-water equations
|
|
|
in a simply connected planar domain defined by a polygonal boundary, where is the water level, and are the velocity fields in the and directions respectively, is the water-body depth, is the Coriolis force and represents the external forces (see figure (1)).
![]() |
Figure 1: Definition of the bottom and free surface. |
In these equations, the change of variables
|
leads to the expressions
|
(1) |
|
(2) |
|
(3) |
here, the unknowns are the conservative quantities , and , that represent the mass and momentum of the physical problem.
After differentiation and some algebra, equations (1) - (3) can be rewritten as
|
(4) |
|
(5) |
|
(6) |
This is the non-conservative form or differential form of the shallow-water equations.
The proposed schemes arise from the integral form of equations (2) and (3), (5) and (6), and a finite difference approximation of equations (1) and (4).
On rectangular regions, the space region can be discretized by taking a grid formed by uniform cells defined as
|
as shown on figure (2).
![]() |
Figure 2: Rectangular meshed region . |
With this discretization of the region, the proposed scheme approximates the value of at the center of each cell by using a classical finite difference approximation. For this case, the partial derivatives can be approximated at a central point as
|
|
and
|
where the subscripts and superscripts represent the spatial position on the grid and the time level, respectively.
Taking into account these approximations, the first equation can be approximated as
|
From here, solving for , the approximation
|
(7) |
can be obtained.
Then, the values of and are approximated with a hybrid Finite Difference-Volume scheme on the edges of each cell. In order to do so the integral form of equations (2) and (3) must be taken into account.
The integral form of equation (2) evaluated over an arbitrary cell is
|
In this expression the partial derivatives can be replaced by their finite difference approximations, calculated on the edges of the cell,
|
|
|
and
|
to obtain
|
|
Since this integral must vanish for every cell , the expression
|
|
can be obtained. Now here, solving for ,
|
|
(8) |
Following a similar logic, the integral form of equation (3)
|
can be treated to obtain an approximation to ,
|
|
(9) |
Now, equations (7), (8) and (9) define a hybrid Finite Difference-Volume scheme for the conservative form of the shallow-water Equations, in rectangular regions.
In the case of the non-conservative form of shallow-water equations, a similar approach, as the one taken for the conservative form, can be chosen. In this occasion we could to take into account a finite difference scheme for equation (4),
|
and the integral form of equations (5) and (6),
|
|
Doing this, a scheme defined by
|
(10) |
|
|
|
(11) |
and
|
|
|
(12) |
can be obtained. This scheme can be applied to the non-conservative form of the shallow-water equations, in rectangular regions.
The hybrid schemes for non-rectangular regions are obtained analogously as the ones for rectangular regions, with the difference that, in this case, instead of replacing the spatial partial derivatives for their finite difference approximation, they are replaced for the generalized finite difference approximations. To learn further about generalized finite difference method see [8,9,10,11,12,13].
For the case of the integral form of (2) and (3),
|
|
the spatial partial derivatives can be replaced by their generalized finite difference approximation, while the temporal partial derivatives will be replaced with their finite difference approximation.
In order to address the definition of generalized finite difference, it is convenient to consider, for each case, the approximation to the first order operator
|
(13) |
where , , , , , , , and are given functions; the operator at some point can be approximated using values of in some neighbor points
|
Thus, a finite difference scheme at is a linear combination
|
where are suitable weights.
Since operator (13) is partially separable, it can be rewritten as
|
where
|
|
and
|
Each operator can be approximated separately, i.e.
|
(14) |
|
(15) |
and
|
(16) |
A finite difference scheme is consistent if the local truncation error satisfies
|
as [14,15].
In this case, this means that
|
(17) |
|
(18) |
and
|
(19) |
Expanding (17), (18) and (19) in Taylor series and regrouping terms, they can be written as
|
(20) |
|
(21) |
and
|
(22) |
where , for .
The consistency condition yields the undetermined systems
|
|
and
|
To find the values that satisfy these conditions requires to solve the systems
|
(23) |
|
(24) |
and
|
(25) |
it has to be noted that each system has equations and unknowns so, in general, in a suitable grid, there is a non trivial kernel. To select a solution, we use a subset of the normal equations of the corresponding least-squares problem; then, considering the last equations in each system
|
|
and
|
which can be solved through a reduced Cholesky factorization, the values of can be obtained.
After this, in order to obtain the values of , y , the first equations in (23), (24) and (25), can be used
|
(26) |
|
(27) |
|
(28) |
The resulting coefficients define the scheme (14 - 16).
This scheme can be used to approximate the operators
|
in order to get, for each equation separately, the approximations
|
(29) |
|
(30) |
Now, for equation (1)
|
a similar path must be taken in order to get a generalized finite difference approximation
|
(31) |
Equations (29), (30) and (31) define a hybrid Generalized Finite Difference-Volume scheme, for the conservative form, for non-rectangular regions.
In an analogous way, for the non-conservative form of the shallow-water equations, if the integral form of equations (5) and (6) is taken into account
|
|
and a generalized finite difference scheme is applied to equation (4),
|
the a scheme defined by
|
(32) |
|
(33) |
|
(34) |
can be found.
Now equations (32), (33) and (34) define a hybrid Generalized Finite Difference-Volume scheme, for the non-conservative form, for non-rectangular regions.
Now, an adequate selection of , the number of points, used by the scheme has to be done in order to represent different characteristics accurately. In this work, the selection of has been done as follows: for , and for and . The chosen stencils are shown in figure (3).![]() |
![]() |
(3) Different stencils used by the scheme. | |
![]() |
For the numerical tests, three different regions were selected: The unit square for the classical finite difference-finite volume hybrid scheme (denoted as QUAD), a widely used geometry denoted as DOME [16], that is a concave region limited by the lines , , and and an approximation to Zirahuen's Lake in Mexico (an endorheic basin), denoted as ZIRA. The corresponding normalized meshes with cells can be seen in figures (4-6).
![]() |
Figure 4: Regular mesh for the unitary square. |
![]() |
Figure 5: Mesh for DOME region. |
![]() |
Figure 6: Mesh for Zirahuen's Lake region. |
For each region 2 different tests were done, the first one using the proposed scheme for the conservative form of the Equations, and the second using the scheme for the non-conservative form.
The initial condition, for both tests, for was chosen to be a droplet with different center, according to the region, as follows:
The initial conditions for and where fixed as for all the regions. Also, reflective boundary conditions were chosen for and .
To produce stable calculations [17], the time discretization was chosen by considering
|
where and are the minimum values of and over the region and . With this, the time interval was divided into steps.
The total amounts of mass and momentum, over all the domain, are given at a time step as
|
and
|
In our case, they were approximated by means of a numerical quadrature.
The set of figures (7) and (8) show the results for the test using the conservative form of the equations in QUAD region, while the set of figures (9) and (10) show the results non-conservative form one for the same region. Following the same idea, the results for the region DOME are shown in the sets of figures (11) to (14), the same for ZIRA region in the sets (15) to (18). These sets of figures are shown as follows: the figures on the left show the behavior of the velocities while the figures on the right show the movement of the water. The figures begin showing the second time step (t = 0.001s) and continue with a plot every of the results every second after (t = 1s, t = 2s, t = 3s, t = 4s, t = 5s).
Figures (19) to (21) show the result of the computed total amount of conservative quantities over all the time steps. In all these figures the blue line represents the total amount of mass while the red-dotted line represents the total amount of momentum.
![]() |
![]() |
(7) Results for QUAD in the conservative case. | |
![]() |
![]() |
![]() |
(8) Results for QUAD in the conservative case. | |
![]() |
![]() |
![]() |
(9) Results for QUAD in the non-conservative case. | |
![]() |
![]() |
![]() |
(10) Results for QUAD in the non-conservative case. | |
![]() |
![]() |
![]() |
(11) Results for DOME in the conservative case. | |
![]() |
![]() |
![]() |
(12) Results for DOME in the conservative case. | |
![]() |
![]() |
![]() |
(13) Results for DOME in the non-conservative case. | |
![]() |
![]() |
![]() |
(14) Results for DOME in the non-conservative case. | |
![]() |
![]() |
![]() |
(15) Results for ZIRA in the conservative case. | |
![]() |
![]() |
![]() |
(16) Results for ZIRA in the conservative case. | |
![]() |
![]() |
![]() |
(17) Results for ZIRA in the non-conservative case. | |
![]() |
![]() |
![]() |
(18) Results for ZIRA in the non-conservative case. | |
![]() |
![]() |
![]() |
(19) and for QUAD. |
![]() |
![]() |
(20) and for DOME. |
![]() |
![]() |
(21) and for ZIRA. |
It must be pointed out that, in figures (19) to (21), the conservation of mass and momentum can be appreciated, since the losses and gains of these conservative quantities, over all the computed time, is small (around ).
Figures (19) to (21) show that, with the criteria used to select , neither spurious oscillations nor instabilities are observed in the tests. Even when these schemes can be applied to non-rectangular regions with ease, the quality of the grid is an important issue that must be taken into account; it is convenient to have quality grids to obtain better results, in all the tests of the present work quality grids were used.
It is important to remark that, in the tests, the boundary conditions where selected to be reflective and the tests show that, both proposed Generalized Finite Difference-Volume hybrid schemes, show a remarkable ability to produce stable results in the selected regions, even in the irregular cases; the numerical tests show that the conservation laws (mass and momentum) are fulfilled, correctly reflecting the expected behavior of a water body. As future work, these schemes need to be implemented in more realistic scenarios, which include inflows, outflows and non-slip boundary conditions.
We want to thank AULA CIMNE-Morelia and Finnish Government Scholarship Pool KM-17-10397 for the financial support for this work. We are grateful to the University of Helsinki for giving all the necessary work materials and the space to work at the Department of Physics in Helsinki, Finland.
[1] C. B. Vreugdenhil. (1994) "Numerical Methods for Shallow-Water Flow". Springer
[2] Po-Wei Li and Chia-Ming Fan. (2017) "Generalized Finite Difference Method for Two-Dimensional Shallow Water Equations", Volume 80. Engineering Analysis with Boundary Elements 58-71
[3] C. K. Chou and C. P. Sun and D. L. Young and J. Sladek and V. Sladek. (2015) "Extrapolated Local Radial Basis Function Collocation Method for Shallow Water Problems", Volume 50. Engineering Analysis with Boundary Elements 275-290
[4] C. P. Sun and D. L. Young and T. F. Chen and C. C. Hsian. (2013) "Application of Localized Meshless Methods to 2D Shallow Water Equation Problems", Volume 37. Engineering Analysis with Boundary Elements 1339-1350
[5] Ulrik Skre Fjordholm and Siddhartha Mishra. (2012) "Accurate numerical discretizations of non-conservative hyperbolic systems", Volume 46. ESAIM: Mathematical Modelling and Numerical Analysis 1 187-206
[6] Bruno Gabutti. (1983) "On Two Upwind Finite-Difference Schemes for Hyperbolic Equations in Non-Conservative Form", Volume 11. Computers & Fluids 3 207 - 230
[7] Carlos Parés. (2006) "Numerical Nethods for Nonconservative Hyperbolic Systems: A Theoretical Framework", Volume 44. SIAM Journal on Numerical Analysis 1 300-321
[8] Francisco Javier Domínguez-Mota and Pablo Michel Fernández-Valdez and Erika Ruiz-Díaz and Gerardo Tinoco-Guerrero and José Gerardo Tinoco-Ruiz. (2014) "An Heuristic Finite Difference Scheme on Irregular Plane Regions", Volume 8. Applied Mathematical Sciences 14 671-683
[9] Francisco Javier Domínguez-Mota and Sanzon Mendoza-Armenta and José Gerardo Tinoco-Ruiz. (2011) "Finite Difference Schemes Satisfying an Optimality Condition", Volume 17. IMACS. IMACS. MASCOT11 Proceedings 195
[10] Francisco Javier Domínguez-Mota and José Gerardo Tinoco-Ruiz and Gerardo Tinoco-Guerrero and Pablo Michel Fernández-Valdez and Erika Ruiz-Díaz. (2012) "A Modified Lax-Wendroff Scheme for Irregular 2D Space Regions", Volume 18. MASCOT12 Proceedings 101-110
[11] Gerardo Tinoco-Guerrero and Francisco Javier Domínguez-Mota and Ariana Gaona-Arias and Martha Leticia Ruiz-Zavala and José Gerardo Tinoco-Ruiz. (2018) "A Stability Analysis for a Generalized Finite-Difference Scheme Applied to the Pure Advection Equation", Volume 147. Mathematics and Computers in Simulation 293-300
[12] Gerardo Tinoco-Guerrero and Francisco Javier Domínguez-Mota and José Gerardo Tinoco-Ruiz and Erika Ruiz-Díaz. (2013) "An Implicit Modified Lax-Wendroff Scheme for Irregular 2D Space Regions", Volume 19. MASCOT13 Proceedings
[13] Gerardo Tinoco-Guerrero and Francisco Javier Domínguez-Mota and José Gerardo Tinoco-Ruiz and Erika Ruiz-Díaz. (2015) "Stability Aspects of a Modified Lax-Wendroff Scheme for Irregular 2D Regions", Volume 20. MASCOT15 Proceedings 161-170
[14] John C. Strikwerda. (2004) "Finite Difference Schemes and Partial Differential Equations". Society for Industrial and Applied Mathematics
[15] J. W. Thomas. (1998) "Numerical Partial Differential Equations: Finite Difference Methods". Springer
[16] Joe F. Thompson and Bharat K. Soni and Nigel P. Weatherill. (1999) "Handbook of Grid Generation". CRC Press 239 - 240
[17] Francisco Alcrudo and Pilar García-Navarro. (1993) "A high-resolution Godunov-type scheme in finite volumes for the 2D shallow-water equations", Volume 16. International Journal of Numerical Methods in Fluids 489 - 505
Published on 09/03/20
Submitted on 14/02/20
Volume 4, 2020
Licence: CC BY-NC-SA license
Are you one of the authors of this document?