Central-upwind scheme for shallow water equations with discontinuous bottom topography

Finite-volume central-upwind schemes for shallow water equations were proposed in [A. Kurganov and G. Petrova, Commun. Math. Sci., 5 (2007), 133–160]. These schemes are capable of maintaining “lake-at-rest” steady states and preserving the positivity of the computed water depth. The well-balanced and positivity preserving features of the central-upwind schemes are achieved, in particular, by using continuous piecewise linear interpolation of the bottom topography function. However, when the bottom function is discontinuous or a model with a moving bottom topography is studied, the continuous piecewise linear approximationmay not be sufficiently accurate and robust.     

A well-balanced van Leer-type numerical scheme for shallow water equations with variable topography

A well-balanced van Leer-type numerical scheme for the shallow water equations with variable topography is presented. The model involves a nonconservative term, which often makes standard schemes difficult to approximate solutions in certain regions. The construction of our scheme is based on exact solutions in computational form of local Riemann problems. Numerical tests are conducted, where comparisons between this van Leer-type scheme and a Godunov-type scheme are provided. Data for the tests are taken in both the subcritical region as well as supercritical region. Especially, tests for resonant cases where the exact solutions contain coinciding waves are also investigated. All numerical tests show that each of these two methods can give a good accuracy, while the van Leer -type scheme gives a better accuracy than the Godunov-type scheme. Furthermore, it is shown that the van Leer-type scheme is also well-balanced in the sense that it can capture exactly stationary contact discontinuity waves.     

Zeroth-order conservation laws of two-dimensional shallow water equations with variable bottom topography

We classify zeroth-order conservation laws of systems from the class of two-dimensional shallow water equations with variable bottom topography using an optimized version of the method of furcate splitting. The classification is carried out up to equivalence generated by the equivalence group of this class. We find additional point equivalences between some of the listed cases of extensions of the space of zeroth-order conservation laws, which are inequivalent up to transformations from the equivalence group. Hamiltonian structures of systems of shallow water equations are used for relating the classification of zeroth-order conservation laws of these systems to the classification of their Lie symmetries. We also construct generating sets of such conservation laws under action of Lie symmetries.     

A parallel well-balanced finite volume method for shallow water equations with topography on the cubed-sphere

A finite volume scheme for the global shallow water model on the cubed-sphere mesh is proposed and studied in this paper. The new cell-centered scheme is based on Osher’s Riemann solver together with a high-order spatial reconstruction. On each patch interface of the cubed-sphere only one layer of ghost cells is needed in the scheme and the numerical flux is calculated symmetrically across the interface to ensure the numerical conservation of total mass. The discretization of the topographic term in the equation is properly modified in a well-balanced manner to suppress spurious oscillations when the bottom topography is non-smooth. Numerical results for several test cases including a steady-state nonlinear geostrophic flow and a zonal flow over an isolated mountain are provided to show the flexibility of the scheme. Some parallel implementation details as well as some performance results on a parallel supercomputer with more than one thousand processor cores are also provided.     

High-order discontinuous element-based schemes for the inviscid shallow water equations: Spectral multidomain penalty and discontinuous Galerkin methods

Two commonly used types of high-order-accuracy element-based schemes, collocation-based spectral multidomain penalty methods (SMPM) and nodal discontinuous Galerkin methods (DGM), are compared in the framework of the inviscid shallow water equations. Differences and similarities in formulation are identified, with the primary difference being the dissipative term in the Rusanov form of the numerical flux for the DGM that provides additional numerical stability; however, it should be emphasized that to arrive at this equivalence between SMPM and DGM requires making specific choices in the construction of both methods; these choices are addressed. In general, both methods offer a multitude of choices in the penalty terms used to introduce boundary conditions and stabilize the numerical solution. The resulting specialized class of SMPM and DGM are then applied to a suite of six commonly considered geophysical flow test cases, three linear and three non-linear; we also include results for a classical continuous Galerkin (i.e., spectral element) method for comparison. Both the analysis and numerical experiments show that the SMPM and DGM are essentially identical; both methods can be shown to be equivalent for very special choices of quadrature rules and Riemann solvers in the DGM along with special choices in the type of penalty term in the SMPM. Although we only focus our studies on the inviscid shallow water equations the results presented should be applicable to other systems of nonlinear hyperbolic equations (such as the compressible Euler equations) and extendable to the compressible and incompressible Navier-Stokes equations, where viscous terms are included.     

Well-posedness for the fifth-order shallow water equations

The Cauchy problems for some kind of fifth-order shallow water equations

     

Meshless method for shallow water equations with free surface flow

Solving problems with free surface often encounters numerical difficulties related to excessive mesh distortion as is the case of dambreak or breaking waves. In this paper the Natural element method (NEM) is used to simulate a 2D shallow water flows in the presence of theses strong gradients. This particle-based method used a fully Lagrangian formulation based on the notion of natural neighbors. In the present study we consider the full non-linear set of Shallow Water Equations, with a transient flow under the Coriolis effect. For the numerical treatment of the nonlinear terms we used a Lagrangian technique based on the method of characteristics. This will allow avoiding divergence of Newton-Raphson scheme, when dealing with the convective terms. We also define a thin area close to the boundaries and a computational domain dedicated for nodal enrichment at each time step. Two numerical test cases were performed to verify the well-founded hopes for the future of this method in real applications.     

