A primitive-variable Riemann method for solution of the shallow water equations with wetting and drying

A Riemann flux that uses primitive variables rather than conserved variables is developed for the shallow water equations with nonuniform bathymetry. This primitive-variable flux is both conservative and well behaved at zero depth. The unstructured finite-volume discretization used is suitable for highly nonuniform grids that provide resolution of complex geometries and localized flow structures. A source-term discretization is derived for nonuniform bottom that balances the discrete flux integral both for still water and in dry regions. This primitive-variable formulation is uniformly valid in wet and dry regions with embedded wetting and drying fronts. A fully nonlinear implicit scheme and both nonlinear and time-linearized explicit schemes are developed for the time integration. The implicit scheme is solved by a parallel Newton-iterative algorithm with numerically computed flux Jacobians. A concise treatment of characteristic-variable boundary conditions with source terms is also given. Computed results obtained for the one-dimensional dam break on wet and dry beds and for normal-mode oscillations in a circular parabolic basin are in very close agreement with the analytical solutions. Other results for a forced breaking wave with friction interacting with a sloped bottom demonstrate a complex wave motion with wetting, drying and multiple interacting wave fronts. Finally, a highly nonuniform, coastline-conforming unstructured grid is used to demonstrate an unsteady simulation that models an artificial coastal flooding due to a forced wave entering the Gulf of Mexico.     

Well-balanced and energy stable schemes for the shallow water equations with discontinuous topography

We consider the shallow water equations with non-flat bottom topography. The smooth solutions of these equations are energy conservative, whereas weak solutions are energy stable. The equations possess interesting steady states of lake at rest as well as moving equilibrium states. We design energy conservative finite volume schemes which preserve (i) the lake at rest steady state in both one and two space dimensions, and (ii) one-dimensional moving equilibrium states. Suitable energy stable numerical diffusion operators, based on energy and equilibrium variables, are designed to preserve these two types of steady states. Several numerical experiments illustrating the robustness of the energy preserving and energy stable well-balanced schemes are presented.     

ADER schemes for the shallow water equations in channel with irregular bottom elevation

This paper deals with the construction of high-order ADER numerical schemes for solving the one-dimensional shallow water equations with variable bed elevation. The non-linear version of the schemes is based on ENO reconstructions. The governing equations are expressed in terms of total water height, instead of total water depth, and discharge. The ENO polynomial interpolation procedure is also applied to represent the variable bottom elevation. ADER schemes of up to fifth order of accuracy in space and time for the advection and source terms are implemented and systematically assessed, with particular attention to their convergence rates. Non-oscillatory results are obtained for discontinuous solutions both for the steady and unsteady cases. The resulting schemes can be applied to solve realistic problems characterized by non-uniform bottom geometries.     

The classical shallow water equations: Symplectic geometry

The classical shallow water equations express the change with time of the height h and the velocity ν of a 1-dimensional fluid:

$\text{νξ}\text{νt}\text{+}\text{νξ}\text{νx}\text{+}\text{νh}\text{νx}\text{=0.}\text{νh}\text{νx}\text{+}\text{νhν}\text{νx}\text{=0}$

. They possess an infinite number of integrals of motion due to Benney [1973] and can be written in Hamiltonian form relative to a symplectic structure introduced by Manin [1978]. The present paper deals with their complete integrability up to the advent of shocks. This is proved in the small under an extra assumption satisfied by most height-velocity pairs: that $\text{h}\text{h′ = ± ν′}$ only at isolated points.     

Augmented Riemann solvers for the shallow water equations over variable topography with steady states and inundation

