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

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.

2015 ◽  
Vol 2015 ◽  
pp. 1-10
Author(s):  
Xinhua Lu ◽  
Bingjiang Dong ◽  
Bing Mao ◽  
Xiaofeng Zhang

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.


1987 ◽  
Vol 183 ◽  
pp. 343-363 ◽  
Author(s):  
P. Ripa

The one-layer reduced gravity (or ‘shallow water’) equations in the f-plane have solutions such that the active layer is horizontally bounded by an ellipse that rotates steadily. In a frame where the height contours are stationary, fluid particles move along similar ellipses with the same revolution period. Both motions (translation along an elliptical path and precession of that orbit) are anticyclonic and their frequencies are not independent; a Rossby number (R0) based on the combination of both of them is bounded by unity. These solutions may be taken, with some optimism, as a model of ocean warm eddies; their stability is studied here for all values of R0 and of the ellipse eccentricity (these two parameters determine uniquely the properties of the solution).Sufficient stability conditions are derived from the integrals of motion; f-plane flows that satisfy them must be either axisymmetric or parallel. For the model vortex, the circular case simply corresponds to a solid-body rotation, and is found to be stable to finite-amplitude perturbations for all values of R0. This includes R0 > ½, which implies an anticyclonic absolute vorticity.The stability of the truly elliptical cases are studied in the normal modes sense. The height perturbation is an n-order polynomial of the horizontal coordinates; the cases for 0 ≤ n ≤ 6 are analysed, for all possible values of the Rossby number and of the eccentricity. All eddies are stable to perturbations with n ≤ 2. (A property of the shallow-water equations, probably related to the last result, is that a general finite-amplitude n-order field is an exact nonlinear solution for n ≤ 2.) Many vortices - noticeably the more eccentric ones - are unstable to perturbations with n ≥ 3; growth rates are O(R02f) where f is the Coriolis parameter.


2009 ◽  
Vol 66 (6) ◽  
pp. 1735-1748 ◽  
Author(s):  
W. T. M. Verkley

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.


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

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


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

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


Sign in / Sign up

Export Citation Format

Share Document