Accessibility / Report Error

Comparison of numerical schemes of river flood routing with an inertial approximation of the Saint Venant equations

Comparação de esquemas numéricos de propagação de vazão com aproximação inercial das equações de Saint Venant

ABSTRACT

The one-dimensional flow routing inertial model, formulated as an explicit solution, has advantages over other explicit models used in hydrological models that simplify the Saint-Venant equations. The main advantage is a simple formulation with good results. However, the inertial model is restricted to a small time step to avoid numerical instability. This paper proposes six numerical schemes that modify the one-dimensional inertial model in order to increase the numerical stability of the solution. The proposed numerical schemes were compared to the original scheme in four situations of river’s slope (normal, low, high and very high) and in two situations where the river is subject to downstream effects (dam backwater and tides). The results are discussed in terms of stability, peak flow, processing time, volume conservation error and RMSE (Root Mean Square Error). In general, the schemes showed improvement relative to each type of application. In particular, the numerical scheme here called Prog Q(k+1)xQ(k+1) stood out presenting advantages with greater numerical stability in relation to the original scheme. However, this scheme was not successful in the tide simulation situation. In addition, it was observed that the inclusion of the hydraulic radius calculation without simplification in the numerical schemes improved the results without increasing the computational time.

Keywords:
Inertial model; Numerical stability; Computational time; HEC-RAS

RESUMO

O modelo inercial unidimensional de propagação de vazão, formulado através de uma solução explícita, possui resultados comparáveis aos modelos que consideram as equações completas de Saint-Venant. A principal vantagem é a simplicidade na formulação com bons resultados. No entanto, o modelo inercial está restrito à um passo de tempo pequeno para evitar instabilidade numérica. Este trabalho propõe seis esquemas numéricos que modificam o modelo inercial unidimensional buscando aumentar a estabilidade numérica da solução. O desempenho dos esquemas numéricos propostos foi comparado em quatro situações de declividade do rio (normal, baixa, alta e muito alta) e em duas situações onde o rio está sujeito ao efeito de jusante (remanso e maré). O modelo hidrodinâmico HEC-RAS foi utilizado para validação das soluções. Os resultados são discutidos em termos de estabilidade, vazão de pico, tempo de processamento, erro de conservação de volume e raiz do erro médio quadrático da vazão (Root Mean Square Error – RMSE). De forma geral os esquemas mostraram melhora relativa a cada tipo de aplicação. Em particular, o esquema numérico aqui denominado Prog Q(k+1)xQ(k+1) se destacou apresentando vantagens com maior estabilidade numérica em relação ao esquema original. Porém este esquema não foi bem sucedido na situação de simulação com maré. Além disso, observou-se que a inclusão do cálculo do raio hidráulico sem simplificação nos esquemas numéricos melhorou os resultados sem aumentar o tempo computacional.

Palavras-chave:
Modelo inercial; Estabilidade numérica; Tempo computacional; HEC-RAS

INTRODUCTION

River’s flow and the propagation of flood waves can be represented by hydrodynamic models that solve numerically the Saint-Venant equations. These are formed by the continuity equation (Equation 1), which represents the mass conservation, and by the dynamic equation (Equation 2), which represents the forces acting on the flow (CUNGE; HOLLY; VERWEY, 1980CUNGE, J. A.; HOLLY, F. M.; VERWEY, A. Practical aspects of computational river hydraulics. Boston: Pitman, 1980. 420 p.; TUCCI, 2005TUCCI, C. E. M. Modelos hidrológicos. 2. ed. Porto Alegre: Editora da UFRGS, 2005.). The two equations are presented considering a one-dimensional flow.

h t + Q B . x = 0 (1)
Q t 1 + ( Q 2 A ) x 2 = g A h x 3 + g A S 0 4 g A S f 5 (2)

where, Q is the flow, h is the water depth, t is the time, B is the cross-section width, x is the longitudinal distance, g is the gravitational acceleration, A is the cross-section area, S0 is the bed slope, and Sf is the friction slope.

The left side of the dynamic equation (Equation 2) refers to the local inertial (1) and advective inertial terms (2). The right side represents the external and internal forces acting on the fluid: pressure (3), weight (4) and friction (5).

The Saint-Venant equations can be simplified by eliminating some terms of the dynamic equation, which result in benefits for the models such as simpler formulation, more accessible programming codes, and gain in computational efficiency. For example, Kinematic Wave, Muskingum and Muskingum-Cunge models (CUNGE, 1969CUNGE, J. A. On the subject of a flood propagation computation method (Muskingum method). Journal of Hydrology Research, v. 7, n. 2, p. 205-230, 1969. http://dx.doi.org/10.1080/00221686909500264.
http://dx.doi.org/10.1080/00221686909500...
; FREAD, 1993FREAD, D. L. Flow routing in handbook of hydrology. In: MAIDMENT, D. R. (Ed.). Handbook of hydrology. New York: McGraw-Hill Inc., 1993. p. 10.1-10.36.; CHANSON, 2004CHANSON, H. The hydraulics of open channel flow: an introduction. 2nd ed. Oxford: Butterworth-Heinemann, 2004. 630 p.) consider only the terms of weight (4) and friction (5). These terms have a greater contribution to the flow acceleration than other (NÉELZ; PENDER, 2009NÉELZ, S.; PENDER, G. (2009). Desktop review of 2D hydraulic modelling packages. Bristol: Defra/Environment Agency. Science Report SC080035.). The Muskingum-Cunge model is often used as a flow propagation module in rainfall-runoff models (TUCCI, 2005TUCCI, C. E. M. Modelos hidrológicos. 2. ed. Porto Alegre: Editora da UFRGS, 2005.; BEVEN, 2011BEVEN, K. J. Rainfall-runoff modelling: the primer. New Jersey: John Wiley & Sons, 2011.). However, the traditional linear Muskingum-Cunge method does not adequately represent the flow in some situations where the flood wave celerity is not constant, or in situations with the backwater effect (PONCE, 1989PONCE, V. M. Engineering hydrology, principles and practices. New Jersey: Prentice Hall, 1989. 640 p.).

Another simplification, considered slightly more complete than Kinematic Wave, Muskingum and Muskingum-Cunge models, is the result of using only the terms of pressure (3), weight (4) and friction (5), called non-inertial (YEN; TSAI, 2001YEN, B. C.; TSAI, C. W.-S. Noninertial wave vs. diffusion wave in flood routing. Journal of Hydrology, v. 244, n. 1-2, p. 97-104, 2001. http://dx.doi.org/10.1016/S0022-1694(00)00422-4.
http://dx.doi.org/10.1016/S0022-1694(00)...
) and diffusion wave models (YEN; TSAI, 2001YEN, B. C.; TSAI, C. W.-S. Noninertial wave vs. diffusion wave in flood routing. Journal of Hydrology, v. 244, n. 1-2, p. 97-104, 2001. http://dx.doi.org/10.1016/S0022-1694(00)00422-4.
http://dx.doi.org/10.1016/S0022-1694(00)...
; CHANSON, 2004CHANSON, H. The hydraulics of open channel flow: an introduction. 2nd ed. Oxford: Butterworth-Heinemann, 2004. 630 p.). The diffusion wave model has the advantage of representing more adequately the attenuation of a flood wave propagating downstream. However, it has limitations in comparison with the complete Saint-Venant equations, as in representation of downstream effects (MONTERO et al., 2013MONTERO, R. A.; SCHWANENBERG, D.; HATZ, M.; BRINKMANN, M. Simplified hydraulic modelling in model predictive control of flood mitigation measures along rivers. Journal of Applied Water Engineering and Research, v. 1, n. 1, p. 17-27, 2013. http://dx.doi.org/10.1080/23249676.2013.827897.
http://dx.doi.org/10.1080/23249676.2013....
).

The inertial model, or local inertial, disregards only the advective inertial term (2). This model, popularized by Bates, Horritt and Fewtrell (2010)BATES, P. D.; HORRITT, M. S.; FEWTRELL, T. J. A simple inertial formulation of the shallow water equations for efficient two-dimensional flood inundation modelling. Journal of Hydrology, v. 387, n. 1, p. 33-45, 2010. http://dx.doi.org/10.1016/j.jhydrol.2010.03.027.
http://dx.doi.org/10.1016/j.jhydrol.2010...
, was widely applied to represent two-dimensional flow in flood areas (NEAL et al., 2011NEAL, J. C.; SCHUMANN, G. J.-P.; FEWTRELL, T. J.; BUDIMIR, M.; BATES, P. D.; MASON, D. Evaluating a new LISFLOOD-FP formulation using data for the summer 2007 floods in Tewkesbury, UK. Journal of Flood Risk Management, v. 4, n. 2, p. 88-95, 2011. http://dx.doi.org/10.1111/j.1753-318X.2011.01093.x.
http://dx.doi.org/10.1111/j.1753-318X.20...
; ALMEIDA et al., 2012ALMEIDA, G. A. M.; BATES, P.; FREER, J. E.; SOUVIGNET, M. Improving the stability of a simple formulation of the shallow water equations for 2-D flood modeling. Water Resources Research, v. 48, n. 5, p. W05528, 2012. http://dx.doi.org/10.1029/2011WR011570.
http://dx.doi.org/10.1029/2011WR011570...
; ALMEIDA; BATES, 2013ALMEIDA, G. A. M.; BATES, P. Applicability of the local inertial approximation of the shallow water equations to flood modeling. Water Resources Research, v. 49, n. 8, p. 4833-4844, 2013. http://dx.doi.org/10.1002/wrcr.20366.
http://dx.doi.org/10.1002/wrcr.20366...
). Montero et al. (2013)MONTERO, R. A.; SCHWANENBERG, D.; HATZ, M.; BRINKMANN, M. Simplified hydraulic modelling in model predictive control of flood mitigation measures along rivers. Journal of Applied Water Engineering and Research, v. 1, n. 1, p. 17-27, 2013. http://dx.doi.org/10.1080/23249676.2013.827897.
http://dx.doi.org/10.1080/23249676.2013....
demonstrated that the inertial model has advantages over the diffusion wave model when compared to complete solutions of the Saint-Venant equation. Fan et al. (2014)FAN, F. M.; PONTES, P. R. M.; PAIVA, R. C. D.; COLLISCHONN, W. Avaliação de um método de propagação de cheias em rios com aproximação inercial das equações de Saint-Venant. Revista Brasileira de Recursos Hídricos, v. 19, n. 4, p. 137-147, 2014. http://dx.doi.org/10.21168/rbrh.v19n4.p137-147.
http://dx.doi.org/10.21168/rbrh.v19n4.p1...
tested the inertial formulation for representing one-dimensional flow in rivers. The authors showed their applicability in scenarios with high and low rivers’ slope and subjected to downstream effects, such as backwater and tide. The results were comparable to the complete Saint-Venant equations. However, it was also observed that the model presents numerical instability for values of the Courant number much smaller than 1 (the stability condition requires that ghΔtΔx must be smaller than 1). Furthermore, Monteiro et al. (2015)MONTEIRO, L. R.; FAN, F. M.; COLLISCHONN, W.; SCHETTINI, E. B. C. Simulação da onda superficial provocada pelo fechamento de comportas utilizando uma aproximação inercial da Equação de Saint-Venant/ Surface waves caused by gates closure simulation with an inertial approximation of Saint-Venant Equation. Revista Brasileira de Recursos Hídricos, v. 20, n. 4, p. 927-936, 2015. http://dx.doi.org/10.21168/rbrh.v20n4.p927-936.
http://dx.doi.org/10.21168/rbrh.v20n4.p9...
tested the inertial solution for the simulation of waves caused by the closure of floodgates, and also obtaining results comparable to the complete solutions of the Saint-Venant equations.

The main advantage of the inertial model is the easy implementation with simple code and good results in relation to the solution of the Saint-Venant equations. As a result, the model has been used in hydrological models and flood models by several authors (DOTTORI; TODINI, 2011DOTTORI, F.; TODINI, E. Developments of a flood inundation model based on the cellular automata approach: testing different methods to improve model performance. Physics and Chemistry of the Earth Parts A/B/C, v. 36, n. 7, p. 266-280, 2011. http://dx.doi.org/10.1016/j.pce.2011.02.004.
http://dx.doi.org/10.1016/j.pce.2011.02....
; NEAL; SCHUMANN; BATES, 2012NEAL, J.; SCHUMANN, G.; BATES, P. A subgrid channel model for simulating river hydraulics and floodplain inundation over large and data sparse areas. Water Resources Research, v. 48, n. 11, 2012. http://dx.doi.org/10.1029/2012WR012514.
http://dx.doi.org/10.1029/2012WR012514...
; ALMEIDA; BATES, 2013ALMEIDA, G. A. M.; BATES, P. Applicability of the local inertial approximation of the shallow water equations to flood modeling. Water Resources Research, v. 49, n. 8, p. 4833-4844, 2013. http://dx.doi.org/10.1002/wrcr.20366.
http://dx.doi.org/10.1002/wrcr.20366...
; COULTHARD et al., 2013COULTHARD, T. J.; NEAL, J. C.; BATES, P. D.; RAMIREZ, J.; ALMEIDA, G. A. M.; HANCOCK, G. R. Integrating the LISFLOOD‐FP 2D hydrodynamic model with the CAESAR model: implications for modelling landscape evolution. Earth Surface Processes and Landforms, v. 38, n. 15, p. 1897-1906, 2013. http://dx.doi.org/10.1002/esp.3478.
http://dx.doi.org/10.1002/esp.3478...
; YAMAZAKI; ALMEIDA; BATES, 2013YAMAZAKI, D.; ALMEIDA, G. A. M.; BATES, P. D. Improving computational efficiency in global river models by implementing the local inertial flow equation and a vector-based river network map. Water Resources Research, v. 49, n. 11, p. 7221-7235, 2013.; SAMPSON et al., 2015SAMPSON, C. C.; SMITH, A. M.; BATES, P. D.; NEAL, J. C.; ALFIERI, L.; FREER, J. E. A high‐resolution global flood hazard model. Water Resources Research, v. 51, n. 9, p. 7358-7381, 2015. PMid:27594719. http://dx.doi.org/10.1002/2015WR016954.
http://dx.doi.org/10.1002/2015WR016954...
). In Brazil, the authors Pontes et al. (2015PONTES, P. R. M.; COLLISCHONN, W.; FAN, F.M.; PAIVA, R.C.D.; BUARQUE, D.C. Modelagem hidrológica e hidráulica de grande escala com propagação inercial de vazões. Revista Brasileira de Recursos Hídricos, vol. 20, n. 4, p. 888-904, 2015. http://dx.doi.org/10.21168/rbrh.v20n4.p888-904.
http://dx.doi.org/10.21168/rbrh.v20n4.p8...
, 2017PONTES, P. R. M.; FAN, F. M.; FLEISCHMANN, A. S.; PAIVA, R. C. D.; BUARQUE, D. C.; SIQUEIRA, V. A.; JARDIM, P. F.; SORRIBAS, M. V.; COLLISCHONN, W. MGB-IPH model for hydrological and hydraulic simulation of large floodplain river systems coupled with open source GIS. Environmental Modelling & Software, v. 94, p. 1-20, 2017. http://dx.doi.org/10.1016/j.envsoft.2017.03.029.
http://dx.doi.org/10.1016/j.envsoft.2017...
) presented coupling of the one-dimensional inertial model in a large-scale hydrological model (MGB-IPH). The authors pointed out that with this model it is possible to perform site simulations, as in floodplains, more adequately than when using other simplified methods. They showed that the scheme has robustness in coupling with hydrological models and, in addition, the parallel processing is easy.

As commented, several studies point to the inertial model as an alternative to the complete solution of the Saint-Venant equations and to other simplified models representing the one-dimensional flow. However, alternatives to improve the model's numerical stability may still be explored. Therefore, this paper proposes and tests some changes in the inertial formulation in order to improve the numerical stability. The numerical scheme of the one-dimensional inertial model, as applied by Fan et al. (2014)FAN, F. M.; PONTES, P. R. M.; PAIVA, R. C. D.; COLLISCHONN, W. Avaliação de um método de propagação de cheias em rios com aproximação inercial das equações de Saint-Venant. Revista Brasileira de Recursos Hídricos, v. 19, n. 4, p. 137-147, 2014. http://dx.doi.org/10.21168/rbrh.v19n4.p137-147.
http://dx.doi.org/10.21168/rbrh.v19n4.p1...
, was altered in relation to the approximation of the numerical derivative, in relation to the numerical method and also in relation to the dynamic equation’s friction term (term 5 in Equation 2). Six new schemes were proposed and tested considering the hydraulic radius without simplification (RhB/(2.h+B)) and six numerical schemes with hydraulic radius simplification (Rh).

The performance of the schemes was compared with original scheme, presented by Bates, Horritt and Fewtrell (2010)BATES, P. D.; HORRITT, M. S.; FEWTRELL, T. J. A simple inertial formulation of the shallow water equations for efficient two-dimensional flood inundation modelling. Journal of Hydrology, v. 387, n. 1, p. 33-45, 2010. http://dx.doi.org/10.1016/j.jhydrol.2010.03.027.
http://dx.doi.org/10.1016/j.jhydrol.2010...
, and with the complete solution of the Saint Venant equations (HEC-RAS hydrodynamic model - USACE, 2016US ARMY CORPS OF ENGINEERS – USACE. HEC-RAS River Analysis System Hydraulic Reference Manual. Version 5.0. Davis: Hydrologic Engineering Center, 2016.), i.e. the HEC-RAS model was considered the “true solution” in tests. This model solves the equations by an implicit finite difference scheme described in USACE (2016)US ARMY CORPS OF ENGINEERS – USACE. HEC-RAS River Analysis System Hydraulic Reference Manual. Version 5.0. Davis: Hydrologic Engineering Center, 2016.. For the sake of isonomy with the knowledge already generated and published on the one-dimensional inertial model application, scenarios identical to those used by Fan et al. (2014)FAN, F. M.; PONTES, P. R. M.; PAIVA, R. C. D.; COLLISCHONN, W. Avaliação de um método de propagação de cheias em rios com aproximação inercial das equações de Saint-Venant. Revista Brasileira de Recursos Hídricos, v. 19, n. 4, p. 137-147, 2014. http://dx.doi.org/10.21168/rbrh.v19n4.p137-147.
http://dx.doi.org/10.21168/rbrh.v19n4.p1...
were considered. The scenarios were: four situations of river’s slope (normal, low, high and very high) and two situations where the river is subject to downstream effect (backwater and tide).

ONE-DIMENSIONAL INERTIAL MODEL

The dynamic equation (Equation 2) in the one-dimensional inertial model is rewritten (Equation 3) considering a rectangular river cross-section, representing the pressure and weight forces in terms of water surface level (y) and estimating the friction force using the Manning equation,

Q t + g B h y x + g | Q | Q n 2 B h R 4 3 = 0 (3)

where y is the water surface level, n is the Manning roughness coefficient, R is the hydraulic radius.

The hydraulic radius is approximated by the depth (Rh), considering that the rivers cross sections have a width much larger than the depth. The derivatives of the dynamic equation are approximated by a numerical scheme of progressive finite differences in space and time, which results in Equation 4, as presented in Fan et al. (2014)FAN, F. M.; PONTES, P. R. M.; PAIVA, R. C. D.; COLLISCHONN, W. Avaliação de um método de propagação de cheias em rios com aproximação inercial das equações de Saint-Venant. Revista Brasileira de Recursos Hídricos, v. 19, n. 4, p. 137-147, 2014. http://dx.doi.org/10.21168/rbrh.v19n4.p137-147.
http://dx.doi.org/10.21168/rbrh.v19n4.p1...
.

Q i + 1 2 k + 1 Q i + 1 2 k Δ t + g B h i + 1 2 k y i + 1 k y i k Δ x + g | Q i + 1 2 k | Q i + 1 2 k + 1 n 2 B ( h i + 1 2 k ) 7 3 = 0 (4)

Where i and k are the indices in space and time, respectively. Δx is the length of the channel section (m).

The hi+12k term is the depth in the contour of the channel stretch (Δx) and can be estimated by the difference of the maximum value of the water level (y) and the maximum value of the bottom level (z) in the center of the channel stretch, according to Equation 5.

h i + 1 2 k = max [ y i k ; y i + 1 k ] max [ z i ; z i + 1 ] (5)

From the explicit solution of Equation 4, the flow is estimated at each time step according to Equation 6. Bates, Horritt and Fewtrell (2010)BATES, P. D.; HORRITT, M. S.; FEWTRELL, T. J. A simple inertial formulation of the shallow water equations for efficient two-dimensional flood inundation modelling. Journal of Hydrology, v. 387, n. 1, p. 33-45, 2010. http://dx.doi.org/10.1016/j.jhydrol.2010.03.027.
http://dx.doi.org/10.1016/j.jhydrol.2010...
considered the multiplication of the flow in time step k by the flow in time step k+1 in friction term (|Qi+12k|Qi+12k+1) to increase the stability of the numerical scheme.

Q i + 1 2 k + 1 = ( Q i + 1 2 k g B Δ t h i + 1 2 k y i + 1 k y i k Δ x ) ( 1 + g Δ t | Q i + 1 2 k | . n 2 B . ( h i + 1 2 k ) 7 3 ) (6)

