No Arabic abstract
We formulate a new approach to solving the initial value problem of the shallow water-wave equations utilizing the famous Carrier-Greenspan transformation [G. Carrier and H. Greenspan, J. Fluid Mech. 01, 97 (1957)]. We use a Taylor series approximation to deal with the difficulty associated with the initial conditions given on a curve in the transformed space. This extends earlier solutions to waves with near shore initial conditions, large initial velocities, and in more complex U-shaped bathymetries; and allows verification of tsunami wave inundation models in a more realistic 2-D setting.
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.
A formulation of the shallow water equations adapted to general complex terrains is proposed. Its derivation starts from the observation that the typical approach of depth integrating the Navier-Stokes equations along the direction of gravity forces is not exact in the general case of a tilted curved bottom. We claim that an integration path that better adapts to the shallow water hypotheses follows the cross-flow surface, i.e., a surface that is normal to the velocity field at any point of the domain. Because of the implicitness of this definition, we approximate this cross-flow path by performing depth integration along a local direction normal to the bottom surface, and propose a rigorous derivation of this approximation and its numerical solution as an essential step for the future development of the full cross-flow integration procedure. We start by defining a local coordinate system, anchored on the bottom surface to derive a covariant form of the Navier-Stokes equations. Depth integration along the local normals yields a covariant version of the shallow water equations, which is characterized by flux functions and source terms that vary in space because of the surface metric coefficients and related derivatives. The proposed model is discretized with a first order FORCE-type Godunov Finite Volume scheme that allows implementation of spatially variable fluxes. We investigate the validity of our SW model and the effects of the bottom geometry by means of three synthetic test cases that exhibit non negligible slopes and surface curvatures. The results show the importance of taking into consideration bottom geometry even for relatively mild and slowly varying curvatures.
We present an extensive numerical comparison of a family of balance models appropriate to the semi-geostrophic limit of the rotating shallow water equations, and derived by variational asymptotics in Oliver (2006) for small Rossby numbers ${mathrm{Ro}}$. This family of generalized large-scale semi-geostrophic (GLSG) models contains the $L_1$-model introduced by Salmon (1983) as a special case. We use these models to produce balanced initial states for the full shallow water equations. We then numerically investigate how well these models capture the dynamics of an initially balanced shallow water flow. It is shown that, whereas the $L_1$-member of the GLSG family is able to reproduce the balanced dynamics of the full shallow water equations on time scales of ${mathcal{O}}(1/{mathrm{Ro}})$ very well, all other members develop significant unphysical high wavenumber contributions in the ageostrophic vorticity which spoil the dynamics.
We present a broadband waveguide for water waves obtained through mere manipulation of seabed properties and without any need for sidewalls. Specifically, we show that a viscoelastic seabed results in a modified effective gravity term in the governing equations of water waves, which provides a generic broadband mechanism to control oceanic wave energy and enables confining surface waves inside a long narrow path without sidewalls. Our findings have promising applications in guiding and steering waves for oceanic wave energy farms or protecting shorelines.
The shallow water equations (SWE) are a widely used model for the propagation of surface waves on the oceans. We consider the problem of optimally determining the initial conditions for the one-dimensional SWE in an unbounded domain from a small set of observations of the sea surface height. In the linear case we prove a theorem that gives sufficient conditions for convergence to the true initial conditions. At least two observation points must be used and at least one pair of observation points must be spaced more closely than half the effective minimum wavelength of the energy spectrum of the initial conditions. This result also applies to the linear wave equation. Our analysis is confirmed by numerical experiments for both the linear and nonlinear SWE data assimilation problems. These results show that convergence rates improve with increasing numbers of observation points and that at least three observation points are required for the practically useful results. Better results are obtained for the nonlinear equations provided more than two observation points are used. This paper is a first step in understanding the conditions for observability of the SWE for small numbers of observation points in more physically realistic settings.