Numerical approximations of one-dimensional linear conservation equations with discontinuous coefficients

Conservative linear equations arise in many areas of application, including continuum mechanics or high-frequency geometrical optics approximations. This kind of equation admits most of the time solutions which are only bounded measures in the space variable known as duality solutions. In this paper, we study the convergence of a class of finite-difference numerical schemes and introduce an appropriate concept of consistency with the continuous problem. Some basic examples including computational results are also supplied.

     

Cauchy problem for viscous rotating shallow water equations

We consider the Cauchy problem for a viscous compressible rotating shallow water system with a third-order surface-tension term involved, derived recently in the modeling of motions for shallow water with free surface in a rotating sub-domain Marche (2007) [19]. The global existence of the solution in the space of Besov type is shown for initial data close to a constant equilibrium state away from the vacuum. Unlike the previous analysis about the compressible fluid model without Coriolis forces, see for instance Danchin (2000) [10], Haspot (2009) [16], the rotating effect causes a coupling between two parts of Hodge's decomposition of the velocity vector field, and additional regularity is required in order to carry out the Friedrichs' regularization and compactness arguments.     

Numerical model for the shallow water equations on a curvilinear grid with the preservation of the Bernoulli integral

Methodological aspects concerning the construction of a two-dimensional numerical model for reservoir flows based on the shallow water equations are considered. A numerical scheme is constructed by applying the control volume method on staggered grids in combination with the Bernoulli integral, which is used to interpolate the desired fields inside a grid cell. The implementation of the method yields a monotone numerical scheme. The results of numerical integration are compared with the exact solution.     

Numerical modelling of free-surface shallow flows over irregular topography with complex geometry

A well-balanced Godunov-type finite volume algorithm is developed for modelling free-surface shallow flows over irregular topography with complex geometry. The algorithm is based on a new formulation of the classical shallow water equations in hyperbolic conservation form. Unstructured triangular grids are used to achieve the adaptability of the grid to the geometry of the problem and to facilitate localised refinement. The numerical fluxes are calculated using HLLC approximate Riemann solver, and the MUSCL-Hancock predictor–corrector scheme is adopted to achieve the second-order accuracy both in space and in time where the solutions are continuous, and to achieve high-resolution results where the solutions are discontinuous. The novelties of the algorithm include preserving well-balanced property without any additional correction terms and the wet/dry front treatments. The good performance of the algorithm is demonstrated by comparing numerical and theoretical results of several benchmark problems, including the preservation of still water over a two-dimensional hump, the idealised dam-break flow over a frictionless flat rectangular channel, the circular dam-break, and the shock wave from oblique wall. Besides, two laboratory dam-break cases are used for model validation. Furthermore, a practical application related to dam-break flood wave propagation over highly irregular topography with complex geometry is presented. The results show that the algorithm can correctly account for free-surface shallow flows with respect to its effectiveness and robustness thus has bright application prospects.     

On the shallow water equations with periodicity in one direction

The linearized 2D shallow water equations, supplemented with suitable boundary conditions in one direction and periodicity in the other direction are considered. The well-posedness of this mixed boundary value problem is proven, using the linear semi-group theory. The well-posedness of the totally periodic problem is also proven.     

Regularized shallow water equations and an efficient method for numerical simulation of shallow water flows

Regularized shallow water equations are derived as based on a regularization of the Navier-Stokes equations in the form of quasi-gasdynamic and quasi-hydrodynamic equations. Efficient finite-difference algorithms based on the regularized shallow water equations are proposed for the numerical simulation of shallow water flows. The capabilities of the model are examined by computing a test Riemann problem, the flow over an obstacle, and asymmetric dam break.     

Space–time discontinuous Galerkin finite element method for shallow water flows

A space–time discontinuous Galerkin (DG) finite element method is presented for the shallow water equations over varying bottom topography. The method results in nonlinear equations per element, which are solved locally by establishing the element communication with a numerical HLLC flux. To deal with spurious oscillations around discontinuities, we employ a dissipation operator only around discontinuities using Krivodonova's discontinuity detector. The numerical scheme is verified by comparing numerical and exact solutions, and validated against a laboratory experiment involving flow through a contraction. We conclude that the method is second order accurate in both space and time for linear polynomials.     

Flatness based trajectory planning for the shallow water equations

A flatness based feed-forward control design approach for an open channel flow modelled by the shallow water equations is discussed. The control input consists of the adjustable height of opening of a delimiting sluice gate. In order to compute the trajectory of the control input, the problem is reduced to an initial value problem w.r.t. the spatial coordinate. The solution to this problem can be obtained with the method of characteristics. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Absorbing boundary conditions for the linearized shallow water equations

A mixed problem imitating the Cauchy problem for the linearized shallow water equations is considered. This problem is also a mixed problem with perfectly absorbing conditions (cp. [1], [3]). An exact formula for the conditions has been given.     

On a numerical flux for the shallow water equations

The finite volume scheme of Vijayasundaram and Osher-Solomon type for shallow water equations are proposed. The numerical results with discontinuous initial condition and the comparison with Lax-Friedrichs numerical flux are presented for homogeneous case. The extension of the method for the inhomogeneous case is described.