Based on the flow values (Qi+12k+1), the continuity equation (Equation 1), discretized by a spatial and time progressive scheme (Equation 7), is applied to estimate the depth at all reaches i at the end of the time interval (k+1). From this value, the water level at the end of the time interval is computed by Equation 8.

h i k + 1 = h i k Δ t B Δ x ( Q i + 1 2 k + 1 Q i 1 2 k + 1 ) (7)
y i k + 1 = z i + h i k + 1 (8)

Equation 6 requires the initial values of the depth (h0) and the water level (y0) in all reaches. For this, the same initial flow is assumed in all reaches and the depth is calculated by considering the continuous and uniform flow, that is, the Manning equation with simplification of the hydraulic radius (Rh). In addition, boundary conditions are required in the first and last channel stretch. The typical upstream condition is a flow hydrograph and the downstream condition can be a rating curve or a normal depth.

The numerical scheme is restricted to the size of time step (Δt) and channel stretch (Δx) to avoid numerical instability. Thus, the choice of Δt and Δx must respect the Courant-Friedrichs-Levy condition (Courant number – C), given by Equation 9.

C = g h Δ t Δ x w h e r e C 1 (9)

The depth (h) used in this calculation corresponds to the largest depth found in the river reach. Although this condition requires a value less than 1, Bates, Horritt and Fewtrell (2010)BATES, P. D.; HORRITT, M. S.; FEWTRELL, T. J. A simple inertial formulation of the shallow water equations for efficient two-dimensional flood inundation modelling. Journal of Hydrology, v. 387, n. 1, p. 33-45, 2010. http://dx.doi.org/10.1016/j.jhydrol.2010.03.027.
http://dx.doi.org/10.1016/j.jhydrol.2010...
observed that the value must be less than 0.7 to ensure the numerical scheme stability in two-dimensional simulations.

After the initial development of the model, Almeida et al. (2012)ALMEIDA, G. A. M.; BATES, P.; FREER, J. E.; SOUVIGNET, M. Improving the stability of a simple formulation of the shallow water equations for 2-D flood modeling. Water Resources Research, v. 48, n. 5, p. W05528, 2012. http://dx.doi.org/10.1029/2011WR011570.
http://dx.doi.org/10.1029/2011WR011570...
proposed a modification in the numerical scheme to improve stability in low friction situations, such as in urban areas. In the new scheme, a weighting (θ) was introduced in the flow derivative in relation to the time. Thus, the flow value at time k and reach i+12 is weighted considering the flow at time k and in the reach i+12, i12, and i+32, according to Equation 10.

Q t = Q i + 1 2 k + 1 ( θ Q i + 1 2 k + ( 1 θ ) 2 ( Q i 1 2 k + Q i + 3 2 k ) ) Δ t (10)

The weighting (θ) varies from 0 to 1. Using Equation 10 to approximate the numerical derivative, the flow in each time step in the new scheme is estimated according to Equation 11.

Q i + 1 2 k + 1 = ( θ Q i + 1 2 k + ( 1 θ ) 2 ( Q i 1 2 k + Q i + 3 2 k ) ) g B Δ t h i + 1 2 k y i + 1 k y i k Δ x ( 1 + g Δ t | Q i + 1 2 k | . n 2 B . ( h i + 1 2 k ) 7 3 ) (11)

According to Almeida et al. (2012)ALMEIDA, G. A. M.; BATES, P.; FREER, J. E.; SOUVIGNET, M. Improving the stability of a simple formulation of the shallow water equations for 2-D flood modeling. Water Resources Research, v. 48, n. 5, p. W05528, 2012. http://dx.doi.org/10.1029/2011WR011570.
http://dx.doi.org/10.1029/2011WR011570...
, the inclusion of the weighting allowed more stability in simulations performed by the authors.

MODIFICATIONS PROPOSED TO THE NUMERICAL SCHEME OF THE INERTIAL MODEL

In the inertial model formulation proposed by Bates, Horritt and Fewtrell (2010)BATES, P. D.; HORRITT, M. S.; FEWTRELL, T. J. A simple inertial formulation of the shallow water equations for efficient two-dimensional flood inundation modelling. Journal of Hydrology, v. 387, n. 1, p. 33-45, 2010. http://dx.doi.org/10.1016/j.jhydrol.2010.03.027.
http://dx.doi.org/10.1016/j.jhydrol.2010...
, the flow at time step k+1 (Equation 6) is solved by the Euler method, i.e. a first order method where the first derivative provides a direct estimate of the slope for entire interval (Δt). In this method, the global truncation error is proportional to the time step size (CHAPRA; CANALE, 2012CHAPRA, S. C.; CANALE, R. P. Numerical methods for engineers. New York: McGraw-Hill, 2012.). The second-order Runge-Kutta method (RK2), equivalent to the Heun method without iteration, and the fourth-order Runge-Kutta method (RK4) are more accurate because the errors are of the order of Δt2 and Δt4, respectively (CHAPRA; CANALE, 2012CHAPRA, S. C.; CANALE, R. P. Numerical methods for engineers. New York: McGraw-Hill, 2012.). Thus, the reduction of the step size decreases the error at a faster rate than in the Euler method (CHAPRA; CANALE, 2012CHAPRA, S. C.; CANALE, R. P. Numerical methods for engineers. New York: McGraw-Hill, 2012.). Another way to improve the accuracy is to use a centered derivative (Cent) instead of the progressive derivative (Prog). This derivative allows a second-order approximation, with an error of the order of Δt2.

In this paper, the modifications of the numerical schemes were done in the derivative approximation (Prog and Cent), in the numerical method (Euler, RK2 ou RK4) and also in the friction term (Equation 2), specifically in the time index (k) of the flow (|Qi+12k|Qi+12k+1).

The properties of the proposed numerical schemes, such as convergence, consistency and stability, were not evaluated. These analyzes indicate how close the numerical scheme is to the analytical solution, and if the errors are not amplified (POPESCU, 2014POPESCU, I. Computational hydraulics. London: IWA Publishing, 2014.). Thus, the numerical schemes were only tested and evaluated in different numerical applications. In each test, 16 simulations were run, each one with the methods listed in Table 1 and described in detail below (“R” means that the hydraulic radius was not simplified).

Table 1
Summary of proposed methods and characteristics.

Scheme Q(k+1) x Q(k-1):

The first modification made in the original numerical scheme was to approximate the flow by a time-centered derivative, which considers the flow in the time steps k+1 and k-1. Thus, the proposed scheme also considered the flow in these time steps in the friction term (|Qi+12k1|Qi+12k+1). Rearranging, the flow Qi+12k+1 is estimated directly by Equation 12. This scheme will be identified by the abbreviation: Cent Q(k+1) x Q(k-1).

Q i + 1 2 k + 1 = ( Q i + 1 2 k 1 2 Δ t g B h i + 1 2 k y i + 1 k y i k Δ x ) ( 1 + g 2 Δ t | Q i + 1 2 k 1 | . n 2 B . ( h i + 1 2 k ) 7 3 ) (12)

Scheme Q(k+1) x Q(k+1):

In this numerical scheme, the friction term was calculated by the flow at time step k+1, (|Qi+12k+1|Qi+12k+1). Equation 13 shows the scheme for a progressive derivative, identified as Prog Q(k+1) x Q(k+1), and Equation 14 for a centered derivative, identified as Cent Q(k+1) x Q(k+1).

g Δ t n 2 B . ( h i + 1 2 k ) 7 3 ( Q i + 1 2 k + 1 ) 2 + Q i + 1 2 k + 1 + Δ t g B h i + 1 2 k y i + 1 k y i k Δ x Q i + 1 2 k = 0 (13)
g 2 Δ t n 2 B . ( h i + 1 2 k ) 7 3 ( Q i + 1 2 k + 1 ) 2 + Q i + 1 2 k + 1 + 2 Δ t g B h i + 1 2 k y i + 1 k y i k Δ x Q i + 1 2 k 1 = 0 (14)

Note that both equations are second-degree polynomial equations: ax2+bx+c=0. Therefore, the flow can be calculated directly by the Bhaskara’s formula, without the need for an iterative method for resolution. The desired estimate is always the largest root given by: x1=b+ b24ac 2a. The reasoning is based on the observation that the root x2 (b b24ac 2a) is always negative, as shown below.

It is known that the roots sum is given by x1+x2=ba. Since b=1 and a>0 (g, n, B, h Δt are always positive), we have x1+x2<0. Moreover, since b24ac2a is always positive, so x1> x2. Hence x2<0, otherwise the sum of the two roots would be positive. Thus, the root x2 does not have physical importance as an always negative flow. On the other hand, the root x1 can be positive or negative depending on the term c. As x1.x2=ca and x2<0, then the root x1 is positive if c<0, and negative if c>0. In fact, if the flow Qi+12k (or Qi+12k1) is negative, as occurs when the river is subject to the downstream effect, the term c will be positive and the flow at time step k+1, Qi+12k+1, is negative.

The scheme is restricted to a new stability condition, in addition to the Courant-Friedrichs-Levy condition, where the discriminant of the Bhaskara’s formula must be positive (b24ac>0) so that the root is not complex. Thus, the condition is given by ac<14. When the flow is positive (x1>0) it was shown that a>0 and c<0, satisfying the condition. However, when the flow is negative (x1<0), the term c is also negative and this condition is mainly dependent on the time step size (Δt), the channel stretch size (Δx) and the level variation between reaches (yi+1kyik).

Scheme Q(k) x Q(k):

Finally, the original numerical scheme was modified by considering the flow at time step k (|Qi+12k|Qi+12k). This scheme allows to calculate the flow directly (Euler) and by the second- and fourth-order Runge-Kutta methods. The flow in time step k+1, calculated by the Euler method, is given by Equation 15 and it is identified as Prog Q(k)xQ(k).

Q i + 1 2 k + 1 = g Δ t n 2 B . ( h i + 1 2 k ) 7 3 | Q i + 1 2 k | Q i + 1 2 k Δ t g B h i + 1 2 k y i + 1 k y i k Δ x + Q i + 1 2 k (15)

In the second-order Runge-Kutta method, the flow determined by the Euler method is an intermediate result, where the slope estimate is calculated at the initial point (predictor equation). In the next step, a new estimate is calculated with the predictor result, i.e., the slope estimate is made at the end point. Then, the mean of the two derivatives is calculated to obtain an improved estimate of the slope over the whole range (broker equation). Therefore, Equation 15 is rearranged to define the derivative function f(Qi+12k ) (Equation 16), the predictor step is given by Equation 17 (equivalent to Equation 15) and the corrector step by Equation 18. This numerical scheme is identified by the abbreviation Prog RK2 Q(k) x Q(k).

f ( Q i + 1 2 k ) = Q i + 1 2 k + 1 Q i + 1 2 k Δ t f ( Q i + 1 2 k ) = g n 2 B . ( h i + 1 2 k ) 7 3 | Q i + 1 2 k | Q i + 1 2 k g B h i + 1 2 k y i + 1 k y i k Δ x (16)
Q i + 1 2 k + 1 p r e d = ( Q i + 1 2 k ) + Δ t f ( Q i + 1 2 k ) (17)
Q i + 1 2 k + 1 = ( Q i + 1 2 k ) + Δ t f ( Q i + 1 2 k + 1 p r e d ) + f ( Q i + 1 2 k ) 2 (18)

Multiple slope estimates are developed in the fourth-order Runge-Kutta method to obtain an improved mean slope (Equations 19 to 22). Each K value represents a slope and the final equation is a weighted average (Equation 23). This scheme is identified by the abbreviation Prog RK4 Q(k)xQ(k). Note that K1 is the slope at the interval beginning and K2 is the slope at the interval final, equivalent to the second-order Runge-Kutta method.

K 1 = f ( Q i + 1 2 k ) (19)
K 2 = f ( Q i + 1 2 k + 1 2 K 1 Δ t ) (20)
K 3 = f ( Q i + 1 2 k + 1 2 K 2 Δ t ) (21)
K 4 = f ( Q i + 1 2 k + K 3 Δ t ) (22)
Q i + 1 2 k + 1 = Q i + 1 2 k + Δ t 6 ( K 1 + 2 K 2 + 2 K 3 + K 4 ) (23)

EXPERIMENTS

The described numerical schemes were evaluated in different applications to compare the performance in relation to the same schemes adapted by Fan et al. (2014)FAN, F. M.; PONTES, P. R. M.; PAIVA, R. C. D.; COLLISCHONN, W. Avaliação de um método de propagação de cheias em rios com aproximação inercial das equações de Saint-Venant. Revista Brasileira de Recursos Hídricos, v. 19, n. 4, p. 137-147, 2014. http://dx.doi.org/10.21168/rbrh.v19n4.p137-147.
http://dx.doi.org/10.21168/rbrh.v19n4.p1...
of the formulation proposed by Bates, Horritt and Fewtrell (2010)BATES, P. D.; HORRITT, M. S.; FEWTRELL, T. J. A simple inertial formulation of the shallow water equations for efficient two-dimensional flood inundation modelling. Journal of Hydrology, v. 387, n. 1, p. 33-45, 2010. http://dx.doi.org/10.1016/j.jhydrol.2010.03.027.
http://dx.doi.org/10.1016/j.jhydrol.2010...
and the formulation proposed by Almeida et al. (2012)ALMEIDA, G. A. M.; BATES, P.; FREER, J. E.; SOUVIGNET, M. Improving the stability of a simple formulation of the shallow water equations for 2-D flood modeling. Water Resources Research, v. 48, n. 5, p. W05528, 2012. http://dx.doi.org/10.1029/2011WR011570.
http://dx.doi.org/10.1029/2011WR011570...
. The schemes were implemented in MATLAB and the simulations were run on eight cores of four Quad-core 3.4 GHz Intel Xeon processors with 8 MB cache memory. The schemes were evaluated for stability (maximum Courant number for which the method remains stable), peak flow at the river’s last reach (Qpeak), volume conservation error (VolErr), and the flow's mean square error (RSME), calculated based on the results of simulations performed in the HEC-RAS hydrodynamic model. In addition, each proposed scheme was solved with and without the hydraulic radius simplification, that is, Rh and R=h.B/(2.h+B).

The applications were equivalent to the scenarios tested by Fan et al. (2014)FAN, F. M.; PONTES, P. R. M.; PAIVA, R. C. D.; COLLISCHONN, W. Avaliação de um método de propagação de cheias em rios com aproximação inercial das equações de Saint-Venant. Revista Brasileira de Recursos Hídricos, v. 19, n. 4, p. 137-147, 2014. http://dx.doi.org/10.21168/rbrh.v19n4.p137-147.
http://dx.doi.org/10.21168/rbrh.v19n4.p1...
. The first test considered a real reach of the São Francisco River, located between Três Marias UHE and the Pirapora city (MG), without lateral contributions. Tests 2, 3 and 4 considered the same reach but with low, high and very high slopes, respectively. In tests 5 and 6, downstream effects were considered: backwater and tide. The hydrograph input was the same as used by Fan et al. (2014)FAN, F. M.; PONTES, P. R. M.; PAIVA, R. C. D.; COLLISCHONN, W. Avaliação de um método de propagação de cheias em rios com aproximação inercial das equações de Saint-Venant. Revista Brasileira de Recursos Hídricos, v. 19, n. 4, p. 137-147, 2014. http://dx.doi.org/10.21168/rbrh.v19n4.p137-147.
http://dx.doi.org/10.21168/rbrh.v19n4.p1...
. The simulated period in all tests was 150 hours. Table 2 presents a summary of the data used in the six tests.

Table 2
Summary of the parameters used in tests.

RESULTS

Test 1

This test considered a real reach of the São Francisco River between the Três Marias UHE and the Pirapora city (MG) with approximate values of 135 km in length, 300 m in width (rectangular cross-section), river’s slope of 29.5 cm/km and Manning’s roughness coefficient of 0.030. The spatial discretization (Δx) was 2 km and time discretization (Δt) was varied between 1 and 4 minutes, assuming C values equal to 0.26, 0.51, 0.77 and 1.03 at the time of greater depth. A constant slope was considered as a downstream contour condition.

The six formulations of the numerical models described (Equations 12, 13, 14, 15, 18 and 23), the formulation proposed by Bates, Horritt and Fewtrell (2010)BATES, P. D.; HORRITT, M. S.; FEWTRELL, T. J. A simple inertial formulation of the shallow water equations for efficient two-dimensional flood inundation modelling. Journal of Hydrology, v. 387, n. 1, p. 33-45, 2010. http://dx.doi.org/10.1016/j.jhydrol.2010.03.027.
http://dx.doi.org/10.1016/j.jhydrol.2010...
(Equation 6) and the formulation proposed by Almeida et al. (2012)ALMEIDA, G. A. M.; BATES, P.; FREER, J. E.; SOUVIGNET, M. Improving the stability of a simple formulation of the shallow water equations for 2-D flood modeling. Water Resources Research, v. 48, n. 5, p. W05528, 2012. http://dx.doi.org/10.1029/2011WR011570.
http://dx.doi.org/10.1029/2011WR011570...
with a weighting factor equal to θ = 0.9 (Equation 11), adopted based on the authors' results, were tested with and without the hydraulic radius simplification.