We present a class of augmented approximate Riemann solvers for the shallow water equations in the presence of a variable bottom surface. These belong to the class of simple approximate solvers that use a set of propagating jump discontinuities, or waves, to approximate the true Riemann solution. Typically, a simple solver for a system of m conservation laws uses m such discontinuities. We present a four wave solver for use with the the shallow water equations—a system of two equations in one dimension. The solver is based on a decomposition of an augmented solution vector—the depth, momentum as well as momentum flux and bottom surface. By decomposing these four variables into four waves the solver is endowed with several desirable properties simultaneously. This solver is well-balanced: it maintains a large class of steady states by the use of a properly defined steady state wave—a stationary jump discontinuity in the Riemann solution that acts as a source term. The form of this wave is introduced and described in detail. The solver also maintains depth non-negativity and extends naturally to Riemann problems with an initial dry state. These are important properties for applications with steady states and inundation, such as tsunami and flood modeling. Implementing the solver with LeVeque’s wave propagation algorithm [R.J. LeVeque, Wave propagation algorithms for multi-dimensional hyperbolic systems, J. Comput. Phys. 131 (1997) 327–335] is also described. Several numerical simulations are shown, including a test problem for tsunami modeling.     

Initial value problem solution of nonlinear shallow water-wave equations

The initial value problem solution of the nonlinear shallow water-wave equations is developed under initial waveforms with and without velocity. We present a solution method based on a hodograph-type transformation to reduce the nonlinear shallow water-wave equations into a second-order linear partial differential equation and we solve its initial value problem. The proposed solution method overcomes earlier limitation of small waveheights when the initial velocity is nonzero, and the definition of the initial conditions in the physical and transform spaces is consistent. Our solution not only allows for evaluation of differences in predictions when specifying an exact initial velocity based on nonlinear theory and its linear approximation, which has been controversial in geophysical practice, but also helps clarify the differences in runup observed during the 2004 and 2005 Sumatran tsunamigenic earthquakes.     

Exact solution of the inverse acoustic scattering problem for the one-dimensional density profile in cylindrical geometry

The inhomogeneous acoustic wave equation in cylindrically symmetric geometry, is reduced to a one-dimensional Schrödinger equation for time harmonic waves. The scattering potential so formed is uniquely reconstructed by the Gelfand-Levitan theory. In physical terms, the scattering potential is a function of the impedance profile, which is then recovered. If the velocity profile is known (or constant) the density profile, as a function of the radial coordinate, can be deduced. A practical experiment is envisaged whereby the back-scattered time sequence is measured in a pulse-echo mode, and the data are used to reconstruct the density profile.     

Parallel multilevel methods for implicit solution of shallow water equations with nonsmooth topography on the cubed-sphere

High resolution and scalable parallel algorithms for the shallow water equations on the sphere are very important for modeling the global climate. In this paper, we introduce and study some highly scalable multilevel domain decomposition methods for the fully implicit solution of the nonlinear shallow water equations discretized with a second-order well-balanced finite volume method on the cubed-sphere. With the fully implicit approach, the time step size is no longer limited by the stability condition, and with the multilevel preconditioners, good scalabilities are obtained on computers with a large number of processors. The investigation focuses on the use of semismooth inexact Newton method for the case with nonsmooth topography and the use of two- and three-level overlapping Schwarz methods with different order of discretizations for the preconditioning of the Jacobian systems. We test the proposed algorithm for several benchmark cases and show numerically that this approach converges well with smooth and nonsmooth bottom topography, and scales perfectly in terms of the strong scalability and reasonably well in terms of the weak scalability on machines with thousands and tens of thousands of processors.     

High-order discontinuous Galerkin schemes on general 2D manifolds applied to the shallow water equations

An innovating approach is proposed to solve vectorial conservation laws on curved manifolds using the discontinuous Galerkin method. This new approach combines the advantages of the usual approaches described in the literature. The vectorial fields are expressed in a unit non-orthogonal local tangent basis derived from the polynomial mapping of curvilinear triangle elements, while the convective flux functions are written is the usual 3D Cartesian coordinate system. The number of vectorial components is therefore minimum and the tangency constraint is naturally ensured, while the method remains robust and general since not relying on a particular parametrization of the manifold. The discontinuous Galerkin method is particularly well suited for this approach since there is no continuity requirement between elements for the tangent basis definition. The possible discontinuities of this basis are then taken into account in the Riemann solver on inter-element interfaces.     

Exact solution of general-relativity equations for vibratory collapse

