On the relevance of the dam break problem in the context of nonlinear shallow water equations

The classical dam break problem has become the de facto standard in validating the Nonlinear Shallow Water Equations (NSWE) solvers. Moreover, the NSWE are widely used for flooding simulations. While applied mathematics community is essentially focused on developing new numerical schemes, we tried to examine the validity of the mathematical model under consideration. The main purpose of this study is to check the pertinence of the NSWE for flooding processes. From the mathematical point of view, the answer is not obvious since all derivation procedures assumes the total water depth positivity. We performed a comparison between the two-fluid Navier-Stokes simulations and the NSWE solved analytically and numerically. Several conclusions are drawn out and perspectives for future research are outlined.Comment: 20 pages, 15 figures. Accepted to Discrete and Continuous Dynamical Systems. Other author's papers can be downloaded at http://www.lama.univ-savoie.fr/~dutyk

Tsunami generation by dynamic displacement of sea bed due to dip-slip faulting

In classical tsunami-generation techniques, one neglects the dynamic sea bed displacement resulting from fracturing of a seismic fault. The present study takes into account these dynamic effects. Earth's crust is assumed to be a Kelvin-Voigt material. The seismic source is assumed to be a dislocation in a viscoelastic medium. The fluid motion is described by the classical nonlinear shallow water equations (NSWE) with time-dependent bathymetry. The viscoelastodynamic equations are solved by a finite-element method and the NSWE by a finite-volume scheme. A comparison between static and dynamic tsunami-generation approaches is performed. The results of the numerical computations show differences between the two approaches and the dynamic effects could explain the complicated shapes of tsunami wave trains.Comment: 16 pages, 10 figures, Accepted to Mathematics and Computers in Simulation. Other author's papers can be downloaded at http://www.cmla.ens-cachan.fr/~dutyk

Tsunami generation by paddle motion and its interaction with a beach: Lagrangian modelling and experiment

A 2D Lagrangian numerical wave model is presented and validated against a set of physical wave-flume experiments on interaction of tsunami waves with a sloping beach. An iterative methodology is proposed and applied for experimental generation of tsunami-like waves using a piston-type wavemaker with spectral control. Three distinct types of wave interaction with the beach are observed with forming of plunging or collapsing breaking waves. The Lagrangian model demonstrates good agreement with experiments. It proves to be efficient in modelling both wave propagation along the flume and initial stages of strongly non-linear wave interaction with a beach involving plunging breaking. Predictions of wave runup are in agreement with both experimental results and the theoretical runup law

Tsunami Propagation from the Open Sea to the Coast

The tsunami generated by the Great East Japan Earthquake (GEJE) caused serious damage to the coastal areas of the Tohoku district. Numerical simulations are used to predict damage caused by tsunamis. Shallow-water equations are generally used in numerical simulations of tsunami propagation from the open sea to the coast. This research focuses on viscous shallow-water equations and attempts to generate a computational method using finite-element techniques based on the previous investigations of Kanayama and Ohtsuka (1978). First, the viscous shallow-water equation system is derived from the Navier-Stokes equations, based on the assumption of hydrostatic pressure in the direction of gravity. The derived equations have the horizontal viscosity term hereditary from the original Navier-Stokes equations. Next, a numerical finite element scheme is shown. Finally, tsunami simulations of Tohoku-Oki are shown using the above-mentioned approach. Our main concern in this chapter is how to set the boundary condition on the open boundary

Spectral/hp element methods: recent developments, applications, and perspectives

The spectral/hp element method combines the geometric flexibility of the classical h-type finite element technique with the desirable numerical properties of spectral methods, employing high-degree piecewise polynomial basis functions on coarse finite element-type meshes. The spatial approximation is based upon orthogonal polynomials, such as Legendre or Chebychev polynomials, modified to accommodate C0-continuous expansions. Computationally and theoretically, by increasing the polynomial order p, high-precision solutions and fast convergence can be obtained and, in particular, under certain regularity assumptions an exponential reduction in approximation error between numerical and exact solutions can be achieved. This method has now been applied in many simulation studies of both fundamental and practical engineering flows. This paper briefly describes the formulation of the spectral/hp element method and provides an overview of its application to computational fluid dynamics. In particular, it focuses on the use the spectral/hp element method in transitional flows and ocean engineering. Finally, some of the major challenges to be overcome in order to use the spectral/hp element method in more complex science and engineering applications are discussed

Operational tsunami modelling with TsunAWI – recent developments and applications

In this article, the tsunami model TsunAWI (Alfred Wegener Institute) and its application for hindcasts, inundation studies, and the operation of the tsunami scenario repository for the Indonesian tsunami early warning system are presented. TsunAWI was developed in the framework of the German-Indonesian Tsunami Early Warning System (GITEWS) and simulates all stages of a tsunami from the origin and the propagation in the ocean to the arrival at the coast and the inundation on land. It solves the non-linear shallow water equations on an unstructured finite element grid that allows to change the resolution seamlessly between a coarse grid in the deep ocean and a fine representation of coastal structures. During the GITEWS project and the following maintenance phase, TsunAWI and a framework of pre- and postprocessing routines was developed step by step to provide fast computation of enhanced model physics and to deliver high quality results

An entropy stable discontinuous Galerkin method for the shallow water equations on curvilinear meshes with wet/dry fronts accelerated by GPUs

We extend the entropy stable high order nodal discontinuous Galerkin spectral element approximation for the non-linear two dimensional shallow water equations presented by Wintermeyer et al. [N. Wintermeyer, A. R. Winters, G. J. Gassner, and D. A. Kopriva. An entropy stable nodal discontinuous Galerkin method for the two dimensional shallow water equations on unstructured curvilinear meshes with discontinuous bathymetry. Journal of Computational Physics, 340:200-242, 2017] with a shock capturing technique and a positivity preservation capability to handle dry areas. The scheme preserves the entropy inequality, is well-balanced and works on unstructured, possibly curved, quadrilateral meshes. For the shock capturing, we introduce an artificial viscosity to the equations and prove that the numerical scheme remains entropy stable. We add a positivity preserving limiter to guarantee non-negative water heights as long as the mean water height is non-negative. We prove that non-negative mean water heights are guaranteed under a certain additional time step restriction for the entropy stable numerical interface flux. We implement the method on GPU architectures using the abstract language OCCA, a unified approach to multi-threading languages. We show that the entropy stable scheme is well suited to GPUs as the necessary extra calculations do not negatively impact the runtime up to reasonably high polynomial degrees (around $N=7$). We provide numerical examples that challenge the shock capturing and positivity properties of our scheme to verify our theoretical findings