All stable schemes presented values similar to the HEC-RAS model since the highest RMSE was 51.48 m3.s-1 for the Prog Q(k+1)xQ(k+1) scheme. Figure 1 shows only the region of the hydrograph peak flow in the river’s last reach for the original scheme, the Cent Q(k+1)xQ(k-1), Prog Q(k+1)xQ(k+1), Prog RK2 Q(k)xQ(k) schemes and the HEC-RAS model, with and without the hydraulic radius simplification. Figure 2 shows the results of the parameters evaluated by varying the time step. In the first column, the results correspond to simulations with the hydraulic radius simplification, and in the second column without. All values are shown in Table A1 of Appendix A Appendix A Supplementary data. Table A1 Results of test 1 (real reach) for different numerical schemes. Method dt = 1 min (C = 0.26) dt = 2 min (C = 0.51) dt = 3 min (C = 0.77) dt = 4 min (C =1.03) Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) HEC-RAS 3399.87 -0.02 3387.35 -0.02 3374.96 -0.02 3362.68 -0.02 Almeida et al. (2012) 3444.76 0.00 1.76 49.49 3444.21 0.00 0.86 48.73 3443.95 0.00 0.59 47.93 -- -- -- -- Bates, Horritt and Fewtrell (2010) 3443.79 0.00 1.27 50.59 3443.73 0.00 0.62 49.27 3443.62 0.00 0.41 48.28 -- -- -- -- Bates, Horritt and Fewtrell (2010) - R 3411.96 0.05 1.25 31.54 3411.91 0.05 0.63 29.84 3411.78 0.05 0.42 28.60 -- -- -- -- Cent Q(k+1) x Q(k-1) 3430.70 0.00 1.60 51.42 -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) - R 3399.17 0.05 1.64 33.06 -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) 3457.21 0.00 1.64 49.85 3470.77 0.00 0.81 48.50 3484.48 0.00 0.58 48.24 3498.38 0.00 0.42 49.08 Prog Q(k+1) x Q(k+1) - R 3425.07 0.05 1.66 30.20 3438.31 0.05 0.85 28.26 3451.66 0.05 0.56 28.13 3465.25 0.05 0.42 29.84 Cent Q(k+1) x Q(k+1) 3457.43 0.00 1.69 49.76 -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k+1) - R 3425.27 0.05 1.71 30.12 -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k) x Q(k) 3430.56 0.00 1.60 51.48 3417.47 0.00 0.79 50.79 3404.51 0.00 0.54 50.15 -- -- -- -- Prog Q(k) x Q(k) - R 3399.04 0.05 1.74 33.11 3386.25 0.05 0.87 32.60 3373.57 0.05 0.57 32.13 -- -- -- -- Prog RK2 Q(k) x Q(k) 3445.48 0.00 3.04 50.47 3451.58 0.00 1.31 48.94 3465.15 0.00 0.80 48.27 3496.15 0.00 0.58 53.43 Prog RK2 Q(k) x Q(k) - R 3413.62 0.05 3.09 31.31 3419.67 0.05 1.30 29.08 3433.14 0.05 0.82 27.96 3464.26 0.05 0.58 35.80 Prog RK4 Q(k) x Q(k) 3444.29 0.00 2.96 50.55 3445.80 0.00 1.48 49.17 3448.54 0.00 1.00 48.17 -- -- -- -- Prog RK4 Q(k) x Q(k) - R 3412.46 0.05 3.09 31.47 3413.95 0.05 1.56 29.63 3416.67 0.05 1.02 28.28 -- -- -- -- Table A2 Results of test 2 (low slope) for different numerical schemes. Method dt = 1 min (C = 0.29) dt = 2 min (C = 0.60) dt = 3 min (C = 0.90) dt = 3,5 min (C = 1.04) Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) HEC-RAS 1562.90 0.64 1561.60 0.64 1560.30 0.64 1559.00 0.64 Almeida et al. (2012) 1617.24 0.57 1.73 19.34 1617.63 0.57 0.89 19.52 1617.76 0.57 0.59 19.58 -- -- -- -- Bates, Horritt and Fewtrell (2010) 1618.03 0.57 1.26 18.51 1618.03 0.57 0.62 18.65 1618.02 0.57 0.43 18.81 -- -- -- -- Bates, Horritt and Fewtrell (2010) - R 1598.35 0.75 1.27 10.99 1598.34 0.75 0.64 11.23 1598.34 0.75 0.42 11.49 -- -- -- -- Cent Q(k+1) x Q(k-1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) 1619.57 0.57 1.51 18.69 1621.12 0.57 0.76 19.06 1622.67 0.56 0.51 19.50 1623.44 0.56 0.43 20.39 Prog Q(k+1) x Q(k+1) - R 1599.83 0.75 1.54 11.25 1601.31 0.74 0.77 11.83 1602.80 0.74 0.51 12.51 1603.55 0.74 0.45 13.12 Cent Q(k+1) x Q(k+1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k+1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k) x Q(k) 1616.50 0.57 1.53 18.36 1614.97 0.57 0.78 18.33 1613.45 0.57 0.51 18.29 -- -- -- -- Prog Q(k) x Q(k) - R 1596.87 0.75 1.53 10.76 1595.41 0.75 0.78 10.74 1593.95 0.75 0.51 10.72 -- -- -- -- Prog RK2 Q(k) x Q(k) 1618.10 0.57 2.93 18.52 1618.34 0.57 1.24 18.68 1618.78 0.57 0.76 18.89 -- -- -- -- Prog RK2 Q(k) x Q(k) - R 1598.42 0.75 3.01 11.00 1598.65 0.75 1.29 11.28 1599.07 0.75 0.78 11.64 -- -- -- -- Prog RK4 Q(k) x Q(k) 1618.05 0.57 2.19 18.52 1618.12 0.57 1.11 18.66 1618.24 0.57 0.73 18.83 -- -- -- -- Prog RK4 Q(k) x Q(k) - R 1598.37 0.75 2.27 10.99 1598.44 0.75 1.17 11.24 1598.55 0.75 0.77 11.53 -- -- -- -- Table A3 Results of test 3 (high slope) for different numerical schemes. Method dt = 1 min (C = 0.18) dt = 2 min (C = 0.37) dt = 3 min (C = 0.55) dt = 4 min (C =0.74) Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) HEC-RAS 4905.42 -0.07 4905.42 -0.07 NaN 4905.42 -0.07 4905.42 -0.07 Almeida et al. (2012) 4906.23 0.00 1.76 57.33 4906.02 0.00 0.88 57.69 -- -- -- -- -- -- -- -- Bates, Horritt and Fewtrell (2010) 4905.97 0.00 1.27 57.74 4905.89 0.00 0.64 57.89 -- -- -- -- -- -- -- -- Bates, Horritt and Fewtrell (2010) - R 4903.83 0.01 1.30 47.64 4903.70 0.01 0.64 47.80 -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) 4889.26 0.00 1.55 61.29 4872.40 0.00 0.76 65.21 -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) - R 4886.97 0.01 1.54 52.09 4869.91 0.01 0.77 56.80 -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) 4923.00 0.00 1.55 54.33 4941.11 0.00 0.77 51.38 4963.35 0.00 0.51 49.14 5000.27 0.00 0.38 48.44 Prog Q(k+1) x Q(k+1) - R 4920.96 0.01 1.54 43.26 4939.08 0.01 0.79 39.31 4960.83 0.01 0.51 36.25 4995.77 0.01 0.38 35.34 Cent Q(k+1) x Q(k+1) 4923.14 0.00 1.53 54.24 4941.25 0.00 0.76 51.17 4961.23 0.00 0.51 48.80 -- -- -- -- Cent Q(k+1) x Q(k+1) - R 4921.11 0.01 1.54 43.19 4939.21 0.01 0.78 39.17 4959.08 0.01 0.51 36.02 -- -- -- -- Prog Q(k) x Q(k) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k) x Q(k) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK2 Q(k) x Q(k) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK2 Q(k) x Q(k) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK4 Q(k) x Q(k) 4908.49 0.00 2.19 57.54 -- -- -- -- -- -- -- -- -- -- -- -- Prog RK4 Q(k) x Q(k) - R 4906.39 0.01 2.26 47.12 -- -- -- -- -- -- -- -- -- -- -- -- Table A4 Results of test 4 (very high slope) for different numerical schemes. Method dt = 80 seg (C = 0.20) dt = 120 seg (C = 0.31) dt = 160 seg (C = 0.41) dt = 200 seg (C ~ 0.53) Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) HEC-RAS 4910.30 -0.09 4910.30 -0.09 4910.30 -0.09 4910.30 -0.09 Almeida et al. (2012) 4942.79 0.00 1.20 30.18 -- -- -- -- -- -- -- -- -- -- -- -- Bates, Horritt and Fewtrell (2010) 4942.78 0.00 1.12 30.19 -- -- -- -- -- -- -- -- -- -- -- -- Bates, Horritt and Fewtrell (2010) - R 4942.18 0.01 1.15 25.75 -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) 4963.23 0.00 1.15 31.24 4978.55 0.00 0.77 32.18 5005.51 0.00 0.58 33.48 -- -- -- -- Prog Q(k+1) x Q(k+1) - R 4962.55 0.01 1.19 26.92 4977.66 0.01 0.78 27.97 5003.56 0.01 0.57 29.41 -- -- -- -- Cent Q(k+1) x Q(k+1) 4963.04 0.00 1.14 31.23 4977.23 0.00 0.78 32.16 4997.92 0.00 0.57 33.41 -- -- -- -- Cent Q(k+1) x Q(k+1) - R 4962.37 0.01 1.16 26.91 4976.51 0.01 0.79 27.94 4996.96 0.01 0.58 29.35 -- -- -- -- Prog Q(k) x Q(k) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k) x Q(k) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK2 Q(k) x Q(k) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK2 Q(k) x Q(k) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK4 Q(k) x Q(k) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK4 Q(k) x Q(k) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Table A5 Results of test 5 (backwater) for different numerical schemes. Method dt = 1 min (C = 0.44) dt = 2 min (C = 0.88) dt = 2,2 min (C = 0.97) dt = 2,5 min (C > 1) Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) HEC-RAS 4011.64 0.00 4011.64 0.00 4011.64 0.00 4011.64 0.00 Almeida et al. (2012) 4092.22 0.00 1.82 21.10 4093.60 0.00 0.91 21.41 4093.67 0.00 0.83 21.44 -- -- -- -- Bates, Horritt and Fewtrell (2010) 4095.07 0.00 1.30 21.82 4095.02 0.00 0.65 21.81 4094.96 0.00 0.59 21.81 -- -- -- -- Bates, Horritt and Fewtrell (2010) - R 4051.15 0.00 1.29 15.76 4051.11 0.00 0.65 15.72 4051.06 0.00 0.59 15.71 -- -- -- -- Cent Q(k+1) x Q(k-1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) 4109.87 0.00 1.69 24.04 4124.91 0.00 0.83 26.80 4127.91 0.00 0.70 27.42 -- -- -- -- Prog Q(k+1) x Q(k+1) - R 4065.47 0.00 1.72 18.30 4079.97 0.00 0.86 21.42 4082.89 0.00 0.72 22.10 -- -- -- -- Cent Q(k+1) x Q(k+1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k+1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k) x Q(k) 4080.54 0.00 1.66 20.29 4066.21 0.00 0.84 19.55 4063.34 0.00 0.77 19.50 -- -- -- -- Prog Q(k) x Q(k) - R 4037.06 0.00 1.76 14.09 4023.14 0.00 0.87 13.49 4020.36 0.00 0.79 13.51 -- -- -- -- Prog RK2 Q(k) x Q(k) 4096.34 0.00 3.07 22.00 4100.76 0.00 1.29 22.73 4102.09 0.00 1.13 23.00 -- -- -- -- Prog RK2 Q(k) x Q(k) - R 4052.42 0.00 3.16 15.90 4056.87 0.00 1.34 16.42 4058.23 0.00 1.14 16.61 -- -- -- -- Prog RK4 Q(k) x Q(k) 4095.46 0.00 2.96 21.87 4096.61 0.00 1.48 22.04 4096.89 0.00 1.00 22.10 -- -- -- -- Prog RK4 Q(k) x Q(k) - R 4051.54 0.00 3.07 15.81 4052.69 0.00 1.52 15.91 4052.99 0.00 1.26 15.94 -- -- -- -- Table A6 Results of test 6 (tide) for different numerical schemes. Method dt = 1 min (C = 0.27) dt = 2 min (C = 0.54) dt = 3,5 min (C = 0.94) dt = 4 min (C = 1.07) Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) HEC-RAS 2125.02 -0.84 2125.02 -0.84 2125.02 -0.84 2125.02 -0.84 Almeida et al. (2012) 2141.89 0.94 1.88 32.73 2140.17 0.95 0.90 31.47 2139.53 1.01 0.52 33.63 -- -- -- -- Bates, Horritt and Fewtrell (2010) 2138.46 0.97 1.33 30.99 2138.44 0.97 0.65 31.20 2138.54 1.02 0.37 33.90 -- -- -- -- Bates, Horritt and Fewtrell (2010) - R 2123.32 1.12 1.32 33.53 2123.32 1.12 0.65 33.44 2123.37 1.17 0.37 35.54 -- -- -- -- Cent Q(k+1) x Q(k-1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k+1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k+1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k) x Q(k) 2138.76 0.96 1.68 31.13 2139.04 0.95 0.85 30.70 2139.45 0.99 0.49 31.17 -- -- -- -- Prog Q(k) x Q(k) - R 2123.59 1.11 1.76 33.59 2123.85 1.10 0.88 32.85 2124.25 1.14 0.51 32.74 -- -- -- -- Prog RK2 Q(k) x Q(k) 2138.45 0.97 3.10 31.01 2138.43 0.97 1.30 31.31 2138.46 1.04 0.70 34.43 2138.49 0.99 0.60 36.46 Prog RK2 Q(k) x Q(k) - R 2123.32 1.12 3.18 33.55 2123.31 1.12 1.32 33.54 2123.29 1.19 0.70 36.03 2123.37 1.14 0.63 37.74 Prog RK4 Q(k) x Q(k) 2138.46 0.97 2.95 31.00 2138.44 0.97 1.46 31.24 2138.51 1.03 0.86 34.07 -- -- -- -- Prog RK4 Q(k) x Q(k) - R 2123.32 1.12 3.03 33.54 2123.31 1.12 1.52 33.48 2123.34 1.18 0.88 35.72 -- -- -- -- .

Figure 1
Region of the hydrograph peak flow in the river's last reach for the Bates, Horritt and Fewtrell (2010)BATES, P. D.; HORRITT, M. S.; FEWTRELL, T. J. A simple inertial formulation of the shallow water equations for efficient two-dimensional flood inundation modelling. Journal of Hydrology, v. 387, n. 1, p. 33-45, 2010. http://dx.doi.org/10.1016/j.jhydrol.2010.03.027.
http://dx.doi.org/10.1016/j.jhydrol.2010...
, Cent Q(k+1)xQ(k-1), Prog Q(k+1)xQ(k+1) and Prog RK2 Q(k)xQ(k) schemes with and without the hydraulic radius simplification considering a river's slope equal to 29.5 cm/km.
Figure 2
Peak flow, volume conservation error, processing time and RMSE in relation to the time step for different formulations of numerical schemes considering a river's slope equal to 29.5 cm/km.

The Prog Q(k+1)xQ(k+1) and Prog RK2 Q(k)xQ(k) schemes were the most stable, without instability until C = 1.03. However, in comparison to the Bates, Horritt and Fewtrell (2010)BATES, P. D.; HORRITT, M. S.; FEWTRELL, T. J. A simple inertial formulation of the shallow water equations for efficient two-dimensional flood inundation modelling. Journal of Hydrology, v. 387, n. 1, p. 33-45, 2010. http://dx.doi.org/10.1016/j.jhydrol.2010.03.027.
http://dx.doi.org/10.1016/j.jhydrol.2010...
model, the processing time was on average 33% higher for the Prog Q(k+1)xQ(k+1) scheme and twice for the Prog RK2 Q(k)xQ(k) scheme.

The hydraulic radius calculation without simplification did not increase the processing time of the schemes and their use reduced RMSE. As an example, the scheme of Bates, Horritt and Fewtrell (2010)BATES, P. D.; HORRITT, M. S.; FEWTRELL, T. J. A simple inertial formulation of the shallow water equations for efficient two-dimensional flood inundation modelling. Journal of Hydrology, v. 387, n. 1, p. 33-45, 2010. http://dx.doi.org/10.1016/j.jhydrol.2010.03.027.
http://dx.doi.org/10.1016/j.jhydrol.2010...
presented a 60% lower RMSE when the hydraulic radius was calculated. Despite this, the volume conservation error increased from 0% to 0.05%, but still remained negligible. This parameter (ErrVol) did not vary with the time step nor with the numerical scheme. It was also observed that the calculation of hydraulic radius did not influence the numerical schemes' stability.

The Almeida et al. (2012)ALMEIDA, G. A. M.; BATES, P.; FREER, J. E.; SOUVIGNET, M. Improving the stability of a simple formulation of the shallow water equations for 2-D flood modeling. Water Resources Research, v. 48, n. 5, p. W05528, 2012. http://dx.doi.org/10.1029/2011WR011570.
http://dx.doi.org/10.1029/2011WR011570...
scheme presented no advantages compared to the Bates, Horritt and Fewtrell (2010)BATES, P. D.; HORRITT, M. S.; FEWTRELL, T. J. A simple inertial formulation of the shallow water equations for efficient two-dimensional flood inundation modelling. Journal of Hydrology, v. 387, n. 1, p. 33-45, 2010. http://dx.doi.org/10.1016/j.jhydrol.2010.03.027.
http://dx.doi.org/10.1016/j.jhydrol.2010...
scheme in terms of stability (see Table A1 of Appendix A Appendix A Supplementary data. Table A1 Results of test 1 (real reach) for different numerical schemes. Method dt = 1 min (C = 0.26) dt = 2 min (C = 0.51) dt = 3 min (C = 0.77) dt = 4 min (C =1.03) Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) HEC-RAS 3399.87 -0.02 3387.35 -0.02 3374.96 -0.02 3362.68 -0.02 Almeida et al. (2012) 3444.76 0.00 1.76 49.49 3444.21 0.00 0.86 48.73 3443.95 0.00 0.59 47.93 -- -- -- -- Bates, Horritt and Fewtrell (2010) 3443.79 0.00 1.27 50.59 3443.73 0.00 0.62 49.27 3443.62 0.00 0.41 48.28 -- -- -- -- Bates, Horritt and Fewtrell (2010) - R 3411.96 0.05 1.25 31.54 3411.91 0.05 0.63 29.84 3411.78 0.05 0.42 28.60 -- -- -- -- Cent Q(k+1) x Q(k-1) 3430.70 0.00 1.60 51.42 -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) - R 3399.17 0.05 1.64 33.06 -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) 3457.21 0.00 1.64 49.85 3470.77 0.00 0.81 48.50 3484.48 0.00 0.58 48.24 3498.38 0.00 0.42 49.08 Prog Q(k+1) x Q(k+1) - R 3425.07 0.05 1.66 30.20 3438.31 0.05 0.85 28.26 3451.66 0.05 0.56 28.13 3465.25 0.05 0.42 29.84 Cent Q(k+1) x Q(k+1) 3457.43 0.00 1.69 49.76 -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k+1) - R 3425.27 0.05 1.71 30.12 -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k) x Q(k) 3430.56 0.00 1.60 51.48 3417.47 0.00 0.79 50.79 3404.51 0.00 0.54 50.15 -- -- -- -- Prog Q(k) x Q(k) - R 3399.04 0.05 1.74 33.11 3386.25 0.05 0.87 32.60 3373.57 0.05 0.57 32.13 -- -- -- -- Prog RK2 Q(k) x Q(k) 3445.48 0.00 3.04 50.47 3451.58 0.00 1.31 48.94 3465.15 0.00 0.80 48.27 3496.15 0.00 0.58 53.43 Prog RK2 Q(k) x Q(k) - R 3413.62 0.05 3.09 31.31 3419.67 0.05 1.30 29.08 3433.14 0.05 0.82 27.96 3464.26 0.05 0.58 35.80 Prog RK4 Q(k) x Q(k) 3444.29 0.00 2.96 50.55 3445.80 0.00 1.48 49.17 3448.54 0.00 1.00 48.17 -- -- -- -- Prog RK4 Q(k) x Q(k) - R 3412.46 0.05 3.09 31.47 3413.95 0.05 1.56 29.63 3416.67 0.05 1.02 28.28 -- -- -- -- Table A2 Results of test 2 (low slope) for different numerical schemes. Method dt = 1 min (C = 0.29) dt = 2 min (C = 0.60) dt = 3 min (C = 0.90) dt = 3,5 min (C = 1.04) Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) HEC-RAS 1562.90 0.64 1561.60 0.64 1560.30 0.64 1559.00 0.64 Almeida et al. (2012) 1617.24 0.57 1.73 19.34 1617.63 0.57 0.89 19.52 1617.76 0.57 0.59 19.58 -- -- -- -- Bates, Horritt and Fewtrell (2010) 1618.03 0.57 1.26 18.51 1618.03 0.57 0.62 18.65 1618.02 0.57 0.43 18.81 -- -- -- -- Bates, Horritt and Fewtrell (2010) - R 1598.35 0.75 1.27 10.99 1598.34 0.75 0.64 11.23 1598.34 0.75 0.42 11.49 -- -- -- -- Cent Q(k+1) x Q(k-1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) 1619.57 0.57 1.51 18.69 1621.12 0.57 0.76 19.06 1622.67 0.56 0.51 19.50 1623.44 0.56 0.43 20.39 Prog Q(k+1) x Q(k+1) - R 1599.83 0.75 1.54 11.25 1601.31 0.74 0.77 11.83 1602.80 0.74 0.51 12.51 1603.55 0.74 0.45 13.12 Cent Q(k+1) x Q(k+1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k+1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k) x Q(k) 1616.50 0.57 1.53 18.36 1614.97 0.57 0.78 18.33 1613.45 0.57 0.51 18.29 -- -- -- -- Prog Q(k) x Q(k) - R 1596.87 0.75 1.53 10.76 1595.41 0.75 0.78 10.74 1593.95 0.75 0.51 10.72 -- -- -- -- Prog RK2 Q(k) x Q(k) 1618.10 0.57 2.93 18.52 1618.34 0.57 1.24 18.68 1618.78 0.57 0.76 18.89 -- -- -- -- Prog RK2 Q(k) x Q(k) - R 1598.42 0.75 3.01 11.00 1598.65 0.75 1.29 11.28 1599.07 0.75 0.78 11.64 -- -- -- -- Prog RK4 Q(k) x Q(k) 1618.05 0.57 2.19 18.52 1618.12 0.57 1.11 18.66 1618.24 0.57 0.73 18.83 -- -- -- -- Prog RK4 Q(k) x Q(k) - R 1598.37 0.75 2.27 10.99 1598.44 0.75 1.17 11.24 1598.55 0.75 0.77 11.53 -- -- -- -- Table A3 Results of test 3 (high slope) for different numerical schemes. Method dt = 1 min (C = 0.18) dt = 2 min (C = 0.37) dt = 3 min (C = 0.55) dt = 4 min (C =0.74) Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) HEC-RAS 4905.42 -0.07 4905.42 -0.07 NaN 4905.42 -0.07 4905.42 -0.07 Almeida et al. (2012) 4906.23 0.00 1.76 57.33 4906.02 0.00 0.88 57.69 -- -- -- -- -- -- -- -- Bates, Horritt and Fewtrell (2010) 4905.97 0.00 1.27 57.74 4905.89 0.00 0.64 57.89 -- -- -- -- -- -- -- -- Bates, Horritt and Fewtrell (2010) - R 4903.83 0.01 1.30 47.64 4903.70 0.01 0.64 47.80 -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) 4889.26 0.00 1.55 61.29 4872.40 0.00 0.76 65.21 -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) - R 4886.97 0.01 1.54 52.09 4869.91 0.01 0.77 56.80 -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) 4923.00 0.00 1.55 54.33 4941.11 0.00 0.77 51.38 4963.35 0.00 0.51 49.14 5000.27 0.00 0.38 48.44 Prog Q(k+1) x Q(k+1) - R 4920.96 0.01 1.54 43.26 4939.08 0.01 0.79 39.31 4960.83 0.01 0.51 36.25 4995.77 0.01 0.38 35.34 Cent Q(k+1) x Q(k+1) 4923.14 0.00 1.53 54.24 4941.25 0.00 0.76 51.17 4961.23 0.00 0.51 48.80 -- -- -- -- Cent Q(k+1) x Q(k+1) - R 4921.11 0.01 1.54 43.19 4939.21 0.01 0.78 39.17 4959.08 0.01 0.51 36.02 -- -- -- -- Prog Q(k) x Q(k) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k) x Q(k) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK2 Q(k) x Q(k) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK2 Q(k) x Q(k) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK4 Q(k) x Q(k) 4908.49 0.00 2.19 57.54 -- -- -- -- -- -- -- -- -- -- -- -- Prog RK4 Q(k) x Q(k) - R 4906.39 0.01 2.26 47.12 -- -- -- -- -- -- -- -- -- -- -- -- Table A4 Results of test 4 (very high slope) for different numerical schemes. Method dt = 80 seg (C = 0.20) dt = 120 seg (C = 0.31) dt = 160 seg (C = 0.41) dt = 200 seg (C ~ 0.53) Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) HEC-RAS 4910.30 -0.09 4910.30 -0.09 4910.30 -0.09 4910.30 -0.09 Almeida et al. (2012) 4942.79 0.00 1.20 30.18 -- -- -- -- -- -- -- -- -- -- -- -- Bates, Horritt and Fewtrell (2010) 4942.78 0.00 1.12 30.19 -- -- -- -- -- -- -- -- -- -- -- -- Bates, Horritt and Fewtrell (2010) - R 4942.18 0.01 1.15 25.75 -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) 4963.23 0.00 1.15 31.24 4978.55 0.00 0.77 32.18 5005.51 0.00 0.58 33.48 -- -- -- -- Prog Q(k+1) x Q(k+1) - R 4962.55 0.01 1.19 26.92 4977.66 0.01 0.78 27.97 5003.56 0.01 0.57 29.41 -- -- -- -- Cent Q(k+1) x Q(k+1) 4963.04 0.00 1.14 31.23 4977.23 0.00 0.78 32.16 4997.92 0.00 0.57 33.41 -- -- -- -- Cent Q(k+1) x Q(k+1) - R 4962.37 0.01 1.16 26.91 4976.51 0.01 0.79 27.94 4996.96 0.01 0.58 29.35 -- -- -- -- Prog Q(k) x Q(k) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k) x Q(k) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK2 Q(k) x Q(k) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK2 Q(k) x Q(k) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK4 Q(k) x Q(k) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK4 Q(k) x Q(k) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Table A5 Results of test 5 (backwater) for different numerical schemes. Method dt = 1 min (C = 0.44) dt = 2 min (C = 0.88) dt = 2,2 min (C = 0.97) dt = 2,5 min (C > 1) Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) HEC-RAS 4011.64 0.00 4011.64 0.00 4011.64 0.00 4011.64 0.00 Almeida et al. (2012) 4092.22 0.00 1.82 21.10 4093.60 0.00 0.91 21.41 4093.67 0.00 0.83 21.44 -- -- -- -- Bates, Horritt and Fewtrell (2010) 4095.07 0.00 1.30 21.82 4095.02 0.00 0.65 21.81 4094.96 0.00 0.59 21.81 -- -- -- -- Bates, Horritt and Fewtrell (2010) - R 4051.15 0.00 1.29 15.76 4051.11 0.00 0.65 15.72 4051.06 0.00 0.59 15.71 -- -- -- -- Cent Q(k+1) x Q(k-1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) 4109.87 0.00 1.69 24.04 4124.91 0.00 0.83 26.80 4127.91 0.00 0.70 27.42 -- -- -- -- Prog Q(k+1) x Q(k+1) - R 4065.47 0.00 1.72 18.30 4079.97 0.00 0.86 21.42 4082.89 0.00 0.72 22.10 -- -- -- -- Cent Q(k+1) x Q(k+1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k+1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k) x Q(k) 4080.54 0.00 1.66 20.29 4066.21 0.00 0.84 19.55 4063.34 0.00 0.77 19.50 -- -- -- -- Prog Q(k) x Q(k) - R 4037.06 0.00 1.76 14.09 4023.14 0.00 0.87 13.49 4020.36 0.00 0.79 13.51 -- -- -- -- Prog RK2 Q(k) x Q(k) 4096.34 0.00 3.07 22.00 4100.76 0.00 1.29 22.73 4102.09 0.00 1.13 23.00 -- -- -- -- Prog RK2 Q(k) x Q(k) - R 4052.42 0.00 3.16 15.90 4056.87 0.00 1.34 16.42 4058.23 0.00 1.14 16.61 -- -- -- -- Prog RK4 Q(k) x Q(k) 4095.46 0.00 2.96 21.87 4096.61 0.00 1.48 22.04 4096.89 0.00 1.00 22.10 -- -- -- -- Prog RK4 Q(k) x Q(k) - R 4051.54 0.00 3.07 15.81 4052.69 0.00 1.52 15.91 4052.99 0.00 1.26 15.94 -- -- -- -- Table A6 Results of test 6 (tide) for different numerical schemes. Method dt = 1 min (C = 0.27) dt = 2 min (C = 0.54) dt = 3,5 min (C = 0.94) dt = 4 min (C = 1.07) Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) HEC-RAS 2125.02 -0.84 2125.02 -0.84 2125.02 -0.84 2125.02 -0.84 Almeida et al. (2012) 2141.89 0.94 1.88 32.73 2140.17 0.95 0.90 31.47 2139.53 1.01 0.52 33.63 -- -- -- -- Bates, Horritt and Fewtrell (2010) 2138.46 0.97 1.33 30.99 2138.44 0.97 0.65 31.20 2138.54 1.02 0.37 33.90 -- -- -- -- Bates, Horritt and Fewtrell (2010) - R 2123.32 1.12 1.32 33.53 2123.32 1.12 0.65 33.44 2123.37 1.17 0.37 35.54 -- -- -- -- Cent Q(k+1) x Q(k-1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k+1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k+1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k) x Q(k) 2138.76 0.96 1.68 31.13 2139.04 0.95 0.85 30.70 2139.45 0.99 0.49 31.17 -- -- -- -- Prog Q(k) x Q(k) - R 2123.59 1.11 1.76 33.59 2123.85 1.10 0.88 32.85 2124.25 1.14 0.51 32.74 -- -- -- -- Prog RK2 Q(k) x Q(k) 2138.45 0.97 3.10 31.01 2138.43 0.97 1.30 31.31 2138.46 1.04 0.70 34.43 2138.49 0.99 0.60 36.46 Prog RK2 Q(k) x Q(k) - R 2123.32 1.12 3.18 33.55 2123.31 1.12 1.32 33.54 2123.29 1.19 0.70 36.03 2123.37 1.14 0.63 37.74 Prog RK4 Q(k) x Q(k) 2138.46 0.97 2.95 31.00 2138.44 0.97 1.46 31.24 2138.51 1.03 0.86 34.07 -- -- -- -- Prog RK4 Q(k) x Q(k) - R 2123.32 1.12 3.03 33.54 2123.31 1.12 1.52 33.48 2123.34 1.18 0.88 35.72 -- -- -- -- ). These schemes were unstable with a time step of 3 min (C = 0.77). All schemes with centered derivative presented instability, with stable results only for the time step of 1 minute (C = 0.26). In addition, these schemes did not improve the RMSE and the volume conservation error, as was expected for a second-order approximation scheme with respect to a first-order scheme.

