Frédéric Dias

Denys Dutykh, Raphaël Poncet, Frédéric Dias

A novel tool for tsunami wave modelling is presented. This tool has the potential of being used for operational purposes: indeed, the numerical code \VOLNA is able to handle the complete life-cycle of a tsunami (generation, propagation and run-up along the coast). The algorithm works on unstructured triangular meshes and thus can be run in arbitrary complex domains. This paper contains the detailed description of the finite volume scheme implemented in the code. The numerical treatment of the wet/dry transition is explained. This point is crucial for accurate run-up/run-down computations. Most existing tsunami codes use semi-empirical techniques at this stage, which are not always sufficient for tsunami hazard mitigation. Indeed the decision to evacuate inhabitants is based on inundation maps which are produced with this type of numerical tools. We present several realistic test cases that partially validate our algorithm. Comparisons with analytical solutions and experimental data are performed. Finally the main conclusions are outlined and the perspectives for future research presented.

Denys Dutykh, Dimitrios Mitsotakis, Xavier Gardeil et al.

The present study is devoted to the problem of tsunami wave generation. The main goal of this work is two-fold. First of all, we propose a simple and computationally inexpensive model for the description of the sea bed displacement during an underwater earthquake, based on the finite fault solution for the slip distribution under some assumptions on the dynamics of the rupturing process. Once the bottom motion is reconstructed, we study waves induced on the free surface of the ocean. For this purpose we consider three different models approximating the Euler equations of the water wave theory. Namely, we use the linearized Euler equations (we are in fact solving the Cauchy-Poisson problem), a Boussinesq system and a novel weakly nonlinear model. An intercomparison of these approaches is performed. The developments of the present study are illustrated on the 17 July 2006 Java event, where an underwater earthquake of magnitude 7.7 generated a tsunami that inundated the southern coast of Java.

Denys Dutykh, Frédéric Dias

The present article is devoted to the influence of sediment layers on the process of tsunami generation. The main scope here is to demonstrate and especially quantify the effect of sedimentation on vertical displacements of the seabed due to an underwater earthquake. The fault is modelled as a Volterra-type dislocation in an elastic half-space. The elastodynamics equations are integrated with a finite element method. A comparison between two cases is performed. The first one corresponds to the classical situation of an elastic homogeneous and isotropic half-space, which is traditionally used for the generation of tsunamis. The second test case takes into account the presence of a sediment layer separating the oceanic column from the hard rock. Some important differences are revealed. We conjecture that deformations in the generation region may be amplified by sedimentary deposits, at least for some parameter values. The mechanism of amplification is studied through careful numerical simulations.

Frédéric Dias, Denys Dutykh, Jean-Michel Ghidaglia

In the study of ocean wave impact on structures, one often uses Froude scaling since the dominant force is gravity. However the presence of trapped or entrained air in the water can significantly modify wave impacts. When air is entrained in water in the form of small bubbles, the acoustic properties in the water change dramatically and for example the speed of sound in the mixture is much smaller than in pure water, and even smaller than in pure air. While some work has been done to study small-amplitude disturbances in such mixtures, little work has been done on large disturbances in air-water mixtures. We propose a basic two-fluid model in which both fluids share the same velocities. It is shown that this model can successfully mimic water wave impacts on coastal structures. Even though this is a model without interface, waves can occur. Their dispersion relation is discussed and the formal limit of pure phases (interfacial waves) is considered. The governing equations are discretized by a second-order finite volume method. Numerical results are presented. It is shown that this basic model can be used to study violent aerated flows, especially by providing fast qualitative estimates.

Frédéric Chardard, Frédéric Dias, Thomas J. Bridges

When solitary waves are characterized as homoclinic orbits of a finite-dimensional Hamiltonian system, they have an integer-valued topological invariant, the Maslov index. We are interested in developing a robust numerical algorithm to compute the Maslov index, to understand its properties, and to study the implications for the stability of solitary waves. The algorithms reported here are developed in the exterior algebra representation, which leads to a robust and fast algorithm with some novel properties. We use two different representations for the Maslov index, one based on an intersection index and one based on approximating the homoclinic orbit by a sequence of periodic orbits. New results on the Maslov index for solitary wave solutions of reaction-diffusion equations, the fifth-order Korteweg-De Vries equation, and the longwave-shortwave resonance equations are presented. Part 1 considers the case of four-dimensional phase space, and Part 2 considers the case of $2n-$dimensional phase space with $n>2$.