Discontinuous boundary implementation for the shallow water equations

Quasi‐bubble finite element approximations to the shallow water equations are investigated focusing on implementations of the surface elevation boundary condition. We first demonstrate by numerical results that the conventional implementation of the boundary condition degrades the accuracy of the velocity solution. It is also shown that the degraded velocity leads to a critical instability if the advection term is present in the momentum equation. Then we propose an alternative implementation for the boundary condition. We refer to this alternative implementation as a discontinuous boundary (DB) implementation because it introduces at each boundary node two independent mass–flux values that result in a discontinuity at the boundary. Numerical results show that the proposed DB implementation is consistent, stabilizes the quasi‐bubble scheme, and leads to second‐order accuracy at the surface elevation specified boundary. Copyright © 2004 John Wiley & Sons, Ltd.     

Nonlinear-dispersive shallow water equations on a rotating sphere and conservation laws

Nonlinear dispersive shallow water equations on a sphere are obtained without using the potential flow assumption. Boussinesq-type equations for weakly nonlinear waves over a moving bottom are derived. It is found that the total energy balance holds for all obtained nonlinear dispersive equations on a sphere.     

Full nonlinear dispersion model of shallow water equations on a rotating sphere

Nonlinear dispersion shallow water equations are derived, which describe propagation of long surface waves on a spherical surface with allowance for rotation of the Earth and mobility of the ocean bottom. Derivation of these equations is based on expanding the solution of hydrodynamic equations on a sphere in small parameters depending on the relative thickness of the water layer and dispersion of surface waves.     

Equations of the shallow water model on a rotating attracting sphere. 1. Derivation and general properties

A shallow water model on a rotating attracting sphere is proposed to describe large-scale motions of the gas in planetary atmospheres and of the liquid in the world ocean. The equations of the model coincide with the equations of gas-dynamic of a polytropic gas in the case of spherical gas motions on the surface of a rotating sphere. The range of applicability of the model is discussed, and the conservation of potential vorticity along the trajectories is proved. The equations of stationary shallow water motions are presented in the form of Bernoulli and potential vorticity integrals, which relate the liquid depth to the stream function. The simplest stationary solutions that describe the equilibrium state differing from the spherically symmetric state and the zonal flows along the parallels are found. It is demonstrated that the stationary equations of the model admit the infinitely dimensional Lie group of equivalence. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 50, No. 2, pp. 24–36, March–April, 2009.     

Numerical solutions of the shallow water equations with discontinuous bed topography

A simple scheme is developed for treatment of vertical bed topography in shallow water flows. The effect of the vertical step on flows is modelled with the shallow water equations including local energy loss terms. The bed elevation is denoted with z_b^‐ for the left and z_b⁺ for the right values at each grid point, hence exactly representing a discontinuity in the bed topography. The surface gradient method (SGM) is generalized to reconstruct water depths at cell interfaces involving a vertical step so that the fluxes at the cell interfaces can accurately be calculated with a Riemann solver. The scheme is verified by predicting a surge crossing a step, a tidal flow over a step and dam‐break flows on wet/dry beds. The results have shown good agreements compared with analytical solutions and available experimental data. The scheme is efficient, robust, and may be used for practical flow calculations. Copyright © 2002 John Wiley & Sons, Ltd.     

Symmetries and exact solutions of the shallow water equations for a two-dimensional shear flow

This paper considers nonlinear equations describing the propagation of long waves in two-dimensional shear flow of a heavy ideal incompressible fluid with a free boundary. A nine-dimensional group of transformations admitted by the equations of motion is found by symmetry methods. Two-dimensional subgroups are used to find simpler integrodifferential submodels which define classes of exact solutions, some of which are integrated. New steady-state and unsteady rotationally symmetric solutions with a nontrivial velocity distribution along the depth are obtained. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 5, pp. 41–54, September–October, 2008.     

Some families of exact solutions of the equations of two-dimensional shallow water theory

Results are presented of a study of the exact solutions of the equations of two-dimensional unsteady and steady shallow water theory, based on the group properties of these equations. The first part presents the group properties of the equations in question; the second part presents the invariant solutions of these equations.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, Vol. 10, No. 6, pp. 62–71, November–December, 1969.The author wishes to thank L. V. Ovsyannikov and N. Kh. Ibragimov for valuable guidance in carrying out this study.     

Discontinuous solutions of the Navier-Stokes equations for compressible flow

We prove local existence and study properties of discontinuous solutions of the Navier-Stokes equations for one-dimensional, compressible, nonisentropic flow. We assume that, modulo a step function, the initial data is in L² and the initial velocity and density are in the space BV. We show that the velocity and the temperature become smoothed out in positive time, and that discontinuities in the density, pressure, and gradients of the velocity and temperature persist for all time. We also show that for stable gases these discontinuities decay exponentially in time, more rapidly for smaller viscosities.     

Backlund transformation and exact solutions for whitham-broer-kaup equations in shallow water

I.IntroductionTodescribethepropagationofshallowwaterwave,manywell-knowncompletelyintegrablemodelsareintroduced,suchasKdVequatioll,Boussinesqequation,K-Pequation,WBKequation,etc.UnderBoussinesqapproximation,Whitham,BroerandKaupl"2'3]obtainednon-lilleurWBKequationwhereu=u(x,l)isthefieldofhorizontalvelocity;v=v(x,t)istheheightthatdeviatefromequilibriumpositionofliquid;a,gareconstantsthatrepresentdift'erentdispersivepower.TheEqs.(1.I),(l.2)areverygoodmodelstodescl.ibedispersivewave.Ifa=0,P…     

Discontinuous solutions of boundary layer equations with a positive pressure gradient

Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 53–62, September–October, 1991.     

Relaxation schemes for the shallow water equations

