No Arabic abstract
In this paper we analyze the stability of equilibrium manifolds of hyperbolic shallow water moment equations. Shallow water moment equations describe shallow flows for complex velocity profiles which vary in vertical direction and the models can be seen as extensions of the standard shallow water equations. Equilibrium stability is an important property of balance laws that determines the linear stability of solutions in the vicinity of equilibrium manifolds and it is seen as a necessary condition for stable numerical solutions. After an analysis of the hyperbolic structure of the models, we identify three different stability manifolds based on three different limits of the right-hand side friction term, which physically correspond to water-at-rest, constant-velocity, and bottom-at-rest velocity profiles. The stability analysis then shows that the structural stability conditions are fulfilled for the water-at-rest equilibrium and the constant-velocity equilibrium. However, the bottom-at-rest equilibrium can lead to instable modes depending on the velocity profile. Relaxation towards the respective equilibrium manifolds is investigated numerically for different models.
The regularisation of nonlinear hyperbolic conservation laws has been a problem of great importance for achieving uniqueness of weak solutions and also for accurate numerical simulations. In a recent work, the first two authors proposed a so-called Hamiltonian regularisation for nonlinear shallow water and isentropic Euler equations. The characteristic property of this method is that the regularisation of solutions is achieved without adding any artificial dissipation or ispersion. The regularised system possesses a Hamiltonian structure and, thus, formally preserves the corresponding energy functional. In the present article we generalise this approach to shallow water waves over general, possibly time-dependent, bottoms. The proposed system is solved numerically with continuous Galerkin method and its solutions are compared with the analogous solutions of the classical shallow water and dispersive Serre-Green-Naghdi equations. The numerical results confirm the absence of dispersive and dissipative effects in presence of bathymetry variations.
We derive boundary conditions and estimates based on the energy and entropy analysis of systems of the nonlinear shallow water equations in two spatial dimensions. It is shown that the energy method provides more details, but is fully consistent with the entropy analysis. The details brought forward by the nonlinear energy analysis allow us to pinpoint where the difference between the linear and nonlinear analysis originate. We find that the result from the linear analysis does not necessarily hold in the nonlinear case. The nonlinear analysis leads in general to a different minimal number of boundary conditions compared with the linear analysis. In particular, and contrary to the linear case, the magnitude of the flow does not influence the number of required boundary conditions.
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.
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 develop an adjoint approach for recovering the topographical function included in the source term of one-dimensional hyperbolic balance laws. We focus on a specific system, namely the shallow water equations, in an effort to recover the riverbed topography. The novelty of this work is the ability to robustly recover the bottom topography using only noisy boundary data from one measurement event and the inclusion of two regularization terms in the iterative update scheme. The adjoint scheme is determined from a linearization of the forward system and is used to compute the gradient of a cost function. The bottom topography function is recovered through an iterative process given by a three-operator splitting method which allows the feasibility to include two regularization terms. Numerous numerical tests demonstrate the robustness of the method regardless of the choice of initial guess and in the presence of discontinuities in the solution of the forward problem.