Weak local residuals as smoothness indicators for the shallow water equations

The system of shallow water equations admits infinitely many conservation laws. This paper investigates weak local residuals as smoothness indicators of numerical solutions to the shallow water equations. To get a weak formulation, a test function and integration are introduced into the shallow water equations. We use a finite volume method to solve the shallow water equations numerically. Based on our numerical simulations, the weak local residual of a simple conservation law with a simple test function is identified to be the best as a smoothness indicator.     

Simulation of the spread of a viscous fluid using a bidimensional shallow water model

In this paper we propose a numerical method to solve the Cauchy problem based on the viscous shallow water equations in an horizontally moving domain. More precisely, we are interested in a flooding and drying model, used to modelize the overflow of a river or the intrusion of a tsunami on ground. We use a nonconservative form of the two-dimensional shallow water equations, in eight velocity formulation and we build a numerical approximation, based on the Arbitrary lagrangian eulerian formulation, in order to compute the solution in the moving domain.     

Simulation of shallow water flows with shoaling areas and bottom discontinuities

A numerical method based on a second-order accurate Godunov-type scheme is described for solving the shallow water equations on unstructured triangular-quadrilateral meshes. The bottom surface is represented by a piecewise linear approximation with discontinuities, and a new approximate Riemann solver is used to treat the bottom jump. Flows with a dry sloping bottom are computed using a simplified method that admits negative depths and preserves the liquid mass and the equilibrium state. The accuracy and performance of the approach proposed for shallow water flow simulation are illustrated by computing one- and two-dimensional problems.     

Group properties and conservation laws for nonlocal shallow water wave equation

Symmetry groups, symmetry reductions, optimal system, conservation laws and invariant solutions of the shallow water wave equation with nonlocal term are studied. First, Lie symmetries based on the invariance criterion for nonlocal equations and the solution approach for nonlocal determining equations are found and then the reduced equations and optimal system are obtained. Finally, new conservation laws are generated and some similarity solutions for symmetry reduction forms are discussed.     

Wave breaking for a periodic shallow water equation

The focus of this paper is on the blow-up of a recently derived one-dimensional shallow water equation which is formally integrable and can be obtained by approximating directly the Hamiltonian for Euler's equations in the shallow water regime. Some new criteria guaranteeing the development of singularities in finite time for strong solutions with regular initial data are obtained for the periodic case.     

Regularized shallow water equations for numerical simulation of flows with a moving shoreline

A numerical algorithm for simulating free-surface flows based on regularized shallow water equations is adapted to flows involving moving dry-bed areas. Well-balanced versions of the algorithm are constructed. Test computations of flows with dry-bed areas in the cases of water runup onto a plane beach and a constant-slope beach are presented. An example of tsunami simulation is given.     

Computations of transport of pollutant in shallow water

The focus of this paper is to simulate the transport of a passive pollutant by a flow modelled by the two-dimensional shallow water equations. Considering the friction terms, new model for simulating the steady and unsteady transport of pollutant is established. Then the adaptive semi-discrete central-upwind scheme based on central weighted essentially non-oscillatory reconstruction is utilized for simulating the two-dimensional steady and unsteady transport of pollutant. The non-oscillatory behavior and accuracy of the scheme are demonstrated by the numerical result.     

Water waves over a rough bottom in the shallow water regime

This is a study of the Euler equations for free surface water waves in the case of varying bathymetry, considering the problem in the shallow water scaling regime. In the case of rapidly varying periodic bottom boundaries this is a problem of homogenization theory. In this setting we derive a new model system of equations, consisting of the classical shallow water equations coupled with nonlocal evolution equations for a periodic corrector term. We also exhibit a new resonance phenomenon between surface waves and a periodic bottom. This resonance, which gives rise to secular growth of surface wave patterns, can be viewed as a nonlinear generalization of the classical Bragg resonance. We justify the derivation of our model with a rigorous mathematical analysis of the scaling limit and the resulting error terms. The principal issue is that the shallow water limit and the homogenization process must be performed simultaneously. Our model equations and the error analysis are valid for both the two- and the three-dimensional physical problems.     

Self-similar solutions and blow-up phenomena for a two-component shallow water system

In this article, we consider a two-component nonlinear shallow water system, which includes the famous 2-component Camassa-Holm and Degasperis-Procesi equations as special cases. The local well-posedess for this equations is established. Some sufficient conditions for blow-up of the solutions in finite time are given. Moreover, by separation method, the self-similar solutions for the nonlinear shallow water equations are obtained, and which local or global behavior can be determined by the corresponding Emden equation.     

Balance-characteristic scheme as applied to the shallow water equations over a rough bottom

The CABARET scheme is used for the numerical solution of the one-dimensional shallow water equations over a rough bottom. The scheme involves conservative and flux variables, whose values at a new time level are calculated by applying the characteristic properties of the shallow water equations. The scheme is verified using a series of test and model problems.     

Two-dimensional numerical modelling of shallow water flows over multilayer movable beds

The two-dimensional modelling of shallow water flows over multi-sediment erodible beds is presented. A novel approach is developed for the treatment of multiple sediment types in morphodynamics. The governing equations include the two-dimensional shallow water equations for hydrodynamics, an Exner-type equation for morphodynamics, a two-dimensional transport equation for the suspended sediments, and a set of empirical equations for entrainment and deposition. Multilayer sedimentary beds are formed of different erodible soils with sediment properties and new exchange conditions between the bed layers are developed for the model. The coupled equations yield a hyperbolic system of balance laws with source terms. As a numerical solver for the system, we implement a fast finite volume characteristics method. The numerical fluxes are reconstructed using the method of characteristics which employs projection techniques. The proposed finite volume solver is simple to implement, satisfies the conservation property and can be used for two-dimensional sediment transport problems in non-homogeneous isotropic beds without need of complicated three-dimensional equations. To assess the performance of the proposed models, we present numerical results for a wide variety of shallow water flows over sedimentary layers. Comparisons to experimental data for dam-break problems over movable beds are also included in this study.     

