ALBERT

Description: Tsunami modeling commonly accepts the shallow-water system as governing equations where the major difficulty is the correct treatment of the non-conservative 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 non-conservative 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. This article is protected by copyright. All rights reserved.