Test 2

This test considered the same conditions of the river reach used in test 1 but it used a river’s low slope of 5 cm/km. The time step was varied in 1, 2, 3 and 3.5 minutes, corresponding to C values of 0.29, 0.60, 1.90 and 1.04, respectively, at the time of greater depth. As in test 1, the numerical schemes' performance was compared. Figure 3 shows the peak flow region in the last reach for the original scheme, the Prog Q(k+1)xQ(k+1), Prog RK2 Q(k)xQ(k) schemes and the HEC-RAS model. Figure 4 shows the performance of all schemes evaluated and Table A2 ( Appendix A Appendix A Supplementary data. Table A1 Results of test 1 (real reach) for different numerical schemes. Method dt = 1 min (C = 0.26) dt = 2 min (C = 0.51) dt = 3 min (C = 0.77) dt = 4 min (C =1.03) Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) HEC-RAS 3399.87 -0.02 3387.35 -0.02 3374.96 -0.02 3362.68 -0.02 Almeida et al. (2012) 3444.76 0.00 1.76 49.49 3444.21 0.00 0.86 48.73 3443.95 0.00 0.59 47.93 -- -- -- -- Bates, Horritt and Fewtrell (2010) 3443.79 0.00 1.27 50.59 3443.73 0.00 0.62 49.27 3443.62 0.00 0.41 48.28 -- -- -- -- Bates, Horritt and Fewtrell (2010) - R 3411.96 0.05 1.25 31.54 3411.91 0.05 0.63 29.84 3411.78 0.05 0.42 28.60 -- -- -- -- Cent Q(k+1) x Q(k-1) 3430.70 0.00 1.60 51.42 -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) - R 3399.17 0.05 1.64 33.06 -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) 3457.21 0.00 1.64 49.85 3470.77 0.00 0.81 48.50 3484.48 0.00 0.58 48.24 3498.38 0.00 0.42 49.08 Prog Q(k+1) x Q(k+1) - R 3425.07 0.05 1.66 30.20 3438.31 0.05 0.85 28.26 3451.66 0.05 0.56 28.13 3465.25 0.05 0.42 29.84 Cent Q(k+1) x Q(k+1) 3457.43 0.00 1.69 49.76 -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k+1) - R 3425.27 0.05 1.71 30.12 -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k) x Q(k) 3430.56 0.00 1.60 51.48 3417.47 0.00 0.79 50.79 3404.51 0.00 0.54 50.15 -- -- -- -- Prog Q(k) x Q(k) - R 3399.04 0.05 1.74 33.11 3386.25 0.05 0.87 32.60 3373.57 0.05 0.57 32.13 -- -- -- -- Prog RK2 Q(k) x Q(k) 3445.48 0.00 3.04 50.47 3451.58 0.00 1.31 48.94 3465.15 0.00 0.80 48.27 3496.15 0.00 0.58 53.43 Prog RK2 Q(k) x Q(k) - R 3413.62 0.05 3.09 31.31 3419.67 0.05 1.30 29.08 3433.14 0.05 0.82 27.96 3464.26 0.05 0.58 35.80 Prog RK4 Q(k) x Q(k) 3444.29 0.00 2.96 50.55 3445.80 0.00 1.48 49.17 3448.54 0.00 1.00 48.17 -- -- -- -- Prog RK4 Q(k) x Q(k) - R 3412.46 0.05 3.09 31.47 3413.95 0.05 1.56 29.63 3416.67 0.05 1.02 28.28 -- -- -- -- Table A2 Results of test 2 (low slope) for different numerical schemes. Method dt = 1 min (C = 0.29) dt = 2 min (C = 0.60) dt = 3 min (C = 0.90) dt = 3,5 min (C = 1.04) Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) HEC-RAS 1562.90 0.64 1561.60 0.64 1560.30 0.64 1559.00 0.64 Almeida et al. (2012) 1617.24 0.57 1.73 19.34 1617.63 0.57 0.89 19.52 1617.76 0.57 0.59 19.58 -- -- -- -- Bates, Horritt and Fewtrell (2010) 1618.03 0.57 1.26 18.51 1618.03 0.57 0.62 18.65 1618.02 0.57 0.43 18.81 -- -- -- -- Bates, Horritt and Fewtrell (2010) - R 1598.35 0.75 1.27 10.99 1598.34 0.75 0.64 11.23 1598.34 0.75 0.42 11.49 -- -- -- -- Cent Q(k+1) x Q(k-1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) 1619.57 0.57 1.51 18.69 1621.12 0.57 0.76 19.06 1622.67 0.56 0.51 19.50 1623.44 0.56 0.43 20.39 Prog Q(k+1) x Q(k+1) - R 1599.83 0.75 1.54 11.25 1601.31 0.74 0.77 11.83 1602.80 0.74 0.51 12.51 1603.55 0.74 0.45 13.12 Cent Q(k+1) x Q(k+1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k+1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k) x Q(k) 1616.50 0.57 1.53 18.36 1614.97 0.57 0.78 18.33 1613.45 0.57 0.51 18.29 -- -- -- -- Prog Q(k) x Q(k) - R 1596.87 0.75 1.53 10.76 1595.41 0.75 0.78 10.74 1593.95 0.75 0.51 10.72 -- -- -- -- Prog RK2 Q(k) x Q(k) 1618.10 0.57 2.93 18.52 1618.34 0.57 1.24 18.68 1618.78 0.57 0.76 18.89 -- -- -- -- Prog RK2 Q(k) x Q(k) - R 1598.42 0.75 3.01 11.00 1598.65 0.75 1.29 11.28 1599.07 0.75 0.78 11.64 -- -- -- -- Prog RK4 Q(k) x Q(k) 1618.05 0.57 2.19 18.52 1618.12 0.57 1.11 18.66 1618.24 0.57 0.73 18.83 -- -- -- -- Prog RK4 Q(k) x Q(k) - R 1598.37 0.75 2.27 10.99 1598.44 0.75 1.17 11.24 1598.55 0.75 0.77 11.53 -- -- -- -- Table A3 Results of test 3 (high slope) for different numerical schemes. Method dt = 1 min (C = 0.18) dt = 2 min (C = 0.37) dt = 3 min (C = 0.55) dt = 4 min (C =0.74) Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) HEC-RAS 4905.42 -0.07 4905.42 -0.07 NaN 4905.42 -0.07 4905.42 -0.07 Almeida et al. (2012) 4906.23 0.00 1.76 57.33 4906.02 0.00 0.88 57.69 -- -- -- -- -- -- -- -- Bates, Horritt and Fewtrell (2010) 4905.97 0.00 1.27 57.74 4905.89 0.00 0.64 57.89 -- -- -- -- -- -- -- -- Bates, Horritt and Fewtrell (2010) - R 4903.83 0.01 1.30 47.64 4903.70 0.01 0.64 47.80 -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) 4889.26 0.00 1.55 61.29 4872.40 0.00 0.76 65.21 -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) - R 4886.97 0.01 1.54 52.09 4869.91 0.01 0.77 56.80 -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) 4923.00 0.00 1.55 54.33 4941.11 0.00 0.77 51.38 4963.35 0.00 0.51 49.14 5000.27 0.00 0.38 48.44 Prog Q(k+1) x Q(k+1) - R 4920.96 0.01 1.54 43.26 4939.08 0.01 0.79 39.31 4960.83 0.01 0.51 36.25 4995.77 0.01 0.38 35.34 Cent Q(k+1) x Q(k+1) 4923.14 0.00 1.53 54.24 4941.25 0.00 0.76 51.17 4961.23 0.00 0.51 48.80 -- -- -- -- Cent Q(k+1) x Q(k+1) - R 4921.11 0.01 1.54 43.19 4939.21 0.01 0.78 39.17 4959.08 0.01 0.51 36.02 -- -- -- -- Prog Q(k) x Q(k) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k) x Q(k) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK2 Q(k) x Q(k) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK2 Q(k) x Q(k) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK4 Q(k) x Q(k) 4908.49 0.00 2.19 57.54 -- -- -- -- -- -- -- -- -- -- -- -- Prog RK4 Q(k) x Q(k) - R 4906.39 0.01 2.26 47.12 -- -- -- -- -- -- -- -- -- -- -- -- Table A4 Results of test 4 (very high slope) for different numerical schemes. Method dt = 80 seg (C = 0.20) dt = 120 seg (C = 0.31) dt = 160 seg (C = 0.41) dt = 200 seg (C ~ 0.53) Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) HEC-RAS 4910.30 -0.09 4910.30 -0.09 4910.30 -0.09 4910.30 -0.09 Almeida et al. (2012) 4942.79 0.00 1.20 30.18 -- -- -- -- -- -- -- -- -- -- -- -- Bates, Horritt and Fewtrell (2010) 4942.78 0.00 1.12 30.19 -- -- -- -- -- -- -- -- -- -- -- -- Bates, Horritt and Fewtrell (2010) - R 4942.18 0.01 1.15 25.75 -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) 4963.23 0.00 1.15 31.24 4978.55 0.00 0.77 32.18 5005.51 0.00 0.58 33.48 -- -- -- -- Prog Q(k+1) x Q(k+1) - R 4962.55 0.01 1.19 26.92 4977.66 0.01 0.78 27.97 5003.56 0.01 0.57 29.41 -- -- -- -- Cent Q(k+1) x Q(k+1) 4963.04 0.00 1.14 31.23 4977.23 0.00 0.78 32.16 4997.92 0.00 0.57 33.41 -- -- -- -- Cent Q(k+1) x Q(k+1) - R 4962.37 0.01 1.16 26.91 4976.51 0.01 0.79 27.94 4996.96 0.01 0.58 29.35 -- -- -- -- Prog Q(k) x Q(k) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k) x Q(k) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK2 Q(k) x Q(k) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK2 Q(k) x Q(k) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK4 Q(k) x Q(k) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK4 Q(k) x Q(k) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Table A5 Results of test 5 (backwater) for different numerical schemes. Method dt = 1 min (C = 0.44) dt = 2 min (C = 0.88) dt = 2,2 min (C = 0.97) dt = 2,5 min (C > 1) Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) HEC-RAS 4011.64 0.00 4011.64 0.00 4011.64 0.00 4011.64 0.00 Almeida et al. (2012) 4092.22 0.00 1.82 21.10 4093.60 0.00 0.91 21.41 4093.67 0.00 0.83 21.44 -- -- -- -- Bates, Horritt and Fewtrell (2010) 4095.07 0.00 1.30 21.82 4095.02 0.00 0.65 21.81 4094.96 0.00 0.59 21.81 -- -- -- -- Bates, Horritt and Fewtrell (2010) - R 4051.15 0.00 1.29 15.76 4051.11 0.00 0.65 15.72 4051.06 0.00 0.59 15.71 -- -- -- -- Cent Q(k+1) x Q(k-1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) 4109.87 0.00 1.69 24.04 4124.91 0.00 0.83 26.80 4127.91 0.00 0.70 27.42 -- -- -- -- Prog Q(k+1) x Q(k+1) - R 4065.47 0.00 1.72 18.30 4079.97 0.00 0.86 21.42 4082.89 0.00 0.72 22.10 -- -- -- -- Cent Q(k+1) x Q(k+1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k+1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k) x Q(k) 4080.54 0.00 1.66 20.29 4066.21 0.00 0.84 19.55 4063.34 0.00 0.77 19.50 -- -- -- -- Prog Q(k) x Q(k) - R 4037.06 0.00 1.76 14.09 4023.14 0.00 0.87 13.49 4020.36 0.00 0.79 13.51 -- -- -- -- Prog RK2 Q(k) x Q(k) 4096.34 0.00 3.07 22.00 4100.76 0.00 1.29 22.73 4102.09 0.00 1.13 23.00 -- -- -- -- Prog RK2 Q(k) x Q(k) - R 4052.42 0.00 3.16 15.90 4056.87 0.00 1.34 16.42 4058.23 0.00 1.14 16.61 -- -- -- -- Prog RK4 Q(k) x Q(k) 4095.46 0.00 2.96 21.87 4096.61 0.00 1.48 22.04 4096.89 0.00 1.00 22.10 -- -- -- -- Prog RK4 Q(k) x Q(k) - R 4051.54 0.00 3.07 15.81 4052.69 0.00 1.52 15.91 4052.99 0.00 1.26 15.94 -- -- -- -- Table A6 Results of test 6 (tide) for different numerical schemes. Method dt = 1 min (C = 0.27) dt = 2 min (C = 0.54) dt = 3,5 min (C = 0.94) dt = 4 min (C = 1.07) Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) HEC-RAS 2125.02 -0.84 2125.02 -0.84 2125.02 -0.84 2125.02 -0.84 Almeida et al. (2012) 2141.89 0.94 1.88 32.73 2140.17 0.95 0.90 31.47 2139.53 1.01 0.52 33.63 -- -- -- -- Bates, Horritt and Fewtrell (2010) 2138.46 0.97 1.33 30.99 2138.44 0.97 0.65 31.20 2138.54 1.02 0.37 33.90 -- -- -- -- Bates, Horritt and Fewtrell (2010) - R 2123.32 1.12 1.32 33.53 2123.32 1.12 0.65 33.44 2123.37 1.17 0.37 35.54 -- -- -- -- Cent Q(k+1) x Q(k-1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k+1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k+1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k) x Q(k) 2138.76 0.96 1.68 31.13 2139.04 0.95 0.85 30.70 2139.45 0.99 0.49 31.17 -- -- -- -- Prog Q(k) x Q(k) - R 2123.59 1.11 1.76 33.59 2123.85 1.10 0.88 32.85 2124.25 1.14 0.51 32.74 -- -- -- -- Prog RK2 Q(k) x Q(k) 2138.45 0.97 3.10 31.01 2138.43 0.97 1.30 31.31 2138.46 1.04 0.70 34.43 2138.49 0.99 0.60 36.46 Prog RK2 Q(k) x Q(k) - R 2123.32 1.12 3.18 33.55 2123.31 1.12 1.32 33.54 2123.29 1.19 0.70 36.03 2123.37 1.14 0.63 37.74 Prog RK4 Q(k) x Q(k) 2138.46 0.97 2.95 31.00 2138.44 0.97 1.46 31.24 2138.51 1.03 0.86 34.07 -- -- -- -- Prog RK4 Q(k) x Q(k) - R 2123.32 1.12 3.03 33.54 2123.31 1.12 1.52 33.48 2123.34 1.18 0.88 35.72 -- -- -- -- ) shows the values of Figure 4.