A simple finite volume method for the shallow water equations

We present a new finite volume method for the numerical solution of shallow water equations for either flat or non-flat topography. The method is simple, accurate and avoids the solution of Riemann problems during the time integration process. The proposed approach consists of a predictor stage and a corrector stage. The predictor stage uses the method of characteristics to reconstruct the numerical fluxes, whereas the corrector stage recovers the conservation equations. The proposed finite volume method is well balanced, conservative, non-oscillatory and suitable for shallow water equations for which Riemann problems are difficult to solve. The proposed finite volume method is verified against several benchmark tests and shows good agreement with analytical solutions.     

Evolution of weak discontinuities in shallow water equations

In this paper, we determine the critical time, when a weak discontinuity in the shallow water equations culminates into a bore. Invariance group properties of the governing system of partial differential equations (PDEs), admitting Lie group of point transformations with commuting infinitesimal operators, are presented. Some appropriate canonical variables are characterized that transform equations at hand to an equivalent form, which admits non-constant solutions. The propagation of weak discontinuities is studied in the medium characterized by the particular solution of the governing system.     

Symmetries of shallow water equations on a rotating plane

We consider a system of nonlinear differential equations which describes the spatial motion of an ideal incompressible fluid on a rotating plane in the shallow water approximation and a more general system of the theory of long waves which takes into account the specifics of shear flows. Using the group analysis methods, we calculate the 9-dimensional Lie algebras of infinitesimal operators admissible by the models. We establish an isomorphism of these Lie algebras with a known Lie algebra of operators admissible by the system of equations for the two-dimensional isentropic motions of a polytropic gas with the adiabatic exponent γ = 2. The nontrivial symmetries of the models under consideration enable us to carry out the group generation of the solutions. The class of stationary solutions to the equations of rotating shallow water transforms into a new class of periodic solutions.     

On curl-preserving finite volume discretizations for shallow water equations

The preservation of intrinsic or inherent constraints, like divergence-conditions, has gained increasing interest in numerical simulations of various physical evolution equations. In Torrilhon and Fey, SIAM J. Numer. Anal. (42/4) 2004, a general framework is presented how to incorporate the preservation of a discrete constraint into upwind finite volume methods. This paper applies this framework to the wave equation system and the system of shallow water equations. For the wave equation a curl-preservation for the momentum variable is present and easily identified. The preservation in case of the shallow water system is more involved due to the presence of convection. It leads to the vorticity evolution as generalized curl-constraint. The mechanisms of vorticity generation are discussed.For the numerical discretization special curl-preserving flux distributions are discussed and their incorporation into a finite volume scheme described. This leads to numerical discretizations which are exactly curl-preserving for a specific class of discrete curl-operators.The numerical experiments for the wave equation show a significant improvement of the new method against classical schemes. The extension of the curl-free numerical discretization to the shallow water case is possible after isolating the pressure flux. Simulation examples demonstrate the influence of the modification. The vortex structure is more clearly resolved.     

A robust central scheme for the shallow water flows with an abrupt topography based on modified hydrostatic reconstructions

We study a second-order central scheme for the shallow water flows with a discontinuous bottom topography based on modified hydrostatic reconstructions (HRs). The first HR scheme was proposed in Audusse et al, which may be missing the effect of the large discontinuous bottom topography. We introduce a modified HR method to cope with this numerical difficulty. The new scheme is well-balanced for still water solutions and can guarantee the positivity of the water depth. Finally, several numerical results of classical problems of the shallow water equations confirmed these properties of the new scheme. Especially, the new scheme yields superior results for the shallow water downhill flow over a step.     

A new highly nonlinear shallow water wave equation

We derive a quasilinear shallow water equation directly from the governing equations for gravity water waves within a certain regime for large-amplitude waves which has not been studied so far. Furthermore, we demonstrate local well-posedness of the corresponding Cauchy problem and finally discuss some aspects of the blowup behavior of solutions.     

Exact solutions of one-dimensional nonlinear shallow water equations over even and sloping bottoms

We establish an equivalence of two systems of equations of one-dimensional shallow water models describing the propagation of surface waves over even and sloping bottoms. For each of these systems, we obtain formulas for the general form of their nondegenerate solutions, which are expressible in terms of solutions of the Darboux equation. The invariant solutions of the Darboux equation that we find are simplest representatives of its essentially different exact solutions (those not related by invertible point transformations). They depend on 21 arbitrary real constants; after “proliferation” formulas derived by methods of group theory analysis are applied, they generate a 27-parameter family of essentially different exact solutions. Subsequently using the derived infinitesimal “proliferation” formulas for the solutions in this family generates a denumerable set of exact solutions, whose linear span constitutes an infinite-dimensional vector space of solutions of the Darboux equation. This vector space of solutions of the Darboux equation and the general formulas for nondegenerate solutions of systems of shallow water equations with even and sloping bottoms give an infinite set of their solutions. The “proliferation” formulas for these systems determine their additional nondegenerate solutions. We also find all degenerate solutions of these systems and thus construct a database of an infinite set of exact solutions of systems of equations of the one-dimensional nonlinear shallow water model with even and sloping bottoms.     

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.