In our previous article (cited as [A]) general requirements were formulated for the solutions of the general-relativity equations that follow from the geometry in the large and from the requirement that the origin be physically realizable [1]. In this article an exact self-consistent equation is described that takes full account of these requirements. The self-consistency is constant with respect to the volume density in the case of a dustlike ball [A17]. Although the self-consistent solution is of measure zero, it enables one to judge the properties of a wide class of solutions as shown by the experience of hydrodynamics.Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Fizika, No. 7, pp. 90–98, July, 1977.     

On numerical solution of the shallow water equations with chemical reactions on icosahedral geodesic grid

The shallow water equations coupled to the set of reaction–advection–diffusion equations are discretized on a geodesic icosahedral mesh using the finite volume technique. The method of solution of this coupled system is based on the principle of semi-discretization. The algorithm is mass conserving and stable for advection with the Courant numbers up to 2.7. The important part of the methodology is the optimization of the node positions in the icosahedral grid. It is shown that a slight adjustment of the mesh is instrumental in improving the accuracy of the numerical approximation. The convergence of the approximation of the differential operators is evaluated and compared to the data published in the literature. Numerical tests performed with the shallow water solver include two advection experiments, steady and unsteady zonal balanced flow, mountain flow, and the Rossby wave. The mountain flow and the Rossby wave cases are used to test the transport properties of the method in the case of both passive and reactive scalar fields. The investigation of essential numerical characteristics of the method is concluded by the simulation of an unstable zonal jet. The numerical simulation is performed using the set of shallow water equations without dissipation as well as with the viscosity term added to the momentum equation. Results show that the behavior of the model is consistent with both the literature published on the subject and the general empirical evidence.     

浅海周期起伏海底环境下的声传播

海底粗糙对水下声传播及水声探测等应用具有重要影响.利用黄海夏季典型海洋环境,分析了同时存在海底周期起伏和强温跃层条件下的声传播特性,结果表明:由于海底周期起伏的存在,对于低频(<1 kHz)、近程(10 km)的声信号,传播损失可增大5—30 d B.总结了声传播损失及脉冲到达结构随声源深度、海底起伏周期及起伏高度等因...     

Exact hierarchical solution of the damped oscillator problem

As a alternative to other exact approaches we present a procedure, in which the eigenvectors of a coupled oscillatory system are calculated in a hierarchically exact (N) manner. In the system a single oscillator is coupled to a dense sequence of background oscillators, which among themselves are not mutually coupled (Oscillatory Fano problem). The eigenvectors beeing known, all decay laws become computable. Two of them are given. It is strongly emphasized thatquite different decay laws evolve fordifferent initial conditions.Dedicated to B. Mühlschlegel on the occasion of his 60th birthday     

Surface-to-volume wave conversion in shallow water with a corrugated bottom

Acoustic transmission between points onshore or in very shallow water and points in deep water is strongly influenced by the shear rigidity of marine sediments, which control the parameters and the very existence of seismoacoustic surface waves. Previously, it was found that coupling between acoustic modes and the seismoacoustic surface waves is normally weak, although not negligible in the case of a gently sloping seafloor and soft sediments. In this paper, the previous work is extended by accounting for the small-scale roughness of the seafloor. The significant role of roughness in coupling between volume and surface waves is demonstrated. The combined effect of bottom topography, roughness, and wave attenuation in soft marine sediments on the sound propagation between points in shallow and deep water is discussed. Published in Russian in Akusticheskiĭ Zhurnal, 2008, Vol. 54, No. 3, pp. 400–407. The article was translated by the author.     

Acoustic tomography of shallow water with unknown relief of hard bottom

The possibility of implementing tomography in shallow water with unknown relief of bottom is considered. It is shown that the use of time delays of different modes at different frequencies allows for reconstruction of bottom relief and the sound speed profile in the water layer without any additional activities for separating the effects of reconstructed characteristics of shallow water on the received data.     

Exact one-body solution of the gauge-covariant gravitational field equations

Dirac has recently introduced a new theory of gravity which includes his large numbers hypothesis. The exact one-body vacuum solution consistent with multiplicative matter creation is obtained and discussed.Supported by a grant from the National Research Council of Canada.     

Exact solution of the AF coherent pulse propagation equations

The exact solution of the adiabatic following pulse propagation equations describing coherent propagation of a single-photon quasi-resonant short pulse of arbitrary shape and phase is reported.