We present a class of first and second order in space and time relaxation schemes for the shallow water (SW) equations. A new approach of incorporating the geometrical source term in the relaxation model is also presented. The schemes are based on classical relaxation models combined with Runge–Kutta time stepping mechanisms. Numerical results are presented for several benchmark test problems with or without the source term present. Copyright © 2003 John Wiley & Sons, Ltd.     

Numerical integration of the shallow water equations over a sloping shelf

A finite-difference method is described for the numerical integration of the one-dimensional shallow water equations over a sloping shelf that allows for a continuously moving shoreline. An application of the technique is made to the propagation of non-breaking waves towards the shoreline. The results of the computation are compared with an evaluation based upon an exact analytical treatment of the non-linear equations. Excellent agreement is found for both tsunami and tidal scale oscillations.     

Performance of finite volume solutions to the shallow water equations with shock‐capturing schemes

Numerical methods have become well established as tools for solving problems in hydraulic engineering. In recent years the finite volume method (FVM) with shock capturing capabilities has come to the fore because of its suitability for modelling a variety of types of flow; subcritical and supercritical; steady and unsteady; continuous and discontinuous and its ability to handle complex topography easily. This paper is an assessment and comparison of the performance of finite volume solutions to the shallow water equations with the Riemann solvers; the Osher, HLL, HLLC, flux difference splitting (Roe) and flux vector splitting. In this paper implementation of the FVM including the Riemann solvers, slope limiters and methods used for achieving second order accuracy are described explicitly step by step. The performance of the numerical methods has been investigated by applying them to a number of examples from the literature, providing both comparison of the schemes with each other and with published results. The assessment of each method is based on five criteria; ease of implementation, accuracy, applicability, numerical stability and simulation time. Finally, results, discussion, conclusions and recommendations for further work are presented. Copyright © 2002 John Wiley & Sons, Ltd.     

Multidimensional upwind schemes for the shallow water equations

A multidimensional discretisation of the shallow water equations governing unsteady free-surface flow is proposed. The method, based on a residual distribution discretisation, relies on a characteristic eigenvector decomposition of each cell residual, and the use of appropriate distribution schemes. For uncoupled equations, multidimensional convection schemes on compact stencils are used, while for coupled equations, either system distribution schemes such as the Lax–Wendroff scheme or scalar schemes may be used. For steady subcritical flows, the equations can be partially diagonalised into a purely convective equation of hyperbolic nature, and a set of coupled equations of elliptic nature. The multidimensional discretisation, which is second-order-accurate at steady state, is shown to be superior to the standard Lax–Wendroff discretisation. For steady supercritical flows, the equations can be fully diagonalised into a set of convective equations corresponding to the steady state characteristics. Discontinuities such as hydraulic jumps, are captured in a sharp and non-oscillatory way. For unsteady flows, the characteristic equations remain coupled. An appropriate treatment of the coupling terms allows the discretisation of these equations at the scalar level. Although presently only first-order-accurate in space and time, the classical dam-break problem demonstrates the validity of the approach. © 1998 John Wiley & Sons, Ltd.     

Singular solutions of equations of shallow shells for a concentrated tangential load

Journal of Applied Mechanics and Technical Physics -     

The influence of normal flow boundary conditions on spurious modes in finite element solutions to the shallow water equations

Finite element solutions of the primitive equation (PE) form of the shallow water equations are notorious for the severe spurious 2Δx modes which appear. Wave equation (WE) solutions do not exhibit these numerical modes. In this paper we show that the severe spurious modes in PE solutions are strongly influenced by essential normal flow boundary conditions in the coupled continuity-momentum system of equations. This is demonstrated through numerical examples that avoid the use of essential normal flow boundary conditions either by specifying elevation values over the entire boundary or by implementing natural flow boundary conditions in the weak weighted residual form of the continuity equation. Results from a series of convergence tests show that PE solutions are of nearly the same quality as WE solutions when spurious modes are suppressed by alternative specification of the boundary conditions. Network intercomparisons indicate that varying nodal support does not excite spurious modes in a solution, although it does enhance the spurious modes introduced when an essential normal flow boundary condition is used. Dispersion analysis of discrete equations for interior and boundary nodes offers an explanation of the observed solution behaviour. For certain PE algorithms a mixed situation can arise where the boundary nodes exhibit a monotonic (noise-free) dispersion relationship and the interior nodes exhibit a folded (noisy) dispersion relationship. We have found that the mixed situation occurs when all boundary nodes are specified elevation nodes (which are enforced as essential conditions in the continuity equation) or when specified flow boundary nodes are treated as natural boundary conditions in the continuity equation. In either case the effect is to generate a solution that is essentially free of noise. Apparently, the monotonic dispersion behaviour at the boundaries suppresses the otherwise noisy behaviour caused by the folded dispersion relation on the interior.     

Bidirectional Whitham equations as models of waves on shallow water

Hammack & Segur (1978) conducted a series of surface water-wave experiments in which the evolution of long waves of depression was measured and studied. This present work compares time series from these experiments with predictions from numerical simulations of the KdV, Serre, and five unidirectional and bidirectional Whitham-type equations. These comparisons show that the most accurate predictions come from models that contain accurate reproductions of the Euler phase velocity, sufficient nonlinearity, and surface tension effects. The main goal of this paper is to determine how accurately the bidirectional Whitham equations can model data from real-world experiments of waves on shallow water. Most interestingly, the unidirectional Whitham equation including surface tension provides the most accurate predictions for these experiments. If the initial horizontal velocities are assumed to be zero (the velocities were not measured in the experiments), the three bidirectional Whitham systems examined herein provide approximations that are significantly more accurate than the KdV and Serre equations. However, they are not as accurate as predictions obtained from the unidirectional Whitham equation.