scholarly journals Convergence Improved Lax-Friedrichs Scheme Based Numerical Schemes and Their Applications in Solving the One-Layer and Two-Layer Shallow-Water Equations

2015 ◽  
Vol 2015 ◽  
pp. 1-10
Author(s):  
Xinhua Lu ◽  
Bingjiang Dong ◽  
Bing Mao ◽  
Xiaofeng Zhang
Keyword(s):  
High Resolution ◽  
Shallow Water ◽  
Shallow Water Equations ◽  
Small Time ◽  
Numerical Schemes ◽  
Time Step ◽  
Numerical Diffusion ◽  
First Order ◽  
Small Time Step ◽  
The One

The first-order Lax-Friedrichs (LF) scheme is commonly used in conjunction with other schemes to achieve monotone and stable properties with lower numerical diffusion. Nevertheless, the LF scheme and the schemes devised based on it, for example, the first-order centered (FORCE) and the slope-limited centered (SLIC) schemes, cannot achieve a time-step-independence solution due to the excessive numerical diffusion at a small time step. In this work, two time-step-convergence improved schemes, the C-FORCE and C-SLIC schemes, are proposed to resolve this problem. The performance of the proposed schemes is validated in solving the one-layer and two-layer shallow-water equations, verifying their capabilities in attaining time-step-independence solutions and showing robustness of them in resolving discontinuities with high-resolution.

scholarly journals A well-balanced spectral volume scheme with the wetting–drying property for the shallow-water equations

2012 ◽  
Vol 14 (3) ◽  
pp. 745-760 ◽  
Author(s):  
Luca Cozzolino ◽  
Renata Della Morte ◽  
Giuseppe Del Giudice ◽  
Anna Palumbo ◽  
Domenico Pianese
Keyword(s):  
Shallow Water ◽  
Shallow Water Equations ◽  
Step Length ◽  
Swash Zone ◽  
Time Step ◽  
Third Order ◽  
Spectral Volume ◽  
The Stability ◽  
The One ◽  
Volume Method

The shallow-water equations are widely used to model surface water bodies, such as lakes, rivers and the swash zone in coastal flows. Physically congruent solutions are characterized by non-negative water depth, and many numerical methods may fail to preserve this property at the discrete level when moving wet–dry transitions are present in the physical domain. In this paper, we present a spectral-volume method for the approximate solution of the one-dimensional shallow-water equations, which is third-order accurate in wet regions, far from discontinuities, and which is well balanced for water at rest states: the stability of the solution is ensured if reconstruction and limitation of variables preserves non-negativity of the depth and a suitable constraint for the time step length is satisfied. A number of numerical experiments are reported, showing the promising capabilities of the model to solve problems with non-trivial topographies and friction.

scholarly journals Flood wave propagation in steep mountain rivers

2012 ◽  
Vol 15 (1) ◽  
pp. 120-137 ◽  
Author(s):  
Gabriella Petaccia ◽  
Luigi Natale ◽  
Fabrizio Savi ◽  
Mirjana Velickovic ◽  
Yves Zech ◽  
...  
Keyword(s):  
Wave Propagation ◽  
Shallow Water ◽  
Shallow Water Equations ◽  
Source Term ◽  
Term Treatment ◽  
Flood Wave ◽  
Numerical Schemes ◽  
Steady State Condition ◽  
One Dimensional ◽  
The One

Most of the recent developments concerning efficient numerical schemes to solve the shallow-water equations in view of real world flood modelling purposes concern the two-dimensional form of the equations or the one-dimensional form written for rectangular, unit-width channels. Extension of these efficient schemes to the one-dimensional cross-sectional averaged shallow-water equations is not straightforward, especially when complex natural topographies are considered. This paper presents different formulations of numerical schemes based on the HLL (Harten–Lax–van Leer) solver, and the adaptation of the topographical source term treatment when cross-sections of arbitrary shape are considered. Coupled and uncoupled formulations of the equations are considered, in combination with centred and lateralised source term treatment. These schemes are compared to a numerical solver of Lax Friedrichs type based on a staggered grid. The proposed schemes are first tested against two theoretical benchmark tests and then applied to the Brembo River, an Italian alpine river, firstly simulating a steady-state condition and secondly reproducing the 2002 flood wave propagation.

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

RBRH ◽  
2018 ◽  
Vol 23 (0) ◽  
Author(s):  
Alice César Fassoni-Andrade ◽  
Fernando Mainardi Fan ◽  
Walter Collischonn ◽  
Artur César Fassoni ◽  
Rodrigo Cauduro Dias de Paiva
Keyword(s):  
Numerical Stability ◽  
Flood Routing ◽  
Computational Time ◽  
Numerical Schemes ◽  
Time Step ◽  
Simple Formulation ◽  
One Dimensional ◽  
Saint Venant Equations ◽  
Original Scheme ◽  
The One

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.

A Balanced Approximation of the One-Layer Shallow-Water Equations on a Sphere

2009 ◽  
Vol 66 (6) ◽  
pp. 1735-1748 ◽  
Author(s):  
W. T. M. Verkley
Keyword(s):  
Shallow Water ◽  
Potential Vorticity ◽  
Shallow Water Equations ◽  
Rossby Waves ◽  
Coriolis Parameter ◽  
Beta Plane ◽  
Propagating Waves ◽  
Barotropic Vorticity Equation ◽  
Barotropic Vorticity ◽  
The One

Abstract A global version of the equivalent barotropic vorticity equation is derived for the one-layer shallow-water equations on a sphere. The equation has the same form as the corresponding beta plane version, but with one important difference: the stretching (Cressman) term in the expression of the potential vorticity retains its full dependence on f 2, where f is the Coriolis parameter. As a check of the resulting system, the dynamics of linear Rossby waves are considered. It is shown that these waves are rather accurate approximations of the westward-propagating waves of the second class of the original shallow-water equations. It is also concluded that for Rossby waves with short meridional wavelengths the factor f 2 in the stretching term can be replaced by the constant value f02, where f0 is the Coriolis parameter at ±45° latitude.

Load Discontinuity and Its Remedy in Step-by-Step Integration

2006 ◽  
Author(s):  
Shuenn-Yih Chang ◽  
Chiu-Li Huang
Keyword(s):  
Numerical Solution ◽  
Small Time ◽  
Step Size ◽  
Time Step ◽  
Integration Procedure ◽  
Amplitude Distortion ◽  
Small Time Step

The discontinuity at the end of an impulse will lead to an extra impulse and thus an extra displacement. Consequently, an amplitude distortion is introduced in the numerical solution. The difficulty arising from the discontinuity at the end of an impulse can be overcome by using a very small time step to perform the step-by-step integration since it reduces the extra impulse and thus extra displacement. However, computational efforts might be significantly increased since the small time step is performed for a complete step-by-step integration procedure. A remedy is devised to computationally efficiently overcome this difficulty by using a very small time step immediately upon termination of the applied impulse. This is because that the extra impulse caused by the discontinuity is almost proportional to the discontinuity value at the end of the impulse and the step size. The feasibility of this proposed remedy is analytically and numerically confirmed herein.

Sign in / Sign up

Export Citation Format

Share Document