Numerical Simulation on a Globally Spherical Geometry

Gayaz Khakimzyanov, Denys Dutykh, Zinaida Fedotova, Oleg Gusev

Research output: Chapter in Book/Report/Conference proceedingChapterpeer-review


In this chapter we focus on numerical aspects, while the model derivation was described in Chap. 3. The algorithm we propose is based on the splitting approach. Namely, equations are decomposed on a uniform elliptic equation for the dispersive pressure component and a hyperbolic part of shallow water equations (on a sphere) with source terms. This algorithm is implemented as a two-step predictor–corrector scheme. On every step we solve separately elliptic and hyperbolic problems. Then, the performance of this algorithm is illustrated on model idealized situations with even bottom, where we estimate the influence of sphericity and rotation effects on dispersive wave propagation. The dispersive effects are quantified depending on the propagation distance over the sphere and on the linear extent of generation region. Finally, the numerical method is applied to a couple of real-world events. Namely, we undertake simulations of the Bulgarian 2007 and Chilean 2010 tsunamis. Whenever the data is available, our computational results are confronted with real measurements.

Original languageBritish English
Title of host publicationLecture Notes in Geosystems Mathematics and Computing
PublisherSpringer Nature
Number of pages47
StatePublished - 2020

Publication series

NameLecture Notes in Geosystems Mathematics and Computing
ISSN (Electronic)2512-3211


Dive into the research topics of 'Numerical Simulation on a Globally Spherical Geometry'. Together they form a unique fingerprint.

Cite this