Figure 3
Region of the hydrograph peak flow in the river's last reach for the Bates, Horritt and Fewtrell (2010)BATES, P. D.; HORRITT, M. S.; FEWTRELL, T. J. A simple inertial formulation of the shallow water equations for efficient two-dimensional flood inundation modelling. Journal of Hydrology, v. 387, n. 1, p. 33-45, 2010. http://dx.doi.org/10.1016/j.jhydrol.2010.03.027.
http://dx.doi.org/10.1016/j.jhydrol.2010...
, Prog Q(k+1)xQ(k+1) and Prog RK2 Q(k)xQ(k) schemes with and without the hydraulic radius simplification considering a river's slope equal to 5 cm/km (the red line is subscribed by the green line).
Figure 4
Peak flow, volume conservation error, processing time and RMSE in relation to the time step for different formulations of numerical schemes considering a river's slope equal to 5 cm/km.

The Prog Q(k+1)xQ(k+1) scheme was the most stable, with stability up to the time step of 3.5 min (C = 1.04), and with an increase in processing time of 20% relative to the original scheme. In addition, it was the scheme that presented the highest RMSE values.

Schemes with centered derivative presented instability at all time steps adopted. The other schemes, such as Bates, Horritt and Fewtrell (2010)BATES, P. D.; HORRITT, M. S.; FEWTRELL, T. J. A simple inertial formulation of the shallow water equations for efficient two-dimensional flood inundation modelling. Journal of Hydrology, v. 387, n. 1, p. 33-45, 2010. http://dx.doi.org/10.1016/j.jhydrol.2010.03.027.
http://dx.doi.org/10.1016/j.jhydrol.2010...
and Almeida et al. (2012)ALMEIDA, G. A. M.; BATES, P.; FREER, J. E.; SOUVIGNET, M. Improving the stability of a simple formulation of the shallow water equations for 2-D flood modeling. Water Resources Research, v. 48, n. 5, p. W05528, 2012. http://dx.doi.org/10.1029/2011WR011570.
http://dx.doi.org/10.1029/2011WR011570...
, presented equal stability up to a time step of 3 min (C = 0.90).

The Prog Q(k)xQ(k) scheme presented the lowest RMSE with practically the same processing time as the original model. The Prog RK2 Q(k)xQ(k) and Prog RK4 Q(k)xQ(k) schemes presented RMSE values closer to the ones in Bates, Horritt and Fewtrell (2010)BATES, P. D.; HORRITT, M. S.; FEWTRELL, T. J. A simple inertial formulation of the shallow water equations for efficient two-dimensional flood inundation modelling. Journal of Hydrology, v. 387, n. 1, p. 33-45, 2010. http://dx.doi.org/10.1016/j.jhydrol.2010.03.027.
http://dx.doi.org/10.1016/j.jhydrol.2010...
, but the processing time was, on average, almost twice.

The volume conservation error did not changed with the time step but increased from 0.57% to 0.75% when the calculation of the hydraulic radius was included in the numerical scheme. This inclusion of this calculation did not increase the processing time, as observed in Figure 4, and led to the results closer to that of the HEC-RAS model, reducing RMSE by 64% compared to the Bates, Horritt and Fewtrell (2010)BATES, P. D.; HORRITT, M. S.; FEWTRELL, T. J. A simple inertial formulation of the shallow water equations for efficient two-dimensional flood inundation modelling. Journal of Hydrology, v. 387, n. 1, p. 33-45, 2010. http://dx.doi.org/10.1016/j.jhydrol.2010.03.027.
http://dx.doi.org/10.1016/j.jhydrol.2010...
scheme.

Test 3

In this test, a high slope of 300 cm/km was considered. The time step was varied between 1 and 4 minutes corresponding to C values of 0.18, 0.37, 0.55, 0.74 at the time of greatest depth. Figure 5 shows a section of the hydrograph at the river's last reach for the original scheme, the Cent Q(k+1)xQ(k-1), Prog Q(k+1)xQ(k+1) schemes and HEC-RAS. Figures 6 shows the main results and Table A3 (Appendix A Appendix A Supplementary data. Table A1 Results of test 1 (real reach) for different numerical schemes. Method dt = 1 min (C = 0.26) dt = 2 min (C = 0.51) dt = 3 min (C = 0.77) dt = 4 min (C =1.03) Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) HEC-RAS 3399.87 -0.02 3387.35 -0.02 3374.96 -0.02 3362.68 -0.02 Almeida et al. (2012) 3444.76 0.00 1.76 49.49 3444.21 0.00 0.86 48.73 3443.95 0.00 0.59 47.93 -- -- -- -- Bates, Horritt and Fewtrell (2010) 3443.79 0.00 1.27 50.59 3443.73 0.00 0.62 49.27 3443.62 0.00 0.41 48.28 -- -- -- -- Bates, Horritt and Fewtrell (2010) - R 3411.96 0.05 1.25 31.54 3411.91 0.05 0.63 29.84 3411.78 0.05 0.42 28.60 -- -- -- -- Cent Q(k+1) x Q(k-1) 3430.70 0.00 1.60 51.42 -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) - R 3399.17 0.05 1.64 33.06 -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) 3457.21 0.00 1.64 49.85 3470.77 0.00 0.81 48.50 3484.48 0.00 0.58 48.24 3498.38 0.00 0.42 49.08 Prog Q(k+1) x Q(k+1) - R 3425.07 0.05 1.66 30.20 3438.31 0.05 0.85 28.26 3451.66 0.05 0.56 28.13 3465.25 0.05 0.42 29.84 Cent Q(k+1) x Q(k+1) 3457.43 0.00 1.69 49.76 -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k+1) - R 3425.27 0.05 1.71 30.12 -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k) x Q(k) 3430.56 0.00 1.60 51.48 3417.47 0.00 0.79 50.79 3404.51 0.00 0.54 50.15 -- -- -- -- Prog Q(k) x Q(k) - R 3399.04 0.05 1.74 33.11 3386.25 0.05 0.87 32.60 3373.57 0.05 0.57 32.13 -- -- -- -- Prog RK2 Q(k) x Q(k) 3445.48 0.00 3.04 50.47 3451.58 0.00 1.31 48.94 3465.15 0.00 0.80 48.27 3496.15 0.00 0.58 53.43 Prog RK2 Q(k) x Q(k) - R 3413.62 0.05 3.09 31.31 3419.67 0.05 1.30 29.08 3433.14 0.05 0.82 27.96 3464.26 0.05 0.58 35.80 Prog RK4 Q(k) x Q(k) 3444.29 0.00 2.96 50.55 3445.80 0.00 1.48 49.17 3448.54 0.00 1.00 48.17 -- -- -- -- Prog RK4 Q(k) x Q(k) - R 3412.46 0.05 3.09 31.47 3413.95 0.05 1.56 29.63 3416.67 0.05 1.02 28.28 -- -- -- -- Table A2 Results of test 2 (low slope) for different numerical schemes. Method dt = 1 min (C = 0.29) dt = 2 min (C = 0.60) dt = 3 min (C = 0.90) dt = 3,5 min (C = 1.04) Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) HEC-RAS 1562.90 0.64 1561.60 0.64 1560.30 0.64 1559.00 0.64 Almeida et al. (2012) 1617.24 0.57 1.73 19.34 1617.63 0.57 0.89 19.52 1617.76 0.57 0.59 19.58 -- -- -- -- Bates, Horritt and Fewtrell (2010) 1618.03 0.57 1.26 18.51 1618.03 0.57 0.62 18.65 1618.02 0.57 0.43 18.81 -- -- -- -- Bates, Horritt and Fewtrell (2010) - R 1598.35 0.75 1.27 10.99 1598.34 0.75 0.64 11.23 1598.34 0.75 0.42 11.49 -- -- -- -- Cent Q(k+1) x Q(k-1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) 1619.57 0.57 1.51 18.69 1621.12 0.57 0.76 19.06 1622.67 0.56 0.51 19.50 1623.44 0.56 0.43 20.39 Prog Q(k+1) x Q(k+1) - R 1599.83 0.75 1.54 11.25 1601.31 0.74 0.77 11.83 1602.80 0.74 0.51 12.51 1603.55 0.74 0.45 13.12 Cent Q(k+1) x Q(k+1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k+1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k) x Q(k) 1616.50 0.57 1.53 18.36 1614.97 0.57 0.78 18.33 1613.45 0.57 0.51 18.29 -- -- -- -- Prog Q(k) x Q(k) - R 1596.87 0.75 1.53 10.76 1595.41 0.75 0.78 10.74 1593.95 0.75 0.51 10.72 -- -- -- -- Prog RK2 Q(k) x Q(k) 1618.10 0.57 2.93 18.52 1618.34 0.57 1.24 18.68 1618.78 0.57 0.76 18.89 -- -- -- -- Prog RK2 Q(k) x Q(k) - R 1598.42 0.75 3.01 11.00 1598.65 0.75 1.29 11.28 1599.07 0.75 0.78 11.64 -- -- -- -- Prog RK4 Q(k) x Q(k) 1618.05 0.57 2.19 18.52 1618.12 0.57 1.11 18.66 1618.24 0.57 0.73 18.83 -- -- -- -- Prog RK4 Q(k) x Q(k) - R 1598.37 0.75 2.27 10.99 1598.44 0.75 1.17 11.24 1598.55 0.75 0.77 11.53 -- -- -- -- Table A3 Results of test 3 (high slope) for different numerical schemes. Method dt = 1 min (C = 0.18) dt = 2 min (C = 0.37) dt = 3 min (C = 0.55) dt = 4 min (C =0.74) Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) HEC-RAS 4905.42 -0.07 4905.42 -0.07 NaN 4905.42 -0.07 4905.42 -0.07 Almeida et al. (2012) 4906.23 0.00 1.76 57.33 4906.02 0.00 0.88 57.69 -- -- -- -- -- -- -- -- Bates, Horritt and Fewtrell (2010) 4905.97 0.00 1.27 57.74 4905.89 0.00 0.64 57.89 -- -- -- -- -- -- -- -- Bates, Horritt and Fewtrell (2010) - R 4903.83 0.01 1.30 47.64 4903.70 0.01 0.64 47.80 -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) 4889.26 0.00 1.55 61.29 4872.40 0.00 0.76 65.21 -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) - R 4886.97 0.01 1.54 52.09 4869.91 0.01 0.77 56.80 -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) 4923.00 0.00 1.55 54.33 4941.11 0.00 0.77 51.38 4963.35 0.00 0.51 49.14 5000.27 0.00 0.38 48.44 Prog Q(k+1) x Q(k+1) - R 4920.96 0.01 1.54 43.26 4939.08 0.01 0.79 39.31 4960.83 0.01 0.51 36.25 4995.77 0.01 0.38 35.34 Cent Q(k+1) x Q(k+1) 4923.14 0.00 1.53 54.24 4941.25 0.00 0.76 51.17 4961.23 0.00 0.51 48.80 -- -- -- -- Cent Q(k+1) x Q(k+1) - R 4921.11 0.01 1.54 43.19 4939.21 0.01 0.78 39.17 4959.08 0.01 0.51 36.02 -- -- -- -- Prog Q(k) x Q(k) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k) x Q(k) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK2 Q(k) x Q(k) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK2 Q(k) x Q(k) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK4 Q(k) x Q(k) 4908.49 0.00 2.19 57.54 -- -- -- -- -- -- -- -- -- -- -- -- Prog RK4 Q(k) x Q(k) - R 4906.39 0.01 2.26 47.12 -- -- -- -- -- -- -- -- -- -- -- -- Table A4 Results of test 4 (very high slope) for different numerical schemes. Method dt = 80 seg (C = 0.20) dt = 120 seg (C = 0.31) dt = 160 seg (C = 0.41) dt = 200 seg (C ~ 0.53) Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) HEC-RAS 4910.30 -0.09 4910.30 -0.09 4910.30 -0.09 4910.30 -0.09 Almeida et al. (2012) 4942.79 0.00 1.20 30.18 -- -- -- -- -- -- -- -- -- -- -- -- Bates, Horritt and Fewtrell (2010) 4942.78 0.00 1.12 30.19 -- -- -- -- -- -- -- -- -- -- -- -- Bates, Horritt and Fewtrell (2010) - R 4942.18 0.01 1.15 25.75 -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) 4963.23 0.00 1.15 31.24 4978.55 0.00 0.77 32.18 5005.51 0.00 0.58 33.48 -- -- -- -- Prog Q(k+1) x Q(k+1) - R 4962.55 0.01 1.19 26.92 4977.66 0.01 0.78 27.97 5003.56 0.01 0.57 29.41 -- -- -- -- Cent Q(k+1) x Q(k+1) 4963.04 0.00 1.14 31.23 4977.23 0.00 0.78 32.16 4997.92 0.00 0.57 33.41 -- -- -- -- Cent Q(k+1) x Q(k+1) - R 4962.37 0.01 1.16 26.91 4976.51 0.01 0.79 27.94 4996.96 0.01 0.58 29.35 -- -- -- -- Prog Q(k) x Q(k) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k) x Q(k) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK2 Q(k) x Q(k) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK2 Q(k) x Q(k) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK4 Q(k) x Q(k) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK4 Q(k) x Q(k) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Table A5 Results of test 5 (backwater) for different numerical schemes. Method dt = 1 min (C = 0.44) dt = 2 min (C = 0.88) dt = 2,2 min (C = 0.97) dt = 2,5 min (C > 1) Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) HEC-RAS 4011.64 0.00 4011.64 0.00 4011.64 0.00 4011.64 0.00 Almeida et al. (2012) 4092.22 0.00 1.82 21.10 4093.60 0.00 0.91 21.41 4093.67 0.00 0.83 21.44 -- -- -- -- Bates, Horritt and Fewtrell (2010) 4095.07 0.00 1.30 21.82 4095.02 0.00 0.65 21.81 4094.96 0.00 0.59 21.81 -- -- -- -- Bates, Horritt and Fewtrell (2010) - R 4051.15 0.00 1.29 15.76 4051.11 0.00 0.65 15.72 4051.06 0.00 0.59 15.71 -- -- -- -- Cent Q(k+1) x Q(k-1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) 4109.87 0.00 1.69 24.04 4124.91 0.00 0.83 26.80 4127.91 0.00 0.70 27.42 -- -- -- -- Prog Q(k+1) x Q(k+1) - R 4065.47 0.00 1.72 18.30 4079.97 0.00 0.86 21.42 4082.89 0.00 0.72 22.10 -- -- -- -- Cent Q(k+1) x Q(k+1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k+1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k) x Q(k) 4080.54 0.00 1.66 20.29 4066.21 0.00 0.84 19.55 4063.34 0.00 0.77 19.50 -- -- -- -- Prog Q(k) x Q(k) - R 4037.06 0.00 1.76 14.09 4023.14 0.00 0.87 13.49 4020.36 0.00 0.79 13.51 -- -- -- -- Prog RK2 Q(k) x Q(k) 4096.34 0.00 3.07 22.00 4100.76 0.00 1.29 22.73 4102.09 0.00 1.13 23.00 -- -- -- -- Prog RK2 Q(k) x Q(k) - R 4052.42 0.00 3.16 15.90 4056.87 0.00 1.34 16.42 4058.23 0.00 1.14 16.61 -- -- -- -- Prog RK4 Q(k) x Q(k) 4095.46 0.00 2.96 21.87 4096.61 0.00 1.48 22.04 4096.89 0.00 1.00 22.10 -- -- -- -- Prog RK4 Q(k) x Q(k) - R 4051.54 0.00 3.07 15.81 4052.69 0.00 1.52 15.91 4052.99 0.00 1.26 15.94 -- -- -- -- Table A6 Results of test 6 (tide) for different numerical schemes. Method dt = 1 min (C = 0.27) dt = 2 min (C = 0.54) dt = 3,5 min (C = 0.94) dt = 4 min (C = 1.07) Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) HEC-RAS 2125.02 -0.84 2125.02 -0.84 2125.02 -0.84 2125.02 -0.84 Almeida et al. (2012) 2141.89 0.94 1.88 32.73 2140.17 0.95 0.90 31.47 2139.53 1.01 0.52 33.63 -- -- -- -- Bates, Horritt and Fewtrell (2010) 2138.46 0.97 1.33 30.99 2138.44 0.97 0.65 31.20 2138.54 1.02 0.37 33.90 -- -- -- -- Bates, Horritt and Fewtrell (2010) - R 2123.32 1.12 1.32 33.53 2123.32 1.12 0.65 33.44 2123.37 1.17 0.37 35.54 -- -- -- -- Cent Q(k+1) x Q(k-1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k+1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k+1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k) x Q(k) 2138.76 0.96 1.68 31.13 2139.04 0.95 0.85 30.70 2139.45 0.99 0.49 31.17 -- -- -- -- Prog Q(k) x Q(k) - R 2123.59 1.11 1.76 33.59 2123.85 1.10 0.88 32.85 2124.25 1.14 0.51 32.74 -- -- -- -- Prog RK2 Q(k) x Q(k) 2138.45 0.97 3.10 31.01 2138.43 0.97 1.30 31.31 2138.46 1.04 0.70 34.43 2138.49 0.99 0.60 36.46 Prog RK2 Q(k) x Q(k) - R 2123.32 1.12 3.18 33.55 2123.31 1.12 1.32 33.54 2123.29 1.19 0.70 36.03 2123.37 1.14 0.63 37.74 Prog RK4 Q(k) x Q(k) 2138.46 0.97 2.95 31.00 2138.44 0.97 1.46 31.24 2138.51 1.03 0.86 34.07 -- -- -- -- Prog RK4 Q(k) x Q(k) - R 2123.32 1.12 3.03 33.54 2123.31 1.12 1.52 33.48 2123.34 1.18 0.88 35.72 -- -- -- -- ) the detailed results.

Figure 5
Region of the hydrograph peak flow in the river's last reach for the Bates, Horritt and Fewtrell (2010)BATES, P. D.; HORRITT, M. S.; FEWTRELL, T. J. A simple inertial formulation of the shallow water equations for efficient two-dimensional flood inundation modelling. Journal of Hydrology, v. 387, n. 1, p. 33-45, 2010. http://dx.doi.org/10.1016/j.jhydrol.2010.03.027.
http://dx.doi.org/10.1016/j.jhydrol.2010...
, Cent Q(k+1)xQ(k-1) and Prog Q(k+1)xQ(k+1) schemes with and without the hydraulic radius simplification considering a river's slope equal to 300 cm/km.
Figure 6
Peak flow, volume conservation error, processing time and RMSE in relation to the time step for different formulations of numerical schemes considering a river's slope equal to 300 cm/km.

As in the previous tests, the Prog Q(k+1)xQ(k+1) scheme was the most stable, reaching the value of C equal to 0.74 compared to 0.37 in the Bates, Horritt and Fewtrell (2010)BATES, P. D.; HORRITT, M. S.; FEWTRELL, T. J. A simple inertial formulation of the shallow water equations for efficient two-dimensional flood inundation modelling. Journal of Hydrology, v. 387, n. 1, p. 33-45, 2010. http://dx.doi.org/10.1016/j.jhydrol.2010.03.027.
http://dx.doi.org/10.1016/j.jhydrol.2010...
and Almeida et al. (2012)ALMEIDA, G. A. M.; BATES, P.; FREER, J. E.; SOUVIGNET, M. Improving the stability of a simple formulation of the shallow water equations for 2-D flood modeling. Water Resources Research, v. 48, n. 5, p. W05528, 2012. http://dx.doi.org/10.1029/2011WR011570.
http://dx.doi.org/10.1029/2011WR011570...
schemes. The increase in processing time compared to the original scheme was 20%.

The Cent Q(k+1)xQ(k+1) scheme also presented better performance than the original model, with stability up to the time step of 3 minutes (C = 0.55). However, the Prog Q(k)xQ(k), Prog RK2 Q(k)xQ(k) and Prog RK4 Q(k)xQ(k) schemes did not show any improvement over the original scheme.

The inclusion of the hydraulic radius calculation in numerical schemes reduced the RMSE by 20% in the original scheme and by 30% in the Prog Q(k+1)xQ(k+1) scheme. Moreover, the increase in the volume conservation error from 0 to 0.01% is negligible considering this calculation.

Test 4

