IFCP Riemann solver: Application to tsunami modelling using GPUs

Abstract

In this work, we present a simplified two-layer model of Savage-Hutter type to simulate tsunamis generated by landslides (see (Fern´andez et al. 2008)). A layer composed of fluidized granular material is assumed to flow within an upper layer composed of an inviscid fluid (e.g. water). The sediment layer ismodelled by a Savage-Hutter type model where buoyancy effects have been considered. The system is discretized using IFCP finite volume scheme. The first order IFCP scheme was introduced in (Fern´andez et al. 2011) and it is constructed by using a suitable decomposition of a Roe matrix by means of a parabolic viscosity matrix, that captures information of the intermediate fields (Intermediate Field Capturing Parabola). Its extension to high order and two-dimensional domains is straightforward. To conclude, some numerical examples are presented