Second‐order finite volume with hydrostatic reconstruction for tsunami simulation

Abstract Tsunami modeling commonly accepts the shallow water system as governing equations where the major difficulty is the correct treatment of the nonconservative term due to bathymetry variations. The finite volume method for solving the shallow water equations with such source terms has receive...

Full description

Bibliographic Details
Published in:Journal of Advances in Modeling Earth Systems
Main Authors: S. Clain, C. Reis, R. Costa, J. Figueiredo, M. A. Baptista, J. M. Miranda
Format: Article in Journal/Newspaper
Language:English
Published: American Geophysical Union (AGU) 2016
Subjects:
Online Access:https://doi.org/10.1002/2015MS000603
https://doaj.org/article/afa2344e867e4b449b878a697d321cda
Description
Summary:Abstract Tsunami modeling commonly accepts the shallow water system as governing equations where the major difficulty is the correct treatment of the nonconservative term due to bathymetry variations. The finite volume method for solving the shallow water equations with such source terms has received great attention in the two last decades. The built‐in conservation property, the capacity to correctly treat discontinuities, and the ability to handle complex bathymetry configurations preserving some steady state configurations (well‐balanced scheme) make the method very efficient. Nevertheless, it is still a challenge to build an efficient numerical scheme, with very few numerical artifacts (e.g., small numerical diffusion, correct propagation of the discontinuities, accuracy, and robustness), to be used in an operational environment, and that is able to better capture the dynamics of the wet‐dry interface and the physical phenomena that occur in the inundation area. In the first part of this paper, we present a new second‐order finite volume code. The code is developed for the shallow water equations with a nonconservative term based on the hydrostatic reconstruction technology to achieve a well‐balanced scheme and an adequate dry/wet interface treatment. A detailed presentation of the numerical method is given. In the second part of the paper, we highlight the advantages of the new numerical technique. We benchmark the numerical code against analytical, experimental, and field results to assess the robustness and the accuracy of the numerical code. Finally, we use the 28 February 1969 North East Atlantic tsunami to check the performance of the code with real data.