The same river’s reach of test 1 with a very high slope of 10 m/km was considered in this test. This slope changes the flow regime to the supercritical regime, where the original inertial model does not show good stability results (FAN et al., 2014FAN, F. M.; PONTES, P. R. M.; PAIVA, R. C. D.; COLLISCHONN, W. Avaliação de um método de propagação de cheias em rios com aproximação inercial das equações de Saint-Venant. Revista Brasileira de Recursos Hídricos, v. 19, n. 4, p. 137-147, 2014. http://dx.doi.org/10.21168/rbrh.v19n4.p137-147.
http://dx.doi.org/10.21168/rbrh.v19n4.p1...
). The time step was varied with values of 80, 120, 160 and 200 seconds, corresponding to the values of C at the time of greatest depth of 0.20, 0.31, 0.41, and 0.51, respectively. Figures 7 and 8 present the main results and Table A4 (Appendix A Appendix A Supplementary data. Table A1 Results of test 1 (real reach) for different numerical schemes. Method dt = 1 min (C = 0.26) dt = 2 min (C = 0.51) dt = 3 min (C = 0.77) dt = 4 min (C =1.03) Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) HEC-RAS 3399.87 -0.02 3387.35 -0.02 3374.96 -0.02 3362.68 -0.02 Almeida et al. (2012) 3444.76 0.00 1.76 49.49 3444.21 0.00 0.86 48.73 3443.95 0.00 0.59 47.93 -- -- -- -- Bates, Horritt and Fewtrell (2010) 3443.79 0.00 1.27 50.59 3443.73 0.00 0.62 49.27 3443.62 0.00 0.41 48.28 -- -- -- -- Bates, Horritt and Fewtrell (2010) - R 3411.96 0.05 1.25 31.54 3411.91 0.05 0.63 29.84 3411.78 0.05 0.42 28.60 -- -- -- -- Cent Q(k+1) x Q(k-1) 3430.70 0.00 1.60 51.42 -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) - R 3399.17 0.05 1.64 33.06 -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) 3457.21 0.00 1.64 49.85 3470.77 0.00 0.81 48.50 3484.48 0.00 0.58 48.24 3498.38 0.00 0.42 49.08 Prog Q(k+1) x Q(k+1) - R 3425.07 0.05 1.66 30.20 3438.31 0.05 0.85 28.26 3451.66 0.05 0.56 28.13 3465.25 0.05 0.42 29.84 Cent Q(k+1) x Q(k+1) 3457.43 0.00 1.69 49.76 -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k+1) - R 3425.27 0.05 1.71 30.12 -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k) x Q(k) 3430.56 0.00 1.60 51.48 3417.47 0.00 0.79 50.79 3404.51 0.00 0.54 50.15 -- -- -- -- Prog Q(k) x Q(k) - R 3399.04 0.05 1.74 33.11 3386.25 0.05 0.87 32.60 3373.57 0.05 0.57 32.13 -- -- -- -- Prog RK2 Q(k) x Q(k) 3445.48 0.00 3.04 50.47 3451.58 0.00 1.31 48.94 3465.15 0.00 0.80 48.27 3496.15 0.00 0.58 53.43 Prog RK2 Q(k) x Q(k) - R 3413.62 0.05 3.09 31.31 3419.67 0.05 1.30 29.08 3433.14 0.05 0.82 27.96 3464.26 0.05 0.58 35.80 Prog RK4 Q(k) x Q(k) 3444.29 0.00 2.96 50.55 3445.80 0.00 1.48 49.17 3448.54 0.00 1.00 48.17 -- -- -- -- Prog RK4 Q(k) x Q(k) - R 3412.46 0.05 3.09 31.47 3413.95 0.05 1.56 29.63 3416.67 0.05 1.02 28.28 -- -- -- -- Table A2 Results of test 2 (low slope) for different numerical schemes. Method dt = 1 min (C = 0.29) dt = 2 min (C = 0.60) dt = 3 min (C = 0.90) dt = 3,5 min (C = 1.04) Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) HEC-RAS 1562.90 0.64 1561.60 0.64 1560.30 0.64 1559.00 0.64 Almeida et al. (2012) 1617.24 0.57 1.73 19.34 1617.63 0.57 0.89 19.52 1617.76 0.57 0.59 19.58 -- -- -- -- Bates, Horritt and Fewtrell (2010) 1618.03 0.57 1.26 18.51 1618.03 0.57 0.62 18.65 1618.02 0.57 0.43 18.81 -- -- -- -- Bates, Horritt and Fewtrell (2010) - R 1598.35 0.75 1.27 10.99 1598.34 0.75 0.64 11.23 1598.34 0.75 0.42 11.49 -- -- -- -- Cent Q(k+1) x Q(k-1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) 1619.57 0.57 1.51 18.69 1621.12 0.57 0.76 19.06 1622.67 0.56 0.51 19.50 1623.44 0.56 0.43 20.39 Prog Q(k+1) x Q(k+1) - R 1599.83 0.75 1.54 11.25 1601.31 0.74 0.77 11.83 1602.80 0.74 0.51 12.51 1603.55 0.74 0.45 13.12 Cent Q(k+1) x Q(k+1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k+1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k) x Q(k) 1616.50 0.57 1.53 18.36 1614.97 0.57 0.78 18.33 1613.45 0.57 0.51 18.29 -- -- -- -- Prog Q(k) x Q(k) - R 1596.87 0.75 1.53 10.76 1595.41 0.75 0.78 10.74 1593.95 0.75 0.51 10.72 -- -- -- -- Prog RK2 Q(k) x Q(k) 1618.10 0.57 2.93 18.52 1618.34 0.57 1.24 18.68 1618.78 0.57 0.76 18.89 -- -- -- -- Prog RK2 Q(k) x Q(k) - R 1598.42 0.75 3.01 11.00 1598.65 0.75 1.29 11.28 1599.07 0.75 0.78 11.64 -- -- -- -- Prog RK4 Q(k) x Q(k) 1618.05 0.57 2.19 18.52 1618.12 0.57 1.11 18.66 1618.24 0.57 0.73 18.83 -- -- -- -- Prog RK4 Q(k) x Q(k) - R 1598.37 0.75 2.27 10.99 1598.44 0.75 1.17 11.24 1598.55 0.75 0.77 11.53 -- -- -- -- Table A3 Results of test 3 (high slope) for different numerical schemes. Method dt = 1 min (C = 0.18) dt = 2 min (C = 0.37) dt = 3 min (C = 0.55) dt = 4 min (C =0.74) Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) HEC-RAS 4905.42 -0.07 4905.42 -0.07 NaN 4905.42 -0.07 4905.42 -0.07 Almeida et al. (2012) 4906.23 0.00 1.76 57.33 4906.02 0.00 0.88 57.69 -- -- -- -- -- -- -- -- Bates, Horritt and Fewtrell (2010) 4905.97 0.00 1.27 57.74 4905.89 0.00 0.64 57.89 -- -- -- -- -- -- -- -- Bates, Horritt and Fewtrell (2010) - R 4903.83 0.01 1.30 47.64 4903.70 0.01 0.64 47.80 -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) 4889.26 0.00 1.55 61.29 4872.40 0.00 0.76 65.21 -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) - R 4886.97 0.01 1.54 52.09 4869.91 0.01 0.77 56.80 -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) 4923.00 0.00 1.55 54.33 4941.11 0.00 0.77 51.38 4963.35 0.00 0.51 49.14 5000.27 0.00 0.38 48.44 Prog Q(k+1) x Q(k+1) - R 4920.96 0.01 1.54 43.26 4939.08 0.01 0.79 39.31 4960.83 0.01 0.51 36.25 4995.77 0.01 0.38 35.34 Cent Q(k+1) x Q(k+1) 4923.14 0.00 1.53 54.24 4941.25 0.00 0.76 51.17 4961.23 0.00 0.51 48.80 -- -- -- -- Cent Q(k+1) x Q(k+1) - R 4921.11 0.01 1.54 43.19 4939.21 0.01 0.78 39.17 4959.08 0.01 0.51 36.02 -- -- -- -- Prog Q(k) x Q(k) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k) x Q(k) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK2 Q(k) x Q(k) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK2 Q(k) x Q(k) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK4 Q(k) x Q(k) 4908.49 0.00 2.19 57.54 -- -- -- -- -- -- -- -- -- -- -- -- Prog RK4 Q(k) x Q(k) - R 4906.39 0.01 2.26 47.12 -- -- -- -- -- -- -- -- -- -- -- -- Table A4 Results of test 4 (very high slope) for different numerical schemes. Method dt = 80 seg (C = 0.20) dt = 120 seg (C = 0.31) dt = 160 seg (C = 0.41) dt = 200 seg (C ~ 0.53) Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) HEC-RAS 4910.30 -0.09 4910.30 -0.09 4910.30 -0.09 4910.30 -0.09 Almeida et al. (2012) 4942.79 0.00 1.20 30.18 -- -- -- -- -- -- -- -- -- -- -- -- Bates, Horritt and Fewtrell (2010) 4942.78 0.00 1.12 30.19 -- -- -- -- -- -- -- -- -- -- -- -- Bates, Horritt and Fewtrell (2010) - R 4942.18 0.01 1.15 25.75 -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) 4963.23 0.00 1.15 31.24 4978.55 0.00 0.77 32.18 5005.51 0.00 0.58 33.48 -- -- -- -- Prog Q(k+1) x Q(k+1) - R 4962.55 0.01 1.19 26.92 4977.66 0.01 0.78 27.97 5003.56 0.01 0.57 29.41 -- -- -- -- Cent Q(k+1) x Q(k+1) 4963.04 0.00 1.14 31.23 4977.23 0.00 0.78 32.16 4997.92 0.00 0.57 33.41 -- -- -- -- Cent Q(k+1) x Q(k+1) - R 4962.37 0.01 1.16 26.91 4976.51 0.01 0.79 27.94 4996.96 0.01 0.58 29.35 -- -- -- -- Prog Q(k) x Q(k) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k) x Q(k) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK2 Q(k) x Q(k) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK2 Q(k) x Q(k) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK4 Q(k) x Q(k) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK4 Q(k) x Q(k) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Table A5 Results of test 5 (backwater) for different numerical schemes. Method dt = 1 min (C = 0.44) dt = 2 min (C = 0.88) dt = 2,2 min (C = 0.97) dt = 2,5 min (C > 1) Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) HEC-RAS 4011.64 0.00 4011.64 0.00 4011.64 0.00 4011.64 0.00 Almeida et al. (2012) 4092.22 0.00 1.82 21.10 4093.60 0.00 0.91 21.41 4093.67 0.00 0.83 21.44 -- -- -- -- Bates, Horritt and Fewtrell (2010) 4095.07 0.00 1.30 21.82 4095.02 0.00 0.65 21.81 4094.96 0.00 0.59 21.81 -- -- -- -- Bates, Horritt and Fewtrell (2010) - R 4051.15 0.00 1.29 15.76 4051.11 0.00 0.65 15.72 4051.06 0.00 0.59 15.71 -- -- -- -- Cent Q(k+1) x Q(k-1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) 4109.87 0.00 1.69 24.04 4124.91 0.00 0.83 26.80 4127.91 0.00 0.70 27.42 -- -- -- -- Prog Q(k+1) x Q(k+1) - R 4065.47 0.00 1.72 18.30 4079.97 0.00 0.86 21.42 4082.89 0.00 0.72 22.10 -- -- -- -- Cent Q(k+1) x Q(k+1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k+1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k) x Q(k) 4080.54 0.00 1.66 20.29 4066.21 0.00 0.84 19.55 4063.34 0.00 0.77 19.50 -- -- -- -- Prog Q(k) x Q(k) - R 4037.06 0.00 1.76 14.09 4023.14 0.00 0.87 13.49 4020.36 0.00 0.79 13.51 -- -- -- -- Prog RK2 Q(k) x Q(k) 4096.34 0.00 3.07 22.00 4100.76 0.00 1.29 22.73 4102.09 0.00 1.13 23.00 -- -- -- -- Prog RK2 Q(k) x Q(k) - R 4052.42 0.00 3.16 15.90 4056.87 0.00 1.34 16.42 4058.23 0.00 1.14 16.61 -- -- -- -- Prog RK4 Q(k) x Q(k) 4095.46 0.00 2.96 21.87 4096.61 0.00 1.48 22.04 4096.89 0.00 1.00 22.10 -- -- -- -- Prog RK4 Q(k) x Q(k) - R 4051.54 0.00 3.07 15.81 4052.69 0.00 1.52 15.91 4052.99 0.00 1.26 15.94 -- -- -- -- Table A6 Results of test 6 (tide) for different numerical schemes. Method dt = 1 min (C = 0.27) dt = 2 min (C = 0.54) dt = 3,5 min (C = 0.94) dt = 4 min (C = 1.07) Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) HEC-RAS 2125.02 -0.84 2125.02 -0.84 2125.02 -0.84 2125.02 -0.84 Almeida et al. (2012) 2141.89 0.94 1.88 32.73 2140.17 0.95 0.90 31.47 2139.53 1.01 0.52 33.63 -- -- -- -- Bates, Horritt and Fewtrell (2010) 2138.46 0.97 1.33 30.99 2138.44 0.97 0.65 31.20 2138.54 1.02 0.37 33.90 -- -- -- -- Bates, Horritt and Fewtrell (2010) - R 2123.32 1.12 1.32 33.53 2123.32 1.12 0.65 33.44 2123.37 1.17 0.37 35.54 -- -- -- -- Cent Q(k+1) x Q(k-1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k+1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k+1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k) x Q(k) 2138.76 0.96 1.68 31.13 2139.04 0.95 0.85 30.70 2139.45 0.99 0.49 31.17 -- -- -- -- Prog Q(k) x Q(k) - R 2123.59 1.11 1.76 33.59 2123.85 1.10 0.88 32.85 2124.25 1.14 0.51 32.74 -- -- -- -- Prog RK2 Q(k) x Q(k) 2138.45 0.97 3.10 31.01 2138.43 0.97 1.30 31.31 2138.46 1.04 0.70 34.43 2138.49 0.99 0.60 36.46 Prog RK2 Q(k) x Q(k) - R 2123.32 1.12 3.18 33.55 2123.31 1.12 1.32 33.54 2123.29 1.19 0.70 36.03 2123.37 1.14 0.63 37.74 Prog RK4 Q(k) x Q(k) 2138.46 0.97 2.95 31.00 2138.44 0.97 1.46 31.24 2138.51 1.03 0.86 34.07 -- -- -- -- Prog RK4 Q(k) x Q(k) - R 2123.32 1.12 3.03 33.54 2123.31 1.12 1.52 33.48 2123.34 1.18 0.88 35.72 -- -- -- -- ) the detailed results.

Figure 7
Region of the hydrograph peak flow in the river's last reach for the Bates, Horritt and Fewtrell (2010)BATES, P. D.; HORRITT, M. S.; FEWTRELL, T. J. A simple inertial formulation of the shallow water equations for efficient two-dimensional flood inundation modelling. Journal of Hydrology, v. 387, n. 1, p. 33-45, 2010. http://dx.doi.org/10.1016/j.jhydrol.2010.03.027.
http://dx.doi.org/10.1016/j.jhydrol.2010...
and Prog Q(k+1)xQ(k+1) schemes with and without the hydraulic radius simplification considering a river's slope equal to 10 cm/km.
Figure 8
Peak flow, volume conservation error, processing time and RMSE in relation to the time step for different formulations of numerical schemes considering a river's slope equal to 10 m/km.

Only the numerical schemes Prog Q(k+1)xQ(k+1), Cent Q(k+1)xQ(k+1),Bates, Horritt and Fewtrell (2010)BATES, P. D.; HORRITT, M. S.; FEWTRELL, T. J. A simple inertial formulation of the shallow water equations for efficient two-dimensional flood inundation modelling. Journal of Hydrology, v. 387, n. 1, p. 33-45, 2010. http://dx.doi.org/10.1016/j.jhydrol.2010.03.027.
http://dx.doi.org/10.1016/j.jhydrol.2010...
and Almeida et al. (2012)ALMEIDA, G. A. M.; BATES, P.; FREER, J. E.; SOUVIGNET, M. Improving the stability of a simple formulation of the shallow water equations for 2-D flood modeling. Water Resources Research, v. 48, n. 5, p. W05528, 2012. http://dx.doi.org/10.1029/2011WR011570.
http://dx.doi.org/10.1029/2011WR011570...
presented stable results. The last two were stable only in the time step adopted in the simulation of 80 sec, corresponding to C = 0.20. The other two schemes were stable up to 160 sec (C = 0.41). Despite the improvement, where the Courant number increased from 0.20 to 0.41, the value is still small compared to other models of flood propagation where C can reach 1.

The Prog Q(k+1)xQ(k+1) and Cent Q(k+1)xQ(k+1) schemes showed practically identical performances with near peak flow, volume conservation error, processing time and RMSE. In addition, the processing time of these schemes was the same as the original scheme with a slightly higher RMSE, 31.24 m3/s compared to 30.19 m3/s.

As in the other tests, the inclusion of the hydraulic radius calculation in the numerical schemes did not increase the processing time, reduced the RMSE from 30.19 to 25.75 in the model of Bates, Horritt and Fewtrell (2010)BATES, P. D.; HORRITT, M. S.; FEWTRELL, T. J. A simple inertial formulation of the shallow water equations for efficient two-dimensional flood inundation modelling. Journal of Hydrology, v. 387, n. 1, p. 33-45, 2010. http://dx.doi.org/10.1016/j.jhydrol.2010.03.027.
http://dx.doi.org/10.1016/j.jhydrol.2010...
while the volume conservation error increased from 0 to 0.01% in all stable numerical schemes, but may be considered negligible.

Test 5

This application was identical to test 1 but considered a constant water level downstream and a constant initial level from the half of the reach as boundary conditions. These conditions represent a dam operating at constant water level in simplified form. Other models of flood propagation, such as Muskingum-Cunge, are invalidated in the presence of dam and backwater (FAN et al., 2014FAN, F. M.; PONTES, P. R. M.; PAIVA, R. C. D.; COLLISCHONN, W. Avaliação de um método de propagação de cheias em rios com aproximação inercial das equações de Saint-Venant. Revista Brasileira de Recursos Hídricos, v. 19, n. 4, p. 137-147, 2014. http://dx.doi.org/10.21168/rbrh.v19n4.p137-147.
http://dx.doi.org/10.21168/rbrh.v19n4.p1...
). The time step was varied with values of 1, 2, 2.2 and 2.5 minutes, corresponding to the maximum C values of 0.44, 0.88, 0.97 and > 1 at the moment of greatest depth (higher values computed during the simulation). Figure 9 shows a hydrograph of the last reach for the original, Prog Q(k+1)xQ(k+1) and Prog RK2 Q(k)xQ(k) schemes. Figure 10 shows the main results and Table A5 (Appendix A Appendix A Supplementary data. Table A1 Results of test 1 (real reach) for different numerical schemes. Method dt = 1 min (C = 0.26) dt = 2 min (C = 0.51) dt = 3 min (C = 0.77) dt = 4 min (C =1.03) Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) HEC-RAS 3399.87 -0.02 3387.35 -0.02 3374.96 -0.02 3362.68 -0.02 Almeida et al. (2012) 3444.76 0.00 1.76 49.49 3444.21 0.00 0.86 48.73 3443.95 0.00 0.59 47.93 -- -- -- -- Bates, Horritt and Fewtrell (2010) 3443.79 0.00 1.27 50.59 3443.73 0.00 0.62 49.27 3443.62 0.00 0.41 48.28 -- -- -- -- Bates, Horritt and Fewtrell (2010) - R 3411.96 0.05 1.25 31.54 3411.91 0.05 0.63 29.84 3411.78 0.05 0.42 28.60 -- -- -- -- Cent Q(k+1) x Q(k-1) 3430.70 0.00 1.60 51.42 -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) - R 3399.17 0.05 1.64 33.06 -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) 3457.21 0.00 1.64 49.85 3470.77 0.00 0.81 48.50 3484.48 0.00 0.58 48.24 3498.38 0.00 0.42 49.08 Prog Q(k+1) x Q(k+1) - R 3425.07 0.05 1.66 30.20 3438.31 0.05 0.85 28.26 3451.66 0.05 0.56 28.13 3465.25 0.05 0.42 29.84 Cent Q(k+1) x Q(k+1) 3457.43 0.00 1.69 49.76 -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k+1) - R 3425.27 0.05 1.71 30.12 -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k) x Q(k) 3430.56 0.00 1.60 51.48 3417.47 0.00 0.79 50.79 3404.51 0.00 0.54 50.15 -- -- -- -- Prog Q(k) x Q(k) - R 3399.04 0.05 1.74 33.11 3386.25 0.05 0.87 32.60 3373.57 0.05 0.57 32.13 -- -- -- -- Prog RK2 Q(k) x Q(k) 3445.48 0.00 3.04 50.47 3451.58 0.00 1.31 48.94 3465.15 0.00 0.80 48.27 3496.15 0.00 0.58 53.43 Prog RK2 Q(k) x Q(k) - R 3413.62 0.05 3.09 31.31 3419.67 0.05 1.30 29.08 3433.14 0.05 0.82 27.96 3464.26 0.05 0.58 35.80 Prog RK4 Q(k) x Q(k) 3444.29 0.00 2.96 50.55 3445.80 0.00 1.48 49.17 3448.54 0.00 1.00 48.17 -- -- -- -- Prog RK4 Q(k) x Q(k) - R 3412.46 0.05 3.09 31.47 3413.95 0.05 1.56 29.63 3416.67 0.05 1.02 28.28 -- -- -- -- Table A2 Results of test 2 (low slope) for different numerical schemes. Method dt = 1 min (C = 0.29) dt = 2 min (C = 0.60) dt = 3 min (C = 0.90) dt = 3,5 min (C = 1.04) Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) HEC-RAS 1562.90 0.64 1561.60 0.64 1560.30 0.64 1559.00 0.64 Almeida et al. (2012) 1617.24 0.57 1.73 19.34 1617.63 0.57 0.89 19.52 1617.76 0.57 0.59 19.58 -- -- -- -- Bates, Horritt and Fewtrell (2010) 1618.03 0.57 1.26 18.51 1618.03 0.57 0.62 18.65 1618.02 0.57 0.43 18.81 -- -- -- -- Bates, Horritt and Fewtrell (2010) - R 1598.35 0.75 1.27 10.99 1598.34 0.75 0.64 11.23 1598.34 0.75 0.42 11.49 -- -- -- -- Cent Q(k+1) x Q(k-1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) 1619.57 0.57 1.51 18.69 1621.12 0.57 0.76 19.06 1622.67 0.56 0.51 19.50 1623.44 0.56 0.43 20.39 Prog Q(k+1) x Q(k+1) - R 1599.83 0.75 1.54 11.25 1601.31 0.74 0.77 11.83 1602.80 0.74 0.51 12.51 1603.55 0.74 0.45 13.12 Cent Q(k+1) x Q(k+1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k+1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k) x Q(k) 1616.50 0.57 1.53 18.36 1614.97 0.57 0.78 18.33 1613.45 0.57 0.51 18.29 -- -- -- -- Prog Q(k) x Q(k) - R 1596.87 0.75 1.53 10.76 1595.41 0.75 0.78 10.74 1593.95 0.75 0.51 10.72 -- -- -- -- Prog RK2 Q(k) x Q(k) 1618.10 0.57 2.93 18.52 1618.34 0.57 1.24 18.68 1618.78 0.57 0.76 18.89 -- -- -- -- Prog RK2 Q(k) x Q(k) - R 1598.42 0.75 3.01 11.00 1598.65 0.75 1.29 11.28 1599.07 0.75 0.78 11.64 -- -- -- -- Prog RK4 Q(k) x Q(k) 1618.05 0.57 2.19 18.52 1618.12 0.57 1.11 18.66 1618.24 0.57 0.73 18.83 -- -- -- -- Prog RK4 Q(k) x Q(k) - R 1598.37 0.75 2.27 10.99 1598.44 0.75 1.17 11.24 1598.55 0.75 0.77 11.53 -- -- -- -- Table A3 Results of test 3 (high slope) for different numerical schemes. Method dt = 1 min (C = 0.18) dt = 2 min (C = 0.37) dt = 3 min (C = 0.55) dt = 4 min (C =0.74) Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) HEC-RAS 4905.42 -0.07 4905.42 -0.07 NaN 4905.42 -0.07 4905.42 -0.07 Almeida et al. (2012) 4906.23 0.00 1.76 57.33 4906.02 0.00 0.88 57.69 -- -- -- -- -- -- -- -- Bates, Horritt and Fewtrell (2010) 4905.97 0.00 1.27 57.74 4905.89 0.00 0.64 57.89 -- -- -- -- -- -- -- -- Bates, Horritt and Fewtrell (2010) - R 4903.83 0.01 1.30 47.64 4903.70 0.01 0.64 47.80 -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) 4889.26 0.00 1.55 61.29 4872.40 0.00 0.76 65.21 -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) - R 4886.97 0.01 1.54 52.09 4869.91 0.01 0.77 56.80 -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) 4923.00 0.00 1.55 54.33 4941.11 0.00 0.77 51.38 4963.35 0.00 0.51 49.14 5000.27 0.00 0.38 48.44 Prog Q(k+1) x Q(k+1) - R 4920.96 0.01 1.54 43.26 4939.08 0.01 0.79 39.31 4960.83 0.01 0.51 36.25 4995.77 0.01 0.38 35.34 Cent Q(k+1) x Q(k+1) 4923.14 0.00 1.53 54.24 4941.25 0.00 0.76 51.17 4961.23 0.00 0.51 48.80 -- -- -- -- Cent Q(k+1) x Q(k+1) - R 4921.11 0.01 1.54 43.19 4939.21 0.01 0.78 39.17 4959.08 0.01 0.51 36.02 -- -- -- -- Prog Q(k) x Q(k) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k) x Q(k) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK2 Q(k) x Q(k) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK2 Q(k) x Q(k) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK4 Q(k) x Q(k) 4908.49 0.00 2.19 57.54 -- -- -- -- -- -- -- -- -- -- -- -- Prog RK4 Q(k) x Q(k) - R 4906.39 0.01 2.26 47.12 -- -- -- -- -- -- -- -- -- -- -- -- Table A4 Results of test 4 (very high slope) for different numerical schemes. Method dt = 80 seg (C = 0.20) dt = 120 seg (C = 0.31) dt = 160 seg (C = 0.41) dt = 200 seg (C ~ 0.53) Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) HEC-RAS 4910.30 -0.09 4910.30 -0.09 4910.30 -0.09 4910.30 -0.09 Almeida et al. (2012) 4942.79 0.00 1.20 30.18 -- -- -- -- -- -- -- -- -- -- -- -- Bates, Horritt and Fewtrell (2010) 4942.78 0.00 1.12 30.19 -- -- -- -- -- -- -- -- -- -- -- -- Bates, Horritt and Fewtrell (2010) - R 4942.18 0.01 1.15 25.75 -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) 4963.23 0.00 1.15 31.24 4978.55 0.00 0.77 32.18 5005.51 0.00 0.58 33.48 -- -- -- -- Prog Q(k+1) x Q(k+1) - R 4962.55 0.01 1.19 26.92 4977.66 0.01 0.78 27.97 5003.56 0.01 0.57 29.41 -- -- -- -- Cent Q(k+1) x Q(k+1) 4963.04 0.00 1.14 31.23 4977.23 0.00 0.78 32.16 4997.92 0.00 0.57 33.41 -- -- -- -- Cent Q(k+1) x Q(k+1) - R 4962.37 0.01 1.16 26.91 4976.51 0.01 0.79 27.94 4996.96 0.01 0.58 29.35 -- -- -- -- Prog Q(k) x Q(k) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k) x Q(k) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK2 Q(k) x Q(k) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK2 Q(k) x Q(k) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK4 Q(k) x Q(k) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK4 Q(k) x Q(k) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Table A5 Results of test 5 (backwater) for different numerical schemes. Method dt = 1 min (C = 0.44) dt = 2 min (C = 0.88) dt = 2,2 min (C = 0.97) dt = 2,5 min (C > 1) Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) HEC-RAS 4011.64 0.00 4011.64 0.00 4011.64 0.00 4011.64 0.00 Almeida et al. (2012) 4092.22 0.00 1.82 21.10 4093.60 0.00 0.91 21.41 4093.67 0.00 0.83 21.44 -- -- -- -- Bates, Horritt and Fewtrell (2010) 4095.07 0.00 1.30 21.82 4095.02 0.00 0.65 21.81 4094.96 0.00 0.59 21.81 -- -- -- -- Bates, Horritt and Fewtrell (2010) - R 4051.15 0.00 1.29 15.76 4051.11 0.00 0.65 15.72 4051.06 0.00 0.59 15.71 -- -- -- -- Cent Q(k+1) x Q(k-1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) 4109.87 0.00 1.69 24.04 4124.91 0.00 0.83 26.80 4127.91 0.00 0.70 27.42 -- -- -- -- Prog Q(k+1) x Q(k+1) - R 4065.47 0.00 1.72 18.30 4079.97 0.00 0.86 21.42 4082.89 0.00 0.72 22.10 -- -- -- -- Cent Q(k+1) x Q(k+1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k+1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k) x Q(k) 4080.54 0.00 1.66 20.29 4066.21 0.00 0.84 19.55 4063.34 0.00 0.77 19.50 -- -- -- -- Prog Q(k) x Q(k) - R 4037.06 0.00 1.76 14.09 4023.14 0.00 0.87 13.49 4020.36 0.00 0.79 13.51 -- -- -- -- Prog RK2 Q(k) x Q(k) 4096.34 0.00 3.07 22.00 4100.76 0.00 1.29 22.73 4102.09 0.00 1.13 23.00 -- -- -- -- Prog RK2 Q(k) x Q(k) - R 4052.42 0.00 3.16 15.90 4056.87 0.00 1.34 16.42 4058.23 0.00 1.14 16.61 -- -- -- -- Prog RK4 Q(k) x Q(k) 4095.46 0.00 2.96 21.87 4096.61 0.00 1.48 22.04 4096.89 0.00 1.00 22.10 -- -- -- -- Prog RK4 Q(k) x Q(k) - R 4051.54 0.00 3.07 15.81 4052.69 0.00 1.52 15.91 4052.99 0.00 1.26 15.94 -- -- -- -- Table A6 Results of test 6 (tide) for different numerical schemes. Method dt = 1 min (C = 0.27) dt = 2 min (C = 0.54) dt = 3,5 min (C = 0.94) dt = 4 min (C = 1.07) Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) HEC-RAS 2125.02 -0.84 2125.02 -0.84 2125.02 -0.84 2125.02 -0.84 Almeida et al. (2012) 2141.89 0.94 1.88 32.73 2140.17 0.95 0.90 31.47 2139.53 1.01 0.52 33.63 -- -- -- -- Bates, Horritt and Fewtrell (2010) 2138.46 0.97 1.33 30.99 2138.44 0.97 0.65 31.20 2138.54 1.02 0.37 33.90 -- -- -- -- Bates, Horritt and Fewtrell (2010) - R 2123.32 1.12 1.32 33.53 2123.32 1.12 0.65 33.44 2123.37 1.17 0.37 35.54 -- -- -- -- Cent Q(k+1) x Q(k-1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k+1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k+1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k) x Q(k) 2138.76 0.96 1.68 31.13 2139.04 0.95 0.85 30.70 2139.45 0.99 0.49 31.17 -- -- -- -- Prog Q(k) x Q(k) - R 2123.59 1.11 1.76 33.59 2123.85 1.10 0.88 32.85 2124.25 1.14 0.51 32.74 -- -- -- -- Prog RK2 Q(k) x Q(k) 2138.45 0.97 3.10 31.01 2138.43 0.97 1.30 31.31 2138.46 1.04 0.70 34.43 2138.49 0.99 0.60 36.46 Prog RK2 Q(k) x Q(k) - R 2123.32 1.12 3.18 33.55 2123.31 1.12 1.32 33.54 2123.29 1.19 0.70 36.03 2123.37 1.14 0.63 37.74 Prog RK4 Q(k) x Q(k) 2138.46 0.97 2.95 31.00 2138.44 0.97 1.46 31.24 2138.51 1.03 0.86 34.07 -- -- -- -- Prog RK4 Q(k) x Q(k) - R 2123.32 1.12 3.03 33.54 2123.31 1.12 1.52 33.48 2123.34 1.18 0.88 35.72 -- -- -- -- ) the detailed results.

Figure 9
Region of the hydrograph peak flow in the river's last reach for the Bates, Horritt and Fewtrell (2010)BATES, P. D.; HORRITT, M. S.; FEWTRELL, T. J. A simple inertial formulation of the shallow water equations for efficient two-dimensional flood inundation modelling. Journal of Hydrology, v. 387, n. 1, p. 33-45, 2010. http://dx.doi.org/10.1016/j.jhydrol.2010.03.027.
http://dx.doi.org/10.1016/j.jhydrol.2010...
, Prog Q(k+1)xQ(k+1) and Prog RK2 Q(k)xQ(k) schemes with and without the hydraulic radius simplification considering a boundary condition downstream of dam and backwater (the red line is subscribed by the green line).
Figure 10
Peak flow, volume conservation error, processing time and RMSE in relation to the time step for different formulations of numerical schemes considering a boundary condition downstream of dam and backwater.

All numerical schemes, with the exception of the centered derivative schemes, showed stable results until the adopted time step of 132 sec (C = 0.97). The processing times of the Prog Q(k+1)xQ(k+1) and Prog Q(k)xQ(k) schemes presented an average increase of 25% in comparison to the original scheme. On the other hand, the Prog RK2 Q(k)xQ(k) and Prog RK4 Q(k)xQ(k) schemes took twice as long as the original scheme.

In this application, the volume conservation error did not change with the inclusion of the hydraulic radius calculation. In addition, the peak flow approached the HEC-RAS model's peak flow and the RMSE decreases.

Test 6

In this test, we considered a situation equal to test 2, low slope, but with different contour conditions. A constant flow of 500 m3/s was adopted as the upstream condition and a variable level as the downstream condition. The water level downstream (y) was given by the equation y= 6+2sen(T π6), which represents a periodic level variation between 4m and 8m with a period of 12 hours. This downstream condition - tidal variation - causes a wave that propagates upstream and which can reverse the flow direction in the river’s last reach during some moments along the tide cycle.

The time step was varied between 1, 2, 3.5 and 4 minutes with maximum values of C equal to 0.27, 0.54, 0.94 and 1.07 (highest values computed during the simulation), respectively. Figure 11 shows a region of the hydrograph peak in the last reach for the original scheme and for the Prog RK2 Q(k)xQ(k) scheme. Figure 12 shows the main results and Table A6 (Appendix A Appendix A Supplementary data. Table A1 Results of test 1 (real reach) for different numerical schemes. Method dt = 1 min (C = 0.26) dt = 2 min (C = 0.51) dt = 3 min (C = 0.77) dt = 4 min (C =1.03) Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) HEC-RAS 3399.87 -0.02 3387.35 -0.02 3374.96 -0.02 3362.68 -0.02 Almeida et al. (2012) 3444.76 0.00 1.76 49.49 3444.21 0.00 0.86 48.73 3443.95 0.00 0.59 47.93 -- -- -- -- Bates, Horritt and Fewtrell (2010) 3443.79 0.00 1.27 50.59 3443.73 0.00 0.62 49.27 3443.62 0.00 0.41 48.28 -- -- -- -- Bates, Horritt and Fewtrell (2010) - R 3411.96 0.05 1.25 31.54 3411.91 0.05 0.63 29.84 3411.78 0.05 0.42 28.60 -- -- -- -- Cent Q(k+1) x Q(k-1) 3430.70 0.00 1.60 51.42 -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) - R 3399.17 0.05 1.64 33.06 -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) 3457.21 0.00 1.64 49.85 3470.77 0.00 0.81 48.50 3484.48 0.00 0.58 48.24 3498.38 0.00 0.42 49.08 Prog Q(k+1) x Q(k+1) - R 3425.07 0.05 1.66 30.20 3438.31 0.05 0.85 28.26 3451.66 0.05 0.56 28.13 3465.25 0.05 0.42 29.84 Cent Q(k+1) x Q(k+1) 3457.43 0.00 1.69 49.76 -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k+1) - R 3425.27 0.05 1.71 30.12 -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k) x Q(k) 3430.56 0.00 1.60 51.48 3417.47 0.00 0.79 50.79 3404.51 0.00 0.54 50.15 -- -- -- -- Prog Q(k) x Q(k) - R 3399.04 0.05 1.74 33.11 3386.25 0.05 0.87 32.60 3373.57 0.05 0.57 32.13 -- -- -- -- Prog RK2 Q(k) x Q(k) 3445.48 0.00 3.04 50.47 3451.58 0.00 1.31 48.94 3465.15 0.00 0.80 48.27 3496.15 0.00 0.58 53.43 Prog RK2 Q(k) x Q(k) - R 3413.62 0.05 3.09 31.31 3419.67 0.05 1.30 29.08 3433.14 0.05 0.82 27.96 3464.26 0.05 0.58 35.80 Prog RK4 Q(k) x Q(k) 3444.29 0.00 2.96 50.55 3445.80 0.00 1.48 49.17 3448.54 0.00 1.00 48.17 -- -- -- -- Prog RK4 Q(k) x Q(k) - R 3412.46 0.05 3.09 31.47 3413.95 0.05 1.56 29.63 3416.67 0.05 1.02 28.28 -- -- -- -- Table A2 Results of test 2 (low slope) for different numerical schemes. Method dt = 1 min (C = 0.29) dt = 2 min (C = 0.60) dt = 3 min (C = 0.90) dt = 3,5 min (C = 1.04) Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) HEC-RAS 1562.90 0.64 1561.60 0.64 1560.30 0.64 1559.00 0.64 Almeida et al. (2012) 1617.24 0.57 1.73 19.34 1617.63 0.57 0.89 19.52 1617.76 0.57 0.59 19.58 -- -- -- -- Bates, Horritt and Fewtrell (2010) 1618.03 0.57 1.26 18.51 1618.03 0.57 0.62 18.65 1618.02 0.57 0.43 18.81 -- -- -- -- Bates, Horritt and Fewtrell (2010) - R 1598.35 0.75 1.27 10.99 1598.34 0.75 0.64 11.23 1598.34 0.75 0.42 11.49 -- -- -- -- Cent Q(k+1) x Q(k-1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) 1619.57 0.57 1.51 18.69 1621.12 0.57 0.76 19.06 1622.67 0.56 0.51 19.50 1623.44 0.56 0.43 20.39 Prog Q(k+1) x Q(k+1) - R 1599.83 0.75 1.54 11.25 1601.31 0.74 0.77 11.83 1602.80 0.74 0.51 12.51 1603.55 0.74 0.45 13.12 Cent Q(k+1) x Q(k+1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k+1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k) x Q(k) 1616.50 0.57 1.53 18.36 1614.97 0.57 0.78 18.33 1613.45 0.57 0.51 18.29 -- -- -- -- Prog Q(k) x Q(k) - R 1596.87 0.75 1.53 10.76 1595.41 0.75 0.78 10.74 1593.95 0.75 0.51 10.72 -- -- -- -- Prog RK2 Q(k) x Q(k) 1618.10 0.57 2.93 18.52 1618.34 0.57 1.24 18.68 1618.78 0.57 0.76 18.89 -- -- -- -- Prog RK2 Q(k) x Q(k) - R 1598.42 0.75 3.01 11.00 1598.65 0.75 1.29 11.28 1599.07 0.75 0.78 11.64 -- -- -- -- Prog RK4 Q(k) x Q(k) 1618.05 0.57 2.19 18.52 1618.12 0.57 1.11 18.66 1618.24 0.57 0.73 18.83 -- -- -- -- Prog RK4 Q(k) x Q(k) - R 1598.37 0.75 2.27 10.99 1598.44 0.75 1.17 11.24 1598.55 0.75 0.77 11.53 -- -- -- -- Table A3 Results of test 3 (high slope) for different numerical schemes. Method dt = 1 min (C = 0.18) dt = 2 min (C = 0.37) dt = 3 min (C = 0.55) dt = 4 min (C =0.74) Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) HEC-RAS 4905.42 -0.07 4905.42 -0.07 NaN 4905.42 -0.07 4905.42 -0.07 Almeida et al. (2012) 4906.23 0.00 1.76 57.33 4906.02 0.00 0.88 57.69 -- -- -- -- -- -- -- -- Bates, Horritt and Fewtrell (2010) 4905.97 0.00 1.27 57.74 4905.89 0.00 0.64 57.89 -- -- -- -- -- -- -- -- Bates, Horritt and Fewtrell (2010) - R 4903.83 0.01 1.30 47.64 4903.70 0.01 0.64 47.80 -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) 4889.26 0.00 1.55 61.29 4872.40 0.00 0.76 65.21 -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) - R 4886.97 0.01 1.54 52.09 4869.91 0.01 0.77 56.80 -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) 4923.00 0.00 1.55 54.33 4941.11 0.00 0.77 51.38 4963.35 0.00 0.51 49.14 5000.27 0.00 0.38 48.44 Prog Q(k+1) x Q(k+1) - R 4920.96 0.01 1.54 43.26 4939.08 0.01 0.79 39.31 4960.83 0.01 0.51 36.25 4995.77 0.01 0.38 35.34 Cent Q(k+1) x Q(k+1) 4923.14 0.00 1.53 54.24 4941.25 0.00 0.76 51.17 4961.23 0.00 0.51 48.80 -- -- -- -- Cent Q(k+1) x Q(k+1) - R 4921.11 0.01 1.54 43.19 4939.21 0.01 0.78 39.17 4959.08 0.01 0.51 36.02 -- -- -- -- Prog Q(k) x Q(k) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k) x Q(k) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK2 Q(k) x Q(k) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK2 Q(k) x Q(k) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK4 Q(k) x Q(k) 4908.49 0.00 2.19 57.54 -- -- -- -- -- -- -- -- -- -- -- -- Prog RK4 Q(k) x Q(k) - R 4906.39 0.01 2.26 47.12 -- -- -- -- -- -- -- -- -- -- -- -- Table A4 Results of test 4 (very high slope) for different numerical schemes. Method dt = 80 seg (C = 0.20) dt = 120 seg (C = 0.31) dt = 160 seg (C = 0.41) dt = 200 seg (C ~ 0.53) Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) HEC-RAS 4910.30 -0.09 4910.30 -0.09 4910.30 -0.09 4910.30 -0.09 Almeida et al. (2012) 4942.79 0.00 1.20 30.18 -- -- -- -- -- -- -- -- -- -- -- -- Bates, Horritt and Fewtrell (2010) 4942.78 0.00 1.12 30.19 -- -- -- -- -- -- -- -- -- -- -- -- Bates, Horritt and Fewtrell (2010) - R 4942.18 0.01 1.15 25.75 -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) 4963.23 0.00 1.15 31.24 4978.55 0.00 0.77 32.18 5005.51 0.00 0.58 33.48 -- -- -- -- Prog Q(k+1) x Q(k+1) - R 4962.55 0.01 1.19 26.92 4977.66 0.01 0.78 27.97 5003.56 0.01 0.57 29.41 -- -- -- -- Cent Q(k+1) x Q(k+1) 4963.04 0.00 1.14 31.23 4977.23 0.00 0.78 32.16 4997.92 0.00 0.57 33.41 -- -- -- -- Cent Q(k+1) x Q(k+1) - R 4962.37 0.01 1.16 26.91 4976.51 0.01 0.79 27.94 4996.96 0.01 0.58 29.35 -- -- -- -- Prog Q(k) x Q(k) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k) x Q(k) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK2 Q(k) x Q(k) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK2 Q(k) x Q(k) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK4 Q(k) x Q(k) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog RK4 Q(k) x Q(k) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Table A5 Results of test 5 (backwater) for different numerical schemes. Method dt = 1 min (C = 0.44) dt = 2 min (C = 0.88) dt = 2,2 min (C = 0.97) dt = 2,5 min (C > 1) Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) HEC-RAS 4011.64 0.00 4011.64 0.00 4011.64 0.00 4011.64 0.00 Almeida et al. (2012) 4092.22 0.00 1.82 21.10 4093.60 0.00 0.91 21.41 4093.67 0.00 0.83 21.44 -- -- -- -- Bates, Horritt and Fewtrell (2010) 4095.07 0.00 1.30 21.82 4095.02 0.00 0.65 21.81 4094.96 0.00 0.59 21.81 -- -- -- -- Bates, Horritt and Fewtrell (2010) - R 4051.15 0.00 1.29 15.76 4051.11 0.00 0.65 15.72 4051.06 0.00 0.59 15.71 -- -- -- -- Cent Q(k+1) x Q(k-1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) 4109.87 0.00 1.69 24.04 4124.91 0.00 0.83 26.80 4127.91 0.00 0.70 27.42 -- -- -- -- Prog Q(k+1) x Q(k+1) - R 4065.47 0.00 1.72 18.30 4079.97 0.00 0.86 21.42 4082.89 0.00 0.72 22.10 -- -- -- -- Cent Q(k+1) x Q(k+1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k+1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k) x Q(k) 4080.54 0.00 1.66 20.29 4066.21 0.00 0.84 19.55 4063.34 0.00 0.77 19.50 -- -- -- -- Prog Q(k) x Q(k) - R 4037.06 0.00 1.76 14.09 4023.14 0.00 0.87 13.49 4020.36 0.00 0.79 13.51 -- -- -- -- Prog RK2 Q(k) x Q(k) 4096.34 0.00 3.07 22.00 4100.76 0.00 1.29 22.73 4102.09 0.00 1.13 23.00 -- -- -- -- Prog RK2 Q(k) x Q(k) - R 4052.42 0.00 3.16 15.90 4056.87 0.00 1.34 16.42 4058.23 0.00 1.14 16.61 -- -- -- -- Prog RK4 Q(k) x Q(k) 4095.46 0.00 2.96 21.87 4096.61 0.00 1.48 22.04 4096.89 0.00 1.00 22.10 -- -- -- -- Prog RK4 Q(k) x Q(k) - R 4051.54 0.00 3.07 15.81 4052.69 0.00 1.52 15.91 4052.99 0.00 1.26 15.94 -- -- -- -- Table A6 Results of test 6 (tide) for different numerical schemes. Method dt = 1 min (C = 0.27) dt = 2 min (C = 0.54) dt = 3,5 min (C = 0.94) dt = 4 min (C = 1.07) Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE Qpeak VolErr Time RMSE (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) (m3/s) (%) (sec) (m3/s) HEC-RAS 2125.02 -0.84 2125.02 -0.84 2125.02 -0.84 2125.02 -0.84 Almeida et al. (2012) 2141.89 0.94 1.88 32.73 2140.17 0.95 0.90 31.47 2139.53 1.01 0.52 33.63 -- -- -- -- Bates, Horritt and Fewtrell (2010) 2138.46 0.97 1.33 30.99 2138.44 0.97 0.65 31.20 2138.54 1.02 0.37 33.90 -- -- -- -- Bates, Horritt and Fewtrell (2010) - R 2123.32 1.12 1.32 33.53 2123.32 1.12 0.65 33.44 2123.37 1.17 0.37 35.54 -- -- -- -- Cent Q(k+1) x Q(k-1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k-1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k+1) x Q(k+1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k+1) -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Cent Q(k+1) x Q(k+1) - R -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- Prog Q(k) x Q(k) 2138.76 0.96 1.68 31.13 2139.04 0.95 0.85 30.70 2139.45 0.99 0.49 31.17 -- -- -- -- Prog Q(k) x Q(k) - R 2123.59 1.11 1.76 33.59 2123.85 1.10 0.88 32.85 2124.25 1.14 0.51 32.74 -- -- -- -- Prog RK2 Q(k) x Q(k) 2138.45 0.97 3.10 31.01 2138.43 0.97 1.30 31.31 2138.46 1.04 0.70 34.43 2138.49 0.99 0.60 36.46 Prog RK2 Q(k) x Q(k) - R 2123.32 1.12 3.18 33.55 2123.31 1.12 1.32 33.54 2123.29 1.19 0.70 36.03 2123.37 1.14 0.63 37.74 Prog RK4 Q(k) x Q(k) 2138.46 0.97 2.95 31.00 2138.44 0.97 1.46 31.24 2138.51 1.03 0.86 34.07 -- -- -- -- Prog RK4 Q(k) x Q(k) - R 2123.32 1.12 3.03 33.54 2123.31 1.12 1.52 33.48 2123.34 1.18 0.88 35.72 -- -- -- -- ) the detailed results.

Figure 11
Region of the hydrograph peak flow in the river's last reach for the Bates, Horritt and Fewtrell (2010)BATES, P. D.; HORRITT, M. S.; FEWTRELL, T. J. A simple inertial formulation of the shallow water equations for efficient two-dimensional flood inundation modelling. Journal of Hydrology, v. 387, n. 1, p. 33-45, 2010. http://dx.doi.org/10.1016/j.jhydrol.2010.03.027.
http://dx.doi.org/10.1016/j.jhydrol.2010...
and Prog RK2 Q(k)xQ(k) schemes with and without the hydraulic radius simplification considering a boundary condition with tide variation (the red line is subscribed by the green line).
Figure 12
Peak flow, volume conservation error, processing time and RMSE in relation to the time step for different formulations of numerical schemes considering a boundary condition with tide variation.

Unlike the other applications, in this test the Prog Q(k+1)xQ(k+1) scheme was not stable in any adopted time step. As shown, the scheme requires that the discriminant of the Bhaskara’s formula be positive and this criterion was not fulfilled even with a very small time step.

The most stable scheme was Prog RK2 Q(k)xQ(k) with C equal to 1.07 (time step of 4 minutes), followed by the other schemes, with C = 0.94 (3.5 min). Despite the improvement in stability, the Prog RK2 Q(k)xQ(k) scheme took, on average, twice the original scheme time and presented higher RMSE.

The volume conservation error and RMSE were equal in all stable schemes, approximately 0.97%, and had a slight increase when the calculation of the hydraulic radius was considered.

CONCLUSIONS

This paper compared the performance of new numerical schemes with inertial approximation of the Saint-Venant equations. The differences between the numerical schemes in terms of stability, peak flow, volume conservation error and RMSE were investigated. All stable schemes presented similar results to the hydrodynamic model HEC-RAS, and can be used with good accuracy in comparion to the model that considers the complete equations of Sant-Venant. The advantages of the proposed schemes are ease of implementation, robustness in the coupling with hydrodynamic models and the possibility of parallelization. However, more comparison of the inertial model with implicit and explicit numerical schemes is recommended for future studies.

With respect to performance measures, Table 3 presents the best scheme(s) in each metric and in each test. The difference of the peak flow and the volume conservation error between schemes can be considered negligible. Furthermore, although the original scheme is faster, the difference in processing time is very small in relation to the Prog Q(k)xQ(k) and the Prog Q(k+1)xQ(k+1) schemes.

Table 3
The best numerical scheme(s) in each test in relation to metrics.

It was observed for all stable numerical schemes that the peak flow approaches the HEC-RAS model's peak flow and that the RMSE decreases when the hydraulic radius calculation is considered without simplification. In this way, the inertial model formulation is more similar to what is done in the HEC-RAS, which does not simplify the hydraulic radius calculation. In addition, the processing time did not increase with this calculation and the volume conservation error showed a negligible increase.

In general, the scheme proposed by Almeida et al. (2012)ALMEIDA, G. A. M.; BATES, P.; FREER, J. E.; SOUVIGNET, M. Improving the stability of a simple formulation of the shallow water equations for 2-D flood modeling. Water Resources Research, v. 48, n. 5, p. W05528, 2012. http://dx.doi.org/10.1029/2011WR011570.
http://dx.doi.org/10.1029/2011WR011570...
presented no advantages in relation to the Bates, Horritt and Fewtrell (2010)BATES, P. D.; HORRITT, M. S.; FEWTRELL, T. J. A simple inertial formulation of the shallow water equations for efficient two-dimensional flood inundation modelling. Journal of Hydrology, v. 387, n. 1, p. 33-45, 2010. http://dx.doi.org/10.1016/j.jhydrol.2010.03.027.
http://dx.doi.org/10.1016/j.jhydrol.2010...
scheme. In all tests, both models presented instability at the same time step and both RMSE and peak flow were equivalent. The processing time was equal or greater in the Almeida et al. (2012)ALMEIDA, G. A. M.; BATES, P.; FREER, J. E.; SOUVIGNET, M. Improving the stability of a simple formulation of the shallow water equations for 2-D flood modeling. Water Resources Research, v. 48, n. 5, p. W05528, 2012. http://dx.doi.org/10.1029/2011WR011570.
http://dx.doi.org/10.1029/2011WR011570...
scheme.

Schemes with centered numerical derivatives showed instability for low C values in tests 1, 2, 5 and 6. Only in the applications of test 3 (high slope) and test 4 (very high slope) the Cent Q(k+1)xQ(k+1) scheme showed an advantage over the others. On the other hand, the second order approximation scheme, Prog RK2 Q(k)xQ(k), did not show stability in these tests (3 and 4). For this scheme, the stability was equivalent to original model in tests 2 and 5 and greater in tests 1 and 6. The disadvantage of the schemes with second and fourth order approximations is the longer processing time.

The parabola scheme, whose solution is a second degree equation (Prog Q(k+1)xQ(k+1)), presented advantages in relation to the Bates, Horritt and Fewtrell (2010)BATES, P. D.; HORRITT, M. S.; FEWTRELL, T. J. A simple inertial formulation of the shallow water equations for efficient two-dimensional flood inundation modelling. Journal of Hydrology, v. 387, n. 1, p. 33-45, 2010. http://dx.doi.org/10.1016/j.jhydrol.2010.03.027.
http://dx.doi.org/10.1016/j.jhydrol.2010...
scheme in most applications. In tests 1, 2, 3 and 4 the proposed model was the most stable, presenting good results with C values of 1.03, 1.04, 0.74 and 0.41, respectively, compared to the maximum values of 0.77, 0.90, 0.37 and 0.20 for the original scheme. In test 5 (dam and backwater), the stability was equal in both schemes (C = 0.97). However, in test 6 (tide), the parabola scheme was not stable at any time step. The proposed scheme is more unstable in negative flow situations.

The parabola scheme's processing time was equal to or slightly higher than the original model. This difference is in the tenths of a second scale and may be considered negligible. Moreover, in applications that do not occur negative flow, the greater stability makes this scheme faster because it allows to adopt a greater time step. As an example, in the supercritical regime application (test 4) the proposed scheme spent 0.58 sec (time step 160 sec) and the original scheme spent 1.12 sec (time step 80 sec). Regarding RMSE, despite the higher values in the parabola scheme in tests 2 (low slope) and 5 (dam and backwater), the difference with the original model is very small. The largest difference found was 2.3 m3/s (RMSE of 22.52 m3/s compared to 24.82 m3/s) in test 5.

Therefore, among all tests performed, we believe that the important information raised in the research is that the greater stability of the parabola scheme makes this scheme more attractive than the original model for use as flow propagation modulus in hydrological models. However this use is limited in some cases where there are negative flows, which reduces one of the main applications of models using the inertial model. This disadvantage could easily be overcome by using an adaptive mixed method which recognizes the possibility of a flow reversal and adapts temporarily to another scheme and then returns to the parabola scheme as soon as the negative flow disappears.

Finally, it is important to emphasize that, according to the results, the inclusion of the hydraulic radius calculation without simplification in the numerical schemes improved the results without increasing the computational time. That is, the inclusion of this consideration approximates the results of physical reality without additional performance costs.

Appendix A Supplementary data.

Table A1
Results of test 1 (real reach) for different numerical schemes.
Table A2
Results of test 2 (low slope) for different numerical schemes.
Table A3
Results of test 3 (high slope) for different numerical schemes.
Table A4
Results of test 4 (very high slope) for different numerical schemes.
Table A5
Results of test 5 (backwater) for different numerical schemes.
Table A6
Results of test 6 (tide) for different numerical schemes.

ACKNOWLEDGEMENTS

The authors would like to thank CNPq for the doctoral fellowship granted to the first author and anonymous reviewers for useful comments on the final draft.

REFERENCES

  • ALMEIDA, G. A. M.; BATES, P. Applicability of the local inertial approximation of the shallow water equations to flood modeling. Water Resources Research, v. 49, n. 8, p. 4833-4844, 2013. http://dx.doi.org/10.1002/wrcr.20366
    » http://dx.doi.org/10.1002/wrcr.20366
  • ALMEIDA, G. A. M.; BATES, P.; FREER, J. E.; SOUVIGNET, M. Improving the stability of a simple formulation of the shallow water equations for 2-D flood modeling. Water Resources Research, v. 48, n. 5, p. W05528, 2012. http://dx.doi.org/10.1029/2011WR011570
    » http://dx.doi.org/10.1029/2011WR011570
  • BATES, P. D.; HORRITT, M. S.; FEWTRELL, T. J. A simple inertial formulation of the shallow water equations for efficient two-dimensional flood inundation modelling. Journal of Hydrology, v. 387, n. 1, p. 33-45, 2010. http://dx.doi.org/10.1016/j.jhydrol.2010.03.027
    » http://dx.doi.org/10.1016/j.jhydrol.2010.03.027
  • BEVEN, K. J. Rainfall-runoff modelling: the primer New Jersey: John Wiley & Sons, 2011.
  • CHANSON, H. The hydraulics of open channel flow: an introduction 2nd ed. Oxford: Butterworth-Heinemann, 2004. 630 p.
  • CHAPRA, S. C.; CANALE, R. P. Numerical methods for engineers New York: McGraw-Hill, 2012.
  • COULTHARD, T. J.; NEAL, J. C.; BATES, P. D.; RAMIREZ, J.; ALMEIDA, G. A. M.; HANCOCK, G. R. Integrating the LISFLOOD‐FP 2D hydrodynamic model with the CAESAR model: implications for modelling landscape evolution. Earth Surface Processes and Landforms, v. 38, n. 15, p. 1897-1906, 2013. http://dx.doi.org/10.1002/esp.3478
    » http://dx.doi.org/10.1002/esp.3478
  • CUNGE, J. A. On the subject of a flood propagation computation method (Muskingum method). Journal of Hydrology Research, v. 7, n. 2, p. 205-230, 1969. http://dx.doi.org/10.1080/00221686909500264
    » http://dx.doi.org/10.1080/00221686909500264
  • CUNGE, J. A.; HOLLY, F. M.; VERWEY, A. Practical aspects of computational river hydraulics Boston: Pitman, 1980. 420 p.
  • DOTTORI, F.; TODINI, E. Developments of a flood inundation model based on the cellular automata approach: testing different methods to improve model performance. Physics and Chemistry of the Earth Parts A/B/C, v. 36, n. 7, p. 266-280, 2011. http://dx.doi.org/10.1016/j.pce.2011.02.004
    » http://dx.doi.org/10.1016/j.pce.2011.02.004
  • FAN, F. M.; PONTES, P. R. M.; PAIVA, R. C. D.; COLLISCHONN, W. Avaliação de um método de propagação de cheias em rios com aproximação inercial das equações de Saint-Venant. Revista Brasileira de Recursos Hídricos, v. 19, n. 4, p. 137-147, 2014. http://dx.doi.org/10.21168/rbrh.v19n4.p137-147
    » http://dx.doi.org/10.21168/rbrh.v19n4.p137-147
  • FREAD, D. L. Flow routing in handbook of hydrology. In: MAIDMENT, D. R. (Ed.). Handbook of hydrology New York: McGraw-Hill Inc., 1993. p. 10.1-10.36.
  • MONTEIRO, L. R.; FAN, F. M.; COLLISCHONN, W.; SCHETTINI, E. B. C. Simulação da onda superficial provocada pelo fechamento de comportas utilizando uma aproximação inercial da Equação de Saint-Venant/ Surface waves caused by gates closure simulation with an inertial approximation of Saint-Venant Equation. Revista Brasileira de Recursos Hídricos, v. 20, n. 4, p. 927-936, 2015. http://dx.doi.org/10.21168/rbrh.v20n4.p927-936
    » http://dx.doi.org/10.21168/rbrh.v20n4.p927-936
  • MONTERO, R. A.; SCHWANENBERG, D.; HATZ, M.; BRINKMANN, M. Simplified hydraulic modelling in model predictive control of flood mitigation measures along rivers. Journal of Applied Water Engineering and Research, v. 1, n. 1, p. 17-27, 2013. http://dx.doi.org/10.1080/23249676.2013.827897
    » http://dx.doi.org/10.1080/23249676.2013.827897
  • NEAL, J. C.; SCHUMANN, G. J.-P.; FEWTRELL, T. J.; BUDIMIR, M.; BATES, P. D.; MASON, D. Evaluating a new LISFLOOD-FP formulation using data for the summer 2007 floods in Tewkesbury, UK. Journal of Flood Risk Management, v. 4, n. 2, p. 88-95, 2011. http://dx.doi.org/10.1111/j.1753-318X.2011.01093.x
    » http://dx.doi.org/10.1111/j.1753-318X.2011.01093.x
  • NEAL, J.; SCHUMANN, G.; BATES, P. A subgrid channel model for simulating river hydraulics and floodplain inundation over large and data sparse areas. Water Resources Research, v. 48, n. 11, 2012. http://dx.doi.org/10.1029/2012WR012514
    » http://dx.doi.org/10.1029/2012WR012514
  • NÉELZ, S.; PENDER, G. (2009). Desktop review of 2D hydraulic modelling packages Bristol: Defra/Environment Agency. Science Report SC080035.
  • PONCE, V. M. Engineering hydrology, principles and practices New Jersey: Prentice Hall, 1989. 640 p.
  • PONTES, P. R. M.; COLLISCHONN, W.; FAN, F.M.; PAIVA, R.C.D.; BUARQUE, D.C. Modelagem hidrológica e hidráulica de grande escala com propagação inercial de vazões. Revista Brasileira de Recursos Hídricos, vol. 20, n. 4, p. 888-904, 2015. http://dx.doi.org/10.21168/rbrh.v20n4.p888-904
    » http://dx.doi.org/10.21168/rbrh.v20n4.p888-904
  • PONTES, P. R. M.; FAN, F. M.; FLEISCHMANN, A. S.; PAIVA, R. C. D.; BUARQUE, D. C.; SIQUEIRA, V. A.; JARDIM, P. F.; SORRIBAS, M. V.; COLLISCHONN, W. MGB-IPH model for hydrological and hydraulic simulation of large floodplain river systems coupled with open source GIS. Environmental Modelling & Software, v. 94, p. 1-20, 2017. http://dx.doi.org/10.1016/j.envsoft.2017.03.029
    » http://dx.doi.org/10.1016/j.envsoft.2017.03.029
  • POPESCU, I. Computational hydraulics London: IWA Publishing, 2014.
  • SAMPSON, C. C.; SMITH, A. M.; BATES, P. D.; NEAL, J. C.; ALFIERI, L.; FREER, J. E. A high‐resolution global flood hazard model. Water Resources Research, v. 51, n. 9, p. 7358-7381, 2015. PMid:27594719. http://dx.doi.org/10.1002/2015WR016954
    » http://dx.doi.org/10.1002/2015WR016954
  • TUCCI, C. E. M. Modelos hidrológicos 2. ed. Porto Alegre: Editora da UFRGS, 2005.
  • US ARMY CORPS OF ENGINEERS – USACE. HEC-RAS River Analysis System Hydraulic Reference Manual. Version 5.0 Davis: Hydrologic Engineering Center, 2016.
  • YAMAZAKI, D.; ALMEIDA, G. A. M.; BATES, P. D. Improving computational efficiency in global river models by implementing the local inertial flow equation and a vector-based river network map. Water Resources Research, v. 49, n. 11, p. 7221-7235, 2013.
  • YEN, B. C.; TSAI, C. W.-S. Noninertial wave vs. diffusion wave in flood routing. Journal of Hydrology, v. 244, n. 1-2, p. 97-104, 2001. http://dx.doi.org/10.1016/S0022-1694(00)00422-4
    » http://dx.doi.org/10.1016/S0022-1694(00)00422-4

Publication Dates

  • Publication in this collection
    2018

History

  • Received
    08 May 2017
  • Reviewed
    15 Oct 2017
  • Accepted
    10 Dec 2017
Associação Brasileira de Recursos Hídricos Av. Bento Gonçalves, 9500, CEP: 91501-970, Tel: (51) 3493 2233, Fax: (51) 3308 6652 - Porto Alegre - RS - Brazil
E-mail: rbrh@abrh.org.br